基于最小方差低通FIR滤波器课程设计说明书

1、中北大学课 程 设 计 说 明 书学生： 杰 学 号：22学 院：信息与通信工程学院 专 业：生物医学工程 题 目： 基于最小方差低通FIR滤波器 指导教师： 权 职称: 副教授 2014 年 1 月 3 日 中北大学课程设计任务书 13/14 学年第 一 学期学 院： 信息与通信工程学院 专 业： 生物医学工程 学 生 姓 名：21学 号： 田野 学 生 姓 名： 22学 号： 杰 学 生 姓 名：45 学 号： 黄志浩 课程设计题目： 基于最小方差低通FIR滤波器 起 迄 日 期： 2013年12 月23 日2014年1月 3 日 课程设计地点： 院楼机房 指 导 教 师： 权 系 主 任

2、：下达任务书日期: 2013 年12月23 日课 程 设 计 任 务 书1设计目的：根据所学的数字信号处理和MATLAB相关知识，用最小方差法设计一个低通FIR滤波器。2设计容和要求（包括原始数据、技术参数、条件、设计要求等）：实现最优化设计，调节这些零点的分布，使得实际滤波器的频率响应Hd(ej)与理想滤波器的频率响应Hd(ej)之间的最大绝对误差最小。3设计工作任务与工作量的要求包括课程设计计算说明书(论文)、图纸、实物样品等：（1）查阅文献资料，了解相关知识；（2）确定设计方法；（3）编写MATLAB程序，并上机调试；（4）撰写课程设计说明书。课 程 设 计 任 务 书4主要参考文献：

3、1 阎石.数字电子技术基础.：高等教育，1998； 2 王远.模拟电子技术.：机械工业，2001； 3 汝全.电子技术常用器件应用手册.：机械工业，2003； 4 毕满清.电子技术实验与课程设计.：机械工业，2006； 5 俊谟.单片机中级教程-原理与应用.：航空航天大学，2002；5设计成果形式与要求：提供课程设计说明书一份；MATLAB程序；6工作计划与进度： 2013年12月23日 12 月25 日：查阅资料；12月26日 12 月 27 日：确定设计方案；12月28 日 12月31 日：编写程序，上机调试；2014年1月 1 日 1 月2 日：完成课程设计说明书；1月 3日：答辩。系主

4、任审查意见：签字：年月日14 / 23目 录1 设计目标·············································&#

5、183;··························12低通FIR滤波器技术指标·····················

6、·····································13低通FIR滤波器的设计···········

7、;·················································13.1低

8、通FIR滤波器阶数的估计················································

9、······13.2最小方差线性相位的误差·········································

10、83;··············13.3 参数、L的确定·································

11、83;···················23.4 参数、的确定····························

12、3;··························33.5参数的确定······················&

13、#183;············································34用直接型结构实现···

14、83;·················································

15、83;··········115 用FDATOOL分析·····································&

16、#183;····························116 误差分析····················

17、;··················································

18、;··126.1 误差产生的原因·············································

19、3;··················126.2 误差的理论计算·····························

20、83;··································126.3用FDATOOL分析不同字长对其幅频响应和相频响应的影响··········

21、··················137 总结·······························

22、·············································148 参考文献···

23、3;·················································

24、3;··················151 设计目标根据所学的数字信号处理和MATLAB相关知识，用最小方差法设计一个低通FIR滤波器。从FIR数字滤波器的系统函数可以看出，极点都是在z平面的原点，而零点的分布是任意的。不同的分布将对应不同的频率响应，最优化设计实际上就是调节这些零点的分布，使得实际滤波器的频率响应Hd(ej)与理想滤波器的频率响应Hd(ej)之间的最大绝对误差最小。2 低通FIR滤波器技术指标（

25、通带截止频率） （阻带截止频率） （通带衰减） （阻带衰减） （通带最大衰减） （阻带最小衰减）3 低通FIR滤波器的设计3.1 低通FIR滤波器阶数的估计46 由于N为偶数，所以可以设计一个1型的低通FIR滤波器。3.2 最小方差线性相位的误差对于基于最小方差的线性相位FIR滤波器的设计下面式子为误差的简化为 其中是低通FIR的振幅响应，是要求的振幅响应，是权重函数。由于所有四种类型的线性相位FIR滤波器的振幅响应可以表示为cos(wk)3.3 参数、L的确定cos(wk) 式中、L的确定a 的确定 由于不同类型也就不尽一样，不同类型时的表达式如下=1 对于1型=cos() 对于2型=sin

26、() 对于3型=sin() 对于4型 由于我们 设计的低通FIR滤波器为1型所以 =1 b 的确定 同样根据不同的类型其的表达式也不一样= 对于1型= 对于2型= 对于3型= 对于4型 我们选择=，对于1型,，C L 的确定 L=M 对于1型 L= 对于2型 L= M-1 对于3型 L= 对于4型根据N与M的关系 M=,表示取不大于的最大整数，所以M=22,L=22。3.4 参数、的确定中、的确定 根据最小方差的相关要求可知 在通带中 在阻带中 在通带中 在阻带中3.5 参数的确定根据上面式子可以确定，L的值和的表达式，由于最小方差是滤波器参数的一个函数。为了得到的最小值，令 由它可生成（L+

27、1）个等式的线性方程组，用来求解。 我们考虑1型线性相位FIR滤波器的设计。在这种情况下，=1，=且L=22。则均方误差的表达式为若有，式中。计算如下1 求H， 1-15时取1，16-22时取0将0到0.35上取均匀的22点最后求的H=1 0.9987 0.9950 0.9888 0.980 0.968 0.955 0.939 0.921 0.900 0.877 0.853 0.825 0.796 0.765 0.732 0.697 0.660 0.622 0.582 0.540 0.498 ;1 0.995 0.980 0.955 0.921 0.877 0.825 0.765 0.697

28、0.622 0.540 0.454 0.363 0.268 0.170 0.07 0.028 0.128 0.226 0.322 0.415 0.588;1 0.988 0.955 0.900 0.825 0.732 0.622 0.498 0.363 0.219 0.071 0.078 0.226 0.369 0.504 0.627 0.737 0.829 0.937 0.655 0.212 0.282 0.707 0.959 0.977 0.755 0.349 0.142 0.599 0.909 0.997 0.841 0.479;1 0.852 0.454 0.078 0.587 0.9

29、23 0.987 0.760 0.309 0.233 0.707 0.972 0.951 0.649 0.156 0.382 0.809 0.997 0.891 0.522 0 0.522 ;1 0.825 0.363 0.226 0.736 0.989 0.897 0.492 0.085 0.633 0.959 0.951 0.610 0.057 0.516 0.909 0.985 0.717 0.199 0.389 0.841 0.999;1 0.796 0.268 0.369 0.856 0.994 0.727 0.163 0.467 0.907 0.977 0.649 0.057 0.

30、558 0.946 0.949 0.564 0.049 0.644 0.975 0.909 0.473;1 0.765 0.170 0.504 0.940 0.937 0.491 0.184 0.774 0.999 0.755 0.156 0.516 0.946 0.932 0.479 0.198 0.783 0.999 0.746 0.142 0.528; 1 0.732 0.071 0.627 0.989 0.821 0.212 0.510 0.959 0.894 0.349 0.349 -0.909 0.948 0.478 0.247 0.841 0.984 0.599 0.107 0.

31、755 0.999;0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;0 0 0 0 0 0 0 0 0 0 0 0 0

32、0 0 0 0 0 0 0 0 0算这个H时算出来的负值取了它的绝对值2 求令d= 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0然后最小均方解可以根据求解线性方程得到。令 a=a=n*pinv(r)用MATLAB可以计算出 ac=a'c =1.0e+005 * Columns 1 through 190.95392.5572-1.4336-0.37351.9754-3.75723.5858-2.6165-2.31652.7508-3.48832.3659-3.81731.90290.5901 0 0 0 0 Columns 20 through

33、 22 0 0 0= 0.95392.5572-1.4336-0.37351.9754-3.75723.5858-2.6165-2.31652.7508-3.48832.3659-3.81731.90290.5901 0 0 0 0 0 0 03 最终的结果由于N的限制所以所设计的滤波器为1型所以=1 =c L=M=22根据表达式 cos(wk)=0.9539+2.5572-1.4336-0.3735+1.9754-3.7572+3.5858-2.6165-2.3165+2.7508-3.4883+2.3659-3.8173+1.9029+0.590116到22项为0 这里我们只是求出了幅频特

34、性，但由于其相频特性是确定的所以在设计中不考虑其相频特性。 求=0.0916()+0.1171()-0.6515-0.1867+0.9877-1.8786+1.7929-1.3082-1.1582+1.3754-1.7441+1.1829-1.9086+0.9514+0.29504 计算误差 由和可求出误差用MATLAB计算误差如下h=1 0.9987 0.9950 0.9888 0.980 0.968 0.955 0.939 0.921 0.900 0.877 0.853 0.825 0.796 0.765 0.732 0.697 0.660 0.622 0.582 0.540 0.498

35、;1 0.995 0.980 0.955 0.921 0.877 0.825 0.765 0.697 0.622 0.540 0.454 0.363 0.268 0.170 0.07 0.028 0.128 0.226 0.322 0.415 0.588;1 0.988 0.955 0.900 0.825 0.732 0.622 0.498 0.363 0.219 0.071 0.078 0.226 0.369 0.504 0.627 0.737 0.829 0.904 0.957 0.989 0.999 ;1 0.980 0.921 0.825 0.697 0.904 0.989 0.988

36、 0.897 0.727 0.492 0.212 0.0085 0.376 0.633 0.833 0.959 0.999 ;1 0.939 0.765 0.498 0.170 0.177 0.504 0.769 0.942 0.999 0.937 0.760 0.277 0.163 0.184 0.510 0.774 0.944 0.999 0.934 0.755 0.485;1 0.921 0.697 0.363 0.028 0.415 0.736 0.942 0.998 0.897 0.655 0.309 0.008 0.467 0.774 0.959 0.993 0.870 0.610

37、 0.254 0.142 0.516; 1 0.900 0.622 0.219 0.226 0.627 0.903 0.999 0.897 0.616 0.212 0.233 0.633 0.906 0.999 0.894 0.610 0.205 0.240 0.638 0.909 0.999 ;1 0.877 0.541 0.071 0.415 0.800 0.989 0.937 0.655 0.212 0.282 0.707 0.959 0.977 0.755 0.349 0.142 0.599 0.909 0.997 0.841 0.479;1 0.852 0.454 0.078 0.5

38、87 0.923 0.987 0.760 0.309 0.233 0.707 0.972 0.951 0.649 0.156 0.382 0.809 0.997 0.891 0.522 0 0.522 ;1 0.825 0.363 0.226 0.736 0.989 0.897 09 -0.909 0.948 0.478 0.247 0.841 0.984 0.599 0.107 0.755 0.999;0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;0 0 0 0

39、0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0a= 0.9539 2.5572 -1.4336 -0.3735 1.9754 -3.7572 3.5858 -2.6165 -2.3165 2.7508 -3.4883 2.3

40、659 -3.8173 1.9029 0.5901 0 0 0 0 0 0 0;s=h*a'd= 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0;c=s-d'e=c*c'e=Columns 1 through 193.6235 2.4192 -0.7942 -1.7213 1.1543 2.2088 1.9128 0.4541 1.2903 2.4710 3.5210 0.6483 -4.8485 -0.3134 -9.8727 0 0 0 02.4192 1.6152 -0.5303 -1.1492 0.7707 1.4747

41、1.2771 0.3032 0.8615 1.6498 2.3508 0.4329 -3.2372 -0.2093 -6.5916 0 0 0 0-0.7942 -0.5303 0.1741 0.3773 -0.2530 -0.4841 -0.4193 -0.0995 -0.2828 -0.5416 -0.7718 -0.1421 1.0628 0.0687 2.1640 0 0 0 0-1.7213 -1.1492 0.3773 0.8177 -0.5483 -1.0493 -0.9086 -0.2157 -0.6129 -1.1738 -1.6726 -0.3080 2.3032 0.14

42、89 4.6899 0 0 0 01.1543 0.7707 -0.2530 -0.5483 0.3677 0.7036 0.6093 0.1447 0.4110 0.7871 1.1216 0.2065 -1.5445 -0.0998 -3.1450 0 0 0 02.2088 1.4747 -0.4841 -1.0493 0.7036 1.3464 1.1660 0.2768 0.7865 1.5063 2.1463 0.3952 -2.9556 -0.1911 -6.0182 0 0 0 01.9128 1.2771 -0.4193 -0.9086 0.6093 1.1660 1.009

43、7 0.2397 0.6811 1.3044 1.8586 0.3422 -2.5594 -0.1655 -5.2116 0 0 0 00.4541 0.3032 -0.0995 -0.2157 0.1447 0.2768 0.2397 0.0569 0.1617 0.3097 0.4413 0.0813 -0.6076 -0.0393 -1.2373 0 0 0 01.2903 0.8615 -0.2828 -0.6129 0.4110 0.7865 0.6811 0.1617 0.4595 0.8799 1.2538 0.2309 -1.7265 -0.1116 -3.5156 0 0 0

44、 02.4710 1.6498 -0.5416 -1.1738 0.7871 1.5063 1.3044 0.3097 0.8799 1.6850 2.4011 0.4421 -3.3064 -0.2137 -6.7325 0 0 0 03.5210 2.3508 -0.7718 -1.6726 1.1216 2.1463 1.8586 0.4413 1.2538 2.4011 3.4214 0.6300 -4.7114 -0.3046 -9.5934 0 0 0 00.6483 0.4329 -0.1421 -0.3080 0.2065 0.3952 0.3422 0.0813 0.2309

45、 0.4421 0.6300 0.1160 -0.8675 -0.0561 -1.7665 0 0 0 0-4.8485 -3.2372 1.0628 2.3032 -1.5445 -2.9556 -2.5594 -0.6076 -1.7265 -3.3064 -4.7114 -0.8675 6.4878 0.4194 13.2106 0 0 0 0-0.3134 -0.2093 0.0687 0.1489 -0.0998 -0.1911 -0.1655 -0.0393 -0.1116 -0.2137 -0.3046 -0.0561 0.4194 0.0271 0.8540 0 0 0 0-9

46、.8727 -6.5916 2.1640 4.6899 -3.1450 -6.0182 -5.2116 -1.2373 -3.5156 -6.7325 -9.5934 -1.7665 13.2106 0.8540 26.8998 0 0 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0Columns 20 through 220 0 00 0 00 0 00 0 00 0 00 0 00 0