误差理论与测量平差实习报告

出单位权中误差,评定各观测值及所求高程精度。四. 计算结果 1. 编写程序代码 n=7,t=3,r=n-t=4 B=[1 0 0;0 1 0;1 0 0;0 1 0;-1 1 0;-1 0 1;0 0 -1] P=[9 00000 0;0 60000 0;0 04000 0;0 00400 0;0 00040 0;0 00007 0;0 00000 4] l=[0;0;4;3;7;2;0] N BB =B'*P*B W=B'*P*l x=inv(N BB )*W V=B*x-l 8 h=[1359;2009;363;1012;657;238;-595] H=h+V X 0 =[6365;7025;6611] X=X 0 +x m 0 =+-sqrt(V'*P*V/r) f =[-1;1;0] Q=f '*inv(Nbb)* f m=m 0 *sqrt(Q) Qxx=inv(N BB )m 1 =m 0 *sqrt(0.0543) m 2 =m 0 *sqrt(0.0759) 2. 计算机计算的结果>> B=[1 0 0;0 1 0;1 0 0;0 1 0;-1 1 0;-1 0 1;0 0 -1] B=100010100010 -110 -101 9 00 -1 >> P=[9 00000 0;0 60000 0;0 04000 0;0 00400 0;0 00040 0;0 00007 0;0 00000 4] P=9000000 0600000 0040000 0004000 0000400 0000070 0000004 >> l=[0;0;4;3;7;2;0] l=0043 10 720 >>N BB =B'*P*B N BB = 24 -4 -7 -4 140 -70 11 >> W=B'*P*l W= -26 40 14 >> x=inv(N BB )*W x=