资源描述:
《数字信号处理实验5有限冲激响应数字滤波器设计》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库。
1、实验5有限冲激响应数字滤波器设计•、实验H的:、加深对数卞滤波器的常川指标现解.2、学习数字滤波器的i5计方法。1+(5
2、^(o)
3、1-Sp—OpOcPassbandStopbandTransitionbandFig1TypicalmagnitudespecificationforadigitalLPF低通滤波器的常用指标:Jl-(5p<
4、7/(Q)
5、6、Q
7、<
8、7/(Q)
9、10、Q
11、12、AB中,熟悉函数Hrl、kaiserord、reraezord^remez的使用:B=firlWn,•high•,’noscale’)设计滤波器;[n,Wn,beta,ftype]=kaiserord(f,a,dev>估计滤波器阶数:[n,fo,ao,w]=roinozord(f,a,dev,fs)计算等波纹诚波器阶数ri和加权函数*(<•>>:B=remez(n.f,a)进行等波纹滤波器的没计,2.阅读扩展练习中的实例,学习FIR滤波器的设计方法及其在MATLAB屮的劣现;3.给出FIK数7•滤波器的冲激响应.绘出它们的幅度和相位频响曲线.讨论它们各自的实现形式和特
13、点,数字滤波器有IIK和FIKM种类型,它们的特点和设计方法不同,四、实验内咨:利用MATLAB编程,分別川说函数法和等波纹滤波器法没H•两种FIR数字滤波器.拊标要求如下:QP1=0.45/T.Qp^=0.65/Tad14、0.1102x(50-8.7)=4.55,□.>□2DJ0A□DA□JDJ02归一化构笵/71□2□JDA□A□ADJ□J归一化相泡/71D□Q□□□□O(3PX吆法(a)(b)上阁给出了以上设计的猢率特性,(a)为N-30迕接截取的频率特性(b〉为凯浆说设计的频率特性,凯堪斑设计对应的MATLAB程序为:wn=kaiser(30,4.55):nn=[0:l:29];alfa=(30-l)/2;hd=sin(0.4*pi*{nn-alfa))./(pi*{nn-alfa)):h=hd.*wn*;[htwl]=freqz(hj);或者:b=firl(29,0.4,kai
15、ser(30,4.55)):[hbwl]=freqz(b>l);plot(wl/pi,20*logl0(abs(hl))):axis([0,1,-80,10]);grid;xlabel(’归一•化頻率/p’);ylabelf幅度/dB’>:还可以使用[n.Wn,beta,ftype]=kaiserord(f,a,dev)函数來估计滤波器阶数等,得到凯寒窗滤波器:fcuts=[0.30.5]:%归—•化顧率omega/pimags=[10];devs=[0.0510"<-2.5)];[n,Wn,beta,ftypo]=kaiserord(fcuts,mags,devs)
16、:ft计算出飢塞窗N,beta的值hh=firl(n,Wn,ftype,kaiser(n+Lbeta),*noscale*);freqz(hh):实呩中,一般词川MATLAB信号处理工具箱函数rernezord來it•算等波纹滤波器阶数N和加权函数W(o),调用函数remez可进行等波纹滤波器的设汁.斑接求出滤波器系数,函数rernezord屮的数组fedge为通带和附带边界频率,数纟flmval足两个边界处的幅偾.而数饥dev足通带和阻带的波动,fs是采样频卒申•位力liz例2:利川盂米兹交矜筇法没计等波纹滤波器,设汁一个线性枏位低通FIR数字滤波器,其指标力:通带
17、辺界頻率fc=800Hz,附带边界fr=10001lz.通带波动<5=0.56®P11带坡小袞减At_40(IB,采样猢率fsMOOOHz解:S1=l-10'J/2°=0.0559,32=10一20=0.01在MATLAB中可以用reroezord和rcmez两个函数设汁,K结果如阁2,MATLAB柺序如下:fedge=[8001000];mval=[l0];dev二[0.05590.01]:fs-4000;[N,fpts,mag,wt]=remezord(fedge,mvahdev,fs):b=remez(N,fpts,mag,wt):[h,w]=