CN108491633B - 一种基于标准椭球牛顿迭代法空间目标/诱饵测量数据产生方法 - Google Patents
一种基于标准椭球牛顿迭代法空间目标/诱饵测量数据产生方法 Download PDFInfo
- Publication number
- CN108491633B CN108491633B CN201810248660.XA CN201810248660A CN108491633B CN 108491633 B CN108491633 B CN 108491633B CN 201810248660 A CN201810248660 A CN 201810248660A CN 108491633 B CN108491633 B CN 108491633B
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- bait
- space target
- under
- target
- 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 27
- 238000005259 measurement Methods 0.000 title claims abstract description 19
- 238000006243 chemical reaction Methods 0.000 claims abstract description 4
- 238000013459 approach Methods 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 2
- 230000005484 gravity Effects 0.000 claims 1
- 238000004422 calculation algorithm Methods 0.000 abstract description 6
- 238000004088 simulation Methods 0.000 abstract description 5
- 238000012545 processing Methods 0.000 abstract description 4
- 238000012795 verification Methods 0.000 description 4
- 238000007689 inspection Methods 0.000 description 2
- 238000009827 uniform distribution Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于标准椭球牛顿迭代法空间目标/诱饵测量数据产生方法,该方法属于空间目标数据处理技术领域。由于空间目标/诱饵的实测数据较少,许多算法的性能就需要利用仿真数据来验证。为了得到更可靠的仿真数据,本发明中的地球模型采用更准确的标准椭球模型,并结合目标空间动力学方程、诱饵释放工程背景等,首先根据空间目标飞行中段动力学模型获得地心惯性坐标系下目标位置和速度信息,然后获得诱饵在地心惯性坐标系下相对空间目标的位置和速度信息,最后由标准椭球模型下地心惯性坐标系和北东上坐标系的转换关系,获得北东上坐标系空间目标/诱饵数据,最后获得极坐标下空间目标/诱饵测量数据。
Description
技术领域
本发明属于空间目标数据处理技术领域,适用于高速空间目标数据仿真***。
背景技术
由于空间目标工程背景比较复杂,涉及到的坐标系就包括地心地固(ECEF)坐标系、地心惯性(ECI)坐标系和北东上(ENU)坐标系等,空间目标在上升飞行到一定预定点之前,还要释放诱饵等虚假目标,空间目标和诱饵的实测数据较少,许多算法的性能就需要利用仿真数据来验证,而传统的处理方法是把地球模型看成是标准圆球,利用标准圆球模型进行简化处理生成相关的空间目标数据,再进行仿真分析验证。标准圆球模型具有算法处理简单等优点,但也具有和地球实际形状差距较大等缺点。如何充分利用更准确的标准椭球模型,再结合空间目标动力学方程、诱饵释放的工程背景等生成更准确的空间目标/诱饵数据,并利用这些数据对相关理论算法和工程模型进行检验验证具有更重要的实际意义。
发明内容
本发明的目的在于提供一种基于标准椭球牛顿迭代法的空间目标/诱饵测量数据产生方法,旨在生成更准确的空间目标数据,为相关理论算法和工程模型的检验验证提供数据基础。
本发明所述的基于标准椭球牛顿迭代法空间目标/诱饵测量数据产生方法,是指根据牛顿迭代法和空间目标飞行中段动力学方程获得ECI坐标系下空间目标位置和速度信息,然后在获得地心惯性坐标系下诱饵相对空间目标位置和速度信息的基础上,获得诱饵在地心惯性坐标系下的位置和速度信息,并由标准椭球模型下ECI坐标系和ENU坐标系的转换关系获得ENU坐标系下空间目标/诱饵位置和速度信息,最后利用ENU坐标系下的位置信息获得极坐标系下空间目标/诱饵径向距离、方位角和俯仰角测量数据。
本发明的方法具体包括以下步骤:
步骤1:利用空间目标动力学方程和牛顿迭代法获得空间目标在地心惯性坐标系下的位置和速度信息;
步骤2:确定地心惯性坐标系下空间目标椭圆轨道参数;
步骤3:利用椭圆轨道参数确定诱饵在地心惯性坐标系下相对空间目标的位置和速度,进而获得诱饵在地心惯性坐标系下的位置和速度信息;
步骤4:地球模型为标准椭球模型时,将地心惯性坐标系下的空间目标/诱饵位置和速度信息转换到雷达站北东上坐标系下,获得该坐标系下的空间目标/诱饵位置和速度信息;
步骤5:由雷达站北东上坐标系下的位置信息获得极坐标系下空间目标/诱饵径向距离、方位角和俯仰角测量数据。
本发明可为空间目标数据处理相关算法的性能验证提供更准确、更贴合工程背景的数据支持。
附图说明
图1是基于标准椭球牛顿迭代法空间目标/诱饵测量数据产生方法流程图。
具体实施方式
本发明的具体实施方式分以下几个步骤:
步骤1:地心惯性(ECI)坐标系下空间目标位置信息获得。
地球模型采用标准椭球体,由关机点时刻空间目标质心在WGS-84坐标系下的纬度、经度和大地高程(Bm、Lm和Hm),可得该时刻空间目标质心在ECEF坐标系下的坐标为
xe=(Nm+Hm)cosBmcosLm (1)
ye=(Nm+Hm)cosBmsinLm (2)
ze=[Nm(1-e2)+Hm]sinBm (3)
其中:e为地球第一偏心率,且
a为参考椭球体的长半轴,b为参考椭球体的短半轴;
Nm为空间目标质心所在点的卯酉圈曲率半径,且
将空间目标关机点时刻地心地固坐标系下的目标位置转换到地心惯性坐标系下,可得
其中:L0为关机点时刻地心地固坐标系的XF轴(OXF)与地心惯性坐标系的XI轴(OXI)相差的角度,由于一般假定空间目标关机点时刻地心地固坐标系和地心惯性坐标系是重合的,所以该时刻L0=0,此时ECEF和ECI坐标系下的目标坐标相同,经过时间t后OXI和OXF的夹角为
Lt=ω·t (7)
这里ω为地球自转角速度。
关机点时刻地心地固直角坐标系下空间目标的x、y、z轴速度可根据常识事先假定,并将其转换到ECI坐标系下可得
由关机点时刻空间目标在地心惯性坐标系下的位置和速度信息,利用牛顿迭代法给出的状态差分方程得到后续时刻地心惯性坐标系下空间目标位置和速度信息:
其中:μ为万有引力常数,rI为ECI坐标系下空间目标质心的地心距,
这里J2为地球二阶带谐系数,re为地球赤道半径;
步骤2:地心惯性坐标系下诱饵相对空间目标位置和速度的产生。
空间目标飞行到一定的预定点时,开始释放诱饵,假设诱饵的释放速度v0在vmin和vmax之间按均匀分布随机取值,诱饵相对空间目标释放速度的x、y、z轴分量分别为vx0、vy0和vz0,
其中:α、β在0到2π之间按均匀分布随机取值。
t时刻地心惯性坐标系下诱饵相对空间目标的位置和速度近似解为:
其中:
(1)空间目标轨道平面的方向数:
(3)空间目标在惯性坐标系下的速度vId及动量矩
(4)椭圆轨道长半轴aE
(5)椭圆轨道偏心率eE
(6)椭圆轨道短半轴bE和焦距cE
cE=aEeE (32)
(8)真近点角fd
在获得真近点角fd后,由式(21)可得ωz,进而由式(15)-式(20)可获得诱饵在地心惯性坐标系下相对空间目标的位置和速度。
步骤3:北东上坐标系(ENU)下空间目标/诱饵位置和速度信息的获得。
利用获得的地心惯性坐标系下诱饵相对空间目标的位置和速度信息,再结合由式(9)得到的空间目标在地心惯性坐标系下的位置和速度,即可获得诱饵在地心惯性坐标系下的位置和速度信息。当地球模型采用标准椭球模型时,将得到的地心惯性坐标系下空间目标/诱饵位置信息(xI,yI,zI)转换到雷达站北东上坐标系下:
从地心惯性坐标系到雷达站北东上坐标系的空间目标/诱饵速度信息转换关系为
其中:
步骤4:极坐标系下空间目标/诱饵径向距离、方位角和俯仰角测量数据的产生。
利用雷达站北东上坐标系下的空间目标/诱饵位置信息(xr,yr,zr)得到极坐标系下的空间目标/诱饵径向距离、方位角和俯仰角测量数据:
Claims (3)
1.一种基于标准椭球牛顿迭代法空间目标/诱饵测量数据产生方法,其特征在于,包括以下步骤:
步骤1:利用空间目标动力学方程和牛顿迭代法获得空间目标在地心惯性坐标系下的位置和速度信息;
步骤2:确定地心惯性坐标系下空间目标椭圆轨道参数,具体包括:空间目标轨道平面的方向数、空间目标质心的地心距及其变化率、空间目标速度、动量矩、椭圆轨道长半轴、椭圆轨道偏心率、椭圆轨道短半轴、焦距、偏近点角和真近点角fd;
步骤3:利用椭圆轨道参数确定诱饵在地心惯性坐标系下相对空间目标的位置和速度,进而获得诱饵在地心惯性坐标系下的位置和速度信息,具体为:
根据真近点角fd和诱饵相对空间目标的释放速度(vx0,vy0,vz0)确定后续t时刻诱饵在地心惯性坐标系下相对空间目标的位置和速度近似解:
其中:eE为椭圆轨道偏心率,
式中μ为万有引力常数,rI为地心惯性坐标系下空间目标质心的地心距;
利用获得的诱饵在地心惯性坐标系下相对空间目标的位置和速度,再结合空间目标在地心惯性坐标系下的位置和速度,获得诱饵在地心惯性坐标系下的位置和速度信息;
步骤4:地球模型为标准椭球模型时,将地心惯性坐标系下的空间目标/诱饵位置和速度信息转换到雷达站北东上坐标系下,获得该坐标系下的空间目标/诱饵位置和速度信息,具体为:
当地球模型采用标准椭球模型时,将获得的地心惯性坐标系下的空间目标/诱饵位置信息(xI,yI,zI)转换到雷达站北东上坐标系下:
式中Lr为WGS-84坐标系中的雷达站心的经度,ω为地球自转角速度;
从地心惯性坐标系到雷达站北东上坐标系的空间目标/诱饵速度信息转换关系为
步骤5:由雷达站北东上坐标系下的位置信息获得极坐标系下空间目标/诱饵径向距离、方位角和俯仰角测量数据。
2.如权利要求1所述的一种基于标准椭球牛顿迭代法空间目标/诱饵测量数据产生方法,其特征在于,步骤1具体为:
由关机点时刻空间目标质心在WGS-84坐标系下的纬度Bm、经度Lm和大地高程Hm确定其在地心地固坐标系下的位置信息(xe,ye,ze),由于关机点时刻地心地固坐标系和地心惯性坐标系是重合的,所以有地心惯性坐标系下的位置信息(xI,yI,zI)和(xe,ye,ze)相等,而该时刻地心惯性坐标系下的空间目标速度信息为
由关机点时刻空间目标在地心惯性坐标系下的位置和速度信息,利用牛顿迭代法给出的状态差分方程得到后续时刻地心惯性坐标系下空间目标位置和速度信息:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810248660.XA CN108491633B (zh) | 2018-03-24 | 2018-03-24 | 一种基于标准椭球牛顿迭代法空间目标/诱饵测量数据产生方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810248660.XA CN108491633B (zh) | 2018-03-24 | 2018-03-24 | 一种基于标准椭球牛顿迭代法空间目标/诱饵测量数据产生方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108491633A CN108491633A (zh) | 2018-09-04 |
CN108491633B true CN108491633B (zh) | 2021-10-01 |
Family
ID=63337519
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810248660.XA Active CN108491633B (zh) | 2018-03-24 | 2018-03-24 | 一种基于标准椭球牛顿迭代法空间目标/诱饵测量数据产生方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108491633B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103941740A (zh) * | 2014-04-15 | 2014-07-23 | 北京控制工程研究所 | 一种考虑地球椭率的多轴机动成像卫星偏航姿态控制方法 |
CN104655135A (zh) * | 2015-02-09 | 2015-05-27 | 南京邮电大学 | 一种基于地标识别的飞行器视觉导航方法 |
CN106197434A (zh) * | 2016-06-30 | 2016-12-07 | 中国电子科技集团公司第五十四研究所 | 基于地面目标点位置的卫星姿态角的计算方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104697900A (zh) * | 2013-12-05 | 2015-06-10 | 上海梭伦信息科技有限公司 | 真空、控温条件下界面张力和接触角值的测量方法及装置 |
KR101723148B1 (ko) * | 2015-10-20 | 2017-04-05 | 임채영 | 태양광 발전용 태양전지 유닛 및 그 제조방법 |
-
2018
- 2018-03-24 CN CN201810248660.XA patent/CN108491633B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103941740A (zh) * | 2014-04-15 | 2014-07-23 | 北京控制工程研究所 | 一种考虑地球椭率的多轴机动成像卫星偏航姿态控制方法 |
CN104655135A (zh) * | 2015-02-09 | 2015-05-27 | 南京邮电大学 | 一种基于地标识别的飞行器视觉导航方法 |
CN106197434A (zh) * | 2016-06-30 | 2016-12-07 | 中国电子科技集团公司第五十四研究所 | 基于地面目标点位置的卫星姿态角的计算方法 |
Non-Patent Citations (2)
Title |
---|
Spatial Clutter Measurement Density Estimation in Nonhomogeneous Measurement Spaces;Kim, Woo Chan等;《2015 18TH INTERNATIONAL CONFERENCE ON INFORMATION FUSION (FUSION)》;20150709;1772-1777 * |
基于飞机弹射座椅姿态求解算法的仿真研究与实现;赵虹;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20120715(第7期);C031-193 * |
Also Published As
Publication number | Publication date |
---|---|
CN108491633A (zh) | 2018-09-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106124170B (zh) | 一种基于高精度姿态信息的相机光轴指向计算方法 | |
CN107544067B (zh) | 一种基于高斯混合近似的高超声速再入飞行器跟踪方法 | |
CN106643804B (zh) | 一种预先确定星载激光测高仪的足印位置的方法 | |
CN104568652B (zh) | 一种高精度测量临近空间大气密度的方法及测量装置 | |
CN107607127B (zh) | 一种基于外场的星敏感器内部参数标定及精度快速验证*** | |
CN107300700B (zh) | 敏捷合成孔径雷达卫星聚束模式姿态机动需求计算方法 | |
CN110646782A (zh) | 一种基于波形匹配的星载激光在轨指向检校方法 | |
CN109932711B (zh) | 一种雷达测量的大气折射修正方法 | |
Kartsan et al. | Applying filtering for determining the angular orientation of spinning objects during interference | |
CN106767928A (zh) | 一种自适应快速传递对准方法 | |
CN109975836B (zh) | 解算ccd图像地面位置方法、装置、电子设备及介质 | |
CN115343744A (zh) | 空中运动目标的光学单双星联合星上定位方法及*** | |
CN107991696A (zh) | 针对终端载体为高轨卫星的卫星导航信号仿真方法 | |
CN108491633B (zh) | 一种基于标准椭球牛顿迭代法空间目标/诱饵测量数据产生方法 | |
CN116992700B (zh) | 一种物流无人机导航精度确定的方法及设备 | |
Huang et al. | Research on UAV flight performance test method based on dual antenna GPS/ins integrated system | |
US11137472B2 (en) | Amplitude based centroiding for angle of arrival estimation | |
CN104931180B (zh) | 一种六自由度气浮台干扰力矩确定方法 | |
CN104391311B (zh) | 基于gps广播数据的星上无源定位方法 | |
Moyano Cano | Quadrotor UAV for wind profile characterization | |
US20210311157A1 (en) | Measuring circuit, measuring device, and program | |
CN109993795B (zh) | 解算sar图像地面位置方法、装置、电子设备及介质 | |
Tian et al. | Launch point estimation with a single passive sensor without trajectory state estimation | |
Lee et al. | A simple prediction method of ballistic missile trajectory to designate search direction and its verification using a testbench | |
Protsenko et al. | Unmanned aerial vehicle positioning by data from pocket device sensors |
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 |