《无限冲激响应滤波器（IIR）算法硬件实现DSP论文》由会员分享，可在线阅读，更多相关《无限冲激响应滤波器（IIR）算法硬件实现DSP论文（22页珍藏版）》请在金锄头文库上搜索。
1、无限冲激响应滤波器（IIR）算法硬件实现摘要20世纪60年代以来，随着大规模集成电路、数字计算机等信息技术的飞速发展，数字信号处理（Digital Signal Processing, DSP）技术应运而生并得到快速发展。在过去的20多年时间里，DSP在理论和应用方面不断地进步和完善，在越来越多的应用领域中迅速取代传统的模拟信号处理方法，并开辟出出许多新的应用领域。目前数字信号处理技术已经在通信、雷达、航空航天、工业控制、生物医学控制、生物医学工程、网络及家电领域得到极为广泛的应用，数字化时代正在到来。1928年，美国德州仪器公司（Texas Instruments Incorporation
2、,简称TI公司）推出该公司的第一款DSPs芯片，很快DSPs芯片就以其数字器件特有的稳定性、可重复性、可大规模集成和易于实现DSP算法等优点，为数字信号处理技术带来了更大的发展和应用前景。采用各种类型DSPs实现系统的数字化处理和控制已经成为未来发展的趋势，并且睡着DSPs运算能力的不断提高，数字信号处理的研究重点由最初的非实用转向高速实时应用。本文主要是应用TI公司生产的一款芯片TMS320F2812 DSPs来实现无限冲击响应滤波器（IIR）算法。本文的主要工作可归结如下：1、 介绍无限冲激响应滤波器（IIR）实现的原理。2、 详细分析无限冲激响应滤波器（IIR）算法。3、 画出无限冲激响
3、应滤波器（IIR）实现的硬件框图及软件流程图。4、 实验调试过程及步骤。5、 分析在CCS软仿真及硬仿真环境下实验结果。6、 反汇编一个小程序及对DSPs实现IIR算法理解。7、 对学习DSPs理论课及做实验的心得体会。通过软件仿真及硬件实现结果表明，用DSPs实现无限冲激响应滤波器（IIR）算法能够很好的达到预期的效果。关键词：IIR, , DSPs, 无限冲激响应，CCS，滤波器引言硬件实现无限冲激响应滤波器（IIR）算法有多种方法,本文介绍用TI公司生产的一款芯片TMS320F2812 DSPs来实现无限冲激响应滤波器（IIR）算法，其中涉及无限冲激响应滤波器（IIR）算法软件仿真和硬件
4、实现。通过实验分析，我们可以看出实验结果和理论结果符合的很好，因此MS320F2812 DSPs可以很好的实现无限冲激响应滤波器的滤波特性。实验原理1、 无限冲激响应数字滤波器的理论基础。2、 模拟滤波器原理（巴特沃斯滤波器、切比雪夫滤波器、椭圆滤波器、贝赛尔滤波器）。3、 数字滤波器系数的确定方法。4、 根据要求设计低通IIR滤波器要求：低通巴特沃斯滤波器在其通带边缘1kHz处的增益为-3dB,12kHz处的阻带衰减为30的dB,采样频率25kHz。设计：-确定待求通带边缘频率fp1Hz、待求阻带边缘频率fs1Hz和待求阻带衰减-20logs dB。模拟边缘频率为：fp1=1000 kHz，
5、fs1=12000kHz，阻带边缘衰减为：-20logs dB=30dB。-用=2f/fs把由Hz表示的待求边缘频率转换成弧度表示的数字频率，得到p1，s1p1=2fp1/fs1 =21000/25000=0.08弧度s1=2fs1/fs =212000/25000=0.96弧度 -计算预扭曲线模拟频率以避免双线性变换带来的失真。由w=2fstan(/2)求得wp1和ws1,单位为弧度/秒wp1=2fstan(p1/2)=6316.5弧度/秒ws1=2fstan(s1/2)=794727.2弧度/秒 -由已给定的阻带衰减-20logs dB确定阻带边缘增益s。 因为-20logs dB=30d
6、B，所以logs=-30/20，s=0.03162-计算所需滤波器阶数：因此，一阶巴特沃斯滤波器就满足要求。-一阶巴特沃斯滤波器的传输函数为：H(s) = wp1/(s+wp1)=6316.5/(s+6316.5)有双线性变换定义s=2fs(z-1)/(z+1)得到的数字滤波器的传输函数为：因此，差分方程为：yn=0.7757yn-1+0.1122xn+0.1122xn-1.硬件框图因此，差分方程为：yn=0.7757yn-1+0.1122xn+0.1122xn-1.硬件框图D/A4路12bit数据地址控制I/OPWMSPISPI4个用户可控开关LEDCPLD扩展DSP引脚RAM:128K*1
7、6bit RAM:128K*16bitDSP总线RS232CAN驱动16路，12bitJATGTMS320F2812PGFAJTAGPORT 片上资源 主频：150MHzA/D RAM:128K*16bit FLASH:128K*16bit CANA/DF2812-A评估板原理图液晶显示及键盘模块信号源模块DSP评估板模块仿真器模块PC机实验中用到的主要硬件模块软件流程图开始初始化工作变量调用波形发生子程序产生混叠的波形（高频+低频）调用IIR滤波器子程序计算当前输出波形发生计算步长用标准C的函数sin函数和cos函数计算当前波形值返回波形值返回计算结果IIR波形用滤波器系数乘以保存的N-1个
8、输入输出值和当前输入值并求和返回计算结果调试过程和步骤1. 实验准备（1） 连接实验设备。（2） 准备信号源进行AD输入。1) 取出2根实验箱附带的信号线。2) 用1根信号线连接实验箱左侧信号源的波形输出A端口和“A/D输入”模块的“ADCIN0”插座，注意插头要插牢、到底。这样，信号源波形输出A的输出波形即可送到 ICETEK-F2812A板的AD输入通道0.3) 用1根信号线连接实验箱左侧信号源的波形输出B端口和“A/D输入”模块的“ADCIN1”插座，注意插头要插牢、到底。这样，信号源波形输出A的输出波形即可送到 ICETEK-F2812A板的AD输入通道1.4) 设置波形输出A:-向内
9、侧按波形频率选择按钮，直到标有正弦波的指示灯亮。-上下调节波形频率选择按钮，直到标有100-1KHz的指示灯亮。-调节幅值调整旋钮，将波形输出A的幅值调到适当位置。 5）设置波形输出B: -向内侧按波形频率选择按钮，直到标有正弦波的指示灯亮。-上下调节波形频率选择按钮，直到标有1K-10KHz的指示灯亮。-调节幅值调整旋钮，将波形输出B的幅值调到适当位置。注意：由于模数输入信号未经任何转换就进入DSP,所以必须保证输入的模拟信号的幅度在03V之间。必须用示波器检测信号范围，保证最小值0V,最大值3V,否则容易损坏DSP芯片的模数采集模块。2. 设置 Code Composer Studio 2
10、.21在硬仿真（Emulator）方式下运行。3. 启动Code Composer Studio 2.21 选择菜单Debug-Reset CPU。4. 打开工程文件 工程目录：D:dspt8mixeiirmixeriir.pjt5. 编译、下载程序，选择菜单Debug-Go Main,使程序运行到 main 函数的入口位置。6. 观察窗口 -打开IIR.c,查看源代码。7. 观察程序运行结果 按CTR控制板的K6键，实现滤波显示，K7键实现混频显示，按K8键实现A、B两信号源分屏显示。8. 观察动态效果，调节信号源输出，观察滤波器输出，改变信号源输入的波形、频率参数，观察动态效果。9. 退出
11、CCS。实验结果一、软仿真环境下的实验结果1. 作如下设置及得到的波形图图12作如下设置得到的波形图图23. 作如下设置得到的波形图图34. 作如下设置得到的波形图图44. 把四幅图片放着一起比较如下图5二、硬仿真环境下的实验结果1. 按CTR控制板的K6键观察到的波形图62. 按K7键观察到的波形图73. 按K8键观察到的波形图8实验结果分析1. 软仿真结果分析-图1是信号经过调制后的模型，毛刺代表其高频分量。-图3是图1中信号的频域图形，从频域可以看出，其频谱既含有低频分量，又含有高频分量。-图2调制信号经过无限冲击响应滤波器（IIR）滤波之后得到的波形。-图4图2中信号的频域图形，从频域
12、可以看出，经过无限冲击响应滤波器（IIR）滤波之后只剩下低频分量，高频分量几乎衰减为0。-图5把四幅图放在一块比较，可以显而易见的看出图形滤波前后的差别，经比较可以看出，在软仿真条件下，用TMS320F2812 DSPs来实现的无限冲击响应滤波器（IIR）具有较好的滤波性能。2.硬仿真结果分析 -图 6 上面为低频信号，下面为高频信号。 -图7为两信号相乘的结果，即调制的结果。 -图8为调制信号经无限冲击响应滤波器（IIR）的结果。比较图6和图8我们可以发现，两图的波形几乎完全相同，由此我们可以得出如下结论： 用TMS320F2812 DSPs来实现的无限冲击响应滤波器（IIR）具有较好的滤波
13、性能。程序算法分析及反汇编1、根据自己的理解，给iir.c程序添加备注#include DSP281x_Device.h / DSP281x Headerfile Include File#include DSP281x_Examples.h / DSP281x Examples Include File#include f2812a.h#includemath.h#define IIRNUMBER 2#define SIGNAL1F 1000#define SIGNAL2F 4500#define SAMPLEF 10000#define PI 3.1415926float InputWave();float IIR();float fBnIIRNUMBER= 0,0.7757 ; /定义差分方程y（n)系数,并赋值float fAnIIRNUMBER= 0.1122,0.1122 ; /定义差分方程x(n)的系数，并赋值float fXnIIRNUMBER= 0.0 ; /定义输入序列float fYnIIRNUMBER= 0.0 ; /定义输出序列float fInput,fOutput;
