代码先锋网 代码片段及技术文章聚合

【学习笔记】GPS测量与数据处理(观测值的线性组合(单差、双差、三差))

技术标签: GPS  单差  双差  三差  GPS求差法

第一次使用播客记录自己的学习,latex公式的确很好使

1. 接收机间单差

信号发射时刻t1t1计算方法

t1=t1+Vti(t1)ρpi(t1))ct1′=t1+Vti(t1)−ρip(t1))c

式中,接收机钟的改正数Vti(t1)Vti(t1)可用t1t1时的伪距观测采用单点定位的方式求得,ρpi(t1)ρip(t1)则需要反复迭代求得
疑问:
1. 接收机中的改正数通过伪距定点定位求得
2. 其距离通过反复迭代求得

书本89页 由于标称时刻t1t1往往是一整数,而发射时刻的t1t1′有一串小数,故计算卫星位置时仍然使用标称时间t1t1来计算卫星位置,加入距离改正数ΔρΔρ:
Δρpi(t1)=ρpi(t1)[Vti(t1)ρpi(t1)c]Δρip(t1)=ρip(t1)⋅[Vti(t1)−ρip(t1)c]

式中 ρpi(t1)ρip(t1)t1t1时刻ii至卫星ρρ的距离变化率????该距离变化率可用多普勒测量的方法求定,一般的测量型GPS接收机均可给出此值。
刚才上文管ρpi(t1)ρip(t1)为卫地距,哎
接收机i,j,t1i,j,t1时刻观测P卫星的载波相位单差观测方程:
Δφpij~=fcΔρptij+fΔVptijΔNpijfcΔ(Vtrop)pijfcΔ(Vion)pijΔφijp~=fcΔρtijp+fΔVtijp−ΔNijp−fcΔ(Vtrop)ijp−fcΔ(Vion)ijp

此处的接受机钟差打错 了应该为ΔVij(t1)ΔVij(t1)接收机钟差不应该挂上卫星的,电离层和对流程延迟是跟卫星接收机两者都有关系
式子改为
Δφpij~=fcΔρptij+fΔVtijΔNpijfcΔ(Vtrop)pijfcΔ(Vion)pijΔφijp~=fcΔρtijp+fΔVtij−ΔNijp−fcΔ(Vtrop)ijp−fcΔ(Vion)ijp

在短基线中,接收机间求差对于卫星星历误差??、电离层延迟、对流层延迟等的影响也可得以削弱,短基线中较为明显。

2. 在接收机和卫星间求二次差

t1t1时刻卫星间求差:qpq−p消除接收机钟差

Δφ~qijΔφ~pij=fc[Δρqij(t1)Δρpij(t1)](ΔNqijΔNpij)fc(Δ(Vtrop)qijΔ(Vtrop)qij)fc(Δ(Vion)qijΔ(Vion)pij)Δφ~ijq−Δφ~ijp=fc[Δρijq(t1)−Δρijp(t1)]−(ΔNijq−ΔNijp)−fc(Δ(Vtrop)ijq−Δ(Vtrop)ijq)−fc(Δ(Vion)ijq−Δ(Vion)ijp)

感叹下,提出GPS这个想法的科学家真的好强
可以用双差公式记为
Δφ~pqij=fcΔρpqij(t1)ΔNpqijfcΔVtroppqijfcΔ(Vion)pqijΔφ~ijpq=fcΔρijpq(t1)−ΔNijpq−fcΔVtropijpq−fcΔ(Vion)ijpq

该观测方程为t1t1时刻在接收机和卫星间求二次差后所得到的双差观测值。
实际工作中,进一步在卫星间求双差是采用下列方式进行的:
选择视场中可观测时间较长、高度角又较大的一颗卫星作为基准星,然后将期许个卫星的单差观测方程分别与基准星的单差观测方程相减,组成双差观测方程。在每个观测历元中,双差观测方程的数量均比单差方程数少一个,但与次同时,该历元接收机相对钟差参数也已被消除。
求差法的优点,在保持原有精度的情况下,可大大减少未知数个数,从而大幅度减少数据处理的工作量。主要是与非差法作为参照对象。

3. 在接收机、卫星和观测历元间求三次差

测站i,jt2pqt2t1:i,j在历元t2对卫星p、q进行同步观测后,t2−t1相减得:

Δφ~pqij(t1,t2)=fcΔρpqij(t1,t2)fcΔ(Vtrop)pqij(t1,t2)fcΔ(Vion)pqij(t1,t2)Δφ~ijpq(t1,t2)=fcΔρijpq(t1,t2)−fcΔ(Vtrop)ijpq(t1,t2)−fcΔ(Vion)ijpq(t1,t2)

在三差观测方程中,整周模糊度参数ΔNpqijΔNijpq已经被消除,因而只剩下三个位置Δx,DeltaYDeltaZΔx,、DeltaY、DeltaZ三个未知参数。
三差法的缺点:
(1)双差观测值比三差观测值只多了7个参数,对于计算机来说相差不多
(2)这一点比较重要,因为由于各种误差的影响,三差解一般为实数解书上讲的是未做取整和回代等工作,故三差解是与浮点解相对的???。 三差方程的几何强度也较差。??
对于上述描述看的不是特别懂,希望学习到后面这些问题能迎刃而解。

4. 求差法的缺点

(1)数据利用的效率低。好的观测值与出了问题的数据配对时会导致数据无法使用。
(2)双差观测值是使用一颗卫星与其他卫星一起求差,故双差观测值之间都包含基准星,故这些双侧观测值为相关观测值。若令单差观测值的方差为σ2σ2,则双差观测值的协方差阵显然为:

Ri=σ22111121111211112Ri=σ2[2111121111211112]

双差观测值的相关性数据处理增加了难度。若不顾及其相关性。仍将其认为独立值(实验结果证明仍能取得较好的结果),但是损害了理论上的严密性。
(3)解的通用性比较差。由于与用户无直接关联的一些参数直接被消去,故无法对其他类型的用户提供使用。非差发能够获得多余参数,供其他类型 的用户使用。

求差法也适用于伪距观测方程。

版权声明:本文为qq_39600809原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/qq_39600809/article/details/79602333

智能推荐

差分数组学习笔记

1.定义 对于已知有n个元素的离线数列d,我们可以建立记录它每项与前一项差值的差分数组f:显然,f[1]=d[1]-0=d[1];对于整数i∈[2,n],我们让f[i]=d[i]-d[i-1]。 2.应用 区间加减操作:O(1) 假如现在对数列中区间[L,R]上的数加上x,令f[L]+=x,f[R+1]-=x即可。 询问每一项的值:O(n) 求差分数组前缀和即可 数组前缀和:O(n) 第...

差分约束 学习笔记

差分约束 差分约束就是:给你 n 个变量,m 个未知数,形成这样几组的不等式:xix_ixi​ - yjy_jyj​ <= ckc_kck​; 然后让你求出一组解,使得所有约束条件都满足; 怎么解呢?可以把上面的式子变形一下:xix_ixi​ <= yjy_jyj​ + ckc_kck​; 这是不是和图论最短路的式子非常相似呢? 所以可以连接 yiy_iyi​ 到 xix_ixi​ 的...

差分数组 学习笔记

我想想我第一次接触差分数组算法,是在洛谷上找树状数组线段树题的时候,遇到的 P1083 借教室。看到题解里有一篇号称更好理解也更好实现的算法,那我为什么不学? 差分数组 首先要说的是,差分数组的思想与前缀和的思想是密不可分的,不如说差分数组就是前缀和数组的逆运算 我们平时是这么思考前缀和数组的? 用sum[i]来存储前i个数的和,然后用sum[r]-sum[l-1]来表示 l ~ r 之间所有数的...

【深度学习笔记】残差

结构 ResNet block有两种,一种两层结构,一种三层结构 我们需要求解的是映射:H(x)。 残差结构将这个问题转换为求解网络的残差映射函数,也就是F(x),其中F(x) = H(x)-x。 从而,问题转换成了 H(x) = F(x)+x。 原因与分析 如果是采用一般的卷积神经网络(没有跳跃连接),原先需要求解的是 H(x) = F(x) 。 那么假设在网络达到某一个深度时,网络已经达到最优...

学习笔记:差分约束

上讲习题 AcWing 1165 从一个字符串的前两个字符向最后两个字符连一条权为这个字符串长度的边,表示用这个字符串对平均值的贡献,一个路径就是一个接龙,一个环就是题目要求的。接下来思考解决这个题,本题又是 a b \cfrac{a}{b} ba​的形式,可以 01 01 01分数规划。设 V V V为用的字符串的集合, ∣ X ∣ |X| ∣X∣表示集合 X X X的大小, l e n X l...

猜你喜欢

R语言学习笔记NO.2 数据处理的基础操作_变量/观测的增改选排

我们先创建一个基础数据框df,用madata记录初始值. 初值备份是一个良好的习惯,这可避免你在误操作后丢失初值. 1.变量的新增 1.1新增变量与已有变量无关 对于上述数据框df,如果想增加一个变量age<-c(32,10,65,150,46): 1.2新增与已有变量相关 如果你想将country和manger变量合并为一个变量:country_manager:那么首先要新增country...

Layout (线性差分)

Like everyone else, cows like to stand close to their friends when queuing for feed. FJ has N (2 <= N <= 1,000) cows numbered 1..N standing along a straight line waiting for feed. The cows are s...

[差分约束系统] 差分约束系统_学习笔记

文章目录 讲解 T1 小k的农场 T2 Candy T3 [HNOI2005]狡猾的商人 讲解 对于一堆不等式 发现多元 于是不能用线性规划解决 考虑元之间的关系: 有向,感觉拓扑 就很像以前做过的拓扑排序的一道题 于是我们寻思能不能转图论 发现可以,就是求最短路 按不等式建DAG 遇到环呢,炸了呗 所以要用SPFA 用小于号 T1 小k的农场 传送门 fa−fb≥c⇒...

测量平差之附有限制条件的条件平差(概括平差模型)

接下来要写的可看成是对前面讲过的四种经典平差方法的概括模型。 从四种基本平差算法的函数模型来看,主要包括如下两种类型的条件方程: {F(L^,X^)=0ϕ(X^)=0{F(L^,X^)=0ϕ(X^)=0 第一种方程中同时含有观测值和未知参数,称为一般条件方程;第二种方程则只含有未知参数而无观测值,称为限制条件方程。 一般而言,对于任何一个平差问题,设想观测值个数为nn,必要观测数为tt,则多余观测...

第三部分 数据结构 第 10 章 基本数据结构

  数学中的集合是不变的,而由算法操作的集合却在整个过程中能增大、缩小或发生其他变化,称这样的集合是动态的。   动态集合中的元素   在动态集合的典型实现中,每个元素都由一个对象来表示,如果有一个指向对象的指针,就能对其各个属性进行检查和操作。一些类型的动态集合假定对象中的一个属性为标识关键字。   动态集合上的操作   分为两类:简单返回有关集合信息的查询操作和改变集合的修改操作。下面列出一些...