DSP_FIR实验报告.docxVIP

实验三 FIR滤波器的DSP实验1. 启动CCS，打开实验工程文件，再编译并装载程序；2. 仿真波形输入与输出波形比较输入信号频谱图输出信号频谱图3.修改代码：滤波器设计要求：线性相位低通滤波器，模拟信号通带为0~1kHz，阻带衰减在5kHz处为-60dB。采样率为20kHz。采用凯塞窗口设计此低通滤波器，设计过程如下：通带边界频率为1KHZ，阻带边界频率是5KHZ有60dB的衰减，可得：过渡带，A=60dB设计凯塞窗口的长度N和参数分别为：使用MATLAB设计此滤波器，过程如下：1）在Matlab的Start 菜单中选择Toolboxes - Filter Design - Filter Design Analysis Tools(fdatool)，或者在命令行中输入fdatool 来启动滤波器设计分析器。2）选择设计fir低通滤波器，阶数为29，采样频率fs=20000Hz，截止频率fc=3000Hz，窗口设置为凯塞窗口，凯塞窗口参数为5.653263）用SPtool 验证所设计滤波器的功能在FDATool 里设计了一个FIR 低通滤波器后，用SPTool 来分析所设计的滤波器.在MATLAB 环境下产生一个混叠波信号，在MATLAB 命令窗口输入：Fs=20000;t=(0:200)/Fs;s=sin(2*pi*t*1000)+sin(2*pi*t*5000);subplot(1,1,1);plot(t,s);axis([0 0.005 -2 2]);grid产生了1000Hz 和5000Hz 的混叠波，并且变量Fs、t、s 将显示在workspace列表中。4)在FDATool 的界面里选择File→Export to →SPTool，然后在将导入SPTool的滤波器命名为filter点击Apply,这样,设计的滤波器就导入了SPTool 中。在弹出的SPTool 主界面里点击菜单File→Import 将信号s 导入,采样频率Fs导入。产生输入输出频谱后，得到5）可以观察到输入信号的波形和频谱为6）可以观察到输出信号的波形和频谱为可以看到，5000HZ频率的信号已被滤波器滤掉。7）将此滤波器的参数输出到C语言头文件中，修改CCS中fir滤波器的代码为float fir[30] ={0.0004001965019842,0.0001796782579809,-0.001610915830384,-0.003928668671765,-0.002911158381528, 0.004441058624642, 0.01416335142997, 0.01444447680484, -0.004466774967774,-0.03524603325736, -0.04900562431329, -0.01223772579527, 0.0834377913716, 0.2040327190773, 0.2883076291491, 0.2883076291491, 0.2040327190773, 0.0834377913716, -0.01223772579527, -0.04900562431329, -0.03524603325736, -0.004466774967774, 0.01444447680484, 0.01416335142997, 0.004441058624642,-0.002911158381528, -0.003928668671765,-0.001610915830384, 0.0001796782579809,0.0004001965019842};8)得到输入波形与频谱如下：输入与输出波形比较输入信号频谱图输出信号频谱图4.通过代码修改、仿真，理解定点运算和浮点运算相比的优缺点。定点与浮点，两者相比较而言，定点DSP处理器速度快，功耗低，价格也便宜；而浮点DSP则计算精度高，动态范围大。5.通过仿真波形的分析、理解滤波器阶数与滤波器延时以及滤波器的阻带衰减之间的关系。 在利用窗函数设计的FIR滤波器中，用同一种窗函数设计的FIR数字滤波器的阶数N不同，所得到的滤波器的过渡带会不一样，N越大，过渡带越窄。滤波器的阻带衰减只有窗形状决定，不受滤波器阶数的影响；但是滤波器的阶数太高，会导致过渡带很窄，暂态增加，延时很长。