程序及運算結果

1、MATLAB程序展示水準網(wǎng)間接平差function v,ch,cx=szw(s,h,B,x0,d,n,t,f)%t匕函數(shù)用于計算改正數(shù)，高差中誤差,高程中誤差 p=diag(1./s);優(yōu)義權陣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高差平差Ifir誤差 Q 0X2層卷平差值中誤差0.COJb 也OKO Tleir fhuiHhunwaEF=0,順1rn00ft0口順蛇UdQ0000.42480000C00.加0口0C0caL 416700Cnc0Q.網(wǎng)Q陽典3,0106L4167-C 4M7L 3753-0.0D030.。此-0. 0342-C.ODQl-C. 0D39H3. 0018言壬二差百中迄差晨在中差也中運甚(,00160.0D2OMA

4、TLAB程序展示一一導線網(wǎng)間接平差loadn.txt%觀測總數(shù)loadt.txt%必要觀測數(shù)loaddws.txt%入點位數(shù)loadcsgs.txt%導入?yún)?shù)個數(shù)loadbs.txt%導入待測邊數(shù)loaddcdh1.txt%1一個待測點號loadx1y1.txt%導入已知點位一的坐標loadxydws.txt溫入末點位的坐標loaddeg.txt%11入由觀測角計算而得的方位角load s.txt%入邊長觀測值loadgcj.txt%各點位觀測角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;%各點位x的近似值y0(dws)=y1;%各點位y的近似值s0=zeros(dws-1,1);%#各點位之間距離的近似值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點坐標平差值=)disp(X(1,1:2)disp( F點坐標平差值=)disp(X(1,3:4)V=B*x-L;uwef=zwc(V,P,n,t,B);%點位精度disp(點位精度 uwe uwf=) disp(uwef) TOC o 1-5 h z disp(各觀測量改正數(shù)Vb1 Vb2 Vb3 Vb4 Vj1 Vj2 Vj3=)disp(V)spcz=s+V(1:3,1);%邊長觀測值平差值disp(邊長觀測平差值 S1 S2 S3=)disp(spcz)gcjpcz=gcj+V(4:7,

8、1);disp(角度觀測平差值J1 J2 J3 J4=)disp(gcjpcz)function A=deg2rad(a)%角度轉弧度A=a*pi/180;pi/10800;pi/648000; Returnfunction x2,y2=zbzs(x1,y1,s,a)% 標正算A=deg2rad(a);x2=x1+s*cos(A);y2=y1+s*sin(A); Returnfunction B1,L1=jsB1L1(x0,y0,s0,s,bs,dcdh1,csgs)%# 算邊長系數(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);%計算角度的%近似方位角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詁長Wil甲*SSI S2 S204.96382D0.1317345. J64T弟房A?差IfiJ】J2 J3 Js230.5406IB0.0138170.6649236,818MATLAB程序展示一一三角網(wǎng)間接平差loadxqs.txt%各起始點x值loadyqs.txt%各起始點y值loadgcj.txt%觀測角loada0.txt%近似坐標方位角loadL1.txtloadL2.txtloads0.txtloa

13、dB1.txtloadB2.txtloadn.txt%觀測總數(shù)loadt.txt%必要觀測數(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(坐標改正數(shù))disp(x)x=x./10;X=x0(1)+x(1);y0(1)+x(2);x0(2)+x(3);y0(2)+x(4);disp(坐標平差值X1 Y1 X2 Y2T=)disp(X)x=x*10;v=B*x-L;S=s0+v(1:7,1)./10;d

14、isp(邊長平差值 P1A P1B P1C P1P2 P2A P2C P2Ddisp(S)GCJ=gcj;GCJ(1:18,3)=gcj(1:18,3)+v(1:18,1);disp(觀測角的平差值,)disp(GCJ)uwx=jszwc(v,n,t,NBB);uwp1=sqrt(uwx(1).A2+uwx(2).A2);disp( P1點的點位中誤差)disp(uwp1)uwp2=sqrt(uwx(3).A2+uwx(4).A2);disp( P2點的點位中誤差)disp(uwp2)function A=deg2rad(a)%角度轉弧度A=a*pi/180;pi/10800;pi/64800

15、0;Returnfunction x0,y0=jszb(xqs,yqs,gcj)%2i 似坐標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)權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運算結果展示中，口 & - PrgE 日、2眼MB * 三片hlDieft書E學。療死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口?？诳贠Gtafioabdfi1 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一標平集值【XI丫】X2 Y2JT=00000I. Ce-b04 京0000000000000001.3180000003.7335000001.557800000000004.43910000000000邊長平差值P1AP1B PIC P1P2 P2A P2c P2D00000i.Ccoa *0000000000000007. 3864000006.91161.000000005.89410 J.OOQO0007. 4d9800

21、1.0000005.9198000 J.CODO000001.00007.43925.9665觀測甬的平差值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點的點位中誤差 0. 2117F2點的點位中誤強 0. 2134控制網(wǎng)概況計算軟件：南方平差易 2005網(wǎng)名：三角網(wǎng)計算日期：2015年11月12日星期四觀測人：記錄人：計算者：檢查者：測量單位：2.50(s)備注：平面控制網(wǎng)等級：國家四等，驗前單位權中誤差:已知坐標點個數(shù)：4未知坐標點個數(shù)：2未知邊數(shù)：0最大點位誤差P

23、2 = 27.4580 (m)最小點位誤差P1 = 27.2219 (m)平均點位誤差 =27.3400 (m)最大點間誤差 =49.3635(m)最大邊長比例誤差=267平面網(wǎng)驗后單位權中誤差=1053.86 (s)閉合差統(tǒng)計報告序號:中點多邊形路徑：極條件閉合差=1764.9，限差=4.1序號:中點多邊形路徑：路徑：極條件閉合差=1196.1，限差=4.3序號:閉合導線路徑：B-P1-A角度閉合差=-2433.62(s),限差=土 8.66(s)序號:閉合導線路徑：P1-P2-A角度閉合差=2439.16(s),限差= 8.66(s)序號:閉合導線路徑：P2-D-A角度閉合差=-36.66

24、(s),限差=土 8.66(s)序號:閉合導線路徑：P1-C-B角度閉合差=2405.56(s),限差=土 8.66(s)序號:閉合導線路徑：C-P2-P1角度閉合差=-2404.72(s),限差= 8.66(s)序號:閉合導線路徑:C-D-P2角度閉合差=-30.28(s),限差=土 8.66(s)起算點數(shù)據(jù)表點名X(m)Y(m)H(m)備注A9684.280043836.8200B10649.550031996.5000C19063.660037818.8600D17814.630049923.1900方向觀測成果表測站照準方向值(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