资源描述：
《FFT频谱分析实验报告.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、实验二：用FFT作谱分析一、实验目的(1)进一步加深DFT算法原理和基本性质的理解(因为FFT只是DFT的一种快速算法，所以FFT的运算结果必然满足DFT的基本性质)。(2)熟悉FFT算法原理和FFT子程序的应用。(3)学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确应用FFT。二、实验原理三、实验内容(1)对2中所给出的信号逐个进行谱分析。 解：（1）n=0:3;xn1=[1111];XK18=fft(xn1,8);XK116=fft(xn1,16);n1=0:7;n2=0:15;subplot(131);s
2、tem(n,xn1);xlabel('n');ylabel('xn1');subplot(132);stem(n1,abs(XK18));xlabel('n1');ylabel('XK18');title('xn的8点');subplot(133);stem(n2,abs(XK116));xlabel('n2');ylabel('XK116');title('xn的16点');（2）n1=0:7;n2=0:15;xn2=[12344321];XK28=fft(xn2,8);XK216=fft(xn2,16);subplot(131);stem(n1,xn2);xl
3、abel('n1');ylabel('xn2');subplot(132);stem(n1,abs(XK28));xlabel('n1');ylabel('XK28');title('xn2的8点');subplot(133);stem(n2,abs(XK216));xlabel('n2');ylabel('XK216');title('xn2的16点');（3）n1=0:7;n2=0:15;xn3=[43211234];XK38=fft(xn3,8);XK316=fft(xn3,16);subplot(131);stem(n1,xn3);xlabel('n1')
4、;ylabel('xn3');subplot(132);stem(n1,abs(XK38));xlabel('n1');ylabel('XK38');title('xn3的8点');subplot(133);stem(n2,abs(XK316));xlabel('n2');ylabel('XK316');title('xn3的16点');（4）n1=0:7;n2=0:15;xn41=cos((pi/4)*n1);xn42=cos((pi/4)*n2);XK48=fft(xn41,8);XK416=fft(xn42,16);subplot(141);stem(n1,
5、xn41);xlabel('n1');ylabel('xn41');subplot(142);stem(n2,xn42);xlabel('n2');ylabel('xn42');subplot(143);stem(n1,abs(XK48));xlabel('n1');ylabel('XK48');title('xn4的8点');subplot(144);stem(n2,abs(XK416));xlabel('n2');ylabel('XK416');title('xn4的16点');（5）n1=0:7;n2=0:15;xn51=sin((pi/8)*n1);xn5
6、2=sin((pi/8)*n2);XK58=fft(xn51,8);XK516=fft(xn52,16);subplot(141);stem(n1,xn51);xlabel('n1');ylabel('xn51');subplot(142);stem(n2,xn52);xlabel('n2');ylabel('xn52');subplot(143);stem(n1,abs(XK58));xlabel('n1');ylabel('XK58');title('xn5的8点');subplot(144);stem(n2,abs(XK516));xlabel('n2');
7、ylabel('XK516');title('xn5的16点');（6）n1=0:7;n2=0:15;xn61=cos(8*pi*n1)+cos(16*pi*n1)+cos(20*pi*n1);xn62=cos(8*pi*n2)+cos(16*pi*n2)+cos(20*pi*n2);XK68=fft(xn61,8);XK616=fft(xn62,16);subplot(141);stem(n1,xn61);xlabel('n1');ylabel('xn61');subplot(142);stem(n2,xn62);xlabel('n2');ylabel('xn6
8、2');s
