导线怎么用计算机平差,在计算机上实现导线网的自动条件平差
《在計(jì)算機(jī)上實(shí)現(xiàn)導(dǎo)線網(wǎng)的自動(dòng)條件平差》由會(huì)員分享,可在線閱讀,更多相關(guān)《在計(jì)算機(jī)上實(shí)現(xiàn)導(dǎo)線網(wǎng)的自動(dòng)條件平差(7頁(yè)珍藏版)》請(qǐng)?jiān)谌巳宋膸?kù)網(wǎng)上搜索。
1、維普資訊 http:/www.cqvip.com計(jì)算機(jī)應(yīng)用; 在計(jì)算機(jī)上實(shí)現(xiàn)導(dǎo)線網(wǎng)自 動(dòng) 條 件 平 差中國(guó)礦大北京研究生部 王紹林 崔明義引喜枝日4余枝。網(wǎng)絡(luò)圖中除樹(shù)枝之外的支、導(dǎo)線網(wǎng)絡(luò)圖很不規(guī)則,幾乎無(wú)規(guī)律可路為余枝。循,且一般導(dǎo)線網(wǎng)的未知數(shù)據(jù)又較多若5、度t 與某結(jié)點(diǎn)連通的支路個(gè)數(shù)稱(chēng)對(duì)導(dǎo)線網(wǎng)進(jìn)行間接平差,需列出許多誤差為該結(jié)點(diǎn)的度方程式 一個(gè)并不復(fù)雜的導(dǎo)線網(wǎng)絡(luò)圖,其 6、虛 擬支路:從導(dǎo)線網(wǎng)的一個(gè)已知誤差方程式繁多而又復(fù)雜,勢(shì)必給導(dǎo)線網(wǎng) 點(diǎn)到另一個(gè)已知點(diǎn)假設(shè)有一條支路連通,平差帶來(lái)困難。相比之下,導(dǎo)線網(wǎng)的條件 這條支路稱(chēng)為虛擬支路。用虛線表示。方程式在其數(shù)量和難易程度上都有一定的如圖 1。
2、(a)為一網(wǎng)絡(luò)圖 I、,優(yōu)越性,而且可在平差過(guò)程中通過(guò)條件不 、V、 為網(wǎng)絡(luò) 圖上的結(jié)點(diǎn),其中 I、符值的大小對(duì)導(dǎo)線網(wǎng)的精度得到直觀的估 為已知點(diǎn) 圖 1(b)為網(wǎng)絡(luò)圖上的一棵價(jià)。但是,由于導(dǎo)線網(wǎng)絡(luò)圖的不規(guī)剛性, 樹(shù), 8、 7、 8、 9、10為樹(shù)枝,則 1,給計(jì)算機(jī)自動(dòng)組成條件方程式帶來(lái)了困 6、 5、 4、 2為余 枝。結(jié)點(diǎn) I、 ,難。在這方面,一些測(cè)繪工作者已作了些 , v、的度 分別為 2, 4, 8、8有益的探索。筆者根據(jù)圖論和離散數(shù)學(xué)的 4、 4。支路 1為虛 擬支路理論,編制了一套計(jì)算機(jī)對(duì)導(dǎo)線網(wǎng)自動(dòng)條 在作了以上的約定后,首先在導(dǎo)線網(wǎng)件平差程序,經(jīng)測(cè)試,取得了較為滿意的 絡(luò)圖。
3、中尋找樹(shù) 。設(shè)這個(gè)導(dǎo)線 網(wǎng)絡(luò)圖中共有結(jié)果。現(xiàn)將其基 本算法介紹如下: 1個(gè)結(jié)點(diǎn),其中有in個(gè)已知點(diǎn),則這個(gè)導(dǎo)線網(wǎng)絡(luò)圖的樹(shù)就有a=n一1條樹(shù)枝網(wǎng)絡(luò)一 、 算法的基本思路 圖有b=nl一1條虛擬 支路 。如 圖 2 (a) 示,該網(wǎng)絡(luò)圖有 6個(gè)結(jié) 點(diǎn),其中 (1)、導(dǎo)線網(wǎng)絡(luò)圖是由導(dǎo)線點(diǎn)和導(dǎo)線邊組成(2) 、 (8) 為 8個(gè)已知點(diǎn), (4),的。為此,先作如下約定:, 5)、 (6)為 8個(gè) 未知 點(diǎn)在 已知 點(diǎn)1、結(jié)點(diǎn)t導(dǎo)線網(wǎng)絡(luò)圖中由兩條以上 1)、(2) 和 (2), (8)之 間,分的導(dǎo)線相交的導(dǎo)線點(diǎn)或已知點(diǎn)為結(jié)點(diǎn)。 別虛擬兩條支路 7, 8。先 在該網(wǎng)絡(luò)圖中2,樹(shù)t包含有導(dǎo)線網(wǎng)絡(luò)圖上的所有。
4、 尋找一棵樹(shù),這棵樹(shù)有a=61=5條結(jié)點(diǎn),且從一個(gè)結(jié)點(diǎn)到另一個(gè)結(jié)點(diǎn)僅有唯 樹(shù)枝 在尋找這棵樹(shù)之前,應(yīng)構(gòu)造一一 的一條通路,結(jié)點(diǎn)之間不能構(gòu)成閉合回 個(gè)初始陣。初始陣以數(shù)據(jù)文件 的形式給路的連通圖。它是導(dǎo)線網(wǎng)絡(luò)圖的子圖,亦 出。初始陣為n (n為結(jié) 點(diǎn)個(gè)數(shù))階矩陣,稱(chēng)為網(wǎng)絡(luò)圖的生成樹(shù)。阻下簡(jiǎn)稱(chēng)樹(shù)。每一行或每一列分別代表一個(gè)結(jié)點(diǎn)矩陣8,樹(shù)枝t構(gòu)成樹(shù)的所有支路稱(chēng)為樹(shù)的元素aIi為結(jié)點(diǎn)i到結(jié)點(diǎn)j的支路編號(hào)。若22維普資訊 http:/www.cqvip.com結(jié)點(diǎn)i到結(jié)點(diǎn) j順導(dǎo)線施測(cè)方向,a I為正 條連通支路,則 Ij為o。以匿 2 (a)l為號(hào),反之為負(fù)號(hào)l若結(jié)點(diǎn)i到結(jié)點(diǎn)j沒(méi)有一 例,該剮絡(luò)圖的初。
5、始陣為1234560701O0l 一70803旺l0800051(1)l 一100026 一 l1030一204l00 56 40I 工(b)圖1(b)圖2初始陣有幾個(gè)明顯的特征:i)主對(duì)角 每一列非零元素的個(gè)數(shù)恰好等于對(duì)應(yīng)結(jié)點(diǎn)線元素均為。ii)元素絕對(duì)值對(duì)稱(chēng)相等, 的度。有 了初始陣,即可求出導(dǎo)線網(wǎng)絡(luò)圖符號(hào)相反,即 lj=一 IIiii)每一行或 的樹(shù)。其通用算法見(jiàn)框圖 323維普資訊 http:/www.cqvip.com在圖 8中,若導(dǎo)線網(wǎng)絡(luò)圈中度數(shù)最大網(wǎng)絡(luò)圈有 T條支路 (包括虛 擬支路) ,則的結(jié)點(diǎn)不止一個(gè),則可取其中任一個(gè)作為關(guān)聯(lián)陣B應(yīng)有T列,每一列代表一條支D點(diǎn)路 B陣前a=n一。
6、1剮為樹(shù) 枝,后Ta歹u在找到樹(shù)以后,即可確定余棱。如翻為余枝。 B陣元素值bIj為十1、一1或o。2 (b) 即為所尋找的樹(shù),由此知道支路+1表示支路j方向遒向結(jié)點(diǎn) i;一1表示1, 5、 6為采枝。找出樹(shù)以后,構(gòu)造一支路 J方 向指向結(jié)點(diǎn)i,o表示支路 j與結(jié)關(guān)聯(lián)矩陣B,B陣有n行 ,每一行代表一點(diǎn)i無(wú)關(guān)。以圖 2 (a)為例,其關(guān)聯(lián) 陣個(gè)結(jié)點(diǎn),行序按結(jié)點(diǎn)編號(hào)排列假設(shè)導(dǎo)線為;圖82d維普資訊 http:/www.cqvip.com78824156+ 10O00+ 100一1+1+10000010 10000+ 10B =ooo+10 1o一1(2)00 1 1 + 10000000 10 。
7、1 + 1陣的每一列元素之和都等于零。這就是說(shuō), 陣是線性相關(guān)的。 陣 何一行元素的絕對(duì)值都等于其余行同列元素代數(shù)和的絕對(duì)值,但符號(hào)相反。因此,在 B陣中去掉任意一行,使其線性無(wú)關(guān),并不影響關(guān)聯(lián)陣所描述的信息不 妨去掉第一行 ,得B陣,78324156,一1+ 1+ 1000O0、1 0 10000+10lB 1 000+10101 I(8)l0011+1000f、0000 10 1+ 1 對(duì) B陣進(jìn)行 分塊,前 5列為 B ,對(duì)應(yīng)于樹(shù)枝部分后 8列為 B :,對(duì)應(yīng)于余枝部分,則有,B =(BlJ Blj(4)B1,陣為a階方陣,BJ 陣為 a(Ta)階 矩陣。由于 B, 陣的行對(duì)應(yīng)于結(jié)點(diǎn),列。
8、對(duì)應(yīng)于樹(shù)枝,所以,B,陣的每一行和每一列都有非零元素存在 ,且 B, 陣線性無(wú) 關(guān)。將 】B-,】的行列式運(yùn)行 (列 )按 代數(shù)余子式展開(kāi),lBl l=1(O)這就是說(shuō),無(wú)論導(dǎo)線網(wǎng)絡(luò)圖形如何 , B,陣都是滿秩陣,有逆陣 Bi;存在。令 A:J=一(BiB1j)(5) Al =E(E為 Ta階單位陣)(6 )A 是單位陣, B,的列數(shù)決定了單位陣的階數(shù),則有 -A =A L| A ,)(7)(7)式即為導(dǎo)線條件方程信息陣。由于余枝個(gè)數(shù)為:i=T a(8)式中T一導(dǎo)線網(wǎng)絡(luò)圖支路個(gè)數(shù),a一樹(shù)枝數(shù)每一個(gè)導(dǎo)線網(wǎng)絡(luò)圖,無(wú)論形狀如何復(fù)雜,其獨(dú)立閉合環(huán)個(gè)數(shù)總是等于余枝個(gè)數(shù)。對(duì) 圖 2(a)來(lái)說(shuō),i=Ta=8。
9、5=3,即圖 2(a)有 8個(gè)獨(dú)立的cI1合環(huán)。A陣為38階陣,每一行代表一個(gè)獨(dú)立的閉合環(huán)信息,即有t25維普資訊 http:/www.cqvip.com7832一 10 11A =1 01100O01在導(dǎo)線網(wǎng)條件平差中,各條件方程式歹立時(shí)的閉合環(huán)方向應(yīng)該一致,這時(shí)的閉合環(huán)方向與導(dǎo)線的施測(cè)方向無(wú)關(guān)。在 A陣中,每一行為一條件方程式的信息,元素 +1表示 閉合環(huán)的方向與對(duì)應(yīng)支路方向一致,元素一 l表示 閉合環(huán)的方向與對(duì)應(yīng)支路的方向相反,元素 0表示 閉合環(huán)不經(jīng)過(guò)該支路。由A陣參照 圖2(a)可知A =;一41G01:I10101 ,第一閉合環(huán)為逆時(shí)針?lè)较?#xff0c;即一7、一3,2、 ,第二閉合環(huán)為順時(shí)針。
10、方向,即一3、8、5、一4,第三閉合環(huán)為順時(shí)針?lè)较?#xff0c;日2、4,6。為解決閉合環(huán)方向的矛盾,程序可判斷各閉合環(huán)的方向,并對(duì)其中個(gè)別閉合環(huán)的方向作調(diào)整。在A陣中,只須將第一行的元素壘部改變符號(hào),即可解決方向的矛盾因此有t4l0 110O O)10A陣即為完整的條件方程式信息陣,且各閉合環(huán)方向一致 由于支路 7、 8為虛 擬支路,所以,第一個(gè)和第二個(gè)條件式均為附臺(tái)條件,第三個(gè)條件式為閉臺(tái)環(huán)條件。 6利用 A陣提供 的信息,程序可自動(dòng)完成對(duì)導(dǎo)線網(wǎng)的條件平差。二、導(dǎo)線網(wǎng)條件平差的數(shù)學(xué)模型OOl本程序采用近似平差方法t 即將角度條件和縱、橫坐標(biāo)條件分別解算平差。首先列出角度條件式,從而求出角度 改正數(shù),利。
11、用改正后的角度求各導(dǎo)線邊的坐標(biāo)方位角,然后分別列出縱、橫坐標(biāo)條件方程式,再分別求出縱、橫坐標(biāo)增量改正數(shù),最后算出各導(dǎo)線點(diǎn)的平差坐標(biāo),并根據(jù)選擇求出各導(dǎo)線邊的坐標(biāo)方位角和邊長(zhǎng),角度和縱 橫坐標(biāo)平差值函數(shù)中誤差,每公里點(diǎn)位中誤差等圖4現(xiàn)以圖 4為例 ,簡(jiǎn)述近似條件平差的數(shù)學(xué)模型。角度條件方程式為26一Vl Vj口一V 4a V!jVj2一Vl3V B。+ V + Vj5十V “ 十W 口一 0V 3+ V jD+V ebV l8一 V 7一V l6一 VV 14一V l5一 V 2+ W nI 0V Il+ V I2+ V 13+V 5a+ V ;b+V 】B+ V l7十V j8+V B4+ V。
12、 2口+V 2l十V 22一V 4b+ W 。j O式中W4l=(Bl+B】0+B +Bjl+Bl4+Bj8+B5 +(380。+(380。一日l(shuí)5)(380。一日2)一10x180。)一( 7一l一 8-2)+I1l380。W tI=Bj+BlB+ j 十B6b+ le+Bl7+ 18+(380。一日6b)+(380。一日1日)(380。 B3)一一10x 180。一 (ds 一9一。)+n x 380。W tI =BlI+Bl2十Bl3+BB +B5b+BlB+Bl7+BI8+Be +B20+B2l+B2j+(360。一B4b)一180。(12+2)+n3360。其法方程式為;N 。k+ 。
13、W = 0式中N = A P AW = gw ql W I W 。I式中可取PI1=E (單位陣)縱坐標(biāo)條件方程式:一Vl V l2+V I。+W Il= 0一V 8一 V d十 Vt5+ W lI= 0V 2+Vi 十 V B+W l= 0式中W ll= 一AxlAX2+x3一(xI X2)W II= 一Ax3一Ax+Ax 5一(Xjx3)W iA x 4+ Ax+ Ax B其法方程式為;Nk +W= 0式中N = A。Pi A:W :CW l W W 1維普資訊 http:/www.cqvip.com(11)(12)( )(14)(15)(16)(17)(18)27維普資訊 http:/w。
14、ww.cqvip.comP。可取橫坐標(biāo)條件方程式,V ,I V ,2+ V ,8PV 3一V ,+ V ,I(19)=V y+ V ,4+ V rB一式中一 Wl= 一 yI一 y2+ y8(yIy 2)W , = 一 y8一 y4+ y5(yj y3)(2O)W ,!= y2+ y+ y日其法方程式為,式中N, Kr w , 0(2王)+Jr,w wN 7w= A , P A!(22)l1=w T=(w= ,I W 7w rIPi可取O00P 一n E IyI精度評(píng)定公式與測(cè)量教科書(shū)中相同,在此不予贅述。源程序清單和平差算側(cè)從略,、 、,、 ,、 ,_(上接第21頁(yè))五、幾點(diǎn)結(jié)論通過(guò)上述分析和試驗(yàn),可以得出以下看法-1、壘盆地求參法比剖面法具有明顯的優(yōu)勢(shì),它不受主斷面和局部地質(zhì)采礦條件的變化影響,放寬了對(duì)地表移動(dòng)觀測(cè)站設(shè)置的要求,2、設(shè)置地表移動(dòng)現(xiàn) 站的兩大類(lèi)五種方法是可行的,而且基本上通過(guò)實(shí)地試驗(yàn)加以證實(shí)。除設(shè)置 在移動(dòng)盆地主斷面謝線外,必須采用壘盆地求參法求得地表移動(dòng)參數(shù),否則,很難達(dá)到設(shè)置觀測(cè)站的目的。8、采用本文提出的設(shè)置地表移動(dòng)觀測(cè)站的方法,不僅可以省勞力、省經(jīng)費(fèi)、省時(shí)間,更重要的是可以緩和觀測(cè)站設(shè)置中的工農(nóng)關(guān)系,為廣泛開(kāi)展地表移動(dòng)基本規(guī)律的研究,安全,合理地解決 “三下采煤?jiǎn)栴}創(chuàng)造更為有利的條件。28。
總結(jié)
以上是生活随笔為你收集整理的导线怎么用计算机平差,在计算机上实现导线网的自动条件平差的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問(wèn)題。
- 上一篇: 公里导线时有没有计算机类计算方法,基于M
- 下一篇: C++/测绘附和导线测量源码