请教关于“实序列快速傅里叶变换”（rfft）程序的问题～

[rabbitjerry]

出自<<数字信号处理C语言程序集>>，但运算结果与书中给的结果不符合，估计书中给的这个程序代码有些错误～哪位高手能帮忙看看错误在哪呢？
实序列快速傅里叶变换：
#include "stdafx.h"
#include "math.h"
void rfft(int n, double x[]);

int main(int argc, char* argv[])
{
int i,n;
double x[64];
n=64;
for(i=0;i<10;i++) x[i]=0;
for(i=0;i<n;i++) x[i]=exp(-(i-10)/15.0)*sin(6.2831853*(i-10)/16.0);
rfft(n,x);
printf("\n DISCRETE FOURIER TRANSFORM\n");
printf(" %10.7f",x[0]);
printf("\n");
printf(" %10.7f + J%10.7f\n",x[1],x[n-1]);
for(i=2;i<n/2;i+=2)
{
printf(" %10.7f + J%10.7f\n",x[i],x[n-i]);
printf(" %10.7f + J%10.7f\n",x[i+1],x[n-i-1]);
//printf("\n");
}
printf(" %10.7f",x[n/2]);
printf("\n");
printf(" %10.7f + J%10.7f\n",x[n/2-1],-x[n/2+1]);
for(i=2;i<n/2;i+=2)
{
printf(" %10.7f + J%10.7f\n",x[n/2-i],-x[n/2+i]);
printf(" %10.7f + J%10.7f\n",x[n/2-i-1],-x[n/2+i+1]);
//printf("\n");
}
return 0;
}

void rfft(int n, double x[]) //实序列快速傅里叶变换
{
int i,j,k,m,i1,i2,i3,i4,n1,n2,n4;
double a,e,cc,ss,xt,t1,t2;
for(j=1,i=1;i<16;i++)
{
m=i;
j= 2*j;
if(j==n) break;
}
n1=n-1;
for(j=0,i=0;i<n1;i++)
{
if(i<j)
{
xt=x[j];
x[j]=x[i];
x[i]=xt;
}
k=n/2;
while (k<(j+1)) {
j=j-k;
k=k/2;
}
j=j+k;
}
for(i=0;i<n;i+=2)
{
xt=x[i];
x[i]=xt+x[i+1];
x[i+1]=xt-x[i+1];
}
n2=1;
for(k=2;k<=m;k++)
{
n4=n2;
n2=2*n4;
n1=2*n2;
e=6.28318530718/n1;
for(i=0;i<n;i+=n1)
{
xt=x[i];
x[i]=xt+x[i+n2];
x[i+n2]=xt-x[i+n2];
x[i+n2+n4]=-x[i+n2+n4];//
a=e;
for(j=1;j<(n4-1);j++)
{
i1=i+j;
i2=i-j+n2;
i3=i+j+n2;
i4=i-j+n1;
cc=cos(a);
ss=sin(a);
a=a+e;
t1=cc*x[i3]+ss*x[i4];
t2=ss*x[i3]-cc*x[i4];
x[i4]=x[i2]-t2;
x[i3]=-x[i2]-t2;//
x[i2]=x[i1]-t1;
x[i1]=x[i1]+t1;
}
}
}
}

[m271262011]

dingdingding
mark

[m271262011]

dingdingding
mark

[大雷子老师]

电路中学过，已经忘记。。。

[o2cola]

看到傅里叶变换 进来了，
曾今极度反感老傅搞出来这么多东西

[guo5115]

HAO

[thrtc]

sinf
ding

[wuhui5125]

真够烦的，下载 东西要积分

[wuhui5125]

真够烦的，下载 东西要积分

[wuhui5125]

真够烦的，下载 东西要积分

[zjweii]

ddddddddddddddddddddddddd

[L_thread]

mask

[HUAXIANGDAGE]

单步跟踪试试！

[chrile]

我是出来打酱油的。飘过

[kingteng]

这么复杂。。。看我写的符合不？http://blog.chinaunix.net/u2/82750/showart_2208849.html

[勇敢的搬砖人]

漂亮美女在这里

[duxing87]

我自己试了一下，当输入1到8计算结果是对的，1到16就有问题，不知道怎么回事

[peijch02]

最好附上数学表达式，

[wangw89]

闲人路过，帮顶一下！

[jardearson]

拿分走人...

[luechen0912]

那本书里我也发现，关于rlfft程序调试的程序结果和书上有不一致的地方。