基于Haar小波變換的碼元速率估計

61251部隊 安佰強

吉林大學計算機科學與技術(shù)學院 61251部隊 鄭 偉

解放軍電子工程學院 羅高健

61251部隊 陳曉輝

碼元速率是數(shù)字通信中最重要的參數(shù)之一，對碼元速率的準確估計在通信對抗中有著重要的意義，它是實施無線電通信干擾時選擇干擾樣式和干擾參數(shù)的重要依據(jù)，也是進行通信信號盲識別和盲解調(diào)的重要前提。本文提出了一種在沒有任何先驗知識的條件下，采用Haar小波對中頻信號進行碼元速率估計的算法。這種算法可以在較低的信噪比下對ASK、PSK和QAM調(diào)制信號進行有效的碼元速率估計，并且容易在FPGA上實現(xiàn)。

1.基本原理

假設(shè)經(jīng)過前端處理后的中頻信號為x(t)，載頻為fc，以采樣頻率fs對x(t)進行采樣，得到離散信號x(k)，設(shè)其數(shù)學模型為：

這里s(k)是調(diào)制信號，n(k)是高斯白噪聲， ωc=2πfc/fs是信號的數(shù)字頻率，θc為載波的初始相位。(k) =Akejφk代表碼元的幅度和相位信息，在同一個碼元內(nèi)，(k)的值不變。根據(jù)幅度Ak和相位φk取值范圍的不同，(k)可以分別表示MASK、MPSK和MQAM信號。

小波變換具有在時間和頻率上突出信號的局部特性的能力，因此可以利用小波變換提取信號碼元間的跳變信息。離散小波變換的公式為：

其中，a是尺度參數(shù)，n是時間平移參數(shù)，*表示復(fù)共軛，ψ是小波函數(shù)，本文選取離散Haar小波函數(shù)，并乘以一個與信號有關(guān)的系數(shù)x(n)。Haar小波函數(shù)為：

增加系數(shù)的Haar小波函數(shù)為：

式中尺度參數(shù)a取偶數(shù)。在不考慮噪聲的情況下，對信號s(k)進行小波變換：

(a,n)的突變特性與尺度參數(shù)a的選擇有著密切的關(guān)系，一般來說，小尺度下的小波變換能精確刻畫信號的邊緣信息，但受噪聲影響較大，不利于碼元速率的精確估計；大尺度下的小波變換對噪聲有較好的平抑作用，但容易淹沒真實局部極值，不能正確反映信號的局部特征。文獻[7][8]研究了如何選擇最優(yōu)尺度和如何進行多尺度小波變換，利用他們提出的方法，可以大幅提高碼元速率估計的準確性。為了便于分析，本文只討論單一尺度的小波變換，并假設(shè)尺度參數(shù)小于一個碼元周期，此時（5）式會出現(xiàn)兩種情況：

1）在尺度參數(shù)范圍內(nèi)無碼元變化，即當 -a/2 ≤k≤a/2- 1時，(k) =Aejφ，此時有：

2）在尺度參數(shù)范圍內(nèi)有碼元變化，即存在i滿足 -a/2 ＜i≤a/2- 1，當 -a/2≤k＜i時，有(k) =A0ejφ0，當i≤k≤a/2- 1時，有(k) =A1ejφ1。為方便計算，不妨設(shè)i=n，此時計算的結(jié)果對分析小波變換的突變特性無太大影響。此時有：

設(shè)PWT為的相位，則：

由（6）和（7）式可以看出，當無碼元變化時，PWT為常數(shù) - (ωc+ π)/2，當碼元發(fā)生變化時（無論是相位變化還是幅度變化），都會使PWT產(chǎn)生突變，因此，PWT具有與碼元相同的周期性。

在實際應(yīng)用中，由于信號的碼元不一定連續(xù)變化，可能會出現(xiàn)連續(xù)相同的碼元，因此直接從PWT中獲得碼元周期信息比較困難?？梢韵葘WT做FFT變換，從頻譜上估計碼元速率，具體方法見下節(jié)內(nèi)容。

2.算法仿真

用MATLAB工具對本文提出的算法進行仿真，以128QAM調(diào)制信號為例。產(chǎn)生載頻fc= 10KHz，采樣率fs= 100KHz，碼元速率fd=4KHz的128QAM數(shù)字通信信號，噪聲為高斯白噪聲，SNR=5dB，取N= 1024點計算小波變換（尺度參數(shù)為2）。

MATLAB產(chǎn)生的128QAM信號如圖1所示。

圖1 128QAM信號（I路和Q路）

利用式（5）獲得信號的小波變換，并計算其相位PWT，如圖2所示。

圖2 128QAM信號的PWT

對相位PWT做FFT變換，獲得相位的頻譜FWT。未經(jīng)處理的頻譜具有較大的零頻。去掉零頻，并設(shè)置門限β，將低于β的頻點濾除。β的取值可根據(jù)頻譜的最大值max(FWT)確定（不含零頻的最大值），在此設(shè) β= 0.8max(FWT)。未經(jīng)過處理的FWT和經(jīng)過處理的FWT見圖3。

圖3 128QAM信號PWT的頻譜

由圖3可以看出，第一個譜峰對應(yīng)的頻率即為碼元速率。也可以先計算出各譜峰間的距離，取其平均值來估計碼元速率，以提高估計的精度。設(shè)譜峰的平均間距為ΔN，則碼元速率可以表示為：

碼元速率的估計精度Δf與頻譜的頻率分辨率相同，即為 Δf=fs/N。

3.仿真結(jié)果

為了分析不同噪聲條件和調(diào)制方式對碼元速率估計性能的影響，我們在不同信噪比下，分別對4ASK、8ASK、BPSK、QPSK、64QAM和128QAM六種調(diào)制信號進行1000次碼元速率估計。設(shè)估計誤差在±Δf內(nèi)為準確估計，則當尺度參數(shù)a=2時，碼元速率估計的準確率如圖4所示。這六種通信信號的參數(shù)為載頻fc= 10KHz，采樣率fs= 100KHz，碼元速率fd=4KHz，碼元個數(shù)N=100。

圖4 碼元速率估計的準確率

從仿真結(jié)果可以看出，在信噪比不低于4dB時，此算法對碼元速率估計的準確率可以達到95%以上。值得指出的是，64QAM和128QAM的估計效果較好，其原因是QAM調(diào)制信號在碼元之間不僅有相位突變，同時還有幅度突變，包含更多的碼元突變信息。

4.結(jié)論

本文提出的利用Haar小波變換估計中頻信號碼元速率的算法，充分利用了碼元變化時的幅度突變和相位突變信息，可以在較低的信噪比下正確獲取MASK、MPSK和MQAM調(diào)制信號的碼元速率。在低尺度參數(shù)下，該算法可以較容易的在FPGA上實現(xiàn)，具有一定的實用價值。

[1] 鄧振淼,劉渝.基于多尺度Haar小波變換的MPSK信號碼速率盲估計[J].系統(tǒng)工程與電子技術(shù),2008,30(1):36-40.

[2] 許小東,路友榮,戴旭初,徐佩霞.一種基于最優(yōu)尺度小波變換的符號率估計算法[J].中國科學技術(shù)大學學報,2008,38(10):1147-1152.

[3] 秦永,萬國金,謝芳娟.基于小波變換的MPSK信號符號率估計方法[J].通信技術(shù),2009,42(1):81-83.

[4] 馮旭哲,楊俊,羅飛路.基于小波變換的通信信號碼元速率估計[J].系統(tǒng)仿真學報,2008,20(5):1259-1261.

[5] 張曉勇,羅來源,羅武忠.一種基于瞬時幅度的PSK信號符號速率估計方法[J].信號處理,2008,24(2):316-319.