数字信号处理.doc

数字信号处理上机作业6.35 利用Matlab设计模拟带阻巴特沃斯滤波器，要求通带下、上截止频率分别为，阻带下、上截止频率分别为，通带最大衰减，阻带最小衰减。解：Matlab源程序及结果如下：Ws1=9;Wsu=11;Wp1=6;Wpu=13;Wp=Wp1,Wpu;Ws=Ws1,Wsu;Rp=1;Rs=20;N,Wc=buttord(Wp,Ws,Rp,Rs,s);B,A=butter(N,Wc,stop,s)W=0:0.001:18;H,W=freqs(B,A,W);H=20*log10(abs(H);plot(W,H),grid onxlabel(频率/(rad/s);ylabel(幅度/dB)系统函数分子、分母多项式系数：幅频响应特性曲线：6.42 利用Matlab工具箱函数，分别设计满足下列指标的数字低通巴特沃斯滤波器和数字低通椭圆滤波器，画出各自的幅频响应特性曲线，并进行比较。设计指标为;通带截止频率，阻带截止频率；通带最大衰减，阻带最小衰减。解：Matlab源程序及结果如下：(1)巴特沃斯数字低通滤波器Wp=0.2*pi;Ws=0.4*pi;Rp=0.5;Rs=30;N,Wc=buttord(Wp,Ws,Rp,Rs,s);B,A=butter(N,Wc)W=0:0.0001:6;H,W=freqs(B,A,W);H=20*log10(abs(H);plot(W,H),grid onxlabel(频率 /(rad/s);ylabel(幅度/dB)系统函数分子、分母多项式系数：幅频响应特性曲线：(2)数字低通椭圆滤波器Wp=0.2*pi;Ws=0.4*pi;Rp=0.5;Rs=30;N,Wpo=ellipord(Wp,Ws,Rp,Rs,s);B,A=ellip(N,Rp,Rs,Wpo);W=0:0.001:6;H,W=freqs(B,A,W);H=20*log10(abs(H);plot(W,H),grid onxlabel(频率/(rad/s);ylabel(幅度 /dB)系统函数分子、分母多项式系数：幅频响应特性曲线：7.22 利用Matlab工具箱函数firl，设计线性相位FIR数字低通滤波器，要求希望逼近的理想低通滤波器的通带截止频率，滤波器长度N=21。试分别选用矩形窗、汉宁窗、汉明窗和布莱克曼窗进行设计，画出每种窗函数所设计的滤波器的单位脉冲响应h(n)波形及幅频响应特性曲线，并进行比较，观察各种窗函数的设计性能。解：Matlab源程序及结果如下：N=21;wc=1/4;n=0:20;fprintf(矩形窗设计结果:)hn=fir1(N-1,wc,boxcar(N);fh=fft(hn,1024);fh=20*log10(abs(fh);wk=0:1023;wk=2*wk/1024;subplot 421,stem(n,hn);grid;subplot 422,plot(wk,fh);grid;title(矩形窗);xlabel(w/pi)矩形窗设计结果:wind=hanning(N);fprintf(汉宁窗设计结果:)hn=hn.*wind;fh=fft(hn,1024);fh=20*log10(abs(fh);subplot 423,stem(n,hn);grid;subplot 424,plot(wk,fh);grid;title(汉宁窗);xlabel(w/pi)axis(0 2 -130 0);汉宁窗设计结果:wind=hamming(N);fprintf(汉明窗设计结果:)hn=hn.*wind;fh=fft(hn,1024);fh=20*log10(abs(fh);subplot 425,stem(n,hn);grid;subplot 426,plot(wk,fh);grid;title(汉明窗);xlabel(w/pi)汉明窗设计结果:axis(0 2 -160 0);wind=blackman(N);fprintf(布莱克曼窗设计结果:)hn=hn.*wind;fh=fft(hn,1024);fh=20*log10(abs(fh);subplot 427,stem(n,hn);grid;subplot 428,plot(wk,fh);grid;title(布莱克曼窗);xlabel(w/pi)布莱克曼窗设计结果:将四种窗函