江濤,錢(qián)富才

(西安理工大學(xué) 自動(dòng)化與信息工程學(xué)院, 陜西 西安 710048)

基于UBB噪聲的高動(dòng)態(tài)GPS載波跟蹤算法研究

江濤,錢(qián)富才

(西安理工大學(xué) 自動(dòng)化與信息工程學(xué)院, 陜西 西安 710048)

GPS接收機(jī)在高動(dòng)態(tài)情況下,載波信號(hào)受到不同類(lèi)型復(fù)雜噪聲的干擾,且很難獲得這些噪聲的統(tǒng)計(jì)特性,而傳統(tǒng)方法對(duì)噪聲的處理都是假定其統(tǒng)計(jì)特性已知,因此,本文提出了一種基于集員濾波(SMF)的解決方案。該方法首先將噪聲定義成一類(lèi)未知但有界(UBB)噪聲,然后利用橢球包含UBB噪聲集合及狀態(tài)集合,進(jìn)而采用集員濾波的方法實(shí)現(xiàn)了載波信號(hào)的多普勒頻率估計(jì),最后通過(guò)Matlab仿真驗(yàn)證該算法的有效性,并與傳統(tǒng)方法進(jìn)行比較。較傳遞函數(shù)設(shè)計(jì)的PLL設(shè)計(jì)方法,提高了環(huán)路的動(dòng)態(tài)適應(yīng)性和精度;較Kalman濾波方法,提高了環(huán)路的魯棒性和快速性。集員濾波為GPS接收機(jī)載波跟蹤環(huán)路的設(shè)計(jì)提供了新的思路。

GPS; 集員濾波; UBB; 多普勒頻率; 載波跟蹤環(huán)路

GPS是由美國(guó)建立的高精度全球?qū)Ш蕉ㄎ幌到y(tǒng),已經(jīng)在陸地、海洋、航空、航天等領(lǐng)域得到廣泛應(yīng)用。一般GPS接收機(jī)在靜態(tài)、中低動(dòng)態(tài)的情況下能夠進(jìn)行準(zhǔn)確的定位及導(dǎo)航,然而在高動(dòng)態(tài)情況下,傳統(tǒng)GPS接收機(jī)無(wú)法準(zhǔn)確的定位導(dǎo)航。目前,國(guó)外對(duì)于高動(dòng)態(tài)GPS接收機(jī)產(chǎn)品及技術(shù)進(jìn)行封鎖,因此這項(xiàng)研究顯得至關(guān)重要。

傳統(tǒng)GPS載波跟蹤環(huán)路設(shè)計(jì)采用的是二階或三階鎖相環(huán)(Phase Locked Loop,簡(jiǎn)稱(chēng)PLL),其中環(huán)路濾波器的作用在于濾除載波環(huán)路中的噪聲,并且使得濾波結(jié)果能夠準(zhǔn)確地反應(yīng)出輸入信號(hào)的相位變化。然而當(dāng)接收機(jī)處在高動(dòng)態(tài)運(yùn)動(dòng)中,環(huán)路中的噪聲特性也隨之發(fā)生了變化。PLL為了提升對(duì)動(dòng)態(tài)的適應(yīng)能力,一般采用的方法是增加環(huán)路帶寬和階數(shù),然而這兩種方法都有各自的局限性。首先,為了適應(yīng)動(dòng)態(tài)應(yīng)力帶來(lái)的變化,增加環(huán)路噪聲帶寬雖然有助于環(huán)路的穩(wěn)定,但同時(shí)引入了更多的環(huán)路熱噪聲,對(duì)于環(huán)路噪聲帶寬的增加，動(dòng)態(tài)應(yīng)力與熱噪聲是一對(duì)矛盾,因此,總體載波的跟蹤性能并不能有效的提升;其次,高動(dòng)態(tài)意味著接收機(jī)接收到的信號(hào)具有較大的多普勒頻率,使得載波相位快速發(fā)生變化,傳統(tǒng)PLL動(dòng)態(tài)適應(yīng)的極限很難達(dá)到高動(dòng)態(tài)的要求。為了提高跟蹤環(huán)路的動(dòng)態(tài)性能,一般的解決方法是增加跟蹤環(huán)路的階數(shù)[1],但同時(shí)也帶來(lái)了設(shè)計(jì)復(fù)雜性和嚴(yán)重的穩(wěn)定性問(wèn)題[2]。因此,由傳遞函數(shù)設(shè)計(jì)的PLL無(wú)法應(yīng)對(duì)高動(dòng)態(tài)等極端條件。

為了解決高動(dòng)態(tài)的跟蹤問(wèn)題,文獻(xiàn)[3]～[7]采用了新的載波跟蹤環(huán)路PLL設(shè)計(jì)結(jié)構(gòu)，通過(guò)建立由載波相位差、多普勒頻率、多普勒頻率變化率組成的狀態(tài)方程,利用Kalman濾波器實(shí)現(xiàn)多普勒頻率估計(jì),取代了傳統(tǒng)的環(huán)路濾波器。由于Kalman濾波在噪聲統(tǒng)計(jì)特性為高斯白噪聲的情況下是一種最優(yōu)估計(jì),GPS接收機(jī)在靜止或者低動(dòng)態(tài)運(yùn)動(dòng)中,載波跟蹤環(huán)路主要受到熱噪聲的影響[8],熱噪聲是由元器件中電子熱振動(dòng)引起的,其統(tǒng)計(jì)特性就是典型的高斯白噪聲,因此,基于Kalman濾波器設(shè)計(jì)的載波跟蹤環(huán)路較傳統(tǒng)方法跟蹤性能有了大幅提升。然而Kalman濾波器在高動(dòng)態(tài)環(huán)境中的局限性也很明顯,首先,在高動(dòng)態(tài)環(huán)境中，動(dòng)態(tài)應(yīng)力成為主要的噪聲源[9],并且GPS接收機(jī)在極限的運(yùn)動(dòng)狀態(tài)下,振動(dòng)引起的時(shí)鐘相位抖動(dòng)、阿倫方差等噪聲已不可忽略[10],這些多源噪聲的疊加如果簡(jiǎn)單的用高斯白噪聲來(lái)描述其統(tǒng)計(jì)特性,顯然是不合理的,當(dāng)這些噪聲統(tǒng)計(jì)特性偏離高斯分布較大時(shí),Kalman濾波就不能得到符合實(shí)際的結(jié)果;其次,Kalman濾波過(guò)程中,過(guò)程噪聲和測(cè)量噪聲的統(tǒng)計(jì)特性是初始給定的,顯然在復(fù)雜的運(yùn)動(dòng)環(huán)境下,噪聲的統(tǒng)計(jì)特性是實(shí)時(shí)變化的,這可能導(dǎo)致濾波適應(yīng)能力變差甚至發(fā)散[11]。

雖然這些噪聲的統(tǒng)計(jì)特性很難通過(guò)數(shù)學(xué)方法得到,但是在特定的環(huán)境下,獲得噪聲邊界是可能的[12],此時(shí)可以采用基于有界集合的算法[13]。20世紀(jì)60年代末,Witsenhausen和Schweppe提出了一類(lèi)基于未知但有界(Unkown But Bounded簡(jiǎn)稱(chēng)UBB)噪聲假設(shè)的估計(jì)方法[14-15],即集員濾波。該方法利用外定界橢球集合包含系統(tǒng)的真實(shí)狀態(tài),只要求系統(tǒng)噪聲有界,且噪聲界已知,而不需要已知噪聲統(tǒng)計(jì)特性。但當(dāng)時(shí)沒(méi)有提出橢球集合的最優(yōu)化問(wèn)題,因而并沒(méi)有得到廣泛的應(yīng)用。直至Fogel和Huang提出了OBE算法[16],得到了最小跡橢球和最小容積橢球集合,以及Maksarov分析了狀態(tài)定界問(wèn)題[17],得到了橢球集合描述,此后學(xué)者們提出了許多相關(guān)的改進(jìn)算法,推動(dòng)了對(duì)于集員濾波的廣泛研究及應(yīng)用。在實(shí)際設(shè)計(jì)過(guò)程中,集員濾波過(guò)程類(lèi)似于Kalman濾波,通過(guò)建立系統(tǒng)狀態(tài)方程,運(yùn)用遞歸的方法進(jìn)行狀態(tài)估計(jì)。由于集員濾波對(duì)噪聲有較強(qiáng)的適應(yīng)性,該方法在高動(dòng)態(tài)環(huán)境下表現(xiàn)出較好的魯棒性和估計(jì)精度。因此,本文將引入集員濾波來(lái)改善高動(dòng)態(tài)GPS載波跟蹤環(huán)路的性能。

1 傳統(tǒng)PLL設(shè)計(jì)的GPS載波跟蹤環(huán)路

1.1 PLL設(shè)計(jì)原理

PLL由鑒相器、環(huán)路濾波器和數(shù)控振蕩器(NCO)組成,時(shí)域結(jié)構(gòu)如圖1所示,鑒相器將輸入信號(hào)相位θi與反饋信號(hào)相位θf(wàn)比較,結(jié)果乘以增益k0,得到相位差θe,經(jīng)過(guò)環(huán)路濾波器得到控制信號(hào)Vo,Vo控制載波NCO改變本地載波的輸出相位。在誤差信號(hào)的控制下,頻率就調(diào)整至與輸入信號(hào)頻率近似,最終使環(huán)路鎖定。

動(dòng)態(tài)環(huán)境的變化引起信號(hào)噪聲特性和多普勒頻率的變化,PLL階數(shù)越高,動(dòng)態(tài)適應(yīng)能力越強(qiáng)。二階PLL對(duì)加速度應(yīng)力敏感,三階PLL對(duì)加加速度應(yīng)力敏感。下面就噪聲對(duì)二階與三階PLL的影響進(jìn)行定量分析。

1.2 噪聲對(duì)動(dòng)態(tài)性能的影響

傳統(tǒng)PLL作為成熟的同步跟蹤技術(shù)已經(jīng)成功的應(yīng)用于GPS載波跟蹤環(huán)路。然而在高動(dòng)態(tài)條件下,環(huán)路中的動(dòng)態(tài)應(yīng)力變化劇烈,噪聲的變化對(duì)PLL有多大的影響呢?為了說(shuō)明問(wèn)題,首先介紹引起載波環(huán)路相位誤差的四種主要噪聲來(lái)源以及總體噪聲。

1) 熱噪聲σt

熱噪聲是PLL中最常見(jiàn)的噪聲源,其分布特性屬于典型的高斯白噪聲。引起鎖相環(huán)相位抖動(dòng)的熱噪聲計(jì)算如下:

(1)

2) 動(dòng)態(tài)應(yīng)力θd

動(dòng)態(tài)應(yīng)力的大小依賴(lài)于信號(hào)所處的動(dòng)態(tài)環(huán)境和帶寬,對(duì)于一個(gè)給定的環(huán)路,高動(dòng)態(tài)生成高動(dòng)態(tài)應(yīng)力,二階環(huán)路動(dòng)態(tài)應(yīng)力誤差定義為:

(2)

三階環(huán)路,動(dòng)態(tài)應(yīng)力誤差定義為:

(3)

式中,dR2/dt2為視線加速度(deg/s2),dR3/dt3為視線加加速度(deg/s3)。

3) 抖動(dòng)誤差συ

振動(dòng)引起的時(shí)鐘相位噪聲可計(jì)算如下:

(4)

4) 阿倫偏差θA

阿倫偏差同樣也能引起時(shí)鐘相位噪聲,對(duì)于二階環(huán)路,阿倫偏差引入的誤差為:

(5)

三階環(huán)路,阿倫偏差引入的誤差為:

(6)

式中,σA(τ)為晶振短期穩(wěn)定度閘門(mén)時(shí)間的方差。

5) 總體噪聲及跟蹤門(mén)限

根據(jù)一般經(jīng)驗(yàn),載波跟蹤環(huán)路門(mén)限定義如下:

(7)

式中,σPLL為均方根誤差為1σ的總體噪聲,該式說(shuō)明總體噪聲引起的1σ誤差必須小于15°,否則整個(gè)環(huán)路的跟蹤性能將不能保證。

文獻(xiàn)[1]給出了GPS載波跟蹤環(huán)路的噪聲及其1σ均方根誤差,本文在此基礎(chǔ)上分析傳統(tǒng)PLL可以承受的動(dòng)態(tài)極限,下面分析噪聲對(duì)動(dòng)態(tài)性能的影響。

忽略振蕩器振動(dòng)引起的抖動(dòng)誤差συ和阿倫偏差θA,式(7)可以表示為:

(8)

則動(dòng)態(tài)應(yīng)力誤差門(mén)限為:

(9)

由式(1)、(2)和(9),進(jìn)而可以確定二階PLL的加速度應(yīng)力門(mén)限:

(10)

由式(1)、(3)和(9),可得三階PLL的加加速度應(yīng)力門(mén)限:

(11)

式中,λL為射頻載波波長(zhǎng)。

如果假定載噪比C/N0=38 dBHz,環(huán)路噪聲帶寬Bn=18 Hz時(shí),T=1 ms。為了保證跟蹤門(mén)限,將參數(shù)代入式(1)、(11)計(jì)算可知,二階PLL必須滿(mǎn)足加速度小于21.56 m/s2,三階PLL必須滿(mǎn)足加加速度小于226.38 m/s3,因此,傳統(tǒng)PLL無(wú)法滿(mǎn)足典型高動(dòng)態(tài)的要求。

從另一視角來(lái)看,熱噪聲不會(huì)隨著動(dòng)態(tài)變化而變化;動(dòng)態(tài)應(yīng)力屬于模型階數(shù)的限制導(dǎo)致的建模誤差，高動(dòng)態(tài)情形下會(huì)引起振蕩器、時(shí)鐘晶振等元器件參數(shù)的變化,也就是除了熱噪聲以外,其余噪聲都會(huì)隨著動(dòng)態(tài)環(huán)境的變化而變化,即噪聲的非高斯部分會(huì)隨著動(dòng)態(tài)的增加而增加?；谶@兩個(gè)方面,我們需要尋找一個(gè)能夠適應(yīng)高動(dòng)態(tài)且噪聲為非高斯的估計(jì)方法來(lái)處理該問(wèn)題。

2 集員魯棒濾波設(shè)計(jì)方法

2.1 狀態(tài)空間模型

傳統(tǒng)PLL是基于傳遞函數(shù)設(shè)計(jì)的濾波器,其本質(zhì)是一種對(duì)系統(tǒng)內(nèi)部結(jié)構(gòu)進(jìn)行不完全描述的頻率分析方法。高動(dòng)態(tài)意味著接收機(jī)存在很大的加速度和加加速度,通過(guò)之前的結(jié)論可知,傳統(tǒng)二階、三階PLL并不能滿(mǎn)足要求,需要高階傳遞函數(shù)設(shè)計(jì)的系統(tǒng)才能降低動(dòng)態(tài)應(yīng)力誤差,然而高階傳遞函數(shù)的設(shè)計(jì)與分析是極其困難的。

集員濾波是一種時(shí)域?yàn)V波方法,它是基于狀態(tài)空間描述的估計(jì)理論,狀態(tài)空間模型能夠描述系統(tǒng)內(nèi)部各個(gè)狀態(tài)以及外部輸入輸出的數(shù)學(xué)關(guān)系,易給出高階系統(tǒng)模型,這就大大提升了系統(tǒng)動(dòng)態(tài)性能。根據(jù)載波相位差、多普勒頻率、多普勒頻率變化率、多普勒頻率變化率的變化率之間的數(shù)學(xué)關(guān)系,建立標(biāo)準(zhǔn)四階狀態(tài)空間[18]。

狀態(tài)方程:

(12)

測(cè)量方程:

(13)

式中,zk是積分時(shí)間T內(nèi),接收載波和本地產(chǎn)生載波的平均相位差,實(shí)際的觀測(cè)量由環(huán)路鑒相器輸出,υk為測(cè)量噪聲序列。

2.2 最優(yōu)濾波過(guò)程與噪聲分析

集員濾波和Kalman濾波相比,首先它們都是基于狀態(tài)空間模型的時(shí)域?yàn)V波方法,其次有著相似的遞歸過(guò)程。在噪聲統(tǒng)計(jì)特性為高斯白噪聲條件下,Kalman濾波是一種最優(yōu)估計(jì)。

Kalman濾波過(guò)程參考文獻(xiàn)[5],本地載波VCO產(chǎn)生的頻率輸入fk更新表示為:

(14)

式中,參數(shù)ka為可調(diào)參數(shù),用于調(diào)節(jié)更新過(guò)程中的振蕩幅度。

Kalman濾波要求ωk和υk的統(tǒng)計(jì)特性為高斯白噪聲且相互獨(dú)立,即:

(15)

(16)

(17)

其中:

設(shè)計(jì)Kalman濾波器過(guò)程中,都是假定系統(tǒng)模型精確已知,并且噪聲為高斯白噪聲,且方差陣Qk、Rk的選取來(lái)自于已知經(jīng)驗(yàn)值或試驗(yàn)測(cè)定值。然而實(shí)際GPS載波跟蹤環(huán)路存在動(dòng)態(tài)模型描述不精確、模型變化及非高斯噪聲干擾,均會(huì)導(dǎo)致Kalman濾波的性能下降甚至發(fā)散。由于真實(shí)系統(tǒng)的不確定性,大量實(shí)際應(yīng)用表明,Kalman濾波并不能保證濾波算法的收斂性和穩(wěn)定性。

2.3 UBB噪聲下的集員濾波

根據(jù)前面的分析可知,傳統(tǒng)PLL存在動(dòng)態(tài)適應(yīng)極限的問(wèn)題,而Kalman濾波存在對(duì)系統(tǒng)模型和噪聲統(tǒng)計(jì)特性過(guò)于依賴(lài)而導(dǎo)致魯棒性差的缺陷,為了解決這個(gè)問(wèn)題,通過(guò)1.2可以看出,除了熱噪聲以外,其余噪聲屬于非高斯,且隨著動(dòng)態(tài)性能的增加，噪聲非高斯部分也隨之增加,因此,UBB噪聲的假設(shè)更符合實(shí)際情況。

另一方面,噪聲引起的相位抖動(dòng)3σ值不超過(guò)π/4,也就是說(shuō),UBB噪聲引起的相位抖動(dòng)設(shè)定的上限就是π/4,因此,需要一種估計(jì)算法來(lái)保證估計(jì)值的邊界。幸運(yùn)的是,集員濾波就是一種基于UBB噪聲的保證估計(jì)算法。

集員濾波的基本思想就是找到由所有與測(cè)量數(shù)據(jù)、模型結(jié)構(gòu)和噪聲的有界假設(shè)相容的狀態(tài)組成的集合,由于此集合中的任何元素都有可能是狀態(tài)真實(shí)值,所以此集合被稱(chēng)為狀態(tài)的可行集。精確描述可行集是非常困難的,所以本文采用集員濾波算法通過(guò)橢球來(lái)包含可行集,從而實(shí)現(xiàn)對(duì)可行集的近似描述。集員濾波在時(shí)間更新和測(cè)量更新過(guò)程中,首先得到包含兩個(gè)橢球向量和的外定界橢球,這時(shí)的外定界橢球是一個(gè)橢球族,接著計(jì)算出該橢球族中的最小容積橢球。

橢球中心為a,P表示橢球形狀大小的正定矩陣,則橢球可表示為如下形式:

(18)

包含已知橢球E(a1,P1)和E(a2,P2)向量和的外定界橢球E(a,P)可表示為:

(19)

橢球E(a,P)的中心:

a=a1+a2

(20)

外定界橢球形狀完全由矩陣P確定,當(dāng)前已知矩陣P1和P2,因此,引入?yún)?shù)p,并且令

(21)

這樣,不同的參數(shù)p可以對(duì)應(yīng)不同的矩陣P,不同的P又可以產(chǎn)生不同的外定界橢球,在所有這些外定界橢球中,尋找容積最小的橢球。

橢球E(a,P)的容積可用P的行列式來(lái)度量,因此,求容積最小的外定界橢球等價(jià)于優(yōu)化問(wèn)題:

(22)

用與文獻(xiàn)[16]相同的方法,可以證明優(yōu)化問(wèn)題(22)等價(jià)于求如下方程:

(23)

由于噪聲以及動(dòng)態(tài)方程的不確定性,假定模型的狀態(tài)值、測(cè)量值、噪聲都是有界的并包含在橢球內(nèi),本文通過(guò)最小容積橢球方法來(lái)更新包含兩個(gè)橢球向量和的外定界橢球,更新橢球的中心,即為狀態(tài)更新值?；谠撍枷?下面給出GPS載波跟蹤環(huán)路模型的集員濾波過(guò)程:

假定ωk∈Rn,υk∈Rm分別為統(tǒng)計(jì)特性未知但有界的過(guò)程噪聲向量和測(cè)量噪聲向量。有界噪聲ωk和υk的橢球集合可以描述為:

(24)

(25)

式中,Qk∈Rn×n及Rk∈Rm×m為先驗(yàn)已知正定矩陣,先驗(yàn)值使得噪聲序列都包含在Qk、Rk設(shè)定的橢球內(nèi)。

假定已知系統(tǒng)初始狀態(tài)x0屬于橢球集合

(26)

k-1時(shí)刻狀態(tài)xk-1的橢球

(27)

根據(jù)最優(yōu)橢球定界算法[19],給出集員濾波過(guò)程。

時(shí)間更新:

(28)

矩陣Pk,k-1為:

(29)

為了使橢球容積達(dá)到最小,參數(shù)pk由以下方程解得,這里要求pk>0。

(30)

測(cè)量更新:

(31)

矩陣Pk為:

(32)

此時(shí),若qk≥0,

(33)

(34)

(35)

為了使橢球Ek的容積達(dá)到最小,參數(shù)qk通過(guò)下式求解:

(36)

(37)

(38)

(39)

3 仿真驗(yàn)證

為了測(cè)試不同方法的跟蹤性能,高動(dòng)態(tài)參照美國(guó)噴氣動(dòng)力實(shí)驗(yàn)室(JPL)的定義,設(shè)GPS接收機(jī)載體在兩個(gè)獨(dú)立的0.5 s內(nèi)分別有±980 m/s3的加加速度,期間被持續(xù)恒定的?245 m/s2加速度所分割,速度的初始值為100 m/s,加速度的初始值設(shè)定為-245 m/s2,整個(gè)運(yùn)動(dòng)過(guò)程8 s,如圖2所示。

GPS測(cè)試信號(hào)由GPS信號(hào)模擬器產(chǎn)生。選取19號(hào)衛(wèi)星,采樣頻率5 MHz,中頻(IF)信號(hào)1.25 MHz,載噪比C/N0=38 dBHz,積分時(shí)間T=1 ms。采集到的中頻信號(hào)通過(guò)Matlab仿真。

仿真首先采用傳統(tǒng)PLL跟蹤多普勒頻率,選取環(huán)路噪聲帶寬Bn=18 Hz,二階PLL仿真結(jié)果見(jiàn)圖3(a)。由于接收機(jī)載體初始加速度為-245 m/s2,絕對(duì)值遠(yuǎn)遠(yuǎn)大于-21.56 m/s2,二階PLL剛鎖定環(huán)路便很快失鎖,跟蹤失敗。三階PLL仿真結(jié)果如圖3(b)所示,三階PLL能夠跟蹤0～3 s的-245 m/s2加速度產(chǎn)生的多普勒頻率,但是3 s后出現(xiàn)了980 m/s3的加加速度,遠(yuǎn)遠(yuǎn)大于226.38 m/s3,因此3 s后三階PLL失鎖,跟蹤失敗。

接著分別采用集員濾波和Kalman濾波(KF)跟蹤多普勒頻率,仿真結(jié)果如圖4所示。集員濾波和Kalman濾波設(shè)計(jì)的環(huán)路都能夠跟蹤整個(gè)運(yùn)動(dòng)過(guò)程引起的多普勒頻率。為了測(cè)試狀態(tài)初值對(duì)兩種算法的影響,設(shè)定初始狀態(tài)多普勒頻率誤差為50 Hz。

圖5給出了集員濾波和Kalman濾波頻率估計(jì)的均方誤差曲線,由仿真結(jié)果可見(jiàn),集員濾波在8 s內(nèi)的誤差小于Kalman濾波,這是因?yàn)榧瘑T濾波具有較快的收斂速度,隨著時(shí)間推移,Kalman濾波精度會(huì)略好于集員濾波,也就是說(shuō)當(dāng)初始誤差較大或者由于外界因素導(dǎo)致?tīng)顟B(tài)值發(fā)生突變時(shí),集員濾波能以更快的收斂速度到達(dá)穩(wěn)態(tài),而Kalman濾波因收斂速度較慢,從而更易在狀態(tài)劇烈變化中導(dǎo)致濾波發(fā)散。

4 結(jié) 論

本文分析了載波信號(hào)的噪聲特性,發(fā)現(xiàn)傳統(tǒng)PLL實(shí)現(xiàn)的載波跟蹤環(huán)路不能適應(yīng)高動(dòng)態(tài)情形，而Kalman濾波雖然能夠適應(yīng)高動(dòng)態(tài),但存在對(duì)非高斯白噪聲適應(yīng)能力差的問(wèn)題。本文采用基于UBB噪聲假設(shè)的集員濾波方法能夠解決高動(dòng)態(tài)、噪聲非高斯情況下的載波跟蹤問(wèn)題,且仿真表明，該方法在保證精度的基礎(chǔ)上具有較強(qiáng)的魯棒性。值得注意的是，集員濾波算法需要設(shè)定噪聲邊界,邊界過(guò)大會(huì)導(dǎo)致濾波精度下降,過(guò)小會(huì)導(dǎo)致估計(jì)值不可靠,如何選取適當(dāng)?shù)倪吔缫蕴岣邽V波性能是未來(lái)研究的方向。

[1]Kaplan E D. Understanding GPS principles and applications [M]. Boston: Artech House, 1996.

[2]Hamm C R, Flenniken W S, Bevly D M. Comparative performance analysis of aided carrier tracking loop algorithms in high noise/high dynamic environments [C]∥Proceedings of International Technical Meeting of the Satellite Division of the Institute of Navigation, Long Beach, 2004: 523-533.

[3]Ziedan N I, Garrison J L. Extended kalman filter-based tracking of weak GPS signals under hign dynamic conditions [C]∥Proceedings of International Technical Meeting of the Satellite Division of the Institute of Navigation, Long Beach, 2004: 20-31.

[4]Psiaki M L. Smoother-Based GPS signal tracking in a software receiver [C]∥ION GNSS 2001, Lake City, 2001: 2900-2913.

[5]Driscoll C O, Petovello M, Lachapelle G. Choosing the coherent integration time for kalman filter-based carrier-phase tracking of GNSS signals [J]. GPS Solutions, 2011, 15(4): 345-356.

[6]Wang J, Liang Q H, Liang K, et al. A new extended kalman filter based carrier tracking loop [C]∥Microwave, Antenna, Propagation and EMC Technologies for Wireless Communications, Guilin, 2009: 1181 -1184.

[7]Yuan G N, Xie Y J, Song Y, et al. Multipath parameters estimation of weak GPS signal based on new colored noise unscented kalman filter [C]∥International Conference on Information and Automation, Harbin, 2010: 1852-1856.

[8]左啟耀, 袁洪, 林寶軍. 高動(dòng)態(tài)環(huán)境下GPS信號(hào)跟蹤環(huán)路優(yōu)化算法研究 [J]. 宇航學(xué)報(bào), 2008, 29(2): 550-555. Zuo Qiyao, Yuan Hong, Lin Baojun. Optimal algorithm research of GPS signal tracking loop under the high dynamic circumstance [J]. Journal of Astronautics, 2008, 29(2): 550-555.

[9]陳竹, 王元?dú)J, 侯孝民. 高動(dòng)態(tài)條件下測(cè)控系統(tǒng)載波跟蹤環(huán)優(yōu)化設(shè)計(jì) [J]. 電訊技術(shù), 2010, 50(8): 57-62. Chen Zhu, Wang Yuanqin, Hou Xiaomin. Optimization of TT and C carrier tracking loop in dynamic environments [J]. Telecommunication Engineering, 2010, 50(8): 57-62.

[10]金天, 王麗華, 黃炳勝. 基于Kalman濾波的GPS跟蹤環(huán)路晶振閃爍噪聲建模方法 [J]. 中國(guó)空間科學(xué)技術(shù), 2012(5): 34-39. Jin Tian, Wang Lihua, Huang Bingsheng. Modeling method for oscillator flicker noise in Kalman-based GPS tracking loop [J]. Chinese Space Science and Technology, 2012(5): 34-39.

[11]Salychev S O. Inertial systems in navigation and geophysics [M]. Moscow: Bauman MSTU Press, 1998.

[12]Walter E, Piet L H. Estimation of parameter bounds from bounded-error data: a survey [J]. Mathematics and Computers Simulation, 1990, 32(3): 449-468.

[13]吳宏鍔, 王國(guó)欣, 牛玉俊. 有界噪聲作用下Duffing系統(tǒng)脈沖控制 [J]. 西安理工大學(xué)學(xué)報(bào), 2013, 29(2): 233-237. Wu Honge, Wang Guoxin, Niu Yujun. Chaos control of Duffing system with bounded noise excitation with periodic impulsive signals [J]. Journal of Xi’an University of Technology, 2013, 29(2): 233-237.

[14]Witsenhausen H S. Sets of possible states of linear systems given perturbed observations [J]. IEEE Transactions on Automatic Control, 1968, 13(5): 556 -558.

[15]Schweppe F C. Recursive state estimation: unknown but bounded errors and system inputs [J]. IEEE Transactions on Automatic Control, 1968, 13(1): 22-28.

[16]Fogel E, Huang Y F. On the value of information in system identification-bounded noise case [J]. Automatica, 1982, 18(2): 229-238.

[17]Maksarov D G, Norton J P. State bounding with ellipsoidal set description of the uncertainty [J]. International Journal of Control, 1996, 65(5): 847-866.

[18]Psiaki M L, Hee J. Extended Kalman filter methods for tracking weak GPS signals [C]∥Proceeding of ION GPS 2002, Portand, 2002: 2539-2553.

[19]柴偉, 孫先仿. 橢球狀態(tài)定界的魯棒算法 [J]. 北京航空航天大學(xué)學(xué)報(bào), 2006, 32(12): 1447-1450. Chai Wei, Sun Xianfang. Robust ellipsoidal state bounding algorithm[J]. Journal of Beijing University of Aeronautics and Astronautics, 2006, 32(12): 1447 -1450.

(責(zé)任編輯 周蓓)

Research on high dynamic GPS carrier tracking algorithm based on UBB noise

JIANG Tao,QIAN Fucai

(Faculty of Automation and Information Engineering, Xi’an University of Technology, Xi’an 710048, China)

In high dynamic condition, the carrier signal of GPS is interfered by different kinds of complex noises, and it is very difficult to obtain the statistic characteristics of the noise, but the traditional method is assumed that the statistic properties of the noise is known, accordingly, this paper proposes a solution based on set membership filter (SMF). Firstly, the noise is defined as an unknown but bounded (UBB) noise, and then the UBB noise set and the state set are contained by the ellipsoid, Doppler frequency estimation of carrier signal is accomplished by using the method of set membership filter next. Finally, the effectiveness of the proposed algorithm is verified by Matlab simulation and compared with the traditional method. Compared with the method based on transfer function of tracking loop, the proposed algorithm improves the dynamic adaptability and tracking precision; both robustness and speediness of the loop are improved compared with the performance of Kalman filter method. Accordingly, set membership filter is a new design method in the field of GPS receiver carrier tracking loop.

GPS; set membership filter; UBB; Doppler frequency; carrier tracking loop

1006-4710(2015)03-0328-07

2015-04-18

國(guó)家自然科學(xué)基金資助項(xiàng)目(61273127,61304204);國(guó)家自然科學(xué)基金重點(diǎn)項(xiàng)目(61533014);陜西省科技創(chuàng)新團(tuán)隊(duì)資助項(xiàng)目(2013KCT-04)。

江濤,男,博士生,研究方向?yàn)樾l(wèi)星導(dǎo)航、系統(tǒng)辨識(shí)與濾波算法。E-mail:jiangtao.xaut@foxmail.com。

錢(qián)富才,男,教授,博導(dǎo),博士,研究方向?yàn)殡S機(jī)控制、系統(tǒng)辨識(shí)、非線性系統(tǒng)等。E-mail:fcqian@xaut.edu.cn。

TN967.1

A