程序及运算结果

1、MATLAB程序展示水准网间接平差function v,ch,cx=szw(s,h,B,x0,d,n,t,f)%t匕函数用于计算改正数，高差中误差,高程中误差 p=diag(1./s);优义权阵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(改正数) v=B*x-l;%改正数disp(v); c0=sqrt(v*p*v)/(n-t);%|位中误差Nbb=B*p*B; Qh=f*(inv(Nbb)*f;%h5 的协因数阵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.网Q阳典3,0106L4167-C 4M7L 3753-0.0D030.。此-0. 0342-C.ODQl-C. 0D39H3. 0018言壬二差百中迄差晨在中差也中运甚(,00160.0D2OMA

4、TLAB程序展示一一导线网间接平差loadn.txt%观测总数loadt.txt%必要观测数loaddws.txt%入点位数loadcsgs.txt%导入参数个数loadbs.txt%导入待测边数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(各观测量改正数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)%# 算边长系数矩阵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)系数矩阵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程序展示一一三角网间接平差loadxqs.txt%各起始点x值loadyqs.txt%各起始点y值loadgcj.txt%观测角loada0.txt%近似坐标方位角loadL1.txtloadL2.txtloads0.txtloa

13、dB1.txtloadB2.txtloadn.txt%观测总数loadt.txt%必要观测数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(坐标改正数)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)优权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口。口口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生林改正数-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控制网概况计算软件：南方平差易 2005网名：三角网计算日期：2015年11月12日星期四观测人：记录人：计算者：检查者：测量单位：2.50(s)备注：平面控制网等级：国家四等，验前单位权中误差:已知坐标点个数：4未知坐标点个数：2未知边数：0最大点位误差P

23、2 = 27.4580 (m)最小点位误差P1 = 27.2219 (m)平均点位误差 =27.3400 (m)最大点间误差 =49.3635(m)最大边长比例误差=267平面网验后单位权中误差=1053.86 (s)闭合差统计报告序号:中点多边形路径：极条件闭合差=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)起算点数据表点名X(m)Y(m)H(m)备注A9684.280043836.8200B10649.550031996.5000C19063.660037818.8600D17814.630049923.1900方向观测成果表测站照准方向值(dms)改正数(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