基于频率抽样法的FIR数字低通滤波器的设计.docVIP

基于频率抽样法的FIR数字低通滤波器的设计 导读：就爱阅读网友为您分享以下“基于频率抽样法的FIR数字低通滤波器的设计”的资讯，希望对您有所帮助，感谢您对92的支持! 基于频率抽样法的FIR数字低通滤波器的设计 1 设计目的 熟悉频率采样法的理论及其应用；掌握频率采样法设计FIR数字滤波器的方法。了解FIR数字滤波器的频率特性和相位特性，观察过渡带取样点对滤波器幅频特性的影响。掌握用频率采样法设计线性相位FIR低通数字滤波器的方法，并掌握该方法的matlab编程和仿真。 2 FIR数字滤波器设计的原理 2.1频率抽样设计法 FIR低通滤波器的设计一般方法有两种，即频率抽样法和窗函数法，频率抽样法设计不同于窗函数法，窗函数是从时域出发，把理想的hd(n)用一定形状得窗函数截取成有限长的h(n)，以此h(n)来近似理想的hd(n)，这样得到的频率响应H(ejw)逼近于所要求的理想的频率响应Hd(ejw)。 频率抽样法则是从频域出发，把给定的理想频率响应Hd(ejw)加以等间隔抽样，即Hd(ejw)|w?2?Nk?Hd(k)然后以此Hd(k)作为实际FIR数字滤波器的频率特性 w?2?Nk的抽样值H(k)，即令H(k)?Hd(k)?Hd(ejw)|知道H(k)k?0,1,?,N?1， 后，由DFT定义，可以用频域的这N个抽样值H(k)来唯一确定有限长序列h(n)，而由X(z)的内插公式知道，利用这N个频域抽样值H(k)同样可求得FIR滤波器的系统函数H(z)及频率响应H(ejw)。这个H(z)或H(ejw)将逼近Hd(z)或 Hd(ejw)，H(z)和H(ejw)的内插公式为 H(z)?1?zN?NN?1?1?Wk?0H(k)?kNz?1 （2.1） N?1H(ejw)??k?0H(k)?(w?2?N k) （2.2） 其中?(w)是内插函数?(w)?1Nsin(wN2w2)e?jw(N?12) （2.3） sin()将式（2.3）代入（2.2）式，化简后可得 H(ejw)?1Ne?j(N?12)wN?1?H(k)ek?0?j?kNsin(sin(w2wN2?w)?kN? （2.4） ))]即 H(e)?ejw?j(N?12)wN?1?H(k)?Nk?01ej?kN(N?1)sin[N(?k?2Nw?ksin(?)2N （2.5） 从内插公式（2.2）看到，在各频率抽样点上，滤波器的实际频率响应是严格地和理想频率响应数值相等，即H(ej2?Nk)?H(k)?Hd(k)?Hd(ej2?Nk)。但是在 抽样点之间的频率响应则是由各抽样点的加权内插函数的延伸叠加而形成的，因而有一定的逼近误差，误差大小取决于理想频率响应曲线形状，理想频率响应特性变化越平缓，则内插值越接近理想值，逼近误差越小，如图2.1梯形理想频率特性所示。反之，如果抽样点之间的理想频率特性变化越陡，则内插值与理想值之误差就越大，因而在理想频率特性的不连续点附近，就会产生肩峰和波纹，如图2.2矩形理想频率特性所示。 图2.1 梯形理想频率特性 图2.2 矩形理想频率特性 2.2线性相位的约束 对于第二类线性相位FIR滤波器，由于h(n)偶对称、N为偶数时, 公式是 H(e)?H(w)ejw?j(N?12)w，其中幅度函数H(w)应为奇对称的，H(w)??H(2??w)， 如果抽样值H(k)?H(ej2?Nk)也用幅值Hk（纯标量）与相角Qk表示，则为： j2?NkH(k)?H(e)?H(2?N2?Nk)ejQK?Hke1NjQK （2.6） 其中Qk必须为： Qk??(HkN?12)k???k(?1 ) （2.7） 必须满足奇对称，即Hk??HN?k。 2.3线性相位第一种频率抽样 N?1 H(k)??n?0*h(n)eN?j2?nk （2.8） 当h(n)为实数时，满足H(k)?H((N?k))NRN(k)?H*(N?k)，由此得出也就是说，H(k)的模H(k)以k?N2为H(k)?H(N?k)，Q(k)??Q(N?k)， 对称中心呈偶对称，H(k)的相角Q(k)以k?N2为对称中心呈奇对称。再利用线性相位的条件Q(ejw)??N?12： w，即可得到（N为偶数） ?2??N?1??N??kk?0,?,?1?????N?