误差理论与测量平差基础..doc

《误差理论与测量平差基础》 实习报告 王驩裕 1420501 201420050135 东华理工大学测绘学院测量系 水准网间接平差 function [V,ZL,SIGMA1,SIGMA2,SIGMA3]=math(B,s,l,L,r) P=diag(1./s); NBB=B*P*B; W=B*P*l; x=inv(NBB)*W; V=B*x-l; ZL=L+V; SIGMA=sqrt(V*P*V/r); E=inv(NBB); SIGMA1=SIGMA*sqrt(E(1)); SIGMA2=SIGMA*sqrt(E(2,2)); D=B*E*B; SIGMA3=SIGMA*sqrt(D(5,5)); end 导线网间接平差 按间接平差法完成一导线网的平差计算。 function [ZX,v,J,H]=nc(s,X,l,beta) L=dms2degrees(beta); alpha0=dms2degrees([226 44 59]); alpha1=alpha0+L(1)-180; alpha2=alpha1+L(2)-180; alpha3=alpha2+L(3)-180; alpha4=alpha3+L(4)-180; alpha=[alpha1;alpha2;alpha3]; e=alpha; c=180*3600*sind(e)/(pi*100); d=-180*3600*cosd(e)/(pi*100); a=c./s; b=d./s; f=cosd(e); g=sind(e); B=[a(1) b(1) 0 0;a(2) b(2) -a(1) -b(1); a(3) b(3) -a(2) -b(2);0 0 -a(3) -b(3); f(1) g(1) 0 0;-f(1) -g(1) f(2) g(2);0 0 -f(3) -g(3)]; P1=diag([1 1 1 1]); SIGMA0=5; SIGMAS=0.5*sqrt(s); D=1./SIGMAS; T=(SIGMA0*D).^2; P2=diag(T); P=blkdiag(P1,P2); NBB=B*P*B; W=B*P*l; x=inv(NBB)*W; ZX=X+x; v=B*x-l; ZL=[L;s]+v/3600; J=degrees2dms(ZL(1:4,:)); H=ZL(5:7,:); end 2.求各导线点的坐标平差值极其点位中误差。 function [ZX,SIGMA1,SIGMA2,SIGMA3,SIGMA4,SIGMA5]=pa(B,S,l,X) P1=eye(10); P2=diag(100000./S); P=blkdiag(P1,P2); NBB=B*P*B; W=B*P*l; x=inv(NBB)*W; ZX=X+x; E=inv(NBB); V=B*x-l; SIGMA0=sqrt(V*P*V/7); G=diag(E); SIGMA1=SIGMA0*sqrt(G(1)+G(2)); SIGMA2=SIGMA0*sqrt(G(3)+G(4)); SIGMA3=SIGMA0*sqrt(G(5)+G(6)); SIGMA4=SIGMA0*sqrt(G(7)+G(8)); SIGMA5=SIGMA0*sqrt(G(9)+G(10)); end 三角网的间接平差 function [ZX,SIGMA1,SIGMA2]=CH(B,l,X) P=eye(18); NBB=B*P*B; W=B*P*l; x=inv(NBB)*W; ZX=X+x; Q=inv(NBB); E=diag(Q); V=B*x-l; SIGMA0=sqrt(V*P*V/12); SIGMA1=SIGMA0*sqrt(E(1)+E(2)); SIGMA2=SIGMA0*sqrt(E(3)+E(4)); end  测量平差软件的使用 TXT导入 控 制 网 平 差 报 告 [控制网概况] 计算软件：南方平差易2005 网名：bjw3 计算日期：1970-01-01 观测人：SOUTH 记录人： 计算者：SOUTH 检查者： 测量单位：SOUTH SURVEY 备注：BJW 平面控制网等级：城市二级，验前单位权中误差：2.50(s) 已知坐标点个数：2 未知坐标点个数：3 未知边数：8 最大点位误差[5] = 0.0080 (m) 最小点位误差[3] = 0.0045 (m) 平均点位误差 = 0.0066 (m) 最大点间误差 = 0.0114(m) 最大边长比例误差 = 693250 平面网验后单位权中误差