燈湖礦區(qū)控制網(wǎng)設計

1、 燈湖礦區(qū)控制網(wǎng)設計一、工程控制網(wǎng)設計任務書燈湖礦區(qū)控制網(wǎng)設計任務通知書XX測繪隊：燈湖礦區(qū)位于XX省XX煤田東北域。本礦區(qū)已列入國家經(jīng)濟建設計劃，準備進行重點開發(fā)。根據(jù)建設計劃，首先需施測80平方公里的1：2000地形圖。圖一為該礦區(qū)1：50000地形圖，其中虛線范圍為1：2000測圖區(qū)。為了加快建井速度，需在豎井間進行對向掘進，最大貫通距離為3km，礦井最大單翼長度5km，井深一般不超過800m。為滿足礦山設計、建井、生產(chǎn)三階段測繪各種比例尺地形圖、井巷貫通以及工業(yè)場地施工測量的需要，應在全礦區(qū)建立統(tǒng)一的具有足夠精度密度的平面控制網(wǎng)與高程控制網(wǎng)。要求你隊在任務通知書下達后一周內(nèi)完成控制網(wǎng)技

2、術(shù)設計，一月內(nèi)完成建網(wǎng)任務。XX省煤炭管理局 年 月 日二、測區(qū)自然地理條件及已有測繪成果資料（一）自然地理條件1地理概況本測區(qū)中心位置為東經(jīng)11713，北緯3412。測區(qū)地面高程為+30+244m。位于XX省XX縣境內(nèi)，屬于華東丘陵地帶，南部及中部為山地。其上林木繁茂，不利通視，其余為平原，密布村莊和廠礦，建筑群遍地，影響觀測。燈湖位于測區(qū)東部，湖面開闊，中有大堤。2交通情況測區(qū)內(nèi)有鐵路支線通過，公路四通八達，村間大道可行汽車，交通方便。3氣候情況測區(qū)地處華東近海地區(qū)，氣候宜人。全年平均降雨量為1040mm，雨量集中于6、7、8三個月份。全年平均氣溫+15，夏季氣溫較高，一般為30左右，冬季

3、有雪，但不寒冷，最低溫度為-5，凍土線深度為0.1m。全年平均風力為23級，夏季略受臺風影響。宜于野外作業(yè)時間為311月份，年平均作業(yè)時間利用率為21天/月。4居民及居民點測區(qū)內(nèi)地少人多，人口稠密，多為漢族，少數(shù)為回族。測量作業(yè)所需人力、物力、財料及食宿均可就地解決。（二）已有測繪成果資料1三角網(wǎng)成果及其精度測區(qū)內(nèi)及附近有國家二等網(wǎng)點三個：玉山、太山、廣具山。系1959年由XX測繪隊施測，作業(yè)所依據(jù)的規(guī)范為一、二、三、四等三角測量細則（1958年）。三點標石保存完好。坐標系統(tǒng)為1954年北京坐標系，三度分帶，中央子午線為117。2高程網(wǎng)成果及其精度國家等水準路線由西向東橫穿測區(qū)北部。根據(jù)二等水

4、準路線圖，本測區(qū)內(nèi)及附近應有三個二等水準點，依點之記只找到其中兩個：-2，-4，標石保存完好。該二等水準路線系國家測繪隊于1979年施測。施測精度及埋石質(zhì)量均符合規(guī)范要求。高程為1956年黃海高程系統(tǒng)。3地形圖資料1984年XX航測隊航測本地區(qū)1：50000地形圖，成圖質(zhì)量良好，可供技術(shù)設計圖上選點之用。二等三角點、水準點坐標和高程：點 名 等級 x（m） y（m） H（m）玉 山 二 3793400 20506000 195.4廣具山 二 3790725 20515900 135.0太 山 二 3786800 20510250 244.0-4 二 35.000-2 二 46.000三、測區(qū)已

5、有地形圖（四）坐標系統(tǒng)的選擇及處理方法的論證、起始數(shù)據(jù)的配置和處理工程平面控制網(wǎng)的坐標系統(tǒng)以與國家大地網(wǎng)的坐標系統(tǒng)取得一致為最理想，此時，工程控制網(wǎng)便能夠成為國家大地網(wǎng)的組成部分。但是要使工程控制網(wǎng)完全采用國家統(tǒng)一坐標系統(tǒng)，必須具備下列條件之一： (1)工程區(qū)域位于高斯投影3帶中央子午線附近，同時測區(qū)平均高程面接近平均海水面； (2)雖然不能符合上述條件，但根據(jù)測區(qū)的y坐標和高程H計算的相對長度綜合變形小于1/40000。事實上很多工礦測區(qū)不能滿足上列條件之一。此時應該結(jié)合測區(qū)已有國家高級點的數(shù)目、位置、精度等情況作下列選擇： 1. 如果測區(qū)內(nèi)國家大地控制點的精度符合現(xiàn)行規(guī)范對起算數(shù)據(jù)的精度要

6、求，且具有必要的或多余的(四個或四個以上)起算數(shù)據(jù)時，可以選擇測區(qū)局部坐標系統(tǒng)(或稱測區(qū)獨立坐標系統(tǒng))，選擇方案可為： (1) 選擇抵償高程面作為投影面，按高斯正形投影3帶計算平面直角坐標； (2) 保持國家統(tǒng)一的橢球面作投影面，選擇抵償投影帶，按高斯正形投影計算平面直角坐標； (3) 選擇測區(qū)平均高程面作為投影面，通過測區(qū)中心的子午線作為中央子午線，按高斯正形投影計算平面直角坐標。2. 測區(qū)內(nèi)雖有多個(兩個及兩個以上)國家大地控制點，但其精度不詳或低于現(xiàn)行規(guī)范對起算數(shù)據(jù)的精度要求。此時可選擇其中一個標石完好，位置適中的國家點作原點，利用其坐標(x,y)和一條邊的坐標方位角T作為起算數(shù)據(jù)，另外

7、重新測定一條或一條以上的起算邊長，將測距結(jié)果化算至測區(qū)平均高程面。按這種方法建立的坐標系統(tǒng)雖屬礦區(qū)局部坐標系統(tǒng)，但其原點的坐標和方位角仍保持了國家統(tǒng)一坐標系統(tǒng)。這種方法已為許多工礦區(qū)和城市所采用，它是彌補國家高級網(wǎng)某些缺陷的好方法。該礦區(qū)坐標系統(tǒng)選擇過程：Ym=(y1+y2)/2;Hm=(H1+H2)/2經(jīng)過計算可得：Ym=-11.5km； Hm=137m所以： =(0.00123*(-11.406)2-15.7*137)/100000=1/503461/40000所以直接采用3度帶國家高斯投影。（五）平面控制網(wǎng)布設方案1）首級GPS網(wǎng)（E級）設計：首級控制網(wǎng)的等級與要求精度及測區(qū)面積有關(guān)。一

8、個點能夠控制多大的面積，或者某一面積范圍內(nèi)應有多少個點，這是點的密度問題。在城市測量規(guī)范和工程測量規(guī)范中已有各等級的邊長系列，實際上已明確了點的密度。 考慮到工程地區(qū)大都采用GPS網(wǎng)作為首級控制網(wǎng)，且GPS網(wǎng)大多布設成類似三角網(wǎng)的形式，。結(jié)合本次任務測區(qū)內(nèi)實際情況，選擇平均邊長為2.5km。按四等GPS網(wǎng)的主要技術(shù)要求布設，在布設過程中，也應注意一點：雖然GPS測量時不要求相鄰點之間通視，布網(wǎng)方式也非常靈活，但考慮到加密導線測量的需要，每個GPS點應與相鄰的12個點通視。綜合考慮以上各種因素，測區(qū)面積約為75平方千米，每個控制點控制的面積約為5.6平方千米(按公式S=0.87D2粗略計算)，所

9、以在測區(qū)內(nèi)布設E級GPS控制點11個，與3個已知點組成由14個點構(gòu)成的GPS平面控制網(wǎng)。設計草圖見下頁：1，設計簡要說明：四臺GPS同時觀測，共8個時段，可獲獨立基線24條，具體可見下頁圖。2，精度評定結(jié)果：運用matlab軟件編程估算其精度，具體代碼見matlab7.0程序附件gps2.m，下表（一）是各控制點點位精度和誤差橢圓參數(shù)。（A表示橢圓方向角，E,F分別表示極大值和極小值誤差，Mx,My分別表示X方向和Y方向的誤差，M表示點位誤差）；表（二）是相對點位精度和誤差橢圓參數(shù)。（T表示橢圓方向角，E,F分別表示極大值和極小值誤差，s是基線長，Ms表示基線誤差，Ms/S表示相對邊長誤差，單

10、位為萬分之一）協(xié)因素陣：由精度估算公式：可得：點位中誤差：表（一） 點號 A(dms) E(mm) F(mm) Mx (mm) My (mm) M(mm)相對點位中誤差：表（二） 基線號 T(dms) E(mm) F(mm) S(km) Ms (mm) Ms/S (1/10000)-1-1 1-22-3 -3-3-3-76-75-64-5-1-41-52-62-73-7-3-98-97-86-86-114-11-2-11-2-108-109-10-2-83，由表（一）可知最弱點中誤差為E09號點，點位精度29.7mm(小于30mm),符合要求，相對點位精度最大方向誤差值為24.6mm（9-10

11、）,相對邊長精度最大值為-3-7號基線(-3號點為廣具山)，其相對精度為1/192382，均符合規(guī)范精度要求。2）一級加密網(wǎng)設計： 根據(jù)一級導線的規(guī)范要求，結(jié)合測區(qū)實際情況，布設成共8個點（包括兩個E級已知點E004,E011）的附合導線網(wǎng)，平均邊長約為0.58km，網(wǎng)型布設情況見前面首級GPS網(wǎng)型設計草圖。已知點和布設的E級GPS控制點的圖上粗略量測坐標如下表： （1954年西安坐標系、中央子午線經(jīng)度117度）序號點號轉(zhuǎn)換后X坐標轉(zhuǎn)換后Y坐標1-1(玉山)3793400205060002-2(E04)3790590205065903F013792880205063304F0237924002

12、05065605F033791910205064506F043791500205066707F05379100020506760精度評定結(jié)果：運用matlab軟件編程估算其精度，具體計算過程及代碼實現(xiàn)見matlab7.0程序附件daoxian2.m，下表（一）是各控制點點位精度和誤差橢圓參數(shù)。（A表示橢圓方向角，E,F分別表示極大值和極小值誤差，Mx,My分別表示X方向和Y方向的誤差，M表示點位誤差）；表（二）是相對點位精度和誤差橢圓參數(shù)。（T表示橢圓方向角，E,F分別表示極大值和極小值誤差，s是基線長，Ms表示基線誤差，Ms/S表示相對邊長誤差，單位為萬分之一）協(xié)因素陣：由精度估算公式：可得

13、：點位中誤差：表（一） 點號 A(dms) E(mm) F(mm) Mx (mm) My (mm) M(mm)相對點位中誤差：表（二） (-1代表E011點位, -2代表E004點位) 基線號 T(dms) E(mm) F(mm) S(km) Ms (mm) Ms/S (1/10000)-1-11-22-33-44-5-2-5由表（一）可知最弱點中誤差為F03號點，點位精度約為13.9mm(小于20mm),符合要求，相對點位精度最大方向誤差值約為9.2mm（1-2）,相對邊長精度最大值為-1-1（即玉山-F01）號基線，其相對精度為1/157673，均符合規(guī)范精度要求。(六)高程控制網(wǎng)設計：

14、根據(jù)兩個個已知的二等高程控制點II-2，II-4，根據(jù)工程實際要求，布設四等水準控制網(wǎng)。1，兩個已知控制點，,水準高程具體如下表：點名H/m備注II-246.000高斯投影L0=1171954國家高程基準II-435.0002，水準網(wǎng)型布設草圖見下頁：點位精度評定： 1，根據(jù)布網(wǎng)，找出高程精度弱點和結(jié)點，用邊長定權(quán)，估算出其點位高程精度，所依據(jù)依據(jù)公式：2，編寫程序計算，具體見附錄matlab7.0程序附件shuizhun2.m，精度評定結(jié)果：協(xié)因素陣：由高差中誤差公式可得：點號 點位精度3，由估算結(jié)果可知，最弱點為F19號點，點位精度約為19.1mm,小于限制誤差20mm，符合規(guī)范精度要求。

15、（七）技術(shù)依據(jù)及作業(yè)方法1. 測量規(guī)范及技術(shù)規(guī)程（1）、全球定位系統(tǒng)（GPS）測量規(guī)范，GB/T18314-2001。（2）、全球定位系統(tǒng)城市測量技術(shù)規(guī)范，CJJ73-97。（3）、工程測量規(guī)范，GB12897-91。（4）、城市測量規(guī)范，CJJ8-99。(5 )、國家三、四等水準測量規(guī)范，GB12898-91。(6 )、1：500、1：1000、1：2000地形圖圖式，GB/T7929-1995。2.標石圖（1）A.2.5 四等導線/三角形網(wǎng)/GPS平面控制點標石埋設規(guī)格應符合圖規(guī)定。圖A.2.5 四等導線/三角形網(wǎng)/GPS平面控制點標石埋設圖（單位：mm）注：1土面；2素土；3凍土線； 4

16、貧混凝土（2）A.2.6 一、二級導線/五等GPS平面控制點標石埋設規(guī)格應符合圖規(guī)定。圖A.2.6 一、二級導線/五等GPS平面控制點標石埋設圖（單位：mm）1 土面；2素土；3貧混凝土（3）A.3.3 四等水準點標石埋設規(guī)格應符合圖的規(guī)定。圖A.3.3 四等水準點標石埋設圖（單位：mm）1地面；2素土；3凍土線；4貧混凝土3，儀器的選擇及檢驗項目要求（1）儀器：GPS單頻接收機，全站儀，DS3水準儀（2）檢驗項目： 水準儀的光學測微器效用正確性的檢驗和分劃值的確定；（3）視準軸與水準軸相互關(guān)系的檢驗；（4）水準尺檢驗：水準尺分劃面彎曲差的檢驗；（5）水準標尺每千米分劃間隔真長的測定；（6）一

17、對水準標尺零點差及基，輔分劃讀數(shù)差常數(shù)的測定；4.觀測方法及各項限差（1）E級GPS觀測方法及各項限差要求：1）觀測組應嚴格按調(diào)度表規(guī)定的時間進行作業(yè)，以保證同步觀測同一衛(wèi)星組。當情況有變化需修改調(diào)度計劃時，應經(jīng)作業(yè)隊負責人同意，觀測組不得擅自更改計劃。2）接收機電源電纜和天線應連接無誤，接收機預置狀態(tài)應正確，然后方可啟動接收機進行觀測。3）各觀測時段的前后各量取天線高一次，兩次量高之差不大于3mm。取平均值作為最后天線高，記錄在手簿。若互差超限，應查明原因，提出處理意見記入手簿備注欄中。天線高是指觀測時天線平均相位中心至測站中心標志面的高度，分為上、下兩段：上段是指相位中心至天線底面的高度，

18、這是常數(shù)hc，由廠家給出；下段是從天線底面至測站中心標志面的高度，由觀測員在現(xiàn)場采用傾斜測量方法直接量取。具體方法是：從三腳架三個空檔（互成120）測量天線底盤下表面至測站中心標志面的距離，互差應小于3mm，取平均值為L，天線底盤半徑為R，再利用廠家提供的hc，按天線高 求出。4）接收機開始記錄數(shù)據(jù)后，作業(yè)人員可使用專用功能鍵選擇菜單，查看測站信息、接收衛(wèi)星數(shù)、衛(wèi)星號、各通道信噪比、實時定位結(jié)果及存貯介質(zhì)記錄情況等。5）儀器工作正常后，作業(yè)員及時（每隔15min）逐項填寫測量手簿中各項內(nèi)容。6）一個時段觀測過程中不得進行以下操作：關(guān)閉接收機以重新啟動；進行自測試（發(fā)現(xiàn)故障除外）；改變衛(wèi)星截止高

19、度角；改變數(shù)據(jù)采樣間隔；改變天線位置；按動關(guān)閉文件和刪除文件等功能。7）觀測員在作業(yè)期間不得擅自離開測站，并防止儀器受震動和被移動，防止人和其它物體靠近天線，遮擋衛(wèi)星信號。8）接收機在觀測過程中不應在接收機近旁使用對講機和手機等通訊設備；雷雨過境時應關(guān)機停測，并卸下天線以防雷擊。9）觀測中應保證接收機工作正常，數(shù)據(jù)記錄正確，每日觀測結(jié)束后，應及時將數(shù)據(jù)下載到計算機硬、軟盤上，確保觀測數(shù)據(jù)不丟失。表1GPS網(wǎng)的主要技術(shù)要求級別平均距離(km)(mm)(110-6)最弱邊相對中誤差E級0.2510201/45000表2E級GPS測量作業(yè)的基本技術(shù)要求級別衛(wèi)星截止高度角()有效觀測衛(wèi)星數(shù)平均重復設站

20、數(shù)時段長度(min)數(shù)據(jù)采樣間隔(s)PDOP值E級1541.6401510表3同步環(huán)坐標分量及環(huán)線全長相對閉合差的規(guī)定（110-6）等級限差類型E 級坐標分量相對閉合差6.0環(huán)線全長相對閉合差10.0（2）一級導線觀測方法及各項限差要求：1，左角采用前-后-后-前，右角采用后-前-前-后，2，水平角采用測回法，豎直角采用中絲法，三絲法3，導線施測采用三聯(lián)腳架全園觀測法施測，水平角觀測的技術(shù)要求按工程測量規(guī)范GB50026-932.3.12.3.10執(zhí)行。4，測角導線轉(zhuǎn)折角有左角和右角之分。在導線前進方向左側(cè)的水平角稱為左角，右側(cè)的水平角稱為右角。閉合導線一般測其內(nèi)角，在公路測量中，附和導線一

21、般測右角，注意全線應統(tǒng)一。各等級的導線測角要求，應滿足規(guī)范。5，測邊光電測距，光電測距導線邊采用單向或往返觀測，導線邊長均觀測2測回，每測回4次讀數(shù)，一測回內(nèi)讀數(shù)較差應小于5mm,單程測回間較差應小于10mm。表一，光電測距導線的水平角技術(shù)要求等級儀器等級測回數(shù)左+右-360之差 /上下半測回之差/測角中誤差/方位角閉合差/一級I級210105表二，光電測距導線的豎直角技術(shù)要求等級儀器等級測回數(shù)豎直角互差/指標差互差/一級I級21010表三，光電測距導線的測距技術(shù)要求 等級儀器等級導線長度/km測回數(shù)測距中誤差/mm相對誤差導線全長相對閉合差一級I級44151/30 0001/14 000（3

22、）四等水準觀測方法及各項限差要求：5，概算內(nèi)容及平差方法：（1）導線測量概算得主要內(nèi)容：作業(yè)測量觀測數(shù)據(jù)的檢查，將觀測值劃算到標石中心，將觀測值劃算到橢球面上和高斯投影平面上。（2）GPS控制網(wǎng)的概算：數(shù)據(jù)處理及基線解算，觀測結(jié)果的外業(yè)的檢核。（3）水準測量的概算得主要內(nèi)容：觀測高差各項改正數(shù)（包括水準標尺每千米長度改正數(shù)，水準面不平行改正數(shù)）的計算和水準點概略高程表的編制。導線控制網(wǎng)，水準網(wǎng)采用間接平差方法，GPS控制網(wǎng)用GPS基線向量網(wǎng)三維無約束平差。6，提供成果資料（1）GPS網(wǎng)概算結(jié)果及平差結(jié)果；（2）一級導線網(wǎng)概算結(jié)果及平差結(jié)果；（3）一級導線網(wǎng)概算結(jié)果及平差結(jié)果。.（八）工作量綜合

23、計算及工作進程計劃表1，查找資料，考察測區(qū)：5天；2，控制網(wǎng)優(yōu)化設計：一周7天；3，準備相關(guān)材料及租借儀器和組織人員：一周7天；4，實地施測：8天；5，數(shù)據(jù)處理及成果提交：3天；共計：5+7+7+8+3=30天，由于天氣等因素影響，再加上彈性時間5天，合計35天。（九）需要設備及經(jīng)費預算1，儀器：GPS單頻接收機四臺，全站儀，DS3水準儀，光電測距儀等。2，粗略經(jīng)費預算：（1）、設計費和施測費：10萬元 （2）、材料費用：2萬 （3）、人員工資：9萬 （4）、儀器租用費：2萬 （5）、交通等其他費用：2萬 工程預計費用合計：25萬元 （十）上交的資料清單1，控制測量技術(shù)設計書；2，外業(yè)觀測原始

24、資料；3，內(nèi)業(yè)計算成果表及控制網(wǎng)圖。附錄matlab程序代碼：1, gps2.mB=zeros(48,22);JK=1,99;1,2;2,3;3,99;7,99;6,7;5,6;4,5;4,99;1,5;2,6;2,7;3,7;9,99;8,9;7,8;6,8;6,11;4,11;11,99;10,99;8,10;9,10;8,99; %基線向量i=1; while 1 %給系數(shù)陣B陣賦值 if i=23; break; endj=JK(i,1);k=JK(i,2); B(i-1)*2+1,(j-1)*2+1)=1; B(i-1)*2+2,(j-1)*2+2)=1; if k=99 a=1;

25、else B(i-1)*2+1,(k-1)*2+1)=-1; B(i-1)*2+2,(k-1)*2+2)=-1; end i=i+1; end S=0.64*4.6,4.6,4.0,4.0,4.2,4.2,3.8,3.8,5.0,5.0,3.7,3.7,2.2,2.2,4.2,4.2,4.3,4.3,3.8,3.8,3.5,3.5,4.0,4.0,3.9,3.9,4.1,4.1,3.6,3.6,3.0,3.0,3.6,3.6,3.9,3.9,4.1,4.1,4.0,4.0,4.0,4.0,2.6,2.6,3.5,3.5,3.7,3.7; %0.64是從圖上距離轉(zhuǎn)換到實際距離的參數(shù) MO=27;

26、 C=(S.*10).2.+100)./(MO2);%以2.5千米為單位權(quán)中誤差(其基線向量中誤差為27mm) Q=diag(C); P=inv(Q);%權(quán)陣 BTPB=B*P*B; QXX=inv(BTPB);%協(xié)因素陣 D=zeros(11,6);%點位精度陣 i=1; while 1 if i=12 disp(D); break; end H=sqrt(QXX(2*(i-1)+1,2*(i-1)+1)-QXX(2*(i-1)+2,2*(i-1)+2)2+4*QXX(2*(i-1)+1,2*(i-1)+2)2); A=90.0;%誤差橢圓的最大值方位角 E=MO*sqrt(0.5*(QXX

27、(2*(i-1)+1,2*(i-1)+1)+QXX(2*(i-1)+2,2*(i-1)+2)+H); F=MO*sqrt(0.5*(QXX(2*(i-1)+1,2*(i-1)+1)+QXX(2*(i-1)+2,2*(i-1)+2)-H); Mx=MO*sqrt(QXX(2*(i-1)+1,2*(i-1)+1); My=MO*sqrt(QXX(2*(i-1)+2,2*(i-1)+2); M=MO*sqrt(QXX(2*(i-1)+1,2*(i-1)+1)+QXX(2*(i-1)+2,2*(i-1)+2);%點位誤差 D(i,:)=A,E,F,Mx,My,M;%誤差橢圓及點位誤差陣 i=i+1;

28、end DX=ones(24,6); i=1; while 1 if i=25 disp(DX); break; end j=JK(i,1);k=JK(i,2); if k=99 Qax=QXX(2*(j-1)+1,2*(j-1)+1); Qay=QXX(2*(j-1)+2,2*(j-1)+2); Qaxy=QXX(2*(j-1)+1,2*(j-1)+2); else Qax=QXX(2*(j-1)+1,2*(j-1)+1)+QXX(2*(k-1)+1,2*(k-1)+1)-2*QXX(2*(j-1)+1,2*(k-1)+1); Qay=QXX(2*(j-1)+2,2*(j-1)+2)+QXX

29、(2*(k-1)+2,2*(k-1)+2)-2*QXX(2*(j-1)+2,2*(k-1)+2); Qaxy=QXX(2*(j-1)+1,2*(j-1)+2)+QXX(2*(k-1)+1,2*(k-1)+2)-QXX(2*(j-1)+1,2*(k-1)+2)-QXX(2*(k-1)+1,2*(j-1)+2); end H=sqrt(Qax-Qay)2+4*Qaxy2); E=MO*sqrt(0.5*(Qax+Qay+H); F=MO*sqrt(0.5*(Qax+Qay-H); T=A; %由于Qaxy=0; s=S(1,2*(i-1)+1);%邊長 Ms=sqrt(cos(T)2*Qax+(s

30、in(T)2*Qay+sin(2*T)*Qaxy);%邊長中誤差 f=(100*Ms)/s;%相對邊長中誤差（單位為萬分之一） DX(i,:)=T,E,F,s,Ms,f;%相對誤差橢圓陣 i=i+1; end 2, daoxian2.mB=zeros(13,10); % 給誤差方程系數(shù)陣的元素賦0,行數(shù)是觀測值個數(shù),列數(shù)是未知數(shù)個數(shù).JKH=-2,-1,1;-1,1,2;1,2,3;2,3,4;3,4,5;4,5,-2;5,-2,-1;JK=1,-1;1,2;2,3;3,4;4,5;5,-2;XYJS=3792880,20506330;3792400,20506560;3791910,2050

31、6450;3791500,20506670;3791000,20506760;3793400,20506000;3790590,20506590;ro=2062.65;i=1; while 1 if i=8;%觀測角的個數(shù)是7,因此i=9時中斷 break; endj=JKH(i,1);k=JKH(i,2);h=JKH(i,3); if j=0 j=5+abs(j);%未知點個數(shù)是6個 end if k=0 k=5+abs(k);%未知點個數(shù)是6個 end if h=0 h=5+abs(h);%未知點個數(shù)是6個 end dxjk=XYJS(k,1)-XYJS(j,1);dyjk=XYJS(k,

32、2)-XYJS(j,2); dxjh=XYJS(h,1)-XYJS(j,1);dyjh=XYJS(h,2)-XYJS(j,2); sjk0=dxjk*dxjk+dyjk*dyjk; sjh0=dxjh*dxjh+dyjh*dyjh; j=JKH(i,1);k=JKH(i,2);h=JKH(i,3); if j=0 & k=0 A=0; end if j=0 B(i,k*2-1)=-ro*dyjk/sjk0;B(i,k*2)=ro*dxjk/sjk0; end if k=0 B(i,j*2-1)=ro*dyjk/sjk0;B(i,j*2)=-ro*dxjk/sjk0; end if k=0 &

33、j=0 B(i,j*2-1)=ro*dyjk/sjk0;B(i,j*2)=-ro*dxjk/sjk0;B(i,k*2-1)=-ro*dyjk/sjk0;B(i,k*2)=ro*dxjk/sjk0; end if j=0 & h=0 A=0; end if j=0 B(i,h*2-1)=ro*dyjh/sjh0;B(i,h*2)=-ro*dxjh/sjh0; end if j=0 & h=0 & h=0 B(i,j*2-1)=B(i,j*2-1)-ro*dyjh/sjh0;B(i,j*2)=B(i,j*2)+ro*dxjh/sjh0;B(i,h*2-1)=ro*dyjh/sjh0;B(i,h*2

34、)=-ro*dxjh/sjh0; end i=i+1; end while 2 if i=14;%觀測角加上觀測邊的總數(shù)是13,因此i=14時中斷 break; endj=JK(i-7,1);k=JK(i-7,2);%減去觀測角的個數(shù)7 if j=0 j=5+abs(j);%未知點個數(shù)是5個 end if k=0 k=5+abs(k);%未知點個數(shù)是5個 end dxjk=XYJS(k,1)-XYJS(j,1);dyjk=XYJS(k,2)-XYJS(j,2); sjk0=sqrt(dxjk*dxjk+dyjk*dyjk); j=JK(i-7,1);k=JK(i-7,2); if j=0 &

35、k=0 A=0; end if j=0 B(i,k*2-1)=dxjk/sjk0;B(i,k*2)=dyjk/sjk0; end if k=0 B(i,j*2-1)=-dxjk/sjk0;B(i,j*2)=-dyjk/sjk0; end if k=0 & j=0 B(i,j*2-1)=-dxjk/sjk0;B(i,j*2)=-dyjk/sjk0;B(i,k*2-1)=dxjk/sjk0;B(i,k*2)=dyjk/sjk0; end i=i+1; end K=1,1,1,1,1,1,1,0.25,0.25,0.25,0.25,0.25,0.25;%粗略定權(quán)，平均邊長約為0.5km P=diag

36、(K);%觀測角和觀測邊的權(quán)陣 BTPB=B*P*B; QXX=inv(BTPB);%點位協(xié)因素陣 D=zeros(5,7);%點位精度陣 i=1; MO=5;%單位權(quán)中誤差為5mm while 1 if i=6 disp(D); break; end H=sqrt(QXX(2*(i-1)+1,2*(i-1)+1)-QXX(2*(i-1)+2,2*(i-1)+2)2+4*QXX(2*(i-1)+1,2*(i-1)+2)2); A=0.5*atan(2*QXX(2*(i-1)+1,2*(i-1)+2)/(QXX(2*(i-1)+1,2*(i-1)+1)-QXX(2*(i-1)+2,2*(i-1)

37、+2)+1.57;%誤差橢圓方向角0- E=MO*sqrt(0.5*(QXX(2*(i-1)+1,2*(i-1)+1)+QXX(2*(i-1)+2,2*(i-1)+2)+H); F=MO*sqrt(0.5*(QXX(2*(i-1)+1,2*(i-1)+1)+QXX(2*(i-1)+2,2*(i-1)+2)-H); Mx=MO*sqrt(QXX(2*(i-1)+1,2*(i-1)+1); My=MO*sqrt(QXX(2*(i-1)+2,2*(i-1)+2); M=MO*sqrt(QXX(2*(i-1)+1,2*(i-1)+1)+QXX(2*(i-1)+2,2*(i-1)+2);%點位誤差 D(i,:)=i,A,E,F,Mx,My,M;%誤差橢圓及點位誤差陣 i=i+1; end DX=ones(6,6); i=1; S=0.54,0.55,0.54,0.57,0.60,0.59,0.65;%邊長 while 1 if i=7 disp(DX); break; end j=JK(i,1);k=JK(i,2); % if k=-1|k=-2|k=-3 if k=0 k=6+abs(k);%未知點個