资源描述：
《fft频谱分析实验报告》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、实验二：用FFT作谱分析一、实验目的(1)进一步加深DFT算法原理和基本性质的理解(为FFT只是DFT的一种快速算法，所以FFT的运算结果必然满足DFT的基本性质)。(2)熟悉FFT算法原理和FFT子程序的应用。(3)学习用FFT对连续信号和时域离散信号进行谱分以便在实析的方法，了解可能出现的分析误差及其原际中正确应用FFTo二、实验原理调用绘图子程序(函绘制
2、雕)
3、曲线实验内容(1)对2中所给出的信号逐个进行谱分析。解：(1)n=0:3;xnl=[l111];XK18=fft(xnl，8);XK116=fft(xnl，16
4、);nl=0:7;n2=0:15;subplot(131);stem(n，xnl);xlabel(.n.);ylabel(’xnr);subplot(132);stem(nl,abs(XK18));xlabel(’nl’)；ylabelCXK18*);titleCxn的8点•);subplot(133);stem(n2,abs(XK116));xlabel(’n2’)；ylabel(,XK116,);title(’xn的16点’)；0.90.80.706c0.5X0.40.30.20-13532521.510.5xn的8点0
5、xn的16点n2(2)nl=0:7;n2=0:15;xn2=[l2344321];XK28=fft(xn2,8);XK216=fft(xn2,16);subplot(131);stem(nl，xn2);xlabel(’nl’)；ylabel(’xn2');subplot(132);stem(nl，abs(XK28));xlabelCnV);ylabel(’XK28’)；title(’xn2的8点’)；subplot(133);stem(n2,abs(XK216));xlabel(’n2’)；ylabel(’XK216’)；t
6、itle(’xn2的16点’)；xn2的8点2045252.1sex02468n120xn2的16点051015n28161421-1089L2X4202468n1(3)nl=0:7;n2=0:15;xn3=[43211234];XK38=ttt(xn3,8);XK316=fft(xn3,16);subplot(131);stem(nl，xn3);xlabel(’nl’)；ylabel(’xn3’);$ubplot(132);stem(nl,abs(XK38));xlabel(’nl’)；ylabel(.XK38’);tit
7、le(’xn3的8点’)；subplot(133);stem(n2,abs(XK316));xlabel(’n2’)；ylabel(’XK316’)；title(’xn3的16点•)；3.532.5空2X1.510.5020©18-16-14-12-00210-X8-6-4-20■8oxn3的16点Q(4)nl=0:7;n2=0:15;xn41=cos((pi/4)*nl);xn42=cos((pi/4)*n2);XK48=fft(xn41,8);XK416=fft(xn42,16);subplot(141);stem(nl
8、,xn41);xlabel(’nl’)；ylabel(’xn4r);subplot(142);stem(n2，xn42);xlabel(’n2’);ylabel(’xn42’)；subplot(143);stem(nl,abs(XK48));xlabel(.nl');ylabel(.XK48’);titleCxn4的8点，);subplot(144);stem(n2,abs(XK416));xlabel(’n2’)；ylabel(’XK416’)；titleCxn4的16点’)；X08060.4020-02•04-064)8
9、)2xn4的16点08060.40.2CM?0M•02-04-06•0.8510n2xn4的8点4r<3525^52x1.505IO150510n2(5)nl=0:7;n2=0:15;xn51=sin((pi/8)*nl);xn52=sin((pi/8)*n2);XK58=fft(xn51,8);XK516=fft(xn52,16);subplot(141);stem(nl,xn51);xlabel('nl');ylabel(’xn5r);subplot(142);stem(n2,xn52);xlabel(’n2’);yla
10、bel(’xn52’)；subplot(143);stem(nl,abs(XK58));xlabel(nl);ylabelCXK58’)；titlefxn5的8点’)；subplot(144);stem(n2,abs(XK516));xlabel(’n2’)；ylabel(.XK516.)
