劉治軍 梁宗闖 邱樂德 周業(yè)軍 齊維孔

(中國空間技術(shù)研究院通信衛(wèi)星事業(yè)部，北京 100094)

?

一種同頻非合作航天器通信信號的盲源分離算法

劉治軍 梁宗闖 邱樂德 周業(yè)軍 齊維孔

(中國空間技術(shù)研究院通信衛(wèi)星事業(yè)部，北京 100094)

針對頻譜混疊的同頻非合作航天器通信信號分離問題，研究了盲源分離技術(shù)，提出一種以立方函數(shù)作為非線性函數(shù)的盲源分離改進(jìn)算法。在介紹盲源分離模型的基礎(chǔ)上，依據(jù)通信信號的亞高斯性，對快速獨立分量分析(FastICA)算法的非線性函數(shù)作適應(yīng)性改進(jìn)，提升了算法的性能。選用8FSK、BPSK、QPSK和頻譜混疊的通信信號與單音干擾信號進(jìn)行混合，仿真并對比分析了4種應(yīng)用不同非線性函數(shù)的FastICA算法。500次的仿真結(jié)果表明：應(yīng)用立方函數(shù)的FastICA算法的分離成功率達(dá)到100%，平均迭代次數(shù)約為25，性能指數(shù)均值為0.115 7，具有更優(yōu)的效能，可很好地解決非合作航天器通信信號的分離問題。

非合作航天器通信信號；立方函數(shù)；盲源分離；頻譜混疊

1 引言

地球同步軌道上衛(wèi)星數(shù)量的增多，以及各種新型通信技術(shù)(如成對載波多址技術(shù))在航天器上的應(yīng)用，使得實際環(huán)境中無線電信號日益密集，自然干擾和人為干擾日益嚴(yán)重。衛(wèi)星數(shù)量日益增多導(dǎo)致衛(wèi)星間距越來越小，衛(wèi)星信號會相互干擾；各種新型通信技術(shù)的應(yīng)用使得下行信號本身就是混合信號，從而使得通信接收、頻譜監(jiān)視等應(yīng)用中遇到的混合信號越來越多。這類信號一般在時域、頻域上易發(fā)生混疊，無法采用傳統(tǒng)的信號分離方法進(jìn)行處理。另外，此類信號多是非合作的，即不知道信號的帶寬、碼速率和載波頻率等基本信息，進(jìn)一步加大了處理難度，因而須要研究同頻非合作航天器通信信號的分離技術(shù)，而盲源分離技術(shù)是一個很好的選擇。

盲源分離技術(shù)的研究起源于20世紀(jì)80年代[1-2]，國內(nèi)的研究始于20世紀(jì)90年代[3]。文獻(xiàn)[4-5]中對盲源分離的理論、方法和應(yīng)用進(jìn)行了概括，文獻(xiàn)[6-12]中將此理論與相應(yīng)的研究領(lǐng)域結(jié)合進(jìn)行了研究。其中：文獻(xiàn)[8-11]中探討了盲源分離在空間電子偵察中的應(yīng)用情況，并進(jìn)行了算法仿真分析。文獻(xiàn)[12]中研究了如何用此技術(shù)解決鄰星干擾抑制與抵消問題，以及在星載船舶自動識別系統(tǒng)(Automatic Identification System，AIS)中的突發(fā)混合信號盲分離問題，并取得了較好的仿真結(jié)果。在算法原理上，文獻(xiàn)[8]中使用峭度作為代價函數(shù)對非高斯性進(jìn)行度量，得到了較好的分離效果；但在實際應(yīng)用中，由于峭度的值只能從測量樣本中估計，峭度方法可能對野值極其敏感，并不是非高斯性的一個魯棒度量[7]。文獻(xiàn)[9-10]中選擇了魯棒性較好的負(fù)熵作為代價函數(shù)，并分別使用并行與串行的FastICA算法，實現(xiàn)了對同頻非合作信號的分離。串行FastICA算法與源數(shù)目一次方成正比，并行算法與源數(shù)目的二次方成正比，因而也不存在不適用高維空間的問題。不過，在非線性函數(shù)的選取上，文獻(xiàn)[9-10]中并未具體針對通信信號進(jìn)行討論，所使用的非線性函數(shù)為雙曲正切函數(shù)，但未能充分考慮通信信號的亞高斯性。

本文在已有研究的基礎(chǔ)上，將盲源分離技術(shù)應(yīng)用于分離頻譜混疊的非合作航天器通信信號。通過結(jié)合通信信號的亞高斯性，對基于負(fù)熵的串行FastICA算法作適應(yīng)性改進(jìn)；對比仿真了使用4種不同非線性函數(shù)的分離算法，驗證了應(yīng)用立方函數(shù)的盲源分離算法的優(yōu)勢。

2 盲源分離模型

最簡單的盲源分離模型是線性混合模型，每個傳感器信號xi(k)都是源信號sj(k)的線性組合。線性瞬時混合模型表示為

(1)

式中：N為源個數(shù)，M為傳感器個數(shù)，簡化模型常常假設(shè)M=N；k為觀測時刻，對于數(shù)字信號來說，表示第k個觀測值；xi(k)為第i個傳感器在第k個時刻接收到的值；sj(k)為第j個源在第k個時刻的觀測值；aij為源信號線性混合的混合系數(shù)。

式(1)寫成矩陣形式為

(2)

式中：x(k)為傳感器接收信號的值，即觀測信號；s(k)為源信號；A為M×N階的混合矩陣，表征源信號的線性混合方式。

盲源分離的任務(wù)是在s和A未知的情況下，僅通過對x的處理，得到分離矩陣W。為了書寫方便，求出源信號s的估計y，對瞬時混合模型的分離可表示為

(3)

當(dāng)W和A的乘積是一個每行和每列只有一個非零元素的矩陣時，就達(dá)到了恢復(fù)信號的目的。在實際應(yīng)用中，對于觀測信號首先需要進(jìn)行預(yù)處理(即乘以預(yù)處理矩陣V)得到信號z，之后再進(jìn)行分離(即乘以分離矩陣W)。

針對含噪聲的情況，若噪聲是加性噪聲，則可以當(dāng)作一個獨立源進(jìn)行分離；若噪聲是非加性的，含在獨立源中，模型中的s就是源和噪聲的混合結(jié)果，分離后得到的y也是帶噪聲的，只要合理設(shè)計分布式天線，即可保證混合矩陣的條件數(shù)較小，噪聲便不會被放大[12]。

圖1是盲源分離系統(tǒng)的示意框圖。前面部分是未知源信號的混合過程；后半部分是盲源信號的分離過程。在盲源信號處理中，一般先要對混合信號進(jìn)行預(yù)處理。預(yù)處理包括中心化和白化：中心化處理是使信號的均值為零；白化處理可去除各觀測信號之間的相關(guān)性，從而簡化后續(xù)獨立分量的提取過程。通常情況下，對數(shù)據(jù)進(jìn)行白化處理[8]以使算法的收斂性較好，工作量少，效率高。

圖1 盲源分離系統(tǒng)模型Fig.1 System model for blind source separation

3 應(yīng)用立方函數(shù)的盲源分離改進(jìn)算法

本文在FastICA算法[7]原理的基礎(chǔ)上，結(jié)合通信信號的亞高斯性，對所選取的非線性函數(shù)進(jìn)行適應(yīng)性改進(jìn)。FastICA算法又稱固定點(Fixed-point)算法，有基于峭度、基于負(fù)熵等形式。峭度和負(fù)熵均可作為非高斯性的度量，但峭度不是一個魯棒的度量，樣本中某個不準(zhǔn)確的值會使峭度發(fā)生非常大的變化；而負(fù)熵則能很好地克服這一缺點[8]。因此，本文采用的是應(yīng)用負(fù)熵最大的FastICA算法。

3.1 FastICA算法基本原理

應(yīng)用負(fù)熵最大的FastICA算法以負(fù)熵最大作為一個搜尋方向，可以實現(xiàn)順序提取獨立源的目的，充分體現(xiàn)探查性投影追蹤(Exploratory Projection Pursuit，EPP)這種傳統(tǒng)線性變換的思想[12]。此外，該 算法采用了定點迭代的優(yōu)化方法，使得收斂更加快速、穩(wěn)健。在該算法中，負(fù)熵可以近似為[7]

(4)

式中：N(Y)為反映Y偏離高斯隨機向量的程度；E[]為均值運算；g()為非線性函數(shù)，其選取規(guī)則詳見第3.2節(jié)；Y為某一隨機向量，YGauss是與Y具有相同協(xié)方差矩陣的高斯隨機向量。

經(jīng)過一系列的數(shù)學(xué)推導(dǎo)后，可以得到FastICA算法如下。

(5)

(6)

式中：W*為W的新值；E{}表示取數(shù)據(jù)的期望；X為觀測矩陣；g′()表示非線性函數(shù)g()的導(dǎo)數(shù)。

3.2 非線性函數(shù)的改進(jìn)

非線性函數(shù)g(y)經(jīng)常取為g1(y)=tanh(a1y)，適用于亞高斯和超高斯信號并存的情況，通用性較強。這里，1≤a1≤2，通常取a1=1。文獻(xiàn)[9-10]中使用此非線性函數(shù)成功分離了幾種頻譜混疊的非合作信號。由于通信信號常常是亞高斯信號，因而選擇專用性更強的非線性函數(shù)會獲得更好的分離結(jié)果。本文列舉了幾個在其他領(lǐng)域使用過的非線性函數(shù)，如表1所示。后文中將使用應(yīng)用這4種非線性函數(shù)的FastICA算法，對頻譜混疊的通信信號進(jìn)行仿真分析，以尋求適用于同頻非合作通信信號的盲源分離技術(shù)。

表1 非線性函數(shù)及其適用范圍

信號是亞高斯的還是超高斯的，可以根據(jù)峭度來判斷。峭度是表征數(shù)據(jù)在均值兩側(cè)集中程度的參數(shù)，定義為

(7)

峭度的值可正可負(fù)，具有負(fù)峭度的隨機變量稱為亞高斯的，峭度為正的隨機變量稱為超高斯的。數(shù)字通信信號的常用調(diào)制方式有BPSK、QPSK和FSK等幾種，一般為亞高斯信號，因而理論上選擇非線性函數(shù)g3和g4所得分離結(jié)果更好。

3.3 改進(jìn)算法

改進(jìn)算法的基本步驟如下[12]。

(1)預(yù)處理：首先對觀測矩陣X進(jìn)行中心化，使它的均值為0；然后對去中心化的數(shù)據(jù)進(jìn)行白化處理，得到新的矩陣Z，X→Z。

(2)設(shè)迭代次數(shù)n=0，迭代精度ε=0.000 01，最大迭代次數(shù)nmax=10 000，待估計分量的個數(shù)m=m0，當(dāng)前提取的源信號序號p=1。

(3)設(shè)定初始迭代矢量(隨機的)Wp=Wp0。

(4)保存上次迭代結(jié)果Wlast，p=Wp；設(shè)定迭代方程Wp=E{Zg(WpTZ)}-E{g′(WpTZ)}W,非線性函數(shù)選用立方函數(shù)。

(6)將數(shù)據(jù)歸一化處理：令Wp=Wp/‖Wp‖。

(7)判斷迭代次數(shù)是否超過限度nmax：若n>nmax，結(jié)束；否則，迭代次數(shù)加1，n=n+1。

(8)若Wp不收斂，即不滿足|Wlast，p-Wp|<ε，返回第4步；否則，p=p+1。

(9)如果已分離信號個數(shù)p大于源個數(shù)m，返回第3步；否則，結(jié)束，最終獲得的分離矩陣為W=[W1W2…Wm]。

3.4 評價指標(biāo)

對盲源分離算法的評價主要包括兩方面：一是算法的收斂速度，二是分離后的信號與源信號的逼近程度。本文采用迭代次數(shù)對收斂速度進(jìn)行衡量；采用相似性系數(shù)與性能指數(shù)兩個指標(biāo)衡量分離結(jié)果的逼近程度。相似性系數(shù)便于直觀理解試驗分離后的結(jié)果，但不便于大量數(shù)據(jù)統(tǒng)計分析；性能指數(shù)則相反。因此，本文先利用相似性系數(shù)進(jìn)行單次試驗，直觀說明分離結(jié)果，再利用性能指數(shù)對大量仿真結(jié)果進(jìn)行統(tǒng)計分析。

1)相似系數(shù)

相似系數(shù)的定義[9]為

(8)

式中：T為觀測樣值個數(shù)。

相似系數(shù)用來度量兩個信號之間的相似程度，它抵消了分離結(jié)果在幅值上存在的差異，從而避免了幅值不確定性的影響。當(dāng)由相似系數(shù)構(gòu)成的相似系數(shù)矩陣的每行和每列僅有一個元素接近于1，則可以認(rèn)為算法分離效果較為理想。

2)性能指數(shù)

性能指數(shù)(PerformanceIndex)的定義[13-14]為

(9)

式中：fij為全局傳輸矩陣F=WA的元素；maxj|Fij|表示全局傳輸矩陣的第i行元素絕對值中的最大值；maxj|fji|表示第i列元素絕對值中的最大值。

若源信號s的估計y與源信號s的波形完全相同，則PI=0，并且式(9)是全局傳輸矩陣為方陣的情況，在其他情況下須對其進(jìn)行修正。

4 仿真算例

4.1 通信信號仿真

在航天數(shù)字通信中，常用調(diào)制方式有FSK、BPSK和QPSK等，因此，本文將8FSK、BPSK、QPSK和單音正弦波4種信號進(jìn)行混合，混合矩陣是一個隨機產(chǎn)生的矩陣，各元素值在(0,1)之內(nèi)，條件數(shù)小于3。設(shè)定仿真條件見表2。為使觀察更形象直觀，在不影響試驗結(jié)果的情況下，將8FSK的調(diào)制偏移量設(shè)置略高。用于仿真的4種信號的時域波形如圖2所示。4種信號的頻譜圖見圖3，可以直觀看出，幾種信號的頻率較為接近，很難用常規(guī)的頻譜分離方法進(jìn)行分離。圖4為4種信號的混合波形，即如圖1所示的觀測向量x(t)，由源信號乘以混合矩陣A所得。后面將以此作為輸入，采用本文的分離算法對x(t)進(jìn)行分離。

表2 仿真條件

圖2 通信信號時域波形Fig.2 Waveforms of communication signals

圖3 通信信號頻譜波形Fig.3 Waveforms of communication signal spectrums

圖4 4種通信信號混合后的波形Fig.4 Waveforms of four mixed communication signals

保持4種信號的信噪比為15 dB，仿真1000次，計算各自峭度并進(jìn)行統(tǒng)計，其結(jié)果如表3所示。3種常用通信信號的峭度范圍均為(-0.5，-0.3)，且都為負(fù)值。亞高斯信號的峭度為負(fù)值，統(tǒng)計結(jié)果驗證了通信信號的亞高斯性。

表3 通信信號的峭度統(tǒng)計

4.2 分離算法仿真結(jié)果與分析

使用相同的混合矩陣對原始信號進(jìn)行混合，對混合信號按第3.3節(jié)的分離算法進(jìn)行分離。4種不同非線性函數(shù)得到4個FastICA算法，使用4種算法分別對混合信號進(jìn)行分離仿真，收斂精度為0.000 01，最大迭代次數(shù)10 000。

根據(jù)不同算法得到的4個相似系數(shù)矩陣依次為

從相似系數(shù)矩陣還可以得到分離后的信號與源信號之間的對應(yīng)關(guān)系。以第1個矩陣為例，其對應(yīng)關(guān)系為s1←y2，s2←y1，s3←y3，s4←y4。分離后的信號與源信號之間的相關(guān)系數(shù)大于0.97，說明分離效果比較理想。

仿真試驗將3種不同調(diào)制方式、中心頻率相同的通信信號與1種在其帶寬內(nèi)的單音干擾信號進(jìn)行混合，目的是仿真真實通信過程中相同中心頻率的不同調(diào)制方式信號混合，并且?guī)拑?nèi)存在單音信號干擾的情況。一般，對于通信信號，信號帶寬在中心頻率的3%以內(nèi)，因而在通信過程中，單音信號會對相同中心頻率的通信信號會產(chǎn)生干擾。幾種相同中心頻率的通信信號混合在一起，是難以通過頻域濾波的方法進(jìn)行分離的，因此必須采用新型的信號分離技術(shù)。

為進(jìn)一步對比分離結(jié)果，分別使用4種非線性函數(shù)進(jìn)行500次分離仿真試驗，結(jié)果如表4所示。4個函數(shù)的結(jié)果的成功率在95%以上，均在可接受范圍內(nèi)。非線性函數(shù)g1與g2的迭代次數(shù)至少60次，而g3與g4的迭代次數(shù)不到30次，并且性能指數(shù)更小，總體而言，使用g3與g4的分離結(jié)果更優(yōu)。

表4 4種非線性函數(shù)仿真所得統(tǒng)計結(jié)果

進(jìn)一步比較各非線性函數(shù)的效能，使用非線性函數(shù)g3所得的分離結(jié)果分別與使用g1，g2，g4所得的結(jié)果進(jìn)行對比，見圖5～7(分布的局部區(qū)域的放大圖，個別分離失敗的點未列入其中)。圖中每個散點的意義是，每次仿真完成后，該次仿真所需要的迭代次數(shù)以及此次仿真結(jié)果所算出的分離性能指數(shù)，因而點越靠近右上方，表明分離結(jié)果越差。在圖5與圖6中，從迭代次數(shù)來看，使用g3要優(yōu)于使用g1與g2，而對于性能指數(shù)，g1與g2的整體分布略大于g3。圖7中，代表g4的散點的整體分布居于代表g3的散點的右上方，因而使用g3的分離結(jié)果略優(yōu)于使用g4的。綜合以上結(jié)果，非線性函數(shù)g3更適合通信信號的盲源分離。

圖5 使用非線性函數(shù)g1與g3分離所得指標(biāo)的分布Fig.5 Distributing of indexes using nonlinear function g1 and g3

圖6 使用非線性函數(shù)g2與g3分離所得指標(biāo)的分布Fig.6 Distributing of indexs using nonlinear function g2 and g3

圖7 使用非線性函數(shù)g3與g4分離所得指標(biāo)的分布Fig.7 Distributing of indexes using nonlinear function g3 and g4

5 結(jié)束語

選擇何種非線性函數(shù)是盲源分離技術(shù)的關(guān)鍵問題。本文針對這一問題的研究結(jié)果表明：在信噪比與混合矩陣條件數(shù)滿足一定條件的情況下，盲源分離技術(shù)可以對頻譜混疊的非合作通信信號進(jìn)行分離；可以完成對通信帶寬內(nèi)單音干擾信號的分離；針對非合作通信信號的盲源分離，采用應(yīng)用立方函數(shù)的FastICA算法優(yōu)于采用雙曲正切函數(shù)的分離算法。本文提出的改進(jìn)算法，可將顯著提高航天器通信信號的抗干擾能力，也能提升針對非合作目標(biāo)信號接收的分離能力。由于本文的算法是在非欠定條件下的分離，而在欠定條件下的盲源分離，還要進(jìn)一步深入研究。

References)

[1]Jutten C，Herault J.Blind separation of sources.Part Ⅰ：an adaptive algorithm based on neuromimetic architecture[J].Signal Processing，1991，24(1)：1-10

[2]Common P，Jutten C，Herault J.Blind separation of sources，Part Ⅱ：problems statement[J].Signal Processing，1991，24(1)：11-20

[3]王衛(wèi)華.盲源分離算法及應(yīng)用研究[D].哈爾濱：哈爾濱工程大學(xué)，2009

Wang Weihua.Application and algorithm of the blind source separation[D].Harbin：Harbin Engineering University，2009 (in Chinese)

[4]張賢達(dá).盲信號分離[J].電子學(xué)報，2001，29(12A)：1766-1771

Zhang Xianda.Blind signals separation[J].Acta Electronica Sinica，2001，29(12A)：1766-1771 (in Chinese)

[5]張賢達(dá)，保錚.通信信號處理[M].北京：清華大學(xué)出版社，2000

Zhang Xianda，Bao Zheng.Communication signals processing[M].Beijing：Tsinghua University Press，2000 (in Chinese)

[6]Cardoso J F.Infomax and maximum likelihood for blind source separation [J].IEEE Signal Proc. Letters，1997,4(4):112-114

[7]Hyv?rinen A.獨立成分分析[M].周宗潭，董國華，徐昕，等譯.北京：電子工業(yè)出版社，2007

Hyv?rinen A.Independent component analysis[M].Zhou Zongtan，Dong Guohua，Xu Xin，et al，translated.Beijing：Publishing House of Electronics Industry，2007 (in Chinese)

[8]李莉.基于峭度的盲分離在通信信號盲偵察中的應(yīng)用[J].通信技術(shù)，2010，43(4)：133-135，138

Li Li.Application of blind source separation based on kurtosis in communication signal blind reconnaissance [J].Communications Technology，2010，43(4)：133-135，138 (in Chinese)

[9]周治宇，陳豪.空間同頻電子偵察信號的盲分離[J].中國空間科學(xué)技術(shù)，2009，29(3)： 10-15，22

Zhou Zhiyu，Chen Hao.Blind separation of space radio reconnaissance signals with same frequency [J].Chinese Space Science and Technology，2009，29(3)： 10-15，22 (in Chinese)

[10]周治宇.基于盲信號分離的同頻信號的串行分離技術(shù)[J].信息與電子工程，2009，7(4)：308-313

Zhou Zhiyu.Serial separation of signals with same frequency based on blind signal separation[J].Information and Electronic Engineering，2009，7(4)：308-313 (in Chinese)

[11]周治宇，陳豪.空間電子偵察信號的盲分離[C]//衛(wèi)星有效載荷技術(shù)學(xué)術(shù)會議論文集.西安：中國空間技術(shù)研究院西安分院，2009：92-98

Zhou Zhiyu，Chen Hao.Blind separation of space radio reconnaissance signals [C]//Proceedings of Satellite Payload Technology Academic Conference.Xi’an：CAST-Xi’an Academy of Space Electronic Information Technology，2009：92-98 (in Chinese)

[12]萬堅，涂世龍，廖燦輝，等.通信混合信號盲分離理論與技術(shù)[M].北京：國防工業(yè)出版社，2012

Wan Jian，Tu Shilong，Liao Canhui，et al.Theory and technology on blind source separation of communication signals [M].Beijing：National Defence Industry Press，2012 (in Chinese)

[13]吳其澤.語音信號的盲源分離算法研究及應(yīng)用[D].廣州：中山大學(xué)，2010

Wu Qize.Researches and application of blind speech signals separation [D].Guangzhou：Sun Yatsen University，2010 (in Chinese)

[14]李舜酩.振動信號的盲源分離技術(shù)及應(yīng)用[M].北京：航空工業(yè)出版社，2011

Li Shunming.Technique and applications of blind vibration signal separation [M].Beijing：Aviation Industry Press，2011 (in Chinese)

(編輯：夏光)

Algorithm of Blind Source Separation for Non-cooperative Spacecraft Communication Signals with Same Carrier Frequencies

LIU Zhijun LIANG Zongchuang QIU Lede ZHOU Yejun QI Weikong

(Institute of Telecommunication Satellite，China Academy of Space Technology，Beijing 100094，China)

To complete the separation of non-cooperative spacecraft communication signal with spectrum aliased，the research of BSS (blind source separation) is carried out.Based on the cubic function as a nonlinear one，the algorithm of BSS is improved.The separation model of BSS is described.The nonlinear function of FastICA (fast independent component analysis) is modified to improve the capability of the algorithm because communication signals are sub-Gaussian signals.Three kinds of common communication signal are mixed with a kind of sine wave interference signal.Several simulations of different algorithms are compared based on different nonlinear functions.The results of 500 simulations show that FastICA algorithm based on cubic function enhances the separation success rate to 100%，and the average number of iteration times is 25,the performance index is 0.1157.The cubic function makes the algorithm get the best effect and completes the separation of non-cooperative spacecraft communication signal well.

non-cooperative spacecraft communication signal；cubic function；blind source separation；aliased spectrum

2014-04-08；

2014-05-08

國家高技術(shù)研究發(fā)展計劃(863計劃)(2012AA01A504)

劉治軍，男，碩士研究生，研究方向為航天器通信技術(shù)。Email：2362172973@qq.com。

TN911.7

A

10.3969/j.issn.1673-8748.2015.06.004