IIR模拟低通滤波器设计.doc

信息科学与工程学院
综合性实验报告

姓 名： 学 号

班 级：

实验项目名称： IIR模拟低通滤波器设计

实验项目性质： 设计性实验

实验所属课程： 数字信号解决

实验室(中心)：

指 导 教 师 ：

实验完毕时间： 年 月 日
教师评阅意见：

署名： 年 月 日
实验成绩：

一、实验目的
1、加深对无限冲激响应IIR 滤波器的常用指标和设计过程的理解。
2、学会用冲激响应不变法把模拟滤波器设计成数字滤波器。
3、进一步掌握matlab在数字信号解决中的应用，以便以后的学习。
二、实验内容及规定
实验内容：规定按照设计指标设计无限冲激响应IIR巴特沃什模拟低通 滤波器。
实验规定：必须掌握IIR巴特沃什模拟低通滤波器的各个指标代表的含义，搞清楚次实验的原理，有也许的话，用冲激响应不变法把模拟滤波器设计成数字滤波器。
三、实验原理
1、低通滤波器的技术规定用图形表达如下：

通带
通带
通带

阻带
过渡带
通带

如上图表达了一个频域设计、一维低通滤波器的技术规定图。 和 分别称为通带截止频率和阻带截止频率。通带频率范围为 ，阻带频率范围为 。从 到 称为过渡带，用 表达，在过渡带里，幅频特性单调下降。在通带和阻带内的衰减一般用数 dB表达。
通带内允许最大衰减是 ，阻带内允许最小衰减是 ，定义分别为：

2、 Matlab信号解决工具箱中提供了设计巴特沃思模拟滤波器的函数buttord、 buttap和butter，格式如下：
l
用于计算巴特沃思模拟低通滤波器的阶N和3dB截止频率Wc (即本书中的符号)。其中，Wp 和Ws分别是滤波器的通带截止频率和阻止截止频率，单位为rad/s；Rp和Rs分别是通带最大衰减系数和阻带最小衰减系数，单位为dB。
l
用于计算N阶巴特沃思归一化（=1）模拟低通滤波器系统函数的零、极点和增益因子，返回长度为N的向量z和p分别给出N个零点和极点，G是滤波器增益。得到的滤波器系统函数形式如下：

其中，和分别是向量z和p的第k个元素。假如要从零、极点得到系统函数的分子和分母多项式系数向量B和A，可以调用结构转换函数。（实验中没有用到）
l ，结构转换后系统函数的形式为

其中，M 是向量B的长度，N是向量A的长度，分别是向量B和A的第k个元素。
（3）
用于计算巴特沃什模拟滤波器系统函数中分子和分母多项式系数向量B和A，其中N和分别是滤波器的阶和3dB截止频率，返回向量B和A中的元素和分别是上面的表达式中的分母和分子系数。ftype缺省时，设计低通滤波器，S缺省时，设计数字滤波器。
四、实验仪器、材料
PC机一台、MATLAB软件实验平台
五、实验过程及原始记录
1、巴特沃什模拟低通 滤波器设计指标如下：
wp=0.4*pi，ws=0.45*pi，Rp=5，Rs=10;
实验所用到的程序如下：
wp=0.4*pi;ws=0.45*pi;Rp=5;Rs=10;
[N,wc]=buttord(wp,ws,Rp,Rs,'s');
[B,A]=butter(N,wc,'s');
w=0:0.001:6;
[H,W]=freqs(B,A,w);
H=20*log10(abs(H));
plot(w,H),grid on;
xlabel(‘频率/(rad/s)');ylabel('幅度/dB')
经MATLAB仿真以后得到的滤波器如下图所示：

运营后在matlab控制窗口上显示如下结果：
N=7
2、为了方便观测故设计了不同阶次N的滤波器，即Butterworth模拟原型低通滤波器，其程序代码如下：
n=0:0.01:2;
for ii=1:4
switch ii
case 1,N=2;
case 2,N=5;
case 3,N=10;
case 4,N=20;
end
[z,p,k]=buttap(N);
[b,a]=zp2tf(z,p,k);
[H,w]=freqs(b,a,n);
magH2=(abs(H)).^2;
hold on;
plot(w,magH2);
end
xlabel('w/wc');
ylabel('|H(jw)|^2');
title('Butterworth模拟原型低通滤波器');
text(1.5,0.18,'n=2')
text(1.3,0.08,'n=5')
text(1.16,0.08,'n=10')
text(0.93,0.98,'n=20')
grid on;
仿真后结果图形如图：

3、 运用冲激响应不变法，把系统函数为的模拟滤波器变换成等价的数字滤波器，采样间隔T=1S。
（此处为了方便计算并没有用环节1中产生的模拟滤波器）
先通过计算得到：，则有如下式子

Matlab实验参考程序如下:
B=[1,1];
A=[1,5,6];
T=1;
Fs=1/T;
[Bz,Az]=impinvar(B,A,Fs); %用冲激响应不变法将模拟滤波器变换成数字滤波器
运营结果：
Bz =1.0000 -0.2209
Az = 1.0000 -0.1851 0.0067
得到的数据和先前计算得到的数据相同，故运用冲激响应不变化法把模拟滤波器设计成数字滤波器对的。
六、实验结果及分析
1、环节1 中产生的低通模拟滤波器的系统函数，经计算为：

通过仿真后的结果，与上式对比，所得结果进分析相同。
此处为N=7的低通模拟滤波器
2、在Butterworth模拟原型低通滤波器中，可以看到，滤波器的幅频特性随着滤波器阶次N的增长而变得越来越好。
3、环节3中，运用冲激响应不变法，把系统函数为的模拟滤波器变换成等价的数字滤波器，在matlab仿真软件中验证了冲激响应不变法的对的性。
七、实验体会
1、在此实验中，把在课堂上一些没搞懂的内容进行了细致的分析，得到冲激响应不变法设计数字滤波器的过程如下：
l 拟定模拟滤波器的系统函数 的技术指标；
l 根据技术指标设计 ，并将其写为 形式；
l 获得冲激响应不变法设计的数字滤波器的系统函数 为

2、根据采样定理可知，只有当模拟滤波器是带限时，即当 ， 时才有
假如模拟信号的频带不是限于 之间，则会在 的奇数倍附近产生频
率混叠，即冲激响应不变法的频率混叠现象。
3、本次对matlab的滤波器设计运用，对以后的通信原理仿真和理论知识有了进一步的了解，对本次通信原理课程设计做了提前预习。在设计滤波器时中，收获最大的是，掌握了冲激响应不变法设计数字滤波器的用法，可以很好的纯熟掌握滤波器的各种设计，对以后的编程有了很大的帮助，同时也对matlab产生了极大的爱好。
4