DSP实验报告-FIR和IIR滤波器设计

....DSP试验一FIR试验目的TMS320VC5416DSPFIRCCSTMS320VC5416试验要求弦波合成文件进展滤波测试，显示出输入和输出信号。试验原理ej,d H(ej)d 0,c

c样响应为h(n)N。h(n)表示截取h(n后的冲激响应，即h(n)(n)h(n)，d d其中(n)即为窗函数，窗长为N。一般的FIRy(n)n1h(k)x(nk)k0ZFIRN1 H(z) h(n)zN1 NN+1FIR11FIR设计思路FIRh(n)matlabmatlabmatlabdatCFIR试验内容滤波器系数的设计matlabfir1fir1的调用格式为〔，〕，6dBN〔单〔n〕=N+1〕的FIR为标量用哈明窗。其中为对归一化的数字频率，其计算公式为，这里计算得出DSPFIR2210FIR对于101134。310FIR

410FIRmatlabwavread去，并编写程序查看其局部信号时域和频域图，并生成两个语音文件的.dat文DSP5〔a〔b〕6〔a〔b。1语音文件1频谱图图5语音文件1〔a〕2

26语音文件2CCSsetup5416CCSproject->new建立fir滤波工程。建一个input.h文件，并将matlab产生的input文件里面的数据放入其中，然后project->addfilestoproject..input.hFIRfird77CCS3.3fir.outfile->loadprogram..->fir.outstartaddressinput和output对应输入和输出的数据。通过转变acquisitionbuffersize和9-12示。8〔a〕时域图〔b〕9输入信号1滤波结果时域图滤波结果频域图10信号1〔a〕〔b〕频域图11输入信号2〔a〕滤波结果时域图

〔b〕滤波结果频域图12信号2configurationDSK5416.cdb，保存至文件的工程名目下。将刚生成的.cdb文件和编写的参与到工程中。将前面的程序文件分别参与工程中，进展编译。链接试验板和电脑。点击debug->connect，连接测试板，file->loadprogram..，从工程文件下的debug文件夹中导入fir.outdebug->run试验结果11Khz、2Khz、3Khz22.5Khz的信号根本滤去。由于滤波器的阶数较小，对于3Khz信号并没有完全滤除，这也是试验预期之中的。对于信号2的复合信号，由频谱图比较可以觉察，大于FIR率为2.5Khz的低通滤波器的功能。试验二IIR试验目的和要求在上面设计了FIR滤波器的根底上，设计一样性能的IIR滤波器。试验原理IIR现，但是IIRIIRFIR馈环节进展了扩展，增加了一个反响环节。IIR〔L-1〕阶前馈环节和〔M-1〕阶反响环节的差分方程表示为：ynL1bxni1ai mi0 m1z进而得到IIR滤波器的传递方程为：

IIRab。试验内容FIRIIRFIRIIRIIRmatlabIIR121314滤波后时域图滤波后频域图13IIR1〔a〕滤波后时域图试验结果

〔b〕滤波后频域图14IIR2IIRFIR2.5KhzIIR2Khz3Khz22.5Khz试验总结10FIRIIRFIRIIRFIR和IIR滤波器各种的优了简洁的滤波器实现算法和程序。由于之前没有较多接触过CCSDSPCCScmd过与同学一起争论和在网上查询才知道解决方法。DSPFIRIIR附录输入信号生成程序〔matlab〕clc;clearall;closeall;FS=8000;%[y,FS,NBITS]=wavread(”C:\Users\Administrator\Desktop\MultiToneTest.wav”);[y,FS,NBITS]=wavread(”C:\Users\Administrator\Desktop\sc03.wav”);figure(1)plot(y(1:500));title(”语言信号时域图”);figure(2)y=abs(fft(y(1:1000)));df=(1:1000)*(FS/1000);plot(df,y(1:1000));title(”语言信号频域图”);y=(y*32767);y=int16(y);fid=fopen(”input.dat”,”w”);form=1:2048fprintf(fid,”%d,\n”,y(m));%输出endfclose(fid);FIR#include“stdio.h“#defineN11 //FIR滤波器的级数+1(由于本试验中滤波器阶数为10)#defineLEN2048 //待滤波的输入数据长度#include“input.h“//输入数据文件longyn;intB[11]={-64,438,-530,-2519,8803,20512,8803,-2519,-530,438,-64};//滤波器系数intinput[LEN]; //输入缓冲，在仿真时将从内存载入intoutput[LEN]; //输出缓冲，直接存放在内存中voidmain{inti,j;int*x;for(j=0;j<LEN-1;j++){x=&input[j];yn=0;for(i=0;i<N-1;i++){yn+=B[i]*(*x++);%滤波}output[j]=yn>>15;}while(1);}IR#include“stdio.h“#include“input.h“#defineN11#defineLEN2048intNUM[10]={69,691,3108,8289,14505,17406,14505,8289,3108,691,69};intDEN[10]={4096,10203,16248,16727,12645,6959,2837,828,166,20,1};longyn,yn1,yn2;intinput[LEN];intinputY[LEN];intoutput[LEN];voidmain{inti,j;int*x,*y;for