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

1、DSP课程设计报告题目：FIR高通滤波器设计题目：FIR高通滤波器设计姓 名学 号教学院系专业年级指导教师目录 TOC o 1-5 h z HYPERLINK l bookmark2 o Current Document 一、设计题目1 HYPERLINK l bookmark4 o Current Document 二、设计目标1 HYPERLINK l bookmark6 o Current Document 三、算法研究与参数计算 1 HYPERLINK l bookmark8 o Current Document 1、FIR 的原理和参数生成公式 12、利用MATLAB计算滤波系数1

2、HYPERLINK l bookmark10 o Current Document 3、输入信号参数计算 2 HYPERLINK l bookmark12 o Current Document 四、编写源程序 3五、调试过程41 、调试前准备 52、MATLAB 的使用5 HYPERLINK l bookmark16 o Current Document 3、编写及编译程序 5 HYPERLINK l bookmark26 o Current Document 4、设置断点和探针 6 HYPERLINK l bookmark28 o Current Document 5、打开观察窗口 6 HY

3、PERLINK l bookmark30 o Current Document 六、实验结果及分析 6 HYPERLINK l bookmark32 o Current Document 1、输入信号的时域波形和频域波形 6 HYPERLINK l bookmark34 o Current Document 2、输出信号的时域波形和频域波形 7 HYPERLINK l bookmark36 o Current Document 七、设计心得8一、设计题目FIR高通滤波器设计二、设计目标设计一个FIR高通滤波器，通带边界频率为6000Hz，采样频 率为20000Hz。FIR滤波器的设计用MATL

4、AB窗函数法进行。三、算法研究与参数计算1、FIR的原理和参数生成公式(1即阶有阻冲徹响应淖波器tFIR)舍戎；尿七)皿 +-舍知1N土(MHFHt设计原理，根据系数h雄偶对称为了简址运畀产生如卜计弊方法期果一个F氓滤液亩一牛冲溅响胸h(Q), h(n n(N-U和K(n)锚绘愉入的时常酒波m输出滤波了5)的n涪血立下方常武：y 5 J=hCO)x(nHh( l)x(n-l )+h(2 )x(n*2)+.+htW)xnn-l)公武 1 2i itaB图3-12、利用MATLAB计算滤波系数在MATLAB界面输入图3-2所示程序，可得到滤波系数并生成INC文件。y=rourLd(firl (6

5、6000/200002, ! high，)+3250): id=f open(J 01 26. inc? J 1/ );fprintf (fid, ? . -word 35dn, 3 yO ;fclose (fid)图3-2输入freqz (y, 1, 512), MATLAB中显示高通滤波器的滤波特 性曲线。如图3-3所示。图3-33、输入信号参数计算MATLAB中输入图3-4中所示程序，包含两种频率成分的正弦 信号，一种信号频率1000Hz，种信号6000Hz。图3-4 其频谱特性曲线如图3-5。0D.10.20.3 OJ 0.50.607030.9Noimailzed Freriuenc

6、y 谀;i 憎誉冏闻刃图3-5四、编写源程序参考资料，编写汇编语言源程序：HIGHPASS .set 1 ;if you want to use ,please set the value to 1.global start,fir .mmregsCOFF_FIR_START: .sect coff_fir.copy0126.incK_FIR_BFFR.set 64d_data_buffer .usect fir_bfr,64FIR_DP.usect fir_vars,0d_filin.usect fir_vars,1output.usect fir_vars,1input.usect fir

7、_vars,1d_filout.usect fir_vars,100hstacksize .set 256 stack .usect fir_vars,stacksize.asgAR4,FIR_DATA_P.asgAR6,INBUF_P.asgAR7,OUTBUF_P.asgAR3,OUTBUF.asgAR2,INBUF.sectfir_prognopstart:stm #stack+stacksize,SPLD #FIR_DP,DPSTM #d_data_buffer,FIR_DATA_PRPTZ A,#K_FIR_BFFR-1STL A,*FIR_DATA_P+STM #d_filin,I

8、NBUF_PSTM #d_filout,OUTBUF_PSTM #output,OUTBUFSTM #input,INBUFSTM #100h,BK fir_loop:NOP ;Add Breakpoint & porbe pointLD *INBUF_P,ASTL A,*INBUFCALL firSTHA,*OUTBUF_P+%STHA,*OUTBUFmain_end:b fir_loopfir:; SSBX SXM; SSBXSTMFRCTSTL; SSBXSTMFRCTSTLSTM#d_data_buffer,FIR_DATA_PA,*FIR_DATA_P#(d_data_buffer+

9、K_FIR_BFFR-1),FIR_DATA_Pfir_task:RPTZA,#K_FIR_BFFR-1RPTZMACD *FIR_DATA_P-,COFF_FIR_START,ARET.end五、调试过程1、调试前准备1）启动SETUP并选择芯片。2）在CCS的安装目录下的myproject的子目录下创建新文件夹, 命名 1207030126。2、MATLAB的使用使用MATLAB生成0126.Inc文件和1207030126.dat文件，保存 在1207030126文件夹中。如图5-1：Cnirent Directory d： .Pr ograni Fi 1 esNJLTLA371 vio

10、rk事 艾a 1 a -All Files VFile TypeSizeLast ModifielJ01207030126.datDAT File2 KB 2015-7-911:07O126.inclINC File1 KB 2015-7-9 11:07*1I2J图5-13、编写及编译程序1）新建 project 并命名 1207030126。2）在project中 添力口 0126.asm 文件和 0126.cmd 文件。如图 5-2。3）编译并链接程序。4 ）设置工程文件，选择菜单“ project”的 “build options 项。图5-25）加载程序。Fileload progra

11、m4、设置断点和探针在标号“fir_loop”下的“NOP”语句处设置软件断点和探针。 选择“file”、“file 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数字高通滤波器，由于 数字信号处理课程是上学期学的，再加上自己对MATLAB软件操作 不熟练。导致前期用了很长一段时间查阅相关书籍，重新学习了数字 滤波器的