哈工程数字信号处理实验报告 谱分析.docVIP

数字信号处理实验 实验五：谱分析 班级： 姓名： 学号： 指导老师： 2010年10月 实验五 谱分析 1.实验原理 2.实验内容 1. w1=boxcar(20) subplot(221),stem(w1); title(boxcar) xlabel(t),ylabel(w1(t)); w2=hanning(20) subplot(222),stem(w2) title(hanning) xlabel(t),ylabel(w2(t)); w3=hamming(20) subplot(223),stem(w3) title(hamming) xlabel(t),ylabel(w3(t)) w4=bartlett(20) subplot(224),stem(w4) title(bartlett) xlabel(t),ylabel(w4(t)) w5=blackman(20) subplot(221),stem(w5); title(blackman) xlabel(t),ylabel(w5(t)); w6=triang(20) subplot(222),stem(w6) title(triang) xlabel(t),ylabel(w6(t)); w7=kaiser(20,80) subplot(223),stem(w7) title(kaiser) xlabel(t),ylabel(w7(t)) w8=chebwin(20,80) subplot(224),stem(w8) title(chebwin) xlabel(t),ylabel(w8(t)) 2. N=20 w1=boxcar(N) [X,W]=dtft(w1,500) subplot(221),plot(W/2/pi,abs(X)); title(boxcar) xlabel(w),ylabel(|W(jw)|); w2=hanning(N) [X,W]=dtft(w2,500) subplot(222),plot(W/2/pi,abs(X)); title(hanning) xlabel(w),ylabel(|W(jw)|); w3=hamming(N) [X,W]=dtft(w3,500) subplot(223),plot(W/2/pi,abs(X)); title(hamming) xlabel(w),ylabel(|W(jw)|); w4=bartlett(N) [X,W]=dtft(w4,500) subplot(224),plot(W/2/pi,abs(X)); title(bartlett) xlabel(w),ylabel(|W(jw)|); w5=blackman(N) [X,W]=dtft(w5,500) subplot(221),plot(W/2/pi,abs(X)); title(blackman) xlabel(w),ylabel(|W(jw)|); w6=triang(N) [X,W]=dtft(w6,500) subplot(222),plot(W/2/pi,abs(X)); title(triang) xlabel(w),ylabel(|W(jw)|); w7=kaiser(N,80) [X,W]=dtft(w7,500) subplot(223),plot(W/2/pi,abs(X)); title(kaiser) xlabel(w),ylabel(|W(jw)|); w8=chebwin(N,80) [X,W]=dtft(w8,500) subplot(224),plot(W/2/pi,abs(X)); title(chebwin) xlabel(w),ylabel(|W(jw)|); 3. w1=boxcar(10) [X,W]=dtft(w1,500) subplot(221),plot(W/2/pi,abs(X)); title(boxcar N=10) xlabel(w),ylabel(|W(jw)|); w2=boxcar(20) [X,W]=dtft(w2,500) subplot(222),plot(W/2/pi,abs(X)); title(boxcar N=20) xlabel(w),ylabel(|W(jw)|); w3=boxcar(50) [X,W]=dtft(w3,500) subplot(223),plot(W/2/pi,abs(X)); title(boxcar N=50) xlabel(w),ylabel(|W(jw)|); w4=boxcar(100) [X,W]=dtft(w