数字信号处理实验三

5IIR一、实验目二、实验计算机、软件环境三、实验基础理1原应h𝑎(𝑡)，h(n)等于h𝑎(𝑡)的取样值。2方a 思路H(s变换h(t采样h(nTh(n换a 将H𝑎(𝑠)𝐻(𝑠)=

对H𝑎(𝑠)

𝑘=1ℎ𝑎(𝑡)= 𝐴𝑘𝑒𝑝𝑘对h𝑎(𝑡)h(n)=

𝐴𝑘𝑒𝑝𝑘𝑛𝑇𝑢(𝑛𝑇)=

H(z)=

3、设计步

𝑘=1确定数字滤波器性能指标pstRpAs

根据指标pstRp和AsHaHa(sNNHa(s)k

A spkkp转换成数字极点epkTkNH(z) N可见H𝑎(𝑠)H(z)

k11epkT1

在中有两种方法可以实现上述变换1residueresiduez函数实现脉冲响应不变法，这两个函数的使用方法如[b,a]=[b,a]=方法2：提供了impinvar函数采用脉冲响应不变法实现模拟滤波器到数字滤波器babzaz，fs（1[bz,az]=impinvar(b,a)1四、实验内实验过程与实验结果调用impinvar函数采用脉冲响应不变法实现模拟滤波器到数字滤波器的变换，并利用freqz函数来计算数字滤波器的频率响应，同时利用grpdelay函数计算群延时响应：a=[12*wc2*wc^2wc^3];gridon;xlabel('\omega(\pi)');ylabel('|H(e^j^\omega)|');gridon;xlabel('\omega(\pi)');ylabel('|H(e^j^\omega)|');gridon;xlabel('\omega(\pi)'),ylabel('Groupdelay');bz az G4.547e-16+0.5813z^-1+0.2114z^-1-0.3984z^-1+0.2475z^-2-0.04321z^-Sampletime:0.1secondsDiscrete-timetransferfunction.1

1PhasePhaseof0

Group1Group0 要求分别采用、切比I型、切比II型和椭圆模拟原型滤波器及脉冲响应不IIR数字滤波器的优缺点及适用范围。实验过程与实验结果①grpdelay函数计算群延时响应：gridon;xlabel('\omega(\pi)');ylabel('|H(e^j^\omega)|');gridon;xlabel('\omega(\pi)');ylabel('|H(e^j^\omega)|');gridon;xlabel('\omega(\pi)'),ylabel('PhaseofH(e^j^\omega)');gridon;xlabel('\omega(\pi)'),ylabel('Groupdelay');bz az G0.000631z^-1+0.0101z^-2+0.01614z^-3+0.004101z^-4+ z^-1-3.364z^-1+5.068z^-2-4.276z^-3+2.107z^-4-0.5706z^-5+0.06607z^-Sampletime:0.1secondsDiscrete-timetransferfunction.Phase0Phase0②切比Ⅰ型滤波器的频率响应，同时利用grpdelay函数计算群延时响应：gridon;xlabel('\omega(\pi)');ylabel('|H(e^j^\omega)|');gridon;xlabel('\omega(\pi)');ylabel('|H(e^j^\omega)|');gridon;xlabel('\omega(\pi)'),ylabel('PhaseofH(e^j^\omega)');gridon;xlabel('\omega(\pi)'),ylabel('Groupdelay');bz az G0.005373z^-1+0.0181z^-2+0.003985z^-1-3.059z^-1+3.832z^-2-2.292z^-3+0.5495z^-Sampletime:0.1secondsDiscrete-timetransferfunction.0

1PhasePhaseof0

GroupGroup50 ③切比Ⅱ型滤波器的频率响应，同时利用grpdelay函数计算群延时响应：gridon;xlabel('\omega(\pi)');ylabel('|H(e^j^\omega)|');gridon;xlabel('\omega(\pi)');ylabel('|H(e^j^\omega)|');gridon;xlabel('\omega(\pi)'),ylabel('Phaseofgridon;xlabel('\omega(\pi)'),ylabel('Groupdelay');bz az G-0.4214+0.9724z^-1-0.8119z^-2+0.4205z^-3+1.663e-06z^-1-1.666z^-1+1.429z^-2-0.5193z^-3+0.09352z^-Sampletime:0.1secondsDiscrete-timetransferfunction.0

1PhasePhaseof0

Group0Group ④波器的频率响应，同时利用grpdelay函数计算群延时响应：gridon;xlabel('\omega(\pi)');ylabel('|H(e^j^\omega)|');gridon;xlabel('\omega(\pi)');ylabel('|H(e^j^\omega)|');gridon;xlabel('\omega(\pi)'),ylabel('PhaseofH(e^j^\omega)');gridon;xlabel('\omega(\pi)'),ylabel('Groupdelay');bz az G0.1855+0.04956z^-1+0.04956z^-2+0.1855z^-1-1.411z^-1+1.236z^-2-0.3549z^-Sampletime:0.1secondsDiscrete-timetransferfunction.10

1PhasePhaseof0

GroupGroup50 freqz函数来计算数字滤波器的频率响应，同时利用grpdelay函数a=[12*wc2*wc^2wc^3];gridon;xlabel('\omega(\pi)');ylabel('|H(e^j^\omega)|');gridon;xlabel('\omega(\pi)');ylabel('|H(e^j^\omega)|');gridon;xlabel('\omega(\pi)'),ylabel('PhaseofH(e^j^\omega)');gridon;xlabel('\omega(\pi)'),ylabel('Groupdelay');bz az G0.1057z^-1+0.06633z^-1-1.649z^-1+1.066z^-2-0.245z^-Sampletime:0.1secondsDiscrete-timetransferfunction.10

1PhasePhaseo