礦大測(cè)繪專業(yè)實(shí)習(xí)報(bào)告

1、大地測(cè)量基礎(chǔ)與衛(wèi)星導(dǎo)航與定位技術(shù)教學(xué)實(shí)習(xí)報(bào)告書(shū) 學(xué)院 環(huán)境與測(cè)繪學(xué)院 班級(jí) 測(cè)繪11-3班 學(xué)號(hào) 07113021 姓名 田孟浩 2013-07-23一、實(shí)習(xí)目的與要求1、進(jìn)一步鞏固、加深大地測(cè)量學(xué)、GPS的有關(guān)理論知識(shí)。2、熟練掌握大地測(cè)量學(xué)、GPS所使用的儀器及觀測(cè)方法，提高控制測(cè)量實(shí)踐能力。3、培養(yǎng)和提高解決實(shí)際問(wèn)題及組織測(cè)繪生產(chǎn)的管理能力。4、進(jìn)行熱愛(ài)測(cè)繪事業(yè)的教育。二、實(shí)習(xí)內(nèi)容與任務(wù)1、進(jìn)行選點(diǎn)、埋石的實(shí)踐訓(xùn)練。2、合作共同完成若干GPS點(diǎn)的靜態(tài)相對(duì)定位觀測(cè)。3、全組完成全站儀、一臺(tái)水準(zhǔn)儀及一對(duì)水準(zhǔn)標(biāo)尺的檢驗(yàn)和記錄。4、全組共同完成一條導(dǎo)線的水平角、垂直角的觀測(cè)工作。每人獨(dú)立完成至少

2、2個(gè)測(cè)站點(diǎn)的水平角和垂直角的觀測(cè)工作，每人獨(dú)立完成至少2個(gè)測(cè)站點(diǎn)的全部記錄工作。5、全組共同完成一條導(dǎo)線的測(cè)距工作。（如用全站儀可同時(shí)進(jìn)行角度與距離的觀測(cè)）6、全組共同完成一條導(dǎo)線上若干測(cè)段的二等水準(zhǔn)觀測(cè)工作。每人都有不少于10站的觀測(cè)和記錄工作。7、每人獨(dú)立完成一份平面控制網(wǎng)和高程控制網(wǎng)的概算及平差計(jì)算。8、每人編寫一份實(shí)習(xí)報(bào)告書(shū)。三. 實(shí)習(xí)的組織和安排：測(cè)繪11-3班共分為五組，每組共6人，然后每組出一人參加由全專業(yè)共同組成的G-S組測(cè)量任務(wù)。每組設(shè)正副組長(zhǎng)各一人，正組長(zhǎng)負(fù)責(zé)分工安排全組的實(shí)習(xí)工作，由副組長(zhǎng)兼任儀器保管員及時(shí)檢查督促安全工作。實(shí)習(xí)宏觀上由教師統(tǒng)一指揮，班干部積極協(xié)助指導(dǎo)教師

3、做好各項(xiàng)工作。我是在測(cè)繪11-3班實(shí)習(xí)三組，兩個(gè)組長(zhǎng)，四個(gè)組員，共六個(gè)人，分別是田孟浩、王名億、蔣路、張學(xué)萬(wàn)、張寶科、張夢(mèng)虹。四、執(zhí)行規(guī)范名稱（1）全球定位系統(tǒng)（GPS）測(cè)量規(guī)范GB/T 18314-2009。（2）全球定位系統(tǒng)城市測(cè)量技術(shù)規(guī)程CJJ73-97。（3）工程測(cè)量規(guī)范GB50026-2007。（4）城市測(cè)量規(guī)范CJJ8-99。（5）國(guó)家一、二等水準(zhǔn)測(cè)量規(guī)范GB/T 12897-91。（6）國(guó)家三、四等水準(zhǔn)測(cè)量規(guī)范GB12898-91。（7）國(guó)家三角測(cè)量規(guī)范GB/T17942-2000。五、實(shí)習(xí)地區(qū)一般情況介紹1. 實(shí)習(xí)地區(qū)所在地理位置實(shí)習(xí)地區(qū)為礦大南湖校區(qū)教學(xué)與生產(chǎn)實(shí)習(xí)基地五號(hào)區(qū)域

4、。GPS點(diǎn)有兩個(gè)，分別為I24和I25點(diǎn)，基本上是繞著圖書(shū)館和鏡湖區(qū)周圍的閉合環(huán)。2. 交通情況我們的測(cè)區(qū)有三邊是公路行走方便，但車輛較多，對(duì)我們的觀測(cè)產(chǎn)生了一定影響。由于我們的測(cè)區(qū)四周與其他測(cè)區(qū)相連，導(dǎo)致我們的控制點(diǎn)經(jīng)常被其他組使用。不利于我們的持續(xù)觀測(cè)3. 地形地物（包括植被）概況測(cè)區(qū)有兩面有著明顯的地形起伏，高差較大，同時(shí)植被茂密，經(jīng)常出現(xiàn)不通視的問(wèn)題。并且測(cè)區(qū)內(nèi)有一條河，河岸上的樹(shù)有時(shí)會(huì)對(duì)測(cè)量造成一些不便。4. 已知控制點(diǎn)利用情況已知點(diǎn)有兩個(gè)I24和I25，這兩個(gè)點(diǎn)通過(guò)GPS組的觀測(cè)，用Compass軟件將基線進(jìn)行解算可以得到這兩個(gè)點(diǎn)的XY坐標(biāo)。六、控制網(wǎng)的布設(shè)、作業(yè)計(jì)劃和作業(yè)組織1.

5、 控制網(wǎng)（GPS網(wǎng)、導(dǎo)線網(wǎng)、水準(zhǔn)網(wǎng)）示意圖（1）GPS網(wǎng)示意圖：2、測(cè)區(qū)分布示意圖3、起算數(shù)據(jù)點(diǎn)號(hào)XYH（56高程）備注泉山（QSGY）3788109.4670514667.667755.950054坐標(biāo)E0883787813.059513642.03439.966454坐標(biāo)3787766.762513589.52039.96980坐標(biāo)S0036.607概略高程七、參加選點(diǎn)在實(shí)地確定控制點(diǎn)的具體位置，力爭(zhēng)做到點(diǎn)位安全可靠，便于觀測(cè)和保存，點(diǎn)間相互通視，控制網(wǎng)的圖形結(jié)構(gòu)良好。在實(shí)地設(shè)立測(cè)量標(biāo)志，埋設(shè)標(biāo)石。其幾何中心表示點(diǎn)位。八.儀器檢驗(yàn)1、全站儀檢驗(yàn)： （1）照準(zhǔn)部水準(zhǔn)軸應(yīng)垂直于豎軸的檢驗(yàn)和校正

6、檢驗(yàn)時(shí)先將儀器大致整平,轉(zhuǎn)動(dòng)照準(zhǔn)部使其水準(zhǔn)管與任意兩個(gè)腳螺旋的連線平行,調(diào)整腳螺旋使氣泡居中,然后將照準(zhǔn)部旋轉(zhuǎn)180度,若氣泡仍然居中則說(shuō)明條件滿足,否則應(yīng)進(jìn)行校正。校正的目的是使水準(zhǔn)管軸垂直于豎軸即用校正針撥動(dòng)水準(zhǔn)管一端的校正螺釘，使氣泡向正中間位置退回一半為使豎軸豎直，再用腳螺旋使氣泡居中即可此項(xiàng)檢驗(yàn)與校正必須反復(fù)進(jìn)行，直到滿足條件為止。 （2）十字絲豎絲應(yīng)垂直于橫軸的檢驗(yàn)和校正檢驗(yàn)時(shí)用十字絲豎絲瞄準(zhǔn)一清晰小點(diǎn)，使望遠(yuǎn)鏡繞橫軸上下轉(zhuǎn)動(dòng)，如果小點(diǎn)始終在豎絲上移動(dòng)則條件滿足否則需要進(jìn)行校正。校正時(shí)松開(kāi)四個(gè)壓環(huán)螺釘（裝有十字絲環(huán)的目鏡用壓環(huán)和四個(gè)壓環(huán)螺釘與望遠(yuǎn)鏡筒相連接。轉(zhuǎn)動(dòng)目鏡筒使小點(diǎn)始終在十

7、字絲豎絲上移動(dòng)，校好后將壓環(huán)螺釘旋緊。2、水準(zhǔn)儀及水準(zhǔn)尺檢驗(yàn)：（1） 光學(xué)測(cè)微器效用正確性的檢驗(yàn)和分劃值的測(cè)定步驟：架設(shè)水準(zhǔn)儀，將儀器視準(zhǔn)軸瞄準(zhǔn)水準(zhǔn)尺，使用傾斜螺旋使水準(zhǔn)氣泡影像精密符合，在一測(cè)回內(nèi)要嚴(yán)格保持傾斜螺旋的位置不變。旋轉(zhuǎn)測(cè)微器螺旋，使測(cè)微器讀數(shù)略大于零，指揮扶尺者將特制分劃尺上下移動(dòng)，直到特制分劃尺上某一分劃線被楔形絲平分。然后將特制分劃尺固定，并在一測(cè)回中保持此位置不變。往測(cè)旋進(jìn)測(cè)微螺旋，用楔形絲分別夾準(zhǔn)特制分劃尺上兩相鄰分劃線，讀取并記錄分劃線編號(hào)和相應(yīng)的測(cè)微器讀數(shù)。返測(cè)按相反的次序，旋出測(cè)微螺旋，使楔形絲分別夾準(zhǔn)特制分劃尺上兩相鄰分劃線，讀取并記錄分劃線編號(hào)和相應(yīng)的測(cè)微器讀數(shù)

8、。至此算作一個(gè)測(cè)回。每測(cè)完兩個(gè)測(cè)回后，應(yīng)略微移動(dòng)分劃尺或變更儀器高度，使每?jī)蓚€(gè)測(cè)回分別觀測(cè)分劃尺上不同的分劃間隔。（2）、視準(zhǔn)軸與水準(zhǔn)管軸相互關(guān)系的檢驗(yàn)校正：、i角檢驗(yàn)：設(shè)水準(zhǔn)管軸不平行于視準(zhǔn)軸，它們之間的交角為i。當(dāng)水準(zhǔn)管氣泡居中時(shí)，視準(zhǔn)軸不在水平線上而傾斜了i角，水準(zhǔn)儀至水準(zhǔn)尺的距離越遠(yuǎn)，由此引起的讀數(shù)偏差也越大。檢驗(yàn)時(shí)，在平坦地面的一直線上選定J1，A，B，J2四點(diǎn)，點(diǎn)間距均為20.6m， 具體檢驗(yàn)結(jié)果如下：i角誤差的檢驗(yàn)對(duì)于精密水準(zhǔn)儀，要求i角小于20。超出限差時(shí)，必須對(duì)其進(jìn)行校正。算出i=12<20,所以檢驗(yàn)合格。九、GPS靜態(tài)相對(duì)定位測(cè)量與計(jì)算成果：1、外業(yè)觀測(cè)：（1) 調(diào)度

9、安排，確定每臺(tái)接收機(jī)觀測(cè)的測(cè)站、開(kāi)機(jī)時(shí)間、搬站交通布置等。 （2) 觀測(cè)組按調(diào)度表規(guī)定的時(shí)間進(jìn)行作業(yè)，保證同步觀測(cè)同一衛(wèi)星組。（3) 每時(shí)段開(kāi)機(jī)前，作業(yè)員量取天線高，并及時(shí)記錄（或輸入）測(cè)站名、年月日、時(shí)段號(hào)、天線高等信息。關(guān)機(jī)后再量取一次天線高作校核，兩次量天線高互差不得大于3mm，取平均值作為最后結(jié)果。若互差超限，應(yīng)查明原因，提出處理意見(jiàn)記入測(cè)量手簿備注欄中。（4) 儀器工作過(guò)程后，作業(yè)員應(yīng)對(duì)照指示燈工作狀態(tài)說(shuō)明，判斷儀器工作是否正常。 （5) 一個(gè)時(shí)段觀測(cè)過(guò)程中不得進(jìn)行以下操作：關(guān)閉接收機(jī)又重新啟動(dòng)；進(jìn)行自測(cè)試(發(fā)現(xiàn)故障除外)；改變衛(wèi)星高度角；改變數(shù)據(jù)采樣間隔；改變天線位置；按動(dòng)關(guān)閉文件

10、和刪除文件等功能鍵。（6) 觀測(cè)員在作業(yè)期間不得擅自離開(kāi)測(cè)站，并應(yīng)防止儀器受震動(dòng)和被移動(dòng)，防止人和其他物體靠近天線，遮擋衛(wèi)星信號(hào)。（7) 接收機(jī)在觀測(cè)過(guò)程中不應(yīng)在接收機(jī)近旁使用對(duì)講機(jī)；雷雨過(guò)境時(shí)應(yīng)關(guān)機(jī)停測(cè)，并卸下天線以防雷擊。（8) 應(yīng)記錄雨、晴、陰、云等天氣狀況。2、內(nèi)業(yè)處理與平差計(jì)算：（1）數(shù)據(jù)傳輸；（2）數(shù)據(jù)預(yù)處理；（3）基線解算；（4）閉合差檢驗(yàn)；（5）平差（自由網(wǎng)平差和約束平差）；（6）網(wǎng)平差成果報(bào)告。3、成果（1）平面坐標(biāo)站點(diǎn)北向(x)/中誤差(m)東向(y)/中誤差(m)中誤差(m)誤差橢圓E(m)F(m)ET(度:分:秒)E0883787813.0590*513642.0340

11、*I0083787764.86480.0025513107.40210.00270.00370.00290.0023121:49:06I0093787637.20490.0047513006.96260.00510.00690.00530.0045121:55:46I0243787894.72700.0062512721.41030.00570.00850.00700.0047141:41:51I0253787992.32660.0070512587.96850.00610.00930.00750.0055148:05:22I0433787313.99260.0092512404.42800.

12、01150.01470.01150.009294:10:01I0453787466.54370.0102512392.40700.02070.02300.02220.0063112:07:59I0683786818.10160.0050512808.93240.00620.00790.00660.0044117:16:56I0743787778.56090.0090513570.51580.00740.01170.00900.0074176:07:29I0753787935.43510.0094513689.11150.00870.01280.00940.00878:08:46I6737866

13、57.22130.0053512678.78920.00610.00800.00670.0045124:14:35泉山點(diǎn)3788109.4670*514667.6677*2 基線向量及改正數(shù) 基線起點(diǎn)->終點(diǎn).時(shí)段DX/改正數(shù)(m)DY/改正數(shù)(m)DZ/改正數(shù)(m)距離/改正數(shù)(m)中誤差/相對(duì)誤差(m)I025->I075.191B-955.5505-434.5615-93.10521053.84490.03240.0088-0.00420.0001-0.00631: 32550I67->I068.191H-74.4072-140.6239132.3218206.9314

14、0.0134-0.00120.0009-0.0005-0.00051: 15429I025->泉山點(diǎn).191B-1788.7431-954.855959.34162028.51500.0171-0.00160.00190.00450.00061: 118946I67->I045.191H459.2072-270.4475673.0827858.51820.03860.0084-0.0091-0.00090.00671: 22243I025->E088.191B-945.3791-350.8361-193.50631026.77760.0174-0.00150.0012-0.0

15、0580.00211: 59125I67->I043.191G409.7618-199.8845546.6210711.79560.02650.0034-0.0129-0.00700.00021: 26888I025->I008.191B-482.2044-82.2523-232.0504541.41840.01730.0025-0.00660.0108-0.00581: 31253I67->I008.191G-97.0160-753.6495912.79751187.68640.01320.02810.0102-0.0257-0.02861: 90314I025->I

16、009.191B-425.472327.6369-337.3635543.69530.01890.0004-0.0032-0.01820.01081: 28816I67->I009.191G-40.2839-643.7602807.48441033.48000.01460.00230.00180.00280.00101: 70927I025->I024.191C-142.2979-14.6311-82.9673165.36730.0201-0.00010.00330.0036-0.00201: 8212I67->泉山點(diǎn).191H-1403.5547-1626.25301204

17、.18942462.66860.01260.0015-0.00720.00910.00841: 195312I025->I074.191B-890.5598-300.8946-221.8900965.85170.0256-0.0018-0.0014-0.00320.00291: 37732I67->E088.191H-560.1907-1022.2333951.34151504.60130.01240.0009-0.00420.00510.00571: 121175I075->泉山點(diǎn).191B-833.1926-520.2943152.4467994.06040.0296-0

18、.01240.01440.01780.00561: 33585I068->I045.191H533.6144-129.8236540.7609770.72740.0390-0.0008-0.00330.01040.00731: 19752I075->E088.191B10.171483.7254-100.4012131.12510.0297-0.01330.01080.00460.00231: 4417I068->I043.191G484.1689-59.2605414.2991639.98060.02680.0072-0.01080.00140.00731: 23838I0

19、75->I008.191B473.3461352.3092-138.9452606.20460.02960.0016-0.00330.0121-0.00351: 20497I068->I008.191G-22.6088-613.0255780.4756992.70020.01340.02070.0099-0.0240-0.02551: 74289I075->I009.191B530.0782462.1985-244.2583744.49480.0306-0.00280.0201-0.00660.01261: 24355I068->I009.191G34.1233-503

20、.1363675.1625842.70690.01500.00380.00260.00410.00191: 56141I075->I024.191C813.2526419.930410.1379915.32740.0313-0.00510.0021-0.0005-0.00361: 29233I068->泉山點(diǎn).191H-1329.1475-1485.62911071.86762263.32210.01280.0026-0.00740.00930.00771: 177442I075->I074.191B64.9907133.6669-128.7848196.66260.0348

21、0.0035-0.01050.0003-0.00621: 5646I068->E088.191H-485.7835-881.6093819.01971297.69560.01260.0030-0.00440.00670.00611: 103286泉山點(diǎn)->E088.191B843.3640604.0197-252.84791067.72410.0086-0.0016-0.0016-0.01170.00061: 124423I045->I043.191G-49.445570.5630-126.4617153.02470.0444-0.0041-0.0032-0.00510.00

22、401: 3450泉山點(diǎn)->I008.191B1306.5387872.6035-291.39191597.93290.00800.0042-0.00760.0070-0.00201: 199783泉山點(diǎn)->I009.191B1363.2708982.4928-396.70501726.60770.01270.00180.0019-0.01600.00621: 135882I045->I009.191G-499.4911-373.3127134.4017637.90090.0389-0.00650.01240.0063-0.00091: 16400泉山點(diǎn)->I024.1

23、91C1646.4452940.2248-142.30881901.33010.01600.0013-0.0010-0.00260.00091: 118969I045->泉山點(diǎn).191H-1862.7619-1355.8055531.10672364.35290.03850.0016-0.00490.00030.00161: 61463泉山點(diǎn)->I074.191B898.1833653.9613-281.23161146.07580.0226-0.0001-0.0029-0.00810.00021: 50678I045->E088.191H-1019.3979-751.785

24、8278.25881296.83530.03820.0024-0.00320.00680.00151: 33925E088->I008.191B463.1747268.5838-38.5440536.79950.00860.0054-0.00410.01880.00131: 62715I043->I008.191G-506.7778-553.7650366.1765835.20330.02620.0172-0.0065-0.0316-0.02001: 31935E088->I009.191B519.9068378.4731-143.8571658.96870.01230.00

25、560.0017-0.00790.00721: 53507I043->I009.191G-450.0456-443.8757260.8634683.82490.0265-0.0044-0.0152-0.00090.01251: 25821E088->I024.191C803.0812336.2050110.5391877.60590.01630.0027-0.00140.00510.00261: 53793I043->泉山點(diǎn).191H-1813.3164-1426.3685657.56842398.96640.0260-0.0050-0.00410.00720.00821:

26、92260E088->I074.191B54.819349.9415-28.383779.40370.0229-0.00190.00350.0079-0.00191: 3474I043->E088.191H-969.9524-822.3488404.72051334.49010.0257-0.00900.01170.01810.00481: 51831I008->I009.191B56.7321109.8893-105.3131162.43470.0128-0.00570.0133-0.02070.02051: 12668I008->I009.191G56.732110

27、9.8893-105.3131162.43470.0128-0.0150-0.01310.0328-0.03541: 12668I008->I024.191C339.906567.6213149.0831377.27290.0163-0.00720.0208-0.0039-0.00431: 23094I008->泉山點(diǎn).191H-1306.5387-872.6035291.39191597.93290.0080-0.0207-0.02150.03120.03431: 199783I008->I074.191B-408.3554-218.642310.1603463.31610

28、.02280.00190.0012-0.0186-0.00261: 20325I008->E088.191H-463.1747-268.583838.5440536.79950.0086-0.0273-0.01260.03330.03221: 62715I009->I024.191C283.1744-42.2680254.3962383.00360.0181-0.0051-0.00050.02100.01031: 21111I009->I074.191B-465.0875-328.5316115.4735581.01080.02420.0057-0.01020.00820.0

29、0281: 23991I009->E088.191H-519.9068-378.4731143.8571658.96870.0123-0.0011-0.00760.00190.00561: 53507I024->I074.191B-748.2619-286.2635-138.9228813.10650.0249-0.00160.00520.0015-0.00061: 32602泉山點(diǎn)->E088.191H843.3640604.0197-252.84791067.72410.00860.00020.0023-0.00280.00211: 1244233、高程擬合成果站點(diǎn)高程(

30、m)高程(m)E08839.9664I00841.14070.0346I00941.38220.0349I02446.51460.0359I02550.26550.0361I04351.59080.0385I04552.11880.0424I06844.90180.0349I07440.03500.0380I07538.62500.0414I6745.93510.0348泉山點(diǎn)55.9500十、水平角觀測(cè)、內(nèi)業(yè)計(jì)算與成果1.觀測(cè)方法。水平角用DJ2經(jīng)緯儀觀測(cè)2個(gè)測(cè)回（注意瞄準(zhǔn)覘板上標(biāo)志中心或地面導(dǎo)線點(diǎn)中心）。每測(cè)回度盤位置分別為：0°00，90°30。每測(cè)回各方向2C互差1

31、3。角度測(cè)回互差限差為9。2.導(dǎo)線邊改正后的斜距化算為高斯平面邊長(zhǎng)；測(cè)角中誤差4；測(cè)距中誤差12mm；計(jì)算導(dǎo)線方位角閉合差（限差為10n），合限并分配角度閉合差；計(jì)算導(dǎo)線坐標(biāo)閉合差（限差為1/15000），合限后轉(zhuǎn)入平差計(jì)算。提供各控制點(diǎn)的54坐標(biāo)。精度評(píng)定：=+20秒<； =0.069m；K=0.0667/1555.9987=1/23328<3.水平角原始數(shù)據(jù)站點(diǎn)覘點(diǎn)讀數(shù)左-右(2C)左減右(2C) °半側(cè)回方向 °盤左°盤右°方向值I24I250 00 00179 59 53-7268 10 47 268 10 50I23 268 10

32、4788 10 45+2268 10 52I24I2590 30 01270 30 05-4268 10 41268 10 44I23358 40 42178 40 52-10268 10 47I17I180 00 00179 59 50+10268 36 36268 36 41I16268 36 3688 36 35+1268 36 45I09I1890 30 00270 29 47+13268 36 40268 36 37 I64359 06 40179 06 2020268 36 334.水平角觀測(cè)成果導(dǎo)線計(jì)算代碼：clear;%清空內(nèi)存yzds=2;wzds=15;%yzds是已知點(diǎn)數(shù)

33、，wzds是未知點(diǎn)數(shù)gcjs=17;gcbs=17;cjzwc=1.5;cbzwc=2.0;cbblwc=0.0;%觀測(cè)角數(shù)=12；觀測(cè)邊數(shù)=6；測(cè)角中誤差=1.5秒；（測(cè)邊固定誤差、比例誤差，以cm為單位）JKH=-1,-2,1;1,-1,2;2,1,3;3,2,4;4,3,5;5,4,6;6,5,7;7,6,8;8,7,9;9,8,10;10,9,11;11,10,12;12,11,13;13,12,14;14,13,15;15,14,-1;-1,15,1;%觀測(cè)角的端點(diǎn)編號(hào),順序是測(cè)站，后視，前視。JK=-2,-1;-1,1;1,2;2,3;3,4;4,5;5,6;6,7;7,8;8,9

34、;9,10;10,11;11,12;12,13;13,14;14,15;15,-1;%觀測(cè)邊的端點(diǎn)編號(hào)YZDZB=3787894.7270,512721.4103;3787992.3266,512587.9685;%已知點(diǎn)坐標(biāo)賦值JSWZD=-1,-2,1;1,-1,2;2,1,3;3,2,4;4,3,5;5,4,6;6,5,7;7,6,8;8,7,9;9,8,10;10,9,11;11,10,12;12,11,13;13,12,14;14,13,15;15,14,-1;-1,15,1;%未知點(diǎn)近似坐標(biāo)推算路線，即測(cè)站、后視、前視；前視點(diǎn)就是需要推算的未知點(diǎn)，此時(shí)測(cè)站和后視點(diǎn)坐標(biāo)已知或已經(jīng)推出

35、。gcjd=268.17972,190.37583,205.09,107.91889,91.81528,258.61972,100.19028,171.92222,141.25639,179.33583,170.96639,100.06417,174.96361,268.61,110.55972,108.42306,139.89083;%觀測(cè)角度gcbc=165.3265,191.6509,113.7752,62.6204,79.7123,174.0970,102.8376,120.5101,54.2983,66.2271,97.5891,36.6804,119.3408,171.8028,7

36、2.4477,50.5438,41.8650;%觀測(cè)邊長(zhǎng)ro=2062.65;pi=3.141592654;%角度化弧度的常數(shù)for i=1:yzds; XYJS(i+wzds,1)=YZDZB(i,1);XYJS(i+wzds,2)=YZDZB(i,2);%將已知點(diǎn)坐標(biāo)值賦給另外一個(gè)數(shù)組，此數(shù)組的前面將存貯未知點(diǎn)坐標(biāo).endfor i=1:wzds;%開(kāi)始計(jì)算未知點(diǎn)坐標(biāo) j=JSWZD(i,1);k=JSWZD(i,2);h=JSWZD(i,3);%從未知點(diǎn)推算路線數(shù)組中讀取推算信息 for jj=1:gcjs; if j=JKH(jj,1) & k=JKH(jj,2) &

37、h=JKH(jj,3) jds=jj; if j<=0 j=abs(j);cjjsjx=YZDZB(j,1);cjjsjy=YZDZB(j,2); else cjjsjx=XYJS(j,1);cjjsjy=XYJS(j,2); end if k<=0 k=abs(k);cjjskx=YZDZB(k,1);cjjsky=YZDZB(k,2); else cjjskx=XYJS(k,1);cjjsky=XYJS(k,2); end dxjk=cjjskx-cjjsjx;dyjk=cjjsky-cjjsjy;%以上為檢索出J、k點(diǎn)的端點(diǎn)坐標(biāo)if abs(dxjk) >= abs(d

38、yjk) afjkt=atan(dyjk/dxjk); if dxjk>=0 & dyjk>=0 afjk=afjkt; end if dxjk<0 afjk=pi+afjkt; end if dxjk>=0 & dyjk<0 afjk=2*pi+afjkt; end end if abs(dxjk) < abs(dyjk) afjkt=atan(dxjk/dyjk); if dyjk>=0 afjk=0.5*pi-afjkt; end if dyjk<0 afjk=1.5*pi-afjkt; end end c=gcjd(1,j

39、ds); chd=(fix(c)+fix(c-fix(c)*100)/60+(c-fix(c)-fix(c-fix(c)*100)/100)*100/36)*pi/180; afjh0=afjk+chd; end j=JSWZD(i,1);h=JSWZD(i,3); for hh=1:gcbs;%此循環(huán)是檢索j、h對(duì)應(yīng)的邊長(zhǎng)序號(hào)，以便調(diào)用該邊長(zhǎng)觀測(cè)值 if j=JK(hh,1) & h=JK(hh,2) bcs=hh; end if j=JK(hh,2) & h=JK(hh,1) bcs=hh; end end end XYJS(h,1)=cjjsjx+gcbc(1,bcs)*

40、cos(afjh0);XYJS(h,2)=cjjsjy+gcbc(1,bcs)*sin(afjh0);%計(jì)算未知點(diǎn)近似坐標(biāo)B=zeros(gcjs+gcbs,wzds*2); % 給誤差方程系數(shù)陣的元素賦0,行數(shù)是觀測(cè)值個(gè)數(shù),列數(shù)是未知數(shù)個(gè)數(shù).i=1; while 1 if i=gcjs+1;%觀測(cè)角的個(gè)數(shù)是gcjs,因此i=gcjs+1時(shí)中斷 break; endj=JKH(i,1);k=JKH(i,2);h=JKH(i,3);%讀取角度信息 if j<=0 j=wzds+abs(j);%未知點(diǎn)個(gè)數(shù)是wzds個(gè) end if k<=0 k=wzds+abs(k);%未知點(diǎn)個(gè)數(shù)是w

41、zds個(gè) end if h<=0 h=wzds+abs(h);%未知點(diǎn)個(gè)數(shù)是wzds個(gè) end dxjk=XYJS(k,1)-XYJS(j,1);dyjk=XYJS(k,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; if abs(dxjk) >= abs(dyjk) afjkt=atan(dyjk/dxjk); if dxjk>=0 & dyjk>=0 afjk=afjkt;

42、end if dxjk<0 afjk=pi+afjkt; end if dxjk>=0 & dyjk<0 afjk=2*pi+afjkt; end end if abs(dxjk) < abs(dyjk) afjkt=atan(dxjk/dyjk); if dyjk>=0 afjk=0.5*pi-afjkt; end if dyjk<0 afjk=1.5*pi-afjkt; end endif abs(dxjh) >= abs(dyjh) afjht=atan(dyjh/dxjh); if dxjh>=0 & dyjh>=0

43、 afjh=afjht; end if dxjh<0 afjh=pi+afjht; end if dxjh>=0 & dyjh<0 afjh=2*pi+afjht; end end if abs(dxjh) < abs(dyjh) afjht=atan(dxjh/dyjh); if dyjh>=0 afjh=0.5*pi-afjht; end if dyjh<0 afjh=1.5*pi-afjht; end endc=gcjd(1,i); chd=(fix(c)+fix(c-fix(c)*100)/60+(c-fix(c)-fix(c-fix(c)*

44、100)/100)*100/36)*pi/180; fwjhc=afjh-afjk; if fwjhc<0 fwjhc=2*pi+fwjhc; end l(i)=(chd-fwjhc)*206265;%此句為計(jì)算角度觀測(cè)值的自由項(xiàng)，單位是秒j=JKH(i,1);k=JKH(i,2);h=JKH(i,3);if j<=0 & k<=0A=0; end if j<=0 & k>=0 B(i,k*2-1)=+ro*dyjk/sjk0;B(i,k*2)=-ro*dxjk/sjk0; end if k<=0 & j>=0 B(i,j*2-

45、1)=-ro*dyjk/sjk0;B(i,j*2)=+ro*dxjk/sjk0; end if k>=0 & 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 & h>=0 B(i,h*2-1)=-ro*dyjh/sjh0;B(i,h*2)=+ro*dxjh/sjh0; end if j>=0 &

46、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; end if j>=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)=+ro*dxjh/sjh0; end i=i+1; end while 2 if i=gcjs+gcbs+1;%觀測(cè)角加上觀測(cè)邊的總數(shù)是gcjs+gcbs,因此i=gcjs+gcbs+1時(shí)中斷

47、 break; endj=JK(i-gcjs,1);k=JK(i-gcjs,2);%減去觀測(cè)邊的個(gè)數(shù),gcjs if j<=0 j=wzds+abs(j);%未知點(diǎn)個(gè)數(shù)是wzds個(gè) end if k<=0 k=wzds+abs(k);%未知點(diǎn)個(gè)數(shù)是wzds個(gè) end dxjk=XYJS(k,1)-XYJS(j,1);dyjk=XYJS(k,2)-XYJS(j,2); sjk0=sqrt(dxjk*dxjk+dyjk*dyjk); l(i)=(gcbc(1,i-gcjs)-sjk0)*100;%此句為計(jì)算邊長(zhǎng)觀測(cè)值的自由項(xiàng)，單位是厘米。 j=JK(i-gcjs,1);k=JK(i-g

48、cjs,2);%減去觀測(cè)邊的個(gè)數(shù),gcjs if j<=0 & k<=0 A=0; end if j<=0 & k>=0 B(i,k*2-1)=dxjk/sjk0;B(i,k*2)=dyjk/sjk0; end if k<=0 & j>=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; e