有限冲激响应数字滤波器设计(docX页)

1、有限冲激响应数字滤波器设计(doc X页)实验5有限冲激响应数字滤波器设计 一、实验目的：、实验原理:1、加深对数字滤波器的常用指标理解。 2、学习数字滤波器的设计方法。OfP issbawiEt 吁bwridTransi tiFig 1 Typical magnitude specification for a digita1 LPF低通滤波器的常用指标：闻上%nr<|Q|<7r四.(1)通带边缘频率；(2)阻带边缘频率;(3)通带起伏；3T91里J1-期(4)通带峰值起伏，/% =701电式为)的(5)阻带起伏，最小阻带衰减。三、预习要求1 .在 MATLAB,熟悉函数 fir

2、l、kaiserord 、remezord、remez的使用；B = fir1(n,Wn,'high','noscale')设计滤波器；n,Wn,beta,ftype = kaiserord(f,a,dev) 估计滤波器数；n,fo,ao,w = remezord (f,a,dev,fs)计算等波纹滤波器阶数n和加权函数w(w); B=remez(n,f,a) 进行等波纹滤波器的设计。2 .阅读扩展练习中的实例，学习FIR滤波器的设计方法及其在MATLA呻的实现；3 .给出FIR数字滤波器的冲激响应，绘出它们的幅度和相位频响曲线，讨论它们各自的实现形式和特点。数

3、字滤波器有IIR和FIR两种类型，它们的特点和设计方法不同。四、实验内容：利用MATLABS程，分别用窗函数法和等波纹滤波器法设计两种FIR数字滤波器，指标要求如下：=0.457 .= 0.65/r通带边缘频率：，通带峰值起伏：门了 V 1口6的品二0.3万二。少=刀。阻带边缘频率：，最小阻带衰减：。as>4d£(1)用窗函数法实现：调用函数n,wn,bta,ftype=kaiserord(f , a, dev, fs) 参数:f=0.3 0.450.65 0.8为对应数字频率Cpi =0.45才。§ =03=0.651=0.75乃a=0 1 0为由f指定的各个频带上

4、的幅值向量，一般只有0和1表示；和f长度关系为(2*a的长度)一2二(f的长度)devs=0.01 0.1087 0.01用于指定各个频带输出滤波器的频率响应与其期望幅值之间的最大输出误差或偏差，长度与 a相等，计算公式：阻带衰减误差=,通 带波动衰减误差=2=10毋加=访=1 10”侬fs缺省值为2HZN0=ceil(4.46*pi/Bt(2);B=fiU(N0-1,wc,'bandpass',kaiser(N0,3.384);figure(5)freqz(B,1)title('kaiser window') figure(6)stem(B)title(

5、9;kaiser window') (2)用等波纹法设计：igurS50 mH13-ril-EFF-r石口TQhum isFf iquemv 产，i W曹 arr- pir ：0.5 Q而 Q； 。将 Q.9BE* Edit Vievj Irispr* Tnsls D esIlzP Wiria, Eel，1Aq.*以岂f）看便口鬲白t50 a 1021.1尸尸J1y一4J- "1 i11111m.L:、 -卜 _、11I11500-1900ShflO2900- 匕* g*Q1 QZ”04050 &0T OS 031ifamna lzfJ Fr&wwncy ra

6、lXsarFplE；调用函数n,fpts,mag,wt=remezord(f ， a， dev) f=0.3 0.45 0.65 0.8a=0 1 0dev=0.01 0.1087 0.01 其含义同函数n,wn,bta,ftype=kaiserord(f ， a，dev， fs) 中的参数相同。f=0.3,0.45,0.65,0.75; m=0 1 0;rp=1;rs=40;dat1=(10A(rp/20)-1)/(10A(rp/20)+1);dat2=10A(-rs/20);rip=dat1,dat2,dat1; M,fo,mo,w=remezord(f,m,rip);M=M+1;hn=remez(M,fo,mo,w);figure(9)freqz(hn,1)title('remez')figure(10)stem(hn)title('remez')山 pdlr-in七 JEnlflFile EdftVie*-"不，?士 Ten tW r>dfWHelp-*口曲。身*r自耳国前即六、实验总结本实验主要是加深对数字滤波器的常用指标理解，学习数字滤波器的设计方法。在MATLAB用函数B =fir1(n,Wn,'high','noscale')来设计滤波器；n,Wn,beta,ftyp