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

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.7S5=2.4 S6=1.4 S7=2.6求解( 1)求个待定点高程 ,H5 的高差中误差； 3、 4号点的高程中误差课程设计内容1、平差程序设计思路：使用间接平差法求解1) 由题意知必要观测数 t=3 ，选取 3、 4、 5 号点高程 X1、 X2、 X3 为参数2) 误差方程：V1=x1v2=x2v3=x1v4=x2 v5=x2-x1+h2-h1-h5v6=x3-x1v7=-x3(3) 取 1M 的观测高程为单位权观测，即 p=1/s ；(4) 求法方程： Nbbx-W=0 Nbb=b ’ pbW=b’ pl⑸求的平差值x=Nbb八-1*W L=l+V V=bx-l( 6)高差权函数式： k=-x1+x2(6)求中误差： 单位权中误差S 0,协因数阵Nbb八-1.求得中误差S2、平差程序流程代码说明 ：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 h7]'s=[1.1 1.7 2.3 2.7 2.4 1.4 2.6]'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;0]W=B'*p*lNbb=B'*p*B x=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 0]'q=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 % 已知点高程，单位 mS 是线路长度h=[h1 h2 h3 h4 h5 h6 h7]'s=[1.1 1.7 2.3 2.7 2.4 1.4 2.6]' %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;0]W=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 00 0 00 0.5882000000000000000 0 00 00.4348 00 0.37040.4167 0 00 0 0 00 0.7143 00 0 00 0 0.3846W =-2.60614.02781.4286Nbb =2.4748 -0.4167 -0.7143-0.4167 1.3753 0-0.7143 0 1.0989 x =-0.24272.85521.1423V =-0.24272.8552-4.2427-0.1448-3.9021-0.6151-1.1423H =1.35882.01190.35881.01190.65310.2374-0.5961成果检核：H1+H5-H2=0 H3+H5-H4=0;H6+H7+H3=(W 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.34500.1608 0.7758 0.10450.3450 0.1045 1.1342j = 19.7994d = 2.2248 f = -10q =0.9850D = 2.2080D1 =1.6208 0.8921 1.30670.8921 1.9597 0.71921.3067 0.7192 2.3694 。