§5.2 点位误差

官方工学老师·20 页·深入(追求细节与边界)·0 次浏览·2 天前
测量误差几何精度数学建模量化分析

点位误差

从测量现象到数学模型,搞清点位误差的来源、量级与极限

按 空格/→ 演示下一步

1 / 20 页

全部页面点击任意一页,跳回舞台从这页播放

测量误差几何精度数学建模量化分析

点位误差

从测量现象到数学模型,搞清点位误差的来源、量级与极限

1第 1 页 · 点位误差

测量误差无处不在

上一页我们定义了什么是点位误差。但误差真的客观存在吗?用定位测量举个例子:同一根标杆,你站过去目测它的位置——每次的结果都一样吗?

重复测量结果不一致
同一观测条件下,对同一目标多次定位,所得坐标不可能完全相同
四大误差来源
仪器不完美、环境在变、人会手抖、方法本身有缺陷——四类因素叠加在一起
随机但有统计规律
看似无序的散布其实服从概率分布,可用数学模型描述其大小和方向
误差会逐级累积
测角、测距的微小误差在坐标推算中不断叠加放大,影响最终点位
射击打靶对应 →定位测量

瞄得再准也会散布在靶心周围——散布客观存在,与定位测量中的误差完全对应

2第 2 页 · 测量误差无处不在

为什么需要点位中误差

上一页我们看到点位在 x、y 两个方向都有误差。那评价一个点的精度,难道要分别说「x 方向几毫米、y 方向几毫米」吗?

误差有多个方向
点位在 x、y 两个方向都存在误差分量,分别报告既繁琐又难以比较
需要综合指标
用单个数字回答「这个点的精度怎么样」,方便横向对比和等级评定
中误差的合成
点位中误差把两方向的分量合成一个数,等价于平面上几何距离的统计形式
力的合成对应 →点位中误差

两个方向的分量像分力,合起来等价于合力,整体描述点的偏离

mp=±mx2+my2m_p = \pm\sqrt{m_x^2 + m_y^2}
3第 3 页 · 为什么需要点位中误差

点位中误差的定义

上一页我们知道了为什么需要点位中误差。这一页,我们从测量学角度严格定义它——它到底在统计什么、几何上又长什么样。

统计定义
点位中误差是 x、y 两方向坐标中误差的平方和再开根号
计算公式
m_p = √(m_x² + m_y²),前提是两方向误差相互独立
几何意义
平面上,它是以观测点为圆心的「误差圆」半径
概率含义
真点位落入该圆内的概率约为 63.2%(二维正态 1σ)
飞镖反复投靶对应 →点位中误差的几何含义

落点呈散点云,云的「特征半径」就是点位中误差

mp=mx2+my2m_p = \sqrt{m_x^2 + m_y^2}
4第 4 页 · 点位中误差的定义

点位中误差的计算方法

三种方法殊途同归,从最直观的分量加和,到矩阵的迹,再到最终点位标准差。

1
分量求和
把 x、y、z 三方向的方差直接相加,得到点位方差
2
矩阵取迹
把分量方差装进协方差矩阵 Q,对角线之和即 trace(Q)
3
开方得σP
对点位方差开根号,得到点位标准差 σ_P
5第 5 页 · 点位中误差的计算方法

方向位差的概念

上一页算出的点位中误差 m_p 是个标量,回答「点整体精度如何」。但工程更常问:沿隧道掘进方向、桥跨方向误差多大?这要把整体精度「拆到」特定方向上——这就是方向位差。

方向位差的定义
位置误差在某特定方向 φ 上的投影分量,衡量该方向上的误差大小
位差计算公式
各方向分量误差的加权和,含 x、y 方向误差及协方差项
与点位中误差的区别
位差随方向 φ 变化;点位中误差是无方向标量。仅各向同性时二者才相等
极值与误差椭圆
存在最大/最小位差方向,对应误差椭圆长、短半轴;各方向位差端点构成误差椭圆
手电筒斜照椭圆石子对应 →方向位差

石子在不同光照角度下投出不同长度的影子——影子长=该方向位差;石子整体大小=点位中误差

mφ2=mx2cos2φ+my2sin2φ+2mxysinφcosφm_\varphi^2 = m_x^2\cos^2\varphi + m_y^2\sin^2\varphi + 2m_{xy}\sin\varphi\cos\varphi
6第 6 页 · 方向位差的概念

协方差矩阵的几何意义

从散点出发,看 Σ 如何一步步决定误差椭圆的形状、大小与朝向。

图解渲染中…
c1二维点位误差的完整描述:σx²、σxy、σy² 三量e1特征值,椭圆两个半轴长 = √λe2特征向量,决定椭圆旋转到哪朝g1|Σ| 是广义方差,反映点位整体离散
7第 7 页 · 协方差矩阵的几何意义

方向位差的几何直观

上一页我们看到协方差矩阵对应一个误差椭圆。实际测量中我们更关心沿某个特定方向的精度——比如沿路线、垂直基线。误差椭圆不是各向同性的圆,沿不同方向切一刀,位差不一样。

误差椭圆
点位不确定度在二维平面上的几何形状,每个方向上的'宽度'并不相等
方向截面
沿指定方向画直线交椭圆,中心到交点的距离即为该方向的位差
长短轴极值
误差最大方向沿椭圆长轴,最小沿短轴,两者互相垂直
手电筒斜照地面对应 →误差椭圆在不同方向的位差

椭圆形光斑沿不同方向量出的'半径'长短不一,对应不同方向的位差

σϕ2=σx2cos2ϕ+σy2sin2ϕ+σxysin2ϕ\sigma_\phi^2 = \sigma_x^2\cos^2\phi + \sigma_y^2\sin^2\phi + \sigma_{xy}\sin 2\phi
8第 8 页 · 方向位差的几何直观

方向位差公式推导

把二维协方差矩阵沿任意方向φ投影,逐步得到该方向上的方差与位差。

1
设方向向量
用单位向量u=(cosφ,sinφ)ᵀ表示任意方向φ
2
建立点位误差
误差向量d=(dx,dy)ᵀ在u方向的投影为δφ=uᵀd
3
写成协方差型
由E(ddᵀ)=Σ,二次型dᵀΣd就是投影方差
4
展开并代入
令c=cosφ、s=sinφ,Q=σx²c²+2σxysc+σy²s²
5
得到方向位差
开方得到σφ=√Q,即该方向一个标准差对应的位差
9第 9 页 · 方向位差公式推导

方向位差公式的物理意义

上一页我们推导出方向位差公式,现在拆开看每项对应协方差矩阵的哪个元素。尤其是 σ_xy 那项——它描述的不是两个方向各自的误差,而是两个方向误差之间是否同步。

σ²_x cos²φ
x方向方差的贡献,φ=0°时该项主导,即纯x方向位差
σ²_y sin²φ
y方向方差的贡献,φ=90°时该项主导,即纯y方向位差
2σ_xy sinφcosφ
协方差的交叉贡献,φ=45°时最大,反映x、y误差的联动
σ_xy=0 的简化
x、y独立时公式退化为σ²_φ=σ²_x cos²φ+σ²_y sin²φ
打靶时枪口在x和y方向的跳动对应 →方向位差公式三项

平方项是各自方向的固有抖动;交叉项是「x偏右时y也偏下」这种联动,对角方向影响最大

σφ2=σx2cos2φ+2σxysinφcosφ+σy2sin2φσ²_φ = σ²_x\cos^2φ + 2σ_{xy}\sinφ\cosφ + σ²_y\sin^2φ
10第 10 页 · 方向位差公式的物理意义

位差的极大值与极小值

上页推导出任意方向 φ 上的位差公式 σφ²,但 φ 还没定——哪个方向误差最大、哪个方向最小?这就是极大值 E 和极小值 F 要回答的问题。

主方向 φ_E
使位差 σφ² 取极值时的方向角,正交方向就是 F 方向
极大值 E
所有方向上位差的最大值,对应误差椭圆的长半轴
极小值 F
所有方向上位差的最小值,对应误差椭圆的短半轴
两个不变量
E+F=σx²+σy²,E·F=σx²σy²-σxy²,坐标轴任意旋转都不变
斜拉一块方形橡皮膜对应 →点位误差在所有方向的分布

被斜拉后变椭圆:长轴方向误差最大,短轴方向误差最小

tan2φE=2σxyσx2σy2E=12(σx2+σy2+(σx2σy2)2+4σxy2)F=12(σx2+σy2(σx2σy2)2+4σxy2)\tan 2\varphi_E=\frac{2\sigma_{xy}}{\sigma_x^2-\sigma_y^2}\\E=\tfrac{1}{2}(\sigma_x^2+\sigma_y^2+\sqrt{(\sigma_x^2-\sigma_y^2)^2+4\sigma_{xy}^2})\\F=\tfrac{1}{2}(\sigma_x^2+\sigma_y^2-\sqrt{(\sigma_x^2-\sigma_y^2)^2+4\sigma_{xy}^2})
11第 11 页 · 位差的极大值与极小值

极值方向的推导

从方向位差公式出发,用求导方法找出主方向角度。

1
列方向位差
设方向角φ,写出位差σφ²作为φ函数的完整三角表达式
2
对φ求导
对σφ²求一阶导数dσφ²/dφ并令其为零,得到极值条件方程
3
二倍角化简
用二倍角公式整理,得到含tan(2φ)的有理方程
4
解出主方向
解方程得两个主方向角,它们之差恰为90°
5
代入求极值
把角度代回位差公式,得到极大值λ₁与极小值λ₂
12第 12 页 · 极值方向的推导

E与F的几何意义

上一页我们解出了位差取极大、极小的方向,那两个方向上的极值 E、F 本身在误差椭圆里对应什么?这是这一页要回答的问题。

E 是长半轴
位差极大值就是误差椭圆的长半轴,对应点位最不确定的主方向
F 是短半轴
位差极小值就是误差椭圆的短半轴,对应与 E 垂直方向上的最小不确定度
两方向正交
E、F 所在的两个极值方向彼此垂直,构成误差椭圆的两条主轴
椭圆即落点区
误差椭圆勾勒出点位真误差的等概率边界,越靠近中心点位越可能落在此处
E²+F² 守恒
无论椭圆怎么旋转,E²+F² 都等于 σ_x²+σ_y²,误差总能量保持不变
晃动小船上打靶对应 →误差椭圆形状

船摇晃最猛的方向就是箭点散布的长边 E,最稳的方向就是短边 F

E2+F2=σx2+σy2E^2+F^2=\sigma_x^2+\sigma_y^2
13第 13 页 · E与F的几何意义

E与F表达任意方向位差

用极大值E和极小值F,把任意方向位差写成一个统一公式。

1
选主轴建系
把坐标轴转φE角,与E、F方向重合
2
交叉项清零
主轴下Qξη=0,公式只剩自协方差项
3
写主轴位差
主轴系下m²φ'=E²cos²φ'+F²sin²φ'
4
回代原方向
用φ'=φ-φE代回,得到任意方向位差
14第 14 页 · E与F表达任意方向位差

误差椭圆与方向位差的关系

从误差椭圆出发,看它如何用 E 和 F 编码任意方向上的位差。

图解渲染中…
a1描述点位不确定性的几何图形a2E和F是椭圆的两个特征参数a5该方向上的位差大小
15第 15 页 · 误差椭圆与方向位差的关系

方向位差的计算示例

python

用一段 NumPy 代码,把方向位差公式落到具体数字上。

代码高亮加载中…

高亮行把 E²、F²、φ₀ 与任意方向位差四步串成一条计算链,每步都对应上一页公式的一项。

16第 16 页 · 方向位差的计算示例

两种位差表达式的对比

同一个方向位差能写出两种形式——协方差矩阵与 E-F 参数。哪种更好用?关键看手头有什么数据。

协方差矩阵形式
  • 表达式:位差² = fᵀ·Q·f(Q 为 2×2 子阵)
  • 适用:已知完整协方差矩阵
  • 优点:保留 x、y 相关性,信息完整
  • 局限:换方向需重做矩阵乘法
E-F 形式
  • 表达式:位差² = E·cos²φ + F·sin²φ
  • 适用:已知 E、F 参数(控制网报告)
  • 优点:一次三角运算即得任意方向
  • 局限:标量组合,掩盖矩阵结构
手头有协方差矩阵用左式;只有 E、F 参数就用右式。两者数学上等价——E、F 就是 Q 的特征值组合。
17第 17 页 · 两种位差表达式的对比

核心概念自测

点击作答

协方差矩阵 Σ = [[E, F], [F, G]] 描述了点位 x、y 方向误差的耦合。关于元素 F = cov(x,y) 的作用,下列说法最准确的是?

18第 18 页 · 核心概念自测

知识要点回顾

  • 点位中误差是误差在二维平面上的标量综合
  • 误差椭圆即协方差矩阵的几何化呈现
  • 方向位差随方位角变化,存在极大极小值
  • E、F 两参数可表达任意方向位差
  • 中误差与位差本质等价,适用场景互补
延伸主题:误差椭圆参数的求解协方差传播律的应用点位精度评定的工程实践
19第 19 页 · 知识要点回顾

课后思考

先自己思考,再看提示——三个问题从回顾、应用到边界,逐层深入。

1点位中误差 E 用一个标量概括点位精度,为什么这种"平均"会掩盖掉方向上的差异?

参考答案E 是各方向位差的某种加权平均,只反映误差的"大小",不反映形状。方向信息藏在协方差矩阵中——F 与极值方向才刻画了几何特征。

2已知某点的协方差矩阵,如何设计附加观测才能最有效地缩小点位误差?说说思路。

参考答案看协方差矩阵的几何——F 越大、椭圆越扁,长轴方向误差越大。在长轴方向增加观测或提高精度,收敛最快;短轴方向收益较小。最优设计即按此分配。

3本章框架建立在二维正态误差假设上。若误差显著偏离正态,误差椭圆还有意义吗?为什么?

参考答案正态假设下误差椭圆是等概率密度边界。非正态时等概率区域不再规则,误差椭圆只能近似描述中心化范围。严格意义上需要用到三阶、四阶矩等更高阶统计量。

20第 20 页 · 课后思考