




版權(quán)說(shuō)明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1、MATLAB程序展示水準(zhǔn)網(wǎng)間接平差function v,ch,cx=szw(s,h,B,x0,d,n,t,f)%t匕函數(shù)用于計(jì)算改正數(shù),高差中誤差,高程中誤差 p=diag(1./s);優(yōu)義權(quán)陣disp( P= ) disp(p) l=h-B*x0-d; W=B*p*l; Nbb=B*p*B; disp( NBB= ) disp(Nbb) x=(inv(Nbb)*W; disp(改正數(shù)) v=B*x-l;%改正數(shù)disp(v); c0=sqrt(v*p*v)/(n-t);%|位中誤差Nbb=B*p*B; Qh=f*(inv(Nbb)*f;%h5 的協(xié)因數(shù)陣disp(高差平差值中誤差) ch=c
2、0*sqrt(Qh);%高差平差值中誤差disp(ch); Nbbn=inv(Nbb);濰逆矩陣disp(高程平差值中誤差) cx=c0*sqrt(diag(Nbbn);%高程平差值中誤差 TOC o 1-5 h z disp(cx); return loadB.txt;loadd.txt;loads.txt;load x0.txt ; loadh.txt;loadn.txt;loadt.txtloadf.txt v,ch,cx=szw(s,h,B,x0,d,n,t,f);G T11 tMJ XO.txtq E F.0. r, n, r, ft。+ U *。: Preqrdm Fik;卜 W、
3、T,B R;015b 卜 bin 卜 同口pptSC)shuiztujlASJlE 的正敢-O.OOT2O. OQ:0|-0. 0042-O.OODJiQ. 003U-C.001B高差平差I(lǐng)fir誤差 Q 0X2層卷平差值中誤差0.COJb 也OKO Tleir fhuiHhunwaEF=0,順1rn00ft0口順蛇UdQ0000.42480000C00.加0口0C0caL 416700Cnc0Q.網(wǎng)Q陽(yáng)典3,0106L4167-C 4M7L 3753-0.0D030.。此-0. 0342-C.ODQl-C. 0D39H3. 0018言壬二差百中迄差晨在中差也中運(yùn)甚(,00160.0D2OMA
4、TLAB程序展示一一導(dǎo)線網(wǎng)間接平差loadn.txt%觀測(cè)總數(shù)loadt.txt%必要觀測(cè)數(shù)loaddws.txt%入點(diǎn)位數(shù)loadcsgs.txt%導(dǎo)入?yún)?shù)個(gè)數(shù)loadbs.txt%導(dǎo)入待測(cè)邊數(shù)loaddcdh1.txt%1一個(gè)待測(cè)點(diǎn)號(hào)loadx1y1.txt%導(dǎo)入已知點(diǎn)位一的坐標(biāo)loadxydws.txt溫入末點(diǎn)位的坐標(biāo)loaddeg.txt%11入由觀測(cè)角計(jì)算而得的方位角load s.txt%入邊長(zhǎng)觀測(cè)值loadgcj.txt%各點(diǎn)位觀測(cè)角digits(12);digits(15);x1=x1y1(1,1);y1=x1y1(1,2);x0=zeros(dws,1);y0=zeros(dw
5、s,1);a=zeros(1,3);x0=x1;y0(1)=y1;for i=2:(dws-1)a=deg(i-1,:);x2,y2=zbzs(x1,y1,s(i-1),a);x1=x2;y1=y2;x0(i)=x1;y0(i)=y1;endx1=xydws(1,1);y1=xydws(1,2);x0(dws)=x1;%各點(diǎn)位x的近似值y0(dws)=y1;%各點(diǎn)位y的近似值s0=zeros(dws-1,1);%#各點(diǎn)位之間距離的近似值for i=1:(dws-1)s0(i)=sqrt(xO(i)-x0(i+1).A2+(y0(i)-y0(i+1).A2);endB1,L1=jsB1L1(x0
6、,yO,s0,s,bs,dcdh1,csgs);B2,L2,ajs=jsB2L2(x0,yO,s0,dws,dcdh1,csgs,gcj,bs);B=B1;B2;L=L1;L2;juw=5; 糊悌中誤差buw=0.5*sqrt(s);%邊中誤差P=zeros(1,length(gcj);for i=1:length(gcj) P(i)=1./(juw.A2);endP=1./(buw.A2),P;P=diag(P)*(juw.A2);NBB=B*P*B;W=B*P*L;x=pinv(NBB)*W;disp( xE yE xF yF=)disp(x)X=zeros(1,length(x);X=x
7、(1)+x0(2) x(2)+y0(2) x(3)+x0(3) x(4)+y0(3);disp( E點(diǎn)坐標(biāo)平差值=)disp(X(1,1:2)disp( F點(diǎn)坐標(biāo)平差值=)disp(X(1,3:4)V=B*x-L;uwef=zwc(V,P,n,t,B);%點(diǎn)位精度disp(點(diǎn)位精度 uwe uwf=) disp(uwef) TOC o 1-5 h z disp(各觀測(cè)量改正數(shù)Vb1 Vb2 Vb3 Vb4 Vj1 Vj2 Vj3=)disp(V)spcz=s+V(1:3,1);%邊長(zhǎng)觀測(cè)值平差值disp(邊長(zhǎng)觀測(cè)平差值 S1 S2 S3=)disp(spcz)gcjpcz=gcj+V(4:7,
8、1);disp(角度觀測(cè)平差值J1 J2 J3 J4=)disp(gcjpcz)function A=deg2rad(a)%角度轉(zhuǎn)弧度A=a*pi/180;pi/10800;pi/648000; Returnfunction x2,y2=zbzs(x1,y1,s,a)% 標(biāo)正算A=deg2rad(a);x2=x1+s*cos(A);y2=y1+s*sin(A); Returnfunction B1,L1=jsB1L1(x0,y0,s0,s,bs,dcdh1,csgs)%# 算邊長(zhǎng)系數(shù)矩陣BB1=zeros(bs,2*(bs+1);L1=s-s0;for i=1:bsk=1;for j=1+2*
9、(i-1):2:(2*csgs+2*(i-1)B1(i,j)=(x0(i+1)-x0(i)*(-1k./s0(i);B1(i,j+1)=(y0(i+1)-y0(i)*(-1)Ak./s0(i);k=k+1;endendB1=B1(1:bs,2*dcdh1-1:2*dcdh1-2+2*csgs);returnfunction B2,L2,ajs=jsB2L2(x0,y0,s0,dws,dcdh1,csgs,gcj,bs)系數(shù)矩陣B2=zeros(dws,2*dws);axs=zeros(dws+1,4);for i=2:dwsaxs(i,1)=3600*(y0(i)-y0(i-1)./(s0(i
10、-1).A2);axs(i,2)=-3600*(x0(i)-x0(i-1)./(s0(i-1).A2);axs(i,3)=3600*(y0(i)-y0(i-1)./(s0(i-1).A2);axs(i,4)=-3600*(x0(i)-x0(i-1)./(s0(i-1).A2);endfor i=1:dwsif i=1B2(i,1:4)=axs(i,:)+axs(i+1,:);elseif i=dwsB2(i,2*dws-3:2*dws)=axs(i,:)+axs(i+1,:);elseB2(i,1+2*(i-2):4+2*(i-2)=axs(i,:);B2(i,1+2*(i-1):4+2*(i
11、-1)=B2(i,1+2*(i-1):4+2*(i-1)+axs(i+1,:);endendendB2=B2(1:dws,2*dcdh1-1:2*dcdh1-2+2*csgs);ajs=zeros(bs+2,1);ajs(1)=46.74972222;for i=2:bs+1ajs(i)=atan(y0(i-1)-y0(i)./(x0(i-1)-x0(i)*180./pi;endajs(bs+2)=144.7675;L2=zeros(bs+1,1);for i=1:bs+1L2(i)=gcj(i)-180+ajs(i)-ajs(i+1);%計(jì)算角度的%近似方位角if L2(i)90L2(i)=
12、L2(i)-180;elseif L2(l)l Vb2 Vb3 TH Vjl Tj2骨3=0.00130.001T o.oin-0.00300.0021-0.0012-0.0002詁長(zhǎng)Wil甲*SSI S2 S204.96382D0.1317345. J64T弟房A?差I(lǐng)fiJ】J2 J3 Js230.5406IB0.0138170.6649236,818MATLAB程序展示一一三角網(wǎng)間接平差loadxqs.txt%各起始點(diǎn)x值loadyqs.txt%各起始點(diǎn)y值loadgcj.txt%觀測(cè)角loada0.txt%近似坐標(biāo)方位角loadL1.txtloadL2.txtloads0.txtloa
13、dB1.txtloadB2.txtloadn.txt%觀測(cè)總數(shù)loadt.txt%必要觀測(cè)數(shù)x0,y0=jszb(xqs,yqs,gcj);P=DP(s0);disp( p=)disp(P)B=B1;B2;L=L1;L2;NBB=B*P*B;disp( NBB=)disp(NBB)W=B*P*L;x=(inv(NBB)*W;disp(坐標(biāo)改正數(shù))disp(x)x=x./10;X=x0(1)+x(1);y0(1)+x(2);x0(2)+x(3);y0(2)+x(4);disp(坐標(biāo)平差值X1 Y1 X2 Y2T=)disp(X)x=x*10;v=B*x-L;S=s0+v(1:7,1)./10;d
14、isp(邊長(zhǎng)平差值 P1A P1B P1C P1P2 P2A P2C P2Ddisp(S)GCJ=gcj;GCJ(1:18,3)=gcj(1:18,3)+v(1:18,1);disp(觀測(cè)角的平差值,)disp(GCJ)uwx=jszwc(v,n,t,NBB);uwp1=sqrt(uwx(1).A2+uwx(2).A2);disp( P1點(diǎn)的點(diǎn)位中誤差)disp(uwp1)uwp2=sqrt(uwx(3).A2+uwx(4).A2);disp( P2點(diǎn)的點(diǎn)位中誤差)disp(uwp2)function A=deg2rad(a)%角度轉(zhuǎn)弧度A=a*pi/180;pi/10800;pi/64800
15、0;Returnfunction x0,y0=jszb(xqs,yqs,gcj)%2i 似坐標(biāo)gcj=deg2rad(gcj);x0=zeros(2,1);y0=zeros(2,1);x0(1)=(xqs(2)./tan(gcj(16)+xqs(3)./tan(gcj(18)-yqs(3)+yqs(2)./(1./ tan(gcj(18)+1./tan(gcj(16);x0(2)=(xqs(4)./tan(gcj(7)+xqs(3)./tan(gcj(9)-yqs(4)+yqs(3)./(1./tan(gcj(9)+1./tan(gcj);y0(1)=(yqs(2)./tan(gcj(16)+
16、yqs(3)./tan(gcj(18)+xqs(3)-xqs(2)./(1./ tan(gcj(18)+1./tan(gcj(16);y0(2)=(yqs(4)./tan(gcj(7)+yqs(3)./tan(gcj(9)+xqs(4)-xqs(3)./(1./ta n(gcj(9)+1./tan(gcj);Returnfunction P=DP(s0)優(yōu)權(quán)P=zeros(1,25);P(1,1:7)=1./s0;P(1,8:25)=1;P=diag(P);Returnfunction uwx=jszwc(v,n,t,NBB)%#算中誤差v=v(8:25,1);uw0=sqrt(v*v./(n
17、-t);QXX=diag(inv(NBB);uwx=uw0*sqrt(QXX);return運(yùn)算結(jié)果展示中,口 & - PrgE 日、2眼MB * 三片hlDieft書E學(xué)。療死rw府.i DittjCm 口 gm,j j E E / 11 nnwE -mid j&4 口 L2Mn bd| BN FJltKW J-L壬1。丸4. POOLD4。494:q niqgtgaogBDQBadd.DDDZddandnd0D1D.DIMLD00D4Bonia aoiDzodd41。q1口 & DOMaDQaDBdDa tomDODadiadaaLnabaqsooiadodni.MI口。pUDgLM 口巾
18、aDDd0OnD0D00i400dD410D400004t。QOOQQg0Q|DODDa0ondoooDa0D14D0000IQQ*QQ電。QQ*CiDDDD4t0D140e0D41t0D1400004IguiqpgQocit ?DBgPJUQQD0DBa口。口口OGtafioabdfi1 1 taai 1小311 y岸山11 90 列0J3000C00Coooooooooo Qdaouoouoo ooaoaqaoaa nooooconaOflOOOOOOOO 0000000000 000000000(1 OOOcOOOodc OQOooooooa h occaiDooooooa0L00093
19、000C03(0b1.00030000000ooa i. oooouoodja00Q0J.OOOUoouooonaao l. 口口口口ooaaonooai. oooonao000000 J. 00001000000000 LOOTO0ooooooooo l.OCOC OflQOOOCOOC oooooooooa ai3oocoocOdOOODCOOC 000000000BB=94.6088 - -22.1152 -11.4010-6.9616生林改正數(shù)-22. 1152 70.5086 -6.8950 -8.4233-11.4010-6. B9E095. 9311-2d 1880-6. 9C4
20、6|-8.4233-20. JE9006.634121至25列-0.10200000o|2.322100000-.2O7200000-0.53420000000000一標(biāo)平集值【XI丫】X2 Y2JT=00000I. Ce-b04 京0000000000000001.3180000003.7335000001.557800000000004.43910000000000邊長(zhǎng)平差值P1AP1B PIC P1P2 P2A P2c P2D00000i.Ccoa *0000000000000007. 3864000006.91161.000000005.89410 J.OOQO0007. 4d9800
21、1.0000005.9198000 J.CODO000001.00007.43925.9665觀測(cè)甬的平差值126.000014.000026.914223. 000039.000050.704530.00005.000038.5663117.000022.000046.563531.000026.000048.544731.000010.OOOO26.2521ZZ. 00002.000039.8315130.00003.000014.990327.000053.000061.514265. 000055.00000.095567. 00002.000060.213247. 00002.000
22、012.355346. 000038.000055.821566. 000034.000055.252166.000046.00009.879429. 000058.000036.7605120.00008.000031.708829. 000052.000054.0305P1點(diǎn)的點(diǎn)位中誤差 0. 2117F2點(diǎn)的點(diǎn)位中誤強(qiáng) 0. 2134控制網(wǎng)概況計(jì)算軟件:南方平差易 2005網(wǎng)名:三角網(wǎng)計(jì)算日期:2015年11月12日星期四觀測(cè)人:記錄人:計(jì)算者:檢查者:測(cè)量單位:2.50(s)備注:平面控制網(wǎng)等級(jí):國(guó)家四等,驗(yàn)前單位權(quán)中誤差:已知坐標(biāo)點(diǎn)個(gè)數(shù):4未知坐標(biāo)點(diǎn)個(gè)數(shù):2未知邊數(shù):0最大點(diǎn)位誤差P
23、2 = 27.4580 (m)最小點(diǎn)位誤差P1 = 27.2219 (m)平均點(diǎn)位誤差 =27.3400 (m)最大點(diǎn)間誤差 =49.3635(m)最大邊長(zhǎng)比例誤差=267平面網(wǎng)驗(yàn)后單位權(quán)中誤差=1053.86 (s)閉合差統(tǒng)計(jì)報(bào)告序號(hào):中點(diǎn)多邊形路徑:極條件閉合差=1764.9,限差=4.1序號(hào):中點(diǎn)多邊形路徑:路徑:極條件閉合差=1196.1,限差=4.3序號(hào):閉合導(dǎo)線路徑:B-P1-A角度閉合差=-2433.62(s),限差=土 8.66(s)序號(hào):閉合導(dǎo)線路徑:P1-P2-A角度閉合差=2439.16(s),限差= 8.66(s)序號(hào):閉合導(dǎo)線路徑:P2-D-A角度閉合差=-36.66
24、(s),限差=土 8.66(s)序號(hào):閉合導(dǎo)線路徑:P1-C-B角度閉合差=2405.56(s),限差=土 8.66(s)序號(hào):閉合導(dǎo)線路徑:C-P2-P1角度閉合差=-2404.72(s),限差= 8.66(s)序號(hào):閉合導(dǎo)線路徑:C-D-P2角度閉合差=-30.28(s),限差=土 8.66(s)起算點(diǎn)數(shù)據(jù)表點(diǎn)名X(m)Y(m)H(m)備注A9684.280043836.8200B10649.550031996.5000C19063.660037818.8600D17814.630049923.1900方向觀測(cè)成果表測(cè)站照準(zhǔn)方向值(dms)改正數(shù)(s)平差后值(dms)備注AB0.000000P123.663028-1611.8923.393839P290.710084-732.7990.584805D122.157306-404.95122.092811BC0.000000P129.882056-2496.8829.464368A59.978362-2381.4659.584216P1A0.000000B126.240028-937.15126.082313C246.382000-475.96246.302404P2312.963861-1847.72313.055089CD0.000000P222.045278-534
溫馨提示
- 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 生物分類與生態(tài)系統(tǒng)研究方法探討試題及答案
- 2024年陪診師考試疾病預(yù)防試題及答案
- 2024陪診師考試心理支持試題及答案
- 去化樓層培訓(xùn)
- 城市污染的成因與防治方法試題及答案
- 電子商務(wù)教師資格證應(yīng)試策略及試題答案
- 深入學(xué)習(xí)監(jiān)理工程師試題及答案
- 黑龍江省七臺(tái)河市勃利縣小學(xué)2024-2025學(xué)年數(shù)學(xué)五下期末經(jīng)典試題含答案
- 黑龍江省佳木斯市同江市2025年四下數(shù)學(xué)期末學(xué)業(yè)質(zhì)量監(jiān)測(cè)試題含解析
- 黑龍江省雙鴨山市尖山區(qū)一中2024-2025學(xué)年高三下學(xué)期調(diào)研物理試題含解析
- 2023機(jī)關(guān)公文寫作與處理PPT模板
- 2022年撫順特殊鋼股份有限公司招聘筆試試題及答案解析
- 兒童抑郁量表CDI
- 幼兒數(shù)字1-100字帖練習(xí)
- 細(xì)胞生物學(xué)-7細(xì)胞信號(hào)轉(zhuǎn)導(dǎo)課件
- 心電監(jiān)護(hù)操作評(píng)分標(biāo)準(zhǔn)
- 攪拌站安全培訓(xùn)試卷
- Q∕SY 02098-2018 施工作業(yè)用野營(yíng)房
- 浙教版勞動(dòng)五年級(jí)下冊(cè) 項(xiàng)目三 任務(wù)三 環(huán)保小車我來(lái)造 教案
- 隔離開(kāi)關(guān)培訓(xùn)課件
- 圖像融合技術(shù)中英文對(duì)照外文翻譯文獻(xiàn)
評(píng)論
0/150
提交評(píng)論