《数字信号处理》课件_第7章.ppt

4.归一化低通滤波器到数字带阻滤波器的变换方法

求解过程与上相同，在此直接给出结果。

(7.4.44)(7.4.43)

式中

求归一化低通滤波器指标时，有

λs=f(ωsl)，-λs=f(ωsu)7.4.4用Matlab设计IIR数字滤波器

Matlab信号处理工具箱提供了丰富的IIR数字滤波器设计函数，分别介绍如下。

IIR数字滤波器设计技术指标如下：

Wp：归一化通带截止频率ωp/π，当滤波器为带通或带阻时，Wp=［ωpl/π，ωpu/π］。

Ws：归一化阻带截止频率ωs/π，当滤波器为带通或带阻时，Ws=［ωsl/π，ωsu/π］。

Rp：通带最大衰减αp。

Rs：阻带最小衰减αs。1.根据滤波器技术指标计算数字滤波器阶数

根据实际数字滤波器技术指标计算数字滤波器阶数。

［N，Wc］=buttord(Wp，Ws，Rp，Rs);

求巴特沃思型滤波器的阶数N和3dB通带截止频率Wc。［N，Wpo］=cheb1ord(Wp，Ws，Rp，Rs)

;求切比雪夫Ⅰ型滤波器的阶数N和通带截止频率Wpo。

［N，Wso］=cheb2ord(Wp，Ws，Rp，Rs);求切比雪夫Ⅱ型滤波器的阶数N和阻带边界频率Wso。

［N，Wpo］=ellipord(Wp，Ws，Rp，Rs);求椭圆型滤波器的阶数N和通带截止频率Wpo。2.根据滤波器阶数和通(阻)带截止频率

求数字滤波器的系统函数

［b，a］=butter(N，Wc，′ftype′);求巴特沃思型数字滤波器系统函数分子、分母多项式的系数向量b和a。向量中的元素bi和ak为下式中分子和分母的系数。

Wc为单项，ftype缺省时，设计低通数字滤波器;ftype=high时，设计高通数字滤波器。

Wc=［Wcl，Wcu］，ftype缺省时，设计带通数字滤波器，通带为［Wcl，Wcu］;ftype=stop时，设计带阻数字滤波器，Wcl、Wcu分别是通带下、上截止频率。

［b，a］=cheby1(N，Rp，Wpo，′ftype′);求切比雪夫Ⅰ型数字滤波器系统函数分子、分母多项式的系数。［b，a］=cheby2(N，Rs，Wso，′ftype′);求切比雪夫Ⅱ型数字滤波器系统函数分子、分母多项式的系数。

［b，a］=ellip(N，Rp，Rs，Wpo，′ftype′);求椭圆型数字滤波器系统函数分子、分母多项式的系数。【例7.4.9】用Matlab设计IIR数字高通滤波器，采用切比雪夫I型低通滤波器原型。

设计指标如下:

解(1)按照首先设计模拟滤波器再用双线性变换法转换为数字滤波器的思路，

Matlab程序如下:

Fs=1000;fp=300;fs=200;

Wp1=1;Ws1=tan(pi*fp/Fs)/tan(fs*pi/Fs);

Rp=1;Rs=40;

%求模拟滤波器的技术指标［N，Wpo］=cheb1ord(Wp1，Ws1，Rp，Rs，′s′); %求归一化低通滤波器的阶数和通带截止频率

［z，p，k］=cheb1ap(N，Rp);

%求归一化低通滤波器的零、极点

［bl，al］=zp2tf(z，p，k);

%求归一化低通滤波器的系统函数

wp=2*tan(pi*fp/Fs)*Fs;

%求模拟高通滤波器的通带截止频率

［b，a］=lp2hp(bl，al，wp);

%求模拟高通滤波器的系统函数

［bz，az］=bilinear(b，a，Fs);

%利用双线性变换法求数字滤波器的系统函数

w=0：0.01*pi：pi;［h，w］=freqz(bz，az，w);

%求数字滤波器的幅频响应程序运行结果如下:

N=5，bz=［0.0079，-0.0397，0.0794，-0.0794，0.0397，-0.0079］;

az=［1.0000，2.2188，3.0019，2.4511，1.2330，0.3109］;以上是首先设计归一化低通滤波器，再通过模拟频率转换求出模拟高通滤波器，最后通过双线性变换库函数求出数字高通滤波器的系统函数。也可以将模拟滤波器指标预畸变

后，直接设计模拟高通滤波器，再用双线性变换