趙 艷，繆善林

(1.通信信息控制和安全技術(shù)重點(diǎn)實(shí)驗室，浙江嘉興 314033;2.中國電子科技集團(tuán)第36研究所，浙江嘉興 314033)

0 引言

盲信號分離是在信號源和傳輸信道完全或部分未知的情況下，只利用傳感器陣列或天線陣列的觀測數(shù)據(jù)來分離、提取源信號的信號處理理論。盲信號分離技術(shù)是近年來信號處理領(lǐng)域的一個研究熱點(diǎn)，主要用于語音信號處理、圖像信號處理、通信信號處理及地震信號處理等。C.Jutten等學(xué)者在20世紀(jì)80年代早期就針對盲信號分離引入了獨(dú)立分量分析(ICA)技術(shù)，ICA技術(shù)包括極大似然估計方法［1］、最大熵方法、最小互信息方法、高階累積量方法［2］、主分量分析(PCA)算法［3，4］、自然梯度算法［5～7］及 FastICA 算法［8，9］等。一般研究數(shù)學(xué)模型中，都假設(shè)源信號的數(shù)目和觀測信號的數(shù)目相等，但實(shí)際電磁環(huán)境較為復(fù)雜，偵收到的源信號數(shù)目未知，即源信號數(shù)目可能比觀測信號數(shù)目少或者多，且觀測到的源信號頻譜可能存在重疊的現(xiàn)象，這就給現(xiàn)有算法的應(yīng)用帶來了一定的局限性。對于源信號數(shù)目多于觀測信號數(shù)目的盲信號分離，雖然已經(jīng)提出了一些解決辦法［10，11］，但在實(shí)際應(yīng)用中效果并不理想，還有待進(jìn)一步的研究。而當(dāng)源信號數(shù)目少于觀測信號數(shù)目時，常規(guī)方法認(rèn)為只要所有的源信號至少被重構(gòu)一次，即源信號可能被多次拷貝輸出，則實(shí)現(xiàn)了對源信號的有效盲分離，但是若要真正地提取源信號，還需要已知源信號數(shù)目及相關(guān)的信號處理技術(shù)才能實(shí)現(xiàn)。針對常規(guī)盲分離方法存在的這些問題，給出了基于廣義特征分解的盲信號分離方法，不需要知道信號數(shù)目就能夠?qū)υ葱盘枌?shí)施有效分離，所得的源信號不會被多次拷貝，并且在實(shí)際的應(yīng)用中取得了較好的效果。

1 基本原理

假設(shè)有M個天線單元接收N個相互獨(dú)立的入射信號(設(shè)為s1，s2，…，sN)，則天線陣列接收到的混疊信號可以表示為

式中，S=［s1，s2，…，sN］T;A 為 M ×N 列滿秩混疊矩陣;X=［x1，x2，…，xM］T為各個天線單元接收到的信號;［·］T為矩陣轉(zhuǎn)置。

為了得到各個輸入信號，需要找到一個去混疊矩陣W，即

式中，Y=［y1，y2，…，yN］T為 S 的估計矩陣。因此，Y中的任意一個元素y都可以認(rèn)為是N個輸入信號的線性混疊信號，即

式中，q=［q1，q2，…，qN］T是 N 維非零向量。

假設(shè)存在一個函數(shù)f使得f(si)(i=1，2，…，N)之間互不相關(guān)，且滿足式(4)和(5)

式中，cov(y)非零;cov(·)為求解協(xié)方差。則定義

由于S及f(S)都是由非相關(guān)信號組成的，其中 f(S)=［f(s1)，f(s2)，…，f(sN)］T，則存在以下關(guān)系

式中，diag［·］為對角陣。于是

即有

根據(jù)式(7)可以得到

另外，根據(jù)式(3)可以得到

根據(jù)前面的分析可以得到

式中，L(sm)=max{L(s1)，…，L(sN)}，則有

當(dāng)且僅當(dāng)y=asm時，式(13)等號才能成立。因此，可以根據(jù)下面的連續(xù)函數(shù)獲得信號sm的最優(yōu)求解

式中，y=wX(w為W的某一行向量)?；謴?fù)sm后即可恢復(fù)出其他信號，并最終獲得解混疊矩陣W。

2 基于廣義特征分解的求解方法

由于 L(s1)，L(s2)，…，L(sN)相互不相等，假設(shè)

根據(jù)前面的分析可以得到

即

式中，wGOP為w1的最優(yōu)求解。為了獲得wGOP，根據(jù)Q(w1)的梯度計算可得

令▽Q(w1)=0可得

則利用廣義特征分解即可獲得w1。

在實(shí)際計算中，需要將原始接收數(shù)據(jù)X進(jìn)行中心化及白化處理，中心化處理即將接收數(shù)據(jù)X變?yōu)榫禐?的信號為

式中，mean(·)為取均值。同時，為了將中心化處理后的信號進(jìn)行白化處理，需要對信號進(jìn)行特征值分解，即

式中，E表示特征向量;D表示特征值。

白化信號可以表示為

根據(jù)式(19)～(22)即可獲得混合前的信號。

3 實(shí)驗結(jié)果

實(shí)驗1:實(shí)驗采用的輸入信號樣式如式(23)表示，其混疊矩陣采用如式(24)所示?；旌锨暗?個原始信號如圖1所示，混合后的3路接收信號如圖2 所示，分離后的3個信號如圖3所示。

實(shí)驗2:某系統(tǒng)接收天線是由5單元天線組成的均勻圓陣，接收來自遠(yuǎn)處的2個信號，信號1的頻率為2 209.075 MHz和碼速率為9.6 ksps的 QPSK信號，信號2的頻率為2 209.045 MHz和碼速率為19.2 ksps的BPSK信號，則分離前后信號的頻譜如圖4所示。

圖4 分離前后信號頻譜圖

根據(jù)分離后的信號進(jìn)行參數(shù)估計，結(jié)果見表1，分離后的兩個信號解調(diào)結(jié)果，如圖5和圖6所示。

表1 分離信號參數(shù)估計

實(shí)驗3:采用實(shí)驗2的系統(tǒng)，接收來自遠(yuǎn)處的兩個信號，信號1的頻率為2 209.075 MHz和碼速率為24.8 ksps的QPSK信號，信號2的頻率為2 209.045 MHz和碼速率為19.2 ksps的BPSK信號，經(jīng)過多普勒頻率漂移之后，頻譜有所重疊，本文盲分離方法與常用的FastICA方法及自然梯度方法的分析處理結(jié)果，如圖7～圖9所示。其中，每幅圖中的第一個為觀測信號的頻譜圖，第二個至第六個為五個分離后的估計信號的頻譜圖。

通過實(shí)驗驗證可以看出，利用盲分離方法，輸出的估計信號中不存在源信號多次拷貝的問題，在源信號數(shù)目未知的情況下即可實(shí)現(xiàn)對頻率相隔很近的源信號進(jìn)行有效分離，與常規(guī)方法相比具有更好的分離效果。

4 結(jié)語

針對當(dāng)前復(fù)雜的電磁環(huán)境，給出了一種基于廣義特征值分解的盲信號分離方法，該方法無需進(jìn)行迭代運(yùn)算即可實(shí)現(xiàn)盲信號分離，從而降低了計算的復(fù)雜度，且對信號的適應(yīng)性具有較好的效果，通過實(shí)驗驗證了該方法的有效性。

［1］LIN S L，TUNG P C.Application of Modified ICA to Secure Communication in Chaotic Systems［C］//Lecture Notes In Computer Science，2007(4707):431-444.

［2］TAOUFIKI M，ADIB A，ABOUTAJDINE D.Blind Separation of Any Source Distributions via High-Order Statistics［J］.Signal Processing，2007，87(8):1 882-1 889.

［3］呂齊，朱孝龍，張賢達(dá).卷積盲信號分離的非線性PCA方法［J］.清華大學(xué)學(xué)報(自然科學(xué)版)，2006，46(4):492-494.

［4］馬麗艷，李宏偉.一種基于非線性PCA的卷積混合盲源分離算法［J］.電子學(xué)報，2008，36(5):1 009-1 012.

［5］付衛(wèi)紅，楊小牛，曾興雯，等.適用于通信偵察的信號盲分離算法［J］.南京理工大學(xué)學(xué)報:自然科學(xué)版，2008，32(2):189-194.

［6］CICHOCKI A，UNBEHAUEN R.Robust Neural Networks with On-line Learning for Blind Identification and Blind Separation of Sources［J］.IEEE Trans.Circuits and Systems I:Fundamentals Theory and Application，1996，43(11):894-906.

［7］AMARI S，CICHOCKI A，YANG H H.A New Learning Algorithm for Blind Signal Separation［C］//Advances in Neural Information Processing Systems，Canbrige，Ma:MIT Press，1995，8:757-763.

［8］郭武，朱長仁，王潤生.一種改進(jìn)的FastICA算法及其應(yīng)用［J］.計算機(jī)應(yīng)用，2008，28(4):960-962.

［9］賈銀潔，許鵬飛.基于FastICA的混合音頻信號盲分離［J］.信息與電子工程，2009，7(4):321-325.

［10］T-W LEE.Independent Component Analysis-Theory and Application［M］.Holland:Kluwer，1998.

［11］GIROLAMI M.Advances in Independent Component Analysis［M］.Berlin:SpringerVerlag，2000.