DSP课程设计---FIR滤波器设计.docxVIP

DSP课程设计题 目: FIR滤波器设计指导老师:姓 名: 学 号: 专业年级: 电 信 08院 系: 电气信息学院西南石油大学2012年01月12日目录一．设计题目：FIR滤波器设计1二．设计要求11.FIR滤波器设计要求12.利用MATLAB计算滤波器系数1三．功能描述11.部分语句说明12.指定存储器13.汇编程序部分说明：1四．算法特点11.FIR的原理及公式12.算法设计1五．相应参数计算11.高通滤波器的滤波系数的计算12.滤波器输入信号参数1六．编写程序11.源程序12..cmd程序13.滤波器系数测试程序14.滤波器输入信号生成程序1七．调试过程11.调试前的准备12.MATLAB的使用13.编写及编译程序14.观察点设置15.运行并观察结果1八．运行结果11.运行结果12.输入信号频谱图1九．设计心得1设计题目：FIR滤波器设计设计要求FIR滤波器设计要求设计一个FIR低通（或高通、带通）滤波器通带边界频率为1500Hz，通带波纹小于1dB；阻带边界频率为2000 Hz，阻带衰减大于40 dB；采样频率为8000 Hz.FIR滤波器的设计可以用MATLAB窗函数法进行。利用MATLAB计算滤波器系数fir(n,wn)例如：round(fir1(16,1500/8000*2)*32768) round(fir1(16,1500/8000*2,’high’)*32768)可用以下语句看频谱特性，如图4—1所示为低通滤波特性图4—1低通滤波特性曲线功能描述部分语句说明1）.global start,fir设定全局变量。2）COFF_FIR_START: .sect”coff_fir”。Include”ZB.inc”(设定系数文件)。3）K_FIR_BFFR .set32(滤波阶数)。4）d_filin（存放输入波形）。5）d_filout（存放输出波形）。指定存储器指定AR4为FIR_DATA_P数据寄存器。指定AR6为INBUF_P输入数据寄存器。指定AR7为OUTBUF_P输出数据寄存器。汇编程序部分说明：start部分：程序初始化部分指定寄存器，清空寄存器。fir_loop部分：循环调入输入数据，并调用子程序fir进行计算。main_end部分：跳转至循环部分。fir部分：子程序部分。其中fir_task部分：进行计算并返回计算结果。算法特点FIR的原理及公式滤波器公式：N=1,2,3…FIR设计原理：根据系数h是偶对称为了简化运算产生如下计算方法如果一个FIR滤波有一个冲击响应，h(0,h(1),…,h(N-1),和想x（n）描 绘的时常滤波你，输出滤波y（n）的n给出以下方程式： 算法设计设N=20,FIR滤波器算法设计如下：1.y（n）=h(0)x(20)+h(1)x(19)+,…,+h(19)x(1)+h(20)x(0)AR4,FIR_DATA_PAR6,INBUF_PAR7,OUTBUF_P程序分别用AR4指向FIR_DATA_PC，AR6指向INBUF_P,AR7指向OUTBUF_P相应参数计算低通滤波器的滤波系数的计算round(fir1(20,1750/9600 )*32768)滤波器输入信号参数编写程序源程序 .globalstart,fir .mmregsCOFF_FIR_START: .sect coff_fir .include.incK_FIR_BFFR .set20d_data_buffer .usectfir_bfr,64FIR_DP .usectfir_vars,0d_filin .usectfir_vars,1output .usectfir_vars,1input .usect fir_vars,1d_filout .usectfir_vars,100hstacksize .set 256stack .usect fir_vars,stacksize.asg AR4,FIR_DATA_P .asg AR6,INBUF_P.asg AR7,OUTBUF_P.asg AR3,OUTBUF.asgAR2,INBUF .sectfir_prognopstart: stm #stack+stacksize,SP LD#FIR_DP,DPSTM#d_data_buffer,FIR_DATA_P; RPTZ A,#K_FIR_BFFR-1;RC=31STLA,*FIR_DATA_P+;STM#d_filin,INBUF_PSTM#d_filout,OUTBUF_P STM#output,OUTBUFSTM#input,INBUF STM#100h,BKfir_loop: NOPLD*INBUF_P,ASTL A,*INBUFCALLfirSTHA,*OUTBUF_P