馬社祥,陳 明**,孟 鑫2,王俊峰
(1.天津理工大學(xué) 電氣電子工程學(xué)院,天津 300384;2.天津理工大學(xué) 海運(yùn)學(xué)院,天津 300384)
星載AIS信號(hào)的頻偏估計(jì)*
馬社祥1,陳 明**1,孟 鑫2,王俊峰1
(1.天津理工大學(xué) 電氣電子工程學(xué)院,天津 300384;2.天津理工大學(xué) 海運(yùn)學(xué)院,天津 300384)
星載自動(dòng)識(shí)別系統(tǒng)(AIS)中信號(hào)存在嚴(yán)重的多普勒頻偏,影響信號(hào)的正確解調(diào),而傳統(tǒng)的基于快速傅里葉變換(FFT)的頻偏估計(jì)算法其估計(jì)范圍無法滿足需要。為此,提出了一種改進(jìn)的自相關(guān)-FFT頻偏估計(jì)算法,通過構(gòu)造輔助函數(shù)以確定信號(hào)頻偏的正負(fù)號(hào),在此基礎(chǔ)上對(duì)基于FFT的頻偏估計(jì)算法進(jìn)行改進(jìn),從而擴(kuò)大頻偏估計(jì)范圍,實(shí)現(xiàn)對(duì)大的多普勒頻移的估計(jì)。仿真和實(shí)驗(yàn)結(jié)果均表明,改進(jìn)算法的頻偏估計(jì)范圍擴(kuò)大為原來的2倍,具有很高的估計(jì)精度,十分接近修正后的克拉美羅界,且算法簡(jiǎn)單、易于實(shí)現(xiàn)。
星載自動(dòng)識(shí)別系統(tǒng);頻偏估計(jì);大多普勒頻移;快速傅里葉變換
自動(dòng)識(shí)別系統(tǒng)(Automatic Identification System,AIS)是一種工作在甚高頻(Very High Frequency,VHF)頻段,采用自組織時(shí)分多址(Self-organized Time Division Multiple Access,SOTDMA)現(xiàn)代通信技術(shù)的廣播式自動(dòng)報(bào)告系統(tǒng)。AIS系統(tǒng)用于岸-船、船-岸及船-船間的通信,其工作半徑約為40 n mile[1]。為了能夠?qū)h(yuǎn)海甚至全球范圍內(nèi)船舶進(jìn)行有效監(jiān)測(cè),以衛(wèi)星平臺(tái)為依托,建立星載AIS系統(tǒng)是一個(gè)很好的選擇。
星載AIS系統(tǒng)由一顆或多顆低軌道衛(wèi)星構(gòu)成[2],AIS接收機(jī)搭載在衛(wèi)星上以實(shí)現(xiàn)對(duì)船舶AIS信號(hào)的接收、解調(diào)和解碼AIS報(bào)文。而衛(wèi)星軌道高度一般為600~1 000 km,且衛(wèi)星的運(yùn)行速度為7.5 km/s[3],因此星載AIS信號(hào)的最大多普勒頻移約為±4 kHz。由于太空環(huán)境復(fù)雜,信號(hào)傳輸衰減較大,使得接收信號(hào)的信噪比較低[4]。因此,為了減少頻偏對(duì)信號(hào)解調(diào)及解碼造成的影響,星載AIS接收機(jī)在進(jìn)行信號(hào)解調(diào)前必須進(jìn)行頻偏的估計(jì)與校正,以便獲得正確的碼元序列,從而保證對(duì)海上船舶的有效監(jiān)控。
對(duì)于頻偏估計(jì),文獻(xiàn)[5]基于最小二乘原理,通過計(jì)算相鄰信號(hào)間的相位差得到頻偏信息。該算法雖具有較大的頻偏估計(jì)范圍,但其估計(jì)性能受噪聲影響較大,在低信噪比下估計(jì)性能很差,無法滿足星載AIS系統(tǒng)的需要。文獻(xiàn)[6]利用AIS信號(hào)中的訓(xùn)練序列等已知信息消除接收信號(hào)的相位信息,在獲得頻偏主值和擴(kuò)展部分的基礎(chǔ)上,通過數(shù)據(jù)擬合的方法消除快速傅里葉變換(Fast Fourier Transform,FFT)運(yùn)算的柵欄效應(yīng),得到準(zhǔn)確的頻偏估計(jì)值,但算法的運(yùn)算量較大,且其估計(jì)精度依賴于較長(zhǎng)的數(shù)據(jù)長(zhǎng)度。1999年,Morelli等[7]利用已估計(jì)出的時(shí)延值確定合適的采樣點(diǎn),然后做該采樣點(diǎn)處的自相關(guān)運(yùn)算得到頻偏估計(jì)值,但該算法的頻偏估計(jì)精度較低,且估計(jì)范圍僅為碼元速率的四分之一。2003年,彭華等人[8]對(duì)接收信號(hào)做自相關(guān)運(yùn)算以獲取頻偏主值,然后利用離散傅里葉變換(Discrete Fourier Transform,DFT)對(duì)頻偏進(jìn)行修正,雖然該算法的頻偏估計(jì)精度較高,但估計(jì)范圍僅為碼元速率的四分之一,無法滿足星載AIS系統(tǒng)的需求。上述算法中,文獻(xiàn)[5-6]是數(shù)據(jù)輔助下的參數(shù)估計(jì)算法,雖然估計(jì)范圍可以滿足星載AIS系統(tǒng)的需要,但需要先驗(yàn)信息,而在實(shí)際應(yīng)用中很難獲得,因此其估計(jì)性能無法保證;文獻(xiàn)[7-8]是非數(shù)據(jù)輔助下的參數(shù)估計(jì)算法,雖不需先驗(yàn)信息,但其估計(jì)范圍十分有限,無法滿足星載AIS系統(tǒng)的需要。因此,研究出一種滿足星載AIS系統(tǒng)頻偏范圍的非數(shù)據(jù)輔助的頻偏估計(jì)算法將會(huì)有很好的應(yīng)用前景。
為了能夠估計(jì)星載AIS信號(hào)的頻偏,本文通過構(gòu)造輔助序列確定接收信號(hào)頻偏符號(hào),并與FFT算法相結(jié)合,進(jìn)而得到改進(jìn)算法。改進(jìn)算法的頻偏估計(jì)范圍從(-1/(4Tb),1/(4Tb))變?yōu)?-1/(2Tb),1/(2Tb)),Tb為碼元周期,使其能夠滿足星載AIS系統(tǒng)的需要,且提高了估計(jì)精度。在實(shí)驗(yàn)室環(huán)境下利用已有實(shí)驗(yàn)儀器對(duì)星載AIS信號(hào)進(jìn)行模擬,實(shí)驗(yàn)結(jié)果表明,改進(jìn)算法的估計(jì)精度較高,與理論仿真基本一致,具有很高的應(yīng)用價(jià)值。
星載AIS系統(tǒng)采用GMSK調(diào)制方式[9],基帶信號(hào)的復(fù)包絡(luò)可表示如下:
s(t)=ejφ(t,α),nTb≤t≤(n+1)Tb。
(1)
式中:φ(t,α)是信號(hào)的相位信息,可表示為
(2)
式中:α={αi}為信息序列,取值為{+1,-1};Tb是碼元周期;h是調(diào)制指數(shù),在GMSK調(diào)制中h=0.5;相位脈沖響應(yīng)
(3)
(4)
(5)
式(4)中:h(t)為高斯濾波器的矩形脈沖響應(yīng),其取值范圍為(0,L1Tb);L1為高斯濾波器持續(xù)碼元個(gè)數(shù);Bb為高斯濾波器的3 dB帶寬。
假設(shè)s(t)在AWGN信道中傳輸,星載AIS接收機(jī)接收到的信號(hào)為
x(t)=ej(2πfet+θ)s(t-τ)+n(t) 。
(6)
式中:fe為頻偏,θ為相移,τ為時(shí)延,η(t)為高斯白噪聲。
對(duì)接收信號(hào)采樣,可得
(7)
式中:T為采樣周期,且T=Tb/Ns,Ns是過采樣因子;將x(nT)簡(jiǎn)記為x(n)。
2003年,彭華等人[8]提出了一種利用相位展開和周期圖進(jìn)行頻偏估計(jì)的方法。該算法是非數(shù)據(jù)輔助下的頻偏估計(jì)算法,是對(duì)接收信號(hào)進(jìn)行處理,通過構(gòu)造輔助序列消除調(diào)制相位信息的影響,得到含有頻偏信息和相移信息的信號(hào)。
文獻(xiàn)[8]給出了FFT算法的頻偏估計(jì)的推導(dǎo)公式,即
(8)
式中:
(9)
(10)
(11)
z(n)=(-1)nx2(n) 。
(12)
式中:x(n)是式(7)得到的采樣信號(hào)(過采樣因子為Ns),m為延遲周期。
因此根據(jù)文獻(xiàn)[8]知,頻偏的估計(jì)范圍為
(13)
當(dāng)m較大時(shí),本算法的估計(jì)范圍約為(-1/(4Tb),1/(4Tb))。由于此算法的估計(jì)范圍無法滿足星載AIS信號(hào)的需要,因此對(duì)該算法進(jìn)行改進(jìn)。
仿真時(shí)均采用AWGN信道,碼元序列長(zhǎng)度N=256 bit,過采樣因子Ns=8,碼元周期Tb=(1/9 600) s,高斯濾波器持續(xù)碼元長(zhǎng)度L1=3,歸一化3 dB帶寬為BT=0.4。進(jìn)行了500次蒙特卡羅實(shí)驗(yàn),且本節(jié)和第4節(jié)中的各圖均是采用此條件仿出的。其中,仿真平臺(tái)選用CPU為AMD Athlon(tm)II X4 640,主頻為3.0 GHz,內(nèi)存為8 GB的電腦;仿真軟件采用MATLAB 7.11,仿真條件均參考ITU-R M.1371-4建議書設(shè)置的,且碼元序列長(zhǎng)度選用一個(gè)標(biāo)準(zhǔn)時(shí)隙的長(zhǎng)度。圖1是FFT算法在[-4 800 Hz,4 800 Hz]內(nèi)的頻偏估計(jì)仿真結(jié)果。從圖1可以看出,此算法的估計(jì)值關(guān)于原點(diǎn)近似中心對(duì)稱,且當(dāng)頻偏大于0而小于2 400 Hz時(shí),其能較準(zhǔn)確地估計(jì);當(dāng)頻偏大于2 400 Hz時(shí),其估計(jì)值小于0,無法估計(jì)出頻偏。而負(fù)頻偏與正頻偏的情況剛好相反,在小于-2 400 Hz時(shí),其估計(jì)值大于0;而大于-2 400 Hz時(shí),能較準(zhǔn)確地估計(jì)。如果我們能夠根據(jù)接收信號(hào)構(gòu)造一個(gè)輔助函數(shù),以明確頻偏的符號(hào),那么就能將原算法的估計(jì)范圍擴(kuò)大1倍。
圖1 FFT算法在[-4 800 Hz,4 800 Hz]內(nèi)的頻偏估計(jì)分布圖Fig.1 Frequency offset estimation of FFT algorithm in [-4 800 Hz,4 800 Hz]
為解決這一困難,由式(7)中的采樣信號(hào)x(n)構(gòu)造輔助序列y(n),即
y(n)=(-1)nx(n),n=0,1,…,L-1,
(14)
則y(n)的DFT為
(15)
式中:M為DFT變換的點(diǎn)數(shù),由頻域采樣定理可知,M≥L。
確定Y(k)的最大譜線的位置,記為最大譜線序列數(shù)[8]。Y(k)的最大譜線序列數(shù)分布圖如圖2所示。
圖2 不同頻偏的最大譜線序列數(shù)分布圖Fig.2 Distribution of the sequence number of the largest spectral line in different frequency offsets