整周未知數(shù)的求解方法

整周未知數(shù)的求解方法孫逸倫09資環(huán)20091303015摘要：整周未知數(shù)確定后，測相偽距與測碼偽距的觀測方程在形式上將一致，此時只要同步觀測的衛(wèi)星數(shù)不少于4，即使觀測一個歷元，也可獲得唯一定位結(jié)果。因此，在載波相位觀測中，如果能預先消去或者快速地解算整周未知數(shù)，將大大縮短必要的觀測時間如果整周未知數(shù)作為待定量，與其它未知參數(shù)一起在數(shù)據(jù)處理中一并求解，則根據(jù)情況，將需要長達1-3小時的觀測時間。因為在同步觀測4顆衛(wèi)星的情況下，為解算整周未知數(shù)，理論上至少觀測3個歷元。但如果同步觀測時間很短，所測衛(wèi)星的幾何分布變化很小，使站星距離變化也很小，將降低不同歷元觀測結(jié)果的作用，在平差計算中，法方程的性質(zhì)將變壞，影響解的可靠性。準確快速地解算整周未知數(shù)，無論對保障相對定位精度，還是開拓高精度動態(tài)定位應用領(lǐng)域，都有重要意義。關(guān)鍵詞：整周未知數(shù)載波相位靜態(tài)相對定位法交換接收天線法搜索法0引言整周未知數(shù)解算方法分類按解算時間長短劃分：經(jīng)典靜態(tài)相對定位法和快速解算法。經(jīng)典靜態(tài)相對定位法：將其作為待定量，在平差計算中求解，為提高解的可靠性，所需觀測時間較長?？焖俳馑惴òǎ航粨Q天線法、P碼雙頻技術(shù)、濾波法、搜索法和模糊函數(shù)法等，所需觀測時間較短，一般為數(shù)分鐘。按接收機狀態(tài)區(qū)分：靜態(tài)法和動態(tài)法。前述的快速算法，雖然觀測時間很短，仍屬靜態(tài)法，動態(tài)法是在接收機載體的運動過程中確定整周未知數(shù)的方法。1確定整周未知數(shù)的經(jīng)典靜態(tài)相對定位法該方法在長距離靜態(tài)相對定位中是一種常用方法，其數(shù)學模型有單差和雙差模型。也可采用三差模型，首先消除整周未知數(shù)，在觀測站坐標確定后，再根據(jù)單差和雙差模型，求解相應的整周未知數(shù)。在平差計算中，整周未知數(shù)的取值分兩種情況：整數(shù)解（固定解）：將平差計算所得的整周未知數(shù)取為相近的整數(shù)，并作為已知數(shù)代入原方程，重新解算其它待定參數(shù)。當觀測誤差和外界誤差（或殘差）對觀測值影響較小時，該方法較有效，一般應用于基線較短的相對定位中。非整數(shù)解（實數(shù)解或浮動解）：如果外界誤差影響較大，求解的整周未知數(shù)精度較低（誤差影響大于半個波長），將其湊成整數(shù)，無助于提高解的精度。此時，不考慮整周未知數(shù)的整數(shù)性質(zhì)，平差計算所得的整周未知數(shù)，不再進行湊整和重新計算。一般用于基線較長相對定位中。2交換接收天線法原理：在觀測之前，先在基準站附近5-10m處選擇一個天線交換點，將兩臺接收機天線分別安置在該基線兩端，同步觀測2-8個歷元后，相互交換天線，并繼續(xù)觀測若干歷元；

最后將兩天線恢復到原來位置。此時固定站與天線交換點之間的基線向量視為起始基線向量，利用天線交換前后的同步觀測量，求解基線向量，進而確定整周未知數(shù)。假設(shè)在固定站1和天線交換點2的接收機，于歷元t1同步觀測了衛(wèi)星j、k,在忽略大氣折射影響的情況下，可得單差觀測方程：j(t)-Pj(tJ-ANj+fAt(t)21111k(k(t)-Pk(t)LANk+fAt(t)21111A9k(t)=-1九相應的雙差觀測方程為1VA9(t)=-1九上式中k(t)-Pk(t)-Pj(t)+Pj(t1VA9(t)=-1九上式中21112111ANj=Nj(t)-Nj(t)20 1 0ANk=Nk(t)-Nk(t)20 1 0當兩接收機交換天線后，于歷元t2同步觀測相同衛(wèi)星j、k,則單差觀測方程為:j(t)-Pj(tJ+ANj+fAt(t)221221A^k(t)=-2九相應的雙差觀測方程為1VA^(t)=-2九1A^k(t)=-2九相應的雙差觀測方程為1VA^(t)=-2九22122k(t)-Pk(t)-Pj(t)+Pj(t)】+ANk—ANj22122212取相應歷元t1、t2的雙差之和，則有=Lpk(t)-Apj(t)+Apk(t)-Apj(t)]九 2211其中ApApk(t)=2Apj(t)=2Apk(t)=1Apj(t)=1Pk(t)-Pk(t)2212Pj(t)-Pj(t)2212Pk(t)-Pk(t)21 11Pj(t)-Pj(t)21 11上述模型與靜態(tài)三差模型相類似，區(qū)別在于上式是根據(jù)不同歷元同步觀測量的雙差之和而建立的。由于所選起始基線很短，此時衛(wèi)星軌道誤差和大氣折射誤差對該模型的影響可忽略不計。上式的求解條件與雙差相同。根據(jù)上式確定起始基線向量后，可根據(jù)雙差模型確定整周未知數(shù)。該方法觀測時間短(數(shù)分鐘)，精度較高，操作方便，在準動態(tài)相對定位中得到應用。3確定整周未知數(shù)的搜索法1990年E.Frei和GBeutler提出了一種快速解算整周未知數(shù)的方法(fastambiguityresolutionapproachFARA)?；舅枷胧牵阂詳?shù)理統(tǒng)計理論的參數(shù)估計和假設(shè)檢驗為基礎(chǔ)，利用初始平差的解向量(點的坐標和整周未知數(shù)的實數(shù)解)及其精度信息(方差與協(xié)方差和單位權(quán)中誤差)，確定在某一置信區(qū)間整周未知數(shù)可能的整數(shù)解的組合，然后將整周未知數(shù)的每一組合作為已知值，重復進行平差計算，其中使估值的驗后方差(或方差和)為最小的

一組整周未知數(shù)就是所搜索的整周未知數(shù)的最佳估值?，F(xiàn)以載波相位觀測值雙差模型為例：假設(shè)在基線兩端對同一組衛(wèi)星(衛(wèi)星數(shù)為nj)進行同步觀測，觀測歷元數(shù)為nt,相應的誤差方程組已知為「8X2N8X2N=8Z8X2N=8Z》21N.」nJ-1LN12經(jīng)過初始平差后，相應整周未知數(shù)解向量的協(xié)因數(shù)陣為QNN,單位權(quán)驗后方差估算式:VtPVm2=0 (n-u)其中n為觀測方程數(shù)，u為未知量個數(shù)，n-u為自由度。則任一整周未知數(shù)經(jīng)初始平差后實數(shù)解的中誤差為m=mQNj 0NiNi在一定置信水平條件下，相應任一整周未知數(shù)的置信區(qū)間為N—mt(a/2)<N<N+mt(a/2)i N i i Niii=l,2,...,nj-l其中t(a/2)為顯著水平a和自由度的函數(shù)。當a和自由度確定后，t(a/2)值可由t值分布表中查得。例如：當取a=0.001，n-u=40時，得t(a/2)=3.55。如果初始平差后得Ni=9.05,mNi=0.78,則Ni的置信區(qū)間為6.28<Ni<11.8。其置信水平為99.9%,在上述區(qū)間整數(shù)Ni的可能取值為6、7、8、9、10、11、12。設(shè)Ci為Ni的可能取值數(shù)，由向量N=(N1,N2,…,Nnj-1),可得整數(shù)組合的總數(shù)nJ—1C=口C=C-C...Ci 1 2 nJ—1i=1如果觀測的衛(wèi)星數(shù)為nj=6,而每個整周未知數(shù)在其置信區(qū)間內(nèi)均有7個可能的整數(shù)取值，按上式可能的組合數(shù)為75=16807，對雙頻接收機則為33614。將上述整周未知數(shù)的各種可能組合，依次作為固定值，代入相應的誤差方程組中，進行平差計算，最終取坐標值的驗后方差為最小的一組平差結(jié)果，作為整周未知數(shù)的最后取值。參考文獻：邱蕾、花向紅等，GPS短基線