fir滤波器在dsp上的实现.docVIP

FＩR滤波器在DＳP上的实现 引言 在信号处理中,滤波占有十分重要的地位。数字滤波是数字信号处理的基本方法。数字滤波与模拟滤波相比有很多优点。它除了可避免模拟滤波器固有的电压漂移、温度漂移和噪声等问题外,还能满足滤波器对幅度和相位的严格要求。数字滤波是语音处理、图像处理、频谱分析等应用中的基本处理算法。ＤSP是一种处理数字信号的专用微处理器,主要应用于实时快速地实现各种信号的数字处理算法。用ＤSP芯片实现数字滤波具有稳定性好、精确度高、不受环境影响等优点。数字滤波器分为有限冲激响应滤波器FIR和无限冲激响应滤波器IＩＲ。 对称FＩＲ滤波器在数字信号处理中应用十分广泛，常用于相位失真要求较高的场合，例如:（１)通信系统：调制解调器、综合业务数据网等,都要求保证数据脉冲的形状和通道中的相关时间。f２)希尔伯特变换器：要求输入输出信号正交。(3)高保真音响系统：音乐的相位失真必须减到最小，尽可能逼真地重现原来的声音等等。由于FIＲ是全零点的滤波器，因而系统总是稳定的。这对于系统综合是很重要的。 FIR滤波器的基本结构 设h(n)（ｎ＝0,1，２,…，Ｎ一1)为滤波器的冲激响应,输入信号为ｘ（ｎ）,则FIＲ滤波器就是要实现下列差分方程: 　 (１) 式(1）就是FIR滤波器的差分方程。ＦIR滤波器的最主要的特点是没有反馈回路，因此它是无条件稳定系统。它的单位脉冲响应h(ｎ)是一个有限长序列。由上面的方程可见，FＩＲ滤波箅法实际上足一种乘法累加运算,它不断地输入样本ｘ（n),经延时做乘法累加，再输出滤波结果ｙ(n）。 ＦＩR滤波器的一个分支的延时线,把每一节的输出加权累加,得到滤波器的输出。结构如图1所示，它由一条均匀间隔的延迟线上对抽失信号进行加权求和构成。 x(n) x(n) h(0) h(1) h(2) h(N-1) y(n) 　 　 　 　 　 图1 ＦIＲ滤波器的直接型结构图 利用ＭAＴLAB确定滤波器的系数 设计一个线性相位FIR低通滤波器，技术指标：通带截止频率fp=1500Hz，阻带起始频率fst=２25０ Hz，通带允许的最大衰减为Ｒp＝0.2５ ｄB，阻带应达到的最小衰减为As=50ｄＢ。滤波器的采样频率为ｆs=16000Hz。 这里采用窗函数设计法设计FIR滤波器,首先根据阻带衰减As＝50dB来选择窗形状，海明窗和布拉克曼窗等窗函数均可提供大于５0dB的衰减。由于海明窗可提供较小的过渡带宽，所以选择海明窗。然后编写ＭAＴLＡＢ程序，求出技术指标中fp、fst相对应的数字频率,求出过渡带宽，由过渡带宽确定窗口长度N,求出低通滤波器的截止频率，由fiｒｌ函数求出滤波器的系数。这部分ＭＡＴLAB程序如下: fp=1５0０；fst=２250;ｆs=16000; ??％输入设计指标 wp=２*ｆp/ｆs; ? ? %求归一化数字通带截止频率 ws=2*fst/ｆs;　?? %求归一化数字阻带起始频率 deｌtaw＝ws-wp;　?? ?? %求过渡带宽 N0=ceil(6.6/ｄelｔaw); ????％求窗口长度 N=N0＋ｍoｄ(N0+ｌ,2);?????%确保窗口长度N为奇数 ｎ＝N-１; ? ??? ％求出滤波器的阶数ｎ wｎ=(wｓ+wp)/2; ? ? %求滤波器的截止频率 b=firl(n,wn);　 ? ％利用ｆｉrｌ函数求出滤波器的系数 bｌ=roｕｎｄ（ｂ*1000０)；? ?%将系数用整数表示 ｆn=fｏpeｎ（‘COＥF_FＩR.inc’，’ｗ’)； %将系数存放在文件COEF_FIR．inc中 fcｌｏｓe(ｆｎ) 运行后得到：N=7１,wｐ=０.187５,ws=0.2813,wｎ=0．2344。 再编制以下MＡＴＬAＢ程序段用于检验Rｐ和As是否满足设计要求: ［H,w]＝freｑｚ（b，1); ? ?%计算频率响应 mag=abｓ(Ｈ）;　 ? ???%求幅频特性 ｄb=20*logｌ０(maｇ/mａx(mag)）; %化为分贝值 ｄw ＝ｐi/512; Rp= －（min（db(１:ｗp*ｐi/ｄw+1)))? %检验通带波动 Ａs= -(mａｘ(dｂ(ｗs*pｉ/dw+１：51２)))??%检验最小阻带衰减 求得滤波器的Rp=０.０404dB,As＝5１．１829dB,满足设计要求。 ＤSP实现FIR滤波器 程序设计的总体思路是:启动ＡDS7864对输入的模拟信号进行模数转换，每采集到一个数据就送人DSP滤波运算,运算结果送ＤＡC７６25转换为模拟量。不断地重复上述过程，在DAＣ7625的输出端就得到滤波后的模拟信号。为了精确地控制ADS7864的采样率,使用TＭS320VC５416内部的定时器控制采样