[发明专利]基于梯度下降法和牛顿法的欠定盲源分离源信号恢复方法

摘要

本发明公开了一种基于梯度下降法和牛顿法的欠定盲源分离源信号恢复方法。本发明首先获取观测信号矩阵，然后对观测信号矩阵中的所有列向量进行聚类得到混合矩阵，根据观测信号矩阵和混合矩阵计算待恢复的源信号列向量，依次利用梯度下降法、牛顿法和投影法更新待恢复的源信号列向量，最终获得恢复的源信号。本发明克服了现有技术存在的源信号恢复精度易受噪声误差影响和计算复杂度偏高的缺点，使得本发明具有能保持较高恢复精度的同时快速的恢复源信号的优点。

主权项

一种基于梯度下降法和牛顿法的欠定盲源分离源信号恢复方法，包括如下步骤：(1)将采集到的通信信号存入观测信号矩阵；(2)对观测信号矩阵进行聚类，得到混合矩阵；(3)将欠定盲源分离的迭代次数初始化为1；(4)按照下式，计算待恢复源信号列向量：s^=AT(AAT)-1x(l)]]>其中，表示待恢复源信号列向量，A表示混合矩阵，T表示转置操作，(·)‑1表示求逆操作，x(l)表示观测信号矩阵中的第l个列向量，l为欠定盲源分离的迭代次数；(5)赋值目标参数：(5a)搜索待恢复源信号列向量中的最大分量；(5b)将搜索到的最大分量的值赋给目标参数；(6)利用梯度下降法更新待恢复源信号列向量：(6a)将梯度下降法的迭代次数初始化为1；(6b)计算梯度下降法的增量；(6c)将待恢复源信号列向量加上梯度下降法的增量，得到利用梯度下降法更新后的待恢复源信号列向量(6d)将梯度下降法的迭代次数加1，得到梯度下降法的当前迭代次数；(6e)判断梯度下降法的当前迭代次数是否小于等于M，若是，执行步骤(6b)，否则，执行步骤(7)，其中，M表示待恢复源信号列向量的个数；(7)利用牛顿法更新待恢复源信号列向量：(7a)计算牛顿方向值；(7b)将利用梯度下降法更新后的待恢复源信号列向量加上牛顿方向值，得到利用牛顿法更新后的待恢复源信号列向量(8)利用投影法更新待恢复源信号列向量：(8a)将投影法的迭代次数初始化为1；(8b)计算投影法的增量；(8c)将利用牛顿法更新后的待恢复源信号列向量加上投影法的增量，得到利用投影法更新后的待恢复源信号列向量(8d)将投影法的迭代次数加1，得到投影法的当前迭代次数；(8e)判断投影法的当前迭代次数是否小于等于M，若是，则执行步骤(8b)，否则，执行步骤(9)，其中，M表示待恢复的源信号列向量的个数；(9)将目标参数乘以0.8，得到更新后的目标参数；(10)判断更新后的目标参数是否大于10‑5，若是，则执行步骤(6)，否则，执行步骤(11)；(11)存储利用投影法更新后的待恢复源信号列向量；(12)将欠定盲源分离的迭代次数加1，得到欠定盲源分离的当前迭代次数；(13)判断欠定盲源分离的当前迭代次数是否小于等于观测信号矩阵的列数，若是，执行步骤(4)，否则，执行步骤(14)；(14)获得恢复源信号矩阵：将欠定盲源分离迭代得到的所有待恢复源信号列向量组成恢复源信号矩阵。