DSP实现FFT的代码.doc

1、DSP 实现的方法/#include “DSP281x_Device.h“ / DSP281x Headerfile Include File/#include “DSP281x_Examples.h“ / DSP281x Examples Include File#include #include #include #define pi atan(1.0)*4 /计算 PI 的值typedef struct /定义表示复数的结构体float real;float imag;FUSHU;FUSHU x1024;int w1024;FUSHU FUCHENG(FUSHU b1,FUSHU b2)

2、 /复数相乘函数FUSHU b3;b3.real=b1.real*b2.real-b1.imag*b2.imag;b3.imag=b1.real*b2.imag+b1.imag*b2.real;return(b3);void FFT(FUSHU *xin,int N) /FFT 运算int m,LH,I,k,J,M,K;float p,ps ;int B,N1;FUSHU w,T;M=10;/对应 N 值为/下面是倒序的程序 LH=N/2;J=LH;N1=N-2; /*变址运算 */for(I=1;I=K)J=J-K;K=K/2;J=J+K;/下面是 DIT-FFT 运算程序for(m=1;m=K)J=J-K;K=K/2;J=J+K;for(m=1;m0.001)break;if(n1024)puts(“FFT 不正确n“);if(n=1024)puts(“FFT 正确n“);return 0;