CN111308432B - 一种利用测速数据评估航天器测距数据精度的方法 - Google Patents
一种利用测速数据评估航天器测距数据精度的方法 Download PDFInfo
- Publication number
- CN111308432B CN111308432B CN201911217022.2A CN201911217022A CN111308432B CN 111308432 B CN111308432 B CN 111308432B CN 201911217022 A CN201911217022 A CN 201911217022A CN 111308432 B CN111308432 B CN 111308432B
- Authority
- CN
- China
- Prior art keywords
- result
- difference
- ranging
- distance measurement
- measurement result
- 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
Images
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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
- G01S13/08—Systems for measuring distance only
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Navigation (AREA)
Abstract
本发明涉及一种利用测速数据评估航天器测距数据精度的方法,包括:S1.获取地面站对航天器的测距结果和测速结果;S2.以相邻时刻的所述测距结果作差获取差分测距结果;S3.基于所述差分测距结果和所述测速结果作差获取双差分结果;S4.基于所述双差分结果计算所述测距结果的精度。本方案适用于在具有航天器高精度测速数据的情况下,对获取的航天器测距数据进行客观、准确的精度评估,用于任务中分析航天器和地面测控设备的测距性能,并确定测距数据用于定轨的权重等。
Description
技术领域
本发明涉及航天器测控领域,尤其涉及一种利用测速数据评估航天器测距数据精度的方法。
背景技术
在航天测控任务中,地面测控设备对航天器的测距数据是航天器导航的重要观测量。对测距数据开展精度评估是分析航天器和地面测控设备的测距性能,并确定测距数据用于定轨的权重的基本前提。对于测距数据的精度评估,可以根据测控链路电平情况利用测距随机误差的理论计算公式予以估算;在了解航天器和地面测控设备性能指标和工作状态的前提下,也可根据经验先进行粗略估计(用于定轨),在定轨后基于测距数据的拟合后残差进一步评估其精度。从原理上讲,理论计算无法顾及未考虑到的模型以外的误差,定轨后的拟合后残差通常会受到定轨误差的附加影响。因此,以上方法均无法客观、准确地反映测距数据的精度。
在目前航天测控的统一载波测控体制下,地面测控设备在获得航天器测距数据的同时通常也获取了高精度的测速数据。从测量模型来看,测速数据是航天器对测站径向距离的时间变化率,测距数据是航天器对测站的径向距离,测速数据是测距数据的时间差分;从信号模型来看,测速数据来自于载波锁相环,测距数据来自于测距通道的锁相环,二者受到的噪声相互独立;此外,在航天测控任务中,测速数据的精度(mm/s级)通常远高于测距数据的精度(m级)。
由于测速数据中包含了测距信息、噪声独立于测距数据、精度远高于测距数据,可望利用测速数据更为客观、准确地评估测距数据的精度。
发明内容
本发明的目的在于提供一种利用测速数据评估航天器测距数据精度的方法,以便评估航天器和地面测控设备的测距性能,确定测距数据用于定轨的权重。
为实现上述发明目的,本发明提供一种利用测速数据评估航天器测距数据精度的方法,包括:
S1.获取地面站对航天器的测距结果和测速结果;
S2.以相邻时刻的所述测距结果作差获取差分测距结果;
S3.基于所述差分测距结果和所述测速结果作差获取双差分结果;
S4.基于所述双差分结果计算所述测距结果的精度。
根据本发明的一个方面,步骤S1中,以预设间隔的第一时刻获取所述测距结果,以具有与所述预设间隔相同时间间隔的第二时刻获取所述测速结果;
所述预设间隔为1s;
所述第二时刻是两个相邻所述第一时刻的中间点时刻。
根据本发明的一个方面,步骤S2中,以与所述第二时刻相同的时刻获取所述差分测距结果。
根据本发明的一个方面,所述差分测距结果为:
dρi=ρi+1-ρi
(i=1,2,...,n-1)
其中,dρ1表示差分测距结果,ρi和ρi+1表示相邻时刻的测距结果。
根据本发明的一个方面,步骤S3中,将对应时刻的所述差分测距结果和所述测速结果作差获取所述双差分结果,且所述双差分结果为:
δi=dρi-vi
(i=1,2,...,n-1)
其中,δi表示双差分结果,vi表示测速结果。
根据本发明的一个方面,所述步骤S4中所述测距结果的精度表示为:
其中,σ表示测距结果的精度。
根据本发明的一种方案,适用于在具有航天器高精度测速数据的情况下,对获取的航天器测距数据进行客观、准确的精度评估,用于任务中分析航天器和地面测控设备的测距性能,并确定测距数据用于定轨的权重等。
根据本发明的一种方案,本发明基于测速数据中包含了测距信息、噪声独立于测距数据、精度远高于测距数据的特点,首先进行相邻两点测距数据差分,然后进行差分测距结果与测速数据差分,通过两次差分后的残差,评估测距数据的精度。该方法充分利用了测距、测速数据的特点,模型、理论清晰,计算方便,且不引入其他无关因素干扰,从而实现测距数据客观、准确的精度评估。
附图说明
图1示意性表示根据本发明的一种实施方式的利用测速数据评估航天器测距数据精度的方法的步骤框图。
具体实施方式
为了更清楚地说明本发明实施方式或现有技术中的技术方案,下面将对实施方式中所需要使用的附图作简单地介绍。显而易见地,下面描述中的附图仅仅是本发明的一些实施方式,对于本领域普通技术人员而言,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
在针对本发明的实施方式进行描述时,术语“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”所表达的方位或位置关系是基于相关附图所示的方位或位置关系,其仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此上述术语不能理解为对本发明的限制。
下面结合附图和具体实施方式对本发明作详细地描述,实施方式不能在此一一赘述,但本发明的实施方式并不因此限定于以下实施方式。
如图1所示,根据本发明的一种实施方式,本发明的一种利用测速数据评估航天器测距数据精度的方法,包括:
S1.获取地面站对航天器的测距结果和测速结果;
S2.以相邻时刻的测距结果作差获取差分测距结果;
S3.基于差分测距结果和测速结果作差获取双差分结果;
S4.基于双差分结果计算测距结果的精度。
根据本发明的一种实施方式,步骤S1中,以预设间隔的第一时刻获取测距结果,以具有与预设间隔相同时间间隔的第二时刻获取测速结果。在本实施方式中,预设间隔为1s,第二时刻是两个相邻第一时刻的中间点时刻。在本实施方式中,设地面测控设备对航天器在t=1,2,…,n(s)的第一时刻的测距结果为ρ1,ρ2,…,ρn,在t=1.5,2.5,…,n-0.5(s)的第二时刻的测速结果为v1,v2,…,vn-1。
根据本发明的一种实施方式,步骤S2中,以与第二时刻相同的时刻获取差分测距结果。在本实施方式中,相邻两点测距数据作差,得到t=1.5,2.5,…,n-0.5(s)时刻的差分测距结果dρ1,dρ2,…,dρn-1。在本实施方式中,差分测距结果表示为:
dρi=ρi+1-ρi
(i=1,2,...,n-1)
其中,dρ1表示差分测距结果,ρi和ρi+1表示相邻时刻的测距结果。
在本实施方式中,通过上述步骤利用相邻时刻的测距数据差分,近似地求解中间时刻(即第二时刻)的速度结果。
根据本发明的一种实施方式,步骤S3中,将对应时刻的差分测距结果和测速结果作差获取双差分结果,且双差分结果为:
δi=dρi-vi
(i=1,2,...,n-1)
其中,δi表示双差分结果,vi表示测速结果。
在本实施方式中,通过上述步骤对比实际测速结果与利用测距数据求得的速度结果的差异(该差异的主要来源是测距数据的随机误差)。
根据本发明的一种实施方式,步骤S4中测距结果的精度表示为:
其中,σ表示测距结果的精度。
在本实施方式中,通过上述步骤求解上述差异的RMS值,并归一化为测距数据的精度,即利用测速数据评估得到的测距数据标准差。
本发明基于测速数据中包含了测距信息、噪声独立于测距数据、精度远高于测距数据的特点,首先进行相邻两点测距数据差分,然后进行差分测距结果与测速数据差分,通过两次差分后的残差,评估测距数据的精度。该方法充分利用了测距、测速数据的特点,模型、理论清晰,计算方便,且不引入其他无关因素干扰,从而实现测距数据客观、准确的精度评估。
以下模拟产生一组测距、测速数据,对本发明的具体实施方式作进一步说明。
设航天器的测距数据对应时刻t=1~10s(每秒一点,共10点),测速数据对应时刻t=1.5~9.5s(每秒一点,共9点),具体数据如下(这里为举例说明实施方式,所取时长较短,实际应用中所取时长应更具统计意义:如5-10分钟,共300-600点):
[ρ1 ρ2 ρ3 ρ4 ρ5 ρ6 ρ7 ρ8 ρ9 ρ10]=[1.538 3.834 0.741 4.862 5.319 4.6926.566 8.343 12.578 12.769];
[v1 v2 v3 v4 v5 v6 v7 v8 v9 v10]=[1.0007 0.9988 1.0007 1.0016 1.0005 1,0010 1.0007 0.9997 1.0003];
基于以上数据,按照如下流程评估航天器测距数据精度。
在本实施方式中,将相邻两点测距数据作差,得到t=1.5,2.5,…,9.5(s)时刻的差分测距结果dρ1,dρ2,…,dρ9:
[dρ1 dρ2 dρ3 dρ4 dρ5 dρ6 dρ7 dρ8 dρ9]
=[ρ2-ρ1 ρ3-ρ2 ρ4-ρ3 ρ5-ρ4 ρ6-ρ5 ρ7-ρ6 ρ8-ρ7 ρ9-ρ8 ρ10-ρ9]
=[1.0007 0.9988 1.0007 1.0016 1.0005 1.0010 1.0007 0.9997 1.0003]
在本实施方式中,将对应时刻的差分测距结果、测速数据进一步作差,得到:
[δ1 δ2 δ3 δ4 δ5 δ6 δ7 δ8 δ9]
=[dρ1-v1 dρ2-v2 dρ3-v3 dρ4-v4 dρ5-v5 dρ6-v6 dρ7-v7 dρ8-v8 dρ9-v9]
=[1.2953 -4.0918 3.1203 -0.5446 -1.6275 0.8730 0.7763 3.2353 -0.8093]
在本实施方式中,计算测距数据精度σ:
其中,所得的值即为利用测速数据评估得到的测距数据标准差。
上述内容仅为本发明的具体方案的例子,对于其中未详尽描述的设备和结构,应当理解为采取本领域已有的通用设备及通用方法来予以实施。
以上所述仅为本发明的一个方案而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
2.根据权利要求1所述的方法,其特征在于,步骤S1中,以预设间隔的第一时刻获取所述测距结果,以具有与所述预设间隔相同时间间隔的第二时刻获取所述测速结果;
所述预设间隔为1s;
所述第二时刻是两个相邻所述第一时刻的中间点时刻。
3.根据权利要求2所述的方法,其特征在于,步骤S2中,以与所述第二时刻相同的时刻获取所述差分测距结果。
4.根据权利要求3所述的方法,其特征在于,所述差分测距结果为:
dρi=ρi+1-ρi
(i=1,2,...,n-1)
其中,dρ1表示差分测距结果,ρi和ρi+1表示相邻时刻的测距结果。
5.根据权利要求4所述的方法,其特征在于,步骤S3中,将对应时刻的所述差分测距结果和所述测速结果作差获取所述双差分结果,且所述双差分结果为:
δi=dρi-vi
(i=1,2,...,n-1)
其中,δi表示双差分结果,vi表示测速结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911217022.2A CN111308432B (zh) | 2019-12-03 | 2019-12-03 | 一种利用测速数据评估航天器测距数据精度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911217022.2A CN111308432B (zh) | 2019-12-03 | 2019-12-03 | 一种利用测速数据评估航天器测距数据精度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111308432A CN111308432A (zh) | 2020-06-19 |
CN111308432B true CN111308432B (zh) | 2022-03-22 |
Family
ID=71157907
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911217022.2A Active CN111308432B (zh) | 2019-12-03 | 2019-12-03 | 一种利用测速数据评估航天器测距数据精度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111308432B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103529482A (zh) * | 2013-10-25 | 2014-01-22 | 中国人民解放军国防科学技术大学 | 一种高精度确定载体动态加速度的方法 |
CN103674034A (zh) * | 2013-12-26 | 2014-03-26 | 北京控制工程研究所 | 多波束测速测距修正的鲁棒导航方法 |
CN105043342A (zh) * | 2015-07-03 | 2015-11-11 | 中国电建集团西北勘测设计研究院有限公司 | 单向精密测距三角高程测量方法 |
CN105974415A (zh) * | 2016-06-24 | 2016-09-28 | 西安电子科技大学 | 一种机载sar方位空变运动误差的高精度补偿方法 |
CN106767824A (zh) * | 2016-12-14 | 2017-05-31 | 中国人民解放军63921部队 | 一种计算双探测器在地外天体表面相对位置的方法 |
CN107782317A (zh) * | 2017-09-20 | 2018-03-09 | 西安空间无线电技术研究所 | 一种倒vlbi深空单向导航定位方法 |
CN108828642A (zh) * | 2018-08-01 | 2018-11-16 | 太原理工大学 | 一种ins辅助bds单频接收机的模糊度快速解算方法 |
CN108957433A (zh) * | 2018-08-14 | 2018-12-07 | 中国电子科技集团公司第三十八研究所 | 一种双基地雷达高精度测距装置及方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE19610970A1 (de) * | 1996-03-20 | 1997-09-25 | Reinhard Dipl Ing Stolle | Radarverfahren und -vorrichtung zur Messung von Entfernungen und Geschwindigkeiten |
ATE299273T1 (de) * | 2002-01-23 | 2005-07-15 | Micro Optronic Messtechnik Gmb | Verfahren und vorrichtung zur optischen entfernungsmessung |
CN102103210A (zh) * | 2009-12-17 | 2011-06-22 | 中国石油大学(北京) | 一种卫星导航系统性能评估系统 |
US8344794B2 (en) * | 2010-06-22 | 2013-01-01 | O2Micro Inc. | Signal monitoring systems |
US9288625B2 (en) * | 2014-06-30 | 2016-03-15 | Qualcomm Technologies International, Ltd. | Method for determining location of wireless devices based on information within messages received from other network devices |
CN104851286B (zh) * | 2014-12-18 | 2017-11-07 | 北京市交通运行监测调度中心 | 一种利用公交车辆gps数据动态评估道路交通状态的方法 |
-
2019
- 2019-12-03 CN CN201911217022.2A patent/CN111308432B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103529482A (zh) * | 2013-10-25 | 2014-01-22 | 中国人民解放军国防科学技术大学 | 一种高精度确定载体动态加速度的方法 |
CN103674034A (zh) * | 2013-12-26 | 2014-03-26 | 北京控制工程研究所 | 多波束测速测距修正的鲁棒导航方法 |
CN105043342A (zh) * | 2015-07-03 | 2015-11-11 | 中国电建集团西北勘测设计研究院有限公司 | 单向精密测距三角高程测量方法 |
CN105974415A (zh) * | 2016-06-24 | 2016-09-28 | 西安电子科技大学 | 一种机载sar方位空变运动误差的高精度补偿方法 |
CN106767824A (zh) * | 2016-12-14 | 2017-05-31 | 中国人民解放军63921部队 | 一种计算双探测器在地外天体表面相对位置的方法 |
CN107782317A (zh) * | 2017-09-20 | 2018-03-09 | 西安空间无线电技术研究所 | 一种倒vlbi深空单向导航定位方法 |
CN108828642A (zh) * | 2018-08-01 | 2018-11-16 | 太原理工大学 | 一种ins辅助bds单频接收机的模糊度快速解算方法 |
CN108957433A (zh) * | 2018-08-14 | 2018-12-07 | 中国电子科技集团公司第三十八研究所 | 一种双基地雷达高精度测距装置及方法 |
Non-Patent Citations (2)
Title |
---|
"Designing optical bandwidth synthesis frequency sequences for VLBI experiments conducted jointly by Chinese deep space stations and VLBI stations";Dezhen Xu等;《2015 Fifth international conference on instrumentation and measurement, computer, communication and control》;20151231;第1374-1378页 * |
"差分VLBI测量的一种实现方案";李金岭等;《测绘科学技术学报》;20070430;第24卷(第2期);第79-82页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111308432A (zh) | 2020-06-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9915550B2 (en) | Method and apparatus for data fusion of a three-axis magnetometer and three axis accelerometer | |
CN109459778A (zh) | 基于抗差方差分量估计的码伪距/多普勒联合测速方法及其应用 | |
EP2480908B1 (en) | Estimating user device location in a wireless network | |
CN106772478A (zh) | 一种基于历元‑星间差分约束的定位方法 | |
US20200103477A1 (en) | Method and system for magnetic fingerprinting | |
CN109490855A (zh) | 一种车载雷达标定方法、装置及车辆 | |
CN107870001A (zh) | 一种基于椭球拟合的磁力计校正方法 | |
US8812225B2 (en) | Electronic navigation device for a human and related methods | |
CN104359496B (zh) | 基于垂线偏差补偿的高精度姿态修正方法 | |
CN110657806B (zh) | 一种基于CKF、chan解算和Savitzky-Golay平滑滤波的位置解算方法 | |
CN103604430A (zh) | 一种基于边缘化ckf重力辅助导航的方法 | |
CN108873029A (zh) | 一种实现导航接收机钟差建模的方法 | |
CN109186591A (zh) | 一种基于系统状态估计的sins/gps高精度重力扰动补偿方法 | |
KR20110085495A (ko) | 센서오차의 작동 중 자동교정 방법과 이를 이용한 관성항법장치 | |
CN109470241A (zh) | 一种具备重力扰动自主补偿功能的惯性导航系统及方法 | |
CN104613966B (zh) | 一种地籍测量事后数据处理方法 | |
CN106468554A (zh) | 一种非接触式的翻滚卫星的惯性参数的测算方法 | |
CN109521444B (zh) | 一种地壳运动gps水平速度场自适应最小二乘拟合推估算法 | |
CN103808349A (zh) | 矢量传感器的误差校正方法和装置 | |
US6298287B1 (en) | System and method of compensating for pressure sensor errors and noise in inertial vertical loop data | |
US8510079B2 (en) | Systems and methods for an advanced pedometer | |
Guo et al. | A robust SINS/USBL integrated navigation algorithm based on earth frame and right group error definition | |
CN107490373A (zh) | 基于编码器和惯性元器件融合的设备位姿估计方法及系统 | |
CN111308432B (zh) | 一种利用测速数据评估航天器测距数据精度的方法 | |
CN111679307B (zh) | 一种卫星定位信号解算方法及装置 |
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 |