本文最后更新于27 天前,其中的信息可能已经过时,如有错误请发送邮件到[email protected]
目录
1. GPS测量中消除误差的方法 原片 @ 00:00*
在GPS定位中,除了利用模型改正误差外,还有其他三种常用的处理误差的方法:求差法、参数法和回避法。
1.1 求差法 (Differencing Method)
Screenshot-[00:46]
- 原理:如果在地面上有两台接收机(计算机),同时对同一颗卫星进行观测(transcript误识别为“关车”)。在特定的瞬间(如0秒、10秒…),卫星钟差、轨道误差等对两台接收机的观测影响是相同的。
- 操作:将两台接收机对同一卫星的观测值相减。
- 优势:
- 不需要建立复杂的物理模型。
- 不需要精确的参数(如星历误差具体数值)。
- 能够直接消除两台接收机共有的误差(Common Errors)。
- 应用:这是GPS测量(特别是相对定位/差分定位)中经常采用的方法。
1.2 参数法 (Parameterization)
Screenshot-[01:24]
- 原理:对于某些误差,虽然有模型可以计算,但模型本身不完美(例如只能改正95%-99%),剩余的残差(Residuals)作为一个未知参数进行估计。
- 典型案例:对流层延迟 (Tropospheric Delay)。
- 对流层模型通常能消除绝大部分误差。
- 剩余的1%-5%误差作为一个位置参数(或常数),在平差计算过程中与坐标参数一起估计出来。
- 注意事项:
- 参数数量要适中:不能将所有未知数都设为参数。
- 参数相关性:如果参数过多,会导致参数之间相关性强,估计结果不稳定,精度下降。
1.3 回避法 (Avoidance Method)
Screenshot-[03:16]
- 定义:通过采用优质的硬件、软件以及选择良好的观测条件来规避误差。
- 适用场景:针对那些模型难以建立、求差法无法消除、且难以作为参数估计的误差。
- 典型案例:多路径误差 (Multipath Error)。
- 难点:多路径效应取决于信号反射点的位置、反射系数、卫星几何关系等,且随卫星运动不断变化,极难建模。
- 措施:
- 选址 (Site Selection):这是最有效的方法。避免将测站设在围墙边、高楼旁或大面积水域(如湖边)附近。
- 有意识地避开信号反射物。
2. 相对论效应概述 原片 @ 06:29*
从本节开始,将详细介绍具体的误差源及其处理方法。首先介绍的是相对论效应 (Relativistic Effects)。
Screenshot-[06:51]
- 背景:GPS精密测量必须考虑相对论效应。
- 产生原因:
- 运动速度不同:卫星高速运动,地面接收机随地球自转,两者速度不同且相对运动(狭义相对论)。
- 引力位不同:卫星处于高空(约20200km),地面处于低空,两者所受地球引力位不同(广义相对论)。
- 影响:导致卫星钟与地面钟的频率(时间流逝速度)产生差异。如果两者效应相同则无影响,但实际上两者处于不同状态,必须进行改正。
3. 狭义相对论效应 原片 @ 08:07*
Screenshot-[08:45]
- 现象:运动的时钟比静止的时钟走得慢(时间膨胀)。
- 公式推导:
设卫星钟在惯性空间中的频率为 f,卫星运行速度为 v_s。根据狭义相对论,卫星上观测到的频率 f_s 为:
f_s = f \sqrt{1 – (\frac{v_s}{c})^2} 其中 c 为光速。 - 近似计算:
由于 v_s \ll c,利用泰勒级数展开并略去高阶项:
\frac{\Delta f}{f} \approx -\frac{1}{2} (\frac{v_s}{c})^2 - 结论:由于卫星高速运动,卫星钟的频率会降低(走慢了)。
4. 广义相对论效应 原片 @ 10:53*
Screenshot-[11:51]
- 现象:处于不同引力势场中的时钟走速不同。引力势越低(越靠近地球),时钟越慢;引力势越高,时钟越快。
- 公式:
广义相对论频移公式为:
\frac{\Delta f}{f} = \frac{\Delta \Phi}{c^2} 其中 \Delta \Phi 为引力位差。 - 模型简化:
将地球视为圆球(或质点),引力位 \Phi = \frac{\mu}{r}。- 地面引力位:\Phi_e = \frac{\mu}{R} (R为地球半径)
- 卫星引力位:\Phi_s = \frac{\mu}{r} (r为卫星轨道半径)
其中 \mu = GM 为地球引力常数。
- 结论:卫星所处位置引力较弱,导致卫星钟频率升高(走快了)。
5. 总的相对论效应与改正措施 原片 @ 13:41*
5.1 综合效应公式
将狭义相对论(速度影响)和广义相对论(引力影响)结合,总的频率漂移为:
\frac{\Delta f}{f} = \frac{1}{c^2} [(\frac{\mu}{R} – \frac{\mu}{r}) – \frac{1}{2}v_s^2]
对于椭圆轨道,经过推导(利用活力公式 v^2 = \mu(\frac{2}{r} – \frac{1}{a})),总的相对论效应可表示为两部分:
- 常数项(固定偏差):假设轨道为圆形 (r=a) 时产生的部分。
- 周期项(变化偏差):由于轨道是椭圆 (e \neq 0) 导致的周期性变化。
5.2 第一部分:固定频率偏差的改正 (Factory Offset)
Screenshot-[16:01]
- 计算数值:
代入常数(\mu, c, R, a),对于GPS卫星(a \approx 26560km):
\Delta f \approx 4.45 \times 10^{-10} \times f 对于基准频率 f = 10.23 \text{ MHz}:
\Delta f \approx 0.00455 \text{ Hz} - 物理意义:如果不做处理,卫星钟上天后,由于相对论效应,频率会变为 10.23 \text{ MHz} + 0.00455 \text{ Hz}。
- 处理方法(调频):
在地面制造卫星钟时,预先将频率调低。- 目标频率:10.23 \text{ MHz} – 0.00455 \text{ Hz} = 10.22999999545 \text{ MHz}。
- 结果:卫星发射入轨后,叠加相对论效应(增加 0.00455 \text{ Hz}),其实际工作频率恰好恢复为标准的 10.23 \text{ MHz}。
- 注意:这一步由卫星制造商完成,用户无需处理。
5.3 第二部分:周期性误差的改正 (Periodic Correction)
Screenshot-[22:43]
- 原因:卫星轨道不是完美的圆,存在偏心率 e(GPS卫星 e \approx 0.001 \sim 0.02)。
- 卫星距离地心的距离 r 和速度 v 随时间周期性变化。
- 变化周期与卫星运行周期一致(约11小时58分)。
- 改正公式:
需要用户根据广播星历中的参数进行计算改正:
\Delta t_r = – \frac{2 \sqrt{\mu a}}{c^2} e \sin E- e:轨道偏心率。
- E:偏近点角 (Eccentric Anomaly)。
- \sqrt{\mu a}:常数项。
- 计算方法:
- e 和 \sqrt{a} 由导航电文直接提供。
- E 是时间的函数,用户在解算开普勒方程(计算卫星位置)时必然会求出 E。
- 量级:
如果 e=0.01,该项引起的误差可达 20 \sim 30 纳秒,相当于 6米 左右的距离误差。 - 结论:这是一个不可忽视的误差,必须在接收机解算定位时通过模型实时改正。
AI 总结
本节课程深入讲解了GPS测量中的误差处理策略及相对论效应的影响。
首先,课程总结了三种非模型化的误差处理手段:求差法通过多台设备同步观测消除共有误差;参数法将模型残余误差作为未知数估计(如对流层延迟);回避法则强调通过选址避开多路径效应等难以建模的干扰。
随后,重点分析了相对论效应。这是由于卫星高速运动(狭义相对论导致钟变慢)和处于高空低引力位(广义相对论导致钟变快)共同作用的结果。
- 固定部分:导致卫星钟频率由于综合效应实际上升约 0.00455 \text{ Hz}。这通过在地面出厂前将原子钟频率预先调低至 10.22999999545 \text{ MHz} 来抵消,用户无需干预。
- 周期性部分:由于轨道偏心率(e \neq 0)引起,最大可造成约6米的距离误差。这部分无法预先消除,必须由用户(接收机)利用广播星历中的参数(e \sin E)进行实时计算和改正。