原创力文档- 1、本文档共25页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
MATLAB条件平差及间接平差
测量平差程序设计;一、条件平差基本原理;1、平差值条件方程：;2、条件方程：;3、改正数方程：;4、法方程：;按条件平差求平差值计算步骤;例 《误差理论与测量平差基础》P77;clc
disp(‘条件平差示例5-1’)
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
;例 《误差理论与测量平差基础》P78;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));作业
文档评论（0）