黎 劍，程 功，陳 占,聶俊堂

(昆明冶金高等專科學(xué)校測繪學(xué)院,云南 昆明 650033)

0 引 言

GNSS高程是基于WGS-84橢球面的大地高H，而實際測量工程則采用以似大地水準面為基準面的正常高Hr，兩者間的差距(即高程異常ξ)導(dǎo)致了GNSS高程在應(yīng)用上的問題。因此，如何高效、精確地實現(xiàn)從大地高H向正常高Hr的轉(zhuǎn)換成為許多測量工作者亟待解決的問題。

GNSS高程擬合模型大體分兩類：線狀擬合模型和面狀擬合模型。面狀擬合模型種類多，各有優(yōu)缺點，曾有同行進行過綜合模型以及分區(qū)擬合的實驗研究，均取得理想效果；線狀模型種類相對較少，最常用的是多項式曲線擬合模型和三次樣條曲線擬合模型,其原理是根據(jù)高程控制點的平面坐標X或Y及其高程異常值ξ，通過構(gòu)造一個插值函數(shù)來擬合線狀工程沿線方向上的似大地水準面曲線，然后據(jù)此內(nèi)插其它點的高程異常。

本文將從三次樣條函數(shù)及其插值的定義入手，根據(jù)其性質(zhì)和嚴密的數(shù)學(xué)推導(dǎo)過程構(gòu)造出分段高程異常函數(shù)的表達式，并將其運用到GNSS高程的線性擬合中。

1 分段高程異常函數(shù)的構(gòu)造

1.1 三次樣條函數(shù)的定義

若函數(shù)S(x)∈C2(a,b)， 且在每個小區(qū)間[xi,xi+1]上是三次多項式，其中a=x0

據(jù)此，可構(gòu)造出分段高程異常函數(shù)ξ(x)。 具體做法是在高程異常擬合中，將三次樣條函數(shù)S(x)用高程異常函數(shù)ξ(x)代替，并以xi(或yi)表示控制點的縱坐標(或橫坐標)，則可將三次樣條函數(shù)的思想應(yīng)用于線狀的高程異常擬合中。

1.2 構(gòu)造分段高程異常函數(shù)ξ(x)

因為ξ(x)在區(qū)間[xi,xi+1]上是三次多項式，故ξ″(x)在[xi,xi+1]上是線性函數(shù)，設(shè)其二階導(dǎo)數(shù)值為ξ″(xi)=Mi(i=0,1,…,n)，則

(1)

式(1)中，hi=xi+1-xi。

對ξ″(x)積分兩次，并根據(jù)三次樣條插值函數(shù)的定義定出積分常數(shù)，得到高程異常函數(shù)的表達式為：

(2)

1.3 確定未知量Mi

為確定未知量Mi，對ξ(x)求導(dǎo)得:

(3)

根據(jù)ξ(x)在[a,b]上二階導(dǎo)數(shù)連續(xù)，在節(jié)點處滿足連續(xù)性條件S′(xi-0)=S′(xi+0)，將(3)式代入該條件中得：

μiMi-1+2Mi+λiMi+1=di

(4)

(5)

要求出ξ(x)，需要在每個小區(qū)間上確定4個待定系數(shù)，通常可在區(qū)間的端點上各加一個條件，稱為邊界條件，分別用m1和mn表示。m1和mn為樣條曲線邊界的一階導(dǎo)數(shù)，即邊界的斜率，可以利用邊界的4個已知點求出邊界的斜率[4]。

常見邊界條件有3種：

3)當f(x)是以xn-x0為周期的函數(shù)時，則ξ(x)也是周期函數(shù)，此時：

ξ(x0+0)=ξ(x0-0)，ξ′(x0+0)=ξ′(x0-0)，ξ″(x0+0)=ξ″(x0-0)。

(6)

(7)

線性方程組(7)式是關(guān)于Mi的三對角線性方程組，此方程組中的系數(shù)矩陣元素λi，μi已完全確定，并且滿足λi≥0，μi≥0，λi+μi=1。因此系數(shù)矩陣為嚴格對角占優(yōu)矩陣，從而方程組(8)有唯一解[4]。

1.4 確定分段高程異常函數(shù)ξ(x)

式(7)求解方法為追趕法，將解得結(jié)果代入式(2)即可確定分段高程異常函數(shù)ξ(x)。

(8)

式(8)即為所構(gòu)造的分段高程異常函數(shù)與控制點的x坐標(或y坐標)相關(guān)的函數(shù)，這一函數(shù)為分段函數(shù)，對應(yīng)每個區(qū)間有相應(yīng)的表達式，代入?yún)^(qū)間內(nèi)的任一x坐標(或y坐標)，即可求出該點的高程異常。

2 實例計算

2.1 實例數(shù)據(jù)

表1為某新修高速公路的GNSS控制網(wǎng)數(shù)據(jù)，公路全長約 35 km，控制點全部實測了四等水準(由于篇幅原因，只給出部分數(shù)據(jù))。

表1 控制點數(shù)據(jù)Tab.1 Data of the control point

2.2 擬合方案

為驗證所述方法的可行性，分別利用5個已知點和10個已知點2種方法進行擬合計算，并運用Matlab軟件按上節(jié)所述方法分別進行擬合計算，擬合曲線見圖1，擬合殘差曲線見圖2。

圖1 Matlab擬合曲線圖 圖2 2種方案精度對比圖Fig.1 Fitting curve of Matlab Fig.2 Precision comparison of the two schemes

圖1中紅色點表示已知點，其余點表示擬合點。將擬合結(jié)果與已知高程異常值進行比較，并求出差值及中誤差，結(jié)果見表2。由圖1可以看出，擬合曲線實際上是由分段的三次多項式曲線拼接而成的連續(xù)曲線，在連接點處，不僅函數(shù)自身連續(xù)，其一階導(dǎo)數(shù)和二階導(dǎo)數(shù)也是連續(xù)的。

2.3 擬合精度

1)由圖2可知,將2種方法對相同點擬合的精度進行對比，用10個已知點擬合的殘差曲線相較于5個已知點來說相對平穩(wěn)，證明穩(wěn)定性較高；2種方案體現(xiàn)出相同的變化趨勢，證明方法可靠。

2)分別計算用5個已知點和10個已知點2種方法擬合的內(nèi)符合精度和外符合精度，結(jié)果由表2、3可知，增加擬合點數(shù)量后，內(nèi)符合精度變化不明顯，說明擬合點數(shù)量增加，可以提高擬合精度。

3)由表3可知，利用構(gòu)造的分段函數(shù)擬合高程異常值與已知的高程異常值差值較小，最大殘差不超過 3 mm，在工程應(yīng)用中，可替代四等水準測量。

表2 不同已知點的擬合精度對比Tab.2 Comparison of fitting accuracy of different known points mm

表3 高程異常值擬合結(jié)果Tab.3 The height anomaly value fitting results mm

3 結(jié) 語

本文引入分段函數(shù)的思想，構(gòu)造分段高程異常擬合函數(shù)，經(jīng)過實驗計算可以看出：利用分段高程異常函數(shù)擬合線狀GNSS高程，確定函數(shù)表達式后只需將坐標值代入對應(yīng)的分段函數(shù)表達式內(nèi)，即可求任一點的高程異常，過程準確嚴密，實驗結(jié)果準確可靠，證明了利用分段函數(shù)擬合線狀GNSS高程的可行性。

此外，筆者同時也用Y坐標進行了擬合實驗，實驗結(jié)果證明:利用X進行擬合的精度高于利用Y的擬合精度，原因有待于進一步研究。為得到良好的擬合結(jié)果，擬合前有必要對異常數(shù)據(jù)進行判定和剔除。