CN105044733B - 一种高精度的导航卫星tgd参数标定方法 - Google Patents
一种高精度的导航卫星tgd参数标定方法 Download PDFInfo
- Publication number
- CN105044733B CN105044733B CN201510523192.9A CN201510523192A CN105044733B CN 105044733 B CN105044733 B CN 105044733B CN 201510523192 A CN201510523192 A CN 201510523192A CN 105044733 B CN105044733 B CN 105044733B
- Authority
- CN
- China
- Prior art keywords
- tgd
- satellite
- ground station
- parameter
- navigation
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 24
- 239000005433 ionosphere Substances 0.000 claims abstract description 13
- 238000004364 calculation method Methods 0.000 claims description 5
- 230000005540 biological transmission Effects 0.000 claims description 3
- 230000001934 delay Effects 0.000 claims description 3
- 238000000926 separation method Methods 0.000 abstract description 13
- 238000005259 measurement Methods 0.000 abstract description 9
- 230000008859 change Effects 0.000 abstract description 3
- 230000008569 process Effects 0.000 abstract description 2
- 238000012545 processing Methods 0.000 abstract description 2
- 238000007796 conventional method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000006185 dispersion Substances 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
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003380 propellant Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012360 testing method Methods 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/03—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
- G01S19/07—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing data for correcting measured positioning data, e.g. DGPS [differential GPS] or ionosphere corrections
Landscapes
- Engineering & Computer Science (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
本发明公开了一种高精度的导航卫星TGD参数标定方法，直接修正导航卫星频间相位中心的垂向误差，对于水平误差，采用多站加权的方式，充分利用地面站的几何布局进行消除；参数标定过程不需要利用卫星质心位置以及质心与相位中心的位置关系，求解相位中心的位置，而且多站加权处理对卫星星下点坐标、电离层穿刺点坐标的精度要求不高，因此TGD参数标定精度不受卫星质心位置的工厂标定误差、在轨卫星的质心在卫星机械坐标系内的变化、卫星姿态测量误差以及卫星轨道测量误差的影响。将导航卫星频间相位中心分离引起的TGD估计误差降低至原来的5％～15％。根据相位中心实际分离的大小，可将相位中心分离引起的TGD估计误差控制在4cm以下甚至9mm以下。
Description
技术领域
本发明涉及导航卫星在轨标定技术领域，尤其涉及一种高精度的导航卫星TGD参数标定方法。
背景技术
电离层延迟是用户接收机测量值的主要误差源之一。消除电离层延迟误差的最有效办法是组合两个频率的测量值，随之而来的问题是引入了导航卫星的通道间时延偏差(TGD)以及用户接收通道间时延偏差(IFB)。导航卫星TGD参数的精确标定对于用户的定位精度和授时精度具有重要意义。另外，卫星导航***可用于测量电离层的总电子含量(TEC)，具有精度高、成本低的优点。该项技术中最重要的影响因素是导航卫星的TGD参数与地面站接收通道的IFB参数。
用户或地面站接收通道的IFB参数可通过共视导航卫星的方法事先标定，或者通过有线连接方式精确标定接收机的频间时延偏差，且通过无线方式标定用户或地面站天线的频间时延偏差。导航卫星在工厂电测阶段都进行信号发射通道时延的标定，由于发射通道所处的卫星内外环境不断变化，并且发射通道自身出现老化，发射通道时延不是一个固定值。部分GPS卫星的L1与L2频点通道时延差的在轨标定结果与地面测试结果相差1～3ns。部分GPS卫星TGD一年内的波动达到0.2～0.7ns，而十年内的变化达到3ns左右。
TGD参数的标定方法得到了大量研究，研究点主要集中在卫星TGD参数的标定方法，以及解决卫星TGD参数与接收机IFB参数的分离问题，例如论文《Variability of GPSSatellite Differential Group Delay Biases》(David S.Coco等，IEEE Transactionson Aerospace and Electronic System，1991年11月，第27卷第6期)、《Compass***TGD与IFB分离解算方法》(吴晓莉等，第一届中国卫星导航学术年会论文集，2010年)、《基于三频数据的北斗卫星导航***DCB参数精度评估方法》(樊家琛等，中国空间科学技术，2013年8月，第4期)。然而目前的TGD参数标定方法忽视了TGD参数的误差源及其对标定结果的影响。人们容易想到和解决伪码测量误差问题。实际上，由于卫星发射链路的色散效应，导航卫星天线在不同频率下的相位中心存在着一定的几何分离，这会影响TGD的标定精确度。
发明内容
有鉴于此，本发明提供了一种导航卫星TGD参数的精确标定方法，能够消除导航卫星天线频间相位中心分离引入的TGD参数标定误差，提高用户接收机的定位与授时精度和电离层总电子含量(TEC)的测量精度。
一种导航卫星TGD参数的精确标定方法，包括如下步骤：
步骤1、测量待测导航卫星的频率为f1和f2的导航天线的相位中心在垂直方向的距离v，其中，垂直方向为导航卫星指向地心的方向；
步骤2、选定n个地面站，其中n不小于3；各地面站所在的经度线和纬度线围成TGD参数标定观测区，其中，在选定的TGD参数标定时间段内，待测导航卫星的星下点轨迹落于所述观测区内；
步骤3、确定导航信号频率为f1和f2的导航天线在各个地面站接收通道内的传输时延之差，标记为τr1-τr2，即实现地面站接收通道IFB的标定；
步骤4、在TGD参数标定时间段内，各个地面站天线分别持续跟踪待测导航卫星，分别通过地面站接收机获得多个观测时刻下导航信号频率为f1和f2的导航天线的伪码测距值ρ1和ρ2；
步骤5、基于任意一个地面站的任意观测时刻建立观测方程：
其中，c为光速，τ1和τ2分别为导航信号频率为f1和f2的导航天线的发射通道时延，s是以导航卫星对地面站天线仰角为变量的倾斜因子，c1,…,c6是电离层模型系数，φ是电离层穿刺点的地理纬度，λcr是电离层穿刺点的太阳时角；
先将同一个地面站在多个观测时刻对应的(1)式表示的观测方程形成方程组，求解τ1-τ2；
然后求得基于该地面站的TGD参数估计值：其中
步骤6、根据星历计算在TGD参数标定时间段的中间时刻该导航卫星的星下点坐标，根据各个地面站地址以及地面站天线参数，求取各地面站天线的相位中心在地球海平面上投影点的坐标，分别计算所述导航卫星的星下点至各地面站天线投影点的距离；其中，用ri表示所述星下点至第i个地面站对应的投影点距离；
步骤7、根据步骤5获得的各个地面站对应的TGD参数估计值以及步骤6获得的各地面站对应的星下点至投影点距离，求出最终的TGD参数估计结果：其中，表示第i个地面站对应的TGD参数估计值；然后用该最终估计结果T′GD对导航卫星进行TGD参数标定。
较佳的，所述TGD参数标定时间段为1小时至4小时。
较佳的，所述倾斜因子其中Re为地球半径，α为导航卫星对地面站天线的仰角，hI为电离层中电子密度最大处的高度，取350km。
较佳的，所述太阳时角λcr的计算公式为：λcr＝λef+Te·ωe，其中λef是电离层穿刺点的地理经度，Te是观测时刻，ωe为地球自转角速度。
本发明具有如下有益效果：
(1)本发明直接修正导航卫星频间相位中心的垂向误差。对于水平误差，则采用多站加权的方式，充分利用地面站的几何布局进行消除。参数标定过程不需要利用卫星质心位置以及质心与相位中心的位置关系，求解相位中心的位置，而且多站加权处理对卫星星下点坐标、电离层穿刺点坐标的精度要求不高，因此TGD参数标定精度不受卫星质心位置的工厂标定误差、在轨卫星的质心在卫星机械坐标系内的变化(由推进剂消耗等因素引起)、卫星姿态测量误差以及卫星轨道测量误差的影响。
(2)本发明将导航卫星频间相位中心分离引起的TGD估计误差降低至原来的5％～15％。根据相位中心实际分离的大小，可将相位中心分离引起的TGD估计误差控制在4cm以下甚至9mm以下。
(3)本发明新增的参数易于获得，且计算复杂度不大。
附图说明
图1为本发明方法的流程框图；
图2为本发明卫星导航信号的生成、传播与接收链路。
具体实施方式
下面结合附图并举实施例，对本发明进行详细描述。
(1)在导航卫星工厂电测阶段，分别标定导航天线f1和f2频率相位中心的坐标，进一步计算二者在垂直方向的距离v和在水平方向的距离h。其中，垂直方向为在轨卫星指向地心的方向，与卫星机械坐标系的Z轴平行，f1和f2频率相位中心的Z坐标之差即为v，如果在轨卫星f2频率相位中心高于f1频率相位中心(即f1频率相位中心的Z坐标大于f2频率相位中心的Z坐标)，则v为正数，否则v为负数，如果二者高度相同(即Z坐标相等)，则v为零；水平方向与垂直方向正交，导航天线f1和f2频率相位中心在卫星机械坐标系XOY平面内的投影之间的距离即为h，h为正数或零。
(2)选定n个地面站参与该导航卫星的TGD参数标定，要求n不小于3。各地面站的最小纬度和最大纬度所在的纬度线之间、以及最小经度和最大经度所在的经度线之间包围的区域定义为TGD参数标定观测区，要求在TGD参数标定时间段(一般为1～4小时)内，卫星的星下点轨迹始终落于观测区内。
(3)事先标定各个地面站接收通道的IFB，即f1和f2频率导航信号在地面站接收通道(含天线)内的传输时延之差，标记为τr1-τr2。
地面站接收机IFB的标定方法比较成熟，参见《Compass***TGD与IFB分离解算方法》、《基于三频数据的北斗卫星导航***DCB参数精度评估方法》等。
(4)在TGD参数标定时间段内，各个地面站天线分别持续跟踪该导航卫星，分别通过地面站接收机获得多个观测时刻下f1和f2频率伪码测距值ρ1和ρ2。
(5)基于某个地面站某时刻的观测量建立方程，
本发明的观测方程中引入v/c，可消除相位中心在垂直方向的偏差对TGD参数标定精度的影响。c为光速，τ1和τ2分别为导航信号频率为f1和f2的导航天线的发射通道时延，c1,…,c6是电离层模型系数，φ是电离层穿刺点的地理纬度，λcr是电离层穿刺点的太阳时角，λcr＝λef+Te·ωe，其中λef是电离层穿刺点的地理经度，Te是观测时刻，ωe地球自转角速度；s是以卫星对地面站天线仰角为变量的倾斜因子，倾斜因子有不同的模型，可以使用其中Re为地球半径，α为卫星对地面站天线的仰角，hI为电离层中电子密度最大处的高度，可取350km。
电离层穿刺点是卫星与地面天线的连线与高度为hI的球面的交点。已知卫星的位置坐标(可通过提取卫星播发的导航电文中星历参数，或者使用地面控制段生成的星历参数，进而计算求得)、地面站天线相位中心的位置坐标、电离层穿刺点的高度hI，容易求得电离层穿刺点的位置坐标，进一步转换为地理纬度φ和经度λef，而后根据λcr＝λef+Te·ωe计算电离层穿刺点的太阳时角λcr。
同一个地面站多个时刻的观测方程形成方程组，求解τ1-τ2，进一步求得基于该地面站的TGD参数估计值，
(6)根据星历参数计算在TGD参数标定时间段的中间时刻该卫星的星下点坐标(即卫星在地球海平面上投影点的坐标)。根据各个地面站地址以及地面站天线参数，求取各地面站天线的相位中心在地球海平面上投影点的坐标。分别计算该星下点至各地面站天线投影点的距离，ri表示第i个地面站对应的星下点至投影点距离；
(7)根据各个地面站估计的TGD参数，求出最终的TGD参数估计结果：其中，表示第i个地面站对应的TGD参数估计值，后用该最终估计结果T′GD对导航卫星进行TGD参数标定。通过采用多站加权的方式，充分利用地面站的几何布局消除相位中心的水平误差对TGD参数标定精度的影响。
下面通过典型的仿真计算案例说明本发明的效果。
MEO卫星的轨道高度为21528km，地球半径取6378km。f1频率值为1575.42MHz，f2频率值为1268.52MHz。选取3个典型的地面监测站，分别为三亚站(位于北纬18°14′、东经109°31′)、喀什站(位于北纬39°28′、东经75°59′)、长春站(位于北纬43°50′、东经125°16′)。监测站围成的TGD参数标定观测区为北纬19°～43°、东经76°～125°。在TGD参数标定观测区内，分别在经度方向和纬度方向上以1°为步进设定卫星f1频点相位中心C1在地表的投影。另外，在0°～360°范围内，以30°为步进设定相对于自C1向东向量的方位角，其中C2为卫星f2频点相位中心。水平偏差h选取了典型的5mm和5cm，垂向偏差v选取了1cm、10cm、20cm、40cm的典型值。
按照传统方法和本发明的新方法，统计计算卫星相位中心分离引起的TGD参数估计RMS误差如下表所示，其中，考虑卫星天线在空间使用过程中出现的形变，分别给出了修正全部垂向偏差、95％的垂向偏差、90％的垂向偏差的结果。
从上表可以看出对于相位中心分离引起的TGD估计误差，传统方法达到厘米级甚至分米级，而本发明的标定方法可将其控制在厘米级甚至毫米级。按照修正95％的垂向偏差，相位中心分离引起的TGD估计误差可达到4cm以下甚至9mm以下，比传统方法减小了85％～95％。
综上所述，以上仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (4)
1.一种导航卫星TGD参数的标定方法，其特征在于，包括如下步骤：
步骤1、测量待测导航卫星的频率为f1和f2的导航天线的相位中心在垂直方向的距离v，其中，垂直方向为导航卫星指向地心的方向；
步骤2、选定n个地面站，其中n不小于3；各地面站所在的经度线和纬度线围成TGD参数标定观测区，其中，在选定的TGD参数标定时间段内，待测导航卫星的星下点轨迹落于所述观测区内；
步骤3、确定导航信号频率为f1和f2的导航天线在各个地面站接收通道内的传输时延之差，标记为τr1-τr2，即实现地面站接收通道IFB的标定；
步骤4、在TGD参数标定时间段内，各个地面站天线分别持续跟踪待测导航卫星，分别通过地面站接收机获得多个观测时刻下导航信号频率为f1和f2的导航天线的伪码测距值ρ1和ρ2；
步骤5、基于任意一个地面站的任意观测时刻建立观测方程：
其中，c为光速，τ1和τ2分别为导航信号频率为f1和f2的导航天线的发射通道时延，s是以导航卫星对地面站天线仰角为变量的倾斜因子，c1,…,c6是电离层模型系数，φ是电离层穿刺点的地理纬度，λcr是电离层穿刺点的太阳时角；
先将同一个地面站在多个观测时刻对应的(1)式表示的观测方程形成方程组，求解τ1-τ2；
然后求得基于该地面站的TGD参数估计值：其中
步骤6、根据星历计算在TGD参数标定时间段的中间时刻该导航卫星的星下点坐标，根据各个地面站地址以及地面站天线参数，求取各地面站天线的相 位中心在地球海平面上投影点的坐标，分别计算所述导航卫星的星下点至各地面站天线投影点的距离；其中，用ri表示所述星下点至第i个地面站对应的投影点距离；
步骤7、根据步骤5获得的各个地面站对应的TGD参数估计值以及步骤6获得的各地面站对应的星下点至投影点距离，求出最终的TGD参数估计结果： 其中，表示第i个地面站对应的TGD参数估计值；然后用该最终估计结果T ′GD 对导航卫星进行TGD参数标定。
2.如权利要求1所述的一种导航卫星TGD参数的标定方法，其特征在于，所述TGD参数标定时间段为1小时至4小时。
3.如权利要求1所述的一种导航卫星TGD参数的标定方法，其特征在于，所述倾斜因子其中Re为地球半径，α为导航卫星对地面站天线的仰角，hI为电离层中电子密度最大处的高度，取350km。
4.如权利要求1所述的一种导航卫星TGD参数的标定方法，其特征在于，所述太阳时角λcr的计算公式为：λcr＝λef+Te·ωe，其中λef是电离层穿刺点的地理经度，Te是观测时刻，ωe为地球自转角速度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510523192.9A CN105044733B (zh) | 2015-08-24 | 2015-08-24 | 一种高精度的导航卫星tgd参数标定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510523192.9A CN105044733B (zh) | 2015-08-24 | 2015-08-24 | 一种高精度的导航卫星tgd参数标定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105044733A CN105044733A (zh) | 2015-11-11 |
| CN105044733B true CN105044733B (zh) | 2017-08-11 |
Family
ID=54451408
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510523192.9A Active CN105044733B (zh) | 2015-08-24 | 2015-08-24 | 一种高精度的导航卫星tgd参数标定方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105044733B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108710144B (zh) * | 2018-04-02 | 2021-02-02 | 武汉大学 | 零/短基线glonass相位通道间偏差校准方法及*** |
| CN110554414B (zh) * | 2019-08-25 | 2022-11-15 | 中国科学院国家授时中心 | 基于气象参数及信号质量的***误差预测方法 |
| US11650327B2 (en) * | 2020-11-20 | 2023-05-16 | Qualcomm Incorporated | Antenna phase center compensation for orbital assistance data |
| CN116208230B (zh) * | 2023-01-19 | 2024-02-13 | 长光卫星技术股份有限公司 | 一种卫星自主数传快速判定及任务参数计算方法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102565814B (zh) * | 2011-12-15 | 2014-10-22 | 北京航空航天大学 | 卫星导航***的信号精度和定位服务可用性的评估方法 |
| CN103777210B (zh) * | 2012-10-22 | 2017-07-11 | 中兴通讯股份有限公司 | 电离层延迟修正参数的使用方法、装置及终端 |
| CN104833993B (zh) * | 2015-05-11 | 2017-07-28 | 中国科学院国家授时中心 | 一种基于卫星和接收机频间偏差之和的北斗定位方法 |
-
2015
- 2015-08-24 CN CN201510523192.9A patent/CN105044733B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105044733A (zh) | 2015-11-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Rama Rao et al. | On the validity of the ionospheric pierce point (IPP) altitude of 350 km in the Indian equatorial and low-latitude sector | |
| Lanyi et al. | A comparison of mapped and measured total ionospheric electron content using global positioning system and beacon satellite observations | |
| CN107290762B (zh) | 基于载波三次差分的室外gnss绝对天线相位中心标校方法 | |
| Yasyukevich et al. | Influence of GPS/GLONASS differential code biases on the determination accuracy of the absolute total electron content in the ionosphere | |
| CN109613585A (zh) | 一种对基站天线超短基线gnss双天线实时测向的方法 | |
| CN104502926A (zh) | 一种基于精密机械臂的室外gnss接收机绝对天线相位中心校正方法 | |
| CN105044733B (zh) | 一种高精度的导航卫星tgd参数标定方法 | |
| CN114488228B (zh) | 一种适用于动态载体平台的gnss多路径误差削弱方法 | |
| CN107607969A (zh) | 一种基于dcb改正的四***伪距定位方法 | |
| JP2010112842A (ja) | Ｔｅｃマップ及び受信機バイアスの作成及び測定方法及び装置 | |
| CN107402394A (zh) | 一种星载测频定位误差源在轨标校方法和装置 | |
| CN107991696B (zh) | 针对终端载体为高轨卫星的卫星导航信号仿真方法 | |
| US8519886B2 (en) | System and method for using GPS interferometry to determine atmospheric conditions | |
| JP5077054B2 (ja) | 移動体用測位システム | |
| CN102023301B (zh) | 一种中轨道卫星搜救***的选星方法 | |
| CN113176596A (zh) | 气压高程约束定位方法 | |
| CN106814376B (zh) | 一种快速精确厘米级单点定位方法 | |
| CN116858290A (zh) | 基于大型无人机的深远海海面高度观测定标方法和*** | |
| Banville et al. | Antenna rotation and its effects on kinematic precise point positioning | |
| CN109471102B (zh) | 一种惯组误差修正方法 | |
| CN109490828A (zh) | 基于同源基线阵列的定位方法 | |
| Xia et al. | Estimation and evaluation of COSMIC radio occultation excess phase using undifferenced measurements | |
| Dong-ju et al. | The application of GIM in precise orbit determination for LEO satellites with single-frequency GPS measurements | |
| Karaush et al. | Calibration of GNSS receiver with GLONASS signals | |
| Zhou et al. | Tropospheric delay correction of VLBI stations for the real-time trajectory determination of the Chang'E-5 spacecraft |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |