魏 榮 寶

(山東省地質(zhì)礦產(chǎn)勘查開發(fā)局第五地質(zhì)大隊(duì)，山東 泰安 271000)

0 引言

在當(dāng)今測(cè)量工作中，GPS越來越發(fā)揮著重要的作用，尤其是在布設(shè)工程控制網(wǎng)方面。利用GPS進(jìn)行工程控制網(wǎng)的布設(shè)與解算，越來越受到測(cè)繪工作者的喜愛。在用GPS進(jìn)行了外業(yè)測(cè)量后，通常用隨機(jī)所帶的解算軟件進(jìn)行解算。在本文中，作者利用條件平差，進(jìn)行GPS測(cè)量數(shù)據(jù)的解算。

1 平差問題

在某次控制測(cè)量中，一組同學(xué)布設(shè)了如圖1所示的GPS網(wǎng)，點(diǎn)1為已知點(diǎn)，點(diǎn)2、點(diǎn)3、點(diǎn)4為待定坐標(biāo)點(diǎn)，現(xiàn)用GPS接收機(jī)觀測(cè)了5條邊的基線向量(ΔXij,ΔYij,ΔZij)。

已知1點(diǎn)的坐標(biāo)為：

X1=-1 054 581.276 1 m，Y1=5 706 987.139 7 m，Z1=2 638 873.815 2 m。

基線向量觀測(cè)值及其協(xié)因數(shù)見表1。

表1 基線向量觀測(cè)值及其協(xié)因數(shù)

設(shè)各基線向量互相獨(dú)立，試用條件平差法求：基線向量改正數(shù)及其平差值。

2 問題解算

2.1 多余觀測(cè)數(shù)的確定

圖1中，測(cè)得5條基線向量，則n=15，t=9，r=6。

2.2 條件方程

1)平差值形式的條件方程。

根據(jù)幾何條件得平差值條件方程：

2)改正數(shù)形式的條件方程。

設(shè)改正數(shù)向量為V=[VX1VY1VZ1VX2VY2VZ2VX3VY3VZ3VX4VY4VZ4VX5VY5VZ5]T，改正數(shù)條件方程形式AV+W=0，則改正數(shù)形式的條件方程為:

-VX1+VX2-VΔX3+(-ΔX1+ΔX2-ΔX3)=0；

-VY1+VY2-VY3+(-ΔY1+ΔY2-ΔY3)=0；

-VZ1+VZ2-VZ3+(-ΔZ1+ΔZ2-ΔZ3)=0；

VX3-VX4+VX5+(ΔX3-ΔX4+ΔX5)=0；

VY3-VY4+VY5+(ΔY3-ΔY4+ΔY5)=0；

VZ3-VZ4+VZ5+(ΔZ3-ΔZ4+ΔZ5)=0。

代入已知數(shù)據(jù)及觀測(cè)值，則得相關(guān)矩陣。

3)條件方程系數(shù)陣A和常數(shù)項(xiàng)矩陣W見表2。

表2 條件方程系數(shù)陣A和常數(shù)項(xiàng)矩陣W

4)權(quán)陣的確定。

依據(jù)平差問題中給定的協(xié)因數(shù)陣，確定權(quán)陣，該權(quán)陣為一個(gè)15行15列的方陣，也是一個(gè)近似對(duì)角的方陣。

2.3 觀測(cè)值的改正數(shù)和平差值

1)改正數(shù)(單位：mm)。

V=[-6.3 -1.9 1.9 6.2 1.8 -1.9 2.6 -7.1 -3.1 -10.3 3.0 4.8 10.2 -3.2 -4.8]T。

2)平差值見表3。

表3 平差值

3 結(jié)語

綜合以上計(jì)算，可以看出，對(duì)于簡單的GPS網(wǎng)，也可以采用條件平差方法進(jìn)行解算，其整個(gè)過程操作間大，容易實(shí)現(xiàn)。

[1] 武漢大學(xué)測(cè)繪學(xué)院測(cè)量平差學(xué)科組.誤差理論與測(cè)量平差基礎(chǔ)[M].第2版.武漢:武漢大學(xué)出版社，2009.

[2] 苗元欣.基于一元線性回歸的變形監(jiān)測(cè)數(shù)據(jù)處理與分析[J].山西建筑，2013，39(33)：206-207.

[3] 泥立麗，王 永.基于Excel的繪制誤差曲線的方法[J].礦山測(cè)量，2010(3)：20-23.

[4] 高騰飛.L波段D-InSAR在礦區(qū)地面沉降監(jiān)測(cè)中的應(yīng)用[J].中國科技論文在線，2016(15)：178-179.