资源描述:
《无限脉冲响应数字滤波器iir设计》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、无限脉冲响应数字滤波器(IIR)设计所谓的数字滤波,指的是输入、输出均为数字信号,通过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分。数字滤波可分为无限脉冲响应(IIR)滤波和有限脉冲响应(FIR)滤波。本次试验利用VISUALDSP++软环境SIMULATOR模拟实现无限脉冲响应(IIR)数字信号处理。无限脉冲响应(IIR)的系统函数为:公式1即如果输入为X(Z),输出为Y(N),则:Y(Z)=X(Z)×H(Z),即公式2本试验中利用的公式是对上面的公式2对了相应的形式变化,利用中间变量公式3那
2、么公式4因此,本次试验设计了一个4阶IIR滤波,其对应的公式3和公式4的时域公式如下:;公式5;公式6在本次设计中由于系数a4,a3,a2,a1,b4,b3,b2,b1都是用户自己初始化的时候给定的且均为常数,所以公式5和公式6也相应可以变为公式7和公式8:公式7公式8其中x(n)是输入的数字序列。从公式5和公式6可知:只要我们设定设计需要的a4,a3,a2,a1,b4,b3,b2,b1,以及初始化w(n-1),w(n-2),w(n-3),w(n-4),就可以得到我们所需要的滤波器了。通过本练习使学生加深对IIR的理
3、解,加深对数字信号处理理论的理解。输入数据放在直接在数组inputs[N]中给出,滤波器系数a4,a3,a2,a1,b4,b3,b2,b1按顺序a2,a4,b2,b4,a1,a3,b1,b3存放在系数缓冲区coeffs[2*SECTIONS]中,w(n-1),w(n-2),w(n-3),w(n-4)初始化存放在delayline[SECTIONS](此处滤波器设计为因果滤波器,n<0时w(n)全为0),程序编译运行时对其进行IIR运算,计算结果存放于输出缓冲区output[N]中。参考程序如下:***********
4、*************************************************************/#defineN_MAX100#defineN30//输入数据的个数N#defineSECTIONS4#definescale0.05078125滤波器.K的K=scale。/************************************************************************/.sectiondata1;.align4;.varinputs[N]=50
5、00.0,4333.0,5465.0,13556.0,7423.0,-5000.0,-4333.0,-5465.0,-13556.0,-7423.0,5000.0,4333.0,5465.0,13556.0,7423.0,-5000.0,-4333.0,-5465.0,-13556.0,-7423.0,5000.0,4333.0,5465.0,13556.0,7423.0,-5000.0,-4333.0,-5465.0,-13556.0,-7423.0;.varcycle_count;//executioncycle
6、counts.align4;.varexpected_output[N]=253.906250,1394.910156,3400.899369,5451.189211,7461.879024,8509.488592,5519.537707,-1306.039745,-6936.575354,-8920.200002,-8397.397895,-4821.647412,1507.253965,6689.519912,8738.028980,8444.656939,4922.543154,-1492.260614,-67
7、30.424864,-8759.801131,-8433.783468,-4908.582993,1492.343760,6724.040670,8757.440081,8435.857839,4910.422599,-1492.627568,-6724.989198,-8757.654792;.align4;.varoutput[N];.sectiondata2;.align4;//a4,a3,a2,a1,b4,b3,b2,b1在varcoeffs[2*SECTIONS]存放的顺序为a2,a4,b2,b4,a1,a
8、3,b1,b3.varcoeffs[2*SECTIONS]=0.1412,0.6272,4.,4.,-0.0255,-0.6108,1.,6.;.vardelayline[SECTIONS]=0.0,0.0,0.0,0.0;//W的初始化值/****************************************************