有没有DFT实数序列变换的源程序啊？

roboy

最好是1024点的，谢谢！

aspen

看看这个是不是你需要的
算法原理与matlab语言的fft和ifft函数同。（可详见matlab的help文件）。

1. #include "stdio.h"
   
2. #include "math.h"
   
3. #define PI 3.1415926
   
4. #define N 7
   
5. void mefft(float pr[],float pi[],int n,float fr[],float fi[])
   
6. {
   
7.         int i,k;
   
8.         for(k=0;k<n;k++)
   
9.         {
   
10.                 fr[k]=0;
    
11.                 fi[k]=0;
    
12.         }
    
13.         for(k=0;k<n;k++)
    
14.         {
    
15.                 for(i=0;i<n;i++)
    
16.                 {
    
17.                         fr[k]=fr[k]+pr[i]*cos(-2*PI*i*k/n)-pi[i]*sin(-2*PI*i*k/n);
    
18.                         fi[k]=fi[k]+pi[i]*cos(-2*PI*i*k/n)+pr[i]*sin(-2*PI*i*k/n);
    
19.                 }
    
20.        
    
21.         }
    
22. }
    
23. 
    
24. void meifft(float fr[],float fi[],int n,float pr[],float pi[])
    
25. {
    
26.         int i,j,k;
    
27.         for(i=0;i<n;i++)
    
28.         {
    
29.                 pr[i]=0;
    
30.                 pi[i]=0;
    
31.         }
    
32.         for(i=0;i<n;i++)
    
33.         {
    
34.                 for(k=0;k<n;k++)
    
35.                 {
    
36.                         pr[i]=pr[i]+(fr[k]*cos(2*PI*i*k/n)-fi[k]*sin(2*PI*i*k/n));
    
37.                         pi[i]=pi[i]+(fi[k]*cos(2*PI*i*k/n)+fr[k]*sin(2*PI*i*k/n));
    
38.                 }
    
39.                 pr[i]=pr[i]/n;
    
40.                 pi[i]=pi[i]/n;
    
41.         }
    
42. }
    
43. 
    
44. main()
    
45. {
    
46.         float fr[N],fi[N];
    
47.         int i,j;
    
48.         static float pr[N]={1,2,3,4,5,6,7};
    
49.         static float pi[N]={0,0,0,0,0,0,0};
    
50.         mefft(pr,pi,7,fr,fi);
    
51.         for(i=0;i<7;i++)
    
52.         {
    
53.                 printf("%f %f\n",fr[i],fi[i]);
    
54.         }
    
55.         meifft(fr,fi,7,pr,pi);
    
56.         for(i=0;i<7;i++)
    
57.         {
    
58.                 printf("%f %f\n",pr[i],pi[i]);
    
59.         }
    
60. }

复制代码

roboy

<P>程序我还没有试，但是大概看了一下就是它了。太谢谢了！</P>

aspen

不用客气，希望以后经常光临，大家共同交流提高。