MATLAB实验iir,用MATLAB设计IIR数字滤波器 实验七.doc

用MATLAB设计IIR数字滤波器 实验七

实验七：用MATLAB设计IIR数字滤波器

1、用双线性变换法设计切比雪夫Ⅱ型数字滤波器，列出传递函数并描绘模拟和数字滤波器的幅频响应曲线。

① 设计一个数字低通滤波器，要求：ωp=0.2П，Rp=1dB；阻带：ωs=0.35П，As=15dB，

滤波器采样频率Fs=10Hz。

程序清单如下：

wp=0.2*pi; %滤波器的通带截止频率

ws=0.35*pi; %滤波器的阻带截止频率

Rp=1;As=15; %滤波器的通阻带衰减指标

ripple=10^(-Rp/20); %滤波器的通带衰减对应的幅度值

Attn=10^(-As/20); %滤波器的阻带衰减对应的幅度值

%转换为模拟滤波器的技术指标

Fs=10;T=1/Fs;

Omgp=(2/T)*tan(wp/2);%原型通带频率的预修正

Omgs=(2/T)*tan(ws/2);%原型阻带频率的预修正

%模拟原型滤波器计算

[n,Omgc]=cheb2ord(Omgp,Omgs,Rp,As,'s')

[z0,p0,k0]=cheb2ap(n,As); %设计归一化的cheb2型模拟滤波器原型

ba1=k0*real(poly(z0)); %求原型滤波器的系数b

aa1=real(poly(p0));

[bb,aa]=lp2lp(ba1,aa1,Omgc); %变换为模拟低通滤波器

%用双线性变换法计算数字滤波器系数

[bd,ad]=bilinear(bb,aa,Fs)

%求数字系统的频率特性

[H,w]=freqz(bd,ad);

dbH=20*log10((abs(H)+eps)/max(abs(H)));

subplot(2,2,1);plot(w/pi,abs(H));

ylabel('|H|');title('幅度响应');axis([0,1,0,1.1]);

set(gca,'XTickMode','manual','XTick',[0,0.2,0.35,1]);

set(gca,'YTickMode','manual','YTick',[0,Attn,ripple,1]);grid

subplot(2,2,2);plot(w/pi,angle(H)/pi);

ylabel('phi');title('相位响应');axis([0,1,-1,1]);

set(gca,'XTickMode','manual','XTick',[0,0.2,0.35,1]);

set(gca,'YTickMode','manual','YTick',[-1,0,1]);grid

subplot(2,2,3);plot(w/pi,dbH);title('幅度响应(dB)');

ylabel('dB');xlabel('频率(pi)');axis([0,1,-40,5]);

set(gca,'XTickMode','manual','XTick',[0,0.2,0.35,1]);

set(gca,'YTickMode','manual','YTick',[-50,-15,-1,0]);grid

subplot(2,2,4);zplane(bd,ad);

axis([-1.1,1.1,-1.1,1.1]);title('零极点图');

n =

3

Omgc =

10.2282

bd =

0.1624 0.0056 0.0056 0.1624

ad =

1.0000 -1.4073 0.9056 -0.1623

由频率特性曲线可知，该设计结果再通阻带截止频率处能满足Rp≤1dB、As≥20dB的设计指标要求，系统的极点全部在单位圆内，是一个稳定系统。由n=3可知，由3阶的模拟低通滤波器原型用双线性变换法设计出来的切比雪夫Ⅱ型数字低通滤波器是一个3阶系统。其传递函数为：

② 设计一个数字高通滤波器，要求：ωp=0.35П，Rp=1dB；阻带：ωs=0.2П，As=15dB，

滤波器采样频率Fs=10Hz。

程序清单如下：

wp=0.35*pi; %滤波器的通带截止频率

ws=0.2*pi; %滤波器的阻带截止频率

Rp=1;As=15; %滤波器的通阻带衰减指标

ripple=10^(-Rp/20); %滤波器的通带衰减对应的幅度值

Attn=10^(-As/20); %滤波器的阻带衰减对应的幅度值

%转换为模拟滤波器的技术指标

Fs=10;T=1/Fs;

Omgp=(2/T)*tan(wp/2);%原型通带频率的预修正

Omgs=(2/T)*tan(ws/2);%原型阻带频率的预修正

%模拟原型滤波器计算

[n,Omgc]=ch