实验6有限冲激响应数字滤波器设计.高毅鹏.doc

实验有限冲激响应数字滤波器设计 实验目的： 二、实验原理： ,则其对应的单位抽样响应为 窗函数设计法的基本原理是用有限长单位抽样响应序列h n 逼近hd n ，由于hd n 往往是无限长序列，且是非因果的，所以用窗函数w n 将hd n 截断，并进行加权处理，得到: h n 就作为实际设计的FIR数字滤波器的单位抽样响应序列，其频率响应函数为： 式中，N为所选窗函数的长度 在MATLAB中，可以用b fir1 ,Wc,’ftype’,taper 等函数辅助设计FIR数字滤波器。归一化频率 ，当设计带通和带阻滤波器时，Wc为双元素相量；ftype代表滤波器类型，如’high’高通，’stop’带阻等；taper为窗函数类型，默认为海明窗；窗函数用blackman, hamming, hanning ， kaiser产生。 三、实验内容 1.设计一FIR低通滤波器，通带截止频率 ，通带衰减不大于0.25dB，阻带截止频率 ，阻带衰减不小于50dB。选择一个合适的窗函数，绘出所设计的滤波器的频率响应图。 wp 0.2*pi;ws 0.3*pi; tr_width ws-wp; M ceil 6.6*pi/tr_width ; n [0:1:M-1]; wc ws+wp /2; hd ideal_lp wc,M ; w_ham hamming M ; h hd.*w_ham; [H,w] freqz h,1,1000,whole ; H H 1:1:501 ; w w 1:1:501 ; mag abs H ; db 20*log10 mag+eps /max mag ; delta_w 2*pi/1000; Rp - min db 1:1:wp/delta_w+1 ; disp [实际通带波动为,num2str Rp ] ; As -round max db ws/delta_w+1:1:501 ; disp [最小阻带衰减为,num2str As ] ; subplot 2,2,1 ;stem n,hd ;title 理想抽样响应 ; axis [0 M-1 -0.1 0.3] ;xlabel n ;ylabel hd n ; subplot 2,2,2 ;stem n,w_ham ;title 海明窗 ; axis [0 M-1 0 1.1] ;xlabel n ;ylabel w n ; subplot 2,2,3 ;stem n,h ;title 实际脉冲响应 ; axis [0 M-1 -0.1 0.3] ;xlabel n ;ylabel h n ; subplot 2,2,4 ;plot w/pi,db ;title 衰减幅度 ; axis [0 1 -100 10] ;xlabel 以pi为单位的频率 ;ylabel 分贝数 ;grid; %子程序ideal_lp function hd ideal_lp wc,M %理想低通滤波器；hd为0到M-1之间的理想脉冲响应；wc为截止频率；M为理想滤波器的长度 alpha M-1 /2; n [0:1:M-1]; m n-alpha+eps; %加一个小数以避免零做除数 hd sin wc*m ./ pi*m ; FIR滤波器的阶次由过渡带宽度决定，窗函数的选择取决于阻带衰减，海明窗和布莱克曼窗均可提供大于50 dB的衰减，由于海明窗提供了较窄的过渡带，因此具有较小的阶数。在设计中没有使用通带波动值Rp 0.25dB，但必须检查设计的实际波动，验证它是否在给定容限内。 2.（选做）?设计一个数字带通滤波器，指标要求如下：通带边缘频率：，，通带衰减。阻带边缘频率：，，最小阻带衰减 。 要求：给出FIR数字滤波器的冲激响应，绘出幅度和相位频响曲线，讨论实现形式和特点。