張立民，譚繼遠(yuǎn)，閆文君

(1.海軍航空大學(xué)信息融合研究所,山東煙臺 264001;2.海軍航空大學(xué)研究生一隊(duì),山東煙臺 264001)

0 引 言

近年來無線通信技術(shù)的迅猛發(fā)展對數(shù)據(jù)傳輸速率、傳輸效率和頻帶利用率提出了更高的要求。選擇高效可行的調(diào)制解調(diào)手段，對提高信號的有效性和可靠性起著至關(guān)重要的作用[1]。隨著通信技術(shù)的發(fā)展，通信信號的調(diào)制方式也多種多樣，當(dāng)對接收信號進(jìn)行解調(diào)或消除未知干擾信號時(shí)，我們需要先對信號進(jìn)行調(diào)制識別。因此調(diào)制識別在許多領(lǐng)域有著廣泛的應(yīng)用,例如電子戰(zhàn),電子干擾和頻譜檢測等領(lǐng)域[2-3]。

目前調(diào)制識別的方法主要分為三種：基于似然的識別方法、基于特征值的識別方法和基于深度學(xué)習(xí)的識別方法[4]。基于似然的識別方法應(yīng)用最廣的是基于最大似然的調(diào)制分類器[5-6]，該方法識別率最高，但是計(jì)算量大，計(jì)算復(fù)雜度高。特征值的識別方法主要是基于高階累積量的識別[6-7]，相比于似然識別，計(jì)算較為簡單，但是識別率比基于似然的低。基于深度學(xué)習(xí)的識別方法還未發(fā)展完善，訓(xùn)練模型和訓(xùn)練樣本的生成與選擇會嚴(yán)重影響識別效果。目前調(diào)制識別研究的難點(diǎn)和熱點(diǎn)主要集中在以下三個(gè)方面[8]：一是如何在非理想條件下進(jìn)行調(diào)制識別[9-11]；二是如何改進(jìn)基于特征值的調(diào)制識別[12-15]；三是如何構(gòu)造一個(gè)簡單有效的神經(jīng)網(wǎng)絡(luò)模型與訓(xùn)練樣本進(jìn)行調(diào)制識別。Aslam等利用遺傳編程和K階最近鄰(KNN)分類器，選擇累積量作為特征進(jìn)行分類[16]，具有良好的分類性能，但是算法復(fù)雜度較高。Dobre等提出了基于秩序統(tǒng)計(jì)的低復(fù)雜度自動調(diào)制分類[17]，但是在分類性能和計(jì)算復(fù)雜度之間需要有一個(gè)權(quán)衡。近年來，基于神經(jīng)網(wǎng)絡(luò)的方法也廣泛應(yīng)用于調(diào)制識別[18][19]，Timothy等提出了利用卷積神經(jīng)網(wǎng)絡(luò)(CNN)進(jìn)行調(diào)制識別[18]，在低信噪比和樣本數(shù)較少的情況下，該算法變現(xiàn)出良好的性能，但是需要提前進(jìn)行神經(jīng)網(wǎng)絡(luò)訓(xùn)練，大量有效的訓(xùn)練樣本的生成還需有所改進(jìn)。王方剛等首次提出了基于經(jīng)驗(yàn)累積分布函數(shù)的KS檢測(Kolmogorov-Smirnov test)方法[13]，并對其進(jìn)行了改進(jìn)，提出了基于折疊經(jīng)驗(yàn)累積分布函數(shù)(cdf)的調(diào)制識別[1]，該方法主要討論了在加性高斯白噪聲條件下對接收信號進(jìn)行識別，識別率得到了一定的提升，但是該論文未對加信道條件下的調(diào)制識別進(jìn)行詳細(xì)的討論。

考慮信道參數(shù)對調(diào)制識別的影響，本文提出了基于折疊經(jīng)驗(yàn)累積分布函數(shù)(cdf)的最小均方差的識別方法，對瑞利時(shí)不變信道下的信號進(jìn)行調(diào)制識別。首先建立了經(jīng)過信道和噪聲的QAM調(diào)制信號模型，然后對接收信號的cdf0和標(biāo)準(zhǔn)調(diào)制信號的cdfk(k=1,2,3分別對應(yīng)4-QAM、 16-QAM和 64-QAM)進(jìn)行折疊，利用最小均方差準(zhǔn)則判斷檢測信號折疊后的cdf0和cdfk的擬合程度，均方差小的，擬合程度高，從而達(dá)到調(diào)制識別的目的。

1 信號模型和方法描述

1.1 信號模型

假設(shè)接收信號的模型為

yn=Hxn+ωn,n=1,2,3,…,N

(1)

其中yn為接收信號，H為衰落信道，xn為復(fù)數(shù)發(fā)射調(diào)制符號，ωn為時(shí)刻n處的噪聲采樣，為加性循環(huán)復(fù)高斯白噪聲。 傳輸?shù)姆杧n是4-QAM、 16-QAM或 64-QAM理想調(diào)制中的星座圖點(diǎn)組成的序列，該星座屬于一組調(diào)制格式{M1,M2,…Mk}。假設(shè)信道H為時(shí)不變?nèi)鹄诺?，循環(huán)復(fù)高斯白噪聲ωn～N(0,σ2)，因此ωn的實(shí)部和虛部獨(dú)立同分布，且分布服從N(0,σ2/2)考慮單位功率星座，因此信噪比SNR=1/σ2。調(diào)制識別就是研究接收信號{yn},n=1,2,3,…,N為何種調(diào)制方式問題。

1.2 基于cdf的最小均方差準(zhǔn)則

當(dāng)對接收信號{yn}進(jìn)行調(diào)制識別時(shí)，我們可以找到關(guān)于{yn}的一個(gè)特征序列{zn}，當(dāng)進(jìn)行PSK調(diào)制識別時(shí)， {zn}可以是信號相位；當(dāng)進(jìn)行QAM調(diào)制識別時(shí)，{zn}可以是信號幅度或者信號的實(shí)部和虛部[10]。然后計(jì)算關(guān)于{zn}的經(jīng)驗(yàn)累積分布函數(shù)F0和理論的理想4-QAM、 16-QAM、 64-QAM調(diào)制{zn}的經(jīng)驗(yàn)累積分布函數(shù)Fk,k=1,2,3，其中k=1,2,3分別對應(yīng)4-QAM、 16-QAM、 64-QAM三種調(diào)制方式。接下來對上述cdf進(jìn)行折疊，使不同調(diào)制的cdf之間有明顯的區(qū)別。最后，利用最小均方差準(zhǔn)則分別計(jì)算折疊后的F0和Fk,k=1,2,3之間的均方差，計(jì)算公式如下，

(2)

其中，Ek代表均方差，F(xiàn)0(zn)和Fk(zn) 分別代表特征序列{zn} 折疊后的cdf和折疊后的理想cdf。

根據(jù)最小均方差準(zhǔn)則

r=argmin(Ek),k=1,2,3

(3)

均方差越小的，表示兩者的cdf擬合程度越高，取三者中均方差最小的即為擬合程度最高的，其所對應(yīng)的調(diào)制方式即為所接收信號的調(diào)制方式。

2 QAM調(diào)制識別

對于QAM調(diào)制識別，4-QAM、 16-QAM和 64-QAM調(diào)制的單位能量星座的信號點(diǎn)集合分別為

-1,1,3,5,7}

(4)

對于QAM調(diào)制信號，信號的虛部和實(shí)部相互獨(dú)立且同分布，因此可以將接收信號yn的虛部和實(shí)部組成特征序列{zn}，z2n-1=R{yn}，z2n=I{yn}(R{·}、I{·}分別代表復(fù)數(shù)的實(shí)部與虛部)，對于上述系統(tǒng)模型，因?yàn)樵肼暒豱為循環(huán)復(fù)高斯白噪聲，其實(shí)部與虛部獨(dú)立同分布于N(0,σ2/2)，因此{(lán)z2n}的實(shí)部與虛部也獨(dú)立同分布。

E(z2n-1)=E(z2n)=E(R(yn))=E(I(yn))=

E(R(Hxn+ωn))=E(I(Hxn+ωn))=

E(R(Hxn))+E(R(ωn))=

R(Hxn)

(5)

其中，E(z2n) 表示求{z2n} 序列的均值，H代表信道參數(shù)。

D(z2n-1)=D(z2n)=D(R(yn))=D(I(yn))=

D(R(Hxn+ωn))=D(I(Hxn+ωn))=

(6)

其中，D(z2n)表示求{z2n}序列的方差。

因此，z2n～N(0,σ2/2)。z2n的pdf為

(7)

其中x∈R{Mk}指的是信號的實(shí)部，H代表信道參數(shù)。將pdf經(jīng)過積分得到關(guān)于z的cdf為

(8)

此時(shí)，未經(jīng)過折疊的pdf和cdf如圖1所示

圖1 未折疊的pdf與cdf

從圖中可以看出，三者的pdf區(qū)分度不是特別高，特征不明顯。接下來，通過折疊，使得三者的pdf更容易區(qū)分。

因?yàn)樯鲜鰌df關(guān)于z=0對稱，接下來將pdf沿著z=0對折，得到一次折疊的pdf[1]如下

z[1]≥0

(9)

其中z[i]代表第i次折疊后的特征序列。

當(dāng)z[1]≤0時(shí)，μk,[1](z[1])=0。將一次折疊的pdf[1]經(jīng)過積分得到關(guān)于z[1]的cdf[1]為

(10)

其中，x∈R+{Mk}指的是星座點(diǎn)的正實(shí)部，經(jīng)過一次折疊后的pdf[1]和cdf[1]如圖2所示。

圖2 經(jīng)過一次折疊的pdf[1]與cdf[1]

其中輸入決策z[i]與s[i]關(guān)系如下

z[i+1]=|z[i]-s[i+1]|+s[i+1]

(11)

經(jīng)過二次折疊的pdf[2]如下

μk,[2](z[2])=μk,[1](z[2])+μk,[1](2s[2]-z[2]),

z[2]≥s[2]

(12)

z[2]≤s[2]時(shí)，μk,[2](z[2])=0。將二次折疊的pdf[2]積分得到cdf[2]如下

(13)

經(jīng)過二次折疊的pdf[2]和cdf[2]如圖3所示。

圖3 經(jīng)過兩次折疊的pdf[2]與cdf[2]

從圖中可以看出經(jīng)過兩次折疊的pdf[2]分為兩個(gè)區(qū)域,其中一個(gè)區(qū)域4-QAM的pdf[2]較大，另一個(gè)區(qū)域16-QAM和64-QAM的pdf[2]較大。積分后兩次折疊的4-QAM的cdf[2]和兩次折疊的16-QAM、64-QAM的cdf[2]的區(qū)分度較大，這時(shí)用到最小均方差準(zhǔn)則計(jì)算出給定樣本的二次折疊cdf[2]與兩次折疊4-QAM的cdf[2]、兩次折疊的16-QAM的cdf[2]的均方差E1,[2]和E2,[2]，若E1,[2]E2,[2]，則信號采用的是16-QAM或者是64-QAM。這時(shí)需要對pdf[2]進(jìn)行第三次折疊，使得折疊后的16-QAM和64-QAM的pdf[2]區(qū)分度明顯，三次折疊的pdf[3]與cdf[3]表達(dá)式如下

μk,[3](z[3])=μk,[2](z[3])+μk,[2](2s[3]-z[3]),

z[3]≥s[3]=μk,[1](2s[2]-2s[3]+z[3])+

μk,[1](2s[3]-z[3])+μk,[1](2s[2]-z[3])+

μk,[1](z[3])

(14)

(15)

經(jīng)過三次折疊的pdf[3]和cdf[3]如圖4所示。

圖4 經(jīng)過三次折疊的pdf[3]和cdf[3]

此時(shí)再一次用到最小均方差準(zhǔn)則判斷擬合程度,若E3,[3]>E2,[3]則為16-QAM調(diào)制，否則，為64-QAM調(diào)制。

檢測算法可以歸納為圖5所示的決策樹。

圖5 最小均方差方法決策樹

3 仿真和結(jié)果

3.1 仿真條件

仿真經(jīng)過1000次蒙特卡洛仿真，接收信號時(shí)為隨機(jī)產(chǎn)生的QAM調(diào)制信號，無特殊說明，仿真條件設(shè)置如下：樣本數(shù)N=4000，調(diào)制信號xn為隨機(jī)產(chǎn)生的序列，H信道模型為時(shí)不變?nèi)鹄诺?，信噪比定義為SNR=10log10(1/σ2)

在實(shí)驗(yàn)中，采用正確的識別概率P衡量算法的性能

p=(λ=ξ|ξ),ξ∈Ω

(16)

Ω={4-QAM,16-QAM,64-QAM}

3.2 AWGN信道下，不同SNR時(shí)該算法對4-QAM和16-QAM調(diào)制的正確識別概率

圖6 AWGN條件下，不同分類器對4-QAM和16-QAM調(diào)制的正確識別概率

由上圖可知，在只加高斯白噪聲的情況下，當(dāng)信噪比小于2 dB時(shí)，該算法性能較差。但當(dāng)信噪比大于2 dB時(shí),算法表現(xiàn)出良好的性能，并且在信噪比達(dá)到3 dB時(shí),正確識別率由于基于折疊KS的識別方法，并且能快速達(dá)到最佳識別效果。

3.3 瑞利信道下，不同SNR時(shí)不同分類器對4-QAM、16-QAM和64-QAM調(diào)制的正確識別概率

對信號在默認(rèn)條件下進(jìn)行仿真,不同識別方法對QAM調(diào)制的正確識別概率如圖7所示，由仿真結(jié)果可以看出相同SNR情況下，文中所提出的基于折疊cdf的最小均方差識別方法識別率絕大多數(shù)最高，且低信噪比情況下，識別概率較其他方法有明顯的提高。隨著SNR增大，識別率增大，當(dāng)SNR=20 dB時(shí)，最大識別率便可以達(dá)到95%，而其他方法在SNR=25 dB時(shí)，識別率才達(dá)到最大。并且在高SNR下比基于累積量的分類器的性能明顯好，較基于累積量的識別算法，最大識別率提高了5%。

圖7 不同SNR時(shí)不同分類器對QAM調(diào)制的正確識別概率

3.4 瑞利信道下，樣本數(shù)對4-QAM和16-QAM調(diào)制的正確識別概率的影響

在不同樣本數(shù)的條件下對4-QAM和16-QAM調(diào)制的正確識別概率如圖8所示，可以看出，算法的正確識別率隨著樣本數(shù)N的增大而增大。并且隨著N的增加，算法的收斂性和正確識別率都得到了提高。

圖8 不同樣本數(shù)時(shí)對4-QAM和16-QAM的正確識別概率

3.5 瑞利信道下，樣本數(shù)對4-QAM、16-QAM和64-QAM調(diào)制的正確識別概率的影響

在不同樣本數(shù)的條件下對4-QAM、16-QAM和64-QAM調(diào)制的正確識別概率如圖9所示，可以看出，算法的正確識別率隨著樣本數(shù)N的增大而增大。這是由于樣本數(shù)N增多，pdf和cdf的統(tǒng)計(jì)特性會更加明顯，更加有利于擬合程度判斷。在樣本數(shù)N=4000時(shí)，算法在低信噪比的識別概率得到了明顯提升。

圖9 不同樣本數(shù)時(shí)的對4-QAM、16-QAM和64-QAM的正確識別概率

3.6 復(fù)雜度分析

3.6.1算法運(yùn)行時(shí)間分析

表1 不同識別方法的識別時(shí)間表

比較了基于累計(jì)量識別、基于KS識別[13]、基于最大似然識別[4]本文算法的復(fù)雜度。從上表中可以看出，基于累積量的識別方法耗時(shí)最短，但是根據(jù)圖6仿真可以看出識別效果較差?；谧畲笏迫坏淖R別耗時(shí)最長，但是該算法的識別效果最佳，本文算法折中了兩者，在保證一定的識別效果的同時(shí)，識別時(shí)間得到了改善。雖然識別時(shí)間并不直接反映算法的復(fù)雜度，但是從一定程度上也反映了算法的復(fù)雜度[17]。

3.6.2算法復(fù)雜度分析

4 結(jié) 語

本文研究了4-QAM、16-QAM和64-QAM的調(diào)制識別問題，考慮了信道以及噪聲的影響。本文方法基于最小均方差擬合，對QAM調(diào)制進(jìn)行識別。建立了QAM調(diào)制的信號模型，推導(dǎo)了接收信號的數(shù)學(xué)分布；計(jì)算了關(guān)于接收信號的pdf和cdf，并對其進(jìn)行折疊，通過折疊能看出不同QAM調(diào)制的cdf有著區(qū)別，通過最小二乘法判定cdf擬合程度達(dá)到QAM調(diào)制識別的目的。仿真結(jié)果表明，方法性能較好，尤其是在低信噪比情況下，正確識別率得到了一定的提高。本文方法的主要缺點(diǎn)時(shí)當(dāng)樣本數(shù)較少時(shí)，識別性能仍舊有待提高。