怎样用窗函数截断IIR滤波器脉冲响应

[wangzhuzhen]

用FIR逼近已已知IIR，采用直接截断法！已经知道了IIR的传递函数，并且使用impz求出脉冲响应，接下来怎么用窗函数截断脉冲响应？MATLAB语言程序如何编写？
请问我下面的做法是否正确：
1。根据传递函数分子分母向量b，a计算传递函数的频率响应；并画出频率响应图。
2。根据IIR的传递函数用impz函数计算其脉冲响应h；
3.对第2步求得的h用kaiser窗函数右乘得到与IIR滤波器传递函数近似的FIR滤波器传递函数的系数向量c。
我用上述方法求得的c，再使用【H,f】=freqz（c，1），求fir的频率响应并绘图plot（fabs（H）），得到的图像与原IIR滤波器传递函数频率响应曲线图像几乎一致，这能说明我的做法是正确的吗？