DSP 第17讲 FIR滤波器及实现.pptVIP

2003.2.16 DSP原理及应用 第6章 应用程序设计 (3) 系数对称FIR滤波器的实现 程序清单： .text start: LD #x_new,DP SSBX FRCT STM #x_new,AR1 STM #x_old+(size-1),AR2 STM #size,BK STM #-1,AR0 PORTR PA1,#x_new ;设置页指针 ;设置小数乘法 ;AR1指向x_new单元 ;AR2指向x_old+3单元 ;设置循环缓冲区长度 ;设置双操作数增量AR0=-1 ;输入数据x(n) 第6章 应用程序设计 程序清单： FIR3: ADD *AR1+0%,*AR2+0%,A RPTZ B,#(size-1) FIRS *AR1+0%,*AR2+0%,COEF STH B,@y PORTW @y,PA0 MAR *+AR1(2)% MAR *AR2+% MVDD *AR1,*AR2+0% BD FIR3 PORTR PA1,*AR1 .end ;完成AH=x(n)+x(n-7) ;B清0，设置重复次数 ;乘法累加、加法操作 B = AH×bi + B AH=(AR1)+(AR2) ;保存y(n)结果 ;输出结果 ;修正AR1，指向New区最老数据 ;修正AR2，指向Old区最老数据 ;New区向Old区传送数据 ;循环 ;输入最新数据至New区 第6章 应用程序设计 6.1 FIR滤波器的DSP实现 6.1.5 FIR滤波器的设计实例 设计一个FIR低通滤波器，其设计参数： 滤波器阶数：N=40； 截止频率：wp=0.35?， ws=0.4?。 第6章 应用程序设计 6.1.5 FIR滤波器的设计实例 根据给定的设计参数，滤波器系数可由MATLAB中的fir2函数产生，函数调用格式： 1.由给定的设计参数确定滤波器的系数 f = [ 0 0.35 0.4 1 ]； m = [ 1 1 0 0 ]； b=fir2( 39, f, m ) 利用freqz函数可绘制滤波器的幅频、相频特性,其格式： freqz(b, 512, 1000) 第6章 应用程序设计 滤波器的频率特性: 1.由给定的设计参数确定滤波器的系数 第6章 应用程序设计 采用循环缓冲区实现FIR滤波器的源程序如下： 2．汇编源程序 .title “FIR.ASM” .mmregs .bss y，1 K_FIR_BFFR .set 40 PA0 .set 0 PA1 .set 1 FIR_COFF_TABLE .usect “FIR_COFF”,40 ;定义数据存储空间 D_DATA_BUFFER .usect “FIR_BFR”，40 ;定义数据存储空间 .data COFF_FIR_START：.word -7*32768/10000,3*32768/10000 ;b0,b1 .word 14*32768/10000,10*32768/10000 ;b2,b3 … … … .word 3*32768/10000,-7*32768/10000 ;b38b39 第6章 应用程序设计 .text .def FIR_INIT .def FIR_TASK FIR_INIT：SSBX FRCT ;设置小数乘法 STM #FIR_COFF，AR5 ;AR1指向b0单元 RPT #K_FIR_BFFR-1 ;设置传输次数 MVPD #COFF_FIR_START,*AR5+ ;系数bi传输至数据区 STM #D