实验一-matlab完成水准网平差

实验一 matlab完成水准网平差实验数据：水准网有2个已知点，3个未知点，7个测段。 已知点高程H1=5.016M H2=6.016h1=1.359; h2=2.009; h3=0.363; h4=1.012;h5=0.657; h6=0.238; h7=-0.595;S1=1.1 S2=1.7 S3=2.3 S4=2.7 S5=2.4 S6=1.4 S7=2.6求解（1）求个待定点高程,H5的高差中误差；3、4号点的高程中误差。课程设计内容1、平差程序设计思路：使用间接平差法求解（1） 由题意知必要观测数t=3，选取3、4、5号点高程X1、X2、X3为参数。（2） 误差方程：V1=x1 v2=x2 v3=x1v4=x2v5=x2-x1+h2-h1-h5v6=x3-x1v7=-x3(3)取1M 的观测高程为单位权观测，即 p=1/s；(4)求法方程：Nbbx-W=0 Nbb=bpb W=bpl(5)求的平差值x=Nbb-1*W L=l+V V=bx-l （6）高差权函数式： k=-x1+x2(6)求中误差： 单位权中误差0，协因数阵Nbb-1.求得中误差2、平差程序流程代码说明 ：h1=1.359;h2=2.009;h3=0.363;h4=1.012;h5=0.657;h6=0.238;h7=-0.595;H1=5.016H2=6.016h=h1 h2 h3 h4 h5 h6 h7s=1.1 1.7 2.3 2.7 2.4 1.4 2.6B=1 0 0 ;0 1 0; 1 0 0;0 1 0 ; -1 1 0 ; -1 0 1 ;0 0 -1 p=diag(1./s)l=0;0;4;3;7;2;0W=B*p*lNbb=B*p*Bx=inv(Nbb)*WV=(B*x-l)H=h+V/1000Q=inv(Nbb) n=7;t=3;j=V*p*Vd= sqrt(j/4)f=-1 1 0q=f*Q*fD=d*sqrt(q)D1=d*sqrt(Q) (3) 平差程序流程代码说明:clccleardisp(观测高差，单位m)h1=1.359;h2=2.009;h3=0.363;h4=1.012;h5=0.657;h6=0.238;h7=-0.595;H1=5.016 % 已知点高程，单位mH2=6.016 % 已知点高程，单位mh=h1 h2 h3 h4 h5 h6 h7s=1.1 1.7 2.3 2.7 2.4 1.4 2.6 % S是线路长度disp(系数矩阵B、l)B=1 0 0 ;0 1 0; 1 0 0;0 1 0 ; -1 1 0 ; -1 0 1 ;0 0 -1 p=diag(1./s) %定义权阵l=0;0;4;3;7;2;0W=B*p*lNbb=B*p*Bdisp(参数的解)x=inv(Nbb)*WV=(B*x-l) % 误差方程(mm)H=h+V/1000 % 观测值的平差值Q=inv(Nbb) %观测值协因数阵n=7 %观测值数t=3 % 必要观测数j=V*p*V %计算单位权中误差的参数d= sqrt(j/4) %单位权中误差f=-1 1 0 %权函数系数阵q=f*Q*f %权函数协因数阵D=d*sqrt(q) %高差中误差D1=d*sqrt(Q) %高程中误差（4）计算结果： p = 0.9091 0 0 0 0 0 0 0 0.5882 0 0 0 0 0 0 0 0.4348 0 0 0 0 0 0 0 0.3704 0 0 0 0 0 0 0 0.4167 0 0 0 0 0 0 0 0.7143 0 0 0 0 0 0 0 0.3846W = -2.6061 4.0278 1.4286Nbb = 2.4748 -0.4167 -0.7143 -0.4167 1.3753 0 -0.7143 0 1.0989x = -0.2427 2.8552 1.1423V = -0.2427 2.8552 -4.2427 -0.1448 -3.9021 -0.6151 -1.1423H =1.3588 2.0119 0.3588 1.0119 0.6531 0.2374 -0.5961成果检核：H1+H5-H2=0；H3+H5-H4=0;H6+H7+H3=0即1.3588+0.6531-2.0119=0; 0.3588+0.6531-1.0119=0；0.2374-0.5961+0.3588=0等式成立 检核通过Q = 0.5307 0.1608 0.3450 0.1608 0.7758 0.1045 0.3450 0.1045 1.1342j = 19.7994d =