FFT实验报告.doc

《FFT实验报告.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、.FFT的C语言编程一、上机目的掌握基2时域抽取FFT算法及其C语言编程。二、上机内容1．根据已给参考程序和说明，自行编写蝶形计算部分的程序；蝶形计算程序如下：temp[re]=w[re]*x[k+B][re]-w[im]*x[k+B][im];temp[im]=x[k+B][re]*w[im]+x[k+B][im]*w[re];x[k+B][re]=x[k][re]-temp[re];x[k+B][im]=x[k][im]-temp[im];x[k][re]=x[k][re]+temp[re];x[k][im]=x[k][im]+temp[im

2、];2．用FFT计算下面8点复数信号的离散傅利叶变换；x(0)=1+5i；x(1)=2+2i；x(2)=5+2i；x(3)=3+7i；x(4)=7+6i；x(5)=5+3i；x(6)=9+1i；x(7)=3+8i；显示结果如下：83152252377653913835.000000i34.000000-8.535534i5.121320-15.999999i7.000000-6.292893i-2.8786809.000001i-5.999998-1.464467i0.8786803.999999i9.000000-7.707107i-7.121

3、320截屏结果如下所示：范文..3．分别计算下面16、32点信号的FFT，fs=16Hz，利用Matlab的做图功能画出两个信号的频谱图(f-

4、X(ejw)

5、)，给出两个信号中所包含的频率值。16点信号：x(0)=0；x(1)=0.7071；x(2)=1；x(3)=0.7071；x(4)=0；x(5)=-0.7071；x(6)=-1；x(7)=-0.7071；x(8)=0；x(9)=0.7071；x(10)=1；x(11)=0.7071；x(12)=0；x(13)=-0.7071；x(14)=-1；x(15)=-0.7071；32点信号：x(0)

6、=0；x(1)=1.0898；x(2)=1.7071；x(3)=1.6310；x(4)=1；x(5)=0.2168；x(6)=-0.2929；x(7)=-0.3244；x(8)=0；x(9)=0.3244；x(10)=0.2929；x(11)=-0.2168；x(12)=-1；x(13)=-1.631；x(14)=-1.7071；x(15)=-1.0898；x(16)=0；x(17)=1.0898；x(18)=1.7071；x(19)=1.631；x(20)=1；x(21)=0.2168；x(22)=-0.2929；x(23)=-0.3244；x

7、(24)=0；x(25)=0.3244；x(26)=0.2929；x(27)=-0.2168；x(28)=-1；x(29)=-1.631；x(30)=-1.7071；x(31)=-1.0898；[00.707110.70710-0.7071-1-0.707100.707110.70710-0.7071-1-0.7071]16点FFT在WIN-TC中的运行结果如下：000.70710100.7071000-0.70710-10-0.70710000.70710100.7071000-0.70710-10范文..-0.707100.000000i0.

8、0000000.000000i0.000000-0.000000i-7.9999600.000000i0.0000000.000000i0.0000000.000000i0.000000-0.000001i0.0000390.000000i0.0000000.000000i0.0000000.000000i0.0000000.000000i-0.0000390.000000i0.0000000.000000i0.0000000.000000i0.0000000.000001i7.9999600.000000i0.00000016点Matlab运行

9、程序如下：x=[00.707110.70710-0.7071-1-0.707100.707110.70710-0.7071-1-0.7071]N=16;y=myfft(x,N);y=abs(y);n=0:N-1;stem(n,y)16点运行结果如下：范文..即y16=[0080000000000080]由图进行分析，已知Fs=16Hz，那么F=Fs/N=1Hz读图知在X(2)有峰值，对应频率为2*1=2Hz即该信号所包含的频率值为2Hz。32点FFT在WIN-TC中的运行结果如下：0.000001i0.0000000.000000i0.00000

10、0-0.000000i-16.0001600.000000i0.000000-0.000001i-15.9999200.000000i0