FIR滤波器在DSP上的实现.docxVIP

精品文档 精品文档 PAGE PAGE12 精品文档 PAGE . FIR滤波器在DSP上的实现 引言 在信号理中，波占有十分重要的地位。数字波是数字信号理的基本 方法。数字波与模波相比有好多点。 它除了可防止模波器固有的 漂移、温度漂移和噪声等外，能足波器幅度和相位的格要求。 数字波是音理、像理、剖析等用中的基本理算法。 DSP是一 种理数字信号的用微理器，主要用于迅速地各样信号的数字 理算法。用DSP芯片数字波具有定性好、精准度高、不受境影响等 点。数字波器分有限冲激响波器 FIR和无限冲激响波器IIR。 称FIR波器在数字信号理中用十分宽泛，常用于相位失真要求高 的合，比如：(1)通信系：制解器、合数据网等，都要求保数据 脉冲的形状和通道中的有关。 f2)希伯特器：要求入出信号正交。 (3)高保真音响系：音的相位失真必减到最小， 尽可能逼真地重原来的声 音等等。由于FIR是全零点的波器，因而系是定的。于系合是 很重要的。 FIR滤波器的基本构造 h(n)(n=0,1,2,?,N一1)波器的冲激响，入信号x(n)，FIR波器 就是要下列差分方程： N 1 (1) y(n) h(k)x(nk) 0 式(1)就是FIR波器的差分方程。FIR波器的最主要的特点是没有反回路，因此它是无条件定系。它的位脉冲响h(n)是一个有限序列。由上面的方程可，FIR波法上足一种乘法累加运算，它不断地入本x(n)， 精品 . 延时做乘法累加，再输 精品 . 出滤波结果y(n)。 FIR滤波器的一个分支的延时线，把每一节的输出加权累加，得到滤波器的 输出。构造如图1所示，它由一条均匀间隔的延迟线上对抽失信号进行加权求和 组成。 x(n) h(0) h(1) h(2) h(N-1) 图1FIR滤波器的直接型构造图 y(n) 利用MATLAB确定滤波器的系数 设计一个线性相位FIR低通滤波器，技术指标：通带截止频次fp=1500Hz，阻 带开端频次fst=2250Hz，通带允许的最大衰减为 Rp=0.25dB，阻带应达到的最小 衰减为As=50dB。滤波器的采样频次为fs=16000Hz。 这里采用窗函数设计法设计 FIR滤波器，首先根据阻带衰减 As=50dB来选择 窗形状，海明窗和布拉克曼窗等窗函数均可提供大于 50dB的衰减。由于海明窗可 提供较小的过渡带宽，所以选择海明窗。然后编写MATLAB程序，求出技术指标 中fp、fst相对应的数字频次，求出过渡带宽，由过渡带宽确定窗口长度N，求出低通滤波器的截止频次，由firl函数求出滤波器的系数。这部分MATLAB程序如下： fp=1500;fst=2250;fs=16000;%输入设计指标 wp=2*fp/fs; %求归一化数字通带截止频次 ws=2*fst/fs; %求归一化数字阻带开端频次 deltaw=ws-wp; %求过渡带宽 精品 . N0=ceil(6.6/deltaw); %求窗口长度 N=N0+mod(N0+l,2); %保证窗口长度N为奇数 n=N-1; %求出滤波器的阶数n wn=(ws+wp)/2; %求滤波器的截止频次 b=firl(n，wn); %利用firl函数求出滤波器的系数 bl=round(b*10000); %将系数用整数表示 fn=fopen(‘COEF_FIR.inc’,’w’); %将系数寄存在文件COEF_FIR．inc 中 fclose(fn) 运行后得到：N=71，wp=0.1875，ws=0.2813，wn=0.2344。 再编制以下MATLAB程序段用于查验Rp和As是否知足设计要求： [H,w]=freqz(b,1); %计算频次响应 mag=abs(H); %求幅频特性 db=20*logl0(mag/max(mag)); %化为分贝值 dw=pi/512; Rp=-(min(db(1:wp*pi/dw+1))) %查验通带波动 As=-(max(db(ws*pi/dw+1:512))) %查验最小阻带衰减 求得滤波器的Rp=0.0404dB，As=51.1829dB，知足设计要求。 4.DSP实现FIR滤波器 程序设计的总体思路是：启动 ADS7864对输入的模拟信号进行模数变换，每 采集到一个数据就送人DSP滤波运算，运算结果送DAC7625变换为模拟量。不断 地重复上述过程，在DAC7625的输出端就得到滤波后的模拟信号。 为了精准地控 精品 . 制 精品 . ADS7864的采样率，使用TMS320VC5416内部的准时器控制采样时间间隔 T， 设置准时器的准时时间等于采样时间间隔 T，并让它工作在中止方式，则准时器 每过T时间就向CPU发出中止恳求，CPU响应中止恳求，转去履行中止服务程序， 在中止服务程序中读取A/D变换结果，对变换