DSP原理及应用大作业.doc

DSP原理及应用大作业 题 目： IIR算法的软件实验 学生姓名： 学 号： 院（系）： 信息工程学院 专 业： 电子信息工程专业 IIR算法实验 摘要：本课题通过软件设计IIR数字滤波器，并对所设计的滤波器进行仿真：应用DSP集成开发环境—CCS调试程序，用TMS320F2812实现IIR数字滤波。具体工作包括：对IIR数字滤波器的基本理论进行分析和探讨。应用DSP集成开发环境调试程序，用TMS320F2812来实现IIR数字滤波。 关键词：数字滤波器；DSP；TMS320F2812；无限冲激响应滤波器（IIR）。 引言：21世纪是数字化的时代，随着信息处理技术的飞速发展，数字信号处理技术逐渐发展成为一门主流技术。相对于模拟滤波器，数字滤波器没有漂移，能够处理低频信号，频率特性可做成非常接近于理想的特性，且精度可以达到很高，容易集成等。这些优势决定数字滤波器的应用越来越广泛。数字滤波器是数字信号处理中最重要的组成部分之一，被广泛应用于语音图像处理、数字通信、谱分析、模式识别、自动控制等领域。 实验原理 1．无限冲激响应数字滤波器的基础理论。 2．模拟滤波器原理（巴特沃斯滤波器、切比雪夫滤波器、椭圆滤波器、贝塞尔滤波器）。 3．数字滤波器系数的确定方法。 4．根据要求设计低通IIR滤波器 要求：低通巴特沃斯滤波器在其通带边缘1kHz处的增益为-3dB，12kHz处的阻带衰减为30dB，采样频率25kHz。设计： -确定待求通带边缘频率fp1Hz、待求阻带边缘频率fs1Hz和待求阻带衰减-20logδsdB。 模拟边缘频率为：fp1=1000Hz，fs1=12000Hz 阻带边缘衰减为：-20logδs=30dB -用Ω=2πf/fs把由Hz表示的待求边缘频率转换成弧度表示的数字频率，得到Ωp1和Ωs1。 Ωp1=2πfp1/fs=2π1000/25000=0.08π弧度 Ωs1=2πfs1/fs=2π12000/25000=0.96π弧度 -计算预扭曲模拟频率以避免双线性变换带来的失真。 由w=2fs tan(Ω/2)求得wp1和ws1，单位为弧度/秒。 wp1=2fs tan(Ωp1/2)=6316.5弧度/秒 ws1=2fs tan(Ωs1/2)=794727.2弧度/秒 -由已给定的阻带衰减-20logδs确定阻带边缘增益δs。 因为-20logδs=30，所以logδs=-30/20，δs=0.03162 -计算所需滤波器的阶数： 因此，一阶巴特沃斯滤波器就足以满足要求。 -一阶模拟巴特沃斯滤波器的传输函数为： H(s)=wp1/(s+wp1)=6316.5/(s+6316.5) 由双线性变换定义s=2fs(z-1)/(z+1)得到数字滤波器的传输函数为： 因此，差分方程为：y[n]=0.3307y[n-1]+0.3346x[n]+0.3346x[n-1] 软件程序流程图 开始 开始 初始化工作变量 调用波形发生子程序产生混叠的波形(高频+低频) 调用IIR滤波子程序计算当前输出 波形发生 计算步长 用标准C的sin函数和cos函数计算当前波形值 返回波形值 IIR滤波 用滤波器系数乘以保存的N-1个输入输出值和当前输入值并求和 返回计算结果 调试过程与步骤： 1．实验准备 -设置软件仿真模式。 -启动CCS 打开工程，浏览程序，工程目录为D:\dsp\t6\iir\iir.pjt 3．编译并下载程序 4．打开观察窗口：*选择菜单View-Graph-Time/Frequency…，进行如下设置： *选择菜单View-Graph-Time/Frequency…，进行如下设置： 5．清除显示：在以上打开的窗口中单击鼠标右键，选择弹出式菜单中“Clear Display”功能。 6．设置断点：在程序iir.c中有注释“/* 请在此句上设置软件断点 */”的语句上 置软件断点。 7.运行并观察结果 ⑴选择“Debug”菜单的“RUN”项，或按F5键运行程序。 ⑵观察“IIR”窗口中时域图形；观察滤波效果。 实验结果 输入波形为一个低频率的正弦波与一个高频的余弦波叠加而成。如图： 通过观察频域和时域图，得知：输入波形中的低频波形通过了滤波器，而高频部分则被衰减。 IIR算法相对应C语言编程及汇编程序 /*=========================================== =====================================