CN107290762A - 基于载波三次差分的室外gnss绝对天线相位中心标校方法 - Google Patents
基于载波三次差分的室外gnss绝对天线相位中心标校方法 Download PDFInfo
- Publication number
- CN107290762A CN107290762A CN201610841265.3A CN201610841265A CN107290762A CN 107290762 A CN107290762 A CN 107290762A CN 201610841265 A CN201610841265 A CN 201610841265A CN 107290762 A CN107290762 A CN 107290762A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- antenna
- gnss
- epsiv
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 230000000694 effects Effects 0.000 claims abstract description 9
- 238000006243 chemical reaction Methods 0.000 claims description 10
- 238000012937 correction Methods 0.000 claims description 8
- 230000009466 transformation Effects 0.000 claims description 8
- 230000008859 change Effects 0.000 claims description 5
- 238000006073 displacement reaction Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000013519 translation Methods 0.000 claims description 3
- 239000005433 ionosphere Substances 0.000 claims 1
- 230000007704 transition Effects 0.000 claims 1
- 239000005436 troposphere Substances 0.000 claims 1
- 238000005259 measurement Methods 0.000 abstract description 5
- 230000001629 suppression Effects 0.000 abstract description 2
- 238000012360 testing method Methods 0.000 description 12
- 238000004364 calculation method Methods 0.000 description 7
- 230000009897 systematic effect Effects 0.000 description 2
- 238000013480 data collection Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
Landscapes
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种基于载波三次差分的室外GNSS绝对天线相位中心标校方法，首先对精密机械臂大角度快速运动时的短基线进行双差解算，然后在模糊度固定之后进行历元间三次差分。在极短的历元间隔可消除参考站天线相位中心误差，并且多路径效应引起的测量误差也得到极大削弱，从而实现GNSS天线绝对相位中心补偿PCO和相位中心变化PCV的精确建模。本发明利用GNSS室外实测数据就能够将待测天线相位中心误差从相对定位残差中分离开来，且对低高度角观测模式下多路径效应抑制效果明显。从而实现1mm精度的PCO和PCV标校，对GPS/BDS/GLONASS/Galileo系统高精度定位有明显的改善效果。
Description
技术领域
本发明属于天线测量技术与卫星导航定位技术领域，具体涉及一种基于载波三次差分的室外GNSS绝对天线相位中心标校方法。
背景技术
天线相位中心即天线接收信号的电气中心，其空间位置在出厂时往往不在天线的几何中心上。一般选取天线底部与天线中轴的交点作为参考点(称天线参考点，英文antenna reference point：ARP)给出其电气中心，这个几何偏差值称为天线相位中心偏差(antenna phase center offset：PCO)。从电气中心出发的信号在理想情况下是一个球形的等相面。然而实际制造的天线不是一个球形等相面，而是一个被扭曲的不规则曲面。在这种情况下，从不同高度角和方位角测得的距离产生系统性的测量偏差，叫天下相位中心变化(antenna phase center variation：PCV)。实践证明，由天线产生的PCO和PCV误差是影响用户定位精度的关键系统误差源。高精度定位需要高精度的天线相位中心改正模型进行精确修正。因此，天线需要做天线的相位中心标校，并将标校值提供给用户使用。
GNSS(Global Navigation Satellite System，全球导航卫星系统)天线相位中心的校正或标定目前主要有微波暗室校正法和利用室外实测GNSS信号的相对定位法两种。微波暗室法采用的是在较为理想的环境下的GNSS模拟信号进行标定，受干扰较少，但这种方法对设备要求较高，测试成本也较高，并且其模拟的观测环境与室外观测环境可能存在一定的差异，导致校正结果可行度不高。而室外实测GNDD信号的相对定位法要求参考站的天线相位中心信息事先已知，得到的待测天线结果是相对于参考天线的相对相位中心模型值，不能得到与参考站无关的绝对天线相位中心模型值。同时，达到1mm标校精度对影响GNSS信号各种信号特别是多路径效应等多种误差的消除提出了非常高的要求。
发明内容
本发明要解决的技术问题是：
(1)常规室外GNSS观测无法获得高精度的接收机绝对天线相位中心校正模型；
(2)低高度角观测数据受多路径效应影响严重，GNSS双差仍然无法消除其影响，超过1mm。
为了解决上述技术问题，本发明提供了一种基于载波三次差分的室外GNSS绝对天线相位中心标校方法。
本发明所采用的技术方案是：一种基于载波三次差分的室外GNSS绝对天线相位中心标校方法，其特征在于，包括以下步骤：
步骤1：利用两个接收机组成超短基线；
所述接收机均配置有GNSS天线，其中一个GNSS天线安装在参考站上，另一个GNSS天线安装在移动站配置的机械臂末端上；调整转动所述机械臂观测得到不同时段的多条基线，通过基线解算得到机械臂的GNSS坐标值，再通过读取对应时段的机械臂坐标值，获得由GNSS坐标值与机械臂坐标值组成的坐标对；通过解算这些坐标对，得到GNSS坐标系与机械臂坐标系之间的转换参数，实现坐标系的统一；
步骤2：固定参考站，将待测天线安装于参考站配置的机械臂末端上，标记待测天线上方位角的起算方向，记录起算方向与机械臂的相对位置；旋转和倾斜机械臂，记录高采样GNSS观测值；
步骤3：将参考站与待测天线组成经典的二次差分(双差)观测模式，利用机械臂提供的位移变化值进行基线解算，同时获得包含已知天线和待测天线相位中心误差在内的双差观测值残差；
步骤4：利用步骤3得到的双差观测值残差，对相邻历元的双差观测残差作历元间三次差，获得三次差分(三差)观测值；以指定的天线北为起算方向，以待测天线参考点为原点，建立天线本体站心坐标系，即北/东/高坐标系；在此坐标系进行绝对天线相位偏差参数建模；通过历元间三差观测值同时解算出北、东、高方向分量PCO(N)、PCO(E)和PCO(U)；
步骤5：将解算的北、东、高方向分量PCO(N)、PCO(E)和PCO(U)回带入步骤4的三差观测值中，得到消除PCO分量的三差观测残差；依据步骤4建立的天线本体站心坐标系计算高度角和方位角，对观测值残差进行分段线性拟合或者球冠谐函数拟合计算得到关于高度角和方位角的相对天线相位中心变化值。
作为优选，步骤1的具体实现包括以下子步骤：
步骤1.1：将参考站与移动站组成基线，移动站标记位置1；通过相对定位得到位置1的GNSS坐标，同时记录机械臂末端天线在机械臂坐标系下的坐标，两者组成坐标对1；
步骤1.2：调整机械臂位置到位置2，将参考站与位置2处的移动站组成基线，通过相对定位得到位置2的GNSS坐标，同时记录机械臂末端天线在机械臂坐标系下的坐标，两者组成坐标对2；
步骤1.3：通过调整机械臂的位置，得到N个坐标对，N≥3；
步骤1.4：采用经典的7参数坐标转换模型得到GNSS坐标系与机械臂坐标系之间的转换参数。
作为优选，步骤1中所述GNSS坐标系与机械臂坐标系之间的转换关系如下：
XGNSS＝X0+(1-k)RR2GXR (1)
其中，XGNSS、XR分别为同一个位置在GNSS坐标系和机械臂坐标系下的坐标值；X0＝(△x，△y，△z)为坐标平移量；(εx,εy,εz)为欧拉角，k为单位尺度因子；RR2G为机械臂坐标系到GNSS坐标系转换矩阵，(Δx,Δy,Δz,εx,εy,εz,k)是待求的7个坐标转换参数。
作为优选，步骤2的具体实现包括以下子步骤：
步骤2.1：调整待测天线，使待测天线与机械臂保持水平，将事先标记的天线指北方向与机械臂的X坐标轴平行，此时待测天线位置标记为位置1；
步骤2.2：水平面内转动机械臂到某一指定的角度，此时待测天线位置标记为位置2；
步骤2.3：竖直面内转动机械臂到某一指定的角度，此时待测天线位置标记为位置3；
步骤2.4：通过机械臂的转动和倾斜，GNSS观测数据采集6至10小时结束。
作为优选，步骤3的具体实现包括以下子步骤：
步骤3.1：在任意时刻t，计算机械臂末端天线在GNSS坐标系下的三维坐标；
步骤3.2：在时刻t，利用公式(1)将机械臂坐标转换为GNSS所在坐标系，对待测天线进行线性化，组成仅包含双差模糊度参数的双差观测方程；
步骤3.3：利用经典模糊度固定方法，固定时刻t对应的双差模糊度，并计算得到双差观测残差；该残差包含参考站天线和待测天线的相位中心误差。
作为优选，步骤3.3的具体实现过程是：
记参考站和移动站的单频载波相位观测方程(这里忽略观测噪声)：
式中，代表单频载波相位观测值，单位为周；p代表卫星号；数字1和2代表测站；ti代表观测历元；ρ代表站星间的几何距离；c代表光速；dtr代表接收机钟差；dt代表卫星钟差；N代表整周模糊度；λ代表对应频率的波长；diono代表电离层延迟；dtrop代表对流层延迟；dmp代表多路径误差；dpcv代表PCV误差；a和e代表接收到的信号的方位角和高度角；
对参考站和移动站的非差观测方程分别进行线性化后，相近历元时刻ti,tj组成双差观测方程，待估参数为模糊度：
对(3)式中的载波相位模糊度采用模糊度固定方法固定，得到的双差观测值残差包含参考站天线和待测天线的相位中心误差、残余多路径误差和双差观测值噪声。
作为优选，步骤4中所述利用步骤3得到的双差观测值残差，对相邻历元的双差观测残差作历元间三次差，获得简称三差观测值；
历元间三次差将消去参考站天线相位中心误差和大大削弱多路径效应：
将所有已知量代入后，可以得到最终的、仅包含待测天线相位中心误差的三差观测残差Tdres：
作为优选，步骤5中球谐函数来拟合PCV的公式如下：
式中，Pnm(cosz)为勒让德函数；n，m，anm，bnm分别为阶数、次数和待估参数；α,z,fi分别为方位角、天顶距(与高度角互余)和频率；
计算PCV时，首先将公式(6)左边ΔPCV用计算PCO三个分量时得到的三差观测残差值代替，Pnm(cosz)勒让德系数、cos z、cos mα、sin mα等通过已知的天顶距和方位角一一计算；然后将公式(6)线性化后通过常规的最小二乘法拟合得到带估参数anm，bnm；最后根据带估参数anm，bnm和公式(6)得到任意高度角和方位角的PCV值。
本发明基于可倾斜和旋转的高精度机械臂采集大量不同方位角和高度角的实测GNSS观测值，不需要双差观测值拟合参考站虚拟观测值，不用考虑参考站与流动站共用一个接收机的复杂工艺，在室外就可以实现GNSS接收机绝对天线相位中心的标定。同时对低高度角观测数据的多路径效应抑制效果明显，标定精度高，操作简单，检校成本低，室外测试环境更接近使用环境，能够实现GNSS天线相位中心标定的精确标定。
具体实施方式
为了便于本领域普通技术人员理解和实施本发明，下面结实施例对本发明作进一步的详细描述，应当理解，此处所描述的实施示例仅用于说明和解释本发明，并不用于限定本发明。
本发明提供的一种基于三次差分法的室外GNSS绝对天线相位中心标校方法，包括以下步骤：
步骤1：选用两个配备GNSS天线(相位中心改正模型已知)的接收机组成超短基线(＜10m)。其中一个接收机天线固定在参考站，另一个接收机天线安装在机械臂末端(移动站)。调整机械臂移动天线，可观测得到不同时段的多条基线(不低于3个位移点)，通过基线解算可得到机械臂末端的GNSS坐标值(即天线的GNSS坐标值)，再通过计算机读取对应时段的机械臂末端的机械臂本体坐标值。通过解算这些坐标对，得到GNSS坐标系与机械臂本体坐标系之间的转换参数，实现坐标系的统一。
其中坐标转换参数计算过程如下：
第一步：将参考站与移动站(标记位置1)组成基线，通过相对定位得到位置1的GNSS坐标，同时记录机械臂末端在机械臂全局坐标系下的坐标，两者组成坐标对1；
第二步：调整机械臂末端的位置(位置2)，将参考站与位置2处的移动站组成基线，通过相对定位得到位置2的GNSS坐标，同时记录机械臂末端在机械臂全局坐标系下的坐标，两者组成坐标对2；
第三步：采用与第一步和第二步类似的方式，通过调整机械臂末端的位置，得到N个坐标对(N≥3)；
第四步：采用经典的7参数坐标转换模型得到GNSS坐标系与机械臂本体坐标系之间的转换参数；
具体转换关系如下：
XGNSS＝X0+(1-k)RR2GXR (1)
其中，XGNSS、XR分别为同一个位置在GNSS坐标系和机械臂坐标系下的坐标值；X0＝(△x，△y，△z)为坐标平移量；(εx,εy,εz)为欧拉角，k为单位尺度因子；RR2G为机械臂坐标系到GNSS坐标系转换矩阵:
式中，(Δx,Δy,Δz,εx,εy,εz,k)是待求的7个坐标转换参数。
步骤2：固定配备任意接收机天线的参考站，将待测天线安装于机械臂末端，并在天线上标记好方位角的起算方向(即“天线指北方向”)并记下这个方向与机械臂的相对位置。快速旋转和倾斜机械臂，同时记录高采样GNSS观测值(采样率≥1Hz)；
其中观测方法如下：
第一步：调整待测天线，使待测天线与机械臂保持水平，将事先标记的天线指北方向与机械臂的X坐标轴平行(待测天线方位角起算方向)，此时待测天线位置标记为位置1；
第二步：水平面内转动机械臂到某一指定的角度(旋转角度任意)，此时待测天线位置标记为位置2；
第三步：竖直面内转动机械臂到某一指定的角度(建议倾斜角度小于等于30°)，此时待测天线位置标记为位置3；
第四步：按第三、四步方式进行机械臂的转动和倾斜，GNSS观测数据采集6至10小时结束；
步骤3：根据步骤2，将参考站与待测天线组成经典的双差观测模式，利用机械臂提供的位移变化值(步骤1实现)进行基线解算。解算参数仅包含双差模糊度参数，采用模糊度固定算法(例如利用LAMBDA)固定模糊度到整数。同时获得包含已知天线和待测天线相位中心误差在内的双差观测值残差。
记参考站和移动站的单频载波相位观测方程：
式中，代表单频载波相位观测值，单位为周；p代表卫星号；数字1和2代表测站；ti代表观测历元；ρ代表站星间的几何距离；c代表光速；dtr代表接收机钟差；dt代表卫星钟差；N代表整周模糊度；λ代表对应频率的波长；diono代表电离层延迟；dtrop代表对流层延迟；dmp代表多路径误差；dpcv代表PCV误差；a和e代表接收到的信号的方位角和高度角。
对参考站和移动站的非差观测方程分别进行线性化后，相近历元时刻ti,tj组成双差观测方程，待估参数为模糊度：
对(3)式中的载波相位模糊度采用LAMBDA等方法固定，得到的双差观测值残差包含参考站和待测天线的相位中心误差、残余多路径误差。
步骤4：利用步骤3得到的同步双差观测值残差，对相邻历元的双差观测残差作历元间三次差。以步骤2指定的天线北为起算方向，以待测天线参考点为原点，建立天线本体站心坐标系(北/东/高坐标系)，在此坐标系进行绝对天线相位偏差(PhaseCenter Offset简称PCO)参数建模。通过历元间三差残差可同时解算出北、东、高方向分量PCO(N)、PCO(E)和PCO(U)。
1)若历元间的时间间隔很短(＜10s)，可以认为在极短的时间内卫星相对于参考站天线的高度角和方位角变化非常小并且两个历元的观测环境是一致的，具有相同的多路径误差。将模糊度固定为整数后，组历元间三差观测方程，将消去参考站天线相位中心误差和大大削弱多路径效应：
将所有已知量代入后，可以得到最终的、仅包含待测天线相位中心误差的三差观测残差Tdres：
2)以步骤2指定的天线北为起算方向，以待测天线参考点为原点，建立天线本体站心坐标系(北/东/高坐标系)，在此坐标系进行绝对天线相位偏差(PhaseCenter Offset简称PCO)参数建模。通过历元间三差残差可同时解算出待测天线北、东、高方向分量PCO(N)、PCO(E)和PCO(U)。
步骤5：根据步骤4，将解算的三个分量的PCO估值带入步骤4的观测方程，得到消除PCO分量新的三差观测残差。依据步骤4建立的坐标系计算高度角和方位角，对观测值残差进行分段线性拟合或者球冠谐函数拟合计算得到关于高度角和方位角的相对天线相位中心变化(Phase Center Variation简称PCV)值。用球谐函数来拟合PCV的公式如下：
式中，Pnm(cosz)为勒让德函数；n，m，anm，bnm分别为阶数、次数和待估参数；α,z,fi分别为方位角、天顶距和频率。实际计算PCV时，首先将公式(6)左边ΔPCV用计算PCO三个分量时得到的三差观测残差值(参考公式(5))代替，Pnm(cosz)勒让德系数、cos z、cos mα、sinmα等通过已知的天顶距和方位角一一计算；然后将公式(6)线性化后通过常规的最小二乘法(拟合)得到带估参数anm，bnm；最后根据带估参数anm，bnm和公式(6)得到任意天顶距(90°-高度角)和方位角的PCV值。
本发明通过超短基线观测消除大部分公共误差，借助高精度的机械臂的快速旋转和倾斜，在极短的时间内对双差观测值进行历元间做差，实现待测接收机天线绝对相位中心补偿(PCO)和相位中心变化(PCV)与参考站天线的分离。最后通过最小二乘算法解算PCO并对观测残差进行拟合计算PCV。本发明能在室外操作，且能达到1mm的PCO和PCV的校正精度，适用于能跟踪GPS/BDS/GLONASS/Galileo等卫星导航系统信号的接收机PCO和PCV的精确标定，从而可以消除由于接收机PCO和PCV的不精确引入的系统性测量误差，对GPS/BDS/GLONASS/Galileo系统高精度定位有明显的改善效果。
应当理解的是，本说明书未详细阐述的部分均属于现有技术。
应当理解的是，上述针对较佳实施例的描述较为详细，并不能因此而认为是对本发明专利保护范围的限制，本领域的普通技术人员在本发明的启示下，在不脱离本发明权利要求所保护的范围情况下，还可以做出替换或变形，均落入本发明的保护范围之内，本发明的请求保护范围应以所附权利要求为准。
Claims (8)
1.一种基于载波三次差分的室外GNSS绝对天线相位中心标校方法，其特征在于，包括以下步骤：
步骤1：利用两个接收机组成超短基线；
所述接收机均配置有GNSS天线，其中一个GNSS天线安装在参考站上，另一个GNSS天线安装在移动站配置的机械臂末端上；调整转动所述机械臂观测得到不同时段的多条基线，通过基线解算得到机械臂的GNSS坐标值，再通过读取对应时段的机械臂坐标值，获得由GNSS坐标值与机械臂坐标值组成的坐标对；通过解算这些坐标对，得到GNSS坐标系与机械臂坐标系之间的转换参数，实现坐标系的统一；
步骤2：固定参考站，将待测天线安装于参考站配置的机械臂末端上，标记待测天线上方位角的起算方向，记录起算方向与机械臂的相对位置；旋转和倾斜机械臂，记录高采样GNSS观测值；
步骤3：将参考站与待测天线组成经典的二次差分观测模式，利用机械臂提供的位移变化值进行基线解算，同时获得包含已知天线和待测天线相位中心误差在内的双差观测值残差；
步骤4：利用步骤3得到的载波双差观测值残差，对相邻历元的载波双差观测残差作历元间三次差，获得载波三次差分观测值，即三差观测值；以指定的天线北为起算方向，以待测天线参考点为原点，建立天线本体站心坐标系，即北/东/高坐标系；在此坐标系进行绝对天线相位偏差参数建模；通过历元间三差观测值同时解算出北、东、高方向分量PCO(N)、PCO(E)和PCO(U)；
步骤5：将解算的北、东、高方向分量PCO(N)、PCO(E)和PCO(U)回带入步骤4的三差观测值中，得到消除PCO分量的三差观测残差；依据步骤4建立的天线本体站心坐标系计算高度角和方位角，对观测值残差进行分段线性拟合或者球冠谐函数拟合计算得到关于高度角和方位角的相对天线相位中心变化值。
2.根据权利要求1所述的基于载波三次差分的室外GNSS绝对天线相位中心标校方法，其特征在于，步骤1的具体实现包括以下子步骤：
步骤1.1：将参考站与移动站组成基线，移动站标记位置1；通过相对定位得到位置1的GNSS坐标，同时记录机械臂末端天线在机械臂坐标系下的坐标，两者组成坐标对1；
步骤1.2：调整机械臂位置到位置2，将参考站与位置2处的移动站组成基线，通过相对定位得到位置2的GNSS坐标，同时记录机械臂末端天线在机械臂坐标系下的坐标，两者组成坐标对2；
步骤1.3：通过调整机械臂的位置，得到N个坐标对，N≥3；
步骤1.4：采用经典的7参数坐标转换模型得到GNSS坐标系与机械臂坐标系之间的转换参数。
3.根据权利要求1或2所述的基于载波三次差分的室外GNSS绝对天线相位中心标校方法，其特征在于，步骤1中所述GNSS坐标系与机械臂坐标系之间的转换关系如下：
XGNSS＝X0+(1-k)RR2GXR (1)
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<msub>
<mi>R</mi>
<mrow>
<mi>R</mi>
<mn>2</mn>
<mi>G</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>cos&epsiv;</mi>
<mi>y</mi>
</msub>
<msub>
<mi>cos&epsiv;</mi>
<mi>z</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>cos&epsiv;</mi>
<mi>y</mi>
</msub>
<msub>
<mi>sin&epsiv;</mi>
<mi>z</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>sin&epsiv;</mi>
<mi>y</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>cos&epsiv;</mi>
<mi>x</mi>
</msub>
<msub>
<mi>sin&epsiv;</mi>
<mi>z</mi>
</msub>
<mo>+</mo>
<msub>
<mi>sin&epsiv;</mi>
<mi>x</mi>
</msub>
<msub>
<mi>sin&epsiv;</mi>
<mi>y</mi>
</msub>
<msub>
<mi>cos&epsiv;</mi>
<mi>z</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>cos&epsiv;</mi>
<mi>x</mi>
</msub>
<msub>
<mi>cos&epsiv;</mi>
<mi>z</mi>
</msub>
<mo>+</mo>
<msub>
<mi>sin&epsiv;</mi>
<mi>x</mi>
</msub>
<msub>
<mi>sin&epsiv;</mi>
<mi>y</mi>
</msub>
<msub>
<mi>sin&epsiv;</mi>
<mi>z</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>sin&epsiv;</mi>
<mi>x</mi>
</msub>
<msub>
<mi>cos&epsiv;</mi>
<mi>y</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>sin&epsiv;</mi>
<mi>x</mi>
</msub>
<msub>
<mi>sin&epsiv;</mi>
<mi>z</mi>
</msub>
<mo>+</mo>
<msub>
<mi>cos&epsiv;</mi>
<mi>x</mi>
</msub>
<msub>
<mi>sin&epsiv;</mi>
<mi>y</mi>
</msub>
<msub>
<mi>cos&epsiv;</mi>
<mi>z</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>sin&epsiv;</mi>
<mi>x</mi>
</msub>
<msub>
<mi>cos&epsiv;</mi>
<mi>z</mi>
</msub>
<mo>+</mo>
<msub>
<mi>cos&epsiv;</mi>
<mi>x</mi>
</msub>
<msub>
<mi>sin&epsiv;</mi>
<mi>y</mi>
</msub>
<msub>
<mi>sin&epsiv;</mi>
<mi>z</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>cos&epsiv;</mi>
<mi>x</mi>
</msub>
<msub>
<mi>cos&epsiv;</mi>
<mi>y</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中，XGNSS、XR分别为同一个位置在GNSS坐标系和机械臂坐标系下的坐标值；X0＝(△x，△y，△z)为坐标平移量；(εx,εy,εz)为欧拉角，k为单位尺度因子；RR2G为机械臂坐标系到GNSS坐标系转换矩阵，(Δx,Δy,Δz,εx,εy,εz,k)是待求的7个坐标转换参数。
4.根据权利要求1所述的基于载波三次差分的室外GNSS绝对天线相位中心标校方法，其特征在于，步骤2的具体实现包括以下子步骤：
步骤2.1：调整待测天线，使待测天线与机械臂保持水平，将事先标记的天线指北方向与机械臂的X坐标轴平行，此时待测天线位置标记为位置1；
步骤2.2：水平面内转动机械臂到某一指定的角度，此时待测天线位置标记为位置2；
步骤2.3：竖直面内转动机械臂到某一指定的角度，此时待测天线位置标记为位置3；
步骤2.4：通过机械臂的转动和倾斜，GNSS观测数据采集6至10小时结束。
5.根据权利要求3所述的基于载波三次差分的室外GNSS绝对天线相位中心标校方法，其特征在于，步骤3的具体实现包括以下子步骤：
步骤3.1：在任意时刻t，计算机械臂末端天线在GNSS坐标系下的三维坐标；
步骤3.2：在时刻t，利用公式(1)将机械臂坐标转换为GNSS所在坐标系，对待测天线进行线性化，组成仅包含双差模糊度参数的双差观测方程；
步骤3.3：利用经典模糊度固定方法，固定时刻t对应的双差模糊度，并计算得到双差观测残差；该残差包含参考站天线和待测天线的相位中心误差。
6.根据权利要求5所述的基于载波三次差分的室外GNSS绝对天线相位中心标校方法，其特征在于，步骤3.3的具体实现过程是：
记参考站和移动站的单频载波相位观测方程：
式中，代表单频载波相位观测值，单位为周；p代表卫星号；数字1和2代表参考站和移动站站；ti代表观测历元；ρ代表站星间的几何距离；c代表光速；dtr代表接收机钟差；dt代表卫星钟差；N代表整周模糊度；λ代表对应频率的波长；diono代表电离层延迟；dtrop代表对流层延迟；dmp代表多路径误差；dpcv代表PCV误差；a和e代表接收到的信号的方位角和高度角；
对参考站和移动站的非差观测方程分别进行线性化后，相近历元时刻ti,tj组成双差观测方程，待估参数为模糊度：
对(3)式中的载波相位模糊度采用模糊度固定方法固定，得到的双差观测值残差包含参考站天线和待测天线的相位中心误差、残余多路径误差和双差观测值噪声。
7.根据权利要求6所述的基于载波三次差分的室外GNSS绝对天线相位中心标校方法，其特征在于：步骤4中所述利用步骤3得到的双差观测值残差，对相邻历元的双差观测残差作历元间三次差，获得简称三差观测值；
历元间三次差将消去参考站天线相位中心误差和大大削弱多路径效应：
将所有已知量代入后，可以得到最终的、仅包含待测天线相位中心误差的三差观测残差Tdres：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>T</mi>
<mi>D</mi>
<mi>r</mi>
<mi>e</mi>
<mi>s</mi>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>d</mi>
<mrow>
<mi>p</mi>
<mi>c</mi>
<mi>v</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>a</mi>
<mn>2</mn>
<mi>q</mi>
</msubsup>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
<mo>,</mo>
<msubsup>
<mi>e</mi>
<mn>2</mn>
<mi>q</mi>
</msubsup>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>d</mi>
<mrow>
<mi>p</mi>
<mi>c</mi>
<mi>v</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>a</mi>
<mn>2</mn>
<mi>p</mi>
</msubsup>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
<mo>,</mo>
<msubsup>
<mi>e</mi>
<mn>2</mn>
<mi>p</mi>
</msubsup>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>d</mi>
<mrow>
<mi>p</mi>
<mi>c</mi>
<mi>v</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>a</mi>
<mn>2</mn>
<mi>q</mi>
</msubsup>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
<mo>,</mo>
<msubsup>
<mi>e</mi>
<mn>2</mn>
<mi>q</mi>
</msubsup>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>-</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>d</mi>
<mrow>
<mi>p</mi>
<mi>c</mi>
<mi>v</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>a</mi>
<mn>2</mn>
<mi>p</mi>
</msubsup>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
<mo>,</mo>
<msubsup>
<mi>e</mi>
<mn>2</mn>
<mi>p</mi>
</msubsup>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
8.根据权利要求1所述的基于载波三次差分的室外GNSS绝对天线相位中心标校方法，其特征在于，步骤5中球谐函数来拟合PCV的公式如下：
<mrow>
<mi>&Delta;</mi>
<mi>P</mi>
<mi>C</mi>
<mi>V</mi>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>,</mo>
<mi>z</mi>
<mo>,</mo>
<msub>
<mi>f</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>N</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>M</mi>
</munderover>
<msub>
<mi>P</mi>
<mrow>
<mi>n</mi>
<mi>m</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>cos</mi>
<mi> </mi>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mrow>
<mi>n</mi>
<mi>m</mi>
</mrow>
</msub>
<mi>cos</mi>
<mi> </mi>
<mi>m</mi>
<mi>&alpha;</mi>
<mo>+</mo>
<msub>
<mi>b</mi>
<mrow>
<mi>n</mi>
<mi>m</mi>
</mrow>
</msub>
<mi>sin</mi>
<mi> </mi>
<mi>m</mi>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
式中，Pnm(cosz)为勒让德函数；n，m，anm，bnm分别为阶数、次数和待估参数；α,z,fi分别为方位角、天顶距和频率；
计算PCV时，首先将公式(6)左边ΔPCV用计算PCO三个分量时得到的三差观测残差值代替，Pnm(cosz)勒让德系数、cos z、cos mα、sin mα等通过已知的天顶距和方位角一一计算；然后将公式(6)线性化后通过最小二乘法拟合得到带估参数anm，bnm；最后根据带估参数anm，bnm和公式(6)得到任意天顶距和方位角的PCV值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610841265.3A CN107290762B (zh) | 2016-09-22 | 2016-09-22 | 基于载波三次差分的室外gnss绝对天线相位中心标校方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610841265.3A CN107290762B (zh) | 2016-09-22 | 2016-09-22 | 基于载波三次差分的室外gnss绝对天线相位中心标校方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107290762A true CN107290762A (zh) | 2017-10-24 |
| CN107290762B CN107290762B (zh) | 2018-09-07 |
Family
ID=60093351
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610841265.3A Active CN107290762B (zh) | 2016-09-22 | 2016-09-22 | 基于载波三次差分的室外gnss绝对天线相位中心标校方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107290762B (zh) |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109031359A (zh) * | 2018-06-11 | 2018-12-18 | 武汉大学 | 一种基于微波暗室的gnss接收天线绝对相位中心校正方法 |
| CN109061559A (zh) * | 2018-06-21 | 2018-12-21 | 同济大学 | 一种uwb基站天线相位中心偏差建模与改正的研究方法 |
| CN110658537A (zh) * | 2019-11-13 | 2020-01-07 | 深圳市鼎耀科技有限公司 | 微波暗室gnss接收机天线绝对相位中心高精度标定的方法 |
| CN112363186A (zh) * | 2019-12-27 | 2021-02-12 | 广东电网有限责任公司电力科学研究院 | 一种卫星天线相位中心参数的计算方法和装置 |
| CN113267794A (zh) * | 2021-07-20 | 2021-08-17 | 杭州中科微电子有限公司 | 一种基线长度约束的天线相位中心校正方法及装置 |
| CN113552599A (zh) * | 2021-07-07 | 2021-10-26 | 杭州中科微电子有限公司 | 基于信息融合的gnss接收机天线相位中心校正方法及装置 |
| CN113568011A (zh) * | 2021-07-16 | 2021-10-29 | 航天科工海鹰集团有限公司 | 一种ppk多径消除技术 |
| CN113703017A (zh) * | 2020-05-22 | 2021-11-26 | 中国人民解放军战略支援部队信息工程大学 | 一种卫星天线相位中心偏差计算方法及装置 |
| CN113721267A (zh) * | 2021-09-01 | 2021-11-30 | 航天科工通信技术研究院有限责任公司 | 一种基于双接收机载波相位差分的gnss欺骗干扰检测方法 |
| CN113885054A (zh) * | 2021-12-07 | 2022-01-04 | 广东皓行科技有限公司 | 天线位置确定方法及装置 |
| CN113899957A (zh) * | 2021-09-30 | 2022-01-07 | 武汉大学 | 基于双差相位观测值的室外天线绝对标定方法及系统 |
| CN115826008A (zh) * | 2023-02-22 | 2023-03-21 | 广州导远电子科技有限公司 | 天线位移距离范围的检测方法、装置、流动站及存储介质 |
| CN116542132A (zh) * | 2023-04-07 | 2023-08-04 | 武汉大学 | 结合球冠谐模型和神经网络的水汽数据校准与优化方法 |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP3667369B1 (en) * | 2018-12-14 | 2023-10-04 | Valeo Comfort and Driving Assistance | Positioning system for a land vehicle and method for computing high-precision gnss positions of a land vehicle |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102830406A (zh) * | 2012-08-30 | 2012-12-19 | 航天东方红卫星有限公司 | 一种gps天线相位中心变化对绝对定位精度的修正方法 |
| CN104502926A (zh) * | 2014-12-17 | 2015-04-08 | 武汉大学 | 一种基于精密机械臂的室外gnss接收机绝对天线相位中心校正方法 |
-
2016
- 2016-09-22 CN CN201610841265.3A patent/CN107290762B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102830406A (zh) * | 2012-08-30 | 2012-12-19 | 航天东方红卫星有限公司 | 一种gps天线相位中心变化对绝对定位精度的修正方法 |
| CN104502926A (zh) * | 2014-12-17 | 2015-04-08 | 武汉大学 | 一种基于精密机械臂的室外gnss接收机绝对天线相位中心校正方法 |
Non-Patent Citations (3)
| Title |
|---|
| 张则宇 等: "三种GPS天线相位中心偏差校准方法比较", 《计测技术-计量、测试与校准》 * |
| 裴超: "《北京交通大学硕士学位论文》", 31 December 2014 * |
| 陈涛 等: "GPS天线相位中心变化精确检测试验研究", 《地震地质》 * |
Cited By (20)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109031359A (zh) * | 2018-06-11 | 2018-12-18 | 武汉大学 | 一种基于微波暗室的gnss接收天线绝对相位中心校正方法 |
| CN109061559A (zh) * | 2018-06-21 | 2018-12-21 | 同济大学 | 一种uwb基站天线相位中心偏差建模与改正的研究方法 |
| CN110658537A (zh) * | 2019-11-13 | 2020-01-07 | 深圳市鼎耀科技有限公司 | 微波暗室gnss接收机天线绝对相位中心高精度标定的方法 |
| CN112363186A (zh) * | 2019-12-27 | 2021-02-12 | 广东电网有限责任公司电力科学研究院 | 一种卫星天线相位中心参数的计算方法和装置 |
| CN112363186B (zh) * | 2019-12-27 | 2023-06-02 | 广东电网有限责任公司电力科学研究院 | 一种卫星天线相位中心参数的计算方法和装置 |
| CN113703017B (zh) * | 2020-05-22 | 2023-06-13 | 中国人民解放军战略支援部队信息工程大学 | 一种卫星天线相位中心偏差计算方法及装置 |
| CN113703017A (zh) * | 2020-05-22 | 2021-11-26 | 中国人民解放军战略支援部队信息工程大学 | 一种卫星天线相位中心偏差计算方法及装置 |
| CN113552599A (zh) * | 2021-07-07 | 2021-10-26 | 杭州中科微电子有限公司 | 基于信息融合的gnss接收机天线相位中心校正方法及装置 |
| CN113552599B (zh) * | 2021-07-07 | 2023-11-21 | 杭州中科微电子有限公司 | 基于信息融合的gnss接收机天线相位中心校正方法及装置 |
| CN113568011A (zh) * | 2021-07-16 | 2021-10-29 | 航天科工海鹰集团有限公司 | 一种ppk多径消除技术 |
| CN113267794A (zh) * | 2021-07-20 | 2021-08-17 | 杭州中科微电子有限公司 | 一种基线长度约束的天线相位中心校正方法及装置 |
| CN113721267B (zh) * | 2021-09-01 | 2024-04-12 | 航天科工通信技术研究院有限责任公司 | 一种基于双接收机载波相位差分的gnss欺骗干扰检测方法 |
| CN113721267A (zh) * | 2021-09-01 | 2021-11-30 | 航天科工通信技术研究院有限责任公司 | 一种基于双接收机载波相位差分的gnss欺骗干扰检测方法 |
| CN113899957A (zh) * | 2021-09-30 | 2022-01-07 | 武汉大学 | 基于双差相位观测值的室外天线绝对标定方法及系统 |
| CN113899957B (zh) * | 2021-09-30 | 2024-03-29 | 武汉大学 | 基于双差相位观测值的室外天线绝对标定方法及系统 |
| CN113885054B (zh) * | 2021-12-07 | 2022-03-25 | 广东皓行科技有限公司 | 天线位置确定方法及装置 |
| CN113885054A (zh) * | 2021-12-07 | 2022-01-04 | 广东皓行科技有限公司 | 天线位置确定方法及装置 |
| CN115826008A (zh) * | 2023-02-22 | 2023-03-21 | 广州导远电子科技有限公司 | 天线位移距离范围的检测方法、装置、流动站及存储介质 |
| CN116542132A (zh) * | 2023-04-07 | 2023-08-04 | 武汉大学 | 结合球冠谐模型和神经网络的水汽数据校准与优化方法 |
| CN116542132B (zh) * | 2023-04-07 | 2024-04-16 | 武汉大学 | 结合球冠谐模型和神经网络的水汽数据校准与优化方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107290762B (zh) | 2018-09-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107290762B (zh) | 基于载波三次差分的室外gnss绝对天线相位中心标校方法 | |
| CN104502926B (zh) | 一种基于精密机械臂的室外gnss接收机绝对天线相位中心校正方法 | |
| CN109031359B (zh) | 一种基于微波暗室的gnss接收天线绝对相位中心校正方法 | |
| Ma et al. | The International Celestial Reference Frame as Realized by VeryLong Baseline Interferometry | |
| CN106873009B (zh) | 利用非组合ppp辅助的长基线模糊度解算方法 | |
| CN105182388B (zh) | 一种快速收敛的精密单点定位方法 | |
| CN107579794A (zh) | 一种基于北斗geo导航卫星的精密共视时间频率传递方法 | |
| CN103323888A (zh) | Gnss大气探测数据中对流层延迟误差的消除方法 | |
| CN104483688A (zh) | 基于北斗卫星导航系统的高精度基线解算方法 | |
| CN109613582B (zh) | 一种车载实时单频米级伪距定位方法 | |
| EP3524931A1 (en) | Method for verifying accuracy of rtk tilt compensation measurement | |
| CN102830406B (zh) | 一种gps天线相位中心变化对绝对定位精度的修正方法 | |
| Geng et al. | Multipath mitigation for GPS/Galileo/BDS-3 precise point positioning with overlap-frequency signals | |
| CN105044733B (zh) | 一种高精度的导航卫星tgd参数标定方法 | |
| CN106814376B (zh) | 一种快速精确厘米级单点定位方法 | |
| Stępniak et al. | Field tests of L1 phase centre variation models of surveying-grade GPS antennas | |
| CN103760582A (zh) | 一种遮挡环境下卫星双差观测结构的优化方法 | |
| CN114296113B (zh) | 一种基于姿态数据的天线相位中心测量装置及方法 | |
| CN113465575B (zh) | 一种基于对流层先验信息约束的高落差山区gnss高精度快速网解方法 | |
| CN113899957B (zh) | 基于双差相位观测值的室外天线绝对标定方法及系统 | |
| Wang et al. | Precise orbit determination of LEO satellite using onboard BDS-3 B1C/B2a observations | |
| Liu et al. | Optimization of Loran-C and Beidou integrated positioning performance | |
| CN110764123A (zh) | 一种基于glonass广播星历的伪距定位改进方法 | |
| CN112363186B (zh) | 一种卫星天线相位中心参数的计算方法和装置 | |
| Xu | A concept of precise kinematic positioning and flight-state monitoring from the AGMASCO practice |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |