高動(dòng)態(tài)環(huán)境中多進(jìn)制擴(kuò)頻系統(tǒng)的載波同步技術(shù)研究

1、高動(dòng)態(tài)環(huán)境中多進(jìn)制擴(kuò)頻系統(tǒng)的載波同步技術(shù)研究杜常林何旭(電子科技大學(xué)通信抗干擾技術(shù)國(guó)家級(jí)重點(diǎn)實(shí)驗(yàn)室,成都 610054摘要:本文針對(duì)高動(dòng)態(tài)環(huán)境中的多進(jìn)制擴(kuò)頻系統(tǒng)的載波同步問(wèn)題進(jìn)行了分析,并在此基礎(chǔ)上提出了一套載波同步算法。計(jì)算機(jī)仿真結(jié)果表明,在一定信噪比下,這種算法可以實(shí)現(xiàn)高動(dòng)態(tài)環(huán)境中載波的有效捕獲和跟蹤和數(shù)據(jù)的相干解調(diào)。關(guān)鍵詞:高動(dòng)態(tài);多進(jìn)制擴(kuò)頻;載波同步Carrier synchronization for M-ary spreaded spectrum system inhigh dynamic environmentDu Changlin He Xu(National Key Lab

2、of Communication, UESTC, Chengdu, 610054Abstract:In this paper, the problem about carrier synchronization for M-ary spreaded spectrum system in high dynamic environment is analysized, and two algorithms are introduced for this difficulty. Its supported by the simulation results that the algorithms c

3、an achieve the carrier acquisition and tracking above certain threshold for SNR, and then coherent de-spreading can be carried out.Keywords:high dynamic; M-ary spreaded spectrum; carrier synchronization1 引言在高動(dòng)態(tài)環(huán)境中,由于機(jī)載和彈載擴(kuò)頻通信機(jī)的動(dòng)態(tài)變化范圍很大,具有很高的速度、加速度和加速度變化率,因此由這些高動(dòng)態(tài)載體發(fā)送的信號(hào)在接收端將產(chǎn)生很大的多普勒頻移、多普勒頻移的一次變化率和二次

4、變化率,對(duì)正確解調(diào)數(shù)據(jù)提出了挑戰(zhàn)12。國(guó)內(nèi)外專家和學(xué)者圍繞載波的捕獲和跟蹤進(jìn)行了大量的研究,但大多數(shù)是圍繞低動(dòng)態(tài)、低碼率的應(yīng)用4展開(kāi)的。本文針對(duì)高動(dòng)態(tài)環(huán)境中的多進(jìn)制擴(kuò)頻信號(hào)的特性,提出了先用FFT1技術(shù)捕獲載波頻率,再用擴(kuò)展卡爾曼濾波(EKF算法6跟蹤載波頻率的高速變化和相位。經(jīng)過(guò)計(jì)算機(jī)仿真,本文的同步算法可以實(shí)現(xiàn)高動(dòng)態(tài)環(huán)境中載波的有效捕獲和跟蹤,以實(shí)現(xiàn)數(shù)據(jù)的相干解調(diào)。2 載波頻率的快速捕獲2.1 FFT快速捕獲原理載波同步包括頻率捕獲和跟蹤兩個(gè)過(guò)程7。本文將載波的頻率捕獲包含在偽碼相位捕獲的過(guò)程中,使捕獲技術(shù)從傳統(tǒng)的偽碼相位和載波頻率的二維搜索過(guò)程變?yōu)閭未a相位的一維搜索過(guò)程。載波快速捕獲環(huán)路

5、的結(jié)構(gòu)框圖如圖1所示。接收機(jī)的中頻接收信號(hào)可以表示為:(cos (I d d r t A PN t t t n t =+(1其中I 表示接收的中頻信號(hào)的標(biāo)稱頻率,d 表示載波多普勒頻偏,為傳輸延時(shí),n (t表示以I 為中心的窄帶噪聲。PN (t表示擴(kuò)頻序列,同步頭所用的擴(kuò)頻序列為M 序列。 部分相關(guān)圖1 基于FFT 的頻率快速捕獲環(huán)路當(dāng)本地偽碼相位和接收信號(hào)偽碼相位誤差|/2c T <=時(shí),可認(rèn)為偽碼相位捕獲成功。在不考慮噪聲時(shí),下變頻至基帶后,經(jīng)累加器對(duì)X 點(diǎn)數(shù)據(jù)進(jìn)行累加,可以得到如下輸出信號(hào):(2(cos (2(I d d k d Q d d k d Y k AX G f t k Y

6、 k AX G f sin t k + (2其中|G (f |=(1/Xsin(X f T c /sin( f T c 是累加器的頻域傳遞函數(shù)。令Z (n =Y I (n +j Y Q (n ,對(duì)Z (n 序列進(jìn)行N 點(diǎn)復(fù)數(shù)FFT ,可以通過(guò)測(cè)量峰值譜線的位置來(lái)完成載波頻偏的初始捕獲。當(dāng)信號(hào)很弱時(shí),載波多普勒頻移的譜線可能被噪聲譜線所湮沒(méi)?？梢圆捎枚嘀芷趫D平均譜估計(jì)法,使有用信號(hào)能量增強(qiáng),以得到正確的譜峰位置。2.2 關(guān)于N 和X 取值的簡(jiǎn)單討論由累加器的頻域特性可知,X 越大,相關(guān)損失越大,頻移的整個(gè)搜索范圍越小。根據(jù)采樣定理,需滿足:X 1/(2T c |f dm |,其中f dm 表示最

7、大多普勒頻偏。增大N 和X 可以提高頻率分辨率,增大X 還有利于降低捕獲門(mén)限,增大N 又意味著FFT 運(yùn)算量的增加,因此存在一個(gè)N 和X 的優(yōu)化取值。3 載波跟蹤算法本次設(shè)計(jì)要求跟蹤算法能實(shí)現(xiàn)高動(dòng)態(tài)環(huán)境中的載波頻率和載波相位的聯(lián)合估計(jì),以實(shí)現(xiàn)相干解擴(kuò),卡爾曼濾波器可以實(shí)現(xiàn)幾個(gè)狀態(tài)的聯(lián)合估計(jì)?？紤]到觀測(cè)變量與狀態(tài)變量的非線性特性,因此本次設(shè)計(jì)在跟蹤時(shí)選擇了擴(kuò)展卡爾曼濾波器(EKF 算法。輸入信號(hào)可以寫(xiě)為如下的形式:cos (n I n Q k k n k I k k sin k k n k Q k +=+r (3其中(k為相位誤差,n (k表示相位噪聲n (t的采樣值,n I (k、n Q (k

8、分別為解擴(kuò)后同相支路和正交支路的輸出噪聲,可以視為相互獨(dú)立的白噪聲。t k 很小時(shí),(k可以近似地展開(kāi)成如下形式:2011(12(1(1(2k k k k f k t f k t + (4其中,t k 為EKF 算法的第k 個(gè)采樣間隔,f 0(k-1、f 1(k-1為接收載波與本地載波第k 個(gè)采樣的頻率偏差及其一次變化率。令狀態(tài)矢量x (k和觀測(cè)矢量h (k分別為:x (k= (k 0(k-1 1(kT , h (k=I(k Q(k T =cos (k sin (k T ,則可以得到如下的狀態(tài)方程和觀測(cè)方程:(k+1=(k+(k(k=(k+(kx x v r h n (5其中v (k、n (k

9、分別為狀態(tài)噪聲矢量和觀測(cè)噪聲矢量,在理想情況下均處理為白噪聲序列。如果采樣間隔取為常數(shù)T ,即t k =T ,則(5式中的狀態(tài)轉(zhuǎn)移矩陣可以寫(xiě)為:21/201001T T T =(6假設(shè)狀態(tài)噪聲為零均值白噪聲,則狀態(tài)噪聲方差矩陣Q 與頻率加速度狀態(tài)分量v 2(k的有如下關(guān)系:4322322286(8221T T T T E k k E k T T T T T =Q (7再令H T (k為觀測(cè)矢量h (k關(guān)于狀態(tài)變量的梯度函數(shù),即(-sin (-1cos (-1T T k k k k k =H ,則可以得到EKF 迭代算法的方程組: T 1T 2T (k k=(k k-1(k(k-(k k -1(

10、k k-1=(k k(k=(k(k k-1(k(k=(k k-1(k(k(k k=(k k-1(k(k(k k-1(k+1k=(k k+xx L r h x xx H P H R L P H P P L H P P P Q(8其中L (k、P (k|k、P (k+1|k通常分別被稱為增益矩陣、估計(jì)誤差相關(guān)矩陣、預(yù)測(cè)誤差相關(guān)矩陣,為略大于1的權(quán)重系數(shù),R 為觀測(cè)噪聲方差矩陣。由于(5式為近似處理,對(duì)t k 的選取應(yīng)充分考慮估計(jì)精度和算法運(yùn)算量的折中。在(4式中,假定有用信號(hào)的幅度為1,因此需要將解擴(kuò)輸出歸一化。4 計(jì)算機(jī)仿真結(jié)果4.1 FFT 捕獲技術(shù)仿真在仿真時(shí),載波標(biāo)稱頻率為f c =3GH

11、z ,R b =7.14Kb/s ,偽碼周期L =128,標(biāo)稱中頻頻率為f I =R c =0.914MHz ,接收信號(hào)r (k為過(guò)采樣中頻信號(hào),過(guò)采樣頻率為f s =6R c =5.484MHz ,后接單邊帶噪聲帶寬為2R c 的低通濾波器。此處SNR 值均指過(guò)采樣信號(hào)r (k上的信號(hào)與噪聲功率比。根據(jù)設(shè)計(jì)指標(biāo),發(fā)射機(jī)與接收機(jī)的相對(duì)最大運(yùn)動(dòng)速度為2Km/s ,則|f dm |=20KHz ,取X =16。FFT 頻率分辨率的選取必須小于后續(xù)EKF 跟蹤算法的同步跟蹤帶寬。高動(dòng)態(tài)環(huán)境的頻偏一次變化率可達(dá)1KHz/s ,因此不宜采用太長(zhǎng)的分析時(shí)間T p 。因此,本次FFT 的分辨率設(shè)計(jì)為452H

12、z ,對(duì)應(yīng)的FFT 點(diǎn)數(shù)為N =128。對(duì)FFT 譜分析輸出采用三周期平均,根據(jù)輸出來(lái)調(diào)整本地載波NCO 。在f dm =20KHz ,SNR =-16dB (解擴(kuò)前,偽碼相位誤差為0、T c /6、T c /3、T c /2時(shí)FFT 的平方幅度輸出如圖2所示。2004006008000501001500200400600800050100150020040060080005010002004006008000204060=Tc/3F F T 輸出F F T 輸出FFT 輸出序號(hào)FFT 輸出序號(hào)=0=Tc/6=Tc/2圖2 SNR =-16dB (解擴(kuò)前時(shí)FFT 平方幅度輸出圖2中V 表示接收

13、端的M 序列與本地M 序列的相位差。第一組FFT 輸出結(jié)果(128個(gè)點(diǎn)均為零,是因?yàn)槠錇槌跏紶顟B(tài)。第二組的譜峰位置均在n =167處,對(duì)應(yīng)頻率估計(jì)值為452Hz×(167-128=20.34KHz ,剩余頻差為-340Hz ,并將此頻率捕獲結(jié)果反饋至載波NCO 。在第三組的譜峰位置均在n =383處,對(duì)應(yīng)的頻率估計(jì)值為452Hz×(383-128×3=-452Hz ,剩余頻差為112Hz 。此后每組FFT 輸出的譜峰位置均在128的整數(shù)倍處,對(duì)應(yīng)與FFT 輸出的頻率零點(diǎn)。同時(shí)可以看出,隨著+減小,輸出譜峰的幅度將增加。在SNR =-16dB (解擴(kuò)前時(shí),可以完成可

14、靠的捕獲。4.2 接收同步頭信息階段的EKF 算法仿真本階段算法的期望跟蹤性能為:頻率估計(jì)誤差小于±15Hz,相位估計(jì)誤差小于±15º,保證在接收數(shù)據(jù)信息時(shí)能采用相干解擴(kuò)。在仿真驗(yàn)證中,取T b =1.3824×10-4s ,采樣間隔T=5t k =6.912×10-4s ,P (0=diag 3 5×103 5×105,E v 22(k=4.386×106,R =0.016I 2×2,算法啟動(dòng)時(shí)的初始頻差設(shè)定為300Hz 。載波NCO 更新周期為T(mén) 0=16T=1.106×10-2s 。根據(jù)S

15、imulink 仿真結(jié)果,SNR =13dB (解擴(kuò)后時(shí)的頻率跟蹤誤差曲線(采樣周期1.3824×10-4s 和相位跟蹤誤差曲線(采樣周期1.106×10-2s 如圖3所示。頻率估計(jì)誤差(H z 0=60º 時(shí)間(s00.511.5-50相位估計(jì)誤差(°時(shí)間(s圖3 SNR =13dB 時(shí)頻率跟蹤誤差和相位跟蹤誤差曲線由圖3可以看出,經(jīng)過(guò)0.1s (約9個(gè)載波NCO 更新周期以后,頻率和相位估計(jì)誤差均趨于穩(wěn)定,均小于±15Hz 和±15º,達(dá)到相位跟蹤的要求。當(dāng)初始相位位差偏離60º時(shí),頻率誤差趨于穩(wěn)定所需的時(shí)間更

16、短。由此可知,采用EKF 算法進(jìn)行相位捕獲所需的時(shí)間至少為0.1s 。4.3 傳數(shù)據(jù)信息時(shí)的載波跟蹤傳數(shù)據(jù)信息時(shí),數(shù)據(jù)率為R b =50Kb/s ,七個(gè)信息比特組成一個(gè)符號(hào)(symbol ,符號(hào)率為R sym =7.14Kb/s ,每個(gè)符號(hào)用128×128Hadamard 矩陣的行(即Walsh 序列作為擴(kuò)頻序列進(jìn)行直接擴(kuò)頻,擴(kuò)頻后的碼片速率為Rc=0.914Mchip/s ,采用相干解擴(kuò)。此處EKF 算法的采樣速率和載波NCO 更新周期不變。Ev 22(k=4.386×106,P (0=diag2.5×10-3 6×102 5×104,R =

17、0.025I 2×2,算法啟動(dòng)時(shí)的初始頻差設(shè)定為10Hz ?？紤]到算法的數(shù)值穩(wěn)定性,每2s 對(duì)估計(jì)誤差協(xié)方差矩陣P 重啟一次,即置為P (0。載體速度和載波頻率變化的動(dòng)態(tài)特性曲線如圖 4所示,SNR =13dB (指解擴(kuò)后時(shí)的頻率跟蹤誤差和相位跟蹤誤差曲線如圖5所示。對(duì)圖4的動(dòng)態(tài)特性進(jìn)行壓縮,可得到較低動(dòng)態(tài)的模型。表1列出了在不同最大加速度下的頻率估計(jì)和相位估計(jì)誤差的均方根,單位分別為赫茲、度,g 為重力加速度。表1 不同SNR 和最大加速度下的估計(jì)誤差均方根0g 4g7g10gSNR(dB頻差相差頻差相差頻差相差頻差相差13 1.5 3.0 3.4 4.4 3.0 4.1 2.6 4.4 14 1.3 2.6 3.2 3.7 2.7 3.6 2.4 4.0 15 1.3 2.5 3.2 3.6 2.7 3.5 2.4 3.8 16 1.32.23.1 3.3 2.6 3.3 2.3 3.7由表1可知,隨著SNR 提高,頻率估計(jì)誤差和相位估計(jì)誤差的均方根有所減小,誤符號(hào)率(SER 也有所降低。當(dāng)動(dòng)態(tài)降低時(shí),仿真結(jié)果表明,頻率估計(jì)誤差和相位估計(jì)誤差的均方根沒(méi)有太大