DNLMS滤波器的FPGA实现

0 引  言

自适应滤波器一直是信号处理领域的研究热点之一，经过多年的发展，已经被广泛应用于数字通信、回声消除、图像处理等领域。自适应滤波算法的研究始于20世纪50年代末，Widrow和Hoff等人最早提出最小均方算法(LMS)。算法由于结构简单，计算量小，易于实时处理，因此在噪声抵消，谱线增强，系统识别等方面得到了广泛的应用。为了克服定步长LMS算法中收敛速度、收敛精度及跟踪速度等对步长大小选取相互矛盾的缺点，人们提出了许多变步长LMS算法，但是，当输入信号具有强相关性时，例如语音信号，LMS算法及NLMS算法的收敛速度将急剧下降。因此，本文对进入自适应滤波器的输入信号首先进行解相关预处理，再对自适应滤波器进行归一化LMS。

很多自适应滤波器是在通用DSP处理器上通过编写软件代码实现。这种滤波器在系统实时性要求不高的情况下，能很好地满足滤波要求，但在实时性要求较高的场合(如3G移动通信，网络会议等)，这种滤波器在处理速度和抗干扰性能等方面已不能满足要求。而用FPGA硬件实现的数字滤波器就能很好的解决并行性和速度问题，而且其具有灵活的可配置特性和优良的抗干扰能力，因而自适应滤波器的FPGA实现成为近年来数字滤波器中研究的热点。Altera公司开发的DSP Builder可以完成FPGA的DSP系统设计，它结合Matlab／Simulink提供的库元件和强大的计算仿真能力，进行图形建模和仿真，还可以完成从Simulink模块文件向VHDL语言的自动转换。

1 算法概述

在标准LMS算法中，其迭代公式如下：

式中：ω(n)为自适应滤波器在时刻n的权矢量；x(n)为时刻n的输入信号矢量；d(n)为参考信号；e(n)是误差信号；μ是步长迭代因子，为了增加LMS算法的收敛速度，文献[1-4]提出了不同的变步长算法代替定步长因子μ，但是很多变步长因子的计算量大，存在开方和指数运算，故在能满足要求的情况下，归一化LMS (NLMS)依然是一个很好的选择，即用μn／[xT(n)x(n)]代替定步长因子μ为了控制失调，在更新公式中引入一个固定的收敛因子μn，这是因为所有的倒数都是基于平方误差的瞬时值而不是均方误差(MSE)。对于NLMS算法，当输入信号具有强相关性时，LMS算法及NLMS算法的收敛速度将急剧下降。因此对进入自适应滤波器的输入信号首先进行解相关预处理，再对自适应滤波器进行归一化LMS，便得解相关归一化LMS(DNLMS)即用式(4)代替式(3)：

2 基于FPGA的系统开发流程

在早期PFGA的DSP系统设计过程中，DSP开发者直接使用VHDL或Verilog HDL语音进行开发。这种开发难度较大，开发效率较低，为了克服这些问题，出现了许多新的基于FPGA的DSP开发工具，如Altera公司的DSP Builder、Xilinx公司的System Generator等。图1是Altera公司的FPGA系统级开发流程。

DSP Builder是Altera公司推出的一个面向DSP开发的系统级工具。它是作为Matlab／Simuli-nk的一个工具箱出现的。它从利用Matlab建立DSP电路模型开始。电路模型可方便地利用Simulink与DSPBuilder中提供的功能模块和IP核进行设计。电路中的功能模块和IP的技术参数、数据格式、数据类型和总线宽度等都能直接设置。电路模型建立后就能进行系统级的仿真。仿真通过后可以运行Signal Compiler，将模型转化成RTL级的VHDL代码和Tcl(工具命令语言)脚本语言。调用QuartusⅡ进行编译，QuartusⅡ根据网表文件及设置的优化约束条件进行布线布局和优化设计的适配，最后生成编程文件和仿真文件。