数字滤波器课程设计

ﻩﻩﻩﻩ课程设计课程设计名称:数字信号处理课程设计专业班级 :电信１203学生姓名：刘海峰学号： 2指导教师:乔丽红课程设计时间:2015/07/01-2015/07/06电子信息工程专业课程设计任务书说明:本表由指导教师填写,由教研室主任审核后下达给选题学生,装订在设计(论文）首页一. 技术要求⏹双线性变换法设计切比雪夫II型数字IIR低通滤波器，⏹要求通带边界频率为40０Hｚ，⏹阻带边界频率分别为500Hｚ，⏹通带最大衰减1dB,⏹阻带最小衰减4０dB，⏹抽样频率为２000Hｚ，二. 设计原理IＩＲ滤波器的设计包括三个步骤：①给出所需要的滤波器的技术指标；②设计一个H（ｚ）使其逼近所需要的技术指标:③实现所设计的Ｈ（z),ＩIR数字滤波器设计的最通用的方法是借助于模拟滤波器的设计方法。

所以ＩＩR数字低通滤波器的设计步骤是：①按一定规则将给出的数字滤波器的技术指标转换为模拟低通滤波器的技术指标；②根据转换后的技术指标设计模拟低通滤波器G（s）:③再按一定规则将G(ｓ）转换成H（z）。

在此过程中,我们用到了很多MＡＴLＡB中的函数，如设计切比雪夫低通滤波器的函数afｄ_ｃhebl、由直接型转换为级联型的函数ｄir2ｃａs、双线性变换的函数bilineaｒ等.其中afｄ _ｃhｅbl用于实现用模拟指标设计一个低通模拟滤波器，ｂilｉnｅar用于利用双线性变换法将模拟低通滤波器转换为数字低通滤波器。

三.程序流程图四：源代码（完美版)％归一化低通滤波器技术指标cｌc；clear all；Ap=１； %最大通带衰减As=40；％最小阻带衰减Ｗ=2０００; ％抽样周期Wp=4０0; ％通带边界频率Ws＝500；％阻带边界频率wp=２*pi*Wp/Ｗ； %归一化通带边界频率ws=2*pi*Ｗｓ/W； %归一化阻带边界频率Wp1=taｎ（wp/2)；%模拟低通滤波器通带边界频率Wｓ1=ｔan(ｗs/2)；％模拟低通滤波器阻带边界频率%归一化切比雪夫II型低通模拟滤波器［N,Wn］=cheｂ2ord（Ｗp1，Ｗｓ1,Aｐ,As,'s'); ％确定滤波器阶数和频率尺度缩放因子［BT,AＴ］＝ｃｈeby２（N,As，Wn，＇ｓ’）;%传输函数的系数[Z,P，K］=cheｂ2aｐ(N，As）；％最小阻带衰减为Ａs(DB）的Ｎ阶归一化模拟切比雪夫２型低通滤波器的零点、极点和增益因子[H，W]=ｚp2tf(Z,Ｐ，K)；%传输函数有理化形式ｆigure；[P,Q］＝fｒeqs(Ｈ，W)；frｅqs(H，W);％模拟滤波器的Ｈ（jｗ）的复频域响应．拉普拉斯格式．(自动挑选２0０个频率点来计算频率响应)fiｇｕrｅ;subｐlot(1,1,１）；fk=0:120０0/512：１2000;wk=2＊pi＊fk；Hk＝frｅqs(BT，ＡＴ,ｗk);％挑选wk个频率点来计算频率响应plot（fk/１00０，２0＊loｇ１０（abs（Ｈk））);gｒid oｎ;ｘlabel('频率’)yｌabel(＇增益衰减’）％用双线性变换法将H(s）转换成数字滤波器H（ｚ）[nｕm,ｄｅｎ］=bｉlｉnear（BT，AＴ,0。

5);%复变量映射［z，p，k］=tｆ2zｐ（ｎum，deｎ）；％显示传输函数disp（＇分子系数：’)；dｉsp（num);ｄisp（’分母系数:’）；disｐ(den)；%计算增益响应w=0:pｉ／255:pi；h=frｅｑz(nuｍ，deｎ,w）;ｇ=２0＊log10（abｓ（h)）；％计算增益衰减fiｇurｅ；ｐｌｏｔ(ｗ/pｉ,ｇ)；grid on；％绘制切比雪夫低通滤波器幅频特性axｉs（［０ 1 －１00 1]）;xlabel（'\oｍegａ/\pi'）;yｌabｅl(＇增益/dB')；tｉtle（'切比雪夫2型低通滤波器幅频响应曲线’)；figurｅ；zｐlaｎe(z,p）；%绘制极零图ａxｉｓ（[－2 2 －２ 2］）；titlｅ（'零极点图’)；figｕｒe；sｕbploｔ（1，１,1）,ploｔ(w/pｉ, angle（h)/pｉ）；grｉd ｏn；xlabel（’归一化角频率＇);yｌabel（'相位响应'）; ％绘制切比雪夫低通滤波器相频特性axis（[0 1 —1 1]）;%输入信号技术指标ｆ1=3００；f２=600;t＝0：0.000５：1；x1=ｓin(2＊pｉ＊ｆ１*t）；％X1信号x2＝sin（2＊pi＊f2*t）；%X2信号ｘ=x１+x2；％X信号ｆｉgurｅ；subploｔ（2，１,1);plot（x１)；ｇrid ｏｎ；％绘制X１波形图axiｓ（[０，5０*pｉ，－3，3］）；xlaｂel(’t'）；yｌabeｌ(＇x1'）；ｔｉtｌe（'x1的波形'）；subpｌot（2,１，2）;plot（ｘ2）;griｄ oｎ； %绘制X2波形图axiｓ([0,５0＊pi,-3，3］);xｌaｂeｌ（＇t＇);ylaｂeｌ('x2＇）；ｔｉｔle(’x2的波形’）；figｕrｅ;ｓubplot（２,１，１);plot（x）；grid ｏn； %绘制X波形图axis(［0，50＊pi,-3，３]）；xlaｂeｌ（＇t');ylaｂel（'x＇）;title（’输入x的波形')；y=fiｌｔｅr(num,deｎ,x)；％Ｘ信号通过切比雪夫ＩI 数字低通滤波器输出Ysubｐlot（2，1,2)pｌｏt（y）；ｇｒｉｄｏn；％绘制输出Ｙ的波形图ａxis（［0,50＊ｐｉ,-3,3］）;xlaｂｅl（'t'）；ylａbｅｌ('y＇）；ｔitle(’滤波器输出y的波形'）；figure；subplot(2,1，1）;plot（abｓ（fft(ｘ））)；grid ｏn； %绘制输入信号频域波形ａxiｓ(［０，８00,0,30０］）;xlabel（'频率'）;ｙlabel(’幅度'）；ｔitｌe（'滤波器输入x的频域波形＇）；sｕｂpｌot（2，1,2);ｐlｏt(abs(ｆｆｔ(y））)；griｄｏn; ％绘制输出信号频域波形axis([0，８０0，０,300］);ｘlabel(’频率');ｙlａbel（’幅度’);tiｔle('滤波器输出y的频域波形’);五.仿真结果图1:模拟滤波器复频域响应(幅频响应—相位响应）1010101-200-1000100200Frequency (rad/s)P h a s e (d e g r e e s )10-11010110-610-410-2100Frequency (rad/s)M a g n i t u d e图2:切比雪夫2数字低通滤波器幅频响应曲线0.10.20.30.40.50.60.70.80.91ω/π增益/d B切比雪夫2型低通滤波器幅频响应曲线图3：传输函数极零图-2-1.5-1-0.500.51 1.52-2-1.5-1-0.500.511.52Real PartI m a g i n a r y P a r t零极点图图４：相频响应曲线0.10.20.30.40.50.60.70.80.91-1-0.8-0.6-0.4-0.200.20.40.60.81归一化角频率相位响应图５：输入信号波形曲线50100150-202tx 1x1的波形50100150-202tx 2x2的波形图6:输入信号经过滤波器输出信号时域波形图50100150-202tx输入x 的波形50100150-202ty 滤波器输出y 的波形图7:输入信号经过滤波器输出信号频域波形图0100200300400500600700800100200300频率幅度滤波器输入x 的频域波形010*******400500600700800100200300频率幅度滤波器输出y 的频域波形图９：传输函数系数如下六。

结论分析1。

通带边界归一化角频率0.4＊ｐi 、阻带边界归一化角频率０.5*pi. 【复合技术要求】0.10.20.30.40.50.60.70.80.91-100-90-80-70-60-50-40-30-20-100ω/π增益/d B切比雪夫2型低通滤波器幅频响应曲线2。

通带边界频率400Hz ，阻带边界频率500Hz,当输入信号频率为300Hz 和600H ｚ时，300Ｈz 可以通过该滤波器，６00H ｚ被滤除。

【复合技术要求】3．对应频域波形010*******40050060070080050100150200250300频率幅度滤波器输入x 的频域波形010*******40050060070080050100150200250300频率幅度滤波器输出y 的频域波形【复合技术要求】 七．参考资料⏹ 程佩青著，《数字信号处理教程》，清华大学出版社，20０1 ⏹ Ｓa ｎjit Ｋ。

Mit ｒa 著,孙洪,余翔宇译,《数字信号处理实验指导书(ＭA ＴLAB 版)》,电子工业出版社，2０05年1月 ⏹ 郭仕剑等，《MAT ＬAB 7．x 数字信号处理》，人民邮电出版社，200６年⏹ 胡广书,《数字信号处理 理论算法与实现》,清华大学出版社，２003年 八.实验心得1．本次课程设计再次让我感受到ＭAT ＬA Ｂ软件功能之强大，今后一定要加强对MA ＴLAB 等工具性软件的学习2.平时只是理论上接触过切比雪夫2数字滤波器，在课程设计一始真正接触到他的时候还真有点手足无措的感觉,由此看出只学习理论知识还是不行的，要注重实践锻炼,从各方面提高自己的能力3。

不过静下心来认真思考,看书,上网查阅相关资料后,就眉目清晰了，而且越来越上手，用了不到两天的时间就完全弄懂了！4.遇到的几个问题经过上机实验时老师的耐心指点很随意就解决了,从而顺利完成了本次课程设计实践，衷心感谢老师的认真倾听和耐心陪伴!。