用MATLAB解决条件平差及间接平差.pptVIP

下载本文档

1508
0
约 24页
2017-04-19 发布于四川
举报
版权申诉

用MATLAB解决条件平差及间接平差.ppt

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

用MATLAB解决条件平差及间接平差

测量程序设计;一、条件平差基本原理;1、平差值条件方程：;2、条件方程：;3、改正数方程：;4、法方程：;按条件平差求平差值计算步骤;例《误差理论与测量平差基础》P74;clc Disp(‘条件平差示例２’) Disp(‘三角形内角观测值’) L1 = [42 12 20] L2 = [78 9 9] L3 = [59 38 40] L = [L1; L2; L3] Disp(‘将角度单位由度分秒转换为弧度’) LL = dms2rad(mat2dms(L)) ; A = [1 1 1] w = sum(LL(:)) - pi w = dms2mat(rad2dms(w)) P = eye(3); Naa = A*inv(P)*A Ka = -inv(Naa)*w V = A*Ka ;例《误差理论与测量平差基础》P75;clc clear h1 = -1.004; h2 = 1.516; h3 = 2.512; h4 = 1.520; HA = 12.013 HB = 10.013 h = [h1 h2 h3 h4] s1 = 2; s2 = 1; s3 = 2; s4 = 1.5; s = [s1 s2 s3 s4]; if H(1,1)+H(2,1)-H(3,1)+HA-HB==0 H(2,1)-H(4,1)==0 disp(‘检核正确) else disp(‘检核错误) end disp(‘平差后的高程值) HC = HA + H(1,1) HD = HA + H(1,1) + H(4,1) ;二、间接???差的基本原理;设误差Δ和参数X的估计值分别为V 和;由最小二乘准则，所求参数的改正数应该满足;设;精度评定;设参数的函数为;中误差为;二、水准网间接平差程序设计;线路编号;程序代码如下：;S = [1.1, 1.7, 2.3, 2.7, 2.4, 4.0] P = C./S % 定义观测值的权， P = diag(P) % 定义权阵 disp(‘参数的解’) x = inv(B’*P*B)*B’*P*l disp(‘误差V(mm), 各待定点的高程平差值L1（m）’) V = B*x - l % 误差方程(mm) L1 = L + V/1000 % 观测值的平差值， disp(‘精度评定’) n = 6; % 观测值的个数 t = 2; % 必要观测数 delta = sqrt(V’*P*V/(n – t))