[交流] 【求助】帮忙看看程序！谢谢了【已完结】 已有2人参与

这是fft的程序 运行中出了问题 不知哪的错误 #include #include #include #define N 1000 typedef struct {         double real;     double img; } complex; void fft();          //快速傅里叶变换 void initw(); void change(); void add(complex,complex,complex *); void mul(complex,complex,complex *); void sub(complex,complex,complex *); void output(); complex x(N),*w; int size_x=0; double PI; int main() {         int i,method;         void fft();         void initw();         void change();     void add(complex,complex,complex *);     void mul(complex,complex,complex *);     void sub(complex,complex,complex *);     void divi(complex,complex,complex *);         PI=atan(1)*4;         printf("please input the size of x:\n";         scanf("%d",&size_x);         printf("please input the date in x(N):\n";         for(i=0;i         scanf("%lf%lf",&x.real,&x.img);         initw();         output();         return 0; } void fft() {         int i=0,j=0,k=0,l=0;         complex up,down,product;         change();         for(i=0;i         {                 l=1<                 for(j=0;j                 {                         for(k=0;k                         {                                 mul(x[j+k+l],w[size_x*k/2/l],&product);                                 add(x[j+k],product,&up);                                 sub(x[j+k],product,&down);                                 x[j+k]=up;                                 x[j+k+l]=down;                         }                 }         } } void initw() {         int i;         w=(complex *)malloc(sizeof(complex)*size_x);         for(i=0;i         {                 w.real=cos(2*PI/size_x*i);                 w.img=-1*sin(2*PI/size_x*i);         } } void change() {         complex temp;         unsigned short i=0,j=0,k=0;         double t;         for(i=0;i         {                 k=1;j=0;                 t=(log(size_x)/log(2));                 while((t--)>0)                 {j=j<<1;j|=(k&1);k=k>>1;}                 if(j>1)                 {temp=x;x=x[j];x[j]=temp;}         } } void output() {         int i;         printf("the result are as follows \n";         for(i=0;i         {                 printf("%.4f",x.real);                 if(x.img>=0.0001)                         printf("+%.4f\n",x.img);                 else if(fabs(x.img)<0.0001)                         printf("\n";                 else printf("%.4f\n",x.img);         } } void add(complex a,complex b,complex *c) {         c->real=a.real+b.real;         c->img=a.img+b.img; } void mul(complex a,complex b,complex *c) {         c->real=a.real*b.real-a.img*b.img;         c->img=a.real*b.img+a.img*b.real; } void sub(complex a,complex b,complex *c) {         c->real=a.real-b.real;         c->img=a.img-b.img; } [ Last edited by 余泽成 on 2010-6-2 at 16:23 ]

回复此楼