2011年DS作业FIR滤波器在DSP上的实现.docVIP

2011年DSP作业FIR滤波器在DSP上的实现 专业：信息与通信系统 引言 在信号处理中，滤波占有十分重要的地位。数字滤波是数字信号处理的基本方法。数字滤波与模拟滤波相比有很多优点。它除了可避免模拟滤波器固有的电压漂移、温度漂移和噪声等问题外，还能满足滤波器对幅度和相位的严格要求。数字滤波是语音处理、图像处理、频谱分析等应用中的基本处理算法。DSP是一种处理数字信号的专用微处理器，主要应用于实时快速地实现各种信号的数字处理算法。用DSP芯片实现数字滤波具有稳定性好、精确度高、不受环境影响等优点。数字滤波器分为有限冲激响应滤波器FIR和无限冲激响应滤波器IIR。 对称FIR滤波器在数字信号处理中应用十分广泛，常用于相位失真要求较高的场合，例如： 1 通信系统：调制解调器、综合业务数据网等，都要求保证数据脉冲的形状和通道中的相关时间。f2 希尔伯特变换器：要求输入输出信号正交。 3 高保真音响系统：音乐的相位失真必须减到最小，尽可能逼真地重现原来的声音等等。由于FIR是全零点的滤波器，因而系统总是稳定的。这对于系统综合是很重要的。 FIR滤波器的基本结构 设h n n 0,1,2,…,N一1 为滤波器的冲激响应，输入信号为x n ，则FIR滤波器就是要实现下列差分方程： 式 1 就是FIR滤波器的差分方程FIR滤波器的最主要的特点是没有反馈回路，因此它是无条件稳定系统。的单位脉冲响应h n 是一个有限长序列。由面的方程可见，FIR滤波箅法实际上足一种乘法累加运算，它不断地输入样本x n ，经延时做乘法累加，再输滤波结果y n 。 FIR滤波器的一个分支的延时线，把每一节的输出加权累加，得到滤波器的输。结构如图1所示，它一条均匀间隔的延迟线上对抽失信号进行加权求和构成。FIR低通滤波器，技术指标：62阶的低通滤波器，采用汉明窗，通带截止频率fp 2000Hz，阻带起始频率fs 3000 Hz，通带允许的最大衰减为Rp 0.25 dB，阻带应达到的最小衰减为As 50dB。滤波器的采样频率为fsample 9000Hz。 这里采用窗函数设计法设计FIR滤波器，首先根据阻带衰减As 50dB来选择窗形状，汉明窗和布拉克曼窗等窗函数均可提供大于50dB的衰减。由于汉明窗可提供较小的过渡带宽，所以选择汉明窗。然后编写MATLAB程序，求出技术指标中fp、fs相对应的数字频率，求出过渡带宽，由过渡带宽确定窗口长度N，求出低通滤波器的截止频率，由firl函数求出滤波器的系数。这部分MATLAB程序如下： n 62 fsample 9000; fp 2000; fs 3000; w hamming n+1 ; wp fp/ fsample/2 ws fs/ fsample/2 h fir1 n,wp,low,w ; m 0:n; figure 1 subplot 2,1,1 ;plot m,w ;grid; axis [0 n 0 1.1*max w ] ; ylabel w n ;xlabel n ; subplot 2,1,2 ;plot m,h ;grid; axis [0 n 1.1*min h 1.18*max h ] ; ylabel h n ;xlabel n ; figure 2 freqz h xlabel \omega/\pi ; ylabel |H e^j\omega |/dB subplot 2,1,2 xlabel \omega/\pi ; ylabel \theta omega / ^o h [h 0]; % 将滤波器的系数输出到firwindow.txt中，以便ＤＳＰ使用 b round h* 2^15-1 ; fid fopen fir1window.txt,w ; for i 1:8 fprintf fid,\r\t%s\t,.word ; for j i-1 *8+1: i-1 *8+7 fprintf fid,%6.0f,b j ; end fprintf fid,%6.0f,b i-1 *8+8 ; end fclose fid ; 运行后得到：n 62，wp 0.4444ws 0.6667。 图２　上图是窗函数参数，下图是滤波器的参数 图３　上图是低通滤波器相对幅值相应，下图是相位相应 由上图３可以看出所设计的滤波器满足要求。 DSP实现FIR滤波器 4.1 主程序firw.asm .mmregs in_buffer .usect in_buf,64 ；输入数据缓冲区 input .usect in_dat,100　　　　；输入数据区 output .usect out_dat,100　　　　；输出数据区 .def start .data co