IIR高通滤波器的设计.docx
- 文档编号:4295680
- 上传时间:2022-11-29
- 格式:DOCX
- 页数:11
- 大小:245.80KB
IIR高通滤波器的设计.docx
《IIR高通滤波器的设计.docx》由会员分享,可在线阅读,更多相关《IIR高通滤波器的设计.docx(11页珍藏版)》请在冰豆网上搜索。
IIR高通滤波器的设计
湖南文理学院实习报告
实习名称:
IIR高通滤波器的设计
教学院部:
电气与信息工程学院
专业班级:
通信工程08101班
学生姓名:
奉繁荣
学生学号:
200816020147
指导教师:
朱明旱
完成时间:
2011年06月20日
报告成绩:
评阅意见:
评阅教师日期
基于双线性变换设计IIR高通滤波器
一、设计目的
(1)加深对滤波器的理解,更好的学习数字信号处理这门课程。
(2)加强对matlab的使用,更好的为以后毕业设计打下基础。
二、设计要求
(1)掌握双线性变换法设计iir高通滤波器的设计。
(2)要有设计的源程序和matlab仿真的结果
(3)要严谨的、要有理有据的编写课程设计报告并及时的提交给老师检查
三、设计原理
3.1数字滤波器的分类
3.1.1从实现的网络结构划分:
无限脉冲响应(IIR)
有限脉冲响应(FIR)
3.1.2数字滤波器的技术要求
是幅频特性,它表示信号通过该滤波器后各频率成分衰减情况,一般选频滤波器的技术要求由它给出。
是相频特性,反映各频率成分通过滤波器后在时间上的延时情况,如果对输出波形有要求,就要考虑相频特性的技术指标。
3.2IIR数字滤波器的设计方法频率变换法:
借助于模拟滤波器的设计
3.2.1设计的框图
3.2.2设计的方法
直接设计法:
直接在频域或时域上进行设计
数字化的方法
冲激响应不变法
阶跃响应不变法
双线性变换法:
数字频率与模拟频率之间非线性,需要在设计之前进行预畸。
3.3双线性变换的原理
3.3.1双线性变换法的映射规则
3.4设计的具体要求
3.4.1设计要求参数和题目
采用双线性变换法设计一巴特沃斯高通滤波器,其通带截止频率为3khz,阻带上限截止频率为1.5khz,通带衰减不大于3db,阻带衰减不小于16db,采样频率为10khz
3.4.2设计的思路
IIR数字滤波器的设计借助模拟滤波器原型,再将模拟滤波器转换成数字滤波器,这些过程已经成为一整套成熟的设计程序。
模拟滤波器的设计已经有了一套相当成熟的方法,它不但有完整的公式,而且还有较为完整的图表查询,因此,充分利用这些已有的资源将会给数字滤波器的设计带来很大的方便。
已知数字低通滤波器的设计要求wp,ws,Rp和As,首先设计一个等效的模拟滤波器,然后再将它映射为所期望的数字滤波器来确定H(z)。
对这个过程所要求的步骤是:
1)选取T并确定模拟频率:
利用设计参数Ωp,Ωs,Rp和As,设计一个模拟滤波器Ha(s);这可以利用模拟滤波器的原型来完成;3)再将
(1)代入Ha(s),求出H(z)。
3.5理论设计的参数计算
3.5.1求对应数字的频率
3.5.2求常数
采用归一化(
)原型低通滤波器作为变换的低通原型,则低通到高通的变换所需的
3.5.3求低通原型
。
设
为满足数字高通滤波器的归一化原型模拟低通滤波器的阻带起始截止频率,则可按
的预畸变换关系来求,得
3.5.4求阶次N。
按阻带衰减求原型归一化模拟低通滤波器的阶次N。
有巴特沃斯低通滤波器频率响应公式去对数,即
其中
解得
取N=3。
3.5.5求归一化巴特沃斯低通原型的
。
去N=3,查表可得
为
3.5.6求数字滤波器的系统函数
最终求得
将
代入,可求得
四、原程序清单
%任务指标
fph=3*1000;%模拟域通带截止频率
fsh=1.5*1000;%模拟域阻带截止频率
fs=3*1000;%采样频率
rp=3;rs=16;%通带衰减,阻带衰减
omegaph=fph/(fs/2);%通带截止频率归一化
omegash=fsh/(fs/2);%阻带截止频率归一化
[N,wn]=cheb1ord(omegaph,omegash,rp,rs)%调用cheb1ord工具函数得到滤波器阶数和截止频率
[BZ,AZ]=cheby1(N,rp,wn,'high')%调用cheby1函数得高通滤波器传递函数系数
ht=impz(BZ,AZ);%调用impz函数画滤波器冲击响应
n=0:
511;%输入信号离散化
t=n/fs;
x=2*sin(2*pi*1000*t)+5*sin(2*pi*3000*t);%输入信号包含f=1000hz,f=3000hz的正弦波
y=filter(BZ,AZ,x);%调用该滤波器得该输出信号
%滤波器频响特性曲线
figure
(1);
freqz(BZ,AZ,512,1.5*10000)
gridon;
%滤波器冲击响应曲线
figure
(2);
plot(ht)
gridon;
%输入信号时域波形图
figure(3);
axis([0,0.01,-7,7]);
plot(t,x);
title('输入信号');
gridon;
%输出信号离散序列
figure(4);
stem(y,'.');
title('输出序列');
gridon;
%输出信号时域波形图
figure(5);
ya=y*sinc(fs*(ones(length(n),1)*t-(n/fs)'*ones(1,length(t))));
plot(t,ya);axis([0,0.01,-7,7]);
title('输出波形');
gridon;
%输入信号频谱图
figure(6);
N1=512;
y=fft(x,N1);
f=(0:
length(y)-1)'*fs/length(y);
mag=abs(y);
stem(f,mag);
title('N1=512点');
%输出信号频谱图
figure(7);
yb=fft(ya,N1);
f=(0:
length(yb)-1)'*fs/length(yb);
mag=abs(yb);
stem(f,mag);
title('N1=512点');
五、仿真结果
N=
3
wn=
0.4333
BZ=
0.0204-0.12250.3063-0.40840.3063-0.12250.0204
AZ=
1.00000.66681.58331.13200.95520.46840.1951
图1滤波器频率响应曲线
图2滤波器冲击响应曲线图
图3输入信号时域波形
图4输出序列图5输出信号时域波形
图6滤波器输入信号频谱图
图7滤波器输出信号频谱图
六.设计的心得和体会
通过这次课程设计,使我掌握了课程设计的基本思路和方法,掌握了课程设计说明书的基本撰写方法,了解了基本的设计思想和设计方法。
同时此次课程设计使我对《数字信号处理原理》这门课程的基础知识和基本理论有了更深的理解和掌握,锻炼了我们综合运用所学知识的能力,并在理论分析设计、计算制图运用标准和规范查阅设计手册与资料以及计算机应用能了等方面得到了初步的训练和提高,培养了我们严谨求实的科学态度。
同时这次课程设计采用分组的形式进行,锻炼了我们的分工合作能力,同时培养了我们的团队意识。
这次课程设计使我进一步熟悉了计算工具软件---MATLAB.并进一步掌握了MATLAB的使用方法.对MATLAB语言的发展和特点有了更深的了解,熟悉其工作环境,并掌握了MATLAB的基本语法,在其应用方面也有了更深的了解,了解了MATLAB程序书写的一般思路及一般步骤,学会了用MATLAB解决复杂信号处理等问题的方法。
使我对MATLAB信号处理工具箱内的模型函数有了深刻的体会和应用。
MATLAB作为一种科学计算软件,有着其独特优势,尤其在科学计算及精确绘图上,与其他计算机语言相比,简洁而且具有智能化,适应我们的思维方式与书写习惯,它使用解析方法工作,直观且简单明了。
人机交互性能好,操作也很简单,在图形描绘上解决手工描绘的各种缺点。
在课程设计的这段时间里,我认为收获还是很多的,不但进一步掌握了数字信号处理的基础知识及一门专业仿真软件的基本操作,还提高了自己的设计能力及动手能力,同时对于模拟滤波器来了个系统的总结。
更多的是让我看清了自己,明白了凡事需要耐心,实践是检验学习的唯一标准。
理论知识的不足在这次课设中表现的很明显。
这将有助于我今后的学习,端正自己的学习态度,从而更加努力的学习。
总之,通过本次课程设计不但让我又学到了一些知识,而且也提高了我的综合能力。
使我在各方面都得到了锻炼,非常感谢我的同组同学,也非常感谢我们的指导老师,使我们这次的课程设计任务圆满完成
七.参考文献
[1].《数字信号处理教程(第三版)》,陈佩青,清华大学出版社。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- IIR 滤波器 设计
![提示](https://static.bdocx.com/images/bang_tan.gif)