实验8FIR数字滤波器的设计实验.docVIP

实验8 FIR数字滤波器设计实验 一．实验目的 1．掌握用窗函数法设计FIR数字滤波器的原理和方法。 2．熟悉线性相位FIR数字滤波器特性。 3．了解各种窗函数对滤波器特性的影响。 二．实验设备 PC兼容机一台，操作系统为WindowsXP，安装CCS软件。 三．实验原理 1．基于窗函数法设计FIR数字滤波器的基础理论。（请参考教材《数字信号处理教程》） 2．根据要求设计FIR数字低通滤波器。 要求：通带边缘频率，阻带边缘频率，阻带衰减70dB，采样频率。 (1) 求各对应的数字频率： 通带截止频率： 阻带截止频率： 阻带衰减相当于： (2) 求。设， 理想数字低通滤波器的截止频率为： 由此可得： 其中， (3) 求窗函数。由阻带衰减确定窗形状，由过渡带宽确定N。 由于，查表可选择布莱克曼窗。 由于布莱克曼窗过渡带宽满足： 故：窗函数长度 所以N=25，窗函数为： (4) 求。 (5) 根据上面计算出的，得到数字滤波器的差分方程为： 3．程序流程图： 四．实验步骤 1．设置软件仿真工作模式。（参看实验5的第四部分第1步）。 2．启动CCS软件。 3．打开工程文件，选择菜单Project→Open。 C:\ICETEK\VC5509Ae\VC5509Ae\Lab0501-FIR\Fir.pjt。 4．打开程序Fir.c，对照实验原理中差分方程的系数修改float fHn[.]中的数值。 5．编译并下载程序。 (1) 选择菜单Project→Rebuild All； (2) 执行File→Load Program，在随后打开的对话框中选择刚刚建立的 C:\ICETEK\VC5509AS60\Lab0501-FIR\Debug\Fir.out。 6．设置断点：在程序Fir.c中有注释“break point”的语句设置软件断点。 7．打开观察窗口： (1) 选择菜单View-Graph-Time/Frequency…，进行如图8-1设置： 图 8-1 (2) 选择菜单View-Graph-Time/Frequency…，进行如图8-2设置： 图 8-2 (3) 在弹出的图形窗口中单击鼠标右键，选择“Clear Display”。 8．运行并观察结果： (1) 选择“Debug”菜单的“Animate”项运行程序。 (2) 观察“Input”、“Output”窗口中时域图形；观察滤波效果。 9．退出CCS。 五．实验报告要求 1.?分别做出输入和输出信号的频率图。 2.?分别做出FIR数字滤波器的时域波形图和幅频特性图，并对照输入信号和输出时域和频率图解释观察到的现象。 3.?逐渐调整fStepSignal2的值，记录并分析实验结果； 4．根据实验原理中的设计要求，试用MATLAB编程实现FIR数字滤波器的设计，并完成与实验结果相同的滤波功能。 开始 初始化工作变量 调用波形发生子程序产生混叠的波形(高频+低频) 调用FIR滤波子程序计算当前输出 波形发生 计算步长 用标准C的sin函数和cos函数计算当前波形值 返回波形值 FIR滤波 用滤波器系数乘以保存的N-1个输入值和当前输入值并求和 返回计算结果