2022年实验一matlab完成水准网平差 .pdf
-
资源ID：27194349
资源大小：52.70KB
全文页数：11页
- 资源格式： PDF
下载积分：4.3金币
快捷下载
会员登录下载
微信登录下载
三方登录下载：
微信扫一扫登录
|
友情提示
2、PDF文件下载后，可能会被浏览器默认打开，此种情况可以点击浏览器菜单，保存网页到桌面，就可以正常下载了。
3、本站不支持迅雷下载，请使用电脑自带的IE浏览器，或者360浏览器、谷歌浏览器下载即可。
4、本站资源下载后的文档和图纸-无水印,预览文档经过压缩，下载后原文更清晰。
5、试题试卷类文档，如果标题没有明确说明有答案则都视为没有答案，请知晓。
|
2022年实验一matlab完成水准网平差 .pdf
实验一matlab 完成水准网平差实验数据：水准网有2 个已知点， 3 个未知点， 7 个测段。已知点高程 H1=5.016M H2=6.016 h1=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 号点名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 11 页 - - - - - - - - - 高程 X1、X2、X3 为参数。（2） 误差方程：V1=x1 v2=x2 v3=x1 v4=x2 v5=x2-x1+h2-h1-h5 v6=x3-x1 v7=-x3 (3)取 1M 的观测高程为单位权观测，即p=1/s；(4)求法方程： Nbbx-W=0 Nbb=b pb W=b pl (5)求的平差值 x=Nbb-1*W L=l+V V=bx-l （6）高差权函数式：k=-x1+x2 (6)求中误差：单位权中误差0，协因数阵Nbb-1.名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 11 页 - - - - - - - - - 求得中误差2、平差程序流程代码说明：h1=1.359; h2=2.009; h3=0.363; h4=1.012; h5=0.657; h6=0.238; h7=-0.595; H1=5.016 H2=6.016 h=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 ; 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 11 页 - - - - - - - - - -1 0 1 ;0 0 -1 p=diag(1./s) l=0;0;4;3;7;2;0 W=B*p*l Nbb=B*p*B x=inv(Nbb)*W V=(B*x-l) 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 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 11 页 - - - - - - - - - D=d*sqrt(q) D1=d*sqrt(Q) (3) 平差程序流程代码说明 : clc clear disp(观测高差，单位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 % 已知点高程，单位m H2=6.016 % 已知点高程，单位m 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 11 页 - - - - - - - - - h=h1 h2 h3 h4 h5 h6 h7 s=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;0 W=B*p*l Nbb=B*p*B disp(参数的解) x=inv(Nbb)*W V=(B*x-l) % 误差方程 (mm) H=h+V/1000 % 观测值的平差值Q=inv(Nbb) %观测值协因数阵名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 11 页 - - - - - - - - - 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 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 7 页，共 11 页 - - - - - - - - - 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.3846 W = -2.6061 4.0278 1.4286 Nbb = 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 8 页，共 11 页 - - - - - - - - - 2.4748 -0.4167 -0.7143 -0.4167 1.3753 0 -0.7143 0 1.0989 x = -0.2427 2.8552 1.1423 V = -0.2427 2.8552 -4.2427 -0.1448 -3.9021 -0.6151 -1.1423 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 9 页，共 11 页 - - - - - - - - - H =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 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 10 页，共 11 页 - - - - - - - - - 0.1608 0.7758 0.1045 0.3450 0.1045 1.1342 j = 19.7994 d = 2.2248 f = -1 1 0 q =0.9850 D = 2.2080 D1 =1.6208 0.8921 1.3067 0.8921 1.9597 0.7192 1.3067 0.7192 2.3694 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 11 页，共 11 页 - - - - - - - - -
