：hFTX
Qhh

FTQXXF

F
T
N
1BB
F
5
f误差理论与测量平差基础课程设计报告
四、程序流程图
根据题目列出条件方程并写成误差方程的形式VBxl
↓
根据误差方程得到矩阵B、l进而写出BT
↓
运用MatLab程序语言定义出权P
↓
用MatLab程序求出参数的改正数xNBB1W，其中NBBBTPB、WBTPl
↓
根据MatLab程序语言求Bx，进而由VBxl写出各观测值的改正数
↓
根据LLV求出各观测值的平差值
↓
检验所求各值是否正确，若无误则往下进行，反之检查各步骤查出错误并改正
↓
由程序计算VTP进而求出VTPV，求单位权中误差，再根据权函数式、协因数传播定律评定各观测值及所求高程的精度
6
f误差理论与测量平差基础课程设计报告
五．程序源代码及说明
121218上午922121218上午926AJAHD121218上午952dispHa5016Hb6016H1359200903631012065702380595B100010100010110101001l0043720S11172327241426Pi
vdiagSx0i
vBPBBPlVBx0lH0HV1000
6t2delsqrtVPV
tBBBPBQ34110i
vBB1103100i
vBB100Q33100i
vBB100Q44010i
vBB010m4delsqrtQ44m3delsqrtQ33m34delsqrtQ34
7
f误差理论与测量平差基础课程设计报告
六．执行调试，得出计算结果
disp水准网间接平差事例水准网间接平差事例Ha50160已知点高程，单位mHa
50160Hb6016已知点高程，单位mHb
60160disp观测高差，单位m观测高差，单位mH1359200903631012065702380595H
13590200900363010120065700238005950disp系数矩阵B系数矩阵BB100010100010110101001B
100010100010110101001l0043720常数项l0043720S11172327241426是线路长度S11000170002300027000240001400026000
8
f误差理论与测量平差基础课程设计报告
Pi
vdiagS定义观测值的权
P
09091
0
0
0
0
0058820
0
0
0
004348
0
0
0
0
003704
0
0
0
0
004167
0
0
0
0
0
0
0
0
0
0
00000071430
00000003846
disp参数的解参数的解x0i
vBPBBPlx0
024272855211423disp误差Vmm各待定点的高程平差值L1（m）
误差Vmm各待定点的高程平差值L1（m）
VBxl误差方程mm
V02427285524242701448390210615111423
H0HV1000观测值的平差值H0
13588201190358810119065310237405961disp精度评定精度评定
6观测值的个数t2必要观测数
9
f误差理论与测量平差基础课程设计报告
delsqrtVPV
t
del
22r