Cycle slip detection and repair principle.ppt

周跳的探测与修复整周模糊度的确定 第一节周跳的探测与修复 1 屏幕扫描法2 高次差法多项式拟合法MW观测值法三差法 1 整周跳变 周跳 CycleSlips 在某一特定时刻的载波相位观测值为如果在观测过程接收机保持对卫星信号的连续跟踪 则整周模糊度将保持不变 整周计数也将保持连续 但当由于某种原因使接收机无法保持对卫星信号的连续跟踪时 在卫星信号重新被锁定后 将发生变化 而也不会与前面的值保持连续 这一现象称为整周跳变 产生周跳的原因 信号被遮挡 导致卫星信号无法被跟踪仪器故障 导致差频信号无法产生卫星信号信噪比过低 导致整周计数错误接收机在高速动态的环境下进行观测 导致接收机无法正确跟踪卫星信号卫星瞬时故障 无法产生信号 周跳的特点 只影响整周计数 周跳为波长的整数倍将影响从周跳发生时刻 历元 之后的所有观测值 周跳将使周跳发生后的所有观测值包含相同的整周计数错误 解决周跳问题的方法 探测与修复设法找出周跳发生的时间和大小参数法将周跳标记出来 引入周跳参数 进行解算 周跳的探测 修复方法 屏幕扫描法方法 人工在屏幕上观察观测值曲线的变化是否连续 特点费时 只能发现大周跳 由于原始的载波观测值变化很快 通常观察的是某种观测值的组合 如 周跳的探测 修复方法 高次差法 周跳的探测 修复方法 续 高次差法的原理由于卫星和接收机间的距离在不断变化 因而载波相位测量的观测值N0 Int Fr 也随时间在不断变化 但这种变化应是有规律的 平滑的 周跳将破坏这种规律性 对于GPS卫星而言 当求至四次差时 其值已趋向于零 残留的四次差主要是由接收机的钟误差等因素引起的 周跳的探测 修复方法 续 高次差法的问题接收机钟差对此方法有效性的影响克服接收机钟差影响的方法 卫星间求差 周跳的探测 修复方法 续 高次差法的问题即使发现相位观测值中存在数周的不规则变化 也很难判断是否存在周跳 所以双差观测值被广泛采用 周跳的探测 修复方法 多项式拟合法 为了便于用计算机计算 常采用多项式拟合的方法 即根据n个相位测量观测值拟合一个n阶多项式 据此多项式来预估下一个观测值并与实测值比较 从而来发现周跳并修正整周计数 这种方法实质上和上面介绍的高次差法是相像的 但便于计算 周跳的探测 修复方法 续 多项式拟合法的应用特点由于四次差或五次差一般巳呈偶然误差特性 无法再用函数来加以拟合 所以用多项式拟合时通常也只需取至4 5阶即可 观测值可以是真正的 非差 相位观测值 也可以是经线性组合后的虚拟观测值 单差观测值和双差观测值 周跳的探测 修复方法 MW观测值法 周跳的探测 修复方法 残差法方法根据平差后的残差 进行周跳的探测与修复特点可以发现小周跳 载波相位双差观测值的残差图 第二节整周模糊度的确定 整周未知数 整周模糊度 Ambiguity 1 静态相对定位中常用的几种方法 待定参数法 经典方法1 取整法2 置信区间法XNi为模糊度的实数解mXNi s0 QNiNi 1 2为该参数的中误差置信区间为 XNi b mXNi XNi b mXNi b xt f 2 根据自由度 f n u 和置信水平 1 从t分布的数值表中查取 如 f 2500 1 99 9 b 3 28整数解在置信区间之内 3 模糊函数法 一 整数解 基本方法1 求初始解确定基线向量的实数解和整周未知数的实数解2 将整周模糊度固定为整数3 求固定解二 实数解 基线较长 误差相关性减弱 初始解的误差将随之增大 从而使模糊度参数很难固定 整数化的意义不大 2 快速定位中常用的方法 走走停停和快速静态定位法是两种具有代表性的快速定位法 确定整周未知数的方法 一 走走停停法 StopandGo 已知基线法交换天线法二 快速静态定位法快速模糊度解算法 FARA 已知基线法 将已修复周跳 剔除粗差后的双差载波相位观测值组成法方程式 然后将已知的基线向量代入法方程式并求解模糊度参数 最后再用取整法或置信区间法将求得的实数模糊度固定为整数 交换天线法 快速模糊度解算法 FARA 由瑞士的E Frei和G Beutler提出过程 快速模糊度解算法 FARA 续 1 搜索候选模糊度 根据P XNi XNAi bmXNi 1 XNi为模糊度的实数解XNAi为相应的候选整数解mXNi s0 qNiNi 1 2为该参数的中误差b xt f 2 根据自由度 f n u 和置信水平 1 从t分布的数值表中查取 这样将 XNi b mXNi XNi b mXNi 中的所有模糊度值挑选出来 构成很多候选模糊度组合 快速模糊度解算法 FARA 续 2 确定最优整数模糊度组合 快速模糊度解算法 FARA 续 3 对备选模糊度组合进行数理统计检验1 互差检验 对XNAik XNAi XNAk进行检核 P XNijk XNAik b mxNik 1 整数模糊度实数差 XNik XNi XNk i k 1 2 r i k 对应的候选整数模糊度差 XNAik XNAi XNAkmXNik s0 qNiNi 2qNiNk qNkNk 1 22 双频检验XNi XNk分别表示对同一卫星的L1 L2载波模糊度的实数解 令 XLik XNi XNk 2 1 XLAik XNAi XNAk 2 1 P XLik XLAik b mXNLik 1 4 确认最优解的三项统计检验 将搜索出来的最优整数模糊度组合 代回原法方程式平差计算 得出基线向量解和方差阵 基线向量的整数解和初始解的一致性检验 整数解和初始解的单位权中误差的一致性检验 整数解中最小单位权中误差与次最小单位权中误差间的显著性检验 快速模糊度解算法 FARA 续 3 动态定位中常用的方法 一 初始化法运动载体处于静止状态时与地面基准