范真真,朱立忠,張笑宇,李 瑩

(沈陽(yáng)理工大學(xué) 信息科學(xué)與工程學(xué)院,遼寧 沈陽(yáng) 110159)

一種小波變換的數(shù)字信號(hào)調(diào)制方式識(shí)別方法研究

范真真,朱立忠,張笑宇,李 瑩

(沈陽(yáng)理工大學(xué) 信息科學(xué)與工程學(xué)院,遼寧 沈陽(yáng) 110159)

數(shù)字通信信號(hào)調(diào)制方式識(shí)別有重要的現(xiàn)實(shí)意義。針對(duì)各數(shù)字信號(hào)的特點(diǎn)，利用Haar小波函數(shù)歸一化前后小波變換模值，實(shí)現(xiàn)對(duì)數(shù)字信號(hào)的類間識(shí)別；根據(jù)數(shù)字信號(hào)進(jìn)制數(shù)，進(jìn)行數(shù)字信號(hào)的類內(nèi)識(shí)別，并通過(guò)軟件仿真驗(yàn)證該方法的有效性。仿真結(jié)果表明，基于小波變換的調(diào)制方式識(shí)別方法可對(duì)數(shù)字信號(hào)進(jìn)行有效的識(shí)別。

數(shù)字信號(hào);小波變換;調(diào)制方式識(shí)別;幅度歸一化

數(shù)字通信信號(hào)調(diào)制方式識(shí)別技術(shù)作為信號(hào)處理領(lǐng)域中的關(guān)鍵技術(shù)被廣泛應(yīng)用于軍事和民用領(lǐng)域[1]。數(shù)字信號(hào)調(diào)制方式識(shí)別研究的主要問(wèn)題包括如FSK、PSK、ASK的類間識(shí)別及如BPSK、QPSK的類內(nèi)識(shí)別[2]。準(zhǔn)確地進(jìn)行數(shù)字信號(hào)調(diào)制方式識(shí)別是后續(xù)數(shù)字信號(hào)處理的前提和基礎(chǔ)。隨著數(shù)字信號(hào)調(diào)制日益廣泛的應(yīng)用，其調(diào)制方式識(shí)別問(wèn)題也越來(lái)越受到重視，針對(duì)各種數(shù)字信號(hào)調(diào)制方式的識(shí)別具有重要的現(xiàn)實(shí)意義。

目前數(shù)字信號(hào)調(diào)制方式識(shí)別算法主要分為基于特征參數(shù)的方法[3]、基于信號(hào)高階累計(jì)量的方法[4]、基于時(shí)頻分析的方法以及基于小波變換的方法[5]等?；谔卣鲄?shù)的方法從時(shí)域提取特征參數(shù)，但受信噪比的影響較大[6]；基于信號(hào)高階累計(jì)量的方法雖對(duì)噪聲不敏感，但對(duì)載波和碼元同步要求較高；基于時(shí)頻分析的方法計(jì)算量較大，分類效果不太明顯；基于小波變換的方法對(duì)初始要求較低，易于實(shí)現(xiàn)，且分類效果較明顯。

本文通過(guò)幅度歸一化前后小波變換模值對(duì)數(shù)字信號(hào)進(jìn)行類間識(shí)別，在此基礎(chǔ)上根據(jù)幅度歸一化前數(shù)字信號(hào)小波變換模值的階梯幅度完成數(shù)字信號(hào)的類內(nèi)識(shí)別，最終完成數(shù)字信號(hào)的調(diào)制方式識(shí)別，并在Matlab7.1仿真平臺(tái)下對(duì)該方法的可行性進(jìn)行驗(yàn)證。

1 數(shù)字信號(hào)機(jī)理分析

常見(jiàn)的數(shù)字信號(hào)主要有振幅鍵控(ASK、QAM)、頻移鍵控(2FSK、4FSK)、相移鍵控(BPSK、QPSK)[7-8]。現(xiàn)對(duì)各類數(shù)字信號(hào)機(jī)理進(jìn)行分析。假設(shè)接收端的信號(hào)r(t)的表達(dá)式為

r(t)=s(t)+n(t)

(1)

式中,n(t)為信道的高斯白噪聲，s(t)是接收到的發(fā)射端產(chǎn)生的調(diào)制信號(hào)，不同調(diào)制方式的調(diào)制信號(hào)表達(dá)式如下：

(2)

(3)

(4)

(5)

式中：Ak為幅值；fk為載波頻率；fc為一固定載波頻率；θk為第k個(gè)碼元的初始相位；θc為一固定初始相位；Ts為碼元寬度。

2 調(diào)制方式識(shí)別

2.1 小波變換基本原理

小波變換是在窗口傅里葉變換基礎(chǔ)上發(fā)展起來(lái)的一種新的時(shí)頻分析方法，其窗口寬度可調(diào)，在時(shí)域、頻域都具有表征信號(hào)局部特征的能力，且具備多分辨分析的特點(diǎn)。一維信號(hào)f(t)的連續(xù)小波變換定義如下：

(6)

式中:a為尺度因子；b為定位參數(shù)；ψ(t)為小波母函數(shù)。在不同的應(yīng)用之中，選擇不同的小波母函數(shù)會(huì)產(chǎn)生不同的結(jié)果。為滿足調(diào)制方式識(shí)別的需求，選擇Haar小波作為小波母函數(shù)。Haar小波不僅簡(jiǎn)單，且對(duì)相位變化的暫態(tài)信號(hào)有較強(qiáng)的檢測(cè)能力。其表達(dá)式如下：

(7)

當(dāng)一維信號(hào)f(t)無(wú)瞬態(tài)變化時(shí)，其小波變換具有恒定的輸出，僅與尺度因子a有關(guān)，而與t無(wú)關(guān)；當(dāng)信號(hào)存在瞬態(tài)變化時(shí)，其小波變換則具有明顯的變化[3]。因此，Haar小波母函數(shù)可應(yīng)用于數(shù)字信號(hào)的調(diào)制方式識(shí)別。

2.2 小波變換模值分析

Haar小波在一個(gè)碼元的范圍之內(nèi)，各數(shù)字信號(hào)的Haar小波變換如下：

(8)

(9)

(10)

(11)

求Haar小波變換后的各數(shù)字信號(hào)模值，其表達(dá)式如下：

(12)

(13)

(14)

(15)

為消除幅度帶來(lái)的影響，對(duì)信號(hào)進(jìn)行幅度歸一化，其表達(dá)式如下：

(16)

對(duì)幅度歸一化后各數(shù)字信號(hào)計(jì)算小波變換模值，其表達(dá)式如下：

(17)

(18)

(19)

(20)

2.3 調(diào)制方式識(shí)別流程

幅度歸一化前2FSK、4FSK、MSK等頻移鍵控信號(hào)，小波變換模值取決于前后碼元的頻率，由于沒(méi)有幅度調(diào)制，幅度歸一化后調(diào)頻信號(hào)小波變換模值與調(diào)相信號(hào)一樣無(wú)明顯的變化，幅度歸一化前后小波變換模值都為多階梯函數(shù)；幅度歸一化前ASK、QAM等振幅鍵控信號(hào)，當(dāng)碼元沒(méi)有變化時(shí)小波變換模值是定值，碼元變化時(shí)小波變換模值取決于前后碼元的幅度且在碼元的交界處前后相位差較大，小波變換模值將產(chǎn)生跳變，其小波變換模值為多階梯函數(shù)，幅度歸一化后調(diào)幅信號(hào)的小波變換模值變?yōu)槎ㄖ?，只是在相位變化處產(chǎn)生跳變，此時(shí)小波變換模值近似為單值函數(shù);幅度歸一化前BPSK等相移鍵控信號(hào)，小波變換模值是定值，僅在相位變化處產(chǎn)生跳變，幅度歸一化后調(diào)相信號(hào)與幅度歸一化前的小波變換模值基本相同，其小波變換的模值均近似為單值函數(shù)。因此可以通過(guò)統(tǒng)計(jì)幅度歸一化前后小波變換模值進(jìn)行調(diào)制方式的類間識(shí)別。

完成類間調(diào)制方式的識(shí)別之后，進(jìn)一步對(duì)各類數(shù)字信號(hào)進(jìn)行類內(nèi)識(shí)別。以頻移鍵控信號(hào)為例進(jìn)行類內(nèi)識(shí)別，對(duì)于多進(jìn)制頻移鍵控信號(hào)MFSK存在M個(gè)載波頻率，幅度歸一化前MFSK信號(hào)小波變換模值將會(huì)產(chǎn)生M個(gè)不同的階梯幅度，可根據(jù)幅度歸一化前小波變換模值不同幅度的階梯數(shù)進(jìn)行MFSK調(diào)制信號(hào)的類內(nèi)識(shí)別。在高斯白噪聲的條件下且信噪比較高時(shí)，幅度歸一化前MFSK信號(hào)小波變換模值可以近似為高斯，其概率密度函數(shù)為M個(gè)高斯函數(shù)的和，因此，對(duì)幅度歸一化前MFSK信號(hào)小波變換模值進(jìn)行幅度直方圖的統(tǒng)計(jì)即可對(duì)MFSK進(jìn)行類內(nèi)識(shí)別，若幅度直方圖中的峰值個(gè)數(shù)在(M/2，M)范圍內(nèi)，可判定輸入的為MFSK信號(hào)。幅度鍵控和相移鍵控的類內(nèi)識(shí)別原理與頻移鍵控基本類似，不再贅述。

基于小波變換的調(diào)制方式識(shí)別流程如圖1所示。

3 仿真分析

為驗(yàn)證基于小波變換的數(shù)字信號(hào)調(diào)制識(shí)別方法的正確性及有效性，采用Matlab平臺(tái)進(jìn)行仿真，Matlab仿真軟件操作界面如圖2所示。仿真引入信號(hào)源采用隨機(jī)二進(jìn)制序列，載波頻率為200Hz，采樣率為2KHz，碼速率為10B的2ASK信號(hào)、QAM信號(hào)、BPSK信號(hào)、2FSK信號(hào)、4FSK信號(hào)及MSK信號(hào)對(duì)該方法的正確性和有效性進(jìn)行仿真分析。

信噪比為10dB條件下各調(diào)制方式幅度歸一化前后的小波變換模值如圖3所示。

由圖3a、3b可以看出，忽略尖峰帶來(lái)的影響，ASK、QAM等幅度鍵控信號(hào)，在幅度歸一化前其小波變換模值變化為多級(jí)函數(shù)，幅度歸一化后近似為單值函數(shù)；由圖3c可以看出，BPSK等相移鍵控信號(hào)，在幅度歸一化前后其小波變換模值均近似為單值函數(shù)；由圖3d、3e、3f可以看出，2FSK、4FSK、MSK等頻移鍵控信號(hào)，在幅度歸一化前后其小波變換模值均為多級(jí)函數(shù)。為進(jìn)行類間信號(hào)的識(shí)別，對(duì)各調(diào)制方式信號(hào)小波變換幅度歸一化前與幅度歸一化后的均值比進(jìn)行統(tǒng)計(jì)分析如表1所示。

圖1 基于小波變換的調(diào)制方式識(shí)別

圖2 Matlab軟件操作界面

由表1可看出，幅度歸一化前后2FSK、4FSK、MSK信號(hào)小波變換模值的均值比在0.95～1.05范圍內(nèi)，幅度歸一化前后2ASK、4ASK、QAM信號(hào)小波變換模值的均值比小于0.9，幅度歸一化前后BPSK、QPSK信號(hào)小波變換模值的均值比均大于1.5。因此,可利用幅度歸一化前后小波變換的模值進(jìn)行有效的類間識(shí)別。

幅度歸一化前2FSK與4FSK信號(hào)小波變換模值的階梯幅度統(tǒng)計(jì)結(jié)果如圖4所示。幅度歸一化前2ASK與4ASK信號(hào)小波變換模值的階梯幅度統(tǒng)計(jì)結(jié)果如圖5所示。幅度歸一化前BPSK與QPSK信號(hào)小波變換相差統(tǒng)計(jì)結(jié)果如圖6所示。

由圖4可以看出，幅度歸一化前2FSK與4FSK信號(hào)小波變換模值的階梯幅度統(tǒng)計(jì)結(jié)果，其中2FSK的階梯幅度近似為2,4FSK的階梯幅度近似為4。由圖5可以看出，幅度歸一化前2ASK與4ASK信號(hào)小波變換模值的階梯幅度統(tǒng)計(jì)結(jié)果，其中2ASK的階梯幅度近似為2,4ASK的階梯幅度近似為4。由圖6可以看出，幅度歸一化前BPSK與QPSK信號(hào)小波變換相差統(tǒng)計(jì)結(jié)果，其中BPSK的相差為1種,QPSK的相差為3種。因此利用此參數(shù)可將類內(nèi)數(shù)字信號(hào)中進(jìn)制數(shù)為2與進(jìn)制數(shù)大于2的數(shù)字信號(hào)有效地區(qū)分開(kāi)。

由上述的仿真結(jié)果可知，基于小波變換調(diào)制方式識(shí)別方法能準(zhǔn)確識(shí)別出BPSK、QPSK、2ASK、4ASK、2FSK、4FSK及 QAM信號(hào)。為進(jìn)一步驗(yàn)證基于小波變換的調(diào)制方式識(shí)別方法的有效性和準(zhǔn)確性，在信噪比為10dB條件下，采用基于小波變換的調(diào)制方式識(shí)別方法分別對(duì)BPSK信號(hào)等數(shù)字信號(hào)的調(diào)制方式進(jìn)行識(shí)別，仿真次數(shù)為100次，仿真結(jié)果如表2所示。

圖3 幅度歸一化前后的小波變換模值

表1 幅度歸一化前后均值對(duì)比

從表2可以看出，當(dāng)信噪比為10dB，仿真次數(shù)為100次時(shí)，對(duì)BPSK、QPSK、2ASK、4ASK、2FSK、4FSK及QAM信號(hào)的調(diào)制方式識(shí)別成功的次數(shù)分別為100、99、100、99、100、100、98次，識(shí)別成功率分別為100%、99%、100%、99%、100%、100%、98%。由此可知，基于小波變換的調(diào)制方式識(shí)別方法能有效識(shí)別出BPSK等數(shù)字調(diào)制信號(hào)，且具有較高的識(shí)別準(zhǔn)確率。

圖4 頻移鍵控信號(hào)小波變換模值幅度直方圖

圖5 幅度鍵控小波變換模值幅度直方圖

圖6 相移鍵控小波變換相差統(tǒng)計(jì)圖

表2 10dB條件下調(diào)制方式識(shí)別準(zhǔn)確率

4 結(jié)束語(yǔ)

針對(duì)數(shù)字信號(hào)的特點(diǎn)，采用基于小波變換的方法對(duì)數(shù)字信號(hào)進(jìn)行調(diào)制方式識(shí)別。首先根據(jù)幅度歸一化前后小波變換模值對(duì)數(shù)字信號(hào)進(jìn)行類間識(shí)別，并在此基礎(chǔ)上根據(jù)幅度歸一化前數(shù)字信號(hào)小波變換模值的階梯幅度判別進(jìn)制數(shù)進(jìn)行類內(nèi)識(shí)別，從而最終確定數(shù)字信號(hào)的調(diào)制方式。通過(guò)仿真結(jié)果驗(yàn)證表明，該方法可對(duì)數(shù)字信號(hào)進(jìn)行有效識(shí)別。

[1]李洋洋,馮永新.基于特征參數(shù)的BPSK信號(hào)調(diào)制識(shí)別方法研究[J].計(jì)算機(jī)仿真,2013,30(3):201-204.

[2]陳健,闊永紅,李建東，等.基于小波變換的數(shù)字調(diào)制信號(hào)識(shí)別方法的研究[J].電子與信息學(xué)報(bào),2008,28(11):2026～2029.

[3]薛磊,劉小秋.基于小波變換的數(shù)字通信信號(hào)識(shí)別[J].電訊技術(shù),2006,(3):52-56.

[4]Qian Bo,Li GuangShen,Feng Yongxin.A Parameters Estimating Method for BOC Modulation Signal[J].Advanced Materials Research,2014,(912-914):1108-1111.

[5]高杰.基于小波變換的MFSK_MQAM信號(hào)類內(nèi)識(shí)別法[J].軟件無(wú)線電,2009,39(8):26-28.

[6]任謙.基于軟件無(wú)線電的調(diào)制方式自動(dòng)識(shí)別[D].北京：北京工業(yè)大學(xué),2013.

[7]于顯耀.數(shù)字調(diào)制信號(hào)調(diào)制方式識(shí)別及參數(shù)估計(jì)研究[D].大連：大連理工大學(xué),2013.

[8]Qian Bo,Zhang Xiaoyu,Sun Shuang.A Parameters Estimation Method for FH Signal Based on SPWVD[J].Advanced Materials Research,2014,(912-914):1112-1115.

(責(zé)任編輯：馬金發(fā))

Research on Modulation Recognition Method of Digital Signals Based on Wavelet Transform

FAN Zhenzhen,ZHU Lizhong,ZHANG Xiaoyu,LI Ying

(Shenyang Ligong University,Shenyang 110159,China)

Modulation recognition of digital signals has much significance.According to characteristics of digital signals,modulation recognition between-classes is completed by the modulus,which is calculated by normalized and non-normalized amplitude through Haar wavelet transform.And modulation recognition within-classes is fulfilled by digital signals ary.And this method is validated through software simulating.Simulation results show that the recognition method based on wavelet transform can recognize digital signals correctly.

digital signals;Wavelet Transform;modulation recognition;normalized amplitude

2014-07-07

范真真(1988—),女，碩士研究生;通訊作者：朱立忠(1967—),男,教授，研究方向：自動(dòng)控制理論及應(yīng)用、智能控制、圖像處理與融合、數(shù)字信號(hào)處理技術(shù)、移動(dòng)通信技術(shù).

1003-1251(2015)03-0010-06

TN911

A