vs做网站的书籍,WordPress大前端DUX5.1,wordpress右键菜单插件,和各大网站做视频的工作第12讲 GNSS授时与PPS
12.1 为什么需要高精度时间
授时的传统理解就是时间传递或者对时。比如以前手机没这么方便时#xff0c;大家还都使用石英钟手表看时间时#xff0c;大家都习惯晚上七点准时对着中央一套的报时#xff0c;来校准你家的机械钟或者挂钟#xff0c;这就…第12讲 GNSS授时与PPS
12.1 为什么需要高精度时间
授时的传统理解就是时间传递或者对时。比如以前手机没这么方便时大家还都使用石英钟手表看时间时大家都习惯晚上七点准时对着中央一套的报时来校准你家的机械钟或者挂钟这就是最简单的时间传递这种时间传递的精度秒的量级对于日常生活已经足够。
但是全世界这么多表如何将大家的时间都统一起来呢。以前的笨方法以及现在还在使用的方法就是大家都把表放到一块一块对比校准一下不就可以做到统一了嘛。但是全世界这么多设备如何大家统一呢。像电视就是时间传递的一种其他的还有电话、电波授时等。
但是此类的授时精度都相对较低但有些场景我们又要比较高的时间精度。比如两个摄像头拍摄的画面如果每个摄像头同1s内拍摄了10000张图片如果我们要将两个设备拍摄的图像拼接起来那么两个设备的时间相差要在1/10000秒以下这样我们才能对齐指定的画面。在我们传统的电力、移动通信等行业有些时间的精度要在纳秒量级也就是一亿分之一秒。
在分布式系统中由于物理上的分散性系统无法为彼此间相互独立的模块提供一个统一的全局时钟而由各个进程或模块各自维护它们的本地时钟。由于这些本地时钟的计时速率、运行环境存在不一致性因此即使所有本地时钟在某一时刻都被校准一段时间后这些本地时钟也会出现不一致。为了这些本地时钟再次达到相同的时间值必须进行时间同步操作。授时技术又是实现时间同步的重要技术手段。
精密时间是科学研究、科学实验和工程技术诸方面的基本物理参量。它为一切动力学系统和时序过程的测量和定量研究提供了必不可少的时基坐标。
精密授时在以通信、电力、控制等工业领域和国防领域有着广泛和重要的应用。现代武器实试验、战争需要它保障智能化交通运输系统的建立和数字化地球的实现需要它支持。现代通信网和电力网建设也越来越增强了对精度时间和频率的依赖。从建立一个现代化国家的大系统工程总体考虑导航定位和授时系统应该说是基础的基础。它对整体社会的支撑几乎是全方位的。
卫星发射
火箭发射升空后不论是发射场、测控站以及测量船都会对火箭进行连续测控。而所有的测控都会根据火箭的飞行轨迹和速度推算出火箭将出现的方位和时间。这时候如果各个测控站的时间不同步就有极大的可能追踪火箭失败。
航天活动
平常在电影中或者新闻里都会看到飞船和目标飞行器交会对接的画面。就像是我国“天舟一号”以及“天宫二号”的交会对接都是需要精密的时间同步对两个飞行器进行姿态的同步观测实时发出测控指令才能保证对接成功否则“太空之吻”的美谈就很有可能变为“太空车祸”。
金融领域
现代金融行业无论银行或证券交易所,都离不开计算机和计算机联网单个计算机业务和计算机联网业务如银行业务往来的发生时刻、金融交易的准确时刻、E-mail 信息和访问时间、数据库处理时间、银联卡/账户的密码识别等等都涉及银行联网计算机之间的时间同步和频率同步都需要精准授时。
电力系统
对于电力系统而言每一条由电厂汇入总干线的电流在同一时刻具有不同的相位。为让各电厂在每一时刻输送的电流按照预定的相位而不至于发生在本该送入波峰而送入反向的波谷导致短路等情况电网整体的工作协调需要一个全局统一的时钟完成。而电网的时钟同步正是由本文的标题”GPS秒脉冲达成时间同步“来完成。
12.2 授时原理
12.2.1 卫星单向法
使用GNSS系统进行高精度授时的原理其实已经呼之欲出我们前面提到的解算结果接收机钟差其实就是本地时钟与GNSS时间框架的差异。
单向法的误差来源是卫星位置与星历的偏差、电离层层扰动以及接收设备误差等。能接收到多颗卫星信号并进行长时间平均的十分理想的情况下该方法可以达到接近1ns左右的时标同步精度。
一般来说对已知精密坐标的固定用户观测1颗卫星就可以实现精密的时间测量或者同步。若观测2颗卫星或者更多卫星则提供了更多的观测量提高了定时的稳健性。
12.2.2 卫星共视法
两个地面站A和B同时观测一枚GNSS卫星通过一定的同时观测时间安排并通过网络传输进行对比也可以实现授时。
12.3 高精度时间传递-PPS
PPS秒脉冲Pulse Per Second
形式化定义GPS接收模块(GPS Receiver)可以产生秒脉冲信号其作用是指示整秒的时间间隔并通过PPS的上升沿标示GPS给出的UTC(Coordinated Universal Time)时间对应的具体时刻。其指示整秒时间精度可达数十纳秒级无累积误差。每秒产生的该PPS与其给出的整秒时刻可以用于晶振精确调频⸺时钟同步。
GPS秒脉冲信号 pps 一秒钟一个它的作用是用来指示整秒的时刻而该时刻通常是用PPS秒脉冲的上升沿来标示。因为GPS可以给出UTC时间但用户收到时是会有延时的那么为了精确授时引入PPS信号上升沿来标示UTC的整秒时刻精度可以到纳秒级而且没有累积误差。
第13讲 仿距差分定位原理-站间单差
从单点定位的定位流程中我们了解到卫星的轨道、钟差以及大气都存在误差这是导致单点定位精度不佳的原因。所以我们再次列出伪距的观测值表达式。 P r , i s ρ r , i s − c δ t s c δ t r , i I r , i s T r , i s d r , i s − d i s ϵ P P_{r,i}^{s} \rho_{r,i}^{s} - c\delta t^{s} c\delta t_{r,i} I_{r,i}^{s} T_{r,i}^{s} d_{r,i}^{s} - d_{i}^{s} \epsilon_{P} Pr,isρr,is−cδtscδtr,iIr,isTr,isdr,is−disϵP
其中 P r , i s P_{r,i}^{s} Pr,is 为卫星 s s s 到接收机 r r r 的第 i i i 频点的伪距观测值 ρ r , i s \rho_{r,i}^{s} ρr,is 为卫星 s s s天线相位中心 到接收机 r r r 天线相位中心的几何距离真实距离 c δ t s c\delta t^{s} cδts 为卫星 s s s 钟误差 c δ t r , i c\delta t_{r,i} cδtr,i 为接收机 r r r 钟误差 I r , i s I_{r,i}^{s} Ir,is 为卫星 s s s 到接收机 r r r 的第 i i i 频点的电离层延迟与频率有关 T r , i s T_{r,i}^{s} Tr,is 为卫星 s s s 到接收机 r r r 的对流层延迟与频率无关 d r , i s d_{r,i}^{s} dr,is 为接收机端第 i i i 频点的伪距硬件延迟与频率有关 d i s d_{i}^{s} dis 为卫星端第 i i i 频点的伪距硬件延迟与频率有关 ϵ P \epsilon_{P} ϵP 为伪距观测值噪声
单位为米
如果我们有两台接收机分别为 r 1 r_1 r1 和 r 2 r_2 r2那么两台接收机的同时观测到卫星 s s s 同一频点伪距观测值如下 P r 1 s ρ r 1 s − c δ t s c δ t r 1 I r 1 s T r 1 s d r 1 − d s ϵ P P_{r_1}^{s} \rho_{r_1}^{s} - c\delta t^{s} c\delta t_{r_1} I_{r_1}^{s} T_{r_1}^{s} d_{r_1} - d^{s} \epsilon_{P} Pr1sρr1s−cδtscδtr1Ir1sTr1sdr1−dsϵP P r 2 s ρ r 2 s − c δ t s c δ t r 2 I r 2 s T r 2 s d r 2 − d s ϵ P P_{r_2}^{s} \rho_{r_2}^{s} - c\delta t^{s} c\delta t_{r_2} I_{r_2}^{s} T_{r_2}^{s} d_{r_2} - d^{s} \epsilon_{P} Pr2sρr2s−cδtscδtr2Ir2sTr2sdr2−dsϵP
其中下标 r 1 r_1 r1 代表1号接收机 r 2 r_2 r2 代表2号接收机。
有关系的项 c δ t s c\delta t^{s} cδts 是卫星钟差对接收机 r 1 r_1 r1 和 r 2 r_2 r2 都有相同的影响。 ϵ P \epsilon_{P} ϵP 伪距测量噪声同样会影响到接收机 r 1 r_1 r1 和 r 2 r_2 r2 的测量结果。 d s d^{s} ds 是卫星硬件延迟对接收机 r 1 r_1 r1 和 r 2 r_2 r2 都有相同的影响。
相互没关系的项 ρ r 1 s \rho_{r_1}^{s} ρr1s 和 ρ r 2 s \rho_{r_2}^{s} ρr2s 表示各自接收机到卫星的真实距离。 c δ t r 1 c\delta t_{r_1} cδtr1 和 c δ t r 2 c\delta t_{r_2} cδtr2 表示各自接收机的钟差。 I r 1 s I_{r_1}^{s} Ir1s 和 I r 2 s I_{r_2}^{s} Ir2s 表示各自接收机的电离层延迟。 T r 1 s T_{r_1}^{s} Tr1s 和 T r 2 s T_{r_2}^{s} Tr2s 表示各自接收机的对流层延迟。 d r 1 d_{r_1} dr1 和 d r 2 d_{r_2} dr2 表示各自接收机的硬件延迟。
13.1 大气误差以及卫星硬件延迟
如果对同一个卫星同一频率的两个接收机观测值进行作差即站间单差 P r 1 − r 2 s ρ r 1 − r 2 s c δ t r 1 − r 2 I r 1 − r 2 s T r 1 − r 2 s d r 1 − r 2 s ϵ P r 1 − r 2 P_{r_1 - r_2}^{s} \rho_{r_1 - r_2}^{s} c\delta t_{r_1 - r_2} I_{r_1 - r_2}^{s} T_{r_1 - r_2}^{s} d_{r_1 - r_2}^{s} \epsilon_{P_{r_1 - r_2}} Pr1−r2sρr1−r2scδtr1−r2Ir1−r2sTr1−r2sdr1−r2sϵPr1−r2
解释这里每一项是什么的差 P r 1 − r 2 s P_{r_1 - r_2}^{s} Pr1−r2s: 第 s s s 频点的站间单差伪距离测量值。表示在两个接收机站点 r 1 r_1 r1 和 r 2 r_2 r2 之间的差分伪距测量值。 ρ r 1 − r 2 s \rho_{r_1 - r_2}^{s} ρr1−r2s: 第 s s s 频点的站间单差真实距离。表示在两个接收机站点 r 1 r_1 r1 和 r 2 r_2 r2 之间的真实距离差。 c δ t r 1 − r 2 c\delta t_{r_1 - r_2} cδtr1−r2: 站间单差钟差项。表示两个接收机站点之间的钟差以光速 c c c 乘以钟差 δ t r 1 − r 2 \delta t_{r_1 - r_2} δtr1−r2 表示。 I r 1 − r 2 s I_{r_1 - r_2}^{s} Ir1−r2s: 第 s s s 频点的站间单差电离层延迟。表示两个接收机站点之间的电离层延迟差。 T r 1 − r 2 s T_{r_1 - r_2}^{s} Tr1−r2s: 第 s s s 频点的站间单差对流层延迟。表示两个接收机站点之间的对流层延迟差。 d r 1 − r 2 s d_{r_1 - r_2}^{s} dr1−r2s: 第 s s s 频点的站间单差硬件延迟。表示两个接收机站点之间的硬件延迟差。 ϵ P r 1 − r 2 \epsilon_{P_{r_1 - r_2}} ϵPr1−r2: 伪距测量噪声。表示伪距测量过程中引入的噪声。站间单差噪声被放大。
公式相减消除的项 c δ t s c\delta t^{s} cδts: 卫星钟差项因为它对两个接收机是相同的在作差时会被消除。 d s d^{s} ds: 卫星硬件延迟项同样在作差时会被消除。
因为我们卫星高度在两万公里左右如果两个接收机距离较近小于20km卫星信号的传播路径十分相似电离层和对流层延迟也基本相等。所以有 I r 1 − r 2 s ≈ 0 I_{r_1 - r_2}^{s} \approx 0 Ir1−r2s≈0 T r 1 − r 2 s ≈ 0 T_{r_1 - r_2}^{s} \approx 0 Tr1−r2s≈0
即有 P r 1 − r 2 s ρ r 1 − r 2 s c δ t r 1 − r 2 d r 1 − r 2 s ϵ P r 1 − r 2 P_{r_1 - r_2}^{s} \rho_{r_1 - r_2}^{s} c\delta t_{r_1 - r_2} d_{r_1 - r_2}^{s} \epsilon_{P_{r_1 - r_2}} Pr1−r2sρr1−r2scδtr1−r2dr1−r2sϵPr1−r2 ρ r 1 − r 2 s \rho_{r_1 - r_2}^{s} ρr1−r2s: 第 s s s 频点的站间单差真实距离。表示两个接收机站点之间的真实距离差。 这是几何距离差分项。 c δ t r 1 − r 2 c\delta t_{r_1 - r_2} cδtr1−r2: 站间单差钟差项。表示两个接收机站点之间的钟差。 d r 1 − r 2 s d_{r_1 - r_2}^{s} dr1−r2s: 第 s s s 频点的站间单差硬件延迟。表示两个接收机站点之间的硬件延迟差。 ϵ P r 1 − r 2 \epsilon_{P_{r_1 - r_2}} ϵPr1−r2: 伪距测量噪声。表示伪距测量过程中引入的噪声。
13.2 轨道误差
13.2.1 单站几何距离公式
几何距离可以表示为卫星轨道和接收机位置的形式。 ρ r s ( r ⃗ s − r ⃗ r ) ⋅ u ⃗ r s \rho_{r}^{s} (\vec{r}^{s} - \vec{r}_{r}) \cdot \vec{u}_{r}^{s} ρrs(r s−r r)⋅u rs
[ 向量*单位向量 模 ]
其中 r ⃗ s ( X s , Y s , Z s ) \vec{r}^{s} (X^{s}, Y^{s}, Z^{s}) r s(Xs,Ys,Zs) 为卫星位置向量 r ⃗ r ( X r , Y r , Z r ) \vec{r}_{r} (X_{r}, Y_{r}, Z_{r}) r r(Xr,Yr,Zr) 为接收机位置向量 u ⃗ r s \vec{u}_{r}^{s} u rs 为视线方向单位向量
从单点定位理论我们如果只知大概路径位置那可以通过将几何距离表示概略路径位置和位置偏差的形式。如下 R ⃗ 0 s ( r ⃗ s − r ⃗ 0 ) ⋅ u ⃗ r s \vec{R}_{0}^{s} (\vec{r}^{s} - \vec{r}_{0}) \cdot \vec{u}_{r}^{s} R 0s(r s−r 0)⋅u rs ρ r s ( r ⃗ s − ( r ⃗ 0 r d r ⃗ ) ) ⋅ u ⃗ r s R 0 s − d r ⃗ ⋅ u ⃗ r s \rho_r^s (\vec{r}^s - (\vec{r}_0^r \vec{dr})) \cdot \vec{u}_r^s R_0^s - \vec{dr} \cdot \vec{u}_r^s ρrs(r s−(r 0rdr ))⋅u rsR0s−dr ⋅u rs
[ 单位向量概略位置和真实位置偏差不大10m以内单位向量忽略]
[ 最终表达是为概略位置位置偏差与单位向量的内积 ]
其中 r ⃗ 0 ( X 0 , Y 0 , Z 0 ) \vec{r}_{0} (X_{0}, Y_{0}, Z_{0}) r 0(X0,Y0,Z0) 为接收机概略位置向量 d r ⃗ ( d x , d y , d z ) \vec{dr} (dx, dy, dz) dr (dx,dy,dz) 为接收机概略位置的偏差向量 u ⃗ r s ( l s , m s , n s ) \vec{u}_{r}^{s} (l^{s}, m^{s}, n^{s}) u rs(ls,ms,ns) 为方向余弦向量其中 l s X s − X 0 R 0 s l^{s} \frac{X^{s} - X_{0}}{R_{0}^{s}} lsR0sXs−X0 m s Y s − Y 0 R 0 s m^{s} \frac{Y^{s} - Y_{0}}{R_{0}^{s}} msR0sYs−Y0 n s Z s − Z 0 R 0 s n^{s} \frac{Z^{s} - Z_{0}}{R_{0}^{s}} nsR0sZs−Z0 R ⃗ 0 s \vec{R}_{0}^{s} R 0s 为概略位置到卫星的几何距离
附加A u ⃗ r s ( l s , m s , n s ) \vec{u}_{r}^{s} (l^{s}, m^{s}, n^{s}) u rs(ls,ms,ns) 为方向余弦向量表示从接收机站点到卫星的单位向量方向。
方向余弦向量的分量公式
方向余弦向量的三个分量 l s , m s , n s l^{s}, m^{s}, n^{s} ls,ms,ns 分别表示方向向量在三个坐标轴上的投影长度与其向量长度之比 l s X s − X 0 R 0 s l^{s} \frac{X^{s} - X_{0}}{R_{0}^{s}} lsR0sXs−X0 m s Y s − Y 0 R 0 s m^{s} \frac{Y^{s} - Y_{0}}{R_{0}^{s}} msR0sYs−Y0 n s Z s − Z 0 R 0 s n^{s} \frac{Z^{s} - Z_{0}}{R_{0}^{s}} nsR0sZs−Z0
其中 X s , Y s , Z s X^{s}, Y^{s}, Z^{s} Xs,Ys,Zs 是卫星在三维空间中的坐标。 X 0 , Y 0 , Z 0 X_{0}, Y_{0}, Z_{0} X0,Y0,Z0 是接收机站点在三维空间中的坐标。 R 0 s R_{0}^{s} R0s 是接收机站点到卫星的距离表示为 R 0 s ( X s − X 0 ) 2 ( Y s − Y 0 ) 2 ( Z s − Z 0 ) 2 R_{0}^{s} \sqrt{(X^{s} - X_{0})^2 (Y^{s} - Y_{0})^2 (Z^{s} - Z_{0})^2} R0s(Xs−X0)2(Ys−Y0)2(Zs−Z0)2 。
通过这些公式可以得到从接收机到卫星的方向余弦向量 u ⃗ r s \vec{u}_{r}^{s} u rs它表示的是单位向量方向。
附加B公式推导 几何距离表示 卫星 s s s 与接收机 r r r 之间的几何距离可以表示为它们位置向量的差 ρ r s ∣ r ⃗ s − r ⃗ r ∣ \rho_r^s \left| \vec{r}^s - \vec{r}_r \right| ρrs∣r s−r r∣ 这里 r ⃗ s \vec{r}^s r s 是卫星的位置向量 r ⃗ r \vec{r}_r r r 是接收机的位置向量。 单点定位近似 如果我们知道接收机的概略位置 r ⃗ 0 r \vec{r}_0^r r 0r 和接收机位置的误差 d r ⃗ \vec{dr} dr 接收机的真实位置 r ⃗ r \vec{r}_r r r 可以表示为 r ⃗ r r ⃗ 0 r d r ⃗ \vec{r}_r \vec{r}_0^r \vec{dr} r rr 0rdr 代入接收机的位置 将接收机的真实位置代入几何距离公式 ρ r s ∣ r ⃗ s − ( r ⃗ 0 r d r ⃗ ) ∣ \rho_r^s \left| \vec{r}^s - (\vec{r}_0^r \vec{dr}) \right| ρrs r s−(r 0rdr ) 点积表示 为了简化我们通常用一个指向卫星的单位方向向量 u ⃗ r s \vec{u}_r^s u rs 表示距离的点积 ρ r s ( r ⃗ s − ( r ⃗ 0 r d r ⃗ ) ) ⋅ u ⃗ r s \rho_r^s (\vec{r}^s - (\vec{r}_0^r \vec{dr})) \cdot \vec{u}_r^s ρrs(r s−(r 0rdr ))⋅u rs 这里 u ⃗ r s \vec{u}_r^s u rs 是从接收机指向卫星的单位方向向量。 近似距离和误差项 使用接收机的概略位置 r ⃗ 0 r \vec{r}_0^r r 0r 计算的近似距离可以定义为 R 0 s ( r ⃗ s − r ⃗ 0 r ) ⋅ u ⃗ r s R_0^s (\vec{r}^s - \vec{r}_0^r) \cdot \vec{u}_r^s R0s(r s−r 0r)⋅u rs 简化表达式 将近似距离代入原始公式并简化 ρ r s ( r ⃗ s − ( r ⃗ 0 r d r ⃗ ) ) ⋅ u ⃗ r s \rho_r^s (\vec{r}^s - (\vec{r}_0^r \vec{dr})) \cdot \vec{u}_r^s ρrs(r s−(r 0rdr ))⋅u rs ρ r s ( r ⃗ s − r ⃗ 0 r − d r ⃗ ) ⋅ u ⃗ r s \rho_r^s (\vec{r}^s - \vec{r}_0^r - \vec{dr}) \cdot \vec{u}_r^s ρrs(r s−r 0r−dr )⋅u rs ρ r s ( r ⃗ s − r ⃗ 0 r ) ⋅ u ⃗ r s − d r ⃗ ⋅ u ⃗ r s \rho_r^s (\vec{r}^s - \vec{r}_0^r) \cdot \vec{u}_r^s - \vec{dr} \cdot \vec{u}_r^s ρrs(r s−r 0r)⋅u rs−dr ⋅u rs ρ r s R 0 s − d r ⃗ ⋅ u ⃗ r s \rho_r^s R_0^s - \vec{dr} \cdot \vec{u}_r^s ρrsR0s−dr ⋅u rs
最终公式
最终的公式是 近似距离 R 0 s ( r ⃗ s − r ⃗ 0 r ) ⋅ u ⃗ r s R_0^s (\vec{r}^s - \vec{r}_0^r) \cdot \vec{u}_r^s R0s(r s−r 0r)⋅u rs 考虑误差修正的真实几何距离 ρ r s ( r ⃗ s − ( r ⃗ 0 r d r ⃗ ) ) ⋅ u ⃗ r s R 0 s − d r ⃗ ⋅ u ⃗ r s \rho_r^s (\vec{r}^s - (\vec{r}_0^r \vec{dr})) \cdot \vec{u}_r^s R_0^s - \vec{dr} \cdot \vec{u}_r^s ρrs(r s−(r 0rdr ))⋅u rsR0s−dr ⋅u rs
13.2.2 站间单差几何距离公式
13.1节中我们已经通过站间单差基本消除了电离层误差和对流层误差以及完全消除了卫星伪距硬件延迟和卫星钟差的影响。
如果我们假设1号接收机为流动站即我们仅能通过单点定位计算其粗略位置而2号接收机为基准站其精确的坐标已知。 ρ r 1 − 2 s ρ r 1 s − ρ r 2 s \rho_{r1-2}^s \rho_{r1}^s - \rho_{r2}^s ρr1−2sρr1s−ρr2s ρ r 1 − 2 s \rho_{r1-2}^s ρr1−2s 表示卫星 s s s 到测站1和测站2之间的几何距离差 ρ r 1 s \rho_{r1}^s ρr1s 表示卫星 s s s 到测站1的几何距离 ρ r 2 s \rho_{r2}^s ρr2s 表示卫星 s s s 到测站2的几何距离
但是我们知道广播星历实际上有误差的如果我们假设其误差是 d r ⃗ s \vec{dr}^s dr s。 d r ⃗ s \vec{dr}^s dr s 表示卫星 s s s 的位置误差向量
所以我们重写两个测站的几何距离公式
流动站
[ 概略位置 流动站位置偏差单位向量 卫星位置偏差单位向量] ρ r 1 s ( r ⃗ s d r ⃗ s − ( r r 1 ⃗ d r r 1 ⃗ ) ) ⋅ u ⃗ r 1 s R 01 s − d r r 1 ⃗ ⋅ u ⃗ r 1 s d r ⃗ s ⋅ u ⃗ r 1 s \rho_{r1}^s (\vec{r}^s \vec{dr}^s - (\vec{r_{r1}} \vec{dr_{r1}})) \cdot \vec{u}_{r1}^s R_{01}^s - \vec{dr_{r1}} \cdot \vec{u}_{r1}^s \vec{dr}^s \cdot \vec{u}_{r1}^s ρr1s(r sdr s−(rr1 drr1 ))⋅u r1sR01s−drr1 ⋅u r1sdr s⋅u r1s r ⃗ s \vec{r}^s r s 表示卫星 s s s 的实际位置向量 r r 1 ⃗ \vec{r_{r1}} rr1 表示测站1的实际位置向量 d r r 1 ⃗ \vec{dr_{r1}} drr1 表示测站1的位置误差向量 u ⃗ r 1 s \vec{u}_{r1}^s u r1s 表示从测站1到卫星 s s s 的单位方向向量 R 01 s R_{01}^s R01s 表示由广播星历计算得到的卫星 s s s 到测站1的几何距离
基准站
[ 概略位置 卫星位置偏差*单位向量] ρ r 2 s ( r ⃗ s d r ⃗ s − ( r r 2 ⃗ ) ) ⋅ u ⃗ r 2 s R 02 s d r ⃗ s ⋅ u ⃗ r 2 s \rho_{r2}^s (\vec{r}^s \vec{dr}^s - (\vec{r_{r2}})) \cdot \vec{u}_{r2}^s R_{02}^s \vec{dr}^s \cdot \vec{u}_{r2}^s ρr2s(r sdr s−(rr2 ))⋅u r2sR02sdr s⋅u r2s r r 2 ⃗ \vec{r_{r2}} rr2 表示测站2的实际位置向量 u ⃗ r 2 s \vec{u}_{r2}^s u r2s 表示从测站2到卫星 s s s 的单位方向向量 R 02 s R_{02}^s R02s 表示由广播星历计算得到的卫星 s s s 到测站2的几何距离
上述中 R 01 s R_{01}^s R01s 代表卫星 s s s 广播星历位置到接收机1的概略位置的几何距离 R 02 s R_{02}^s R02s 代表卫星 s s s 广播星历位置到接收机2的真实几何距离其他符号相似。
那么有 ρ r 1 − 2 s ρ r 1 s − ρ r 2 s ( R 01 s − d r r 1 ⃗ ⋅ u ⃗ r 1 s d r ⃗ s ⋅ u ⃗ r 1 s ) − ( R 02 s d r ⃗ s ⋅ u ⃗ r 2 s ) \rho_{r1-2}^s \rho_{r1}^s - \rho_{r2}^s (R_{01}^s - \vec{dr_{r1}} \cdot \vec{u}_{r1}^s \vec{dr}^s \cdot \vec{u}_{r1}^s) - (R_{02}^s \vec{dr}^s \cdot \vec{u}_{r2}^s) ρr1−2sρr1s−ρr2s(R01s−drr1 ⋅u r1sdr s⋅u r1s)−(R02sdr s⋅u r2s)
那么有 ρ r 1 − 2 s ( R 01 s − R 02 s ) − d r r 1 ⃗ ⋅ u ⃗ r 1 s δ ρ s \rho_{r1-2}^s (R_{01}^s - R_{02}^s) - \vec{dr_{r1}} \cdot \vec{u}_{r1}^s \delta \rho^s ρr1−2s(R01s−R02s)−drr1 ⋅u r1sδρs
[流动站概略位置到广播星历的距离 - 基准站位置到广播星历的距离- 流动站接收机位置误差*单位向量轨道误差在视线方向上投影的差值] δ ρ s \delta \rho^s δρs 为轨道误差在视线方向投影的差值
令 δ ρ s \delta \rho^s δρs 为轨道误差在视线方向投影的差值 δ ρ s d r ⃗ s ⋅ u ⃗ r 1 s − d r ⃗ s ⋅ u ⃗ r 2 s d r ⃗ s ⋅ ( u ⃗ r 1 s − u ⃗ r 2 s ) \delta \rho^s \vec{dr}^s \cdot \vec{u}_{r1}^s - \vec{dr}^s \cdot \vec{u}_{r2}^s \vec{dr}^s \cdot (\vec{u}_{r1}^s - \vec{u}_{r2}^s) δρsdr s⋅u r1s−dr s⋅u r2sdr s⋅(u r1s−u r2s)
[ 轨道误差*两个视线单位向量的差 ] 都是极小量轨道误差可以忽略。下面给出实际计算证明
所以上述中我们主要考虑是 s s s 的轨道误差在两个测站观测方向上的投影差。
如果两个测站的方向余弦向量相差足够小那么轨道误差通过星间单差后就可以忽略。
我们来简单假设一下如果轨道误差在视线投影误差为20米当测站相距20km时那站间单差后的残余误差为 20 20000 × 20 1 20000 \frac{20}{20000} \times 20 \frac{1}{20000} 2000020×20200001米。
所以通过站间单差后轨道的误差我们也可以忽略。 13.3 差分定位流程
文章转载自： http://www.morning.dcccl.cn.gov.cn.dcccl.cn http://www.morning.rhkgz.cn.gov.cn.rhkgz.cn http://www.morning.mwbqk.cn.gov.cn.mwbqk.cn http://www.morning.kxnjg.cn.gov.cn.kxnjg.cn http://www.morning.lgsqy.cn.gov.cn.lgsqy.cn http://www.morning.smggx.cn.gov.cn.smggx.cn http://www.morning.ckwrn.cn.gov.cn.ckwrn.cn http://www.morning.zdkzj.cn.gov.cn.zdkzj.cn http://www.morning.mldrd.cn.gov.cn.mldrd.cn http://www.morning.htbgz.cn.gov.cn.htbgz.cn http://www.morning.muzishu.com.gov.cn.muzishu.com http://www.morning.ypdmr.cn.gov.cn.ypdmr.cn http://www.morning.rkfxc.cn.gov.cn.rkfxc.cn http://www.morning.zlnmm.cn.gov.cn.zlnmm.cn http://www.morning.jpnw.cn.gov.cn.jpnw.cn http://www.morning.rjqtq.cn.gov.cn.rjqtq.cn http://www.morning.hkysq.cn.gov.cn.hkysq.cn http://www.morning.mqfw.cn.gov.cn.mqfw.cn http://www.morning.lwrcg.cn.gov.cn.lwrcg.cn http://www.morning.cpljq.cn.gov.cn.cpljq.cn http://www.morning.bqmdl.cn.gov.cn.bqmdl.cn http://www.morning.ysybx.cn.gov.cn.ysybx.cn http://www.morning.lznqb.cn.gov.cn.lznqb.cn http://www.morning.rpjr.cn.gov.cn.rpjr.cn http://www.morning.dkgtr.cn.gov.cn.dkgtr.cn http://www.morning.kqbzy.cn.gov.cn.kqbzy.cn http://www.morning.xhqr.cn.gov.cn.xhqr.cn http://www.morning.xdmsq.cn.gov.cn.xdmsq.cn http://www.morning.tntqr.cn.gov.cn.tntqr.cn http://www.morning.sfrw.cn.gov.cn.sfrw.cn http://www.morning.rxxdk.cn.gov.cn.rxxdk.cn http://www.morning.yfnjk.cn.gov.cn.yfnjk.cn http://www.morning.rcrnw.cn.gov.cn.rcrnw.cn http://www.morning.hxxyp.cn.gov.cn.hxxyp.cn http://www.morning.zdsdn.cn.gov.cn.zdsdn.cn http://www.morning.pflpb.cn.gov.cn.pflpb.cn http://www.morning.wdhlc.cn.gov.cn.wdhlc.cn http://www.morning.bbtn.cn.gov.cn.bbtn.cn http://www.morning.trnl.cn.gov.cn.trnl.cn http://www.morning.qkkmd.cn.gov.cn.qkkmd.cn http://www.morning.gywfp.cn.gov.cn.gywfp.cn http://www.morning.fzlk.cn.gov.cn.fzlk.cn http://www.morning.fjscr.cn.gov.cn.fjscr.cn http://www.morning.dxpqd.cn.gov.cn.dxpqd.cn http://www.morning.tklqs.cn.gov.cn.tklqs.cn http://www.morning.wwthz.cn.gov.cn.wwthz.cn http://www.morning.lqypx.cn.gov.cn.lqypx.cn http://www.morning.kmbgl.cn.gov.cn.kmbgl.cn http://www.morning.rkxdp.cn.gov.cn.rkxdp.cn http://www.morning.brzlp.cn.gov.cn.brzlp.cn http://www.morning.zybdj.cn.gov.cn.zybdj.cn http://www.morning.rnytd.cn.gov.cn.rnytd.cn http://www.morning.rgxll.cn.gov.cn.rgxll.cn http://www.morning.xcbnc.cn.gov.cn.xcbnc.cn http://www.morning.qmzhy.cn.gov.cn.qmzhy.cn http://www.morning.wchsx.cn.gov.cn.wchsx.cn http://www.morning.rqxtb.cn.gov.cn.rqxtb.cn http://www.morning.zrbpx.cn.gov.cn.zrbpx.cn http://www.morning.nlrp.cn.gov.cn.nlrp.cn http://www.morning.ybyln.cn.gov.cn.ybyln.cn http://www.morning.mtmph.cn.gov.cn.mtmph.cn http://www.morning.ubpsa.cn.gov.cn.ubpsa.cn http://www.morning.rkdw.cn.gov.cn.rkdw.cn http://www.morning.xhklb.cn.gov.cn.xhklb.cn http://www.morning.tqrjj.cn.gov.cn.tqrjj.cn http://www.morning.frcxx.cn.gov.cn.frcxx.cn http://www.morning.pxtgf.cn.gov.cn.pxtgf.cn http://www.morning.qqpg.cn.gov.cn.qqpg.cn http://www.morning.jxhlx.cn.gov.cn.jxhlx.cn http://www.morning.qnypp.cn.gov.cn.qnypp.cn http://www.morning.smfbw.cn.gov.cn.smfbw.cn http://www.morning.c7501.cn.gov.cn.c7501.cn http://www.morning.gslz.com.cn.gov.cn.gslz.com.cn http://www.morning.lgxzj.cn.gov.cn.lgxzj.cn http://www.morning.jtwck.cn.gov.cn.jtwck.cn http://www.morning.jcffp.cn.gov.cn.jcffp.cn http://www.morning.zlbjx.cn.gov.cn.zlbjx.cn http://www.morning.dzyxr.cn.gov.cn.dzyxr.cn http://www.morning.kcwkt.cn.gov.cn.kcwkt.cn http://www.morning.tbkqs.cn.gov.cn.tbkqs.cn
