数字信号处理课程设计-基于MATLAB的巴特沃斯低通滤波器以切贝雪夫低通滤波器的设计.doc

各专业全套优秀毕业设计图纸吉首大学信息科学与工程学院课程设计报告书课 程 ： 数字信号处理教程 课 题： 基于MATLAB的巴特沃斯低通滤波器以 切贝雪夫低通滤波器的设计 姓 名： 学 号： 专 业： 通信工程 年 级： 2012级 指导教师： 基地指导教师： 2014年 12 月一、项目介绍与设计目的1. 通过实验加深对巴特沃斯低通滤波器以及切贝雪夫低通滤波器基本原理的理解。2. 学习编写巴特沃斯低通滤波器和切比雪夫的MATLAB仿真程序。3.滤波器的性能指标如下：通带截止频率fp=3kHz，通带最大衰减Rp=2dB，阻带截止 频率fst=6kHz，阻带最小衰减As=30db二、设计方案1项目环境要求 MATLAB 软件2 设计内容一理论设计： 模拟巴特沃思低通滤波器的设计 1.有技术指标可求的设计参数 p=6000rad/s; st=12000rad/s; Rp=2db; As=30 db 2.求N; =5.369 取N=63. 确定参数 =1971.114. 求系统函数Han（s）=1/(1+3.8637033S+7.4641016S2+9.141620S3+7.4641016S4+3.8637033S5+S6)5.去归一化 H（s）=Han（s/c） =0.7648/(0.7648+3.0910s5+6.2424s4+7.9947S3+6.6300S2+0.6949)化简可得 Ha(s)=5.8650*1025/(s6+7.6158*104s5+2.9000*109s4+7.0010*1013s3+1.1267*10*18s2+1.1496*1032s+5.8650*1025)切比雪夫低通滤波器1. 由技术指标可得e =1/(100.1As)=0.031N=arcch/arch()取N=4系统函数为：Ha(s)=0.0316s4+3.5954*108S2+5.1099*1017/(S4+6.6533*104S3+2.2247*109S2+4.3659*1013S+5.1099*1017)二程序设计(1) 巴特沃斯低通滤波器程序：clc,clear all;OmegaP=2*pi*3000;OmegaS=2*pi*6000;Rp=2;As=30;N=ceil(log10(10(As/10)-1)/(10(Rp/10)-1)/(2*log10(OmegaS/OmegaP)OmegaC=OmegaP/(10(Rp/10)-1)(1/(2*N);z0,p0,k0=buttap(N);p=p0*OmegaC;a=real(poly(p)k=k0*OmegaCN;b0=real(poly(z0);b=k*b0w0=OmegaP,OmegaS;H,w=freqs(b,a);Hx=freqs(b,a,w0);dbHx=-20*log10(abs(Hx)/max(abs(H)plot(w/(2*pi)/1000,20*log10(abs(H);xlabel(f(kHz);ylabel(dB);axis(-1,12,-55,1)set(gca,xtickmode,manual,xtick,0,1,2,3,4,5,6,7,8,9,10);Set(gca,ytickmode,maunal,ytick,-50,-40,-30,-2,-10,0);grid;输出结果N=6b=5.8650e+25a=1 7.5158e+04 2.9000e+09 7.0010e+13 1.1267e+18 1.1496e+22 5.8650e+25dbHx=2.0000 33.7692在matlab中显示如下：程序运行结果如下故系统函数为：Ha(s)=5.8650*1025/(s6+7.6158*104s5+2.9000*109s4+7.0010*1013s3+1.1267*10*18s2+1.1496*1032s+5.8650*1025)(2) 切贝雪夫型低通滤波器程序：clc，clear all;OmegaR=2*pi*3000;OmegaS=2*pi*6000;Rp=2;As=30;g=sqrt(10(As/10)-1)/(10(Rp/10)-1);OmegaR=OmegaS/OmegaPN=ceil(log10(g+sqrt(g*g-1)/log10(OmegaR+sqrt(OmegaR*OmegaR-1)OmegaC=OmegaS;z0,p0,k0=cheb2ap(N,As);a0=real(poly(p0)aNn=a0(N+1);p=p0*OmegaC;a=real(poly(p)aNu=a(N+1);b0=real(poly(z0);M=length(b0);bNn=b0(M);z=z0*OmegaC;b=real(poly(z);bNu=b(M);k=k0*(aNu*bNn)/(aNn*bNu);b=k*bw0=OmegaP,OmegaS;H,w=freqs(b,a);Hx=freqs(b,a,w0);dbHx=-20*log10(abs(Hx)/max(abs(H)plot(w/(2*pi)/1000,20*log10(abs(H);xlabel(f(kHz);ylabel(dB);axis(-1,12,-55,1)set(gca,xtickmode,manual,xtick,0,1,2,3,4,5,6,7,8,9,10);Set(gca,ytickmode,maunal,ytick,-50,-40,-30,-20,-10,0);grid;输出结果：N=4b=0.0316 0 3.5954e+08 0 5.1099e+17A=1 6.6533e+04 2.2247e+09 4.3659e+13 5.1099e+17程序运行结果如下图此低通滤波器的系统函数为：Ha(s)=0.0316s4+3.5954*108S2+5.1099*1017/(S4+6.6533*104S3+2.2247*109S2+4.3659*1013S+5.1099*1017)三、总结和分析 通过本次butterworth低通滤波器的设计，使我对低通数字滤波器的工作原理和特性有了深刻的认识和了解，实验主要用到了matlab软件。对matlab的一些运用也有了一定的了解,对以后的继续学习有了一定的基础。在这次实验中通过查阅资料和结合平时所学的知识才得以完成。四 参考资料网页：关于matla R2010b的安装数字信号处理教程第四版 程佩青指导老师评语：成绩评定： 指导老师签名： 基地指导老师签名： 年 月 日教研室意见教研室主任签章年 月 日学院意见分管院长签章年 月 日9