资源描述：
《误差理论与测量平差基础》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、《误差理论与测量平差基础》实习报告王驩裕1420501201420050135东华理工大学测绘学院测量系12水准网间接平差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)
2、);D=B*E*B';SIGMA3=SIGMA*sqrt(D(5,5));end12导线网间接平差1.按间接平差法完成一导线网的平差计算。function[ZX,v,J,H]=nc(s,X,l,beta)L=dms2degrees(beta);alpha0=dms2degrees([2264459]);alpha1=alpha0+L(1)-180;alpha2=alpha1+L(2)-180;alpha3=alpha2+L(3)-180;alpha4=alpha3+L(4)-180;alpha=[al
3、pha1;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)00;a(2)b(2)-a(1)-b(1);a(3)b(3)-a(2)-b(2);00-a(3)-b(3);f(1)g(1)00;-f(1)-g(1)f(2)g(2);00-f(3)-g(3)];P1=diag([1111]);SIGMA0=5;
4、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,:);end122.求各导线点的坐标平差值极其点位中误差。function[ZX,SIGMA1,SIGMA2,SIGMA3,SIGMA4,SIGMA5]=pa(
5、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));12SIGMA2=SIGMA0*sqrt(G(3)+G(4));SIGMA3=SIGMA0*sqrt(G(5)+G(6));SIGMA4=SIGMA0*s
6、qrt(G(7)+G(8));SIGMA5=SIGMA0*sqrt(G(9)+G(10));end12三角网的间接平差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(
7、4));end121212测量平差软件的使用1.TXT导入控制网平差报告12[控制网概况]计算软件：南方平差易2005网名：bjw3计算日期：1970-01-01观测人：SOUTH记录人：计算者：SOUTH检查者：测量单位：SOUTHSURVEY备注：BJW平面控制网等级：城市二级，验前单位权中误差：2.50(s)已知坐标点个数：2未知坐标点个数：3未知边数：8最大点位误差[5]=0.0080(m)最小点位误差[3]=0.0045(m)平均点位误差=0.0066(m)最大点间误差=0.0114(m)最
8、大边长比例误差=693250平面网验后单位权中误差=0.57(s)[边长统计]总边长：24226.019(m),平均边长：3028.252(m),最小边长：2115.919(m),最大边长：4363.611(m)[闭合差统计报告]序号:<1>:中点多边形路径：[]极条件闭合差=0.9,限差=±5.3序号:<2>:大地四边形序号:<3>:闭合导线路径:2-3-4-1极条件闭合差=1.4,限差=±6.8路径：[2-3-1]角度闭合差=0.30(s),限差=±
