《基于最小方差法低通fir的设计大学毕设论文》由会员分享，可在线阅读，更多相关《基于最小方差法低通fir的设计大学毕设论文（14页珍藏版）》请在金锄头文库上搜索。
1、基于最小方差低通FIR滤波器设计说明书 （一）设计目标 根据所学的数字信号处理和MATLAB相关知识，用最小方差法设计一个低通FIR滤波器。从FIR数字滤波器的系统函数可以看出，极点都是在z平面的原点，而零点的分布是任意的。不同的分布将对应不同的频率响应，最优化设计实际上就是调节这些零点的分布，使得实际滤波器的频率响应Hd(ej)与理想滤波器的频率响应Hd(ej)之间的最大绝对误差最小。（二）低通FIR滤波器技术指标1 （通带截止频率） （阻带截止频率） （通带衰减） （阻带衰减） （通带最大衰减） （阻带最小衰减）(三)低通FIR滤波器的设计3.1 低通FIR滤波器阶数的估计 46 由于N为
2、偶数，所以可以设计一个1型的低通FIR滤波器。3.2 对于基于最小方差的线性相位FIR滤波器的设计下面式子为误差的简化为 其中是低通FIR的振幅响应，是要求的振幅响应，是权重函数。由于所有四种类型的线性相位FIR滤波器的振幅响应可以表示为 cos(wk)3.3 cos(wk) 式中、L的确定a 的确定 由于不同类型也就不尽相同，不同类型时的表达式如下 =1 对于1型 =cos() 对于2型 =sin() 对于3型 =sin() 对于4型 由于我们 设计的低通FIR滤波器为1型所以 =1 b 的确定 同样根据不同的类型其的表达式也不一样 = 对于1型 = 对于2型 = 对于3型 = 对于4型 我
3、们选择=，对于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+1）个等式的线性方程组，用来求解。 我们考虑1型线性相位FIR滤波器的设计。在这种情况下，=1，=且L=22。则均方误差的表达式为 若有 ，式中。计算如下1 求H， 1-15时取1，16-22时取0将0到0.3
4、5上取均匀的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 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
5、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.541 0.363 0.170 0.028 0.226 0.415 0.588 0.737 0.856 0.942 0.989 0.998 0.967 0.897 0.792 0.654 0.492 ;1 0.968 0.877 0.732 0.540 0.315 0.071 0.177 0.415 0.627 0.800 0.924
6、0.989 0.99 0.937 0.821 0.655 0.312 0.212 0.035 0.282 0.510 ; 1 0.955 0.825 0.622 0.363 0.071 0.226 0.504 0.737 0.904 0.989 0.988 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
7、 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 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.4
8、15 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.587 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 0.492 0.085 0.633 0.959 0.951 0.
9、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.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.74
10、6 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.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
11、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算这个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=ac = 1.0e+005 * Columns 1 thro
12、ugh 19 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.5901 0 0 0 0 Columns 20 through 22 0 0 0= 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.5901 0 0 0 0 0 0 03 最终的结果由于N的限制所以所设计的
13、滤波器为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.5901 16到22项为0 这里我们只是求出了幅频特性，但由于其相频特性是确定的所以在设计中不考虑其相频特性。 求=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
