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

DSP课程设计报告题目：FIR高通滤波器设计姓名学号教学院系专业年级指导教师目录一、设计题目1二、设计目的1三、算法研究与参数计算11、FIR的原理和参数生成公式12、运用MATLAB计算滤波系数13、输入信号参数计算2四、编写源程序3五、调试过程41、调试前准备52、MATLAB的使用53、编写及编译程序54、设立断点和探针65、打开观测窗口6六、实验结果及分析61、输入信号的时域波形和频域波形62、输出信号的时域波形和频域波形7七、设计心得8一、设计题目FIR高通滤波器设计二、设计目的设计一个FIR高通滤波器，通带边界频率为6000Hz，采样频率为20230Hz。FIR滤波器的设计用MATLAB窗函数法进行。三、算法研究与参数计算1、FIR的原理和参数生成公式图3-12、运用MATLAB计算滤波系数在MATLAB界面输入图3-2所示程序，可得到滤波系数并生成INC文献。图3-2输入freqz（y，1，512），MATLAB中显示高通滤波器的滤波特性曲线。如图3-3所示。图3-33、输入信号参数计算MATLAB中输入图3-4中所示程序，包含两种频率成分的正弦信号，一种信号频率1000Hz，一种信号6000Hz。图3-4其频谱特性曲线如图3-5。图3-5四、编写源程序参考资料，编写汇编语言源程序：HIGHPASS.set1;ifyouwanttouse,pleasesetthevalueto1.globalstart,fir.mmregsCOFF_FIR_START:.sect"coff_fir".copy"0126.inc"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.asgAR4,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;AddBreakpoint&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）启动SETUP并选择芯片。2）在CCS的安装目录下的myproject的子目录下创建新文献夹，命名。2、MATLAB的使用使用MATLAB生成0126.Inc文献和.dat文献，保存在文献夹中。如图5-1：图5-13、编写及编译程序1）新建project并命名。2）在project中添加0126.asm文献和0126.cmd文献。如图5-2。3）编译并链接程序。4）设立工程文献，选择菜单“project”的“buildoptions…”项。图5-25）加载程序。File—loadprogram4、设立断点和探针在标号“fir_loop”下的“NOP”语句处设立软件断点和探针。选择“file”、“fileI/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从输入信号波形和输出