matlab平差结果转换,用MATLAB解决条件平差及间接平差.ppt

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

测量程序设计;一、条件平差基本原理;1、平差值条件方程：;2、条件方程：;3、改正数方程：;4、法方程：;按条件平差求平差值计算步骤;例 《误差理论与测量平差基础》P74;clcDisp(‘条件平差示例２’)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(:)) - piw = dms2mat(rad2dms(w))P = eye(3);

Naa = A*inv(P)*A'Ka = -inv(Naa)*w

V = A'*Ka;例 《误差理论与测量平差基础》P75;clcclearh1 = -1.004;h2 = 1.516;h3 = 2.512;h4 = 1.520;HA = 12.013HB = 10.013h = [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(‘检核错误')enddisp(‘平差后的高程值')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*ldisp(‘误差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))