滤波器实验5基于matlab的数字滤波器设计

滤波器实验5基于matlab的数字滤波器设计内容摘要：

1,2) plot(w1/pi,angle(hh1)) grid xlabel(39。 归一化频率 w39。 )。 ylabel(39。 相位 /rad39。 )。 h1 = Columns 1 through 8 Columns 9 through 16 Columns 17 through 24 Columns 25 through 31 图 41 如果直接用 freqz(h1,1,256)。 画图得： 12）用等波纹法设计： 调用函数 [n,fpts,mag,wt]=remezord(f， a， dev) f=[ ] a=[0 1 0] dev=[ ] 其含义同 函数 [n,wn,bta,ftype]=kaiserord（ f， a， dev， fs）中的参数相同。 [n,fpts,mag,wt]=remezord([ ],[0 1 0],[ ])。 %用 remezord 函数估算出 remez 函数要用到的阶 n、归一化频带边缘矢量 fpts、频带内幅值响应矢量 mag 及加权矢量 w，使 remez函数设计出的滤波器满足 f、 a 及 dev 指定的性能要求。 h2=remez(n,fpts,mag,wt)。 %设计出等波纹滤波器 [hh2,w2]=freqz(h2,1,256)。 figure(2) subplot(2,1,1) plot(w2/pi,20*log10(abs(hh2))) grid xlabel(39。 归一化频率 w39。 )。 ylabel(39。 幅度 /db39。 )。 subplot(2,1,2) plot(w2/pi,angle(hh2)) grid xlabel(39。 归一化频率 w39。 )。 ylabel(39。 相位 /rad39。 )。 h2 h2 = Columns 1 through 9 Columns 10 through 17 图 42 用 freqz(h2,1,256)。 直接得图： 二、 对课本作业。 21）用汉宁窗实现： 1）手动计算 n=:。 hd=sin(*pi*(n29))./(pi*(n29))。 win=+*cos(2*pi*(n29)/58)。 h1=2*cos(pi*(n29)/2).*hd.*win。 [hh1,w1]=freqz(h1,1,256)。 figure(1) subplot(2,1,1) plot(w1,20*log10(abs(hh1))) grid xlabel(39。 数字频率 w/rad39。 )。 ylabel(39。 幅度 /db39。 )。 subplot(2,1,2) plot(w1,angle(hh1)) grid xlabel(39。 数字频率 w/rad39。 )。 ylabel(39。 相位 /rad39。 )。 图 43 22）用自带的 fir1 函数： n=59。 wn=[3/8 5/8]。 h2=fir1(n,wn,39。 bandpass39。 ,hann(n+1))。 [hh2,w2]=freqz(h2,1,256)。 figure(2) subplot(2,1,1) plot(w2,20*log10(abs(hh2))) xlabel(39。 数字频率 w/rad39。 )。 ylabel(39。 幅度 /db39。 )。 grid su。