测量平差教案
第一章绪论
第一节观测误差
一、观测值中为什么存在观测误差?
观测条件对观测成果产生影响,不可避免产生观测误差。有观测就有误差的结论。
二、观测误差的计算
给出观测误差计算的纯量表达式和矩阵表达式。
三、观测误差的分类及其处理
1、分类
给出误差分类的表达式,粗差、系统误差和偶然误差的定义。结合测角、测距和水准测量的全过程,让学生分析哪些因素引起的误差属于粗差,那些哪些因素引起的误差属于系统误差,那些哪些因素引起的误差属于偶然误差。
2、处理
总结粗差、系统误差和偶然误差的处理方法,让学生举例说明测量上哪些操作是为了消除系统误差影响的,那些计算改正为了消除系统误差影响的。
四、测量平差的任务
根据一系列含有观测误差的观测值求待定量的最佳估值。
第二节测量平差学科的研究对象
研究对象为含有观测误差的各类观测值。举例说明。
第三节测量平差的简史和发展
一、测量平差理论的发展
1、经典平差理论的发展
主要介绍高斯创立最小二乘原理和马尔可夫创立高斯-马尔可夫平差模型的历史背景和过程。
2、近代平差理论的发展
主要介绍二十世纪四十年代以后出现的近代平差理论,结合导线网平差和我国南极考察、建站,重点介绍方差分量估计和秩亏网平差的理论、方法及其用途。
二、平差计算方法的发展
1、手算阶段
2、半自动平差阶段
3、全自动平差阶段
第四节测量平差的任务和内容
一、任务
讲授测量平差的基本理论和基本方法,为进一步学习和研究测量平差打下深入的基础。
二、内容
课本各章的内容。
小结:本节介绍了观测条件的定义,观测条件与观测误差的关系,观测误差的定义、处理,以及测量平差的发展概况。
第二章误差分布与精度指标
第一节正态分布
一、一维正态分布
绘一维正态分布图,列出分布函数,讲解,强调两个分布参数的含义。
二、n维正态分布讲解
绘n维正态分布图,列出分布函数,讲解,强调两个分布参数的含义。
第二节偶然误差的规律性
一、偶然误差分布
1、描述误差分布的三种方法
(1)列表法(通过实例列表讲解)
(2)绘图法(通过实例绘图讲解)
(3)密度函数法(通过实例绘图讲解)
二、偶然误差的分布特性
(1) 在一定的观测条件下,误差的绝对值不会超过一定的限值。(界限性)(2) 绝对值较小的误差比绝对值较大的误差出现的概率要大。(小误差占优性)(3) 绝对值相等的正负误差出现的概率相等。(对称性)三、两个重要概念
(1) 由偶然误差的界限性,可以依据观测条件来确定误差限值(2) 由偶然误差的对称性知观测量的期望值就是其真值。
小结:偶然误差有其统计规律,研究偶然误差的分布规律是为了更好的研究偶然误差的处理问题。
第三节衡量精度的指标;
第四节精度、准确度与精确度;第五节测量不确定度
一、精密度指标
(一)观测量的精密度指标
1、观测条件与精密度
配合误差分布曲线讲解精密度的定义和观测条件与精密度的关系。
2、几种常用的精密度指标
(1)方差与标准差
推导相应公式,给出其估值公式,讲解应用实例
(2) 极限误差
分析误差出现在某一范围内的概率的大小,给出极限误差定义公式
(3) 相对误差
给出相对精度的定义,用实例讲解其应用范围。
(4) 平均误差与或然误差
给出平均误差和或然误差的定义,讲解其在国际上应用的范围和地区,以及其与中误差的关系。
(二)观测向量的精度指标
1、n维随机向量的方差阵
导出n维随机向量的方差阵表达形式,指出该阵是对称矩阵,并讲解矩阵中各元素的含义,同时给出当n维随机向量中各随机变量不相关时的矩阵形式。
2、两随机向量的互协方差阵
导出两个随机向量互协方差阵表达形式,并讲解矩阵中各元素的含义,同时给出当维随机向量不相关时的矩阵形式。
二、准确度和精确度指标
分别给出准确度和精确度的定义,及其数值指标,绘图讲解其几何意义。三、测量不确定度
给出测量数据的不确定性、不确定度的概念,可测不确定度的计算方法,不可测不确定度的估计方法。
小结:精度指标分为精密度指标、准确度指标和精确度指标三种,观测成果的质量应用精确度指标衡量,精密度指标中的方差、极限误差、相对误差几个指标应重点掌握。
第三章协方差传播律及权
第一节数学期望的传播律
;)(C C E =
);()(X CE CX E =
);()()()(2121n n X E X E X E X X X E +++=+++ 当X i 相互独立时(i =1,2, …,n ),
)()()(),,,(2121n n X E X E X E X X X E =
第二节协方差传播律
协方差传播律是观测值(向量)与其函数(向量)之间精度传递的规律。
一误差的传递
1、线性函数误差的传递
02211...f x f x f x f Y n n ++++= n x n x x Y f f f ?++?+?=? (2121)
推导上述公式,讲解式中符号的含义 2、非线性函数误差的传递
()n x x x f Y ...
2
1
=
n x n x x Y f f f ?++?+?=? (2121)
推导上述公式,讲解式中符号的含义 3、函数向量误差的传递
Y=FX+F 0 Y=F(X) ΔY =F ΔX
讲解式中符号的含义,强调矩阵表达式与纯量表达式之间的相互表式 二、协方差的传递 1、基本公式 函数向量
Y=F(X) Z=K(X)
其误差向量为
ΔY =F ΔX ΔZ =K ΔX
则随机向量与其函数向量间的方差传递公式为
????
?????====
F
D K D
K D F D K D K D F
D F D T X
ZY
T X YZ T X Z
T
X Y 证明第一、第三式,并说明同理可证二、四式。
2独立观测量函数的方差传递
2
21222221212...n
T
X
Y f f f F
D F
σσσσ+++== 讲解式中符号的含义,说明公式应用的条件,强调公式的重要性。 3、分块向量函数向量的方差传递
?????
???????
=+Y
X Z r t r t 1
,1,1, ??????
?
???????=X
r
r YX t
r XY
r
t X t
t D D D
D D Z ,,,, 证明上式,对阵中元素加以说明,给出两向量不相关时该矩阵的形式。
通过五个典型例题的讲解说明方差-协方差传播公式的应用方法和计算中需注意的问题。
小结:协方差传播律是观测值(向量)与其函数(向量)之间精度传递的规律,用其解决观测值函数(向量)的精度评定问题。本节重点是利用协方差传播律解题的方法和步骤,以及只有一个观测值函数,且观测值之间不相关时的协方差传播公式的应用。
第三节协方差传播律的应用
1、水准测量的精度
绘制具有N 个测站的水准高差示意图,应用协方差传播公式导出高差中误差计算公式:
站σσ
N h =? 进一步导出S 公里观测高差的中误差计算公式:
k m ??σσ
S h = 举例说明公式的应用。
2、同精度独立观测值的算数平均值的精度
由算术平均值公式,应用协方差传播公式导出其中误差计算公式
N
x σ
σ=
举例说明公式的应用。 3、若干独立误差的联合影响
n z ?++?+?=? 21
222212n Z σσσσ+++=
即观测结果的方差,等于各独立误差所对应的方差之和。
4、平面控制点的点位精度
绘支导线略图,求未知点点位中误差,用两种方法求解。 解法一: (1)、列函数式 (2)线性化
(3)应用协方差传播公式计算坐标方差 (4)计算点位方差
解法二:利用纵向方差和横向方差进行计算。
小结:本节的重点内容为水准测量高差和同精度独立观测算数平均值的精度计算问题,应熟记计算公式,能熟练应用公式进行相关计算。
第四节权与定权的常用方法
一、权的定义
权是衡量各观测值在平差结果中应起作用大小的数值。
()n i i
i
P ,...,2,1220==
σ
σ
P i 为观测值L i 的权,2
0σ是可以任意选定的比例常数。 观测值的权与观测值的方差成反比。
二、单位权方差
权的作用是衡量观测值的相对精度,称其为相对精度指标。确定一组权时,只能用同一个σ0, 令σi =σ0,则得:
i
i
P ===022020
21σσσσ
上式说明20σ是单位权(权为1)观测值的方差,简称为单位权方差。凡是方差等于2
0σ的观测值,其权必等于1。权为1的观测值,称为单位权观测值。无论2
0σ取何值,权之间的比例关系不变。
举例(例〔1〕、例〔2〕)讲解。
三、测量中常用的定权方法
1、水准测量的权
(1)、用测站数定权(山地、起伏较大的丘陵)
利用用测站数计算高差中误差的公式和权的定义式导出利用测站数定权的公式。
N
C
P h =
解释式中符号的含义。 (2)、用路线长度定权(平地)
利用用路线长度计算高差中误差的公式和权的定义式导出利用路线长度定权的公式。
S
C
P h =
解释式中符号的含义。
举例(例〔3〕、例〔4〕、例〔5〕)讲解。 2、距离量测的权
距离长度可通过钢尺丈量或测距仪测距得到。下面分别讨论两种情况下的定权方法。 1) 钢尺量距的权
S
C P
S
=
解释式中符号的含义。 2) 测距仪测距的权
S
S
P =022σσ
610-标称标称+?=S S σσσ
解释式中符号的含义。 3、等精度观测算术平均值的权
利用等精度独立观测值算术平均值的方差计算公式和权的定义式导出利用观测次数定权的公式
C
P
n S
=
说明公式中符号的含义。
小结:权是用来衡量观测成果的相对精度的,单位权方差可以根据计算方便任意选定,但观测值之间的比例关系不变。水准测量的权与测站数或路线长度成反比;钢尺量测的权与距离长度成反比,光电测距的权用定义式计算,其中测距方差由固定误差和比例误差两项组成;等精度算术平均值的权与观测次数成正比。应熟记定权公式,明确公式中各符号的含义,掌握利用公式解题的方法。
第五节协因数和协因数传播律
一、协因数
定义协因数
n i p i
i
ii
Q ,...,2,1/120
2
===σ
σ
权可表示为
i
ii
P Q =1
方差和标准差可表式为
Q
Q ii
i
ii i
σσσσ
2
02,
=?=
二、协因数阵
1、n 维随机向量X 的协因数阵 定义互协因数:
ij
ij
Q =σσ
02
利用方差协方差与协因数弧协因数的关系导出协因数阵
??
???
????
????
?==
Q Q
Q
Q Q Q D
Q
nn n n
X
X
对称............
1
222
1121120
σ
上式矩阵中,ij
ji
Q Q
=。当Q ij =0(i ≠j)时,则X i 和X j 互不相关。
2、分块向量的协因数阵
??
????=Y X Z
???
?????=Q
Q Q
Q Q
Y
YX
XY X X
式中,Q X 、Q Y 分别为X 、Y 向量的自协因数阵,而Q XY 、Q YX 分别为X 向量关于Y 向量的互协因数阵,Q XY 与Q YX 互为转置。当Q XY 等于零时,表示X 、Y 互不相关。
三、权阵
??????
???
???
??=??
??
??????????==--P P
P P
P P P P
P Q Q
Q Q Q Q Q Q
Q
Q P
nn n n n n
nn m n n n ...
......
.....................
.........
2
1222
21112111
2
1
222
21
112111 观测值的权一般要通过对权阵求逆得到协因数阵,再利用权与协因数的倒数关系求权。当权阵为对角阵时,
P
P ii
i
=。
举例(例〔1〕、例〔2〕)讲解、分析
四协因数传播律
将协方差传播公式乘以
20
1
σ,并顾及Y Y D Q 20
1
σ=
,X X D Q 20
1σ=,即可得到观测向量X 与其函数向量
Y 、Z 之间的协因数传播公式。
列出相应公式,以及只有一个函数,且观测值之间不相关时的协因数传播公式。 举例(例〔3〕、例〔4〕)讲解、分析
小结:权与协因数互为倒数关系,权阵与协因数阵互为逆阵关系,一般情况下给了观测值的权阵求观测值的权要先求权阵的逆阵得到其协因数阵,再利用权与协因数的关系求权;协因数传播律与协方差传播律公式相仿,只记住其中一套公式,再记住协因数阵与协方差阵的关系即可。
第六节由真误差计算中误差及其实际应用
一、利用不同精度的真误差计算单位权中误差的基本公式
利用协因数传播律导出利用不同精度的真误差计算单位权中误差的基本公式
σ=+
P n
??
二、由真误差计算中误差的实际应用 1、由三角形闭合差计算测角中误差
利用协因数传播律导出由三角形闭合差计算测角中误差的公式
[]n
WW 3?+=σ
说明公式的不严密性。 2、利用双观测列之差求中误差 (1)求单位权中误差
利用协因数传播律导出利用双观测列之差求单位权中误差的公式 不等精度观测
2σ=+
P dd n
等精度观测
2σ=+
dd n
说明公式中符号的含义。
(2)求双观测列单次观测的中误差
i
L L P i
i
10
σσσ==''' (3)求双观测列平均值的中误差 利用协因数传播律导出相应公式 不等精度观测
i
L P i
210
σσ= 等精度观测
i
L
_σσ=+0
12
水准测量双观测平差应用例题
小结:本节重点是利用双观测之差计算中误差的公式及其应用,该公式在测量中应用广泛,应重点掌握。
第七节系统误差的传播
一、观测值的系统误差与综合误差的方差 1、观测值的系统误差——偏差 导出偏差表达公式
2、观测值的综合误差方差——可靠性
如果系统误差部分是偶然中误差部分的三分之一或更小时,则可将系统误差的影响忽略不计。 二、系统误差的传播 导出传播公式
三、系统误差与偶然误差的联合传播 导出传播公式
小结:了解系统误差的传播规律。
第四章 平差数学模型与最小二乘原理
第一节测量平差概述
一、测量控制网简介
1.高程控制网(水准网或三角高程网) 包括闭合水准网和符合水准网。
绘出三组不同网形的水准网。
网中元素:已知高程点,未知高程点和观测高差。
2.平面控制网
(1)三角网
根据观测量的不同,三角网分为测角三角网、测边三角网和边角同测三角网。
1)测角三角网
包括独立三角网和符合三角网。
绘出一组不同网形的三角网。
网中元素:已知点,未知点和观测角度。
2)测边三角网
包括独立测边网和符合测边网。
绘出一组不同网形的测边网。
网中元素:已知点,未知点和观测边长。
3)边角三角网
包括独立边角网和符合边角网。
绘出一组不同网形的边角网。
网中元素:已知点,未知点,观测角度和边长。
(2)导线网
包括独立导线网和符合导线网。
绘出一组不同网形的边角网。
网中元素:已知点,未知点,观测角度和边长。
还有三维网、GPS控制网、航测控制网、工程专用网等将在后续相应课程中介绍。
二、必要起算数据
确定几何(物理)图形的位置所必须具有的已知数据
①水准网(三角高程网):一个已知点高程
②测站平差:一个已知方位
③测角网:一个已知点坐标,一个相邻已知方位,一个相邻已知边长或两个相邻点坐标。
④测边网和边角网:一个已知点坐标,一个相邻已知方位。
各种控制网中少于等于必要起算数据的控制网成为独立网,多于必要起算数据的控制网成为非独立网或附合网。
三、必要观测及其数目的确定
确定几何、物理模型的形状、大小所必须进行的观测称为必要观测,其符号用符号t表示。
高程网: t=p-q-1
测站平差: t=p-q-1必要起算数据
测角网: t=2p-q-4
测边网和边角网:t=2p-q-3
P:总点数或总方向数(测站平差);q
必要起算数据之外的起算数据
四、多余观测及其数目的确定
必要观测之外的观测称为多余观测,其数目用符号r 表示 多余观测数=观测总数-必要观测数(r=n-t )
五、必要观测和多余观测数目计算练习
计算图3-1至图3-7的必要观测数和多余观测数。
小结:本节介绍了测量控制网的类型,和各类控制网中应具备的必要起算元素,必要观测元素,应重点掌握必要观测元素数和多余观测元素数的计算。
第二节函数模型
1、条件平差法
()
0~
1
,=L F r
2、间接(参数)平差法
??
? ??=1,1,~~t n X F L 3、附有参数的条件平差法
()
0~
,~1
,=X L F c
4、附有限制条件的间接(参数)平差法
()
??
???=??? ??=0~
~~1,1,1,X X F L s u n φ 用简单控制网图形举例说明。
第三节函数模型的线性化
设
??
? ??=1,1,1,~,~u n c X L F F 用泰勒公式导出F 的线性形式为
()
x B A X L F F ~,0+?+=
根据上述函数模型线性化过程,可将各种平差方法的函数模型线性化
1、 条件平差法
0=-?W A 式中 L
L F A ~??=,()L F W -=
2、间接平差法
l x B -=?~
式中 0
~X X F B ??=,()0
X
F L l -=,
3、附有参数的条件平差法
0~=-+?W x B A
式中 (
)0
,X L F W -=,
4、附有限制条件的间接平差法
~~=--=?x W x c l
x B
式中
~
,X
X
u
s c ?Φ
?=,()
X W x Φ-=
第四节测量平差的数学模型
1、各种平差方法的随机模型
n
n n
n n
n P Q D ,12
0,20,-==σσ
2、各种平差方法的数学模型
各种平差方法函数模型的线性形式分别与平差的随机模型联立,即为相应平差方法的数学模型。 小结:本次课所讲的各种平差方法的函数模型均能建立各观测值之间的函数关系式,正确建立这种关系式,是正确求得观测值最可靠结果的前提。
第五节 参数估计与最小二乘原理
一、引例
已知平面三角形三内角应满足
1
2
3
1800~
~~L L L ++-=
或
03
2
1
=-++?
??W
式中
()L L L W 1110180++-=
上方程中有三个未知数,是相容方程,只能在某一准则下求得式中未知数的估值。 二、最小二乘准则:
min 1=??=Φ-D T
顾及方差阵与权阵的关系,并用Δ的估值V 代替Δ又可得
Φ==T
V PV min
观测量真值向量的估值公式为:
V L L
+=? 式中L
?称为观测向量的“最或然值”向量或“观测值的平差值”向量;V 称为改正数向量。 三、最小二乘估计
根据最小二乘准则进行的估计称为最小二乘估计,按此准则求得一组估值的过程,称为最小二乘平差,由此而得到的一组估值是满足方程的唯一解。
如果方差阵D 和权阵P 是非对角阵,则表示观测值是相关的,按此准则进行的平差即称为相关观测平差。如果是对角阵,则表示观测值是彼此不相关的,此时称为独立观测平差。
当观测值不相关,即P 为对角阵时,则有
Φ==∑=+++==T
i
i n
i
n
n
V PV P V P V P V P V 1
21
1
22
2
22 (i)
当观测值不相关,并为等精度,即P=I 时,则有:
Φ==∑=+++==T
i n
i
n
V V V V V V 1
21
22
22 (i)
小结:最小二乘原理是测量平差的基本原理,按最小二乘准则求得的观测量及其函数的结果是最可靠的结果,后续所讲所有平差方法均按此准则求解。
第五章 条件平差
第一节条件平差原理
一、条件方程和改正数条件方程
列出用观测值真值和真误差表示的条件平差函数模型
导出用按最小二乘准则求得的观测值平差值和观测值改正数表示的条件平差的函数模型
()
1
,1
,0?r r L
F =——条件方程 01
,1
,,=-r n n r W V A ——改正数条件方程
()L F W -=——改正数条件方程常数项(闭合差)计算式
举例(单三角形函数模型的建立)
二、条件方程的纯量表达式和矩阵表达式
r 个条件方程的纯量表达式:
()()(
)
???
???
?===0?,,?,?0?,,?,?0?,,?,?21212211n r n
n L L L F L L L F L L L F 线性化后得改正数条件方程
???
???
?=-+++=-+++=-+++000221122112211r n n b n n a n n w v r v r v r w v b v b v b w v a v a v a
其中
()()()???
???
?-=-=-=n r r n b n a L L L F w L L L F w L L L F w ,,,,,,,,,21212211
令
????????????=n n n r r r b b b a a a A 2121
21
, ?
?
????
??????=r r w w w W 211,, ?????
???????=n n v v v V 211, 则改正数条件方程及其闭合差计算的矩阵表达式分别为
0=-W AV
()L F W -=
三、基础方程
按求函数极值的拉格朗日乘数法,设其乘数为()T r b a r k k k K =1
,,称为联系数向量。组成函数
()W AV K PV V T T --=Φ2,
对其求导整理得改正数V 的计算公式
K QA K A P V T T ==-1——改正数方程
当P 为对角阵时,改正数方程的纯量形式为
()n i k r k b k a v r i b i a i p i i
,,2,1,
1
=+++=
改正数条件方程与改正数方程联立,称为条件平差的基础方程。此时,方程的个数与未知数的个数相
同,方程有唯一解。
四、基础方程的解
将改正数方程代入改正数条件方程,得
0=-W K AQA T ,
令
T T aa A AP AQA N 1-==,
得
0=-W K N aa ——联系数法方程
秩()()
()r A R AQA R N R T aa ===,即aa N 是个r 阶的满秩方阵,由此解出
W N K aa 1
-=
当P 为对角阵时,法方程的纯量形式为
[][][][][][][][][]
???
?
??
?=-+++=-+++=-+++000r r p rr b p br a p ar b
r
p
br b
p
bb a
p
ab a
r
p
ar b
p
ab a
p
aa w k k k w k k k w k k k 解出K ,将其代入改正数方程,求出改正数V ,在按V L L
+=?可求得平差值。 五、条件平差步骤及示例
用具有两个条件的符合水准网为例讲解。
小结:本节应熟记条件方程,改正数条件方程,改正数条件方程闭合差计算式,法方程,改正数方程的表达形式,掌握用条件平差法平差的方法、步骤。
第二节条件方程
一、水准网(同§5.1中所述,略) 二、测角网
1.单三角形(同§5.1中所述,略)
2.中心多边形
以中心三边形为例,画出示意图,列出其条件方程和改正数条件方程的一般表达式。 重点讲解极条件的列立方法和规律。
举例(中心三边形实例)列条件方程和改正数条件方程。 3、大地四边形
画出示意图,列出其条件方程和改正数条件方程的表达式。 重点讲解极条件的列立方法和规律。
举例
上图中,若以对角线交点为极列极条件,其极条件闭合差超限,说明角度观测存在问题,如何返工?先让让学生回答,然后教师讲解。
三、测边网
1.中心多边形
画出测边中心三边形示意图。
(1)列出以反算角表示的条件方程和改正数条件方程
(2)建立反算角改正数与边改正数之间的关系
(3)导出以边改正数表示的条件方程
2.大地四边形
画出测边大地四边形示意图。
(1)列出以反算角表示的条件方程和改正数条件方程
(2)建立反算角改正数与边改正数之间的关系
(3)导出以边改正数表示的条件方程
四、边角网
如图,t=2p-q-3=8-3-3=2,r=n-t=8-2=6
应列出6个条件方程
条件分析:内角和条件 2个
正弦条件 2个
固定角条件 1个
规定边条件 1个
边角网条件方程列立例题讲解分析。
小结:条件方程列立,首先应能正确确定应列的条件数目,保证方程之间不相关,其次应能分析条件类型,最后应掌握各类方程的列立规律,正确列出条件方程。
第三节 精度评定
一、单位权方差估值计算
r
?2
0PV
V T =ο
PV V T 的计算:
1、()
权阵为对角阵时=2n 222211P P P n T V V V PV V +++ 2、K N K K AQPQA K K QA P K QA PV V aa T T T T T T T ===)()( 3、K W K A V K QA P V PV V T T T T T T ===)( 二、协因数阵 设
[
]
T T
T T T
T L
V K W L Z ?= 列出各分块向量解的表达式及其微分式,利用协因数传播律导出各量的协因数阵和各量之间的互协因数阵的结果列于相应表中
表中L
?与V 、W 、K 的互协因数阵为零,说明L ?与V 、W 、K 统计不相关 证明:表中VV Q 、L L Q ??的计算表达式。 三、观测值平差值的精度评定
L L L L Q D ?
?2
0???σ= 四、平差值函数的精度评定
1.平差值函数表达式及其协因数计算 列出平差值函数表达式
)
,,(=n 21?,???L L L f ? 按泰勒公式展开,并按协因数传播律导出平差值函数协因数的计算公式
()AQf N AQf Qf f Q aa T
T 1
--=??
f i (i =1,2,…,n)为偏导数值。
2.权函数式
n L L
n L L L L L d
L f L d L f L d L f d ????????2?21?1===?
??? ????++???? ????+???? ???? =?
n
n L d f L d f L d f ???2211+++ =——权函数式 3.平差值函数的方差
??
??σ??2
0???Q D = 小结:本节主要介绍了利用改正数计算单位权中误差的公式,各种平差量协因数和互协因数及方差协方差的计算,平差值函数式和权函数式的列立方法,平差值函数协因数和互协因数及方差协方差的计算方法,应重点掌握。
第六章 附有参数的条件平差
一、概述
设ABD X ∠=~
,又可列出1个极条件和一个固定边条件
极条件为(以A 点为极):(
)(
)(
)
1~sin ~~sin ~~sin ~
sin ~sin ~~sin 5
8696
75=+-+L L L X L L X L L
固定边条件为(由AC 边推算到AB 边):()
X
L L L L S S AC AB ~
sin ~sin ~
~sin ~sin 38
62+= 或
()
1~
sin ~sin ~
~sin ~sin 38
62=+X
L S L L L S AB AC 由于选了一个参数,增加了一个条件,一般情况下,若选了u 个参数,则条件方程的数目为c=r+u.从以上5 个方程出发进行平差,就是附有参数的条件平差方法。
二、基础方程
观测量L ~和X ~的最佳估值V L L +=?,x
X X ??0+=,用奇表示的附有参数的条件平差函数模型为
()
1
,1
,0?,?r r X L
F =——条件方程 或
0?1
,1
,,1
,,=-+c u u c n n c W x
B V A ——改正数条件方程 ()
0,X L F W -=——改正数条件方程常数项(闭合差)计算式
按求函数极值的拉格朗日乘数法,设其乘数为()T
r b a r k k k K =1
,,称为联系数向量。组成函数
()W x
B AV K PV V T T -+-=Φ?2, 将Φ对V 和x
?分别求一阶导数,并令其为零,导出改正数V 的计算公式 K QA K A P V T T ==-1——改正数方程
附有参数的条件平差的基础方程为:
0?1
,1
,,1
,,=-+c u u c n n c W x
B V A K QA K A P V T T ==-1
0=K B T
方程的个数与未知数的个数相同,方程有唯一解。
三、基础方程的解
将改正数方程代入改正数条件方程,并令T
T aa A AP AQA N 1-==,则得法方程
0?=-+W x
B K N aa 0=K B T ——法方程
秩()()
()c A R AQA R N R T aa ===,即aa N 是个c 阶的满秩方阵,顾及B
N B N aa T bb 1
-=,
由法方程可解出
,?1
1W N B N x aa T bb --= ()x B W N QA V aa
T ?1
-=-, 四、精度评定
(一)、单位权方差估值计算
u
-c r ?2
0PV
V PV V T T ==ο