一种变步长lp范数LMS算法的稀疏系统辨识方法与流程

本发明属于稀疏信号处理领域，涉及一种变步长lp范数lms算法的稀疏系统辨识方法，是在含噪声环境下的稀疏系统辨识，可用于宽带无线通信系统的信道估计等。

背景技术：

自适应滤波具有较强的适应性和更好的过滤性能，该算法已经在信道均衡、线性预测、谱分析和系统识别、雷达、回波消除等领域得到了广泛应用。自适应滤波的目的是对不确定的系统或信息进行处理。这里的“不确定性”意味着所研究的信息处理过程及其环境的数学模型尚未完全确定，它包含一些未知和随机的因素，例如，一些噪音信号会以不同方式影响信息处理，这些干扰和噪声信号的统计特性常常是未知的。面对客观存在的各种不确定性，如何综合处理信息，使某些指定的性能指标达到最优或近似最优，是自适应滤波算法需要解决的问题。

widrow和hoff提出的最小均方误差(lms)自适应算法是最流行的自适应算法之一，该算法的权值更新函数采用固定步长，而步长的大小则反映了稳态误差和收敛速度之间的折衷，小的步长可给出小的稳态误差，但收敛速度较慢；大的步长可给出快的收敛速度却也导致了稳态误差较大，收敛速度和稳态误差性能不可同时兼得。为解决该问题，raymondh.kwong等研究者提出了可变步长lms算法(visualstepsizeleastmeansquare，vsslms)，该算法的步长调整由预测误差的平方值进行控制，大的预测误差值会导致步长增加以提供更快的跟踪，而小的预测误差会导致步长减小以产生更小的稳态误差。lms/f算法是由sjlim，jgharris提出的一种由lms算法和lmf算法组合而成的自适应滤波算法，也可以看成是一种步长可变的自适应滤波算法，该算法在权值向量距离最优值较远时，收敛速度和稳定性比lms算法好，在权值向量接近最优解时，稳定性比lmf算法好，但由于lms/f算法在权值向量更新中采用误差的三阶函数,因而会对附加噪声和信噪比都非常敏感。

对于稀疏未知系统，上述自适应算法并不能有效利用系统中的稀疏结构信息，参数估计性能效率不高。

技术实现要素：

针对上述现有技术中描述的不足，本发明提供一种变步长lp范数lms算法的稀疏系统辨识方法，本发明是为了提高自适应滤波算法的收敛速度和稀疏未知系统的估计精度，将系统迭代过程中产生的误差平方根引入到步长控制中，并由正常数vth平衡系统的收敛速度和稳态误差，为了挖掘未知信道的稀疏性，在代价函数中引入了lp范数，本发明在保持传统lms算法优良的收敛速度和稳态性能的基础上，进一步提高了系统的收敛速度和估计精度。

为解决上述技术问题，本发明所采用的技术方案是：

一种变步长lp范数lms算法的稀疏系统辨识方法，步骤如下：

a01.训练输入信号x(n)：

x(n)＝[x0x1…xn-1]^t；

其中，n为当前时刻，n为输入序列的信号长度，t为矩阵转置符号。

a02.初始化自适应滤波系统的抽头权重向量h(n)。

对抽头权重向量h(n)进行初始化是将所有抽头权重向量h(n)置为零。

a03.发送输入信号x(n)至待识别稀疏系统和自适应滤波系统。

a04.获取待识别稀疏系统的输出信号d(n)，d(n)为系统的期望输出信号：

d(n)＝x^t(n)h0+v(n)；

其中，h0为稀疏信道矢量，v(n)为零均值高斯白噪声。

a05.计算自适应滤波系统的输出信号y(n)，y(n)为系统的实际输出信号：

y(n)＝h^t(n)x(n)；

其中，h(n)为抽头权重向量。

a06.计算待识别稀疏系统与自适应滤波系统之间的输出误差e(n)：

e(n)＝d(n)-h^t(n)x(n)。

a07.构造代价函数l(n)：

其中，γ为稀疏约束加权值，是p范数，且0＜p＜1。

a08.计算代价函数l(n)的梯度

其中，p为大于0小于1的数值，|h(n)|^1-p为抽头权重向量|h(n)|的“1-p”次幂。

a09.构造动态步长μ(n)：

其中，μ为固定步长，vth是调节因子，且0＜vth＜1，为输出误差e(n)的平方根。

a10.计算自适应滤波器的抽头权重向量的更新函数：

h(n+1)＝h(n)+μ(n)e(n)x(n)-k(n)psgn[h(n)]/(ε+|h(n)|^1-p)；

k(n)＝γμ(n)；

其中，k(n)为正则化参数，ε为很小的正数。

a11.判断[e(n)-e(n-1)]＞ε0是否成立，若成立，则返回步骤a01；如果不等式不成立，则系统收敛，迭代结束，得到连续的稀疏系统的辨识值。

而系统收敛的判定方法为：

且

其中，λmax是输入信号的自相关矩阵的最大特征值，r是输入信号的自协方差矩阵。

本发明的优点在于：将误差的平方根引入到系统迭代过程时的步长控制中，并引入一个正常数以平衡系统的收敛速度和稳态误差，使得在迭代过程中能够实时调整系统的步长；引入的lρ范数自适应零吸引子可以加快小系数的收敛速度，消除大系数引起的系统辨识估计偏差。本发明具有较快的收敛速度和较高的估计精度、系统的稀疏度变化时，本发明的跟踪能力强，计算复杂度较低。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1是本发明的系统结构框图。

图2是系统信噪比(snr)为10db时，传统的lms、lms/f、vss-lms算法与本发明的性能仿真对比图。

图3是系统信噪比(snr)为20db时，传统的lms、lms/f、vss-lms算法与本发明的性能仿真对比图。

图4是信噪比为10db时，当信道的稀疏度产生突变的情况下，传统的lms、lms/f、vss-lms算法与本发明的跟踪情况的对比图。

图5是信噪比为20db时，当信道的稀疏度产生突变的情况下，传统的lms、lms/f、vss-lms算法与本发明的跟踪情况的对比图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有付出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本发明提出一种变步长lp范数lms算法的稀疏系统辨识方法。该方法在固定步长μ的基础上引入了误差函数的平方根来调节步长，利用平方根对小数(小于1)的放大作用对误差进行非线性放大，当误差较大时，该动态步长能提供一个较大值以促进均方误差的收敛；当误差较小时，动态步长能提供一个较小值以降低稳态误差，为了利用未知信道的的稀疏性，在代价函数中引入了lρ范数。

本实施例的一种变步长lp范数lms算法的稀疏系统辨识方法，系统信道长度设为128。自适应滤波系统如图1所示，方法的流程为：首先进行滤波器权值的初始化设置，设定初始化条件；训练输入信号并获取期望信号的值，将输入信号经过自适应滤波器获取输出信号、计算滤波器输出信号和期望信号的误差、构建变步长函数、代价函数中引入p范数稀疏惩罚项、求出权值更新向量，并判断其是否收敛，若不收敛，则返回对输出进行滤波处理，每次返回后需要更新输入信号x(n)；若收敛，则迭代结束。

具体步骤为：

a01.训练输入信号x(n)：

x(n)＝[x0x1…xn-1]^t；

其中，n为当前时刻，n为输入序列的信号长度，t为矩阵转置符号。

a02.初始化自适应滤波系统的抽头权重向量h(n)。

h(n)＝[h0h1…hn-1]^t；

h(n)是长为n的自适应滤波器信道估计值。

而对抽头权重向量h(n)进行初始化是将所有抽头权重向量h(n)置为零。

a03.发送输入信号x(n)至待识别稀疏系统和自适应滤波系统。

a04.获取待识别稀疏系统的输出信号d(n)，d(n)为系统的期望输出信号：

d(n)＝x^t(n)h0+v(n)；

其中，h0为稀疏信道矢量，是长为n(n＝128)的最优(维纳)解，仅由k个非零信道抽头组成，即稀疏度为k，分别取2和4；v(n)为零均值高斯白噪声且假设与x(n)统计独立。

a05.计算自适应滤波系统的输出信号y(n)，y(n)为系统的实际输出信号：

y(n)＝h^t(n)x(n)；

其中，h(n)为抽头权重向量。

a06.计算待识别稀疏系统与自适应滤波系统之间的输出误差e(n)：

e(n)＝d(n)-h^t(n)x(n)。

e(n)是第n次迭代得到的信道系数估计误差。

a07.构造代价函数l(n)，为了对稀疏系统进行辨识，在代价函数中引入能够表征稀疏特性的lp范数：

其中，γ为稀疏约束加权值，γ＝0.1；是p范数，且0＜p＜1。

当p趋近于零时，lρ范数近似于l0范数：

当p趋近于1时，lρ范数近似于l1范数：

a08.计算代价函数l(n)的梯度

其中，p为大于0小于1的数，|h(n)|^1-p为抽头权重向量|h(n)|的“1-p”次幂。

a09.构造动态步长μ(n)：

其中，μ为固定步长，在本实施例中μ＝0.012，vth是调节因子，vth＝0.8，为输出误差e(n)的平方根。

该步长利用平方根对小数(随着迭代的进行，误差一般都小于1)的放大挖掘作用，将误差放大，从而加速了误差的收敛。vth是一个常数调节因子，0＜vth＜1，通过和vth的调节可以使自适应滤波器在迭代初期误差较大时，步长μ(n)大于标准lms的步长μ，收敛速度加快；随着迭代进行，误差减小，μ(n)将小于标准lms的步长μ，稳态误差降低。

a10.计算自适应滤波器的抽头权重向量的更新函数：

h(n+1)＝h(n)+μ(n)e(n)x(n)-k(n)psgn[h(n)]/(ε+|h(n)|^1-p)；

k(n)＝γμ(n)；

其中，k(n)为正则化参数，ε＝0.05。

a11.判断[e(n)-e(n-1)]＞ε0是否成立，若成立，则返回步骤a01；如果不等式不成立，则系统收敛，迭代结束，得到连续的稀疏系统的辨识值。

而系统收敛的判定方法为：

且

其中，λmax是输入信号的自相关矩阵的最大特征值，r是输入信号的自协方差矩阵。

而且，本发明将传统算法与本发明算法进行了仿真实验，如图2和3，分别是系统信噪比(snr)为10db和20db时，传统的lms、lms/f、vss-lms算法与本发明的性能仿真对比图。

图4和5是信噪比为10db和20db时，当信道的稀疏度产生突变的情况下，传统的lms、lms/f、vss-lms算法与本发明的跟踪情况的对比图。

从对比图可知，本发明具有较快的收敛速度和较高的估计精度、系统的稀疏度变化时，本发明的跟踪能力强，计算复杂度较低。

上面所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。