实验一matlab完成水准网平差

实验一matlab完成水准网平差
实验一 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号点高程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、平差程序流程代码说明：
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 ; -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
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
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)。