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

测量程序设计一、条件平差基本原理1、平差值条件方程：2、条件方程：3、改正数方程：4、法方程：按条件平差求平差值计算步骤例《误差理论与测量平差基础》P74clc
Disp(‘条件平差示例２’)
Disp(‘三角形内角观测值’)
L1=[421220]
L2=[7899]
L3=[593840]
L=[L1;L2;L3]
Disp(‘将角度单位由度分秒转换为弧度’)
LL=dms2rad(mat2dms(L))

A=[111]
w=sum(LL(:))-pi
w=dms2mat(rad2dms(w))
P=eye(3);

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

V=A'*Ka
例《误差理论与测量平差基础》P75clc
clear
h1=-1.004;
h2=1.516;
h3=2.512;
h4=1.520;
HA=12.013
HB=10.013
h=[h1h2h3h4]'
s1=2;
s2=1;
s3=2;
s4=1.5;
s=[s1s2s3s4]'

ifH(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))