【正文】 (j,k)=sqrt((X(j)X(k))^2+(Y(j)Y(k))^2) endendfor j=2 for k=3:6 DX(j,k)=X(k)X(j)。)。b=DY(j,k)/S(j,k) B1=[a b a b]endfor i=2 j=input(39。k=39。j=39。) a=DX(j,k)/S(j,k)。)。b=DY(j,k)/S(j,k)。k=input(39。 B5=[0 0 a b]endfor i=6 j=input(39。k=39。j=39。) a=DX(j,k)/S(j,k)。)。b=DY(j,k)/S(j,k)。B2。B6。)Nbb=B39。精度估算39。P1P4邊權(quán)函數(shù)式的矩陣形式39。k=input(39。c=r*b/S(j,k)。j=39。)a=DX(j,k)/S(j,k)。P2P4邊權(quán)函數(shù)式的矩陣形式39。k=input(39。c=r*b/S(j,k)。*Nbb1*F1Q23=F239。)P=[ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ]l=[ ]for j=1 for k=2:6 DX(j,k)=X(k)X(j)。 S(j,k)=sqrt((X(j)X(k))^2+(Y(j)Y(k))^2) endendfor i=1 j=input(39。k=39。)。b=DY(j,k)/S(j,k)。k=input(39。 B3=[a b 0 0]endfor i=4 j=input(39。k=39。j=39。) a=DX(j,k)/S(j,k)。)。b=DY(j,k)/S(j,k)。k=input(39。 B7=[0 0 a b]endfor i=8 j=input(39。k=39。誤差方程系數(shù)39。B4。B8]disp(39。法方程逆矩陣39。計(jì)算邊長(zhǎng)和方位精度39。j=39。)a=DX(j,k)/S(j,k)。P2P3邊權(quán)函數(shù)式的矩陣形式39。k=input(39。c=r*b/S(j,k)。j=39。)a=DX(j,k)/S(j,k)。計(jì)算協(xié)因數(shù)39。*Nbb2*F340。*Nbb1*F1Q23=F239。c=r*b/S(j,k)。k=input(39。P2P4邊權(quán)函數(shù)式的矩陣形式39。)a=DX(j,k)/S(j,k)。j=39。c=r*b/S(j,k)。k=input(39。P1P4邊權(quán)函數(shù)式的矩陣形式39。精度估算39。)Nbb=B39。B6。B2。b=DY(j,k)/S(j,k)。)。) a=DX(j,k)/S(j,k)。j=39。k=39。 B5=[0 0 a b]endfor i=6 j=input(39。k=input(39。b=DY(j,k)/S(j,k)。)。) a=DX(j,k)/S(j,k)。j=39。k=39。b=DY(j,k)/S(j,k) B1=[a b a b]endfor i=2 j=input(39。)。 S(j,k)=sqrt((X(j)X(k))^2+(Y(j)Y(k))^2) endendfor j=2 for k=3:6 DX(j,k)=X(k)X(j)。*Nbb2*F3附錄B當(dāng)進(jìn)行兩次觀測(cè)時(shí)的程序：X=[ ]Y=[ ]disp(39。計(jì)算協(xié)因數(shù)39。)a=DX(j,k)/S(j,k)。j=39。c=r*b/S(j,k)。k=input(39。P2P3邊權(quán)函數(shù)式的矩陣形式39。)a=DX(j,k)/S(j,k)。j=39。計(jì)算邊長(zhǎng)和方位精度39。法方程逆矩陣39。B8]disp(39。B4。誤差方程系數(shù)39。k=39。 B7=[0 0 a b]endfor i=8 j=input(39。k=input(39。b=DY(j,k)/S(j,k)。)。) a=DX(j,k)/S(j,k)。j=39。k=39。 B3=[a b 0 0]endfor i=4 j=input(39。k=input(39。b=DY(j,k)/S(j,k)。)。k=39。 S(j,k)=sqrt((X(j)X(k))^2+(Y(j)Y(k))^2) endendfor i=1 j=input(39。)P=[ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ]l=[ ]for j=1 for k=2:6 DX(j,k)=X(k)X(j)。在此向他們表示我由衷的謝意，并祝所有的老師培養(yǎng)出越來(lái)越多的優(yōu)秀人才，桃李滿天下，祝愿各位同學(xué)前程似錦！感謝在寫(xiě)論文階段給予我莫大幫助的同學(xué)。從最初的選題，到收集資料，到課題研究、編程序、撰寫(xiě)、修改，到論文最終定稿，他給了我細(xì)心的指導(dǎo)和幫助。5．本文進(jìn)一步還需要做的就是多進(jìn)行幾次觀測(cè)，看觀測(cè)到幾次的時(shí)候?qū)鹊挠绊憰?huì)很小，考慮一些外界條件和天氣對(duì)精度的影響，并將程序進(jìn)一步完善。精度評(píng)估作為一個(gè)重要的技術(shù)指標(biāo)，是控制網(wǎng)優(yōu)化設(shè)計(jì)的必不可少的一個(gè)環(huán)節(jié)。在收集資料的時(shí)候，我也收集了很多以前學(xué)者對(duì)GPS控制網(wǎng)精度估算的方法，比如通過(guò)VC++編寫(xiě)程序?qū)刂凭W(wǎng)的精度進(jìn)行估算，南方平差易軟件對(duì)精度的估算。采用不同次數(shù)觀測(cè)時(shí)段的方法，計(jì)算出各個(gè)點(diǎn)位精度、邊長(zhǎng)精度和方位精度，并對(duì)估算結(jié)果做出了對(duì)比和分析得出結(jié)論。1．進(jìn)行一次觀測(cè)時(shí)段的精度估算程序編制（1）按如下公式求得待定點(diǎn)近似坐標(biāo)： （） （）， （）， （）（2）按如下公式求得、并編程得到運(yùn)行結(jié)果如圖： （） （） （） 的運(yùn)行結(jié)果 的運(yùn)行結(jié)果 的運(yùn)行結(jié)果（3）計(jì)算求得系數(shù)矩陣B： 系數(shù)矩陣B的運(yùn)行結(jié)果（4）： 法方程系數(shù)矩陣和法方程逆矩陣的運(yùn)行結(jié)果（5）精度估算? 計(jì)算坐標(biāo)點(diǎn)位中誤差： ? ： 邊長(zhǎng)中誤差和方位中誤差的運(yùn)行結(jié)果2．進(jìn)行兩個(gè)時(shí)段觀測(cè)的精度估算程序編制 法方程系數(shù)矩陣和法方程逆矩陣的運(yùn)行結(jié)果（1）計(jì)算坐標(biāo)點(diǎn)位中誤差： （2）計(jì)算邊長(zhǎng)中誤差與方位中誤差： 邊長(zhǎng)中誤差與方位中誤差的運(yùn)行結(jié)果 本章小結(jié)本章主要介紹了MATLAB語(yǔ)言的發(fā)展、特點(diǎn)以及用MATLAB對(duì)平面控制網(wǎng)精度進(jìn)行估算。6．眾多面向應(yīng)用的工具箱和模塊集MATLAB的工具箱加強(qiáng)了對(duì)工程及科學(xué)中特殊應(yīng)用的支持。4．可擴(kuò)展性、開(kāi)放性強(qiáng)M文件是可見(jiàn)的MATLAB程序，所以我們可以查看源代碼。使用MATLAB可事半功倍地開(kāi)發(fā)自己的程序。這一點(diǎn)是MATLAB成為了數(shù)學(xué)分析、算法開(kāi)發(fā)及應(yīng)用程序開(kāi)發(fā)的良好環(huán)境。隨著MATLAB版本的提高，MATLAB的特點(diǎn)更為突出，數(shù)據(jù)類型更加豐富，人機(jī)交互界面更加方便，功能也更加強(qiáng)大。第5章 基于MATLAB的GNSS控制網(wǎng)精度估算程序設(shè)計(jì)與實(shí)現(xiàn)根據(jù)以上原理和思路，我將采用MATLAB軟件編寫(xiě)程序，并用算例進(jìn)行了驗(yàn)證。公式為[17~18] （） 本章小結(jié)本章介紹了平面控制網(wǎng)精度估算的原理，控制網(wǎng)的精度估算與控制網(wǎng)的優(yōu)化設(shè)計(jì)一樣，都是控制測(cè)量中較重要的一項(xiàng)工作。設(shè)網(wǎng)中未知點(diǎn)總數(shù)為，未知數(shù)的編序?yàn)?，則的形式為 （）當(dāng)然，由于N是對(duì)稱矩陣，也是對(duì)稱矩陣。 邊長(zhǎng)觀測(cè)示意圖 角度觀測(cè)示意圖1．，其中為邊的邊長(zhǎng)觀測(cè)值、二點(diǎn)坐標(biāo)值為參數(shù)，可以列出邊長(zhǎng)誤差方程式 （） 式中，若兩點(diǎn)中有一點(diǎn)為已知點(diǎn)或兩點(diǎn)都為已知點(diǎn)時(shí)，則誤差方程的系數(shù)為“0”，這是不用再進(jìn)行求解的，因此，在編輯設(shè)計(jì)數(shù)據(jù)時(shí)可以采用數(shù)字或字母加以區(qū)別，以節(jié)省計(jì)算機(jī)編程和運(yùn)行的時(shí)間。（3）聯(lián)合平差網(wǎng)的聯(lián)合平差是一個(gè)嚴(yán)格而有效的平差方法，所采用的觀測(cè)值不僅包括GNSS觀測(cè)值，除此以外還采用了含邊長(zhǎng)、角度、方向和高差等地面常規(guī)觀測(cè)值[15]。（2）二維平差所謂二維平差是指在平面二維坐標(biāo)系下進(jìn)行平差，所采用的GNSS觀測(cè)值為在一個(gè)二維坐標(biāo)下的GNSS基線向量觀測(cè)值，所確定出的結(jié)果為點(diǎn)的二維平面坐標(biāo)。第4章 平面控制網(wǎng)精度估算的原理平面控制網(wǎng)形一旦設(shè)計(jì)好后，按照已定的設(shè)計(jì)路線施工，在沒(méi)有測(cè)量以外誤差的影響下，控制網(wǎng)的點(diǎn)位精度是可估計(jì)的。全面系統(tǒng)的研究了GPS工程控制網(wǎng)的布設(shè)，探討工程控制網(wǎng)的設(shè)計(jì)方法。對(duì)收集的資料加以分析和研究，選取可靠和有價(jià)值的部分作為設(shè)計(jì)時(shí)參考。5．時(shí)限要求提交成果的有效期。1．測(cè)區(qū)概況測(cè)區(qū)的概況包括測(cè)區(qū)的地理位置、測(cè)區(qū)范圍、控制網(wǎng)的控制面積有多大。用計(jì)算機(jī)進(jìn)行設(shè)計(jì)計(jì)算的過(guò)程中，人機(jī)對(duì)話，使設(shè)計(jì)過(guò)程高效、靈活。與預(yù)定的精度要求、成本和可靠性要求等相比較，根據(jù)計(jì)算所提供的信息及設(shè)計(jì)者的經(jīng)驗(yàn)，對(duì)控制網(wǎng)的基準(zhǔn)、網(wǎng)形、觀測(cè)精度等進(jìn)行修正。 控制網(wǎng)設(shè)計(jì)階段的劃分情況表分 類固定參數(shù)待定參數(shù)含 義零類設(shè)計(jì)（ZOD）,基準(zhǔn)設(shè)計(jì)一類設(shè)計(jì)（FOD）圖形設(shè)計(jì)二類設(shè)計(jì)（SOD）觀測(cè)精度設(shè)計(jì)三類設(shè)計(jì)（THOD），部分已有網(wǎng)的改進(jìn) 優(yōu)化設(shè)計(jì)的方法 控制網(wǎng)優(yōu)化設(shè)計(jì)的方法有：解析法和模擬法。固定參數(shù)是和，待定參數(shù)為。2．一類設(shè)計(jì)（圖形設(shè)計(jì)）。由測(cè)量平差可知，對(duì)于間接平差有 （）1．零類設(shè)計(jì)（基準(zhǔn)設(shè)計(jì)）。由于觀測(cè)量不同，采用儀器也不同，其計(jì)算過(guò)程就都不一樣，很難用一個(gè)完整表達(dá)式表達(dá)出來(lái)，只能視情況而定。具體地說(shuō)就是采用下列的某一原則： 最大原則：在費(fèi)用一定的情況下，使控制網(wǎng)的質(zhì)量標(biāo)準(zhǔn)最大或者可靠性指標(biāo)能滿足一定要求時(shí)使精度指標(biāo)達(dá)到最高。如果在控制網(wǎng)的優(yōu)化設(shè)計(jì)中只用可靠性標(biāo)準(zhǔn)，而不用精度指標(biāo)、費(fèi)用指標(biāo)等作為目標(biāo)進(jìn)行設(shè)計(jì)，獲得合理的觀測(cè)方案就比較困難，這樣做會(huì)出現(xiàn)費(fèi)用高，優(yōu)化解不穩(wěn)定