張廣吉，袁仕繼，劉志華，黃文晶

(解放軍63888 部隊，河南 濟(jì)源 454650)

1 問題描述

人們對信號進(jìn)行處理，其最終目的就是要對單通道或多通道的有限長觀測數(shù)據(jù)進(jìn)行必要的處理，從而檢測出隱含于這些觀測數(shù)據(jù)中的各種重要信息并加以利用。高階統(tǒng)計量在信號處理與系統(tǒng)分析中扮演著一個極為重要的角色。根據(jù)最近的資料顯示，在通信、生物醫(yī)學(xué)工程、語音處理、地震信號分析、圖象處理、雷達(dá)、聲納等領(lǐng)域都進(jìn)行了有關(guān)高階統(tǒng)計量處理的研究，具體應(yīng)用于時延估計、系統(tǒng)辨識、自適應(yīng)濾波及陣處理等方面。高階統(tǒng)計量能提供比功率譜更多的有用信息，能夠有效地檢測信號幅度以外的其他信息，這就具有明顯的優(yōu)點。高階統(tǒng)計量作為非高斯信號處理的主要分析工具，不僅提供高階相關(guān)信息，而且能夠衡量隨機序列偏離正態(tài)的程度，并對加性高斯噪聲(白色或有色)不敏感。這些特點都非常有利于特征參數(shù)的提取。

2 基于高階統(tǒng)計量的信號特征分析

所謂高階統(tǒng)計量，通常理解為高階矩、高階累積量以及它們的譜。本文主要討論高階矩和高階累積量，對它們的譜暫不作分析。首先將利用特征函數(shù)來定義高階矩和高階累積量，然后引出它們的有關(guān)性質(zhì)，最后推導(dǎo)出高階矩和高階累積量之間的變換關(guān)系。

2.1 通信信號高階統(tǒng)計量的算法實現(xiàn)

通過接收機，在不同的時間里對幾個不同電臺所發(fā)射的不同調(diào)制樣式的短波通信信號進(jìn)行了不同采樣速率的接收。對接收到的信號進(jìn)行去均值處理，去其直流成分，然后再對信號的高階統(tǒng)計量(二階、三階、四階)進(jìn)行計算，并畫出曲線圖進(jìn)行比較。

如果，x={xi}，i =1，…，N 表示隨機序列或信號的采樣，則序列的矩為:，其中E{}表示期望值的運算符。其中p(x)表示xn的概率密度函數(shù)。當(dāng)采集得到的信號時間長度比較長時，認(rèn)為是各態(tài)歷經(jīng)和平穩(wěn)的，同時進(jìn)行了去均值處理。

這樣則采樣信號的矩可由下式估計:

同樣它的n 階累計量:

這樣，前4 個累計量可以用下式計算:

2.2 信號的高階統(tǒng)計量實驗計算

對同一電臺、同一種調(diào)制樣式信號進(jìn)行分析，由于比較得到AM 信號的高階統(tǒng)計量參數(shù)較穩(wěn)定，分別將兩個時間段內(nèi)對同一電臺AM 信號的采樣進(jìn)行了分析畫圖比較發(fā)現(xiàn)，測量的時間不同對同一電臺AM 信號的高階統(tǒng)計量參數(shù)影響不大，它們的各階統(tǒng)計量的波動范圍都幾乎在同一個很小的區(qū)間內(nèi)，這反應(yīng)了該電臺的同一調(diào)制樣式的信號的相對穩(wěn)定性。現(xiàn)以在不同時間內(nèi)對AM 666 kHz 電臺的二階矩和二階累計量為例，如圖1 所示。

對4 個不同電臺的AM 信號進(jìn)行了采樣，它們是800 kHz、864 kHz、936 kHz、1 592.5 kHz。發(fā)現(xiàn)它們各自的高階統(tǒng)計量參數(shù)仍是較平穩(wěn)，不同的是各自的中心值和波動范圍不同，以此可以區(qū)分各自電臺。現(xiàn)將它們的各階統(tǒng)計量畫出，如圖2 所示。

對所有電臺的高階統(tǒng)計量參數(shù)計算、畫曲線圖后發(fā)現(xiàn)一個共同特點:對一定采樣速率得到的采樣信號，當(dāng)對它們合并數(shù)目越多時(即新的采樣信號的時間長度越長)，信號的高階統(tǒng)計量參數(shù)曲線越平穩(wěn)，即參數(shù)越集中;而對它們截的越短時(即新的采樣信號的時間長度越短)，信號的高階統(tǒng)計量參數(shù)曲線越雜亂，抖動越大，即參數(shù)越分散?，F(xiàn)以電臺800 kHz 的三階統(tǒng)計量為例，如圖3 所示。

圖1 AM 666 kHz 電臺的二階矩和二階累計量

圖2 不同信號的高階統(tǒng)計量

以電臺864 kHz 的三階統(tǒng)計量為例，如圖4 所示。通過實驗結(jié)果的比較，得出了上述結(jié)論。在對同一電臺的同一種調(diào)制樣式(AM)信號在不同時間內(nèi)的比較中得到，高階統(tǒng)計量的值相對穩(wěn)定，區(qū)別也只是隨機噪聲帶來的。這就為利用高階統(tǒng)計量去識別這一信號的電臺提供了實驗依據(jù)。

同樣，在對不同采樣時間長度的比較中得到:采樣時間越長，高階統(tǒng)計量越穩(wěn)定。對于這種現(xiàn)象，可以這樣理解:因為在計算信號高階統(tǒng)計量的時候，所采用的算法的前提條件是信號為各態(tài)歷經(jīng)和平穩(wěn)的，而實際試驗中采樣的信號，由于噪聲的影響和采樣時間較短等因素，分析的信號是不太滿足上述條件。但是如果采樣信號時間越長，它就是越接近各態(tài)歷經(jīng)和平穩(wěn)的，也就越滿足前提條件。當(dāng)合并信號，實際上是在加長采樣信號的時間長度，而截短信號，是在減少采樣信號時間。而同一電臺發(fā)射的同一調(diào)制樣式的信號基本是穩(wěn)定的，這樣當(dāng)越接近條件時，對實驗的效果越好，高階統(tǒng)計量參數(shù)曲線也就越顯得平穩(wěn)。

圖3 800 kHz 的三階統(tǒng)計量

圖4 864 kHz 的三階統(tǒng)計量

3 結(jié)束語

由于諸多因素的影響，對通信信號高階統(tǒng)計量特征的研究還沒有深入下去。這篇論文僅是對該方面的一個淺顯的探索，只是其中的一小部分。對于下一步工作，可以從以下幾個方面開展:

1)解決算法。由于計算信號高階統(tǒng)計量的計算量太大，總是尋找簡便的算法，這樣就會出現(xiàn)不少的前提條件，如何使采樣得到的信號更滿足前提條件，這是需要考慮的;同時，對由于前提條件的不同而產(chǎn)生的高階統(tǒng)計量的誤差，還要進(jìn)行分析校正。

2)信號樣本。由于本文只著重對AM 信號進(jìn)行了分析，而沒對其他的信號作進(jìn)一步研究，這樣所得到的參數(shù)還不充足?？梢栽趯Ω鞣N調(diào)制樣式的信號采樣時，不斷改變采樣信號的時間長度，看看在什么樣的一個時間范圍內(nèi)，信號的高階統(tǒng)計量可以進(jìn)一步的穩(wěn)定，這對信號識別很有借鑒價值。

3)高階矩譜和高階累積量譜。在這次分析中，只是針對于信號的時域分析，在對信號歸一化(去均值)后利用算法進(jìn)行的計算。而作為工程上所研究的高階統(tǒng)計量還包括高階矩譜和高階累積量譜，即站在頻譜的角度上對高階統(tǒng)計量進(jìn)行分析。高階矩譜和高階累積量譜會含有對提取信號特征更為有用的信息和規(guī)律。

［1］樊昌信. 通信原理［M］. 5 版. 北京: 國防工業(yè)出版社，1999.

［2］Richard A.通信電子戰(zhàn)系統(tǒng)導(dǎo)論［M］.北京:電子工業(yè)出版社，2003.

［3］王炳武.MATLAB5.3 教程［M］.北京:中國水利水電出版社，2000.

［4］張賢達(dá). 現(xiàn)代信號處理［M］. 北京: 清華大學(xué)出版社，1995.

［5］姚文冰.基于高階累積量的抗噪語音識別［D］.武漢:華中科技大學(xué)，2001.

［6］梁清泉.基于高階統(tǒng)計量的波達(dá)方向估計技術(shù)的研究［D］.西安:西北工業(yè)大學(xué)，2002.

［7］宋文武，丁志堯，張凱，等. 短波通信系統(tǒng)阻塞干擾控制方法研究［J］.裝備環(huán)境工程，2010（2）:68-70