基于最小方差带阻FIR的设计.docVIP

基于最小方差带阻（模拟频变）FIR的设计 电气工程学院 电子071班 20074470137 翁强 设计目的 （1）、掌握用最小方差优化设计法设计FIR滤波器的原理及方法，熟悉相应的MATLAB编程。 （2）、熟悉线性相位FIR滤波器的幅频特性和相频特性。 带阻滤波器的技术指标 1、 2、（通带衰减） （阻带衰减） 设计流程图 模拟低通滤波器的最小方差实现 （1）、由带阻滤波器的技术指标写出相应的低通滤波器的技术指标如下： （通带截止频率） （阻带截止频率） （通带衰减） （阻带衰减） （通带最大衰减） （阻带最小衰减） （2）、低通FIR滤波器的设计 1 低通FIR滤波器阶数的估计 46 由于N为偶数，所以可以设计一个1型的低通FIR滤波器。 2对于基于最小方差的线性相位FIR滤波器的设计下面式子为误差的简化为 其中是低通FIR的振幅响应，是要求的振幅响应，是权重函数。由于所有四种类型的线性相位FIR滤波器的振幅响应可以表示为 cos(wk) 3 的确定 由于不同类型也就不尽相同，不同类型时的表达式如下 =1 对于1型 =cos() 对于2型 =sin() 对于3型 =sin() 对于4型 由于我们设计的低通FIR滤波器为1型所以 =1 4 的确定 同样根据不同的类型其的表达式也不一样 = 对于1型 = 对于2型 = 对于3型 = 对于4型 我们选择=，对于1型,， 5 L 的确定 L=M 对于1型 L= 对于2型 L= M-1 对于3型 L= 对于4型 根据N与M的关系 M=,表示取不大于的最大整数，所以M=22,L=22。 6 根据上面式子可以确定，L的值和的表达式，由于最小方差是滤波器参数的一个函数。为了得到的最小值，令 由它可生成（L+1）个等式的线性方程组，用来求解。 我们考虑1型线性相位FIR滤波器的设计。在这种情况下，=1，= 且L=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 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.62