信号分析与处理,离散傅里叶变换与反变换

《信号分析与处理,离散傅里叶变换与反变换》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、#include"math.h"#include"stdio.h"#defineswap(a,b);T=(a);(a)=(b);(b)=TvoidIfft(floatA[],floatB[],unsignedM);main(){floatA[1024],B[1024];unsignedlongN,I;unsignedM;printf("inputM=");scanf("%d",&M);N=1<

2、"inputB[]:");for(I=0;I

3、Tr,WTi,theta,Tr,Ti;N=1<I){floatT;swap(A[I],A[J]);swap(B[I],B[J]);}K=N>>1;while(K>=2&&J>=K){J-=K;K>>=1;}J+=K;}for(L=1;L<=M;L++){LE=1<

4、;P+=LE){Q=P+LE1;Tr=Wr*A[Q]-Wi*B[Q];Ti=Wr*B[Q]+Wi*A[Q];A[Q]=A[P]-Tr;B[Q]=B[P]-Ti;A[P]+=Tr;B[P]+=Ti;}WTr=Wr;WTi=Wi;Wr=WTr*W1r-WTi*W1i;Wi=WTr*W1i+WTi*W1r;}}return;}#include"math.h"#include"stdio.h"#defineswap(a,b);T=(a);(a)=(b);(b)=TvoidIfft(floatA[],floatB[],unsignedM);main()

5、{floatA[1024],B[1024];unsignedlongN,I;unsignedM;printf("inputM=");scanf("%d",&M);N=1<

6、for(I=0;II){floatT;swap(A[I],A[J]);swap(B[I],B[J]);}K=N>>1;while(K>=2&&J

7、>=K){J-=K;K>>=1;}J+=K;}for(L=1;L<=M;L++){LE=1<

8、r*W1r-WTi*W1i;Wi=WTr*W1i+WTi*W1r;}}return;}#include#include