載波相偏估計算法研究

金慧琴，王正磊，朱安周

（1海軍航空工程學(xué)院電子信息工程系，山東煙臺264001；2.91697部隊山東青島266000）

根據(jù)實現(xiàn)結(jié)構(gòu)的不同，載波同步算法有兩種實現(xiàn)方式：閉環(huán)結(jié)構(gòu)和開環(huán)結(jié)構(gòu)[1]。閉環(huán)結(jié)構(gòu)精度高，實現(xiàn)簡單，但收斂慢；開環(huán)結(jié)構(gòu)雖然實現(xiàn)復(fù)雜，但其捕獲時間短。本文主要考慮基于開環(huán)結(jié)構(gòu)的載波同步算法。載波同步包括載波頻率同步和相位同步兩部分[2]，分別通過載波頻偏估計和校正、載波相偏估計和校正實現(xiàn)同頻和同相。短波信道屬于頻率選擇性衰落信道，關(guān)于此類衰落信道條件下的頻偏估計的研究很少，大部分研究都是針對高斯白噪聲信道的。這些研究的前提是，經(jīng)過均衡的信號已經(jīng)去除了碼間干擾，使得均衡后的信號近似符合高斯白噪聲條件[3]，此前提成立的條件是均衡算法不受頻偏的影響。否則，由于頻偏的影響，均衡效果很差，本文在研究載波相偏估計算法基礎(chǔ)上，給出改進的載波相偏估計算法。

1 載波相偏估計算法

利用載波頻偏估計算法估計出頻偏后，對接收信號進行頻偏補償，以消除頻偏的影響：

由于頻偏估計算法本身具有一定的誤差，使得估計頻偏與實際頻偏存在一定的誤差，稱為剩余頻偏[4]，它會產(chǎn)生一定的相位誤差；同時，發(fā)送端與接收端之間的不同步使得接收信號本身存在相位誤差，基于以上兩個原因經(jīng)過頻偏校正后的信號，會存在相位旋轉(zhuǎn)，需要載波相偏估計與校正技術(shù)來消除相位誤差的影響。

1.1 最大似然估計

進入載波相偏估計器之前，假定已經(jīng)完成定時同步、頻偏估計及校正。將總通信時間分成若干個短時間段，在每個短時間段內(nèi)，剩余頻偏相對于符號速率已經(jīng)非常小，由它引起的相位變化非常緩慢[5]，此時載波相位可表示為：

其中表示第i個時間段內(nèi)的平均相位，n(k)為噪聲，此時可以認為在短時間段內(nèi)不存在剩余頻偏，只有相偏。信號模型表示如下：

其中sk是發(fā)送符號，φ是需要估計的載波相偏；wk是均值為零、方差為2σ2的高斯白噪聲，實部和虛部相互獨立；N為觀察數(shù)據(jù)長度。觀察數(shù)據(jù)的聯(lián)合概率密度函數(shù)為：

簡化式（4），去除無關(guān)項：

MPSK調(diào)制信號的載波相偏盲估計式為：

若發(fā)送符號sk已知，則數(shù)據(jù)輔助的載波相偏估計式為：

1.2 相位跳變

當信號的傳輸時間較短，只進行一次相偏估計就可以。當信號傳輸時間較長時，再將剩余頻偏等效于相偏進行相偏估計是不合理的，誤差很大。此時可以將整個傳輸時間分成若干個短時間段，在每個短時間段內(nèi)，剩余頻偏的影響很小，對每段數(shù)據(jù)進行相偏估計[6]。但是這種處理方式會引起段與段之間的相位跳變，導(dǎo)致誤碼率的增大，影響整個系統(tǒng)的通信性能。

針對相位跳變問題，F(xiàn)itz給出了一種簡單方法檢測并消除這種跳變，稱為Fitz算法[7-8]，假定第i段傳輸信號內(nèi)的相位估計為φi，消除相位跳變后的相位為：

其中α≤1。當信噪比較高時，此算法能比較有效的消除相位跳變；但是當信噪比較低時，會頻繁地出現(xiàn)相位跳變[9-10]，此時該算法消除相位跳變的能力較差。

文獻[11]給出了另一種解決方案，稱為Zheng算法，此時表示為：

該算法同樣存在信噪比較小時會頻繁出現(xiàn)相位跳變的問題[12-13]。

1.3 改進的載波相偏估計算法

為了推導(dǎo)改進的載波相偏估計算法，首先對剩余頻偏導(dǎo)致的相偏及利用Fitz算法估計的相偏進行仿真，觀察相偏的特點。發(fā)送信號采用8PSK調(diào)制，信噪比為15 dB，剩余歸一化頻偏為0.000 2，數(shù)據(jù)長度為40段，每段由32個觀察數(shù)據(jù)組成。

圖1 每段估計的相偏曲線圖

從圖1中可以看出，相偏估計值隨機的分布在真實相偏值附近，它的回歸曲線是一條直線，這意味著相偏估計值可以近似為：

式中，φ(k)為根據(jù)第k段觀察數(shù)據(jù)估計的相偏，L為時間段數(shù)，φ(1),φ(2),…,φ(L)相互獨立，ξk是均值為零、方差為σ2的高斯噪聲，a和b為未知常數(shù)。如果可以得到a和b的值，就可以通過式（11）預(yù)計后面時間段的載波相偏值，以糾正噪聲引起的相位跳變[14]。

由于ξk滿足高斯分布，且φ(1),φ(2),…,φ(L)間相互獨立，得到φ(1),φ(2),…,φ(L)的聯(lián)合概率密度函數(shù)為：

通過最大似然法估計a和b，要使P(a,b)取最大值，只要下式最?。?/p>

將式（13）分別對a和b求偏導(dǎo)，并令其為零：

得方程組：

方程組的系數(shù)行列式為：

則方程組（16）存在唯一解：

將L段數(shù)據(jù)作為一組，根據(jù)每組數(shù)據(jù)都可以計算出一條回歸曲線[15-16]，設(shè)每組得到的回歸曲線為：

其中M表示分組數(shù)。

假設(shè)已經(jīng)得到第m組的回歸曲線，則可以根據(jù)回歸曲線預(yù)測第m+1組的載波相偏值：

另外根據(jù)第m+1組的數(shù)據(jù)計算得到的回歸曲線為：

這樣對于第m+1組數(shù)據(jù)來說，根據(jù)第m組和第m+1組的回歸曲線，對每一段數(shù)據(jù)，都可以得到這 2個相偏估計值，根據(jù)這2個相偏值可以對相位跳變進行修正：

綜上所述，本算法的基本流程為：

1）將信號分段，每N個信號序列組成1段，每L段信號組成1組，共有M組；

2）用Fitz算法計算第1組的載波相偏，根據(jù)式（18）得到回歸曲線，然后根據(jù)式（20）得到第1組的最終相偏估計值。

3）假設(shè)已得到第i組的最終相偏估計值，利用第i組的回歸曲線根據(jù)式（21）預(yù)測第i+1組的相偏估計值φ′m+1(k)；然后根據(jù)第i+1組的數(shù)據(jù)得到第i+1組的回歸曲線，根據(jù)式（22）得到另外一組相偏值m+1(k)；最后根據(jù)式（23）對相偏進行修正以消除相位跳變得到第i+1組的最終相偏值。

4）重復(fù)步驟3）直到信號結(jié)束。

2 性能仿真分析

下面通過仿真比較提出的新算法和Zheng算法、Fitz算法的性能。采用8PSK調(diào)制信號，總的觀察數(shù)據(jù)長度為1 280個符號，每個時間段包含32個符號，每20段為一組，符號速率fs為240 0baud，剩余載波頻偏為0.000 2fs。當SNR分別為15 dB、10 dB時3種算法的估計性能分別如圖2、圖3所示；當SNR為10 dB時存在相偏及相偏校正后的信號星座圖如圖4所示；利用3種相偏估計算法校正相偏后系統(tǒng)的誤碼率性能曲線如圖5所示。

圖2 SNR=15dB時性能比較

圖3 SNR=10dB時性能比較

圖4 信號星座圖

在圖5中的New指的是本文提出的新算法，可以看出，當信噪比較大時Zheng算法和Fitz算法都能夠有效地消除相位跳變的影響，但是信噪比較小時，這兩種算法都會出現(xiàn)相位跳變現(xiàn)象，嚴重影響信號解調(diào)結(jié)果，導(dǎo)致誤碼率的增高；而新算法不管信噪比大小，都能有效地消除相位跳變的影響，利用此算法校正相偏后的性能和無相偏時的性能基本接近。

3 結(jié)束語

考察了頻率選擇性衰落信道條件下的載波頻偏估計算法，由于估計頻偏和實際頻偏總會存在誤差，也就是存在剩余頻偏，當時間較短、剩余頻偏相對符號速率較小時，可以將剩余頻偏處理為相偏，通過考察估計相偏的特點提出了一種新的相偏估計算法，能夠有效的實現(xiàn)對相偏的估計與校正。

參考文獻：

[1]謝麗君，譚立志.基于FPGA的MSK調(diào)制器設(shè)計與實現(xiàn)[J].電子設(shè)計工程，2012，20（6）:16-18.

[2]孫錦華，韓會梅.低信噪比下時頻聯(lián)合的載波同步算法[J].西安交通大學(xué)學(xué)報，2015，49（2）:62-67.

[3]劉曉東，陳衛(wèi)東.一種寬范圍高精度的載波頻偏估計算法[J].無線電工程，2014，44（3）:43-45.

[4]聶少軍，何兵哲，王宏卓.一種新的大載波頻偏估計算法[J].電子設(shè)計工程，2014，22（1）:168-170.

[5]張帥，張曉林.數(shù)字電視大范圍載波頻偏快逮捕獲算法[J].北京航空航天大學(xué)學(xué)報，2013，39（2）:225-229.

[6]張毅，歐陽志新，鄧云凱，等.高性能無數(shù)據(jù)輔助QPSK頻偏估計新算法[J].西安電子科技大學(xué)學(xué)報，2013，40（2）:187-193.

[7]包超峰，王瑞峰，孔德龍.基于改進Fitz算法的多普勒頻移估計方法的研究[J].計算機測量與控制，2014，22（2）:611-613.

[8]張媛.突發(fā)數(shù)據(jù)傳輸中的同步技術(shù)研究[D].南京:南京理工大學(xué)，2012.

[9]楊晗竹，李廣俠，趙陸文，等.突發(fā)信號前向頻偏估計方法綜述[J].無線電通信技術(shù)，2015，41（6）:10-15.

[10]王芳，陳勇，葉志清.一種快速高精度的改進Fitz頻率估計算法[J].江西師范大學(xué)學(xué)報，2014，38（6）:605-609.

圖5 存在相偏及相偏校正后的系統(tǒng)誤碼率

[11]曹亞桃.頻偏估計算法研究及FPGA實現(xiàn)[D].天津:天津大學(xué)，2014.

[12]孟慶萍，周新力，田偉.基于FFT和長時延自相關(guān)函數(shù)的頻偏估計方法[J].計算機工程與設(shè)計，2013，34（3）:799-803.

[13]俞靚.基于相位差的頻偏估計技術(shù)研究[D].天津:天津理工大學(xué)，2013.

[14]孟慶萍，周新力，田偉.一種改進的Fitz載波頻偏估計方法及其仿真[J].電訊技術(shù)，2012，52（10）:130-132.

[15]譚曉衡，張毛.一種高精度的改進FFT頻偏估計算法[J].重慶理工大學(xué)學(xué)報，2010，1（7）:71-75.

[16]姚國義，李鑫，蘭瑞田.一種突發(fā)通信信號檢測與頻偏估計算法[J].無線電通信技術(shù)，2014，40（4）:22-23.