李曉光，潘克剛，齊 帥，潘 晨

（中國(guó)人民解放軍陸軍工程大學(xué) 通信工程學(xué)院，江蘇 南京 210007）

0 引 言

自動(dòng)調(diào)制分類是一種識(shí)別接收信號(hào)的調(diào)制類型技術(shù)。它在諸如智能調(diào)制解調(diào)器、軟件定義無線電、頻譜監(jiān)測(cè)、電子監(jiān)控和威脅分析、信號(hào)確認(rèn)、干擾識(shí)別以及可編程和可重構(gòu)無線電系統(tǒng)[1-4]中有著廣泛應(yīng)用。其中，關(guān)鍵技術(shù)是設(shè)計(jì)能夠識(shí)別使用的調(diào)制類型的智能分類器。

根據(jù)調(diào)制識(shí)別分類研究情況，在分類過程中分別使用多個(gè)特征提取，一般涉及二階、四階和六階的矩和累積量[5-7]。為了有效識(shí)別調(diào)制信號(hào)，基于計(jì)算的八階累計(jì)量，開發(fā)和驗(yàn)證用于MPSK相位多級(jí)數(shù)字調(diào)制的自動(dòng)分類技術(shù)。相比于前面的研究，在特征提取過程中選取與3個(gè)八階累積量關(guān)系相關(guān)的特征參數(shù)。

所提出的分類算法將使用Matlab進(jìn)行仿真驗(yàn)證，并使用GNU Radio和Boost庫(kù)在C++中實(shí)現(xiàn)該識(shí)別模塊。GNU Radio和USRP軟件無線電是一種開放的可重新編程的無線電通信系統(tǒng)，根據(jù)軟件改變便攜式USRP硬件設(shè)置，結(jié)合GNU Radio軟件就可設(shè)計(jì)所需的通信系統(tǒng)。

1 多相數(shù)字調(diào)制（MPSK）

M進(jìn)制相移鍵控信號(hào)載波相位有M個(gè)值，已調(diào)信號(hào)的時(shí)域表達(dá)式為：

MPSK調(diào)制的星座平面中，點(diǎn)分為兩部分：相位分量和正交分量，分別稱為描述星座Xi的表達(dá)式的實(shí)部和虛部：

其中A是信號(hào)幅度值，φ是星座的相位角。

2 高階矩和高階累積量

X=cul(x1,x2,…xk)是一向量，Ix={1,2,…,k}是其指示符集。mx( I )和cx( I )是X的子向量X1的矩和累積量，可以寫出：

于是，得到矩-累積量公式（簡(jiǎn)稱M-C公式）：

以一種類似的方式，累積量被定義為一個(gè)復(fù)雜的、零均值和階數(shù)為p的隨機(jī)過程Cpq,x=cum[Xp-qX—p]。

對(duì)于MPSK調(diào)制復(fù)數(shù)隨機(jī)過程，定義了算子cum[·]。對(duì)于k=8，l=0、1、4，根據(jù)M-C公式有：

3 分類器設(shè)計(jì)

根據(jù)以前的結(jié)果，在MPSK星座計(jì)算八階累積量變量，信號(hào)模型為加性高斯白噪聲數(shù)字調(diào)制信號(hào)。假設(shè)在載波、定時(shí)同步的前提下，經(jīng)過下變頻，匹配濾波器輸出的復(fù)基帶序列表達(dá)式為：

其中，p為信號(hào)的平均功率，an為發(fā)送碼元序列，Ts是碼元周期，θc為相位偏差，ωc為載頻，g(·)為發(fā)送的碼元波形，n(i)為零均值的加性復(fù)高斯白噪聲序列，N為觀測(cè)數(shù)據(jù)的長(zhǎng)度。

接收端已知載頻信息且達(dá)到定時(shí)同步的情況下，對(duì)待識(shí)別信號(hào)進(jìn)行下變頻處理，得到了采樣復(fù)信號(hào)序列的表達(dá)式為：

通過對(duì)高階累積量理論的研究可以知道，對(duì)于信號(hào)，階數(shù)k＞2的累積量對(duì)噪聲nk具有抑制作用。因此，這里只分析ak的各累積量。對(duì)于MPSK信號(hào)而言，它的碼元同步采樣復(fù)信號(hào)的表達(dá)式為：

在假定發(fā)送的信息序列為獨(dú)立同分布的情況下，計(jì)算2PSK信號(hào)的高階矩，得到如下八階累積量：

在MPSK（2PSK、4PSK、8PSK和 16PSK） 計(jì)算八階累積量結(jié)果，如表1所示。

從表1可以看出，累積量結(jié)果取決于信號(hào)的平均功率和相位偏差。

表1 MPSK八階累積量

根據(jù)表1中提到的理論結(jié)果，提出將其定義為用于識(shí)別調(diào)制的特征向量，具有以下關(guān)系：

根據(jù)式（15）得到FM的確定性值。考慮表1的結(jié)果，對(duì)于每個(gè)要分類的調(diào)制類型，有：

可以采用定義新的特征參數(shù)：

來實(shí)現(xiàn)MPSK信號(hào)的分類[8]。

可見，DΩ結(jié)果與信號(hào)的平均功率和相位偏差無關(guān)。所以，分類器邏輯設(shè)計(jì)如下：

假設(shè)觀測(cè)值為FR=[+]，實(shí)測(cè)值D=+對(duì)不同的MPSK調(diào)制類型實(shí)現(xiàn)分類。

4 仿真分析

仿真分析MPSK信號(hào)分類識(shí)別性能。通過Matlab產(chǎn)生仿真MPSK信號(hào)，仿真條件為：采樣頻率60 kHz，載波頻率10 kHz，碼元速率為2 500 Hz，符號(hào)數(shù)為1 000，噪聲為加性高斯窄帶噪聲，SNR=-5～25 dB。

圖1為根據(jù)這一特征參數(shù)進(jìn)行500次蒙特卡羅試驗(yàn)對(duì)信號(hào)進(jìn)行分類識(shí)別得到的結(jié)果。該參數(shù)在SNR＞2.5 dB的情況下，可以對(duì)MPSK信號(hào)分類。

圖2為輸入符號(hào)數(shù)N分別為1 000和500時(shí)，進(jìn)行1 500次蒙特卡洛試驗(yàn)的出的的識(shí)別率。圖2（a）中，在SNR≥6 dB時(shí)，識(shí)別率可達(dá)到100%；而N為500的圖2（b）中顯示，在SNR=7時(shí)，識(shí)別率還沒達(dá)到100%，尤其是16PSK表現(xiàn)更差?？梢姡?hào)數(shù)越多，識(shí)別率越高。

圖1 分類識(shí)別結(jié)果

圖2 識(shí)別率

綜合識(shí)別結(jié)果和識(shí)別率情況，當(dāng)信號(hào)SNR＞5 dB時(shí)，可以實(shí)現(xiàn)正確的分類，驗(yàn)證了分類方法的有效性和準(zhǔn)確性。

5 軟件無線電平臺(tái)設(shè)計(jì)

5.1 軟件無線電平臺(tái)簡(jiǎn)介

GNU Radio和USRP軟件無線電平臺(tái)，軟件部分主要基于Linux操作系統(tǒng)，通信系統(tǒng)模型由C++和Python語(yǔ)言編程構(gòu)成。C++語(yǔ)言用來編寫各種信號(hào)處理模塊，而這些信號(hào)處理模塊在GNU Radio中被稱為“Block”；Python用來連接各個(gè)Block，使之成為一個(gè)腳本文件，從而實(shí)現(xiàn)通信系統(tǒng)的某些功能，而這個(gè)腳本文件被稱為“Flow graph”。USRP是硬件部分，作用是收發(fā)射頻信號(hào)并將其轉(zhuǎn)換為基帶信號(hào)連接PC機(jī)[9]。

通過軟件開發(fā)通信系統(tǒng)具有良好的性能，只需利用系統(tǒng)自帶的通信模塊功能和開發(fā)的Block，就可改變系統(tǒng)的功能。因此，利用該系統(tǒng)進(jìn)行仿真實(shí)現(xiàn)和做成工程項(xiàng)目也已經(jīng)普遍。此外，該平臺(tái)軟件是一個(gè)開源項(xiàng)目，硬件部分的價(jià)格不高，技術(shù)要求較低，用戶均可參與該平臺(tái)的建設(shè)[10-11]。

5.2 分類器設(shè)計(jì)

5.2.1 C++程序設(shè)計(jì)流程圖

C++程序設(shè)計(jì)流程圖，如圖3所示[12-13]。

將輸入分成M段長(zhǎng)度為L(zhǎng)的數(shù)據(jù)，根據(jù)式（15）、式（16）和式（17）計(jì)算D(i),i=1,2,…,M，求出平均值D。根據(jù)比較設(shè)置的閥值和參照分類器邏輯設(shè)計(jì)（式（18））輸出調(diào)制類型識(shí)別結(jié)果，并將數(shù)據(jù)推送到下一個(gè)信號(hào)處理模塊。

5.2.2 grc/myblk_CUMA.xml設(shè)計(jì)

主要設(shè)計(jì)Block名稱、控制參數(shù)和輸入輸出口，如表2所示。

圖3 流程圖

表2 相關(guān)設(shè)計(jì)參數(shù)

其中，在程序運(yùn)行過程中，實(shí)時(shí)改變參數(shù)的值callback關(guān)鍵字，調(diào)用C++的成員函數(shù)，修改BLOCKSIZE的值，用來控制每次截取數(shù)據(jù)塊大小。

5.2.3 Block設(shè)計(jì)

分類器為同步類型模塊名稱為CUMA，輸入輸出端數(shù)據(jù)類型為復(fù)數(shù)（Complex）。設(shè)計(jì)2個(gè)輸出口連接不同的解調(diào)模塊，配置參數(shù)為處理數(shù)據(jù)模塊大小，默認(rèn)24 000。

圖4為3部分的拼合。右上方為Block已經(jīng)注冊(cè)到系統(tǒng)內(nèi)，左上方為創(chuàng)建的Block調(diào)用情況，下方則為Block設(shè)置參數(shù)。

圖4 分類器

5.3 軟件無線電平臺(tái)設(shè)計(jì)和結(jié)果

發(fā)射機(jī)（如圖5所示）隨機(jī)信號(hào)生成模塊，調(diào)制模塊生成MPSK復(fù)包絡(luò)信號(hào)，使用USRP設(shè)置中心頻率為1 GHz的信號(hào)進(jìn)行無線發(fā)射。為了監(jiān)控星座圖，可在發(fā)射端構(gòu)建信號(hào)監(jiān)測(cè)模塊。接收機(jī)如圖6所示，UHD中心頻率1 GHz，包括分類器識(shí)別模塊、可控信道信道模塊、CMA盲均衡模塊、鎖相環(huán)模塊和信號(hào)監(jiān)測(cè)模塊。

圖5 發(fā)射圖

圖6 接收?qǐng)D

圖7 接收端信號(hào)星座圖

圖7 是模塊搭建后測(cè)試2PSK、4PSK、8PSK、16PSK信號(hào)經(jīng)過可變信道進(jìn)行參數(shù)設(shè)置，經(jīng)過分類器模塊處理后的星座圖?？梢?，分類器能準(zhǔn)確識(shí)別輸入信號(hào)，模塊性能良好。

6 結(jié) 語(yǔ)

在計(jì)算八階累積量的基礎(chǔ)上，定義新的特征向量和特征參數(shù)，而這些特征不受信號(hào)的幅度和相位偏差影響，可利用特征值設(shè)計(jì)MPSK信號(hào)調(diào)制識(shí)別分類器。通過Matlab仿真進(jìn)行對(duì)比論證，結(jié)果證明了該分類器的有效性和正確性。此外，基于GNU Radio和URSP平臺(tái)實(shí)現(xiàn)了模塊化設(shè)計(jì)，通過MPSK的調(diào)制信號(hào)驗(yàn)證，證明分類器模塊具有良好的性能，為該方向工程實(shí)現(xiàn)奠定了良好的試驗(yàn)基礎(chǔ)。