巴特沃斯Ⅱ型低通滤波器和切比雪夫Ⅱ型低通滤波器IIR低通数字滤波器设计文档格式.docx
- 文档编号:22512553
- 上传时间:2023-02-04
- 格式:DOCX
- 页数:18
- 大小:34.60KB
巴特沃斯Ⅱ型低通滤波器和切比雪夫Ⅱ型低通滤波器IIR低通数字滤波器设计文档格式.docx
《巴特沃斯Ⅱ型低通滤波器和切比雪夫Ⅱ型低通滤波器IIR低通数字滤波器设计文档格式.docx》由会员分享,可在线阅读,更多相关《巴特沃斯Ⅱ型低通滤波器和切比雪夫Ⅱ型低通滤波器IIR低通数字滤波器设计文档格式.docx(18页珍藏版)》请在冰豆网上搜索。
IIR低通数字滤波器设计
指导教师:
XX
设计要求:
用冲激响应不变法和双线性变换法分别以巴特沃斯滤波器、切比雪夫I型滤波器、切比雪夫U型滤波器为原型设计IIR低通数字滤波器。
1、IIR低通数字滤波器指标为:
■•p=0.2二数字通带截止频率(弧度)
■•s=0.3二数字阻带截止频率(弧度)
Rp=1dB通带衰减(dB);
As=15dB阻带衰减(dB)
2、采用冲激响应不变法分别设计上述三种数字滤波器;
3、分别绘制这三种数字滤波器的幅度响应曲线和相位响应曲线;
4、根据设计结果,对三种滤波器的性能进行比较和分析;
5、采用双线性变换法重新设计上述三种数字滤波器;
6分别绘制这三种数字滤波器的幅度响应曲线和相位响应曲线;
7、对三种滤波器的性能进行比较;
8、将其与使用冲激响应不变法的设计结果进行比较和分析
时间
设计内容
2011628
查阅资料,确定方案
2011629
设计总体方案
2011630
设计程序
2011.7.1
调试程序
2011.7.2
撰写报告
2011.7.5
答辩
周
星期一
星期二
星期三
星期四
星期五
课程设计评语表
指导教师评语:
成绩:
年月日
IIR低通数字滤波器设
设计目的和意义
目的:
1.深入理解数字信号处理基础知识的理解;
2.加深对MATLA基础知识的理解;
3.掌握低通数字滤波器的设计方法;
4.了解冲激响应不变法的基本原理和特点;
5.了解双线性变换法的基本原理和特点;
意义:
通过课程设计设计可以加深我们对课本基础知识的理解,对已经学习的知识进行实践训练,起到了理论联系实践的作用。
在设计过程中,一定会遇到很多的困难和问题,在解决问题的过程中,不仅锻炼了我解决实际问题的能力,而且也培养了我设计的综合能力。
总之,理论联系实践,对我来说是非常的重要。
IIR低通数字滤波器设计是滤波器设计中很经典的问题,而滤波器设计则是是《数字信号处理》的核心内容。
所以,IIR低通数字滤波器设计是数字信号处理的经典内容。
1.数字滤波器原理
与模拟滤波器类似,数字滤波器按频率特性划分为低通、高通、带通、带阻、全通等类型。
由于频率响应的周期性,频率变量以数字频率w来表示(w=QT=Q/fs,Q为模拟角频率,T为抽样时间间隔,fs为抽样频率),所以数字滤波器设计中必须给出抽样频率。
一般情况下,数字滤波器是一个线性移不变离散时间系统,利用有限精度算法来实现。
具体的实现方法有很多,不过主流的方法是:
先设计出对应的模拟滤波器,再将模拟滤波器数字化为数字滤波器。
2.脉冲响应不变法工作原理
冲激响应不变法遵循的准则是使数字滤波器的单位取样响应与参照的模拟滤波器的脉冲响应的取样值完全一样,即h(n)=ha(nT),其中T为取样周期。
从而使数字
滤波器在时域上模仿模拟滤波器,但是缺点是产生频率响应的混叠失真。
实际是由模拟滤波器转换成为数字滤波器,就是要建立模拟系统函数Ha(S)与
数字系统函数H(z)之间的关系,即S域到Z域的映射。
3.双线性变换法工作原理
使数字滤波器的频率响应与模拟滤波器的频率响应相似。
双线性变换法也是一种由S平面到Z平面的映射过程,双线性变换法与脉冲响应不变法不同,它是一种从S平面到z平面简单映射。
双线性变换中数字频率w和模拟频率Q之间的非线性关系使其只能应用于片段长度的函数,通常低通、高通、带通和带阻等滤波器等具有分段恒定的频率特性,可以采用预畸变的方法来补偿频率畸变,因此可以采用双线性变换设计方法。
三、详细设计步骤
(一)巴特沃斯滤波器
(I)冲级响应不变法
1.设计流程
(1)BUTTORD函数求滤波器的最小阶数N和截止频率OmgS;
(2)BUTTER函数求模拟滤波器响应函数的分子多项式系数B和分母多项式系数A;
(3)IMPINVAR函数求数字滤波器响应函数的分子多项式系数BZ和分母
多项式系数AZ;
(4)FREQZ函数求数字滤波器响应函数表达式;
2.实现过程及代码
%给出数字滤波器参数
Wp=0.2*pi;
Ws=0.3*pi;
Rp=1;
%单位是分贝,转换为10A(-Rp/20)
Rs=15;
%单位是分贝,转换为10A(-Rs/20)
T=1;
Fs=1/T;
%转换为模拟滤波器参数
OmgP=Wp*T;
%模拟通带截止频率
OmgS=Ws*T;
%模拟阻带截止频率
%模拟低通滤波器设计
%获得模拟滤波器阶数和通带截止频率
[N,OmgS]=buttord(OmgP,OmgS,Rp,Rs,'
s'
);
%设计模拟低通滤波器原型,得到模拟滤波器的零极点
[Z,P,K]=buttap(N);
%零极点转化为传递函数表示的滤波器
[B0,A0]=zp2tf(Z,P,K);
%将模拟低通滤波器原形转化为要求设计的低通滤波器
[B,A]=lp2lp(B0,A0,OmgS);
%数字滤波器设计
[Bz,Az]=impinvar(B,A,Fs);
%计算频率响应
[H,W]=freqz(Bz,Az);
%绘制图线
subplot(211)
plot(W/pi,abs(H));
xlabel('
频率'
ylabel('
幅度'
axis([0,1,0,1.1]);
set(gca,'
XTickMode'
'
manual'
XTick'
[00.20.30.51]);
YTickMode'
YTick'
[010A(-Rs/20)10A(-Rp/20)1]);
gridon
subplot(212)
plot(W/pi,angle(H)/pi);
title('
相频响应'
相位'
axis([0,1,-1,1]);
(U)双线性变换法
(1)BUTTORD函数求滤波器的最小阶数N和截止频率OmgS
(2)BUTTAP函数求模拟滤波器的零点和极点
(3)BILINEAR函数求数字滤波器零点和极点
OmgP=(2/T)*tan(Wp/2);
OmgS=(2/T)*tan(Ws/2);
[N,OmgS]=buttord(OmgP,OmgS,Rp,Rs,'
[Bz,Az]=bilinear(B,A,Fs);
幅频响应'
(二)切比雪夫I型滤波器
(1)CHEB1ORD函数求滤波器的最小阶数N和截止频率OmgS;
(2)CHEB1AP函数求模拟滤波器响应函数的分子多项式系数B和分母多项式系数A;
[N,OmgS]=cheb1ord(OmgP,OmgS,Rp,Rs,'
[Z,P,K]=cheb1ap(N,Rp);
[Bz,Az]=impinvar(B,A,Fs);
[010(Rs/20)10(Rp/20)1]);
(3)BILINEAR函数求数字滤波器响应函数的分子多项式系数BZ和分母
set(gca,'
gridon
(三)切比雪夫U型滤波器
(1)CHEB20RD函数求滤波器的最小阶数N和截止频率OmgS;
(2)CHEB2AP函数求模拟滤波器响应函数的分子多项式系数B和分母多项式系数A;
[N,OmgS]=cheb2ord(OmgP,OmgS,Rp,Rs,'
%设计模拟低通滤波器原型,得到模拟滤波器的零极点[Z,P,K]=cheb2ap(N,Rs);
[BO,AO]=zp2tf(Z,P,K);
[B,A]=lp2lp(BO,AO,OmgS);
subplot(212)plot(W/pi,angle(H)/pi);
(1)CHEB2ORD函数求滤波器的最小阶数N和截止频率OmgS;
[Z,P,K]=cheb2ap(N,Rs);
%绘制图线subplot(211)plot(W/pi,abs(H));
四、设计结果及分析
锁4-
t.-i
频率
相频响应
:
•;
五、体会
掌握了现代数字滤波器的基本设计理念,这是完成这个设计最有效的方法。
这个课程设计没有太大的难度,6个滤波器虽然不同,但是编程的方法都差不多,只是调用的函数和使用的方法有略微的不同。
通过这个课程设计,给我感触最深的是,设计的思想和方法是最重要的,至于具体的滤波器设计,知道了方法,就只需要查阅帮助文档,寻找相关的函数。
我总结了如下方法:
1.数字滤波器的参数
2.模拟滤波器的参数
3.模拟滤波器的设计
4.模拟滤波器转换为数字滤波器
六、参考文献
【1】《离散信号处理》奥本海姆,清华大学出版社
【2】《信号与系统基础》陈从颜、翟军永编著,机械工业出版社
【3】《电子滤波器设计》Arthur.B.Williams编著,科学出版社
【4】《数字信号处理教程》程佩青编著,清华大学出版社
【5】《数字信号处理及MATLA实现》恩德编著,机械工业出版社
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 巴特沃斯 型低通 滤波器 IIR 数字滤波器 设计