数字滤波器的MATLAB设计与DSP上的实现(2)

假设实际工程需要设计一个线性相位带通FIR滤波器指标如下：

fn=[1000,1375,3625,4000];

a=[0,1,0];

dev=[0.0005,0.05,0.0005];

阻带最小衰减60dB，再设采样频率fs=10kHz，则根据阻带最小衰减来选择凯塞（Kaiser）窗，利用MATLAB的求阶函数和FIR滤波器的设计函数，可以快速地设计出所需的数字滤波器。这两个设计函数如下：

[N,Wn,beta,ftype]=kaiserord(fn,a,dev,fs)；

b=fir1(n,Wn,ftype,Kaiser(n 1,beta)；

最后，利用的滤波器分析函数freqz分析所设计出的滤波器的幅频特性和相频特性，并用图形显示函数plot将它们显示出来，如图1所示。由图1可见，设计结果满足指标要求。

2 数字滤波器的实现方法

数字滤波器的实现方法一般有以下几种。

①采用加法器、乘法器、延时器设计专用的滤波电路。

②在通用计算机系统中加上专用的加速处理机设计实现。

③用通用的可编程DSP芯片实现。

④用专用的DSP芯片实现。在一些特殊的场合，要求的信号处理速度极高，用通用DSP芯片很难实现。这种芯片将相应的滤波算法在芯片内部用硬件实现，无需进行编程。

⑤采用FPGA/CPLD设计实现。

在上述几种方法中，第②种方法的缺点是速度较慢，一般可用于DSP算法的模拟。第①、④种方法专用性强，应用受到很大的限制。第③、⑤种方法都可以通过编程来实现各种数字滤波，但是，第③种因有专用的指令来实现滤波运算编程实现容易，而第⑤种方法编程实现较为困难。

3 数字滤波器的DSP实现

DSP是一种实时、快速、特别适合于实现各种数字信号处理运算的微处理器。由于它由具有丰富的硬件资源、改进的哈佛结构、高速数据处理能力和强大的指令系统，而在通信、航空、航天、雷达、工业控制、网络及家用电器等各个领域得到广泛应用。DSP分为定点和浮