CN103439704B - 一种基于虚拟站的三程观测数据处理方法 - Google Patents
一种基于虚拟站的三程观测数据处理方法 Download PDFInfo
- Publication number
- CN103439704B CN103439704B CN201310314991.6A CN201310314991A CN103439704B CN 103439704 B CN103439704 B CN 103439704B CN 201310314991 A CN201310314991 A CN 201310314991A CN 103439704 B CN103439704 B CN 103439704B
- Authority
- CN
- China
- Prior art keywords
- station
- observed quantity
- cell site
- satellite
- phantom
- 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
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种基于虚拟站的三程观测数据处理方法,属于航天测量与控制领域,该方法首先设置虚拟站M,所述虚拟站M位于发射站S与接收站R的中点;计算基于虚拟站的理论测距观测量ρM;使用ρM代替基于发射站S与接收站R的理论测距观测量ρC;计算基于虚拟站的理论测速观测量使用代替基于发射站S与接收站R的理论测速观测量本发明的有益效果是:基于虚拟站的理论观测量计算简单快捷,满足计算精度要求,可以用来替代计算复杂的基于发射站S与接收站R的理论观测量,在航天测控领域可以进行推广。
Description
技术领域
本发明专利属于航天测量与控制领域,涉及建立基于虚拟站的测量模型代替双站收发USB测量模型的数据处理方法。
背景技术
多年来我国的航天测控网普遍采用上下行统一S波段载波对航天器进行跟踪测量,其中上下行多采用同一个天线收发。近年来随着任务类型的多样化、国际合作和设备升级改造等带来的现实问题,导致部分任务中采用不同的USB上下行天线(称三程测距测速)。不同的天线存在站址坐标不同、设备状况、积分时间和测量精度不一致等问题,从而使得在定轨数据处理中不可能套用已有的测量模型和软件。例如,在某航天任务中,某站对某卫星的测量设备为USB,其发射天线和接收天线站址不同;测速的积分间隔为2秒,远大于一般USB采用的0.4秒。如果将该站测量模型等同于同一站址的双程USB测量模型,将导致测距最大有250米的误差;测速最大有1m/s的误差,最终给定轨带来较大的误差。
目前,国内在针对不同站址上下行USB测量数据处理方面的研究较少,需要针对实际问题建立一套测量模型,并且在工程中加以验证。
发明内容
本发明的目的是:为了解决某些情况下卫星测控时USB站发射天线和接收天线不在一起,导致的无法使用传统USB双程测距测速模型进行轨道确定的问题,给出了一套新颖的、效率较高的基于虚拟站的三程观测模型数据处理方法。
本发明的技术方案是:基于虚拟站的三程观测模型数据处理方法,包括如下步骤:
步骤一:参阅图1,设置虚拟站M,所述虚拟站M位于发射站S与接收站R的中点;发射站S到卫星V的距离为ρS,实际观测量为ρS′,卫星V到接收站R的距离为ρR,实际观测量为ρR′;发射站S到卫星V之间连线与卫星V到接收站R之间连线的夹角为θ;基于发射站S与接收站R的实际测距观测量为基于发射站S与接收站R的理论测距观测量为实际测速观测量为其中ρOb和ρOe分别为计算差分测速的开始时刻和结束时刻对应的实际测距观测量,T为计算差分测速的开始时刻到结束时刻的时间间隔;理论测速观测量为其中ρCb和ρCe分别为计算差分测速的开始时刻和结束时刻对应的理论测距观测量,T为计算差分测速的开始时刻到结束时刻的时间间隔;
步骤二:计算基于虚拟站的理论测距观测量其中,为虚拟站M发送信号时刻对应的位置矢量,为卫星V收到信号时刻对应的位置矢量;为虚拟站M收到信号时刻对应的位置矢量,为卫星V转发信号时刻对应的位置矢量;
步骤三:使用ρM代替基于发射站S与接收站R的理论测距观测量ρC;
根据几何关系可以得到 由于SR<<ρS,ρR,即sinθ≈θ,代入上式得该量为小量,所以ρM与ρC之间误差可以忽略;
步骤四:计算基于虚拟站的理论测速观测量其中ρMe和ρMb分别为计算差分测速的开始时刻和结束时刻对应的基于虚拟站的理论测距观测量,T为开始时刻到结束时刻的时间间隔。
步骤五:使用代替基于发射站S与接收站R的理论测速观测量
计算使用代替的误差利用步骤三中的误差结果,容易得到其中θb和θe分别为计算差分测速的开始时刻和结束时刻对应的SV与VR的夹角,由于SR<<ρS,ρR,故为小量,该误差可以忽略。
至此,得到了基于虚拟站的理论测距观测量ρM和理论测速观测量本数据处理方法完毕。
本发明的有益效果是:基于虚拟站的理论观测量计算简单快捷,满足计算精度要求,可以用来替代计算复杂的基于发射站S与接收站R的理论测速观测量,在航天测控领域可以进行推广。
附图说明
图1是卫星测站几何关系示意图;
图2是卫星第8圈某站测距修正误差示意图;
图3是卫星第9圈某站测距修正误差示意图;
图4是卫星第13圈某站测距修正误差示意图;
具体实施方式
本实施例中,采用实测的嫦娥二号(CE-2)USB测距测速数据,采用宇航动力学国家重点实验室自主研发的多星并行高精度轨道计算软件系统(PASAX)进行,具体实施方式如下:
步骤一:参阅图1,以CE-2工程中某上行测控站作为发射站R,某下行测控站作为接收站S,将RS连线的中点设置为虚拟站M,计算M站坐标。实际观测量为2011年5月1日10时至当日12时,R站上行发射,经CE-2转发,S站下行接收的USB测距测速数据,采样率为5s;准备好该时间段内CE-2的精密轨道作为下一步计算的基础。本例中,两站的距离RS约500m,ρS,ρR约38万km。
步骤二:利用虚拟站M的坐标和CE-2精密轨道,计算各对应时刻M站和CE-2的位置矢量,进而计算各时刻基于虚拟站的理论测距观测量
步骤三:利用发射站S、接收站R的坐标和CE-2精密轨道,计算各对应时刻S、R站的位置矢量和CE-2的位置矢量,进而计算各时刻基于发射站S与接收站R的理论测距观测量计算各个时刻两种理论测距观测量的差别结果显示两者差别的最大值小于0.1m;可以使用ρM代替基于发射站S与接收站R的理论测距观测量ρC。
步骤四:利用虚拟站M的坐标和CE-2精密轨道,按照平均测速公式计算各对应时刻基于虚拟站的理论测速观测量
步骤五:利用发射站S、接收站R的坐标和CE-2精密轨道,按照平均测速公式计算各对应时刻基于发射站S与接收站R的理论测速观测量计算各个时刻两种理论测速观测量的差别结果显示两者差别的最大值小于0.05mm/s,;故可以使用代替基于发射站S与接收站R的理论测速观测量
根据实际数据的验证,本发明的观测模型数据处理方法准确可信,满足工程计算要求。
在该实施例中,图2是卫星第8圈某站测距修正误差示意图;图3是卫星第9圈某站测距修正误差示意图;图4是卫星第13圈某站测距修正误差示意图。图2-4表明使用本数据处理方法,这3圈测站测距修正误差的最大值约0.07m,小于所规定的误差限。故本发明的观测模型数据处理方法准确可信,满足计算精度要求。
Claims (1)
1.基于虚拟站的三程观测模型数据处理方法,其特征在于:包括如下步骤:
步骤一:设置虚拟站M,所述虚拟站M位于发射站S与接收站R的中点;发射站S到卫星V的距离为ρS,实际观测量为ρ′S,卫星V到接收站R的距离为ρR,实际观测量为ρ′R;发射站S到卫星V之间连线与卫星V到接收站R之间连线的夹角为θ;基于发射站S与接收站R的实际测距观测量为基于发射站S与接收站R的理论测距观测量为实际测速观测量为其中ρOb和ρOe分别为计算差分测速的开始时刻和结束时刻对应的实际测距观测量,T为计算差分测速的开始时刻到结束时刻的时间间隔;理论测速观测量为其中ρCb和ρCe分别为计算差分测速的开始时刻和结束时刻对应的理论测距观测量,T为计算差分测速的开始时刻到结束时刻的时间间隔;
步骤二:计算基于虚拟站的理论测距观测量ρM:其中, 为虚拟站M发送信号时刻对应的位置矢量,为卫星V收到信号时刻对应的位置矢量; 为虚拟站M收到信号时刻对应的位置矢量,为卫星V转发信号时刻对应的位置矢量;
步骤三:使用ρM代替基于发射站S与接收站R的理论测距观测量ρC;
步骤四:计算基于虚拟站的理论测速观测量 其中ρMb和ρMe分别为计算差分测速的开始时刻和结束时刻对应的基于虚拟站的理论测距观测量,T为开始时刻到结束时刻的时间间隔;
步骤五:使用代替基于发射站S与接收站R的理论测速观测量
至此,得到了基于虚拟站的理论测距观测量ρM和理论测速观测量
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310314991.6A CN103439704B (zh) | 2013-07-24 | 2013-07-24 | 一种基于虚拟站的三程观测数据处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310314991.6A CN103439704B (zh) | 2013-07-24 | 2013-07-24 | 一种基于虚拟站的三程观测数据处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103439704A CN103439704A (zh) | 2013-12-11 |
CN103439704B true CN103439704B (zh) | 2015-05-13 |
Family
ID=49693409
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310314991.6A Active CN103439704B (zh) | 2013-07-24 | 2013-07-24 | 一种基于虚拟站的三程观测数据处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103439704B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101881820A (zh) * | 2010-07-09 | 2010-11-10 | 中国科学院测量与地球物理研究所 | 用于地球静止轨道卫星精密定轨的镜面投影方法 |
CN102243314A (zh) * | 2011-04-14 | 2011-11-16 | 山东大学 | 一种规律移动的船载基准站实时差分测量装置及其工作方法 |
EP2397868A1 (en) * | 2010-06-15 | 2011-12-21 | The European Union, represented by the European Commission | Method of providing an authenticable time-and-location indication |
CN103064097A (zh) * | 2011-02-14 | 2013-04-24 | 天宝导航有限公司 | 使用用于合成参考数据的电离层模型的gnss信号处理 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6992617B2 (en) * | 2003-11-13 | 2006-01-31 | Global Locate, Inc. | Method and apparatus for monitoring the integrity of satellite tracking data used by a remote receiver |
-
2013
- 2013-07-24 CN CN201310314991.6A patent/CN103439704B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2397868A1 (en) * | 2010-06-15 | 2011-12-21 | The European Union, represented by the European Commission | Method of providing an authenticable time-and-location indication |
CN101881820A (zh) * | 2010-07-09 | 2010-11-10 | 中国科学院测量与地球物理研究所 | 用于地球静止轨道卫星精密定轨的镜面投影方法 |
CN103064097A (zh) * | 2011-02-14 | 2013-04-24 | 天宝导航有限公司 | 使用用于合成参考数据的电离层模型的gnss信号处理 |
CN102243314A (zh) * | 2011-04-14 | 2011-11-16 | 山东大学 | 一种规律移动的船载基准站实时差分测量装置及其工作方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103439704A (zh) | 2013-12-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103675861B (zh) | 一种基于星载gnss多天线的卫星自主定轨方法 | |
CN101403790B (zh) | 单频gps接收机的精密单点定位方法 | |
CN103471544B (zh) | 一种基于北斗的高精度位移形变监测应用系统 | |
CN101738598B (zh) | 对卫星或飞行器轨道和位置实时精密测量方法 | |
CN101156080B (zh) | 一种基于gps导航卫星系统发射信号主要用于地形运动、大型设施以及民用建筑工程监控和测量的系统和方法 | |
CN102998690B (zh) | 一种基于gps载波双差方程的姿态角直接求解方法 | |
CN103674034B (zh) | 多波束测速测距修正的鲁棒导航方法 | |
CN101446634A (zh) | 一种高精度位置、方位角和俯仰角的组合测量方法及装置 | |
CN102981144A (zh) | 空中运动平台对目标的三维无源定位方法 | |
CN102608642A (zh) | 北斗/惯性组合导航系统 | |
CN102565834A (zh) | 一种单频gps测向系统及其测向定位方法 | |
CN105425261A (zh) | 基于GPS/Beidou2/INS的组合导航与定位方法 | |
CN107462905A (zh) | 一种北斗gps双模定位方法、定位终端及定位系统 | |
CN102591343A (zh) | 基于两行根数的卫星轨道维持控制方法 | |
CN102323597A (zh) | 编队绕飞卫星基于gps的星间基线测量方法 | |
CN101833103A (zh) | 一种基于伪卫星的深凹露天矿定位方法 | |
CN201266237Y (zh) | 高精度gps差分系统 | |
CN106773783B (zh) | 一种连续运行参考站系统的半物理仿真测试系统及方法 | |
CN103344977A (zh) | 一种精度可切换的定位跟踪装置及定位跟踪方法 | |
CN103968836A (zh) | 一种基于时序伪距差分推算移动目标位置的方法及装置 | |
CN103235321A (zh) | Gps单站伪距定位精密授时方法 | |
CN107782317B (zh) | 一种倒vlbi深空单向导航定位方法 | |
CN107728125A (zh) | 一种电离层闪烁的卫星信号模拟系统及使用方法 | |
CN106840113A (zh) | 一种基于星基差分增强技术的深远海波浪和潮位测量方法 | |
CN104913743A (zh) | 基于惯性测量的电力铁塔变形监测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |