史建青，董春來(lái)

(淮海工學(xué)院，江蘇連云港222005)

一、平差原理及數(shù)學(xué)模型

導(dǎo)線網(wǎng)平差采用適宜電算的簡(jiǎn)接平差法，導(dǎo)線網(wǎng)誤差方程包括角度觀測(cè)誤差方程和邊長(zhǎng)觀測(cè)誤差方程。

間接平差基礎(chǔ)方程及其解為

按函數(shù)極值的求法，極值函數(shù)為

求其一階偏導(dǎo)數(shù)，并令其為0，得

代入誤差方程得

即為法方程式

1．邊長(zhǎng)觀測(cè)誤差方程

設(shè)導(dǎo)線邊jk 的兩端點(diǎn)坐標(biāo)(Xj，Yj)、(Xk，Yk)為未知參數(shù)。則 jk 的距離可表示為其為非線性函數(shù)，應(yīng)先進(jìn)行線性化。

于取初始近似值 X°j、Y°j、X°k、Y°k處進(jìn)行 Taylor級(jí)數(shù)展開(kāi)，略去二次以及二次以上項(xiàng)，整理得

其中

當(dāng)j點(diǎn)已知時(shí)

當(dāng)k點(diǎn)已知時(shí)

2．角度觀測(cè)誤差方程

當(dāng)一個(gè)測(cè)站上有兩個(gè)以上方向觀測(cè)時(shí)，由方向觀測(cè)值求得的角值之間是相關(guān)的。如圖1所示。

圖1 角度觀測(cè)示例圖

設(shè)j、k、h均為待定點(diǎn)，Li為測(cè)角網(wǎng)中的一個(gè)觀測(cè)值，得

以坐標(biāo)方位角改正數(shù)方程代入誤差方程，則角度的誤差方程為

二、導(dǎo)線網(wǎng)平差程序設(shè)計(jì)流程及關(guān)鍵程序

1．程序設(shè)計(jì)流程

程序設(shè)計(jì)具體流程示意圖如圖2所示。

圖2 導(dǎo)線網(wǎng)平差程序設(shè)計(jì)流程圖

2．關(guān)鍵程序示例

導(dǎo)線網(wǎng)程序坐標(biāo)估算主要程序如下：

三、導(dǎo)線網(wǎng)平差計(jì)算示例

1)導(dǎo)線網(wǎng)平差過(guò)程主要包括：角度閉合差的計(jì)算與調(diào)整;導(dǎo)線邊坐標(biāo)方位角的計(jì)算;相鄰導(dǎo)線點(diǎn)之間的坐標(biāo)增量計(jì)算;坐標(biāo)增量閉合差的計(jì)算與調(diào)整;導(dǎo)線坐標(biāo)點(diǎn)的計(jì)算。

2)首先把已知點(diǎn)的坐標(biāo)和觀測(cè)數(shù)據(jù)(距離和角度)輸入到Excel中，并保存，如圖3所示。

圖3 導(dǎo)線網(wǎng)平差的已知點(diǎn)坐標(biāo)和觀測(cè)數(shù)據(jù)

3)單擊“數(shù)據(jù)導(dǎo)入”，選擇框中選擇編輯好的數(shù)據(jù)，再點(diǎn)擊計(jì)算數(shù)據(jù)，即可計(jì)算出待求坐標(biāo)，如圖4所示。

圖4 導(dǎo)線網(wǎng)平差的計(jì)算結(jié)果

4)單擊顯示導(dǎo)線，即可顯示。如圖5所示。

圖5 導(dǎo)線網(wǎng)圖形顯示

四、結(jié) 論

測(cè)繪科學(xué)是一門以大規(guī)模數(shù)據(jù)甚至是海量數(shù)據(jù)處理、分析與應(yīng)用為基礎(chǔ)的學(xué)科，其各項(xiàng)具體工作都涉及大量的計(jì)算。Matlab在大規(guī)模數(shù)據(jù)處理特別是矩陣運(yùn)算方面具有其他程序設(shè)計(jì)語(yǔ)言難以比擬的優(yōu)越性。通過(guò)本文數(shù)據(jù)處理與分析，筆者得出以下結(jié)論：

1)利用Matlab強(qiáng)大的矩陣運(yùn)算功能，可以簡(jiǎn)單快捷地實(shí)現(xiàn)導(dǎo)線網(wǎng)的間接平差計(jì)算，相信對(duì)于其他的平差計(jì)算，Matlab一樣可以簡(jiǎn)單快捷的實(shí)現(xiàn)。

2)Matlab具有數(shù)據(jù)可視化的功能，利用強(qiáng)大的繪圖工具及繪圖函數(shù)，可以計(jì)算和繪制復(fù)雜的多維曲線。

3)Matlab是適合個(gè)人應(yīng)用的強(qiáng)有力的面向矩陣的高級(jí)程序設(shè)計(jì)語(yǔ)言，具有與其他語(yǔ)言編寫(xiě)的程序結(jié)合和輸入輸出格式化數(shù)據(jù)的能力。隨著Matlab的普及應(yīng)用，它在測(cè)繪領(lǐng)域中發(fā)揮更大的作用。

[1]周建興，豈興明，矯津毅，等．MATLAB從入門到精通[M]．北京：人民郵電出版社，2008．

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

[3]崔利，武文波．測(cè)繪領(lǐng)域中MATLAB的應(yīng)用[J]．遼寧工程技術(shù)大學(xué)學(xué)報(bào)：自然科學(xué)版，2005(Z1)：47-49．

[4]高寧，高彩云．MATLAB在測(cè)繪領(lǐng)域中的應(yīng)用[J]．平頂山工學(xué)院學(xué)報(bào)，2008，17(1)：28-31．

[5]王繼剛．關(guān)于在測(cè)繪工程專業(yè)中設(shè)置Matlab課程的探討[J]．北京測(cè)繪，2009(2)：76-78．

[6]趙亞紅，郝延錦，徐明俠．Matlab在測(cè)繪工程領(lǐng)域中的幾點(diǎn)應(yīng)用[J]．地礦測(cè)繪，2011(1)：38-40．