二维自适应滤波器设计

1、苏州大学本科生毕业设计（论文）学院(部)电子信息学院题 目二维自适应滤波器设计目录摘要1Abstract2第一章 绪论31.1背景和意义31.2自适应滤波器的原理及应用31.2.1自适应滤波器的基本原理41.2.2自适应滤波器的应用51.3自适应滤波算法51.4论文结构5第二章2D-LMS自适应滤波算法72.1 LMS算法简介72.2 2D-LMS算法82.3 MATLAB仿真92.3.1仿真条件92.3.2仿真结果与分析10第三章2D-NLMS自适应滤波算法123.1 NLMS简介123.2 2D-NLMS算法123.3 MATLAB仿真123.3.1仿真条件123.3.2仿真结果与分析13

2、第四章2D-APA自适应滤波算法164.1 APA简介164.2 2D-APA算法164.3 MATLAB仿真174.3.1仿真条件174.3.2仿真结果与分析18第五章 二维自适应噪声消除（2D-ANC）215.1 简介215.2 MATLAB仿真215.2.1仿真条件215.2.2仿真结果22第六章 总结与展望246.1总结246.2展望24参考文献25致谢27摘要自适应滤波器通过使用递归算法来估计未知系统，具有较强的适应性和跟踪能力。而将一维自适应滤波算法扩展到二维结构便得到二维自适应滤波器。二维结构对自适应滤波器的图像数据处理能力有很大提高，二维自适应滤波器能够考虑到二维固有的非平稳统

3、计特性数据以及二维统计相关性。目前，二维自适应滤波器已应用于各种图像处理应用，如图像恢复，图像增强，自适应噪声消除和二维系统识别等。本文首先对自适应滤波器的原理进行了讨论，详细介绍了三种二维自适应滤波器算法，即二维最小均方(2D-LMS)算法、二维归一化最小均方(2D-NLMS)算法和二维反射投影(2D-APA)算法，并通过使用MATLAB软件对其仿真得到MSD曲线，比较不同算法的性能。然后将三种算法用于图像处理中的图像噪声消除，对图像的去噪性能进行比较。关键词：二维自适应滤波器，噪声消除，MATLAB仿真作 者：王映涛指导老师：倪锦根AbstractThe adaptive filter e

4、stimates the unknown system by using a recursive algorithm and has strong adaptability and tracking ability. A one-dimensional adaptive filter algorithm is extended to a two-dimensional structure to obtain a two-dimensional adaptive filter. The two-dimensional structure greatly improves the image da

5、ta processing capability of the adaptive filter. The two-dimensional adaptive filter can take into account inherent non-stationary statistical characteristic data and two-dimensional statistical correlation. Currently, 2D adaptive filters have been applied to various image processing applications su

6、ch as image restoration, image enhancement, adaptive noise cancellation, and 2D system identification. This paper first discusses the principle of adaptive filter, and introduces three two-dimensional adaptive filter algorithms in detail, namely two-dimensional least mean square (2D-LMS) algorithm,

7、two-dimensional normalized least mean square (2D-NLMS) algorithm and two-dimensional reflection projection (2D-APA) algorithm are used to simulate the MSD curve by using MATLAB to compare the performance of different algorithms. Then the three algorithms are used for image with noise elimination in

8、image processing, and the image cancellation performance is compared.Keywords: Two-dimensional adaptive filters, Noise cancellation, MATLAB simulation Written by Wang ying-tao Supervised by Ni Jin-gen第一章 绪论1.1背景和意义自适应滤波器的开发与应用是在现今自适应信号处理领域中非常活跃的课题之一。由于其具备较强的自学习、自跟踪的能力以及算法的易实现性等特点，其理论及方法在通信、雷达、声呐、控制、

9、导航、地震学以及生物医学等领域有着广泛应用1。传统的滤波器算法，如维纳滤波器是一种基于最小均方误差准则建立的最优线性滤波器，它利用信号的统计特性，对含噪声信号进行滤波。当输入过程是广义平稳的，且输入信号的统计特性已知的情况下，维纳滤波器的输出与期望输出误差的均方值达到最小。但在实际应用中，信号的输入统计特性可能是未知的，也可能是变化的。在这种情况下，自适应滤波器有较好的滤波性能，不需要知道所处理信号的所有统计特性，且系数也不是固定不变的，而是能通过多次迭代的运算来实现系数的更新变化，完成滤波运算。即自适应滤波器能够根据环境情况的变化来改变滤波器的参数和结构。而在滤波器算法中，滤波器系数是根据运

10、算过程中的误差函数或均方误差来进行调整的2。目前的自适应滤波算法主要有最小均方算法(LMS)、仿射投影算法(APA)和最小二乘法(RLS)以及一些改进的算法，如归一化LMS算法(NLMS)，变步长LMS算法(VSS-LMS)等。二维自适应滤波器因为能够考虑到二维固有的非平稳统计特性以及二维的统计相关性，在图像处理中有着广泛应用，如图像恢复，图像增强，自适应噪声消除，二维自适应线性增强器以及二维系统识别3。而将一维自适应滤波算法扩展到二维结构便得到二维自适应滤波器，并由此提出了2D-LMS算法(二维最小均方算法)，2D-NLMS算法(归一化的二维最小均方算法)，2D-APA算法(基于仿射投影的二

11、维自适应算法)。1.2自适应滤波器的原理及应用自适应滤波器一般包括两个基本过程，其中滤波过程通过对一系列的输入数据产生输出响应，而自适应过程就是为了给滤波过程提供一种可调参数的自适应算法，这两个过程是相互影响的。1.2.1自适应滤波器的基本原理自适应滤波器的一般结构如图1-1所示，设为输入信号，为输出信号，为期望信号，而误差信号。通过利用误差信号构建目标函数，并通过使误差信号尽可能的趋近于0，实现自适应滤波器的输入信号与输出信号的最佳匹配。图1-1自适应滤波器一般结构图自适应滤波器主要包括三个部分：(1)自适应滤波器的结构：自适应滤波器可以通过多种不同结构实现，并且不同结构还会对计算复杂度、收

12、敛速度等性能造成影响。如FIR滤波器又称为横向滤波器；格型预测器，具有模块结构；脉动阵列，即并行计算网络。本文主要参考横向滤波器，以此为基础建立的自适应横向滤波器是自适应滤波器的基本结构。(2)性能判据：自适应滤波过程是一个迭代收敛过程，本文采用稳态均方误差(Mean Square Deviation，MSD)来分析自适应滤波器的算法性能。根据图1-1，将一维扩展到二维，则均方差可以定义为：(1.1)其中，表示未知系统，表示网络中第节点的权值向量。(3)算法：自适应算法使用不同的目标函数，通过输入信号和期望信号的迭代来确定滤波器参数。不同的算法对计算复杂度，最优解的存在性等自适应特性有直接影响

13、。1.2.2自适应滤波器的应用通过对自适应滤波器的期望信号和输入信号使用不同的选取方式，就可以得到自适应滤波器的不同应用。自适应滤波器的应用主要分为四种基本类型，包括系统的辨识(图1-1就是一个简单的一维自适应滤波器的辨识模型)；还包括干扰消除(例如噪声消除)；逆模型；预测。本次毕业设计主要涉及二维自适应滤波器在图像噪声消除的应用以及二维自适应滤波器的辨识的应用。1.3自适应滤波算法本次毕业设计涉及的自适应滤波算法目前主要分为以下几种：(1)最小均方算法(Least Mean Square, LMS)：是目前使用广泛的滤波器算法，具有计算复杂度低，稳定性强，易于收敛等特点。而二维最小均方算法(

14、2D-Least Mean Square, 2D-LMS)实际上是LMS算法的在二维上的扩展，用于非平稳图像的估计，但其特征值差异具有高度敏感的特性，并且其收敛速度较慢，在许多应用中不适用。(2)归一化最小均方(Normalized Least Mean Square, NLMS)算法：通过将LMS算法进行归一化处理就能得到NLMS算法，NLMS算法考虑了滤波器输入幅度对其收敛速度的影响。而且可以防止输入信号的突变带来的自适应系数的突变，使算法性能得到提升。将其扩展到二维便得到二维归一化最小均方(2D-Normalized Least Mean Square, 2D-NLMS)算法4。(3)仿

15、射投影(Affine Project, APA)算法：当输入数据高度相关时，使用仿射投影算法能够提升收敛性能，但计算复杂度大大提高。将其扩展到二维得到基于反射投影的二维自适应算法(2D-Affine Project, 2D-APA)算法5。该算法收敛速度快，跟踪能力强，缺点就是计算复杂度过高。1.4论文结构本文的主要内容如下：第一章介绍了自适应滤波器的背景及意义，自适应滤波器的原理与应用，本次毕业设计使用到的几种自适应算法，并说明了本文的内容及主要工作。第二章介绍了二维最小均方算法(2D-Least Mean Square, 2D-LMS)，并进行仿真，分析其稳态性能。第三章介绍了二维归一化最

16、小均方(2D-Normalized Least Mean Square, 2D-NLMS)算法，并与2D-LMS算法进行仿真比较6。第四章介绍了基于仿射投影的二维自适应算法(2D-Affine Project, 2D-APA)，并与2D-LMS，2D-NLMS两种算法进行仿真比较。第五章将2D-LMS，2D-NLMS，2D-APA三种算法应用于图1-3所示系统，对同一含噪图像进行自适应噪声消除，并通过去噪结果比较几种算法的性能。第六章对全文进行简单总结和展望。第二章2D-LMS自适应滤波算法2.1 LMS算法简介最小均方算法(Least Mean Square, LMS)由B.维德罗于上世纪7

17、0年代提出7。LMS算法利用线性组合器实现滤波器，具有结构简单，性能稳定等特点，得到了广泛应用。该滤波器的输出为 (2.1.1)误差信号可以表示为 (2.1.2)其中为输出信号，为期望信号。该算法的滤波器系数由输入信号与误差信号来决定，从而实现最小均方误差准则8。对于，求估计误差与抽头输入的内积，得到自适应滤波系数相对于维纳解的误差。可以表示为 (2.1.3)其中，利用最陡下降法可以得到权值向量 (2.1.4)其中，为代价函数的梯度向量。 表示迭代次数，为步长参数且，因子的引入可以使计算方便。由式（2.1.4）可知，在到次迭代时，校正量 (2.1.5)设和是联合广义平稳随机过程，可以通过矩阵和

18、向量来估计。则 (2.1.6)其中 (2.1.7)由此得到的梯度算法即为LMS算法。将(2.1.6)，(2.1.7)代入(2.1.4)中可以得到其迭代方程为 (2.1.8)2.2 2D-LMS算法将最小均方算法(Least Mean Square, LMS)扩展到二维，便得到二维最小均方算法(2D-Least Mean Square, 2D-LMS)。该算法采用有限脉冲响应（FIR）线性模型。为了简化估计任务并减少实时应用中的计算量，提出了这种方法。让它作为一个线性二维FIR模型的输入，每个点定义在规定间隔的区间上，即 , 其中M1和M2指定输入数据的顺序。为二维有限脉冲响应（FIR）数字滤波

19、器的输出信号，表达式如下： (2.2.1)其中是输入信号，为滤波器权矢量， 和 指定了FIR滤波器的阶数。该二维未知系统是由满足正态分布的随机数组成的的矩阵，为了比较方便，将其转换为长度为的一维列向量。在迭代期间，这些列向量可以表示为 (2.2.2) (2.2.3)其中是迭代次数且。从(2.2.2)和(2.2.3)可知，二维FIR滤波器的输出可以表示为 (2.2.4)2D-LMS算法同LMS算法一样基于最陡下降法，并且根据该原理，二维滤波器权矢量可以表示为 (2.2.5)其中为步长，为迭代k次的误差信号，可以表示为 (2.2.6)另外，为期望信号。2D-LMS算法是实现二维维纳滤波器的一种实用

20、方案，没有明确求解维纳-霍夫方程。 因此，2D-LMS自适应滤波器的迭代公式可以表示为 (2.2.7)2.3 MATLAB仿真2.3.1仿真条件本节基于MATLAB软件，使用2D-LMS算法实现了自适应滤波器的系统辨识结构的应用9。图2-1二维自适应滤波器辨识结构图如图2-1二维自适应滤波器辨识结构图所示，仿真条件如下10：(1) 未知系统：在公式(2.2.2)和(2.2.3)中已经知道，未知二维系统是一个的矩阵，为了与自适应滤波器方便比较，将其转换为长度为的一维列向量。本次实验中，通过randn函数产生正态分布的矩阵，均值为0，方差为1，该一维列向量为。(2) 输入信号：首先由randn函数

21、产生一个二维高斯白噪声矩阵，大小为，将其通过如公式(2.3.1)所示的二维低通滤波器，就能得到有色输入信号。(2.3.1)其中滤波器模型系数，。注：将矩阵大小设为主要考虑到第五章二维自适应噪声消除应用中处理的图片大小。(3)误差信号：误差信号可由公式(2.1.2)扩展到二维得到，公式如下： (2.3.2)(4)测量噪声：是一个均值为0，噪声方差为0.001的二维高斯白噪声，并且与输入信号互相独立11。(5)性能判据：采用稳态均方误差(Mean Square Deviation，MSD)来分析二维最小均方算法(2D-Least Mean Square, 2D-LMS)的性能，单位为分贝(dB)。

22、仿真结果都是经过30次独立实验平滑后得到的结果12。2.3.2仿真结果与分析如图2-2所示，实现了基于2D-LMS算法的辨识系统的仿真。图2-2 2D-LMS算法仿真图由于步长因子的选取会对2D-LMS算法的性能产生影响，所以本次实验分别采用了不同步长因子进行仿真，步长因子取值，和。从图中可以知道，当步长因子时，收敛曲线大约在迭代3000次达到稳定，收敛值为；当步长因子时，收敛曲线大约在迭代2000次达到稳定，收敛值为；当步长因子时，收敛曲线大约在迭代1000次达到稳定，收敛值为。通过以上结果可以得出，步长因子影响算法性能。当步长较大时(如0.012)， 2D-LMS算法收敛速度较快，但稳态偏

23、差大，即对未知向量的估计时间较快，但是估计的准确度不高；而步长较小时(如0.006)，稳态偏差会相对较小，但收敛速度缓慢，即估计时间较慢但准确度较高。第三章2D-NLMS自适应滤波算法3.1 NLMS简介归一化最小均方(Normalized Least Mean Square, NLMS)算法相较于LMS算法，具有更快的收敛速度13。其迭代公式可以表示为 (3.1.1)其中的选取必须实现快速收敛，所以为了提高收敛速度，可以通过选取合适的值来尽可能地减小瞬时平方误差14。值的表达式如下 (3.1.2)将（3.1.2）代入（3.1.1）可得 (3.1.3)为了控制失调量，需要引入固定收敛因子，同时

24、为了避免较小时，步长太大，引入常数。所以归一化NLMS算法的迭代公式为 (3.1.4)3.2 2D-NLMS算法将一维NLMS算法扩展到二维，便得到了二维归一化最小均方(2D-Normalized Least Mean Square, 2D-NLMS)算法。和一维NMLS算法一样，为了克服2D-LMS算法收敛速度慢的缺点，提出了2D-NLMS算法 ，考虑了滤波器输入幅度的影响。2D-NLMS算法的迭代公式为 (3.2.1)其中为了防止较小，引入的常数。3.3 MATLAB仿真3.3.1仿真条件本节基于MATLAB软件，使用2D-NLMS算法实现了自适应滤波器的系统辨识结构的应用，如图2-1二维

25、自适应滤波器辨识结构图所示，仿真条件如下15：(1)未知系统：同2.3.1中的(1)，未知二维系统是一个的矩阵，为了与自适应滤波器方便比较，将其转换为长度为的一维列向量。本次实验中，通过randn函数产生正态分布的矩阵，均值为0，方差为1，该一维列向量为。(2)输入信号：首先由randn函数产生一个二维高斯白噪声矩阵，大小为，将其通过如公式(2.3.1)所示的二维低通滤波器，就能得到有色输入信号。其中滤波器模型系数，。注：将矩阵大小设为主要考虑到第五章二维自适应噪声消除应用中处理的图片大小。(3)误差信号：误差信号为。 (4)测量噪声：是一个均值为0，噪声方差为0.001的二维高斯白噪声，并且

26、与输入信号互相独立。(5)性能判据：采用稳态均方误差(Mean Square Deviation，MSD)来分析二维归一化最小均方(2D-Normalized Least Mean Square, 2D-NLMS)算法的性能，单位为分贝(dB)。仿真结果都是经过30次独立实验平滑后得到的结果。3.3.2仿真结果与分析(1)实现基于2D-NLMS算法的辨识系统的仿真。图3-1 2D-NLMS算法仿真图如图3-1所示，由于步长因子的选取会对2D-NLMS算法的性能产生影响，所以本次实验分别采用了不同步长因子进行仿真，步长因子分别取值，和。从图中可以知道，当步长因子时，收敛曲线大约在迭代18000次

27、达到稳定，收敛值为；当步长因子时，收敛曲线大约在迭代7000次达到稳定，收敛值为；当步长因子时，收敛曲线大约在迭代4000次达到稳定，收敛值为。通过以上结果可以得出，步长因子影响算法性能。当步长较大时(如0.35)，2D-NLMS算法收敛速度较快，但稳态偏差大，即对未知向量的估计时间较快，但是估计的准确度不高；而步长较小时(如0.10)，稳态偏差会相对较小，但收敛速度缓慢，即估计时间较慢但准确度较高16。(2)将2D-NLMS算法与2D-LMS算法进行比较。在图3-2中，将2D-LMS算法的步长因子设置为0.001，2D-NLMS算法的步长因子设置为0.09，使得两种算法在的收敛值下收敛。在相

28、同的收敛值下，2D-LMS算法大约在25000次达到稳定，而2D-NLMS算法大约在21000次达到稳定。图3-2 2D-NLMS算法与2D-LMS算法对比图一 在图3-3中，将2D-LMS算法的步长因子设置为0.006，2D-NLMS算法的步长因子设置为0.45，使得两种算法在迭代3000次时达到稳定。在相同的迭代次数下，2D-NLMS算法的收敛值，而2D-LMS算法的收敛值。图3-3 2D-NLMS算法与2D-LMS算法对比图二 通过以上仿真结果可以得出，在达到相同的收敛值时，2D-NLMS算法的迭代次数更少，在迭代次数相同时，2D-LMS算法的稳态偏差更低，即2D-NLMS算法的收敛速度

29、更快，对未知系统的估计更加准确17。第四章2D-APA自适应滤波算法4.1 APA简介仿射投影算法(Affine Project, APA)是基于归一化最小均方误差算法(NLMS)提出的，是NLMS算法的一种多维推广(在投影阶数的值为1时，该算法等价于NLMS算法)18。仿射投影算法提高了收敛速度，改善了跟踪性能，但计算复杂度大大增加，相对的抗干扰能力较弱19。仿射投影算法的主要思想是重复利用过去的数据信号，以提高自适应滤波算法的收敛速度。它将当前的输入与前一时刻的输入组成一个的多维仿射投影矩阵，同时使用在同一次迭代中，即 (4.1.1)在这个条件下，每次迭代时的误差已经不再是一个数值，而是一

30、个长为的向量，并且期望响应，自适应滤波器输出也已经成为一个长度为的向量。根据NLMS算法的权值向量迭代公式（3.1.4），可以整理出APA算法的迭代公式： (4.1.2)公式中的也由一个较小的常数变成一个矩阵，本文中 ,其中为阶的单位矩阵。4.2 2D-APA算法将仿射投影算法(APA)扩展到二维，便得到二维仿射投影算法(2D-Affine Project, 2D-APA)。2D-APA算法与2D-LMS算法相比，具有更好的收敛性能，然而计算复杂度也大大增加20。该算法可以表示为使下列目标函数最小化的二维自适应滤波器： (4.2.1)期望信号为 (4.2.2)在2D-APA算法中，使用的多维仿

31、射投影矩阵来引入参数和,同时应用于同一迭代中，两个参数的取值范围 和。2D-APA算法的输入信号为矩阵，它的维度为，可以表示为 (4.2.3)其中 (4.2.4)这是维数为的矩阵。 此外，（4.2.2）中的是维数为的期望信号： (4.2.5)从上面可以看出，2D-APA算法的自适应滤波器系数的迭代公式可以表示为 (4.2.6)其中误差信号可以表示为 (4.2.7)4.3 MATLAB仿真4.3.1仿真条件本节基于MATLAB软件，使用2D-NLMS算法实现了自适应滤波器的系统辨识结构的应用，如图2-1二维自适应滤波器辨识结构图所示，仿真条件如下：(1)未知系统：同2.3.1中的(1)，未知二维

32、系统是一个的矩阵，为了与自适应滤波器方便比较，将其转换为长度为的一维列向量。本次实验中，通过randn函数产生正态分布的矩阵，均值为0，方差为1，该一维列向量为。(2)输入信号：首先由randn函数产生一个二维高斯白噪声矩阵，大小为，将其通过如公式(2.3.1)所示的二维低通滤波器，就能得到有色输入信号。其中滤波器模型系数，。注：将矩阵大小设为主要考虑到第五章二维自适应噪声消除应用中处理的图片大小。(3)误差信号：误差信号为。(4)正则化因子为0.001。(5)值为投影阶数，当值为1时，该算法等价于2D-NLMS算法。(6)测量噪声：是一个均值为0，噪声方差为0.001的二维高斯白噪声，并且与

33、输入信号互相独立。(7)性能判据：采用稳态均方误差(Mean Square Deviation，MSD)来分析二维仿射投影算法(2D-Affine Project, 2D-APA)的性能，单位为分贝(dB)。仿真结果都是经过30次独立实验平滑后得到的结果。4.3.2仿真结果与分析(1)实现基于2D-APA算法的辨识系统的仿真(投影阶数)。图4-1 2D-APA算法仿真图如图4-1所示，由于步长因子的选取会对2D-APA算法的性能产生影响，所以本次实验分别采用了不同步长因子进行仿真，步长因子分别取值，和。从图中可以知道，当步长因子时，收敛曲线大约在迭代5000次达到稳定，收敛值为；当步长因子时，

34、收敛曲线大约在迭代2000次达到稳定，收敛值为；当步长因子时，收敛曲线大约在迭代1000次达到稳定，收敛值为。通过以上结果可以得出，步长因子影响算法性能。当步长较大时(如0.9)， 2D-APA算法收敛速度较快，但稳态偏差大，即对未知向量的估计时间较快，但是估计的准确度不高；而步长较小时(如0.05)，稳态偏差会相对较小，但收敛速度缓慢，即估计时间较慢但准确度较高。(2)投影阶数的不同取值的2D-APA算法对比仿真。对不同投影阶数的2D-APA算法进行仿真，如图4-2，本次仿真的步长因子均为0.05。当投影阶数时，收敛曲线在迭代41000次时达到稳定，收敛值；当投影阶数时，收敛曲线在迭代150

35、00次时达到稳定，收敛值；当投影阶数时，收敛曲线在迭代4000次时达到稳定，收敛值；当投影阶数时，收敛曲线在迭代2000次时达到稳定，收敛值。通过本次仿真结果，可知当较小时，算法失调量较小，但算法的收敛速度较慢。反之较大时算法收敛较快，但失调量较大。其中当值为1时，该算法等价于2D-NLMS算法。图4-2 不同投影阶数的2D-APA算法仿真(3)将2D-APA算法与2D-NLMS算法，2D-LMS算法进行比较。在图4-3中，将2D-LMS算法的步长因子设置为0.001，2D-NLMS算法的步长因子设置为0.09，2D-APA算法的步长因子设置为0.015。使得三种算法在的收敛值处收敛。在相同的

36、收敛值下，2D-LMS算法大约在25000次达到稳定，2D-NLMS算法大约在21000次达到稳定，2D-APA算法大约在12000次达到稳定。图4-3 2D-APA算法与2D-LMS，2D-NLMS算法对比图一在图4-4中，将2D-LMS算法的步长因子设置为0.006，2D-NLMS算法的步长因子设置为0.45，将2D-APA算法的步长因子设置为0.05，使得三种算法在迭代3000次时达到稳定。在相同的迭代次数下，2D-NLMS算法的收敛值，2D-LMS算法的收敛值，2D-APA算法的收敛值。图4-4 2D-APA算法与2D-LMS，2D-NLMS算法对比图二通过以上仿真结果可以得出，在达到

37、相同收敛值时，2D-APA算法的迭代次数最少；在迭代次数相同时，2D-LMS算法的稳态偏差最低。即收敛速度2D-APA算法2D-NLMS算法2D-LMS算法 。2D-APA算法对未知系统的估计更加准确，但计算复杂度较大。第五章 二维自适应噪声消除（2D-ANC）5.1 简介在自适应滤波器的应用中，自适应噪声消除的目的是从接受信号中减去噪声，从而改善信噪比21。二维自适应噪声消除（2D-ANC）可用于二维图像的处理，如图片去噪，图像回复，图像增强。5.2 MATLAB仿真5.2.1仿真条件如图5-1二维自适应噪声消除结构图所示：期望信号为原图信号(即输入信号)与噪声信号的叠加，噪声信号由高斯白噪

38、声信号通过一个二维的低通滤波器得到。仿真条件如下：图5-1二维自适应噪声消除(2D-ANC)结构图(1)仿真了两组图片：输入信号是一个像素的原始图像；输入信号是一个像素的原始图像。(2)由randn函数产生两个二维高斯白噪声矩阵和，大小分别为和，一般将其通过二维低通滤波器就能得到噪声信号和22。(3)本次毕业设计通过将输入信号与噪声信号转换为长度为90000的两个列向量；将输入信号与噪声信号转换为长度为65536的两个列向量。其中噪声信号由转换为列向量的高斯噪声信号通过一个一阶系统得到。使用了2D-LMS算法，2D-NLMS算法和2D-APA算法，对两个图像进行了去噪(4)在二维自适应噪声消除

39、的应用中，2D-LMS算法的步长因子设置为0.000004，2D-NLMS算法的步长因子设置为0.004，2D-APA算法的步长因子设置为0.00423。5.2.2仿真结果图5-2a 原始图像 图5-2b 原始图像图5-3a噪声污染图像 图5-3b噪声污染图像图5-4a 2D-LMS去噪图像 图5-4b 2D-LMS去噪图像图5-5a 2D-NLMS去噪图像 图5-5b 2D-NLMS去噪图像图5-6a 2D-APA去噪图像 图5-6b 2D-APA去噪图像通过对三组仿真结果图进行对比，2D-APA算法的图像去噪性能比2D-LMS算法和2D-NLMS算法更好。第六章 总结与展望6.1总结本文通

40、过对自适应滤波器算法的阐述，介绍了自适应滤波器结构，滤波器的应用，包括辨识结构的应用与自适应噪声消除的应用。简单介绍了常用的自适应滤波算法：LMS算法，NLMS算法，APA算法。将一维自适应滤波算法扩展到二维结构得到二维自适应滤波器 ：2D-LMS算法(二维最小均方算法)，2D-NLMS算法(归一化的二维最小均方算法)，2D-APA算法(基于仿射投影的二维自适应算法)。介绍了其推导公式及迭代算法。通过仿真，比较了三种算法的性能，得到了以下结论：算法的性能取决于步长因子，如果步长因子大，收敛较快，但稳态偏差大；如果步长因子小，稳态偏差会降低，但收敛缓慢。2D-NLMS算法的收敛速度比2D-LMS

41、算法快，而2D-APA算法的整体性能要优于前两种算法。2D-APA算法对未知系统的估计更加准确，在二维噪声消除的应用中，去噪效果较好，得到的去噪图片较为清晰。唯一的缺点就是计算复杂度较大。6.2展望本文虽然实现了三种算法2D-LMS算法(二维最小均方算法)，2D-NLMS算法(二维归一化的最小均方算法)，2D-APA算法(二维仿射投影算法)。但存在计算复杂度较高的问题。这也是二维自适应滤波器最重要的问题之一，所以可以提出一种方法来降低计算复杂度，如选择部分更新自适应滤波器系数。因此如何优化二维自适应滤波器算法，提收敛性能可以进一步研究。参考文献1 赫金（Haykin S）著, 郑宝玉等译. 自

42、适应滤波器原理(第五版)M. 北京电子工业出版社, 2016.2 迪尼兹（Diniz P S R）著, 刘郁林等译. 自适应滤波算法与实现（第四版）M. 北京: 电子工业出版社, 2014.3 Abadin M S E, Aali S N. The novel two-dimensional adaptive filter algorithms with the performance analysisJ. Signal Processing, 2014, 103: 348-366.4 伍星, 陈朝阳, 陈新武. 二维 NLMS 自适应算法的结构设计J. 计算机与数字工程2005, 12: 8

43、-10.5 倪元敏, 王军. 基于投影算法的二维自适应滤波器J. 重庆工学院学报(自然科学), 2008, 08: 112-118.6 伍星. 用于小目标检测的二维NLMS自适应滤波器的设计与实现D. 华中科技大学, 2005.7 胡伟. 子带自适应滤波器及其应用D. 苏州大学, 2016.8 李宁. LMS自适应滤波算法的收敛性能研究与应用D. 哈尔滨工程大学, 2009.9 贺宽, 黄涛. 基于Matlab的自适应滤波器设计J. 武汉理工大学学报, 2008, 02: 71-73.10 Ohki M, Hashiguchi S. Two-dimensional LMS adaptive f

44、iltersJ. IEEE Transactions on Consumer Electronics, 1991, 37: 66-73.11 Azimi-Sadjadi M R, Pan H. Two-dimensional block diagonal LMS adaptive filteringJ. IEEE Transactions on signal processing, 42(9): 2420-2429.12 Glentis G O. An efficient affine projection algorithm for 2D FIR adaptive filtering and linear predictionJ. Signal Processing, 2005, 86: 98-116.13 Douglas S C. Analysis and implementation of the max-NLMS adaptive filterJ. Proceedings of 29th