胡 輝，雷明東，李 林，方 玲，徐 健

（華東交通大學(xué)信息工程學(xué)院，江西 南昌 330013）

0 引言

人們常將應(yīng)用于導(dǎo)彈、火箭、衛(wèi)星、航天飛船等航天軍工領(lǐng)域的接收機定義為高動態(tài)接收機。在高動態(tài)環(huán)境下，接收機接收到的信號會具有較高的多普勒及多普勒頻率變化率，這會使GPS接收機對信號的跟蹤產(chǎn)生較大的困難［1］。由于載體的動態(tài)在偽碼跟蹤環(huán)路上引起的多普勒可以通過載波環(huán)的輔助來進行消除，所以對于高動態(tài)接收機來說，設(shè)計一種性能優(yōu)越的載波跟蹤環(huán)是至關(guān)重要的。

常用的載波跟蹤環(huán)有鎖頻環(huán)（FLL）和鎖相環(huán)（PLL），其中FLL 跟蹤的動態(tài)范圍較大，但跟蹤精確度不高；而PLL，跟蹤精度較高，但跟蹤的動態(tài)范圍較?。?］。針對上述問題，文獻［3-6］提出了先利用FLL將頻率牽引到PLL 的跟蹤門限內(nèi)，然后再切換到PLL跟蹤的模式，該類復(fù)合環(huán)較好地解決了對動態(tài)的適應(yīng)性問題，但存在牽引時間不易確定、環(huán)路切換頻繁、強制切換造成濾波不連續(xù)等缺陷。為了兼顧PLL 和FLL 的特性，文獻［2］提出了FLL輔助PLL 的復(fù)合環(huán)路結(jié)構(gòu)。在此基礎(chǔ)上，文獻［7-8］為了更好地適應(yīng)動態(tài)性和加快收斂速度，在文獻［2］的基礎(chǔ)上增加了FLL 環(huán)路來進行頻率牽引，但是該復(fù)合環(huán)路也存在切換門限不易確定、環(huán)路切換頻繁、濾波不連續(xù)等問題。

針對上述問題，本文提出了一種基于頻率牽引和自適應(yīng)加權(quán)技術(shù)的載波跟蹤組合環(huán)路控制策略。

1 傳統(tǒng)GPS信號跟蹤技術(shù)

GPS接收機一般載波跟蹤環(huán)如圖1所示，在完成對信號的下變頻、濾波、數(shù)字量化后，將中頻數(shù)字信號與本地復(fù)現(xiàn)的C／A 碼、載波做相關(guān)，然后經(jīng)過積分-清零器形成Ip（k）、Qp（k）兩路值。環(huán)路鑒別器利用Ip（k）、Qp（k）得到載波誤差值，將該誤差值濾波去噪后，加載到本地載波NCO 上，產(chǎn)生新的載波信號，從而保持對輸入信號的跟蹤。

圖1 一般GPS載波跟蹤環(huán)路Fig.1 General GPS carrier tracking loop

文中的Ip（k）和Qp（k）兩路信號模型的表達式如下所示：

其中A 為信號的幅度，D（k）為導(dǎo)航數(shù)據(jù)位，Δfd（k）為多普勒殘差，ε（k）為碼相位偏差，T 為預(yù)檢積分時間，）和）分別為引進的同相和正交噪聲分量。

而載波跟蹤環(huán)路根據(jù)選取的鑒別器的不同，常用的載波跟蹤環(huán)路有FLL和PLL。當(dāng)環(huán)路的鑒別器被確定后，主要分析的是環(huán)路參數(shù)對于跟蹤環(huán)路的影響。

1.1 鎖相環(huán)跟蹤技術(shù)

由于導(dǎo)航數(shù)據(jù)位中存在數(shù)據(jù)位的跳變PLL 應(yīng)選取對導(dǎo)航數(shù)據(jù)位不敏感的Costas環(huán)。為了更好地適應(yīng)動態(tài)性，本文特選取鑒相范圍廣的二象限反正切鑒相器，其相關(guān)公式和特性曲線參見文獻［10］。PLL環(huán)路經(jīng)驗門限公式可以表示如下：

式中：

其中σtPLL表示由熱噪聲引起的相位抖動；θerror表示動態(tài)應(yīng)力引起的相位抖動；σv表示由振動引起的振動器顫動；θA表示由阿倫偏差引起的振蕩器顫動；Bn為PLL環(huán)路的噪聲帶寬，ω0為環(huán)路濾波器的自然圓頻率；T 為預(yù)檢積分時間；C／N0為載噪比，由式（4）可知載噪比越高，熱噪聲誤差越小，那么環(huán)路能適應(yīng)的動態(tài)應(yīng)力也相應(yīng)的提高；表示為視線方向上的動態(tài)應(yīng)力，則d代表最大視線方向上的加加速度應(yīng)力，其計算公式如下所示：

其中a′為加加速度，fL1為載波的頻率，c為光速。為了更好地適應(yīng)高動態(tài)場景，本文預(yù)檢積分時間取值1ms。在對式（3）求導(dǎo)后，可以得到3階PLL環(huán)路最佳噪聲帶寬滿足的表達式為［6］：

從式（7）可知，在取環(huán)路帶寬的時候，要折衷地考慮動態(tài)應(yīng)力、熱噪聲誤差、預(yù)檢積分時間之間的關(guān)系。

1.2 鎖頻環(huán)跟蹤技術(shù)

在GPS跟蹤環(huán)中常用的鑒頻器主要有3種，其鑒頻公式及特性曲線見文獻［10］。而為了有效地應(yīng)對加加速度帶來的沖擊，本文復(fù)合環(huán)特選取鑒頻范圍廣的四象限反正切鑒頻器。FLL 環(huán)跟蹤誤差的一般經(jīng)驗門限取值為：

式中3σtFLL表示FLL熱噪聲頻率顫動，ferror表示FLL環(huán)的動態(tài)應(yīng)力誤差，相關(guān)的表達式參見文獻［10］。從式（8）可知，減小預(yù)檢積分時間T，可以相應(yīng)的增加動態(tài)應(yīng)力的門限，提高環(huán)路對動態(tài)的適應(yīng)能力。在T 取值1ms的情況下，2階FLL環(huán)路噪聲帶寬、載噪比、動態(tài)應(yīng)力三者滿足的表達式如下所示：

從上面的分析可以看出，對于FLL和PLL來說，影響其跟蹤性能的主要因素為熱噪聲和載體的動態(tài)性，所以對于不同的動態(tài)場景，要綜合熱噪聲和動態(tài)應(yīng)力因素的影響，來設(shè)計環(huán)路中的預(yù)檢積分時間T、環(huán)路噪聲帶寬Bn。而在本文中，環(huán)路使用的為有源比例濾波器，具體的參數(shù)設(shè)置，可以參見文獻［10］。

2 高動態(tài)GPS載波跟蹤組合環(huán)路控制策略

2.1 高動態(tài)GPS載波跟蹤組合環(huán)路控制策略原理

本文提出的高動態(tài)GPS載波跟蹤組合環(huán)路控制策略原理（后簡稱載波跟蹤組合環(huán)路）如圖2所示。載波跟蹤組合環(huán)路主要由兩個模塊組成：一個是頻率牽引模塊，另一個為復(fù)合環(huán)路模塊。為了較好地適應(yīng)動態(tài)性并快速捕獲信號，捕獲算法中預(yù)檢積分時間取1ms，頻率井取值為666.67Hz，即捕獲后的最大多普勒頻差將達到333.33Hz，由于頻差較大，先由頻率牽引模塊將頻差牽引到較小的范圍內(nèi)，然后再利用復(fù)合環(huán)路實現(xiàn)對信號的精確跟蹤。

而復(fù)合環(huán)路模塊先利用I、Q 兩路值進行環(huán)路的跟蹤狀態(tài)判決，輸出控制信號對復(fù)合環(huán)路中鑒頻器和鑒相器輸出值進行自適應(yīng)加權(quán)，來調(diào)節(jié)復(fù)合環(huán)路中FLL和PLL的主導(dǎo)作用，以此形成一個動態(tài)載波跟蹤組合環(huán)路控制策略。

圖2 載波跟蹤組合環(huán)路框圖Fig.2 Block diagram of combination of carrier tracking loop

2.2 頻率牽引模塊

根據(jù)鑒頻器的鑒頻特性，本文選擇鑒頻范圍廣的四象限反正切鑒頻器作為頻率牽引模塊中的鑒頻器，頻率牽引模塊的原理圖如圖3所示。

圖3 頻率牽引模塊原理框圖Fig.3 Block diagram of frequency pulling module

由于接收機動態(tài)具有隨機性，估計值fe可視為一個變量，而頻率牽引模塊輸出的最終頻率誤差控制量fe＿out是對所有觀察量fe的最大似然估計，其計算公式［9］如式（10）：

其 中N 為 頻 率 牽 引 階 段 的 觀 測 值 次 數(shù)，femax和femin分別是N 次觀測值中的最大值和最小值，在頻率牽引模塊中，為了和捕獲環(huán)路相一致，將頻率牽引模塊中的預(yù)檢積分時間T 也設(shè)置為1ms，那么根據(jù)四象限反正切鑒頻器的工作特點，其鑒頻的范圍為±500Hz，可完全覆蓋捕獲后的頻率偏差。

而根據(jù)動態(tài)場景的不同也應(yīng)設(shè)置觀測次數(shù)N的值，其計算公式如式（11）［9］：

其中：Δf=500-333.33≈167Hz為四象限反正切鑒頻器的工作頻差冗余量。fL1為衛(wèi)星L1波段的載波頻率，a為接收機運動的加速度，C為光速。在本文試驗中的最大加速度為a=150 g，則t牽引≤27.4ms。由于導(dǎo)航數(shù)據(jù)位的限制，選取t牽引=20 ms，則N =20。

2.3 復(fù)合環(huán)路模塊

由于在高動態(tài)場景下，輸入信號與本地信號的頻差容易超出PLL 環(huán)的同步帶，造成失鎖顯現(xiàn)；而FLL環(huán)的同步帶和捕獲帶比PLL更寬，但是其精度不如PLL環(huán)，所以利用FLL 和PLL 相結(jié)合的組合環(huán)路是一種較好的解決跟蹤環(huán)路對動態(tài)適應(yīng)性、跟蹤精度要求的方法，所以本文復(fù)合環(huán)路采用2 階FLL輔助3階PLL的結(jié)構(gòu)，復(fù)合環(huán)路模塊的原理框圖如圖4所示。

圖4 復(fù)合環(huán)路模塊原理框圖Fig.4 Block diagram of compound loop module

相比傳統(tǒng)復(fù)合環(huán)路，文中復(fù)合環(huán)路先利用判決算法對跟蹤狀態(tài)進行粗略的預(yù)測（門限值判決），然后利用自適應(yīng)加權(quán)算法對FLL和PLL進行調(diào)節(jié)。在接收機具有較大動態(tài)性的時候，可以通過降低PLL、增加FLL在復(fù)合環(huán)路中的作用，使得FLL 能快速將輸入信號和本地信號之間的頻差牽引到PLL的快捕帶內(nèi)，加快復(fù)合環(huán)路的收斂速度；而接收機在動態(tài)性較弱的時候，可以通過增加PLL、降低FLL在復(fù)合環(huán)路中的作用，來提高環(huán)路的跟蹤精度。復(fù)合環(huán)路模塊中的算法具體如下：

1）門限值判決算法

門限值判決公式如下：

式中Ip（k）、Qp（k）的表達式如式（1）、（2）所示。在環(huán)路基本鎖定后，殘余多普勒Δfd趨近于零，此時碼 相 位 誤 差 也 基 本 為 零， 那 么R（ε（k））sinc（Δfd（k）πT）的值趨近于1，在忽略噪聲的影響后，判決式子簡化為：

其中θk為殘余載波相位，如果環(huán)路跟蹤穩(wěn)定后，則θk趨近于零，那么判決式子的值趨近于1。對一般高動態(tài)條件下，在忽略弱信號影響的情況下，該算法能粗略地反映出跟蹤環(huán)路的跟蹤狀態(tài)，從而為后面的自適應(yīng)權(quán)值算法提供加權(quán)控制信息。

2）自適應(yīng)加權(quán)算法

①權(quán)值因子的選取

在復(fù)合環(huán)模塊中，設(shè)FLL 中鑒頻器輸出值為Δf-d，由分析可知，該值能時刻地反應(yīng)跟蹤環(huán)的跟蹤狀態(tài)，如果輸入信號剛結(jié)束頻率牽引而進入復(fù)合環(huán)路時，輸入信號和本地信號還會具有較大的頻差，此時估計值Δf-d會偏大；當(dāng)環(huán)路處于鎖定狀態(tài)時，載體突然遭到較大的動態(tài)性沖擊，此時輸入信號具有較大的多普勒及多普勒頻率變化率，也會使得估計值Δf-d偏大。而當(dāng)環(huán)路趨于鎖定狀態(tài)或者環(huán)路已經(jīng)鎖定，但載體的運動緩和時，此時估計值Δf-d會較小。所以為了對FLL和PLL在不同狀態(tài)下的作用進行自適應(yīng)調(diào)節(jié)，特選取α、β兩對權(quán)值系數(shù)如下：

在本文中，接收機的最大加速度a=150 g、a′=±50 g／s，當(dāng)T 取值為1 ms時，（Hz）。由于sin（x）函數(shù)在這個區(qū)間上為遞增函數(shù)，當(dāng)Δf-d增加時，α1會隨著增加，而β1 會相應(yīng)的減小，式中T 為環(huán)路的預(yù)檢積分時間。

②當(dāng)判決值Vjudgment小于門限值H 時，說明輸入信號與本地信號頻差還較大，此時復(fù)合環(huán)需以FLL為主、PLL為輔。對FLL和PLL輸出值進行加權(quán)：

如果當(dāng)判決值Vjudgment大于門限值H 時，說明環(huán)路的輸入信號與本地信號的頻差較小，此時主要以PLL工作為主、FLL為輔。對FLL和PLL輸出值進行加權(quán)：

③將自適應(yīng)加權(quán)后的ΔF、ΔP 通過濾波器濾波融合后來調(diào)整NCO產(chǎn)生新的本地載波信號，從而保持環(huán)路對于輸入信號的跟蹤。

3 試驗驗證

為了驗證載波跟蹤組合環(huán)路控制策略的性能，本文在環(huán)路噪聲帶寬和預(yù)檢積分時間設(shè)置都相同的情況下，將傳統(tǒng)的2 階FLL 輔助3 階PLL（簡稱FLL／PLL）的復(fù)合環(huán)與載波跟蹤組合環(huán)路控制策略（簡稱F-PLL）在相同動態(tài)場景下對GPS信號進行了跟蹤試驗。

1）載體以最大加速度a=-150 g，a′=±50 g／s的加加速度運動的狀態(tài)圖如圖5所示。

在綜合考慮了熱噪聲誤差和動態(tài)應(yīng)力的情況下，根據(jù)式（7）、式（9），F(xiàn)LL／PLL、F-PLL 兩種復(fù)合環(huán)路中FLL的帶寬設(shè)置為10Hz，PLL的帶寬設(shè)置為18Hz，環(huán)路對載體運動狀態(tài)（1）跟蹤的相關(guān)結(jié)果圖如圖6-圖11所示。

從圖6、圖7中可以看出，由于F-PLL使用了自適應(yīng)加權(quán)技術(shù)，使得環(huán)路對于突發(fā)的加加速度對環(huán)路帶來的沖擊比FLL／PLL能更快地收斂。從圖8、圖9可以看出，F(xiàn)-PLL中鑒相器估計出的相位誤差均值比FLL／PLL小一個數(shù)量級，所以相比于FLL／PLL來說精度更高。

從圖10、11 可以看出，F(xiàn)LL／PLL 的I、Q 支路輸出值的離散圖比F-PLL的輸出值更加分散，前者的Q 路值集中于±600之間，而后者的Q 路值集中于±400之間，這說明前者的Q 路信號中除了含有噪聲外還含有大量的信號信息。因此F-PLL 相比于FLL／PLL來說具有更佳的跟蹤性能。

圖5 載體的運動狀態(tài)1Fig.5 Carrier state of motion 1

圖6 FLL／PLL中鑒相器輸出值Fig.6 Phase detector output value（without control strategies）

圖8 FLL／PLL中3階PLL環(huán)鑒相器輸出均值Fig.8 Output mean of 3-order PLL loop phase detector（without control strategies）

圖9 3階PLL環(huán)鑒相器輸出均值Fig.9 3-order PLL loop phase detector（with control strategies）

圖10 I、Q 兩路輸出值離散圖Fig.10 I，Q branch discrete figure（without control strategies）

圖11 I、Q 兩路輸出值離散圖Fig.11 I，Q branch discrete figure（with control strategies）

載體以a=150 g 勻加速運動的運動狀態(tài)圖如圖12所示。

根據(jù)式（7）、式（9），F(xiàn)LL／PLL、F-PLL兩個復(fù)合環(huán)路中FLL的帶寬設(shè)置為10 Hz，PLL 的帶寬設(shè)置為15Hz，環(huán)路對載體運動狀態(tài)（2）跟蹤的相關(guān)結(jié)果如圖13、圖14所示。從圖13、14可以看出，由于在初始階段，150 g／s的勻加速運動會短暫的帶來較大的沖擊，對于FLL／PLL來說未能跟蹤該信號，而F-PLL在大約200ms后進行收斂，能正確地解算出導(dǎo)航電文，說明F-PLL具有更佳的動態(tài)適應(yīng)性能。

圖12 載體的運動狀態(tài)2Fig.12 Carrier state of motion 2

圖13 FLL／PLL環(huán)跟蹤后的導(dǎo)航電文信息Fig.13 The output value of the navigation message（without control strategies）

圖14 F-PLL環(huán)路跟蹤后的導(dǎo)航電文信息Fig.14 The output value of the navigation message（with control strategies）

從上面的兩組實驗可以看出，基于頻率牽引和自適應(yīng)加權(quán)技術(shù)的高動態(tài)GPS載波跟蹤組合環(huán)路控制策略，能使得復(fù)合環(huán)路在遇到較大的動態(tài)沖擊下，能較快地收斂，具有更佳的跟蹤性能。

4 結(jié)論

本文在分析了FLL和PLL跟蹤技術(shù)的基本原理及特點的基礎(chǔ)上，提出了一種基于頻率牽引和自適應(yīng)加權(quán)技術(shù)的載波跟蹤組合環(huán)路控制策略。該策略先通過頻率牽引算法將捕獲后的頻差牽引到較小的范圍內(nèi)，然后利用帶有跟蹤狀態(tài)判決算法來自適應(yīng)調(diào)節(jié)FLL 和PLL 主導(dǎo)作用的復(fù)合環(huán)路對信號進行跟蹤。在對不同高動態(tài)運動場景下的GPS信號進行跟蹤仿真試驗后表明，該策略較好地解決了傳統(tǒng)復(fù)合環(huán)路存在的切換門限不易確定、環(huán)路切換頻繁、濾波不連續(xù)等問題，具有更好的跟蹤精度和對動態(tài)的適應(yīng)性，載體能在最大加速度a=150 g，最大加加速度a′=50 g／s的動態(tài)性指標(biāo)下正常工作，這為后期開發(fā)高性能GPS接收機奠定了基礎(chǔ)。

［1］王偉，張廷新，史平彥.高動態(tài)環(huán)境下GPS信號跟蹤算法綜述［J］.空間電子技術(shù)，2000（1）：1-9.

［2］Ward P W.Performance comparisons between FLL，PLL and a novel FLL-assisted-PLL carrier tracking loop under RF interference conditions［C］／／Proceedings of the 11th International Technical Meeting of the Satellite Division of The Institute of Navigation（ION GPS 1998）.US：ION，1998：783-795.

［3］Stefan K，Christian A，Daniel G，et al.GNSS receiver with vector based FLL-assisted PLL carrier tracking loop［C］／／ION GNSS 21st international Technical Meeting of the Satellite Division.Savannah CA：The institute of Navigation，2008：197-203.

［4］梁丹丹，張一，張中兆.高動態(tài)直擴接收機載波跟蹤技術(shù)研究［J］.電子技術(shù)應(yīng)用，2005（9）：51-53.

［5］田明坤，邵定蓉.高動態(tài)GPS接收機的一種設(shè)計方案［J］.遙測遙控，2002（3）：15-20.

［6］沈鋒，徐定杰，薛冰.基于FLL／PLL 相結(jié)合的載波跟蹤在導(dǎo)航接收機中的應(yīng)用［J］.中國航海，2004（2）：51-54.

［7］羅大成，王仕成，劉志國.一種基于FQFD／FLL／PLL的混合載波跟蹤算法［J］.航天控制，2009，27（1）：10-14.

［8］羅大成，王仕成，曾洪貴，等.一種GPS 軟件接收機的設(shè)計［J］.儀器儀表學(xué)報，2008，29（9）：1856-1861.

［9］胡輝，孫函子，紀(jì)兆云.高動態(tài)GPS載波跟蹤算法和環(huán)路控 制 策 略 研 究［J］.宇 航 學(xué) 報，2011，32（8）：1805-1812.

［10］Kaplan E D，Hegarty C.Understanding GPS：principles and applications［M］.London：Artech House Publishers，2005.