DSP高通滤波器课程设计报告

1、而船石油火毒DSP课程设计报告题目：FIR高通滤波器设计姓名学号教学院系专业年级指导教师目录一、设计题目1二、设计目标1三、算法研究与参数计算11、FIR的原理和参数生成公式12、利用MATLA时算滤波系数13、输入信号参数计算2四、编写源程序3五、调试过程41、调试前准备42、MATLAB勺使用53、编写及编译程序54、设置断点和探针55、打开观察窗口6六、实验结果及分析61、输入信号的时域波形和频域波形62、输出信号的时域波形和频域波形6七、设计心得7一、设计题目FIR高通滤波器设计二、设计目标设计一个FIR高通滤波器，通带边界频率为6000Hz,采样频率为20000HzFIR滤波器的设计

2、用MATLA函函数法进行。三、算法研究与参数计算1、FIR的原理和参数生成公式ojn阶有限冲激啊而通波器cm)公式；/(，)=(无*)+-(H-1+灯)公式】TNR12.FIR仪计原理；根据索嬖h是偶对称为T褐化运算产生如F计算方法如果一个网代滤波有一个冲激唁皮.h,h(1)和苒描说地人的时需滤波m辎出渔波式n)的n给出以卜.方程式；y(n)-h(O)x(n)+h(l)x(n-lhh(2)x(a-2>.+h(N1闻Xm切公式17图3-12、利用MATLA升算滤波系数在MATLABF面输入图3-2所示程序，可得到滤波系数并生成INC文件。y=roond(firl(64j6000/20000

3、*2,!high1八3#焚)；,Eid=fopen(T0126.inc1jTwr，);fprintf(fid,',-wordayJ;fclose(fid)图3-2输入freqz(y,1,512),MATLA中显示局通滤波器的滤波特性曲线。如图3-3所示图3-33、输入信号参数计算MATLA计输入图3-4中所示程序，包含两种频率成分的正弦信号，一种信号频率1000Hz,一种信号6000H4CoiriVTi&rLdlWizidow»i=D：l:255;y=romd(sin(2pi*i*1000/20000+sin(Z*pi*i*6000/20000)*32768/2);f

4、idfopenC1207030126.daf/v)；fprintf(fid/1651工DI0n);Eptirrtf(£l4jy);fclose(fid)图3-4其频谱特性曲线如图3-5。五三名mu>3一落*ufc附而丫H2印系> >lr«杳：口I图3-5四、编写源程序参考资料，编写汇编语言源程序：HIGHPASS.set1;ifyouwanttouse,pleasesetthevalueto1.globalstart,fir.mmregsCOFF_FIR_START:.sect"coff_fir".copy"0126.inc&

5、quot;K_FIR_BFFR.set64d_data_buffer.usect"fir_bfr",64FIR_DP.usect"fir_vars",0d_filin.usect"fir_vars",1output.usect"fir_vars",1input.usect"fir_vars",1d_filout.usect"fir_vars",100hstacksize.set256stack.usect"fir_vars",stacksize.asgAR

6、4,FIR_DATA_P.asgAR6,INBUF_P.asgAR7,OUTBUF_P.asgAR3,OUTBUF.asgAR2,INBUF.sect"fir_prog"nopstart:stm#stack+stacksize,SPLD#FIR_DP,DPSTM#d_data_buffer,FIR_DATA_PRPTZA,#K_FIR_BFFR-1STLA,*FIR_DATA_P+STM#d_filin,INBUF_PSTM#d_filout,OUTBUF_PSTM#output,OUTBUFSTM#input,INBUFSTM#100h,BKfir_loop:NOP;Ad

7、dBreakpoint&porbepointLD*INBUF_P,ASTLA,*INBUFCALLfirSTHA,*OUTBUF_P+%STHA,*OUTBUFmain_end:bfir_loopfir:;SSBXSXM;SSBXFRCTSTM#d_data_buffer,FIR_DATA_PSTLA,*FIR_DATA_PSTM#(d_data_buffer+K_FIR_BFFR-1),FIR_DATA_Pfir_task:RPTZA,#K_FIR_BFFR-1MACD*FIR_DATA_P-,COFF_FIR_START,ARET.end五、调试过程1、调试前准备1）启动SETU所

8、选择芯片。2）在CCS勺安装目录下的myproject的子目录下创建新文件夹，命名120703012&2、MATLAB勺使用使用MATLABk成0126.Inc文件和1207030126.dat文件，保存在1207030126文件夹中。如图5-1:irectorvLV-Pron-amFilesKMATLrBTlVworl:囱于地I国0,AlLFil电工L|Fi"Typ.SimL”ModiMalSl207030126.datDATFile2KB201S7-91107Ol26jncjINCFile1KB20157-911:072J2J_/_图5-13、编写及编译程序1）新建pro

9、ject并命名120703012&堂 Pj1«CEL £门”-Projects日图 1207Q012C,vjt国 oia »d：I nsMioe j皿 Lj Gf ma叵 |Ifl 0L26, inrI Libruri«sE L_ Source2）在project中添加0126.asm文件和0126.cmd文件。如图5-2。3）编译并链接程序。4）设置工程文件，选择菜单“project”的“buildoptions”项。图5-25）加载程序。File-loadprogram4、设置断点和探针在标号"fir_loop”下的“NOP语句处设

10、置软件断点和探针选择“filefile I/O;按照实验程序操作5、打开观察窗口选择菜单“view”、“Graph"、“time/Freauency”。按要求设置。六、实验结果及分析1、输入信号的时域波形和频域波形输入信号的时域波形如图6-1所示，频域波形如图6-2所示图6-1图6-2这个信号由1000Hz和6000Hz的信号叠加而得2、输出信号的时域波形和频域波形输出信号的时域波形如图6-3所示，频域波形如图6-4所示图6-3图6-4从输入信号波形和输出信号波形对比看，低频信号（1000Hz）不能顺利通过滤波器，被削减到几乎为0,而高频信号（6000HZ）可以顺利通过。因此这是一个性能良好的高通滤波器。七、设计心得本次我的课程设计是基于窗函数法的FIR数字高通滤波器，由于数字信号处理课程是上学期学的，再加上自己对MATLA歆件操作不熟练。导致前期用了很长一段时间查阅相关书籍，重新学习了数字滤波器的相关知识和MATLA歆件