求c++傅立叶变换的源程序

t12211017

不知能否给我传一份,谢谢大家

francexuqu

1. #define SWAP(a,b) tempr=(a);(a)=(b);(b)=tempr
   
2. 
   
3. fft( double* data, int nn, int isign)
   
4. // numerical recipies routine "four1.c"
   
5. {
   
6.         int n,mmax,m,j,istep,i;
   
7.         double wtemp,wr,wpr,wpi,wi,theta;
   
8.         double tempr,tempi;
   
9. 
   
10.         n=nn << 1;
    
11.         j=1;
    
12.         for (i=1;i<n;i+=2) {
    
13.                 if (j > i) {
    
14.                         SWAP(data[j],data);
    
15.                         SWAP(data[j+1],data[i+1]);
    
16.                 }
    
17.                 m=n >> 1;
    
18.                 while (m >= 2 && j > m) {
    
19.                         j -= m;
    
20.                         m >>= 1;
    
21.                 }
    
22.                 j += m;
    
23.         }
    
24.         mmax=2;
    
25.         while (n > mmax) {
    
26.                 istep=2*mmax;
    
27.                 theta=6.28318530717959/(isign*mmax);
    
28.                 wtemp=sin(0.5*theta);
    
29.                 wpr = -2.0*wtemp*wtemp;
    
30.                 wpi=sin(theta);
    
31.                 wr=1.0;
    
32.                 wi=0.0;
    
33.                 for (m=1;m<mmax;m+=2) {
    
34.                         for (i=m;i<=n;i+=istep) {
    
35.                                 j=i+mmax;
    
36.                                 tempr=wr*data[j]-wi*data[j+1];
    
37.                                 tempi=wr*data[j+1]+wi*data[j];
    
38.                                 data[j]=data-tempr;
    
39.                                 data[j+1]=data[i+1]-tempi;
    
40.                                 data += tempr;
    
41.                                 data[i+1] += tempi;
    
42.                         }
    
43.                         wr=(wtemp=wr)*wpr-wi*wpi+wr;
    
44.                         wi=wi*wpr+wtemp*wpi+wi;
    
45.                 }
    
46.                 mmax=istep;
    
47.         }
    
48. }
    
49. 
    
50. #undef SWAP

复制代码

[ 本帖最后由 风花雪月 于 2007-12-24 11:14 编辑 ]

hyl2323

我在matlab版发了一个帖子，本人拙作。

风花雪月

原帖由 t12211017 于 2007-12-15 21:37 发表

                               
登录/注册后可看大图

不知能否给我传一份,谢谢大家

搜索本站，我记得在信号处理版有好几个

zz80

学习了