t12211017 发表于 2007-12-15 21:37

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

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

francexuqu 发表于 2007-12-16 14:11

回复 #1 t12211017 的帖子

#define SWAP(a,b) tempr=(a);(a)=(b);(b)=tempr

fft( double* data, int nn, int isign)
// numerical recipies routine "four1.c"
{
      int n,mmax,m,j,istep,i;
      double wtemp,wr,wpr,wpi,wi,theta;
      double tempr,tempi;

      n=nn << 1;
      j=1;
      for (i=1;i<n;i+=2) {
                if (j > i) {
                        SWAP(data,data);
                        SWAP(data,data);
                }
                m=n >> 1;
                while (m >= 2 && j > m) {
                        j -= m;
                        m >>= 1;
                }
                j += m;
      }
      mmax=2;
      while (n > mmax) {
                istep=2*mmax;
                theta=6.28318530717959/(isign*mmax);
                wtemp=sin(0.5*theta);
                wpr = -2.0*wtemp*wtemp;
                wpi=sin(theta);
                wr=1.0;
                wi=0.0;
                for (m=1;m<mmax;m+=2) {
                        for (i=m;i<=n;i+=istep) {
                              j=i+mmax;
                              tempr=wr*data-wi*data;
                              tempi=wr*data+wi*data;
                              data=data-tempr;
                              data=data-tempi;
                              data += tempr;
                              data += tempi;
                        }
                        wr=(wtemp=wr)*wpr-wi*wpi+wr;
                        wi=wi*wpr+wtemp*wpi+wi;
                }
                mmax=istep;
      }
}

#undef SWAP

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

hyl2323 发表于 2007-12-19 16:21

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

风花雪月 发表于 2007-12-24 11:15

原帖由 t12211017 于 2007-12-15 21:37 发表 http://www.chinavib.com/forum/images/common/back.gif
不知能否给我传一份,谢谢大家

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

zz80 发表于 2008-1-26 15:03

学习了
页: [1]
查看完整版本: 求c++傅立叶变换的源程序