[发明专利]基于改进时频单源区的欠定混合矩阵盲估计方法

摘要

本发明公开了一种基于改进时频单源区的欠定混合矩阵盲估计方法。本发明首先获取观测信号矩阵，对观测信号矩阵进行伽柏Gabor变换得到时频信号矩阵，然后构建时频单源区信号矩阵，剔除时频单源区信号矩阵中的低能量值的列向量，最终对时频单源区信号矩阵进行聚类得到欠定混合矩阵。本发明克服了现有技术存在的混合矩阵估计误差较高和容易受到源信号稀疏性的影响的缺点，使得本发明具有能保持较低混合矩阵估计误差的同时降低对源信号稀疏性的要求。

主权项

一种基于改进时频单源区的欠定混合矩阵盲估计方法，包括如下步骤：(1)将采集到的通信信号存入观测信号矩阵中；(2)对观测信号矩阵进行伽柏Gabor变换，得到时频信号矩阵；(3)构建时频单源区信号矩阵：(3a)按照下式，计算每一个时频信号矩阵中列向量的角度系数：ρ=|(R{X(t,k)})TI{X(t,k)}||R{X(t,k)}||||I{X(t,k)}|||]]>其中，ρ表示时频信号矩阵中列向量的角度系数，R{·}表示取实部操作，X(t,k)表示时频信号矩阵中采集时刻t和频率点k对应的列向量，(·)T表示转置操作，I{·}表示取虚部操作，||·||表示取2范数操作，|·|表示取绝对值操作。(3b)判断每一个时频信号矩阵中列向量的角度系数是否大于0.95，若是，则执行步骤(3c)，否则，则执行步骤(3d)；(3c)将时频信号矩阵中角度系数大于0.95的列向量存入时频单源区信号矩阵；(3d)剔除时频信号矩阵中角度系数小于0.95的列向量；(4)剔除时频单源区信号矩阵中的低能量值的列向量：(4a)按照下式，计算门限值：ε＝0.2max(||X(t,k)||2)其中，ε表示门限值，max(·)表示取最大值操作；(4b)按照下式，计算时频单源区信号矩阵中列向量的能量值：E=||X^(t,k)||2]]>其中，E表示时频单源区信号矩阵中列向量的能量值，表示时频单源区信号矩阵中采集时刻t和频率点k对应的列向量；(4c)判断时频单源区信号矩阵中列向量的能量值是否大于门限值，若是，则执行步骤(4d)，否则，则执行步骤(4e)；(4d)保留时频单源区信号矩阵中能量值大于门限值的列向量；(4e)剔除时频单源区信号矩阵中能量值大于门限值列向量；(5)对时频单源区信号矩阵进行聚类，得到欠定混合矩阵。