王思秀 張蕾 任艷 馮長(zhǎng)征

摘 要：在高速移動(dòng)通信系統(tǒng)中，收發(fā)雙方往往面臨較大多普勒頻移和有限導(dǎo)頻開銷兩個(gè)問題，從而嚴(yán)重影響傳統(tǒng)載波同步模式（TCSP）的整體性能。為此，提出了一種基于自相關(guān)增量的載波參數(shù)解耦合技術(shù)（CPEDT），并應(yīng)用到傳統(tǒng)載波同步模式（CPEDTTCSP）中。首先在接收端任取一段導(dǎo)頻信號(hào)進(jìn)行去調(diào)制操作，然后對(duì)去調(diào)制信號(hào)進(jìn)行有效延遲長(zhǎng)度為α的相關(guān)運(yùn)算，一方面利用該相關(guān)運(yùn)算的結(jié)果作頻偏估計(jì)，另一方面利用取α為導(dǎo)頻長(zhǎng)度一半的相關(guān)運(yùn)算結(jié)果的共軛形式與去調(diào)制信號(hào)一同作最大似然相偏估計(jì)。理論分析和仿真結(jié)果表明，當(dāng)導(dǎo)頻初始位置為零時(shí)，提出的CPEDTTCSP可以實(shí)現(xiàn)TCSP中頻偏估計(jì)和相偏估計(jì)的解耦合，同時(shí)能夠使最大似然相偏估計(jì)的復(fù)乘運(yùn)算量由L降為1，從而更加適用于高速移動(dòng)通信。

關(guān)鍵詞：多普勒效應(yīng);數(shù)據(jù)輔助;最大似然準(zhǔn)則;載波同步;高速移動(dòng)通信

中圖分類號(hào)：TN911.23

文獻(xiàn)標(biāo)志碼：A

Carrier parameter decoupling technique based on autocorrelation increment

WANG Sixiu*， ZHANG Lei， REN Yan， FENG Changzheng

College of Computer Science and Engineering， Xinjiang University of Finance and Economics， Urumqi Xinjiang 830012， China

Abstract：

In the highspeed mobile communications， transceivers always face large Doppler shift and limited pilot overhead， which severely influence the overall performance of the Traditional Carrier Synchronization Pattern （TCSP）. Thus， an autocorrelation increment based Carrier Parameter Estimation Decoupling Technique （CPEDT） was proposed and was applied to the TCSP （CPEDTTCSP）. Firstly， a pilot signal with certain length was selected at the receiving end to perform the operation of modulation removal， and then the correlation operation with an effective delay length α was performed on the signal with modulation removal. The frequency offset was estimated by the result of the correlation operation， and the conjugate form of the correlation operation result with α as half of the pilot length was used to make the maximum likelihood phase offset estimation with the signal with modulation removal. Theoretical analysis and simulation results show that with pilot starting location of zero， the CPEDT TCSP can implement the decoupling between the frequency offset estimation and the phase offset estimation in the TCSP， and can reduce the computational complexity of complex multiplication from L to 1 in the maximum likelihood phase offset estimation， therefore is more suitable for highspeed mobile communications.

Key words：

Doppler effect; dataaided; maximum likelihood criterion; carrier synchronization; highspeed mobile communications

0?引言

在諸如衛(wèi)星通信和無人機(jī)通信等高速移動(dòng)通信中，收發(fā)端往往會(huì)面臨兩個(gè)不利因素：較大的多普勒頻移和有限的導(dǎo)頻資源，從而會(huì)嚴(yán)重影響接收端相干解調(diào)性能[1-4]， 這也正是引入載波同步技術(shù)的原因。對(duì)于傳統(tǒng)載波同步模式（Traditional Carrier Synchronization Pattern，TCSP），其處理流程可以等效成一個(gè)頻偏估計(jì)器級(jí)聯(lián)一個(gè)相偏估計(jì)器，如圖1所示， 即頻偏估計(jì)器決定了后續(xù)相偏估計(jì)器的性能。這是因?yàn)榻?jīng)補(bǔ)償器后的殘留頻偏就是一個(gè)時(shí)變的累積相偏。

針對(duì)上述通信場(chǎng)景存在的大多普勒頻移，很多文獻(xiàn)提出了不同的載波同步算法[5-8]。針對(duì)衛(wèi)星通信數(shù)據(jù)鏈下的低信噪比環(huán)境，文獻(xiàn)[5]提出了一種Turbo碼輔助時(shí)頻域聯(lián)合的載波同步算法;針對(duì)光通信，文獻(xiàn)[6]提出了一種基于期望最大準(zhǔn)則的載波同步算法。在導(dǎo)頻開銷足夠大的情況下，這兩種載波同步算法都可以獲得優(yōu)異的估計(jì)性能。但在諸如衛(wèi)星通信、無人機(jī)通信等高速移動(dòng)通信中，導(dǎo)頻資源是非常有限的， 這樣就很可能導(dǎo)致頻偏估計(jì)算法失準(zhǔn)，從而影響到后續(xù)相偏估計(jì)算法的性能。為此，從載波同步算法的估計(jì)性能限——克拉美羅界出發(fā)，文獻(xiàn)[7-8]將接收端采樣零時(shí)刻設(shè)置在數(shù)據(jù)幀結(jié)構(gòu)的正中間，發(fā)現(xiàn)了載波同步算法中頻偏估計(jì)和相偏估計(jì)的克拉美羅界的解耦合現(xiàn)象; 然而文獻(xiàn)中并沒有進(jìn)一步討論此時(shí)載波同步算法的性能。另外從實(shí)際應(yīng)用考慮，人為地設(shè)置采樣零時(shí)刻的位置是不可取的。

另一方面，載波同步算法的性能與是否使用數(shù)據(jù)輔助技術(shù)密切相關(guān)[9]。由此，載波同步算法可分為數(shù)據(jù)輔助型和非數(shù)據(jù)輔助型兩類，其中，前者利用了稱為導(dǎo)頻序列的數(shù)據(jù)信號(hào)通過一種線性去調(diào)制方式來獲得僅包含載波頻偏和相偏的單音信號(hào)[10-11];后者則利用了接收到的未知數(shù)據(jù)信息通過諸如冪次運(yùn)算等非線性去調(diào)制方式或者利用解調(diào)/譯碼軟信息的方法來獲得上述的單音信號(hào)[12-14]。顯然，由于采用了非線性變換，后者會(huì)具有較高的信噪比門限和復(fù)雜度。鑒于此，數(shù)據(jù)輔助型的載波同步算法被普遍應(yīng)用到高速移動(dòng)通信中。

考慮高速移動(dòng)通信存在的上述兩個(gè)不利因素，基于文獻(xiàn)[7-8]的啟發(fā)，本文提出了一種基于自相關(guān)增量的載波參數(shù)解耦合技術(shù)（Carrier Parameter Estimation Decoupling Technique，CPEDT），并將其應(yīng)用到了數(shù)據(jù)輔助型的傳統(tǒng)載波同步模式（CPEDTTCSP）中。該技術(shù)的基本原理如下：在接收端，首先利用接收到的導(dǎo)頻信號(hào)進(jìn)行線性去調(diào)制操作，得到一種重要的去調(diào)制信號(hào);然后利用自相關(guān)技術(shù)的思想對(duì)得到的去調(diào)制信號(hào)實(shí)現(xiàn)有效延遲長(zhǎng)度為α的自相關(guān)操作;最后利用得到的自相關(guān)增量進(jìn)行頻偏估計(jì)，同時(shí)利用設(shè)定α為去調(diào)制序列長(zhǎng)度一半的自相關(guān)增量的共軛形式連同去調(diào)制信號(hào)一起進(jìn)行最大似然相偏估計(jì)，這樣就實(shí)現(xiàn)了傳統(tǒng)載波同步模式中頻偏估計(jì)與相偏估計(jì)的解耦合。

1?系統(tǒng)模型

考慮加性高斯白噪聲信道下的具有理想定時(shí)[15]的單載波傳輸系統(tǒng)。接收端經(jīng)過匹配濾波和波特采樣后，經(jīng)過能量歸一化的等效基帶離散信號(hào)可以表示為：

r（k）=s（k）exp（j（2πfdTsk+θ））+n（k）; k∈κ（1）

其中： fd為由收發(fā)雙方相對(duì)運(yùn)動(dòng)而產(chǎn)生的最大多普勒頻移，即fd=fc·（v/c）（fc為載波頻率，v為視線方向上的相對(duì)速度，c為光速）;θ為由傳輸信道引入的隨機(jī)相偏，且在[-π， π]區(qū)間內(nèi)均勻分布;Ts為符號(hào)周期;s（k）表示能量歸一化的調(diào)制信號(hào);n（k）表示均值為0、實(shí)部和虛部方差均為N0/2的復(fù)高斯隨機(jī)變量;另外，κ{N，N+1，…，N+L-1}表示圖2所示的數(shù)據(jù)幀格式中導(dǎo)頻符號(hào)對(duì)應(yīng)的采樣時(shí)刻集，且令κ=L， 這里N和分別表示前導(dǎo)數(shù)據(jù)長(zhǎng)度（后文稱之為導(dǎo)頻起始位置）和包含諸如“數(shù)據(jù)導(dǎo)頻”結(jié)構(gòu)的總長(zhǎng)度。

為了實(shí)現(xiàn)數(shù)據(jù)輔助的估計(jì)模式，需要獲得無調(diào)制信息的信號(hào)，即對(duì)式（1）的兩邊同時(shí)乘以s（k）（表示取復(fù)數(shù)共軛），可得一個(gè)線性去調(diào)制信號(hào)z（k）：

z（k）r（k）s（k）=exp（j（2πfdTsk+θ））+d（k）（2）

其中d（k）n（k）s（k）仍為白噪聲。

2?載波參數(shù)估計(jì)解耦合技術(shù)的提出

為了引入載波參數(shù)估計(jì)解耦合技術(shù)（CPEDT），首先給出基于去調(diào)制信號(hào)的自相關(guān)公式，即：

R（α）=1L-α∑N+L-α-1k=Nz（k）z（k+α）=

1L-α∑L-α-1k=0z（k+N）z（k+α+N）=

exp（j2πfdTsα）+1L-α ψ（α）; 1≤α

其中，α為有效延遲長(zhǎng)度，ψ（α）為噪聲累加項(xiàng)，其具體表達(dá)形式如下：

ψ（α）∑N+L-α-1k=Nd（k+α）exp（-j2πfdTs（k））+d（k）·exp（j2πfdTs（k+α））+d（k）d（k+α））

接下來，基于自相關(guān)增量R（α），圖3給出了CPEDT的原理框圖。

下面給出CPEDT的實(shí)現(xiàn)過程：將得到的自相關(guān)增量R（α）與去調(diào)制信號(hào)z（k）一起送到補(bǔ)償器中，并且令α=（L-1）/2[7]，從而可以得到載波參數(shù)估計(jì)解耦合因子D（）：

D（）z（k）·（R（α）α=L-12）*=z（k）·R*（）=z（k）·[e-j2πfdTs× + ψ*（）]（4）

最后將得到的自相關(guān)增量R（α）和載波參數(shù)估計(jì)解耦合因子D（）一同送至TCSP的頻偏估計(jì)器和相偏估計(jì)器中。顯然，如圖3所示，這種設(shè)計(jì)可以實(shí)現(xiàn)整個(gè)載波同步的過程，且不需要使用其他頻偏估計(jì)算法。

3?載波參數(shù)估計(jì)解耦合技術(shù)的應(yīng)用

下面本文將提出的CPEDT應(yīng)用到TCSP中，便可得到載波參數(shù)估計(jì)解耦合基的傳統(tǒng)載波同步模式（記作CPEDTTCSP），如圖4所示，實(shí)現(xiàn)了串行處理向并行處理“轉(zhuǎn)換”。

這里以最大似然準(zhǔn)則為例來闡述CPEDT的有效性。假設(shè)TCSP采用了基于自相關(guān)增量的頻偏估計(jì)算法和基于最大似然準(zhǔn)則的相偏估計(jì)算法。

3.1?TCSP

首先利用基于自相關(guān)增量的頻偏估計(jì)算法獲得一個(gè)頻偏的估計(jì)值，記為d;然后將該估計(jì)值d與去調(diào)制信號(hào)z（k）一起送到補(bǔ)償器中，并將補(bǔ)償器的輸出信號(hào)送到基于最大似然準(zhǔn)則的相偏估計(jì)器中，可得：

s=arg∑N+L-1k=Nz（k）·exp（-j2πdTsk）=

argN+L-1k=N[z（k）·exp（-j2πdTsk）]L次復(fù)乘（5）

s=argexp（jθ）∑N+L-1k=Nexp（j2π（fd-d）Tsk）+=arg{exp（jθ）exp（j2π（fd-d）Ts（N+L-1））·

sin[π（fd-d）TsL]sin[π（fd-d）Ts]+}≈

θ，（fd-d）Ts≈0

H（θ， fd， d，N，L），（fd-d）Ts0 （6）

其中：H（θ， fd， d，N，L）為一個(gè)表征相位模糊的函數(shù)，其取值與實(shí)際相偏、頻偏、頻偏估計(jì)值和導(dǎo)頻起始位置以及導(dǎo)頻長(zhǎng)度有關(guān);亦為噪聲累加項(xiàng)，具有不同于的表達(dá)形式：

∑N+L-1k=Nd（k）exp（-j2πdTsk）

同理，當(dāng)信噪比很高時(shí)，有||≈0，那么式（6）也成立。從式（6）的結(jié)果可以看出，當(dāng)且僅當(dāng)剩余頻偏|（fd-d）Ts|≈0，相偏估計(jì)值才近似等于其真實(shí)值。但在諸如衛(wèi)星通信等通信系統(tǒng)中導(dǎo)頻資源是非常有限的，那么頻偏估計(jì)值就很可能會(huì)遠(yuǎn)離其真實(shí)值，從而導(dǎo)致后續(xù)的相位估計(jì)產(chǎn)生相位模糊問題。

EDTTCSP

將載波參數(shù)估計(jì)解耦合因子D（）送到基于最大似然準(zhǔn)則的相偏估計(jì)器中，可得

p=arg∑N+L-1k=ND=arg{∑N+L-1k=Nz（k）·

exp-j2πfdTs×L-12+ψL-12}=

argN+L-1k=Nz（k）·[exp（-jπfdTs（L-1））+′]1次復(fù)乘（7）

p=argexp（jθ）∑N+L-1k=Nexp（-jπfdTs（L-1））+′=

argexp（jθ）exp（j2πfdTsN）sin（πfdTsL）sin（πfdTs）+′≈

θ，N=0且fdTs≤1/L

H（θ， fd，N，L），N≠0或fdTs>1/L （8）

其中：H（θ， fd，N，L）也是一個(gè)表征相位模糊的函數(shù)，其取值與實(shí)際相偏、頻偏和導(dǎo)頻初始位置以及導(dǎo)頻長(zhǎng)度有關(guān);′亦為噪聲累加項(xiàng)，且具有如下的形式：

′∑N+L-1k=Nd（k）exp-j2πfdTs×L-12+

d（k）ψ（L-12）+ψL-12exp（j（2πfdTsk+θ））

當(dāng)信噪比較高時(shí)，有ψ′≈0，此時(shí)式（8）成立。由式（8）可以發(fā)現(xiàn)，當(dāng)N=0且|fdTs|≤1/L時(shí)，即使存在較大的頻偏，相偏估計(jì)值仍近似等于其真實(shí)值。換句話說，傳統(tǒng)載波同步模式中存在的相位模糊問題就可以一定程度上避免了。另外，比較式（5）和式（7）可知，經(jīng)過CPEDT處理后，TCSP中最大似然相偏估計(jì)的復(fù)乘運(yùn)算量由需要L次降至1次，顯然得到了較大的復(fù)雜度改善。

4?仿真結(jié)果與分析

經(jīng)過CPEDT處理后，解決了TCSP中存在的“頻偏估計(jì)直接影響相偏估計(jì)”問題，實(shí)現(xiàn)了TCSP的并行處理，同時(shí)還降低了TCSP中最大似然相偏估計(jì)的計(jì)算復(fù)雜度。值得一提的是：在有無使用CPEDTTCSP中，頻偏估計(jì)算法可以是一樣的， 因此，評(píng)估有無使用CPEDTTCSP的性能可等效為評(píng)估兩者最大似然相偏估計(jì)的性能（即式（6）和式（8））。不失一般性，仿真中調(diào)制方式為最小頻移鍵控（MinimumShift Keying， MSK）調(diào)制，使用導(dǎo)頻長(zhǎng)度L=10，導(dǎo)頻起始位置N≥0（需要說明的是N并不是導(dǎo)頻長(zhǎng)度，而是數(shù)據(jù)長(zhǎng)度）。由式（8）可知，最大似然相偏估計(jì)器可以抵抗的歸一化頻偏范圍為fdTs≤1/10。

1）導(dǎo)頻起始位置選?。杭僭O(shè)工作信噪比Eb/N0分別為4dB、6dB和8dB，歸一化頻偏fdTs=0.08<0.1、相偏θ=π/2=90°。圖5給出了不同導(dǎo)頻起始位置N下的相偏估計(jì)均方誤差曲線。

由圖5的仿真結(jié)果可以看出，不同信噪比下，即使設(shè)置不同的導(dǎo)頻初始位置N，都可以得到類似的結(jié)果，即當(dāng)導(dǎo)頻初始長(zhǎng)度N=0時(shí)，相偏估計(jì)性能最好;但隨著導(dǎo)頻初始位置N的增大，其性能會(huì)急劇惡化。為了在大頻偏下獲得盡可能好的估計(jì)性能，導(dǎo)頻初始位置N應(yīng)設(shè)置為0。此時(shí)圖2所示的數(shù)據(jù)幀格式可以看成衛(wèi)星通信等通信系統(tǒng)所用的DVBT2數(shù)據(jù)幀格式[16]。在后文的仿真參數(shù)中，設(shè)置N=0。

2）相偏估計(jì)期望性能：假設(shè)工作信噪比Eb/N0=8dB，歸一化頻偏fdTs=0.08，0.008，0.0008，相偏θ∈[-π，π]≈[-3.14，3.14]rad意味著調(diào)制信號(hào)按逆時(shí)針旋轉(zhuǎn)和順時(shí)針旋轉(zhuǎn)的最大角度為π。圖6給出了不同相偏下，有無使用CPEDTTCSP對(duì)應(yīng)的相偏估計(jì)期望性能。

從圖6的仿真結(jié)果可以發(fā)現(xiàn)，在不同的相偏范圍內(nèi)，即使所加歸一化頻偏較大時(shí)，使用了CPEDTTCSP中相偏估計(jì)期望值與其真實(shí)值幾乎完全重合，而TCSP中相偏估計(jì)期望值相差很大。但隨著所加頻偏的減小，比如當(dāng)歸一化頻偏為0.000-8 時(shí)，TCSP中的相偏估計(jì)期望值與其真實(shí)值也相差無幾。因此可以得到如下猜想：對(duì)于較大的頻偏，CPEDT就會(huì)顯著改善TCSP中的相偏估計(jì)性能（具體仿真結(jié)果見圖7）。

3）相偏估計(jì)的抗頻偏性能：假設(shè)工作信噪比Eb/N0=8dB，相偏θ=π/2=90°，歸一化頻偏fdTs∈[-0.1，0.1]。圖7給出了有無使用CPEDTTCSP中相偏估計(jì)的抗頻偏能力曲線。

從圖7的仿真結(jié)果可知，對(duì)于一定的頻偏范圍，TCSP中相偏估計(jì)性能會(huì)受到所加頻偏的嚴(yán)重影響，但使用了CPEDTTCSP中相偏估計(jì)可以在較大的頻偏下仍能正常工作即完成對(duì)相偏的準(zhǔn)確估計(jì)。這是因?yàn)镃PEDT能夠一定程度上補(bǔ)償由殘留頻偏引起的時(shí)變累積相偏，從而顯著降低頻偏估計(jì)對(duì)相偏估計(jì)性能的影響。

5?結(jié)語(yǔ)

針對(duì)高速移動(dòng)通信的特點(diǎn)及其存在的載波同步問題，提出了一種基于自相關(guān)增量的載波參數(shù)解耦合技術(shù)（CPEDT），并將其應(yīng)用到傳統(tǒng)載波同步模式（TCSP）中，得到了一種載波參數(shù)解耦合技術(shù)基的傳統(tǒng)載波同步模式（CPEDTTCSP）。理論分析和仿真結(jié)果都表明：將導(dǎo)頻初始位置設(shè)置為零，可以降低頻偏估計(jì)對(duì)相偏估計(jì)的初始影響。在此基礎(chǔ)上，CPEDTTCSP可以實(shí)現(xiàn)相偏估計(jì)與頻偏估計(jì)的解耦合;且與TCSP中的相偏估計(jì)相比，CPEDTTCSP中的相偏估計(jì)具有更好的估計(jì)性能、更低的計(jì)算復(fù)雜蘇和更高的抗頻偏能力。此外，所提的CPEDT還適用于其他的載波同步方算法，具有一定的普適性。

參考文獻(xiàn) （References）

[1]ANTN A， GARC?AROJO I， GIRN A， et al. Analysis of a distributed array system for satellite acquisition[J]. IEEE Transactions on Aerospace and Electronic Systems， 2017， 53（3）： 1158-1168.

[2]郭業(yè)才，趙衛(wèi)娟，張秀再. 衛(wèi)星信道三狀態(tài)Markov模型設(shè)計(jì)與仿真[J]. 計(jì)算機(jī)應(yīng)用， 2016， 36（1）： 91-95. （GUO Y C， ZHAO W J， ZHANG X Z. Design and simulation of satellite channel threestate Markov model[J]. Journal of Computer Applications， 2016， 36（1）： 91-95.）

[3]XIAO Z， XIA P， XIA X. Enabling UAV cellular with millimeterwave communication： potentials and approaches[J]. IEEE Communication Magazine， 2016， 54（5）： 66-73.

[4]徐進(jìn)，劉重陽(yáng)，廖文. 無人機(jī)測(cè)控系統(tǒng)抗干擾性的優(yōu)化[J]. 國(guó)外電子測(cè)量技術(shù)， 2017， 36（8）： 26-29. （XU J， LIU C Y， LIAO W. Optimization of UAV control system antijamming[J]. Foreign Electronic Measurement Technology， 2017， 36（8）： 26-29.）

[5]SUN J， YU Z. Joint timefrequency domain and codeaided carrier synchronization algorithm at low SNR[C]// Proceedings of the IEEE 7th Joint International Information Technology and Artificial Intelligence Conference. Piscataway： IEEE， 2014： 545-549.

[6]ZIBAR D， de CARVALHO L H， ESTARAN J， et al. Joint iterative carrier synchronization and signal detection employing expectation maximization[J]. Journal of Lightwave Technology， 2014， 32（8）： 1608-1615.

[7]RICE F. Carrierphase and frequencyestimation bounds for transmissions with embedded reference symbols[J]. IEEE Transactions on Communications， 2006， 54（2）： 221-225.

[8]YU Z， SUN J， BAI B， et al. A phase incrementbased frequency estimator for general PSAM in burst communications[C]// Proceedings of the IEEE 83rd Vehicular Technology Conference. Piscataway： IEEE， 2016： 1-5.

[9]凌琪琪，羅志年. 基于數(shù)據(jù)相關(guān)疊加訓(xùn)練序列的載波頻偏與信道聯(lián)合估計(jì)[J]. 計(jì)算機(jī)工程， 2017， 43（9）： 138-143. （LING Q Q， LUO Z N. Joint estimation of carrier frequency offset and channel based on data dependent superimposed training sequence [J]. Computer Engineering， 2017， 43（9）： 138-143.）

[10]ZHANG W， YANG L， CHENG X， et al. On the optimality of dataaided coarse timing with dirty templates[J]. IEEE Transactions on Vehicular Technology， 2014， 63（4）： 1759-1769.

[11]GOMAA A， ALDHAHIR N. Phase noise in asynchronous SCFDMA systems： performance analysis and dataaided compensation[J]. IEEE Transactions on Vehicular Technology， 2014， 63（6）： 2642-2652.

[12]SINGH P， KIM B W， JUNG S Y. Preamblebased synchronization scheme for electromagnetic wireless nanocommunications[J]. IET Communications， 2017， 11（7）： 1097-1105.

[13]WANG Q， HUANG B， XU Z. Nondataaided timing recovery algorithm for MILSTD SOQPSK[J]. Electronics Letters， 2015， 51（5）： 423-425.

[14]宋靖，張劍云，鄭志東，等. 分布式全相參雷達(dá)相干參數(shù)估計(jì)性[J]. 電子與信息學(xué)報(bào)， 2014， 36（8）： 1926-1931. （SONG J， ZHANG J Y， ZHENG Z D， et al. Coherent parameters estimation performance for distributed aperture coherent radar[J]. Journal of Electronics and Information Technology， 2014， 36（8）： 1926-1931.）