一种非负自适应滤波器的制造方法
【专利摘要】本发明公开了一种非负自适应滤波器，属于数字滤波器设计领域。该滤波器由最小化系数向量的l0范数和对数绝对误差函数相结合的代价函数获得。l0范数的最小化加快了非负自适应滤波器估计未知系统的收敛速度，而对数绝对误差函数最小化提高了非负自适应滤波器的鲁棒性。该非负自适应滤波器可以应用于电子、通信系统受到大脉冲噪声干扰的场合。
【专利说明】
-种非负自适应滤波器
技术领域
[0001] 本发明公开了一种非负自适应滤波器，属于数字滤波器设计领域。
【背景技术】
[0002] 系统辨识是自适应信号处理的一个重要分支，传统的自适应信道均衡、自适应噪 声消除、自适应回声抵消、主动噪声控制等诸多问题都可W归结为系统辨识问题。在一些应 用中，由于受到系统内在的物理特性的限制，需要对待估计的系统的系数的最优值进行非 负性约束。运种非负性约束条件下的系统辨识问题是非负性约束条件下的最优化问题的一 种具体表现形式。非负性约束条件下的最优化问题在理论和工程实践中经常设及，是目前 的研究热点之一，例如非负最小二乘、非负矩阵分解等理论及其在天体物理图像去模糊、计 量化学光谱发散反卷积、高光谱图像分析等问题中的应用。然而，传统的非负最小二乘法需 要进行批处理，因此不适合用来在线处理非负性约束条件下的系统辨识问题。
[000引为了克服非负最小二乘自适应滤波器实时性差的问题，Jie Chen等人 [Nonnegative least-mean-square algorithm. IEEE Transactions on Signal Processing，2011，59( 11): 5225-5235.]提出了一种在非负性约束条件下的系统辨识方法。 由于该自适应滤波器的迭代公式的表现形式和迭代特征与最小均方(LMS)自适应滤波器具 有相似性，将其称为非负最小均方(NNLMS)自适应滤波器。NNLMS自适应滤波器结构简单且 易于实现。
[0004] 收敛速度、稳态失调和鲁棒性是自适应滤波器的Ξ个重要性能指标。收敛速度的 快慢决定了自适应滤波器逼近未知系统需要花费的时间，而稳态失调的高低决定了自适应 滤波器逼近未知系统能够达到的精度，鲁棒性决定了自适应滤波器能否收敛，运Ξ个指标 同时影响着信号处理的质量。一方面稀疏系统存在于多个领域，如视频会议和免提电话中 的回声路径。当未知系统的最优估计向量*°为稀疏时，自适应滤波器的收敛速度较慢。另一 方面，在有些情况下，自适应滤波器所处的环境非常恶劣，其采集到的信号很可能受到脉冲 噪声的干扰，脉冲噪声的幅值有时是输出误差信号中所包含的有用信号幅值的几百甚至几 千倍，基于最小均方误差准则的NNLMS自适应滤波器在该环境中稳态失调较大。
【发明内容】
[0005] 为了解决上述问题，本发明的目的是提出了一种非负自适应滤波器。用于解决稀 疏系统辨识中NNLMS自适应滤波器收敛速度慢、鲁棒性差的缺点。该非负自适应滤波器采用 抗脉冲干扰的方法和增加零吸引子的方法来更新其系数向量，从而提高非负系统辨识的性 能。
[0006] 为实现上述的方案，本发明采用如下技术特征：
[0007] -种非负自适应滤波器，其特征在于:所述非负自适应滤波器采用抗脉冲干扰的 方法和增加零吸引子的方法来更新其系数向量。
[000引较佳的，该非负自适应滤波器更新其系数向量包含W下步骤：
[0009] 1)通过η时刻的输入信号x(n)和期望信号d(n)计算误差信号e(n)，即e(n)=d(n)- wT(n)x(n)，其中，x(n) = [x(n)，x(n-l)，···，χ(η-Μ+1)]τ为由输入信号的前Μ个样值{x(n)，x (n-l)，…，x(n-M+l)}构成的输入向量，w(n) = [wl(n)，W2(n)，···，WM(n)]τ为非负自适应滤波 器的Μ个抽头系数构成的系数向量，Τ表示转置运算；
[0010] 2)将输入向量χ(η)的元素作为对角元素，生成对角矩阵Dx(n)，即Dx(n) =diag{x (n)，x(n-l)，…，x(n-M+l)}，再由滤波器系数向量w(n)、误差信号e(n)及误差信号的绝对值 e(n) I计算可W降低对脉冲干扰敏感性的分量
[OOW 3)计算非负自适应滤波器在η时刻的Μ个零吸引子邑1(11)，巧{1，2^-，1}，即邑1(11) =-曲iax{l-iki(n)，0}，其中β为较小的正常数，再由Μ个零吸引子形成加快非负自适应滤波 器收敛速度的分量f2(n)=G(n)w(n)，其中，G(n)为个零吸引子为元素的对角矩阵，即G (n)=diag{gi(n)，g2(n)，...，卵（η)};
[0012] 4)采用迭代公式w(n+l)=w(n)+yfi(n)+kf2(n)来更新非负自适应滤波器的系数向 量，其中，μ为非负自适应滤波器的步长参数，k为决定零吸引子强度的权重系数。
[0013] 本发明技术方案的原理在于：
[0014] 对一个向量取1〇范数，就是计算该向量中非零元素的个数，该值是系统稀疏特性 的典型表征。本发明将1〇范数引入建立非负自适应滤波器的代价函数，运种扩展相当于在 自适应滤波器更新中加入了零吸引子，从而加快了自适应滤波器在估计稀疏系统时的收敛 速度。此外，本发明还将一种可W降低自适应滤波器对脉冲干扰敏感性的对数函数作为建 立非负自适应滤波器的代价函数，从而提高自适应滤波器的抗脉冲干扰能力。
[0015] 与现有技术相比，本发明的有益效果是：
[0016] 既能加快对稀疏系统辨识的收敛速度，又具有较强的鲁棒性。同时实验结果表明， 本发明提出的自适应滤波器能够加快NNLMS自适应滤波器估计稀疏系统的收敛速度，并提 高其抗脉冲噪声干扰的能力。
【附图说明】
[0017] 下面结合附图及实施例对本发明作进一步描述：
[0018] 图1为非负自适应滤波器结构原理图；
[0019] 图2为自适应滤波器在实施例所述条件下归一化均方偏差的比较。
【具体实施方式】
[0020] 实施例
[0021] 本实施例采用计算机实验的方法验证非负自适应滤波器的性能。实验中使用本发 明公开的非负自适应滤波器在脉冲噪声干扰的环境下对未知稀疏系统进行辨识，并将其性 能与NNLMS自适应滤波器性能进行对比。如图1所示，本发明公开的非负自适应滤波器辨识 该未知稀疏系统包含W下步骤：
[0022] 1)通过η时刻的输入信号x(n)和期望信号d(n)计算误差信号e(n)，即e(n)=d(n)- wT(n)x(n)，其中，x(n) = [x(n)，x(n-l)，···，χ(η-Μ+1)]τ为由输入信号的前Μ个样值{x(n)，x (n-l)，…，x(n-M+l)}构成的输入向量，w(n) = [wl(n)，W2(n)，···，WM(n)]τ为非负自适应滤波 器的Μ个抽头系数构成的系数向量，Τ表示转置运算；
[0023] 2)将输入向量χ(η)的元素作为对角元素，生成对角矩阵Dx(n)，即Dx(n) =diag{x (n)，x(n-l)，…，x(n-M+l)}，再由滤波器系数向量w(n)、误差信号e(n)及误差信号的绝对值 e(n) I计算可W降低对脉冲干扰敏感性的分量
[0024] 3)计算非负自适应滤波器在η时刻的Μ个零吸引子邑1(11)，巧{1，2，。'，1}，即邑1(11) =-曲iax{l-iki(n)，0}，其中β为较小的正常数，再由Μ个零吸引子形成加快非负自适应滤波 器收敛速度的分量f2(n)=G(n)w(n)，其中，G(n)为个零吸引子为元素的对角矩阵，即G (n)=diag{gi(n)，g2(n)，...，卵（η)};
[002引 4)采用迭代公式w(n+l)=w(n)+yfi(n)+kf2(n)来更新非负自适应滤波器的系数向 量，其中，μ为非负自适应滤波器的步长参数，k为决定零吸引子强度的权重系数。
[0026] 实验中输入信号x(n)和加性噪声v(n)为零均值的高斯白噪声序列，其方差分别为 曰x2=l和〇/ = 0.001。在噪声中叠加了干扰脉冲，该脉冲噪声由伯努利过程和高斯过程的乘 积产生，即卓)=7(?)如0,其中η(η)为高斯白噪声序列，为伯努利序列。伯努利序列的 概率分布满化巧抑Η) = 1) = 0 01，巧例")=0) = 0.Q9，信号干扰比取为-lOdB。未知系统的 系数向量取为[0.8,0,0.6,0,0.5,0,0.2，0,0，-0.!，0，-0.3]τ。在非负性约束条件下，使 用高斯白噪声作为系统输入信号时，未知稀疏系统对应的最优非负系数向量为*°= [0.8， 0,0.6，0,0.5,0,0.2，0,0,0，0,0]τ。自适应滤波器的初始系数向量使用均匀分布随机函数 产生。采用归一化均方偏差（NMSD)相对于迭代次数的函数作为性能指标，其表达式为 201ogi〇( ||w°-w(n) M/lkl I )，单位为分贝（地）。所有的NMSD曲线为200次独立实验取平均 的结果。
[0027] 如图2所示，分别采用NNLMS自适应滤波器和本发明公开的非负自适应滤波器对脉 冲噪声环境中稀疏系统进行估计。NNLMS自适应滤波器步长取为μ = 0.001，本发明公开的非 负自适应滤波器步长取为μ = 0.01，其他参数为0=10、k=0.001。
[0028] 由实验结果可知：本发明公开的非负自适应滤波器比順LMS自适应滤波器具有更 快的收敛速度和更小的稳态失调。
[0029] 上述实施例只为说明本发明的技术构思及特点，其目的在于让熟悉此项技术的人 是能够了解本发明的内容并据W实施，并不能W此限制本发明的保护范围。凡根据本发明 精神实质所做的等效变换或修饰，都应涵盖在本发明的保护范围之内。
【主权项】
1. 一种非负自适应滤波器，其特征在于:所述非负自适应滤波器采用抗脉冲干扰的方 法和增加零吸引子的方法来更新其系数向量。2. 根据权利要求1所述的非负自适应滤波器，其特征在于:所述非负自适应滤波器更新 其系数向量包含以下步骤： 1) 通过η时刻的输入信号x(n)和期望信号d(n)计算误差信号e(n)，即e(n) = d(n)-wT(n) x(n)，其中，x(n) = [x(n)，x(n-l)，…，x(n-M+l)]T 为由输入信号的前 Μ 个样值{x(n)，x(n-1)，…，x(n-M+l)}构成的输入向量，w(n) = [wi(n)，W2(n)，···，WM(n)]T为非负自适应滤波器 的Μ个抽头系数构成的系数向量，T表示转置运算； 2) 将输入向量χ(η)的元素作为对角元素，生成对角矩阵Dx(n)，即Dx(n) =diag{x(n)，χ (η-1)，…，χ(η-Μ+1)}，再由滤波器系数向量w(n)、误差信号e(n)及误差信号的绝对值|e(n) 计算出用于降低对脉冲干扰敏感性的分量3) 计算非负自适应滤波器在η时刻的Μ个零吸引子gi(n)，i e {1，2，···，M}，即gi(n)=-0 max{l-i3Wl(n)，0}，其中β为较小的正常数，再由M个零吸引子形成可以加快非负自适应滤波 器收敛速度的分量f 2(n)=G(nMn)，其中，G(n)为以Μ个零吸引子为元素的对角矩阵，即G (n)=diag{gi(n)，g2(n)，···，gM(n)}; 4) 采用迭代公式<11+1)='\￥(11)+以；^(11)+1^2(11)来更新非负自适应滤波器的系数向量， 其中，μ为非负自适应滤波器的步长参数，k为决定零吸引子强度的权重系数。
【文档编号】H03H21/00GK106059531SQ201610347326
【公开日】2016年10月26日
【申请日】2016年5月24日
【发明人】倪锦根, 赵凯
【申请人】苏州大学