IIR滤波器设计与验证:冲激响应不变法实现

3星 · 超过75%的资源 需积分: 9 8 下载量 117 浏览量 更新于2024-09-13 收藏 80KB DOC 举报
本文档主要介绍了IIR(无限 impulse response)数字滤波器在数字信号处理系统中的实现方法,特别是使用冲击响应不变法来设计滤波器。首先,作者详细解释了冲激响应不变法的概念,这种方法是通过保持系统的冲激响应不变,从而将连续域的滤波器转换到离散域。在这个过程中,参数如通带最大衰减Rp、阻带最小衰减Rs以及采样频率FS都起到了关键作用。 计算巴特沃夫滤波器的阶数n和截止频率wn,利用`buttord`函数进行,根据指定的通带和阻带要求确定了滤波器的设计参数。接着,通过`butter`函数生成系统函数的分子向量Bs和分母向量As。为了进行离散化,使用`impinvar`函数进行脉冲不变型变换,得到新的系统函数Bz和Az。`filt`函数被用来实际应用这个滤波器于信号上,比如对于不同的正弦信号x1、x2和x3,分别进行滤波操作。 幅频特性和相频特性是评估滤波器性能的重要指标。通过`freqz`函数计算出滤波器的频率响应,然后用`plot`函数绘制幅度响应曲线和相位响应曲线,以可视化地展示滤波器的特性。最后,文档展示了实际应用滤波器对输入信号x1和x2进行滤波的结果,通过`stem`函数显示输入信号和经过滤波后的输出信号,直观地展示了滤波器在实际信号处理中的效果。 这篇文档涵盖了IIR数字滤波器设计的基本原理、具体实现步骤,以及如何运用它来处理模拟信号,并通过实例验证了滤波器的性能。这对于理解和使用数字信号处理技术,特别是在 DSP(数字信号处理器)中设计和应用滤波器是非常有帮助的。
2020-07-17 上传
巴特沃兹滤波器的vb.net程序。 ' 使用双线性变换法的 Butterworth 型 IIR 数字滤波器设计程序 ' ' 形参说明如下 : ' ' PbType ----------- 输入整型量 ,滤波器通带类型 : ' PbType = 0 : 低通滤波器 ; ' PbType = 1 : 高通滤波器 ; 2 ' PbType = 2 : 带通滤波器 ; ' PbType = 3 : 带阻滤波器 . ' fp1 ----------- 输入双精度量 , 低通或高通滤波器的通带边界频率 ( Hz ); 带通或带阻滤波器的通带低端边 ‘ 界频率 ( Hz ). ' fp2 ----------- 输入双精度量 , 带通或带阻滤波器的通带低端边界频率 ( Hz ). ' Apass ----------- 输入双精度量 , 通带衰减 ( dB ). ' fs1 ----------- 输入双精度量 , 低通或高通滤波器的阻带边界频率 ( Hz ); 带通或带阻滤波器的阻带高端边 ‘ 界频率 ( Hz ). ' fs2 ----------- 输入双精度量 , 带通或带阻滤波器的阻带高端边界频率 ( Hz ). ' Astop ----------- 输入双精度量 , 阻带衰减 ( dB ). ' fsamp ----------- 输入双精度量 , 采样频率 ( Hz ). ' points ----------- 输入整型量 , 幅频特性计算点数 . ' ord ----------- 输入整型量 , 滤波器阶数 . ' NumSec( ) -------- 输出双精度量 , 转移函数二阶节的分子多项式系数二维数组 . ' 元素 NumSec( k, i ) 中 , ' k : 二阶节序号 ; ' i : 多项式系数 , i = 0 相应于常数项 . ' DenSec( ) -------- 输出双精度量 转移函数二阶节的分母多项式系数二维数组 . ' 元素 DenSec( k, i ) 中 , ' k : 二阶节序号 ; ' i : 多项式系数 , i = 0 相应于常数项 . ' NumSec_Z( ) ------ 输出双精度量 系统函数二阶节的分子多项式系数二维数组 . ' 元素 NumSec_Z( k, i ) 中 , ' k : 二阶节序号 ; ' i : 多项式系数 , i = 0 相应于常数项 . ' DenSec_Z( ) ------ 输出双精度量 系统函数二阶节的分母多项式系数二维数组 . ' 元素 DenSec_Z( k, i ) 中, ' k : 二阶节序号 ; ' i : 多项式系数 , i = 0 相应于常数项 . ' AR( ) ------------ 输出双精度量 ,滤波器的幅频特性数组 . ' Sub Butterworth(PbType As Integer, fp1 As Double, fp2 As Double, Apass As Double, fs1 As Double, fs2 As Double, Astop As Double, fsamp As Double, points As Integer, ord As Integer, NumSec() As Double, DenSec() As Double, NumSec_Z() As Double, DenSec_Z() As Double, AR() As Double) Dim i%, j%, k%, ord_t% Dim angle#, emp1#, temp2#, temp3# Dim ratio(0 To 50) As Double '''''''''''''''''''' If PbType = 0 Then ' 低通滤波器 ; wpass = 2# * Pi * fpass / fsamp: wstop = 2# * Pi * fstop / fsamp ' 通带、阻带边界频率 omikaP = Tan(wpass / 2#): omikaS = Tan(wstop / 2#) epass = epson(Apass): estop = epson(Astop) ' 根据对幅频特性的技术要求 ,计算模拟滤波器的阶数 orde = Ne_B(estop, epass, omikaS, omikaP) ord = Fix(orde) + 1 omk0 = omika0(omikaP, epass, ord) ' 调用 Fz_LP 子程序, 将低通模拟滤波器的转移函数变量 s 映射为低通数字滤波器的系统函数变量 z Call Fz_LP(F1(), F2(), ord_t) ' End If '''''''''''''''''''' If PbType = 1 Then ' 高通滤波器 ; wpass = 2# * Pi * fpass / fsamp: wstop = 2# * Pi * fstop / fsamp ' 通带、阻带边界频率 3 omikaP = 1# / Tan(wpass / 2#): omikaS = 1# / Tan(wstop / 2#) epass = epson(Apass): estop = epson(Astop) ' 根据对幅频特性的技术要求 ,计算模拟滤波器的阶数 orde = Ne_B(estop, epass, omikaS, omikaP) ord = Fix(orde) + 1 omk0 = omika0(omikaP, epass, ord) ' 调用 Fz_HP 子程序,将高通模拟滤波器的转移函数变量 s 映射为高通数字滤波器的系统函数变量 z Call Fz_HP(F1(), F2(), ord_t) End If '''''''''''''''''''' If PbType = 2 Then ' 带通滤波器 ; wp1 = 2# * Pi * fp1 / fsamp: wp2 = 2# * Pi * fp2 / fsamp ' 通带上下边界频率 ws1 = 2# * Pi * fs1 / fsamp: ws2 = 2# * Pi * fs2 / fsamp ' 阻带上下边界频率 Ci = BpC(wp1, wp2) omikaP = Abs((Ci - Cos(wp2)) / Sin(wp2)) omikaS1 = Abs((Ci - Cos(ws1)) / Sin(ws1)) omikaS2 = Abs((Ci - Cos(ws2)) / Sin(ws2)) If omikaS1 <= omikaS2 Then omikaS = omikaS1 Else omikaS = omikaS2 End If epass = epson(Apass): estop = epson(Astop) ' 根据对幅频特性的技术要求 ,计算模拟滤波器的阶数 orde = Ne_B(estop, epass, omikaS, omikaP) ord = Fix(orde) + 1 omk0 = omika0(omikaP, epass, ord) ' 调用 Fz_BP 子程序,将带通模拟滤波器的转移函数变量 s 映射为带通数字滤波器的系统函数变量 z Call Fz_BP(fp1, fp2, fsamp, F1(), F2(), ord_t) End If