数字信号课程设计郑大.docVIP

PAGE  PAGE 23 实验一：用DFT作谱分析 一、实验目的 1、进一步加深对DFT的基本性质的理解。 2、熟悉FFT算法原理和FFT子程序的应用。 3、学习用DFT对连续信号和离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确应用DFT。 二、上级程序： x1=[1 1 1 1]; x2=[1 2 3 4 4 3 2 1]; x3=[4,3,2,1,1,2,3,4]; figure subplot(1,3,1);stem(x1); title(序列x1) xlabel(n); ylabel(x1); subplot(1,3,2);stem(x2); title(序列x2) xlabel(n); ylabel(x2); subplot(1,3,3);stem(x3); title(序列x3) xlabel(n); ylabel(x3); N1=8; F1x1=fft(x1,N1); F1x2=fft(x2,N1); F1x3=fft(x3,N1); figure subplot(1,3,1);stem(abs(F1x1)); title(序列x1的8点谱分析) xlabel(n); ylabel(F1x1); subplot(1,3,2);stem(abs(F1x2)); title(序列x2的8点谱分析) xlabel(n); ylabel(F1x2); subplot(1,3,3);stem(abs(F1x3)); title(序列x3的8点谱分析) xlabel(n); ylabel(F1x3); N2=16; F2x1=fft(x1,N2); F2x2=fft(x2,N2); F2x3=fft(x3,N2); figure subplot(1,3,1);stem(abs(F2x1)); title(序列x1的16点谱分析); xlabel(n); ylabel(F2x1); subplot(1,3,2);stem(abs(F2x2)); title(序列x2的16点谱分析); xlabel(n); ylabel(F2x2); subplot(1,3,3);stem(abs(F2x3)); title(序列x3的16点谱分析); xlabel(n); ylabel(F2x3); N3=16; n=0:N3-1; x4=cos(n*pi/4); x5=sin(n*pi/8); figure; subplot(1,2,1);stem(n,x4); title(序列x4) xlabel(n); ylabel(x4); subplot(1,2,2);stem(n,x5); title(序列x5) xlabel(n); ylabel(x5); N4=8; F3x4=fft(x4,N4); F3x5=fft(x5,N4); figure subplot(1,2,1);stem(abs(F3x4)); title(序列x4的8点谱分析); xlabel(n); ylabel(F3x4); subplot(1,2,2);stem(abs(F3x5)); title(序列x5的8点谱分析); xlabel(n); ylabel(F3x5); N5=16; F4x4=fft(x4,N5); F4x5=fft(x5,N5); figure subplot(1,2,1);stem(abs(F4x4)); title(序列x4的16点谱分析); xlabel(n); ylabel(F4x4); subplot(1,2,2);stem(abs(F4x5)); title(序列x5的16点谱分析); xlabel(n); ylabel(F4x5); N6=64; F5x4=fft(x4,N6); F5x5=fft(x5,N6); figure w=2/N6*[0:N6-1]; subplot(1,2,1);plot(w,abs(F5x4)); title(序列x4的64点谱分析); xlabel(w); ylabel(F5x4); subplot(1,2,2);plot(w,abs(F5x5)); title(序列x5的64点谱分析); xlabel(w); ylabel(F5x5); N7=256; F6x4=fft(x4,N7); F6x5=fft(x5,N7); figure w=2/N7*[0:N7-1]; subplot(1,2,1);plot(w,abs(F6x4)); title(序列x4的256点谱分析); xlabel(w); ylabel(F6x4); subplot(1,2,2);plot(w,abs(F6x5)); title(序列x5的256点谱分析); xlabel(w);