DSP实验报告解析.docVIP

DSP实验报告  软件实验 １无限冲激响应滤波器(IIR) 算法 一．实验目的 1 ．掌握设计 IIR 数字滤波器的原理和方法。 2 ．熟悉 IIR 数字滤波器特性。 3 ．了解 IIR 数字滤波器的设计方法。 二．实验设备 PC 兼容机一台，操作系统为 Windows2000( 或 Windows98 ，WindowsXP ，以下默认为Windows2000) ，安装 Code Composer Studio 2.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.7757y[n-1]+0.1122x[n]+0.1122x[n-1] 。 ５．程序流程图： 四．实验步骤 1 ．实验准备： - 设置软件仿真模式 - 启动 CCS 2 ．打开工程，浏览程序，工程目录为 C:\ICETEK\VC5416AES61\VC5416AES61\Lab0502-IIR\IIR.pjt 。 3 ．编译并下载程序。 4 ．打开观察窗口： * 选择菜单 View-Graph-Time/Frequency… 进行如下图所示设置。 * 选择菜单 View-Graph-Time/Frequency… ，进行如下设置： 5 ．清除显示：在以上打开的窗口中单击鼠标右键，选择弹出式菜单中“ CleAe Display ”功能。 6 ．设置断点：在程序 iir.c 中有注释“ break point ”的语句上设置软件断点。 7 ．运行并观察结果： ⑴ 选择“ Debug ”菜单的“ Animate ”项，或按 F12 键运行程序。 ⑵ 观察“ IIR ”窗口中时域图形；观察滤波效果。 8 ．退出 CCS 五．实验结果 输入波形为一个低频率的正弦波与一个高频的余弦波叠加而成。如图： 通过观察频域和时域图， 得知： 输入波形中的低频波形通过了滤波器， 而高频部分则被衰减。 实验代码： #include math.h #define IIRNUMBER 2 #define SIGNAL1F 1000 #define SIGNAL2F 4500 #define SAMPLEF 10000 #define PI 3.1415926 float InputWave(); float IIR(); float fBn[IIRNUMBER]={ 0.0,0.7757 }; float fAn[IIRNUMBER]={ 0.1122,0.1122 }; float fXn[IIRNUMBER]={ 0.0 }; float fYn[IIRNUMBER]={ 0.0 }; float fInput,fOutput; float fSignal1,fSignal2; float fStepSignal1,fStepSignal2; float f2PI; in