- 1、本文档共32页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
实验十一 用双线性变换法设计 IIR数字滤波器;掌握用双线性变换法设计IIR DF的原理及具体设计方法,熟悉用双线性变换法设计IIR DF的计算机编程。
观察用双线性变换法设计的DF的频响特性,了解双线性变换法的特点。
熟悉用双线性变换法设计BW和CB型DF的全过程。;二、实验设备; 为了克服冲激响应不变法产生的频率混叠现象,这是从S平面到Z平面的标准变换z=esT的多值对应关系导致的,为了克服这一缺点,产生了双线性变换法。 ;双线性变换法的映射函数:;(1)把s=σ+jΩ带入上式得: ;由于幅度为1,把ω看做Ω的函数求解,得到: ;求最小阶数N,ωc
Buttord,cheblord
Cheb2ord,ellipord;利用模拟滤波器设计IIR数字滤波器的步骤;4) 设计归一化低通滤波器Ga(p) 。;巴特沃斯(Butterworth)滤波器、切比雪夫(Chebyshev)滤波器、椭圆(Cauer)滤波器、贝塞尔(Bessel)滤波器;;[n,wn]=buttord(wp,ws,Rp,Rs,’s’);;2)模拟低通原型;4)模拟→数字 ;[n,wc]=buttord(wp,ws,rp,rs);设计任意选频的巴特沃斯滤波器
对于低通滤波器,wpws;
对于高通滤波器,wpws;
对于带通滤波器,wp和ws分别为具有两个元素的矢量wp=[wp1, wp2]和ws=[ws1, ws2],并且ws1wp1 wp2ws2;
对于带阻滤波器,wp和ws分别为具有两个元素的矢量,wp=[wp1, wp2]和ws=[ws1, ws2],并且wp1ws1 ws2wp2;;[b,a]=butter(n,wc):设计2N阶带通数字滤波器
Wn是两个元素的向量Wn=[w1 w2],通频带为w1ww2
[b,a]=butter(n,wc, ‘high’):设计N阶高通数字滤波器
归一化截止频率为Wn, 以π为单位。
[b,a]=butter(n,wc, ‘low’): 设计N阶低通数字滤波器
归一化截止频率为Wn,以π为单位。
[b,a]=butter(n,wc, ‘stop’): 2N阶带阻滤波器,
Wn是两个元素的向量wc=[w1, w2],设计阻带为w1 ww2 ,以π为单位。;FS=1000;
Fl=200;Fh=300; %通带、阻带截止频率
Rp=1;Rs=25;
wp=Fl*2*pi/FS; ws=Fh*2*pi/FS; % 边界频率? 数字频率
OmegaP=2*FS*tan(wp/2); %频率预畸
OmegaS=2*FS*tan(ws/2);
%选择滤波器的最小阶数
[n,Wn]=buttord(OmegaP,OmegaS,Rp,Rs,s); %此处是代入经预畸变后获得的归一化模拟频率参数
[bt,at]=butter(n,Wn,s); % 设计n阶的巴特沃思模拟滤波器
[bz,az]=bilinear(bt,at,FS); %双线性变换为数字滤波器
[H,W] = freqz(bz,az); %求解数字滤波器的频率响应
plot( W*FS/(2*pi) , abs(H)); grid;
xlabel(频率/Hz); ylabel(幅值);;三、实验原理;例2 设计一个Butterworth型高通数字滤波器,在通带中允许的最大衰减为0.5dB,阻带内的最小衰减为40dB,通带上限临界频率为40Hz,阻带下限临界频率为30Hz。它的采样率为fs=100Hz;三、实验原理;chebyshevII型滤波器(阻带等波纹)设计
数字域:
[b,a]=cheby2(n,Rs,Wn) n阶chebyshevII DF
截止频率由Wn确定,阻带内的波纹由Rs确定
[b,a]=cheby2(n,Rs,Wn,’ftype’)
ftype=high: 设计截止频率为Wn的高通滤波器;
ftype=stop: 可设计出带阻滤波器
[z,p,k]=cheby2(n,Rs,Wn)
[zp,k]= cheby2(n,Rs,Wn,’ftype’)
[A,B,C,D]= cheby2 (n,Rs,Wn)
[A,B,C,D]= cheby2(n,Rs,Wn,’ftype’)
模拟域:
[b,a]=cheby2(n,Rs,Wn,’s’) n阶chebyshevII型 AF
截止频率为Wn
其余形式类似于数字域的;例3 设计一个CB-
文档评论(0)