哈工程数字信号处理实验报告5.doc

数字信号处理实验 实验五：谱分析班级： 姓名： 学号： 指导老师： 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=20w1=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(w4,500)subplot(224),plot(W/2/pi,abs(X);title(boxcar N=100)xlabel(w),ylabel(|W(jw)|);4. f=25/16;dt=0.01;N=230;n=0:(N-1);x=0.75+3.4*cos(2*pi*f*dt*n)+2.7*cos(4*pi*f*dt*n)+1.5*sin(3.5*pi*f*dt*n)+2.5*sin(7*pi*f*dt*n);w1=boxcar(N);y1=x.*w1;Y1,W=dtft(y1,1000);subplot(221),plot(W/2/pi,abs(Y1);grid,title(信号周期0.9的矩形窗（幅频）);xlabel(f);ylabel(|Y1|);subplot(222),plot(W/2/pi,angle(Y1);grid,title(信号周期0.9的矩形窗（相频）);xlabel(f);ylabel(Y1);w2=hanning(N);y2=x.*w2;Y2,W=dtft(y2,1000);subplot(223),plot(W/2/pi,abs(Y2);grid,title(信号周期0.9的汉宁窗（幅频）);xlabel(f);ylabel(|Y2|);subplot(224),plot(W/2/pi,angle(Y2);grid,title(信号周期0.9的汉宁窗（相频）);xlabel(f);ylabel(Y2);f=25/16;dt=0.01;N=230;n=0:(N-1);x=0.75+3.4*cos(2*pi*f*dt*n)+2.7*cos(4*pi*f*dt*n)+1.5*sin(3.5*pi*f*dt*n)+2.5*sin(7*pi*f*dt*n);w1=hamming(N);y1=x.*w1;Y1,W=dtft(y1,1000);subplot(221),plot(W/2/pi,abs(Y1);grid,title(信号周期0.9的汉明窗（幅频）);xlabel(f);ylabel(|Y1|);subplot(222),plot(W/2/pi,angle(Y1);grid,title(信号周期0.9的汉明窗（相频）);xlabel(f);ylabel(Y1);w2=bartlett(N);y2=x.*w2;Y2,W=dtft(y2,1000);subplot(223),plot(W/2/pi,abs(Y2);grid,title(信号周期0.9的巴特利特窗（幅频）);xlabel(f);ylabel(|Y2|);subplot(224),plot(W/2/pi,angle(Y2);grid,title(信号周期0.9的巴特利特窗（相频）);xlabel(f);ylabel(Y2);f=25/16;dt=0.01;N=230;n=0:(N-1);x=0.75+3.4*cos(2*pi*f*dt*n)+2.7*cos(4*pi*f*dt*n)+1.5*sin(3.5*pi*f*dt*n)+2.5*sin(7*pi*f*dt*n);w1=blackman(N);y1=x.*w1;Y1,W=dtft(y1,1000);subplot(221),plot(W/2/pi,abs(Y1);grid,title(信号周期0.9的布莱克曼窗（幅频）);xlabel(f);ylabel(|Y1|);subplot(222),plot(W/2/pi,angle(Y1);grid,title(信号周期0.9的布莱克曼窗（相频）);xlabel(f);ylabel(Y1);w2=triang(N);y2=x.*w2;Y2,W=dtft(y2,1000);subplot(223),plot(W/2/pi,abs(Y2);grid,title(信号周期0.9的triang窗（幅频）);xlabel(f);ylabel(|Y2|);subplot(224),plot(W/2/pi,angle(Y2);grid,title(信号周期0.9的triang窗（相频）);xlabel(f);ylabel(Y2);f=25/16;dt=0.01;N=230;n=0:(N-1);x=0.75+3.4*cos(2*pi*f*dt*n)+2.7*cos(4*pi*f*dt*n)+1.5*sin(3.5*pi*f*dt*n)+2.5*sin(7*pi*f*dt*n);w1=kaiser(N,2);y1=x.*w1;Y1,W=dtft(y1,1000);subplot(221),plot(W/2/pi,abs(Y1);grid,title(信号周期0.9的kaiser窗（幅频）);xlabel(f);ylabel(|Y1|);subplot(222),plot(W/2/pi,angle(Y1);grid,title(信号周期0.9的kaiser窗（相频）);xlabel(f);ylabel(Y1);w2=chebwin(N,10);y2=x.*w2;Y2,W=dtft(y2,1000);subplot(223),plot(W/2/pi,abs(Y2);grid,title(信号周期0.9的切比雪夫窗（幅频）);xlabel(f);ylabel(|Y2|);subplot(224),plot(W/2/pi,angle(Y2);grid,title(信号周期0.9的切比雪夫窗（相频）);xlabel(f);ylabel(Y2);f=25/16;dt=0.01;N=282;n=0:(N-1);x=0.75+3.4*cos(2*pi*f*dt*n)+2.7*cos(4*pi*f*dt*n)+1.5*sin(3.5*pi*f*dt*n)+2.5*sin(7*pi*f*dt*n);w1=boxcar(N);y1=x.*w1;Y1,W=dtft(y1,1000);subplot(221),plot(W/2/pi,abs(Y1);grid,title(信号周期0.9的矩形窗（幅频）);xlabel(f);ylabel(|Y1|);subplot(222),plot(W/2/pi,angle(Y1);grid,title(信号周期0.9的矩形窗（相频）);xlabel(f);ylabel(Y1);w2=hanning(N);y2=x.*w2;Y2,W=dtft(y2,1000);subplot(223),plot(W/2/pi,abs(Y2);grid,title(信号周期0.9的汉宁窗（幅频）);xlabel(f);ylabel(|Y2|);subplot(224),plot(W/2/pi,angle(Y2);grid,title(信号周期0.9的汉宁窗（相频）);xlabel(f);ylabel(Y2);f=25/16;dt=0.01;N=282;n=0:(N-1);x=0.75+3.4*cos(2*pi*f*dt*n)+2.7*cos(4*pi*f*dt*n)+1.5*sin(3.5*pi*f*dt*n)+2.5*sin(7*pi*f*dt*n);w1=hamming(N);y1=x.*w1;Y1,W=dtft(y1,1000);subplot(221),plot(W/2/pi,abs(Y1);grid,title(信号周期0.9的汉明窗（幅频）);xlabel(f);ylabel(|Y1|);subplot(222),plot(W/2/pi,angle(Y1);grid,title(信号周期0.9的汉明窗（相频）);xlabel(f);ylabel(Y1);w2=bartlett(N);y2=x.*w2;Y2,W=dtft(y2,1000);subplot(223),plot(W/2/pi,abs(Y2);grid,title(信号周期0.9的巴特利特窗（幅频）);xlabel(f);ylabel(|Y2|);subplot(224),plot(W/2/pi,angle(Y2);grid,title(信号周期0.9的巴特利特窗（相频）);xlabel(f);ylabel(Y2);f=25/16;dt=0.01;N=282;n=0:(N-1);x=0.75+3.4*cos(2*pi*f*dt*n)+2.7*cos(4*pi*f*dt*n)+1.5*sin(3.5*pi*f*dt*n)+2.5*sin(7*pi*f*dt*n);w1=blackman(N);y1=x.*w1;Y1,W=dtft(y1,1000);subplot(221),plot(W/2/pi,abs(Y1);grid,title(信号周期0.9的布莱克曼窗（幅频）);xlabel(f);ylabel(|Y1|);subplot(222),plot(W/2/pi,angle(Y1);grid,title(信号周期0.9的布莱克曼窗（相频）);xlabel(f);ylabel(Y1);w2=triang(N);y2=x.*w2;Y2,W=dtft(y2,1000);subplot(223),plot(W/2/pi,abs(Y2);grid,title(信号周期0.9的triang窗（幅频）);xlabel(f);ylabel(|Y2|);subplot(224),plot(W/2/pi,angle(Y2);grid,title(信号周期0.9的triang