测量数据处理理论与方法-0

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
武汉大学测绘学院 黄海兰
1
2015-10-25
武汉大学测绘学院 孙 海燕
课程教学大纲
9 教学目的与要求 本课程通过对现代测量数据处理理论和方法的系统介绍,使
学生了解目前的研究方向,明确目标,通过该课程,进一步掌 握专业知识,逐步培养学生进行数据处理的能力,为进一步后 续专业课程的学习和科研项目的研究打下良好的基础。 9 课程目标 使学生不仅具有扎实的理论基础,而且具有开阔的思路和较强 的解决实际问题的能力。
Φ( X~ ) = 0
Φ( Xˆ ) = 0
武汉大学测绘学院
黄海兰
绪论
3、经典测量平差的数学模型(间接平差)
1)数学模型
L~ = L + Δ E(Δ) = 0 D(Δ) = σ 02Q
P = Q−1
V T PV = min V = Bxˆ − L
2)解算
B T PB xˆ = B T PL
xˆ = ( B T PB ) −1 B T PL
观测值个数 n = 5
必要观测数 t = 3
多余观测数 r = 2
参数个数 u = t = 3
H B = Xˆ1 H C = Xˆ 2
H D = Xˆ 3
方程个数 c
c=u+r =n=5
Lˆ1 =Xˆ 1 − H A Lˆ2 =Xˆ 2 − H A
Lˆ3 =Xˆ 3 − H A Lˆ4 = − Xˆ1 + Xˆ 2 Lˆ5 =Xˆ 2 − Xˆ 3
必要观测数
t ==44 − 1 = 3
S P1
(
Xˆ 1
,
Yˆ 1
)
L1 L6
参数个数 u = 4 u > t
非独立参数个数 s = u − t = 1
限制条件方程:
( Xˆ 2 − Xˆ1)2 +(Yˆ2 −Yˆ1)2 − S = 0
L2 A
方程个数 c = u + r = s + t + r = s + n
c=u+r =t+r =n
Lˆ = F ( Xˆ )
Lˆ n×1
=
nB×t
Xˆ + t×1
d
n×1
V
n×1
=
nB×t
xˆ + t×1
l
n×1
(l
=
L

0 ( BX
+
d ))
待求量: n + t
Lˆ Xˆ
-17-
武汉大学测绘学院 黄海兰
附有限制条件的间接平差
已知点:A、B 坐标及边长S, P1、 P2 待定点:
P2 ( Xˆ 2 , Yˆ2 ) L5
L3 L4 B
观测方程个数:n
参数间的限制条件方程数:s
Lˆ = F ( Xˆ ) Φ( Xˆ ) = 0
附有限制条件的间接平差
-18-
武汉大学测绘学院 黄海兰
9
2015-10-25
平差函数模型的比较
函数模型
观测数 必要观测数 多余观测数 参数个数
方程数 待求量数 方程形式
响,以保证测量成果的精度。
三、 平差方法
经典平差:条件平- 差-、-1间5 接-9平-差3 0、+附+有参+15数+的21条+2件7 平差、附有限 制条件的间接平2差7 和21附-18有-1限2 制-6条件3的+ 条9+1件2 +平18差+2等4 五种平差方法。
6武汉大学测绘学院 黄海兰
7
2015-10-25
估计方法,并讨论其实际应用;
¾ 研究讨论动态线性系统,建立动态系统的卡尔曼滤波算法;
武汉大学测绘学院 孙
海燕 武汉大学测绘学院 黄海兰
课程教学大纲
9 教学内容 ¾ 研究稳健估计的基本理论,讨论常用的几种抗差最小二乘法
及其应用;介绍几种特殊问题的估计方法,如有偏估计、整 体最小二乘估计、半参数估计以及非线性最小二乘估计等基 本思想。
⎡ ⎢ ⎢
Lˆ1 Lˆ2
⎢ ⎢
Lˆ3
⎢ ⎢
Lˆ4
⎤ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥
=
⎡1 ⎢0 ⎢0 ⎢−1 ⎢⎣ 0
0 1 0 1 1
0⎤ 0⎥ 1⎥ 0⎥ −1⎥⎦
⎡ ⎢ ⎢
Xˆ1
⎤ ⎥ ⎥
⎢ ⎢ ⎢

2
⎥ ⎥ ⎥


⎢ ⎢⎣

3
⎥ ⎥⎦
+
⎡−H ⎢−H ⎢−H ⎢0 ⎢⎣ 0
A A A
⎤ ⎥ ⎥ ⎥ ⎥⎦
f (X) = 0
X =[x1 x2 ] " xm T
m>6
武汉大学测绘学院 黄海兰
绪论
平差问题中涉及的量分为:
起算数据
观测值: 参数:
L~ = L + Δ X~
Lˆ = L + V Xˆ = X 0 + xˆ
函数模型: f ( L~ ) = 0 f ( L~ , X~ ) = 0 Φ( X~) = 0
f (Lˆ) = 0
f (Lˆ, Xˆ ) = 0 Φ ( Xˆ ) = 0
观测值个数 n ,必要观测数 t ,参数个数 u
武汉大学测绘学院 黄海兰
11
2015-10-25
绪论
2、函数模型的分类
1)条件平差(u = 0 )
f (L~) = f (L + Δ) = 0 f (Lˆ) = f (L +V ) = 0
1、广义测量平差,崔希璋等,2000 2、测量数据测量理论与方法,邱卫宁等,2008 3、Parameter estimation and hypothesis testing in linear models, Koch, 1999 4、Adjustment theory: an introduction, Teunissen, 2000
武汉大学测绘学院 孙
海燕 武汉大学测绘学院 黄海兰 3
2015-10-25
课程教学大纲
9 教学内容
第一章 估计方法和广义测量平差原理
1.1 概述 1.2 多维正态分布 1.3 极大似然估计 1.4 最小二乘估计 1.5 极大验后估计 1.6 最小方差估计 1.7 线性最小方差估计 1.9 广义测量平差原理
⎢⎣ Lˆ5 ⎥⎦

B Xˆ d
观测方程:
Lˆ = B Xˆ + d
5×1 5×3 3×1 5×1
将Lˆ =
L + V 和Xˆ
=
0 X
+ xˆ代入:
误差方程:
V = B xˆ + l
5×1 5×3 3×1 5×1
(l = L − ( BX 0 + d ))
间接平差
参数个数 u = t
参数独立 观测方程个数:
9 观测误差:指待观测量的真值与其观测值之差。 不论观测条件如何,测量误差总是不可避免的。
9 多余观测:为了检查观测值中是否存在错误,并提高观测成果 的精度。
9 测量平差:依据某种最优化准则,由一系列带有观测误差的测 量数据,求定未知量的最佳估值及精度的理论和方法。
武汉大学测绘学院 黄海兰
回顾:测量平差的基本概念
条件平差
n t r=n-t 0
附有参数的 条件平差
n t r=n-t 0<u<t且独立
间接平差 附有限制条件
的间ቤተ መጻሕፍቲ ባይዱ平差
n t r=n-t u=t且独立
n
t
r=n-t U>t且包含t个
独立
r n F (Lˆ ) = 0
AΔ +W = 0
r+u
n+u
F (Lˆ Xˆ ) = 0
AΔ+BX~+W =0
r+u=n n+u Lˆ = F ( Xˆ ) Δ = BX~ − l
武汉大学测绘学院 孙
海燕 武汉大学测绘学院
黄海兰
课程教学大纲
9 教学内容 第二章 最小二乘平差的统一理论和方法 2.1 概述 2.2 秩亏自由网平差 2.4 极大验后滤波与推估 2.5 最小二乘配置 2.7 随机模型具有奇异协因数阵的平差
武汉大学测绘学院 孙
海燕 武汉大学测绘学院 黄海兰
4
课程教学大纲
武汉大学测绘学院 孙
海燕 武汉大学测绘学院 黄海兰
2
2015-10-25
9 教学内容
课程教学大纲
¾ 针对经典测量平差的局限性,系统地研究极大验后估计等若
干估计方法,构造广义测量平差原理;
¾ 讨论并突破经典平差理论及应用上限制,完成最小二乘平差
的统一理论与方法建立;
¾ 讨论测量平差的随机模型验后估计的赫尔默特法及二次无偏
sin sin(Lˆ1
Xˆ sin(Lˆ3 + Lˆ5 ) sin Lˆ1 + Lˆ4 ) sin(Lˆ6 − Xˆ ) sin Lˆ3
−1 =
0
Lˆ7 + Xˆ + α BD − α BC = 0
武汉大学测绘学院 黄海兰
8
2015-10-25
间接平差
1
A
B Xˆ1
4
2
Xˆ 2
C
3
5
D Xˆ 3
9 教学内容 第三章 平差随机模型的验后估计 3.1概述 3.2赫尔默特方差分量估计 3.3 方差-协方差分量估计 3.4 二次无偏估计法
2015-10-25
武汉大学测绘学院 孙
海燕 武汉大学测绘学院 黄海兰
课程教学大纲
9 教学内容 第四章 动态线性系统的卡尔曼滤波 4.1连续线性系统的数学模型 4.2离散线性系统的数学模型 4.3 离散线性系统的卡尔曼滤波 4.4 动态测量系统的卡尔曼滤波
二、 误差理论
研究内容包括:误差分布、精度指标、误差估计、误差检验、 误差分析以及误差预测和控制。 9 误差分布:服从正态分布(假定系统误差已经通过某种手段
得以消除,而且不存在粗差,其数学期望(真值)为零)。 9 精度指标:方差 9 误差估计:总是与平差参数估计同时进行,而且依附于平差
参数估计之中。 9 误差检验:目的是要在平差问题中排除系统误差和粗差的影
条件平差
例1水准网如右图,D为已知点,观测值及其权阵如下:
L = (0.023 1.114 1.142 0.078 0.099 1.216)T P = diag (1 1 1 2 .5 2 .5 2 .5 )
解:(1)列出条件方程
AV +W = 0,
3,6 6,1 3,1 3,1
W = AL + A0
隐含的假设: | Q |≠ 0 、| B T PB |≠ 0 ⇒ rank ( B ) = t
在什么情况上述假设无法满足?
武汉大学测绘学院 黄海兰
12
2015-10-25
绪论
4、针对经典平差模型应想到的问题
1)为什么取 P = Q−1
2)随机模型
⎧ E(Δ) = 0
⎨ ⎩
D(Δ)
=
σ
02Q
⇒ Δ 为偶然误差 ,即
2015-10-25
测量数据处理理论与方法
(现代测量数据处理理论)
黄海兰
武汉大学测绘学院
2015年1武0汉月大学1测1绘日学院 孙
海燕
课程教学大纲
9 学时安排 36学时共12次课 (6-17周,7-9节,1-328) 9 前导课程 误差理论与平差基础、高等测量平差、概率论、线性代数 (矩阵论)、高等数学 9 运用 大地测量、导航、 GNSS、摄影测量与遥感等领域 9 参考文献
φ = fL + f0
Qφφ = fQf T
σˆφ2φ = σ 02Qφφ
武汉大学测绘学院 黄海兰
10
2015-10-25
绪论
一、经典测量平差
1、函数模型
例:如图,与三角形有关的量有
1)三个内角(3)
2)三条边(3)
3)三个点的坐标(6)
4)三条边的正、反方位角(6) 合计 :18 个量 确定三角形形状、大小、位置需要6个量 即18个量中只有6个函数独立的量
9 教学内容
课程教学大纲
第六章 几种特殊问题的估计方法
6.1 有偏估计
6.2 整体最小二乘估计
6.3 半参数估计
6.4 非线性最小二乘估计 9 考核方式
闭卷考试
武汉大学测绘学院 孙
海燕 武汉大学测绘学院
黄海兰
6
2015-10-25
回顾:测量平差的基本概念
一、 测量平差问题
9 观测:指用一定的仪器、工具、传感器或其他手段在一定空间 获取反映地球及其他实体与空间分布有关信息的过程和数据。
2)附有参数的条件 平差 (0 < u < t 参数独立)
f ( L~ , X~ ) = 0
f ( Lˆ , Xˆ ) = 0
3)间接平差(u = t 参数独立)
L~ = f ( X~ )
Lˆ = f ( Xˆ )
4)附有限制条件的间接平差(u > t )
L~ = f ( X~ )
Lˆ = f ( Xˆ )
⎡1 0 0 1 −1 0 ⎤
⎡2⎤
A = ⎢⎢0 1 0 0 1 −1⎥⎥ , W = ⎢⎢− 3⎥⎥ ,
⎢⎣0 0 −1 −1 0 1 ⎥⎦
⎢⎣− 4⎥⎦
-15-
武汉大学测绘学院 黄海兰
附有参数的条件平差
F (L~ X~) = 0
Lˆ1 + Lˆ2 + Lˆ3 − 180 = 0
x
Lˆ4 + Lˆ5 + Lˆ6 −180 = 0
r+u=n+s
n+u Lˆ = F ( Xˆ ) Δ = BX~ − l
Φ( Xˆ ) = 0 CX~ +WX = 0
武汉大学测绘学院 黄海兰
回顾:测量平差的基本概念
四、平差结果的精度评定
(1) 根据平差后求得的改正数来估计单位权中误差,即
σˆ0 =
V T PV r
(2) 应用协因数传播律,计算观测值函数的协因数以及 观测值的平差值函数的协因数。
武汉大学测绘学院 孙
海燕 武汉大学测绘学院 黄海兰
5
2015-10-25
课程教学大纲
9 教学内容 第五章 平差模型的稳健估计 5.1 概述 5.2 稳健估计原理 5.3 基于选权迭代法的稳健估计方法 5.4 几种常用的抗差最小二乘法 5.5 相关观测的稳健估计方法
武汉大学测绘学院 孙
海燕 武汉大学测绘学院 黄海兰
相关文档
最新文档