FFT频谱分析实验报告.docVIP

实验二： 用FFT作谱分析 实验目的 (1) 进一步加深DFT算法原理和基本性质的理解(因为FFT只是DFT的一种快速算法， 所以FFT的运算结果必然满足DFT的基本性质)。 (2) 熟悉FFT算法原理和FFT子程序的应用。 (3) 学习用FFT对连续信号和时域离散信号进行谱分析的方法， 了解可能出现的分析误差及其原因， 以便在实际中正确应用FFT。 二、实验原理 实验内容 (1) 对 2 中所给出的信号逐个进行谱分析。 ? 解：（1） n=0:3; xn1=[1 1 1 1]; XK18=fft(xn1,8); XK116=fft(xn1,16); n1=0:7; n2=0:15; subplot(131); stem(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=[1 2 3 4 4 3 2 1]; XK28=fft(xn2,8); XK216=fft(xn2,16); subplot(131); stem(n1,xn2); xlabel(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=[4 3 2 1 1 2 3 4]; XK38=fft(xn3,8); XK316=fft(xn3,16); subplot(131); stem(n1,xn3); xlabel(n1); 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,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); xn52=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); 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(