切比雪夫I型数字带通滤波器

课程设计课程名称：数字信号处理题目编号：0801题目名称：切比雪夫I型带通IIR数字滤波器专业名称：电子信息工程班级：1203班学号：20124470323学生姓名：段超任课教师：陈忠泽2015年08月30日1目录1.数字滤波器的设计任务及要求（编号0801）...................22.数字滤波器的设计及仿真...................................22.1数字滤波器（编号0801）的设计...................................22.2数字滤波器（编号0801）的性能分析...............................43.数字滤波器的实现结构对其性能影响的分析...................63.1数字滤波器的实现结构一（0801）及其幅频响应......................73.2数字滤波器的实现结构二（0801）及其幅频响应......................93.3数字滤波器的实现结构对其性能影响的小结........................124.数字滤波器的参数字长对其性能影响的分析..................124.1数字滤波器的实现结构一（0801）参数字长及幅频响应特性变化......144.2数字滤波器的实现结构二（0801）参数字长及幅频响应特性变化.......174.3数字滤波器的参数字长对其性能影响的小结........................205.结论及体会..............................................205.1滤波器设计、分析结论..........................................215.2我的体会......................................................215.3展望..........................................................2121.数字滤波器的设计任务及要求（1）切比雪夫I型带通IIR数字滤波器各项指标如下：阻带下截止频率1s=radid32)ln(;通带下截止频率1p=radiidd20)(32)ln(log10;通带上截止频率2p=radiidd20)(32)ln(-1log10;阻带上截止频率2s=radid32)ln(-1;通带最大衰减p=1dB;阻带最小衰减s=60dB.其中的id为抽到的题目的四位数编号，我的题目编号是0801，所以取id为801.（2）滤波器的初始设计通过手工计算完成；（3）在计算机辅助计算基础上分析滤波器结构对其性能指标的影响（至少两种合适的滤波器结构进行分析）；（4）在计算机辅助计算基础上分析滤波器参数的字长对其性能指标的影响；（5）以上各项要有理论分析和推导、原程序以及表示计算结果的图表；（6）课程设计结束时提交说明书。2.数字滤波器的设计及仿真2.1数字滤波器的设计(1)根据任务,确定性能指标:在设计带通滤波器之前,首先根据工程实际的需要确定滤波器的技术指标:（手工计算完成）3通带下截止频率1p=0.354π通带上截止频率2p=0.646π阻带下截止频率1s=0.208π阻带上截止频率2s=0.792π阻带最小衰减s=60dB通带最大衰减p=1dB抽样频率sF=8kHz（2）有所给的技术参数可直接得到各模拟频率。1p=1p*sF,2p=2p*sF,1s=1s*sF,2s=2s*sF1p=2832,2p=5168,1s=1664,2s=6336(3)根据教程表7.8的变换关系，设归一化（_p）低通滤波器的阻带衰减用表示，则有pB=2p-1p=2336，0p=21pp=3825.67322st=pstpstB20222=1.7231st=pstpstB10122=-3.053取st_=min（|1st|,|2stst2|）=1.723,可满足阻带衰减要求。（4）求波纹系数=1101.0p=0.50885（5）根据st_及给定的p，s，利用教程（7.5.65）式，可得切比雪夫I型低通滤波器的阶次为4pstarccharcchs/1102.01=)723.1()2147.1965(arccharcch=1398.12765.8=26114.7取N=8.（6）查教程表7.5，可得N=8，p=1dB时的归一化原型切比雪夫I型低通滤波器的系统函数)(sHan)(sHan=01723.01073.04478.08468.08369.16552.14230.29198.0123456782118ssssssss(7)按教表7.8的相应变换关系，求出模拟带通滤波器系统函数为：ppsBssanbpsHsH202)()(65614441108955.1100.1101722.7100.8ss（8）按冲激响应不变法的修正公式（7.78）式，求所需数字滤波器的系统函数)(zH87654321765435847.282727.80350.269831.57331.156014.27911.55316.010013.00015.00015.00003.00001.0zzzzzzzzzzzzz161514131211109121110984301.02275.01159.34457.14690.102324.42727.213787.70004.00012.00005.00024.00038.0zzzzzzzzzzzzz2.2数字滤波器的性能分析MATLAB程序如下:5Rs=1;As=60;Fs=8000;T=1/Fs;wp1=0.354*pi;wp2=0.646*pi;ws1=0.208*pi;ws2=0.792*pi;Omegap1=wp1*Fs;Omegap2=wp2*Fs;Omegap=[Omegap1,Omegap2];Omegas1=ws1*Fs;Omegas2=ws2*Fs;Omegas=[Omegas1,Omegas2];bw=Omegap2-Omegap1;w0=sqrt(Omegap1*Omegap2);[N,OmegaC]=cheb1ord(Omegap1,Omegas1,Rp,As,'s')N=8OmegaC=8.8970e+03[z0,p0,k0]=cheb1ap(N,Rp);ba1=k0*real(poly(z0));aa1=real(poly(p0));[ba,aa]=lp2bp(ba1,aa1,w0,bw)；[bd,ad]=impinvar(ba,aa,Fs)；%(C,B,A)=dir2par(bd,ad)[H,w]=freqz(bd,ad);dbH=20*log10((abs(H)+eps)/max(abs(H)));subplot(2,2,1),plot(w/pi,abs(H));subplot(2,2,2),plot(w/pi,angle(H)/pi);subplot(2,2,3),plot(w/pi,dbH);subplot(2,2,4),zplane(ba,ad);相应的图形6图1频率特性曲线和零极点图3.数字滤波器的实现结构对其性能影响的分析在理想状态下，对于同一个传递函数几乎对应着无数种等效结构，然而这些结构却并不一定都能实现。在无限参数字长的情况下，所有能实现传递函数的结构之间，其表现完全相同。然而在实际中，由于参数字长有限的限制，各实现结构的表现不同。下面我们就对比直接型和级联型两种结构在本设计中对性能指标的影响。在MATLAB中可以利用FDATOOL工具箱建构不同类型的数字滤波器。再次为了使对比效果明显，将上述初步设计的切比雪夫带通IRR数字滤波器的设计参数字长（即转移函数中分子分母各项前的系数）进行保留小数点后3位的进一步缩减。缩减参数如下：7Nun=[000000.002-0.002-0.0010.004-0.002-0.0010.0010000]Den=[1-0.5325.791-2.60115.733-5.98326.035-8.27328.585-7.37921.273-4.23210.469-1.4463.116-0.2280.430]将上述参数输入FDATOOL中的filtercoefficients工具中3.1数字滤波器的实现结构一及其幅频响应直接型结构滤波器的转移函数为：)(zH87654321765585.28273.8035.26983.5733.15601.2791.5532.01001.0002.0002.0zzzzzzzzzzz161514131211109111098430.0228.0116.3446.1469.10232.4273.21379.7001.0001.0002.0004.0zzzzzzzzzzzz差分方程：)16(430.0)15(228.0)14(116.3)13(446.1)12(469.10)11(232.4)10(273.21)9(379.7)8(585.28)7(273.8)6(035.26)5(983.5)4(733.15)3(601.2)2(791.5)1(532.0)11(001.0)10(001.0)9(002.0)8(004.0)7(001.0)6(002.0)5(002.0)(nynynynynynynynynynynynynynynynynxnxnxnxnxnxnxny8直接一型的结构流图如下：x(n)y(n)5z0.002-0.5321z1z-0.0025.7911z…………1z0.0010.4301z图2Direct-FormI滤波器结构流程图选择filterstructure选项框中的Direct-FormI选项，点击窗口下方的ImportFilter按钮，构建直接1型结构的椭圆带通IRR数字滤波器，结果如图所示。图3Direct-FormI型结构滤波器幅频响应图9表1Direct-FormI结构滤波器性能指标影响(1p，2p，1s，2s单位为rad/sample；s,p,单位为dB）由图3和表1可以看出：（1）滤波器幅频曲线在通带和阻带内波动幅度不均匀；（2）阻带最小衰减与通带最大衰减都存在误差；（3）wp和ws分别较初始设计变化了。3.2数字滤波器的实现结构一及其幅频响应级联型结构滤波器的转移函数：性能指标初始设计DirectFormI相对误差（%）1p0.3540.3513184-0.0026818-0.762p0.6460.6492920.0032920.511s0.2080.29748540.0894854432s0.7920.7630615-0.0289385-3.65p10.9859966-0.01400341.4s6061.991721.991723.310)9643.08642.011108.13048.11)(9629.07747.016328.04550.11()9608.08774.011)(8924.03369.011)(8809.04108.016304.01()8759.06715.011283.1)(8638.05222.01)(8103.00328