利用matlab设计巴特沃斯低通滤波器.doc

三峡大学 课 程 设 计 报 告 专业班级 课 程 数字信号处理课程设计 学 号 2009142116 学生姓名 姜祥奔 指导教师 王露 2012年 5 月 平时成绩（20%） 报告成绩（40%） 答辩成绩（40%） 总成绩 数字信号处理课程设计 实验一：用双线性变换法和脉冲响应不变法设计IIR数字滤波器 采用双线性变换法设计一个巴特沃斯数字低通滤波器，要求：通带截止频率100Hz，阻带截止频率200Hz，通带衰减指标Rp小于2dB，阻带衰减Rs大于15dB，滤波器采样频率Fs=500Hz。绘制频率响应曲线。 理论部分：（原理及设计过程） 第一步：论的关系及数字域性能的公式表示。模拟频率与数字频率之间为线性关系，为抽样周期，故 修正后 修正后 按衰减的定义 设处频率响应幅度归一化为1,即，则上两式变成 ………………（1） ………………（2） 这就是数字滤波器的性能指标的表达式。 2 下面把数字低通滤波器的性能要求转变为“样本”模拟低通滤波器的性能要求。由，按修正式，设没有混叠效应（即混叠效应设计完成后再进行校验） 则有 ……………（3） 利用（3）式，由（1）、（2）式可写出模拟低通滤波器的指标为 …….（4） …（5） 3 计算“样本”模拟低通滤波器所需的阶数N及3dB截止频率。巴特沃思低通滤波器的幅度平方函数是 以分贝形式表示上式，即 ……………（6） 把求出的性能指标关系（4）式、(5)式代入(6)式得 先用等号来满足指标，可得 ＝ ＝ 解此两方程，得N=1.3709, N是滤波器阶次，必须取整数，为了满足或超过给定指标，故应选取比求出的N大一点的整数，故取N＝2,代入通带条件（6-118a）式，得＝1.038,即 N＝2, ＝1.038 ……………（7） 显然，利用（7）式的N和，阻带指标正好满足，而通带指标则低于2dB，这正好对减少冲激响应不变法造成的频谱混叠效应是有利的。 程序部分： wp=100*2*pi; %数字滤波器的通带截止频率 ws=200*2*pi; %数字滤波器的阻带截止频率 Rp=2;As=15; %输入滤波器的通阻带衰减指标 %转换为模拟滤波器指标 Fs=500;T=1/Fs; wp1=wp*T; ws1=ws*T; Omgp=(2/T)*tan(wp1/2); %原型通带频率预修正 Omgs=(2/T)*tan(ws1/2); %原型通阻频率预修正 %模拟原型滤波器计算 [n, Omgc]=buttord(Omgp, Omgs,Rp,As,s); %计算阶数n和截止频率 [z0,p0,k0]=buttap(n); %归一化原型设计 [ba,aa]=zp2tf(z0,p0,k0); %把滤波器零极点转化为传递函数模型 [ba1,aa1]=lp2lp(ba,aa, Omgc); %变换为模拟低通滤波器系数b,a %或者[ba1,aa1]=butter(n, Omgc,’s’); %用双线性变换法计算数字滤波器系数 [bd,ad]=bilinear(ba1,aa1,Fs) ; %用双线性变换法求数字滤波器系数b,a %求数字系统的频率特性 [H,w]=freqz(bd,ad); plot(w*Fs/(2*pi),20*log10(abs(H))); ylabel(增益/dB);xlabel(频率（Hz）);title(数字滤波器幅度响应);grid on; 数字滤波器的频率响应图： 实验二：用窗函数法设计FIR数字滤波器