1.一種基于韋羅內(nèi)塞映射的欠定盲分離混合矩陣估計方法,包括如下步驟:
(1)構(gòu)建觀測信號矩陣:
(1a)構(gòu)建一個行數(shù)與接收端天線個數(shù)相等、且元素全為0的觀測信號矩陣;
(1b)將接收端天線采集到的通信信號,按列存入到觀測信號矩陣中;
(2)對觀測信號矩陣進行韋羅內(nèi)塞映射:
(2a)按照下式,計算韋羅內(nèi)塞映射向量的維數(shù):
其中,m表示韋羅內(nèi)塞映射向量的維數(shù),!表示階乘操作,i表示韋羅內(nèi)塞映射的階數(shù),M表示接收端天線的總數(shù);
(2b)按照下式,計算每一個韋羅內(nèi)塞映射向量:
Vj=[v(x1),v(x2),…,v(xT)]
其中,Vj表示第j個韋羅內(nèi)塞映射向量,v(·)表示取韋羅內(nèi)塞系數(shù)操作,x1表示觀測信號矩陣的第1個列向量,x2表示觀測信號矩陣的第2個列向量,xT表示觀測信號矩陣的最后一個列向量;
(2c)將所有的韋羅內(nèi)塞映射向量按列組成韋羅內(nèi)塞映射矩陣;
(3)構(gòu)造韋羅內(nèi)塞映射的法向量:
(3a)對韋羅內(nèi)塞映射矩陣進行奇異值分解,得到韋羅內(nèi)塞映射的左奇異矩陣;
(3b)取出韋羅內(nèi)塞映射的左奇異矩陣中最后一個列向量,作為韋羅內(nèi)塞映射的法向量;
(4)按照下式,計算韋羅內(nèi)塞映射的一階導(dǎo)數(shù):
Dp=nonzero(cL⊙P)v
其中,Dp表示韋羅內(nèi)塞映射的一階導(dǎo)數(shù),nonzero(·)表示剔除向量中的零元素操作,c表示韋羅內(nèi)塞映射的法向量,L表示元素全為1的行向量,⊙表示點乘操作,P表示韋羅內(nèi)塞映射的冪次,v表示韋羅內(nèi)塞映射矩陣中的最后一個列向量;
(5)對韋羅內(nèi)塞映射的一階導(dǎo)數(shù)進行譜圖劃分,得到子空間矩陣;
(6)對子空間矩陣的所有交線列向量進行聚類,得到欠定混合矩陣。
2.根據(jù)權(quán)利要求1所述的基于韋羅內(nèi)塞映射的欠定盲分離混合矩陣估計方法,其特征在于:步驟(5)中所述對韋羅內(nèi)塞映射的一階導(dǎo)數(shù)進行譜圖劃分的具體步驟如下:
第1步,對韋羅內(nèi)塞映射的一階導(dǎo)數(shù)矩陣的拉普拉斯系數(shù)進行特征譜分析,得到特征值矩陣;
第2步,對特征值矩陣的所有行向量進行聚類,得到子空間矩陣。
3.根據(jù)權(quán)利要求1所述的基于韋羅內(nèi)塞映射的欠定盲分離混合矩陣估計方法,其特征在于:步驟(6)中所述對子空間矩陣的所有列向量進行聚類的具體步驟如下:
第1步,從子空間矩陣的所有列向量中選取第一個分量為負(fù)數(shù)的列向量,對該列向量的所有元素乘以-1得到翻轉(zhuǎn)列向量;
第2步,將翻轉(zhuǎn)列向量和未被選取的列向量按列組成翻轉(zhuǎn)矩陣;
第3步,對翻轉(zhuǎn)矩陣進行歸一化處理,得到歸一化信號矩陣;
第4步,利用K-均值聚類方法,對歸一化信號矩陣中的所有列向量進行聚類,得到所有的聚類中心;
第5步,將所有的聚類中心按列組成欠定混合矩陣。