【正文】 ，附合線路所包含的基線邊數(shù)，也不能超過一定的數(shù)量。它廣泛用于工程放樣、邊界測(cè)量、地籍測(cè)量和碎部測(cè)量等。 環(huán)形網(wǎng) 星形網(wǎng) GNSS網(wǎng)的布設(shè)原則網(wǎng)的圖形設(shè)計(jì)雖然主要決定于用戶的要求，但是經(jīng)費(fèi)、時(shí)間和人為的消耗以及所需接收設(shè)備的類型、數(shù)量和后勤保障條件等，也都與網(wǎng)的設(shè)計(jì)有關(guān)。2．至少有3個(gè)GPS控制網(wǎng)點(diǎn)與地面控制網(wǎng)點(diǎn)重合，才能求定GNSS網(wǎng)坐標(biāo)與原有地面控制網(wǎng)坐標(biāo)之間的坐標(biāo)轉(zhuǎn)換參數(shù)。并要遠(yuǎn)離變電所，高壓線及微波輻射干擾源。因此，各個(gè)點(diǎn)之間不需要通視的。因此，經(jīng)典控制網(wǎng)的一類優(yōu)化設(shè)計(jì)在GPS網(wǎng)中稱為網(wǎng)形的結(jié)構(gòu)設(shè)計(jì)。經(jīng)典控制網(wǎng)中，誤差具有累積性，網(wǎng)中各邊的相對(duì)精度和方位精度不均勻。控制網(wǎng)設(shè)計(jì)的目標(biāo)，指的是控制網(wǎng)應(yīng)達(dá)到的質(zhì)量標(biāo)準(zhǔn)，它是設(shè)計(jì)的依據(jù)和目的，同時(shí)又是評(píng)定網(wǎng)質(zhì)量的指標(biāo)。設(shè)坐標(biāo)未知參數(shù)的方差—協(xié)方差陣 （） 則作為整體精度標(biāo)準(zhǔn)的指標(biāo)有：A最優(yōu)：若 （）N最優(yōu)：即DXX的范數(shù)║DXX║滿足： ║DXX║=min （）D最優(yōu)：若 （）成立，則稱為D最優(yōu)。未知數(shù)某些函數(shù)的精度：比如控制網(wǎng)中推算邊長(zhǎng)、方位角的精度等，設(shè)有： （） 則有： （） （2）可靠性標(biāo)準(zhǔn)可靠性理論的前提是考慮觀測(cè)值中含有隨機(jī)誤差和粗差，并把粗差歸入函數(shù)模型之中來評(píng)價(jià)網(wǎng)的質(zhì)量。定義 （） 為第個(gè)觀測(cè)值的多余觀測(cè)分量，且 （為多余觀測(cè)數(shù)） （）內(nèi)部可靠性指標(biāo)：在顯著水平下，以檢驗(yàn)功效發(fā)現(xiàn)粗差的下界為 （） （）外部可靠性指標(biāo)：表示發(fā)現(xiàn)不了的粗差對(duì)平差結(jié)果的影響的指標(biāo)。所以多余觀測(cè)分量與可靠性標(biāo)準(zhǔn)有直接聯(lián)系，若要求可靠性指標(biāo)在一定的范圍內(nèi)，就相當(dāng)于對(duì)多余觀測(cè)分量和總的多余觀測(cè)提出制約。網(wǎng)的優(yōu)化設(shè)計(jì)，就是得出在費(fèi)用使用最少的情況下，使其精度指標(biāo)和可靠性指標(biāo)能滿足布網(wǎng)方案的要求。網(wǎng)的優(yōu)化設(shè)計(jì)中，觀測(cè)費(fèi)用為主要考慮的對(duì)象。可分為零、一、二、三類優(yōu)化設(shè)計(jì)。主要采用自由網(wǎng)平差和進(jìn)行s變形，得到位置、定向和尺度參數(shù)等的一組基準(zhǔn)數(shù)值。3．二類設(shè)計(jì)（權(quán)設(shè)計(jì)）。固定參數(shù)是和部分，待定參數(shù)為部分和，它是在舊網(wǎng)的基礎(chǔ)上，對(duì)現(xiàn)有網(wǎng)和現(xiàn)有設(shè)計(jì)進(jìn)行改進(jìn)，通過設(shè)計(jì)新點(diǎn)位、引入附加觀測(cè)值和觀測(cè)值權(quán)的改變，導(dǎo)致點(diǎn)位增刪或移動(dòng)，觀測(cè)值的增刪或精度改變，已達(dá)到全網(wǎng)點(diǎn)位精度的要求。（2）模擬法模擬法是根據(jù)經(jīng)驗(yàn)設(shè)計(jì)的初步控制網(wǎng)的網(wǎng)形和設(shè)計(jì)精度要求的準(zhǔn)則，模擬一組起始數(shù)據(jù)與觀測(cè)值（在設(shè)計(jì)圖上量?。┹斎胗?jì)算機(jī)，按參數(shù)平差法，組成誤差方程式、法方程，并對(duì)法方程系數(shù)矩陣求逆，從而得到未知參數(shù)的協(xié)因數(shù)陣以及方差陣，計(jì)算未知參數(shù)及其函數(shù)的精度，估算成本，或進(jìn)一步計(jì)算可靠性數(shù)值等信息。模擬法的優(yōu)點(diǎn)是計(jì)算簡(jiǎn)單，設(shè)計(jì)程序易于編制，優(yōu)化過程可據(jù)作業(yè)員經(jīng)驗(yàn)進(jìn)行人工干預(yù)。對(duì)于一項(xiàng)GNSS測(cè)量工程項(xiàng)目，一般有如下一些要求。4．提交成果的內(nèi)容用戶需要提交哪些成果，除了提交最終的結(jié)果外，是否還需要提交原始數(shù)據(jù)或中間數(shù)據(jù)等。這些資料包括測(cè)區(qū)的各種比例尺地形圖，交通圖，氣象資料；已有的大地測(cè)量成果資料，例如點(diǎn)之記、成果表、技術(shù)總結(jié)、三角網(wǎng)略圖及水準(zhǔn)路線圖。GPS工程控制網(wǎng)優(yōu)化設(shè)計(jì)，是實(shí)施其控制測(cè)量工作的第一步，是一項(xiàng)基礎(chǔ)性的工作。控制網(wǎng)的布設(shè)在精度估算方面是必不可少的一個(gè)環(huán)節(jié)。網(wǎng)的三維平差一般是在三維空間大地坐標(biāo)系或三維直角大地坐標(biāo)系或三維空間大地坐標(biāo)系下進(jìn)行的。（2）約束平差網(wǎng)的約束平差中所采用的觀測(cè)值是GNSS基線向量，但與無約束平差不同的是，在平差過程中引入了會(huì)使網(wǎng)的尺度和方位產(chǎn)生由非基線向量所引起的外部起算數(shù)據(jù)。 平面控制網(wǎng)精度評(píng)估的理論依據(jù)——參數(shù)平差為了表示直觀，繪出邊長(zhǎng)、。 平面控制網(wǎng)的精度評(píng)定1．未知點(diǎn)坐標(biāo)的協(xié)方差因數(shù)陣平面控制網(wǎng)以未知點(diǎn)坐標(biāo)，為參數(shù)的參數(shù)平差的優(yōu)點(diǎn)之一，就是法方程系數(shù)陣的逆矩陣即為未知點(diǎn)坐標(biāo)的協(xié)方差因數(shù)陣，即 （）式中，N為法方程系數(shù)矩陣。2．參數(shù)中誤差第號(hào)未知點(diǎn)縱橫坐標(biāo)方位的中誤差為 （）點(diǎn)位中誤差為 （）邊長(zhǎng)中誤差為 （）邊長(zhǎng)是兩未知點(diǎn)，間的長(zhǎng)度。研究了利用平面控制網(wǎng)精度估算的原理，采用定權(quán)的方式，對(duì)控制網(wǎng)進(jìn)行整體平差，來估算平面控制網(wǎng)的精度。自1984年推向市場(chǎng)以來，經(jīng)過不斷的完善和發(fā)展，現(xiàn)已成為國(guó)際上科技與工程應(yīng)用領(lǐng)域最具有影響力的科學(xué)計(jì)算軟件。在MATLAB中，無論是問題還是結(jié)果的表達(dá)都采用我們熟悉的數(shù)學(xué)表示方法，而不需要用傳統(tǒng)的編程語言進(jìn)行前后處理。2．直觀靈活的語言MATLAB不僅是一套打包好的函數(shù)庫(kù)，同時(shí)也是一種高級(jí)的、面向?qū)ο蟮木幊陶Z言?？梢暬ぞ甙ㄇ驿秩?、線框圖、偽彩圖、光源、三維等高圖、圖像顯示、動(dòng)畫、體積可視化等。只需調(diào)用不同的繪圖函數(shù)，即可在圖上標(biāo)出圖題、標(biāo)注XY軸、繪制格，簡(jiǎn)單易行。 GNSS網(wǎng)精度估算程序設(shè)計(jì)一側(cè)邊控制網(wǎng)，網(wǎng)中P1，P2為未知點(diǎn)，P3，P4，P5，P6為已知點(diǎn)，選用某測(cè)距儀觀測(cè)8條邊[23]。本章根據(jù)已有GNSS網(wǎng)圖件、已知起算數(shù)據(jù)和觀測(cè)邊長(zhǎng)，利用參數(shù)平差的方法列出誤差方程和系數(shù)陣的逆陣，然后再求得參數(shù)函數(shù)的中誤差進(jìn)行平面控制網(wǎng)的精度估算。GNSS作業(yè)精度高、受環(huán)境的和距離影響限制小。GPS控制網(wǎng)優(yōu)化設(shè)計(jì)的目的，是通過對(duì)不同方案布網(wǎng)的精度、可靠性、經(jīng)濟(jì)指標(biāo)等方面綜合比較，挑選出技術(shù)上可行，勞動(dòng)力較低、工作量較少、作業(yè)效率較高的最優(yōu)方案。隨著觀測(cè)時(shí)段的增長(zhǎng)，GPS工程控制網(wǎng)中的點(diǎn)位精度、邊長(zhǎng)精度和方位精度都在逐漸的提高。參考文獻(xiàn)參考文獻(xiàn)在文中用右上角標(biāo)進(jìn)行標(biāo)注[1] 孔祥元,郭際明,[M]參考專著的標(biāo)準(zhǔn)引用格式.武漢:武漢大學(xué)出版社,2003:222223.[2] [J]參考期刊論文的引用格式.現(xiàn)代測(cè)繪,2003,26(5):68.[3] 劉立龍,[J].桂林工學(xué)院學(xué)報(bào),2002,(3): 225227.[4] [M].鄭州:黃河水利出版社,2005:46.[5] 黃子發(fā),熊永良,[M].北京:科學(xué)出版社,2009: 239251.[6] 程鵬飛,蔡艷輝,[M].北京:測(cè)繪出版社,2009:257304.[7] 李征航,[M].武漢:武漢大學(xué)出版社,2005:134135.[8] 張勤,[M].北京:科學(xué)出版社,2006:144147.[9] 孔祥元,郭際明,[M].武漢:武漢大學(xué)出版社,2011:222226.[10] 鮑建寬,李秀海,朱繼文,[M].哈爾濱:哈爾濱地圖出版社,2008:3236.[11] 徐愛功,韓曉東,[M].徐州:中國(guó)礦業(yè)大學(xué)出版社,2009:134138,158.[12] 李國(guó)波,[J].湘潭師范學(xué)院學(xué)報(bào)（自然科學(xué)版）,2004,26(3):1214.[13] [M].北京:測(cè)繪出版社,1992.[14] 李青岳,[M].北京:測(cè)繪出版社, 2008.[15] 吳學(xué)偉,[M].北京:科學(xué)出版社,2006:129130.[16] [M].武漢大學(xué)出版社,2009:111125.[17] [J].工程勘察,1992,(2):5862.[18] [M].北京:測(cè)繪出版社,2009:205212.[19] 蔡旭輝,劉衛(wèi)國(guó),[M].北京:人民郵電出版社,2009:12.[20] 孫曉,王勇,[J].焦作工學(xué)院學(xué)報(bào)(自然科學(xué)版),2002,21(5):393395.[21] [M].北京:清華大學(xué)出版社,2011:14.[22] 褚洪生,杜增吉,[M].北京:機(jī)械工業(yè)出版社,2007:14.[23] 隋立芬,[M].北京:測(cè)繪出版社,2010:7779.[24] ，：The .[25] Wolfgang Torge. Geodesy. Second Edition. Berlin. New York, 1991.[26] Baarda, Concepts in Commision ,New Series , , Delft,1967.致 謝在這篇論文完成之際，我衷心地感謝我的指導(dǎo)老師，李秀海老師在我大學(xué)的最后學(xué)習(xí)階段——畢業(yè)設(shè)計(jì)階段給自己的指導(dǎo)，李老師對(duì)我們認(rèn)真負(fù)責(zé)，及時(shí)傳遞相關(guān)信息，并在我們做課題期間給予了極大的幫助和支持。我還要感謝我的各位同學(xué)，很榮幸和他們成為同學(xué)和朋友，大家在一起互相幫助，互相鼓勵(lì)，共同進(jìn)步，這段美好的學(xué)習(xí)和生活經(jīng)歷將成為我一生最珍貴的回憶。當(dāng)觀測(cè)一次時(shí)的權(quán)為39。 DY(j,k)=Y(k)Y(j)。k=input(39。j=39。) a=DX(j,k)/S(j,k)。)。b=DY(j,k)/S(j,k)。k=input(39。 B4=[a b 0 0]endfor i=5 j=input(39。k=39。j=39。) a=DX(j,k)/S(j,k)。)。b=DY(j,k)/S(j,k)。k=input(39。 B8=[0 0 a b]enddisp(39。B3。B7。*P*Bdisp(39。)disp(39。)j=input(39。k=39。d=r*a/S(j,k)F1=[c a d b]disp(39。)。b=DY(j,k)/S(j,k)。)j=input(39。k=39。d=r*a/S(j,k)F3=[c a d b]disp(39。*Nbb2*F2Q24=F339。 DY(j,k)=Y(k)Y(j)。j=39。) a=DX(j,k)/S(j,k)。k=input(39。 B2=[a b 0 0]endfor i=3 j=input(39。k=39。j=39。) a=DX(j,k)/S(j,k)。)。b=DY(j,k)/S(j,k)。k=input(39。 B6=[0 0 a b]endfor i=7 j=input(39。k=39。j=39。) a=DX(j,k)/S(j,k)。)B=[B1。B5。法方程系數(shù)矩陣39。)A=pinv(Nbb)disp(39。)Nbb1=[ ]Nbb2=[ ]r=206265disp(39。)。b=DY(j,k)/S(j,k)。)j=input(39。k=39。d=r*a/S(j,k)F2=[c a d b]disp(39。)。b=DY(j,k)/S(j,k)。)Q14=F139