资源描述：
《实验一_matlab完成水准网平差.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、.word格式,实验一matlab完成水准网平差实验数据：水准网有2个已知点，3个未知点，7个测段。已知点高程H1=5.016MH2=6.016h1=1.359;h2=2.009;h3=0.363;h4=1.012;h5=0.657;h6=0.238;h7=-0.595;S1=1.1S2=1.7S3=2.3S4=2.7S5=2.4S6=1.4S7=2.6求解（1）求个待定点高程,H5的高差中误差；3、4号点的高程中误差。课程设计内容1、平差程序设计思路：使用间接平差法求解（1）由题意知必要观测数t=3，选取3、4、5号点高程X
2、1、X2、X3为参数。（2）误差方程：V1=x1v2=x2v3=x1v4=x2v5=x2-x1+h2-h1-h5,专业.专注..word格式,v6=x3-x1v7=-x3(3)取1M的观测高程为单位权观测，即p=1/s；(4)求法方程：Nbbx-W=0Nbb=b’pbW=b’pl(5)求的平差值x=Nbb^-1*WL=l+VV=bx-l（6）高差权函数式：k=-x1+x2(6)求中误差：单位权中误差δ0，协因数阵Nbb^-1.求得中误差δ2、平差程序流程代码说明：h1=1.359;h2=2.009;h3=0.363;h4=1.
3、012;h5=0.657;h6=0.238;h7=-0.595;H1=5.016H2=6.016h=[h1h2h3h4h5h6h7]'s=[1.11.72.32.72.41.42.6]'B=[100;010;,专业.专注..word格式,100;010;-110;-101;00-1]p=diag(1./s)l=[0;0;4;3;7;2;0]W=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=[-110]'
4、q=f'*Q*fD=d*sqrt(q)D1=d*sqrt(Q)(3)平差程序流程代码说明:clccleardisp(‘观测高差，单位m’),专业.专注..word格式,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=[h1h2h3h4h5h6h7]'s=[1.11.72.32.72.41.42.6]'%S是线路长度disp(‘系数矩阵B、l’)B=[100;010;100;
5、010;-110;-101;00-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),专业.专注..word格式,H=h+V/1000%观测值的平差值Q=inv(Nbb)%观测值协因数阵n=7%观测值数t=3%必要观测数j=V'*p*V%计算单位权中误差的参数d=sqrt(j/4)%单位权中误差f=[-110]'%权函数系数阵q=f'*Q*f%权函数协因数阵D=d*sqrt(q)%高差中
6、误差D1=d*sqrt(Q)%高程中误差（4）计算结果：p=0.909100000000.588200000000.434800000000.37040000000,专业.专注..word格式,0.416700000000.714300000000.3846W=-2.60614.02781.4286Nbb=2.4748-0.4167-0.7143-0.41671.37530-0.714301.0989x=-0.24272.85521.1423V=-0.24272.8552-4.2427-0.1448,专业.专注..word格式
7、,-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=0即1.3588+0.6531-2.0119=0;0.3588+0.6531-1.0119=0；0.2374-0.5961+0.3588=0等式成立检核通过Q=0.53070.16080.34500.16080.77580.10450.34500.10451.1342j=19.7994d=2.2248f=-1,专业.专注..
8、word格式,10q=0.9850D=2.2080D1=1.62080.89211.30670.89211.95970.71921.30670.71922.3694,专业.专注.
