基于matlab利用双线性变换法设计iir带阻滤波器

基于matlab利用双线性变换法设计iir带阻滤波器内容摘要：

拓分量的频谱交叠，即产生频率响应的混叠失真，如图74所示。 这时数字滤波器的频响就不同于原模拟滤波器的频响，而带有一定的失真。 当模拟滤波器的频率响应在折叠频率以上处衰减越大、越快时，变换后频率响应混叠失真就越小。 这时，采用脉冲响应不变法设计的数字滤波器才能得到良好的效果。 图12脉冲响应不变法中的频响混叠现象 对某一模拟滤波器的单位冲激响应ha(t)进行采样，采样频率为fs，若使fs增加，即令采样时间间隔（T=1/fs）减小，则系统频率响应各周期延拓分量之间相距更远，因而可减小频率响应的混叠效应。 脉冲响应不变法优缺点：从以上讨论可以看出，脉冲响应不变法使得数字滤波器的单位脉冲响应完全模仿模拟滤波器的单位冲激响应，也就是时域逼近良好，而且模拟频率Ω和数字频率ω之间呈线性关系ω=ΩT。 因而，一个线性相位的模拟滤波器（例如贝塞尔滤波器）通过脉冲响应不变法得到的仍然是一个线性相位的数字滤波器。 脉冲响应不变法的最大缺点是有频率响应的混叠效应。 所以，脉冲响应不变法只适用于限带的模拟滤波器(例如，衰减特性很好的低通或带通滤波器)，而且高频衰减越快，混叠效应越小。 至于高通和带阻滤波器，由于它们在高频部分不衰减，因此将完全混淆在低频响应中。 如果要对高通和带阻滤波器采用脉冲响应不变法，就必须先对高通和带阻滤波器加一保护滤波器，滤掉高于折叠频率以上的频率，然后再使用脉冲响应不变法转换为数字滤波器。 当然这样会进一步增加设计复杂性和滤波器的阶数。 IIR数字滤波器是一种离散时间系统，其系统函数为假设M≤N，当M＞N时,系统函数可以看作一个IIR的子系统和一个(MN)的FIR子系统的级联。 IIR数字滤波器的设计实际上是求解滤波器的系数和 ，它是数学上的一种逼近问题，即在规定意义上（通常采用最小均方误差准则）去逼近系统的特性。 如果在S平面上去逼近，就得到模拟滤波器；如果在z平面上去逼近，就得到数字滤波器。 ⑵用双线性变换法设计IIR数字滤波器脉冲响应不变法的主要缺点是产生频率响应的混叠失真。 这是因为从S平面到Ｚ平面是多值的映射关系所造成的。 为了克服这一缺点，可以采用非线性频率压缩方法，将整个频率轴上的频率范围压缩到π/T～π/T之间，再用z=esT转换到Z平面上。 也就是说，第一步先将整个S平面压缩映射到S1平面的π/T～π/T一条横带里；第二步再通过标准变换关系z=es1T将此横带变换到整个Z平面上去。 这样就使S平面与Z平面建立了一一对应的单值关系，消除了多值变换性，也就消除了频谱混叠现象，映射关系如图13所示。 图13双线性变换的映射关系为了将S平面的整个虚轴jΩ压缩到S1平面jΩ1轴上的π/T到π/T段上，可以通过以下的正切变换实现（15）式中,T仍是采样间隔。 当Ω1由π/T经过0变化到π/T时，Ω由∞经过0变化到+∞，也即映射了整个jΩ轴。 将式（15）写成将此关系解析延拓到整个S平面和S1平面，令jΩ=s，jΩ1=s1，则得再将S1平面通过以下标准变换关系映射到Z平面z=es1T从而得到S平面和Z平面的单值映射关系为：（16）（17）式（16）与式（17）是S平面与Z平面之间的单值。