基于DSP芯片实现数字信号的相关、卷积、滤波运算.pdfVIP

基于DsP芯片实现 数字信号的相关、卷积、滤波运算 高小玲王擎天 总装备部指挥技术学院 摘要 本文介绍了I骼A100开发板的结构，并针对该结构提出了对数字信号进行相 关、卷积、有限长数字滤波的运算方法．最后对无限长数字信号的滤波运算进行了探讨。 关键词 相关卷积有限长数字滤波无限跃数字滤波 1 IMsAl00芯片介绍 IMsAloO是英国lnmos公司生产的高速、高精度的32级可级联的数字式横向滤波器 它的主要部分是乘法累加阵列，其原理如图1所示： X 在这种实现方式中，输入信号同时送到所有32个乘法器中，并以分散的形式完成相乘 和相加。每个输入样点数据都是同时送到A100的32级的每一级。在每一级中，当前输入 的样点与存储在存储器中的系数相乘，并与前一级中的经一个时钟周期延时后的输出值相 加，在t=kT时滤波器的输出为： 其中，x(kT)代表第k个输入样点数据，c(o)到c(N—1)是N阶的系数。 y(七r)=c(O)x(七71)+c(1)z【(七一1)r+·-·+c(Ⅳ一1)z【七一Ⅳ+1)丁】 AloO有四个接口，可通过它们进行数据交换，其中存储器接口允许访问下列寄存器： 系数寄存器，配置和状态寄存器，用丁乘法——累加器组的数据输入输出寄存器，一般的 微机可通过存储器接口与AloO配合。 2用A100实现相关、卷积 2 1相关 数字信号的相关运算可表示为： Ⅳ一1 —f￡砂(，71)=1／Ⅳ≥，x【(女一，”)71]y(￡7’) ±=0 由丁Aloo可以实现32阶的乘法和加法，冈I比可川来实现相关。在廊用A100芯片实现 两个信号或序列x(k)和y(n)的相关时，要把两个信号中的一个(比如x(k))存在AlOO 的一套系数积存器中．这些样点应该从左到打存入相应的系数寄存器中：如果基准波形x(k) 小丁-32个样点K，则最左边的未使用的系数寄存器应置为O。例如对一个30个样点的基准 信号．其分配如图2： 0 l 2 3 … 29 30 3l 图230点相关器系数寄存器的放置 另一个信号y(n)的样点要加到A100的输入上，则相应序列应该是两个信号的互相关 函数。 2．2卷积 两个数字信号的卷积可表示为： Ⅳ一I c砂(珊丁)=1／Ⅳ∑x(￡r)y【(Ⅲ一k)丁】 卷积与相关的区别在于先要将“k)反褶，然后移位m个单位，在与x(k)相乘求和，用 A100实现时仅表现在存放系数时顺序不同，如上例，存放时为： 0 l 2 3 ．．． 29 30 3l 图3 30点卷积器系数寄存器的放置 同样将另一个信号y(n)的样点要加到Aloo的输入上，则输出序列应该是相应于两个信 号的卷积函数。 3用A100构造滤波器 3．1实现有限长数宇滤波器(FlR) 用A100构造FIR滤波器时，首先要根据设计要求采用适当的方法得到滤波器系数，根 据不同精度要求，将系数表示为4，8．12或16bits格式，滤波器可通过将系数装入A100 存储体而构成。 31．1不加窗求滤波器系数 数字滤波器是处理以频率f。抽样的信号，因此它的频响是以f。为周期的周期函数，利 用Fourier来表示其冲击响应和频率响应： H(∞)=∑自(月弘”“^(”)=1／c队C，’iH@弦”。d珊 1 … ∞。是以rad／s为单位的采样频率，对于一个滤波器，利用已定义的理想频响H(m)，由 (2)式求得冲击响应h(n)。作为一个例子，讨论截止频率为∞。的理想低通滤波器H佃) 在—∞，∞((o，内为1．O，其余为O时的冲击响应h(n)： 所求32个系数h(O)…h(31