精品文档-数字信号处理实验(MATLAB版)刘舒帆-第21章.pptVIP

由频率特性曲线可知，该设计结果在通阻带截止频率处能满足Rp≤1dB.As≥20dB的设计指标要求。由n＝3可知，由3阶的模拟低通原型用双线性变换法设计出来的切比雪夫Ⅰ型数字带通滤波器是一个6阶的系统，极点全部在Z平面的单位圆内，是一个稳定的系统。这个滤波器的传递函数应为注意:在使用［z0，p0，k0］＝cheb1ap(n，Rp)设计模拟低通原型时，需要输入通带衰减Rp，即切比雪夫Ⅰ型模拟低通原型是以通带衰减Rp为主要设计指标的。因此，由模拟低通原型变为数字带通(或带阻)滤波器时，使用lp2bp(或lp2bs)语句要求输入模拟通带带宽W0和中心频率BW，应采用通带截止频率来计算，即

bw＝Omgp2－Omgp1；w0＝sqrt(Omgp1*Omgp2)；%[ZK(]模拟滤波器通带带宽和中心频率

如果将例21-3改为:采用双线性变换法设计一个切比雪夫Ⅱ型数字带通滤波器，其它条件不变，则需要修改下面几句程序:bw＝Omgs2－Omgs1；w0＝sqrt(Omgs1*Omgs2)；%[ZK(]模拟滤波器阻带带宽和中心频率

［n，Omgn］＝cheb2ord(Omgp，Omgs，Rp，As，￠s￠)%计算阶数n和截止频率

［z0，p0，k0］＝cheb2ap(n，As)；%设计归一化的模拟原型滤波器

采用阻带截止频率来计算W0和BW，是因为切比雪夫Ⅱ型模拟低通原型是以阻带衰减As为主要设计指标的。程序运行结果如下:

n＝3

Omgn＝6.640760.2343

bd＝0.2537－0.0000－0.4733－0.00000.4733－0.0000－0.2537

ad＝1.00000.00000.00080.00000.42060.0000－0.0343

频率特性及零极点图形如图21-4所示。图21-4用双线性变换法设计切比雪夫Ⅱ型带通数字滤波器由程序数据和曲线可知，该设计结果在通阻带截止频率处能满足Rp≤1dB.As≥20dB的设计指标要求。由n＝3可知，由3阶的模拟低通原型用双线性变换法设计出来的切比雪夫Ⅱ型数字带通滤波器是一个6阶的系统，极点全部在z平面的单位圆内，是一个稳定的系统。这个滤波器的传递函数应为5.用双线性变换法设计IIR数字带阻滤波器

例21-4采用双线性变换法设计一个切比雪夫Ⅰ型数字带阻滤波器，要求:下通带wp1＝0.2p，上通带wp2＝0.8p，Rp＝1dB；阻带下限ws1＝0.3p，阻带上限ws2＝0.7p，As＝20dB，滤波器采样频率Fs＝1000Hz。

解由题目可知，本例只是将例21-3的条件改为相反，即将原带通滤波器通带的频率区域改为带阻滤波器阻带的频率区域，将原带通滤波器阻带的频率区域改为带阻滤波器通带的频率区域。程序只需作5句修改:ws1＝0.3*pi；ws2＝0.7*pi；%数字滤波器的阻带截 止频率

wp1＝0.2*pi；wp2＝0.8*pi；%数字滤波器的通带截止 频率

［ba，aa］＝lp2bs(ba1，aa1，w0，bw)；%变换为模拟 带阻滤波器

程序运行结果如下:

n＝3

*1.1市场与市场营销1.2我国汽车市场的发展与现状复习思考题实验21用双线性变换法设计IIR数字滤波器一、实验目的

(1)加深对双线性变换法设计IIR数字滤波器基本方法的了解。

(2)掌握用双线性变换法设计数字低通、高通、带通、带阻滤波器的方法。

(3)了解MATLAB有关双线性变换法的子函数。二、实验涉及的MATLAB子函数

bilinear

功能:双线性变换——将s域(模拟域)映射到z域(数字域)的标准方法，将模拟滤波器变换成离散等效滤波器。

调用格式:

［numd，dend］＝bilinear(num，den，Fs)；将模拟域传递函数变换为数字域传递函数，Fs为取样频率。

［numd，dend］＝bilinear(num，den，Fs，Fp)；将模拟域传递函数变换为数字域传递函数，Fs为取样频率，Fp为通带截止频率。［zd,pd,kd］＝bilinear(z,p,k,Fs)；将模拟域零极点增益系数变换到数字域,Fs为取样频率。

［zd,pd,kd］＝bilinear(z,p,k,Fs,Fp)；将模拟域零极点增益系数变换到数字域,Fs为取样频率,Fp为通带截止频率。

［Ad,Bd,Cd,Dd］＝bilinear(A,B,C,D,Fs)；将模拟域状态变量系数变换到数字域,F