CN112540367A - 空间目标雷达定轨实时识别方法、设备和存储介质 - Google Patents
空间目标雷达定轨实时识别方法、设备和存储介质 Download PDFInfo
- Publication number
- CN112540367A CN112540367A CN202011390432.XA CN202011390432A CN112540367A CN 112540367 A CN112540367 A CN 112540367A CN 202011390432 A CN202011390432 A CN 202011390432A CN 112540367 A CN112540367 A CN 112540367A
- Authority
- CN
- China
- Prior art keywords
- target
- trace data
- track
- orbit determination
- space
- 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 39
- 230000006872 improvement Effects 0.000 claims abstract description 12
- 238000012937 correction Methods 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 9
- 230000009466 transformation Effects 0.000 claims description 7
- 239000013598 vector Substances 0.000 claims description 6
- 239000000126 substance Substances 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 abstract description 6
- 238000012544 monitoring process Methods 0.000 abstract description 6
- 230000035484 reaction time Effects 0.000 abstract description 5
- 238000005516 engineering process Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 239000003795 chemical substances by application Substances 0.000 description 2
- 238000002592 echocardiography Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000000630 rising effect Effects 0.000 description 2
- 230000001174 ascending effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
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
- 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
- G01S13/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本申请公开了空间目标雷达定轨实时识别方法、设备和存储介质,方法包括:获取目标点迹数据;对目标点迹数据进行坐标转换;判断目标点迹数据中是否存在轨道根数;根据目标点迹数据执行初轨确定,得到初始轨道根数;根据目标点迹数据执行轨道改进,得到精密轨道根数;根据目标点迹数据执行轨道匹配。通过初轨确定、轨道改进、轨道匹配运算,能够快速有效对雷达观测到的空间目标进行识别,并成功和目标数据库中对应目标匹配,有效缩短了空间目标从发现到识别所需要的反应时间,有效增强了预警雷达对空间目标的预警监视能力。
Description
技术领域
本申请涉及信息技术领域,尤其涉及空间目标雷达定轨实时识别方法、装置、设备和存储介质。
背景技术
对空间目标的快速捕获识别与高精度跟踪监视是利用并控制空间资源的前提和基础。自20世纪50年代开始雷达目标识别的相关研究就开始出现,其主要原理是通过对雷达观测信号回波进行分析和处理,从中提炼典型数据或者特征标识,从而对目标属性进行识别和判断,早期的雷达目标识别研究是基于目标RCS的测量来进行,随着技术发展,越来越多可用于目标识别的信息被不断挖掘出来,相关研究和应用技术也不断出现,例如基于回波调制谱特征的目标识别技术、基于目标极点分布特征的识别技术、基于目标回波极化特征的目标识别技术、高分辨力雷达目标成像识别等等。
然而,现有的空间目标雷达定轨实时识别方法,空间目标从发现到识别所需要的反应时间长,导致预警雷达对空间目标的预警监视能力差。
发明内容
本申请的主要目的在于解决现有的空间目标雷达定轨实时识别方法反应时间长、预警监视能力差的技术问题。
一种空间目标雷达定轨实时识别方法,其特征在于,依次包括以下步骤:
S1、获取目标点迹数据;
S2、对目标点迹数据进行坐标转换;
S3、判断目标点迹数据的数量是否大于或等于N0;若为是,则执行S4;若为否,则执行S1;
S4、判断目标点迹数据中是否存在轨道根数;若为否,则执行S5;若为是,则执行S6;
S5、判断目标点迹数据的数量是否大于或等于N1;若为是,则执行S13;若为否,则执行S7;
S6、判断目标点迹数据的数量是否大于或等于N2;若为是,则执行S8;若为否,则执行S7;
S7、根据目标点迹数据执行初轨确定,得到初始轨道根数,执行S1;
S8、判断目标点迹数据的数量是否是N3的倍数;若为是,则执行S9;若为否,则执行S1;
S9、根据目标点迹数据执行轨道改进,得到精密轨道根数;
S10、判断目标点迹数据的数量是否大于或等于N1;若为是,则执行S11;若为否,则执行S1;
S11、根据目标点迹数据执行轨道匹配;
S12、判断轨道匹配是否成功;若为是,则输出空间目标信息;若为否,则执行S13;
S13、提示识别失败;
其中,N0、N1、N2、N3分别为正整数,且N1>N2>N3>N0,且N1和N2分别是N3的整数倍。
可选地,所述轨道根数包含半长轴、偏心率、倾角、升交点赤经、近地点幅角和升交角距。
可选地,所述根据目标点迹数据执行初轨确定,得到初始轨道根数,包括:根据目标点迹数据执行初轨确定,得到空间目标的运行速度。
可选地,所述根据目标点迹数据执行初轨确定,得到空间目标的运行速度,包括以下步骤:
获取3个连续时刻t1、t2和t3的目标点迹数据;其中,t1<t2<t3;
从3个连续时刻的目标点迹数据中确定3个时刻的地心位置矢量r1、r2和r3;
根据3个时刻的地心位置矢量r1、r2和r3,确定空间目标在t2时刻的所述运行速度。
可选地,所述运行速度的表达式为:
v2=-d1r1+d2r2+d3r3
其中,
v2为所述运行速度。
可选地,所述根据目标点迹数据执行轨道改进,得到精密轨道根数,包括:
获取航迹数据;
以初始轨道根数为初始值,进行有摄星历计算,得到与实际观测时刻对应的预报数据;
根据所述航迹数据和所述预报数据,得到当前的观测残差;
根据当前的观测残差得到修正的轨道根数,并且计算最近两次观测残差的均方根;
根据最近两次观测残差的均方根,得到迭代参数;若所述迭代参数小于预设阈值,则将修正的轨道根数设定为所述精密轨道根数。
可选地,所述观测残差的表达式为:
ΔYj=Y-Yj
其中,j表示迭代次数,ΔYj为第j轮迭代对应的观测残差,Y为航迹数据,Yj为第j轮迭代中计算得到的预报数据。
可选地,所述根据所述观测残差得到修正的轨道根数,包括:
根据所述观测残差,得到轨道根数修正值;
所述轨道根数修正值的表达式为:
ΔX=(BTWB)-1BTWΔY
其中,B为观测残差表达式的系数矩阵,它的求解和实际选取的摄动力计算模型有关,W为权矩阵,ΔY为观测残差;
根据所述轨道根数修正值,得到修正的轨道根数。
基于相同的技术构思,本申请还提供了一种计算机设备,包括输入输出单元、存储器和处理器,所述存储器中存储有计算机可读指令,所述计算机可读指令被所述处理器执行时,使得所述处理器执行如上述的空间目标雷达定轨实时识别方法中的步骤。
基于相同的技术构思,本申请还提供了一种存储有计算机可读指令的存储介质,所述计算机可读指令被一个或多个处理器执行时,使得一个或多个处理器执行如上述的空间目标雷达定轨实时识别方法中的步骤。
本申请的有益效果:通过初轨确定、轨道改进、轨道匹配运算,能够快速有效对雷达观测到的空间目标进行识别,并成功和目标数据库中对应目标匹配,有效缩短了空间目标从发现到识别所需要的反应时间,有效增强了预警雷达对空间目标的预警监视能力。
附图说明
图1为本申请实施例中空间目标雷达定轨实时识别方法的流程示意图。
图2为本申请实施例中计算机设备的结构示意图。
具体实施方式
应当理解,此处所描述的具体实施例仅用以解释本申请,并不用于限定本申请。
本技术领域技术人员可以理解,除非特意声明,这里使用的单数形式“一”、“一个”、“所述”和“该”也可以包括复数形式。应该进一步理解的是,本申请的说明书中使用的措辞“包括”是指存在所述特征、程序、步骤、操作、元件和/或组件,但是并不排除存在或添加一个或多个其他特征、程序、步骤、操作、元件、组件和/或它们的组。
图1为本申请一些实施方式中一种空间目标雷达定轨实时识别方法的流程图,如图1所示,依次包括以下步骤S1-S13。
S1、获取目标点迹数据,执行S2;
一些实施方式中,在雷达任务规划席位通过内网UDP协议接收来自雷达显控席的目标点迹数据。
S2、对目标点迹数据进行坐标转换,执行S3。
雷达观测基于地球坐标系,空间目标轨道是基于天球坐标系,这就涉及到地球坐标系与天球坐标系之间的转换。
协议天球坐标系由国际天文联合会(IAU)和国际地球自转和参考系服务(IERS)发布,目前采用的是国际天球参考系(ICRS)。依据坐标原点的不同,ICRS可分为太阳系质心天球参考系(BCRS)和地球质心天球参考系(GCRS)。BCRS用于计算行星的运动轨道,编制星表;GCRS用于计算卫星轨道,编制卫星星历。ICRS由国际天球参考框架(ICRF)来实现。1997年IAU第23届大会上,通过并决定自1998年1月1日起,在天文研究、空间探测、大地测量以及地球动力学等领域中采用ICRS。
协议地球坐标系由国际地球参考系(ITRS)实现。GCRS是一个相当好的准惯性系,卫星的轨道计算一般都是在GCRS中进行。这就必须涉及到GCRS与ITRS间的坐标转换问题。
一些实施方式中,S2包含以下步骤S21:将目标点迹数据由站心地平坐标系Xh转换为ITRS坐标系XGO。
可选地,转换表达式为:
其中,Rn(θ)表示绕n轴(n=1,2,3分别对应x,y,z三个轴)转动角θ的坐标变换矩阵,P2=diag[1,-1,1],λ、为测站的天文经纬度,而测站在ITRS坐标系中的直角坐标XGOC由测站的地理经纬高坐标(B,L,H)计算。站心地平坐标系与ITRS坐标系的定义如表1所示。
表1
一些实施方式中,S2包含以下步骤S21:将目标点迹数据由ITRS坐标系XGO转换为GCRS坐标系XGCRS。
ITRS与GCRS的转换早期是基于春分点的。目前,IERS(2010)建议使用基于无旋转原点(NRO)的转换方法。转换过程中涉及到两个中间坐标系:地球中间坐标系(TIRS)和天球中间坐标系(CIRS),它们的定义见表2。
表2
在t时刻,ITRS和GCRS的转换是两个三维直角坐标系间的转换,可选地,转换表达式为:
XGCRS=M(t)RCIO(t)W(t)XITRS (2)
式中,M(t)、RCIO(t)和W(t)分别表示由于CIP在GCRS中的运动(岁差章动)、地球的自转以及CIP在ITRS中的运动(极移)引起的旋转矩阵。
S3、判断目标点迹数据的数量是否大于或等于N0;若为是,则执行S4;若为否,则执行S1。
S4、判断目标点迹数据中是否存在轨道根数;若为否,则执行S5;若为是,则执行S6。
一些实施方式中,轨道根数包含半长轴a、偏心率e、倾角i、升交点赤经Ω、近地点幅角ω和升交角距L。
轨道根数,或称轨道要素或轨道参数,是用来描述天体在其轨道运行状态的一组参数。通常情况下指的是用经典万有引力定律描述天体按圆锥曲线运动时所必需的6个参数。
S5、判断目标点迹数据的数量是否大于或等于N1;若为是,则执行S13;若为否,则执行S7。
S6、判断目标点迹数据的数量是否大于或等于N2;若为是,则执行S8;若为否,则执行S7。
S7、根据目标点迹数据执行初轨确定,得到初始轨道根数,执行S1。
初轨确定是用少量观测数据粗略确定空间目标轨道的过程,主要用于刚入轨或初次捕获的空间目标的定轨。通过初轨确定可以迅速提供空间目标运动的粗略情况,为轨道改进提供初值。
一些实施方式中,所述根据目标点迹数据执行初轨确定,得到初始轨道根数,包含以下步骤:根据目标点迹数据执行初轨确定,得到空间目标的运行速度。具体地,包含以下步骤。
S71、获取3个连续时刻t1、t2和t3的目标点迹数据;其中,t1<t2<t3;
S72、从3个连续时刻的目标点迹数据中确定3个时刻的地心位置矢量r1、r2和r3;
S73、根据3个时刻的地心位置矢量r1、r2和r3,确定空间目标在t2时刻的运行速度V2。
具体地,运行速度V2的表达式为:
v2=-d1r1+d2r2+d3r3 (3)
其中,
S8、判断目标点迹数据的数量是否是N3的倍数;若为是,则执行S9;若为否,则执行S1。
S9、根据初始轨道根数执行轨道改进,得到精密轨道根数,执行S10。
轨道改进是一种精密测定天体轨道的方法。这种方法以天体的某一初始轨道为依据,利用尽可能多的观测资料,逐次改进轨道要素,最后求出天体的精密轨道。
一些实施方式中,步骤S9包含以下步骤S91-S95:
S91、获取航迹数据。
S92、以初始轨道根数为初始值,进行有摄星历计算,得到与实际观测时刻对应的预报数据。
S93、根据所述航迹数据和所述预报数据,得到当前的观测残差。
可选地,观测残差的表达式为:
ΔYj=Y-Yj (4)
其中,j表示迭代次数,ΔYj为第j轮迭代对应的观测残差,Y为航迹数据,Yj为第j轮迭代中计算得到的预报数据。
S94、根据当前的观测残差得到修正的轨道根数,并且计算最近两次观测残差的均方根。
可选地,所述根据所述观测残差得到修正的轨道根数,包含以下步骤S941-S942:
S941、根据所述观测残差,得到轨道根数修正值。
进一步地,轨道根数修正值的表达式为:
ΔX=(BTWB)-1BTWΔY (5)
其中,B为观测残差表达式的系数矩阵,它的求解和实际选取的摄动力计算模型有关,W为权矩阵,ΔY为观测残差。
S942、根据所述轨道根数修正值,得到修正的轨道根数。
进一步地,修正的轨道根数的表达式为:
Xj+1=Xj-ΔXj (6)
可选地,观测残差的均方根的表达式为:
其中,m、n为预设常数。
S95、根据最近两次观测残差的均方根,得到迭代参数;若所述迭代参数小于预设阈值,则将修正的轨道根数设定为所述精密轨道根数。
可选地,迭代参数的表达式为:
比较两次残差的均方根,如果δ<μ则迭代终止。其中μ为小量,一般可取为0.01。如果δ>μ则将该轮修正的轨道根数作为输入再次带入有摄星历计算,继续迭代直到最终满足δ<μ为止。
S10、判断目标点迹数据的数量是否大于或等于N1;若为是,则执行S11;若为否,则执行S1。
S11、根据所述精密轨道根数执行轨道匹配,执行S12。
空间目标是沿着固有轨道运动的,它的六个轨道根数是决定其运行规律的主要参数,也是进行空间目标识别的基础。因此,空间目标的识别可以通过将观测数据与已知目标的轨道根数样本进行匹配来实现。
S12、判断轨道匹配是否成功;若为是,则输出空间目标信息;若为否,则执行S13。
一些实施方式中,轨道匹配的表达式为:
其中,δi为待识别目标相对样本目标的轨道根数误差,i=1,…,6分别表示半长轴a、偏心率e、倾角i、升交点赤经Ω、近地点幅角ω和升交角距L(可与真近点角进行互算)相对误差。wi(i=1,…,6)为六个加权系数,满足w1+w2+…+w6=1,T为比较门限。如果上式满足,则目标匹配成功。
首先,定轨相对误差δi由下式计算:
权重系数wi和门限T的设置是轨道匹配的关键,本申请利用NASA公布的空间目标数据库,根据某型雷达的威力范围,从中选出N个近地点高度小于2000km的目标作为样本库。设置雷达的站址和跟踪空域范围,依据雷达的测量精度,对这N个目标按3σ准则仿真得到带观测误差的距离、方位和俯仰数据(以下记为RAE数据)。由RAE数据定轨获得这N个目标的轨道根数,然后按表达式(5)计算相对误差,这样对每个轨道根数都得到一个相对误差序列接着计算相对误差序列的标准差σi:
则各轨道根数的权重系数wi可由下式计算获得:
有了wi后,可以对这N个目标分别计算参数Tj(j=1,…,N)。
然后计算它的均值mT和标准差σT,则门限T可由下式计算获得:
T=mT+(3~5)σT (14)
从上式也可以看出,轨道改进精度越高(即δi越小),进行轨道匹配的门限T就越小。
S13、提示识别失败。
其中,N0、N1、N2、N3分别为正整数,且N1>N2>N3>N0,且N1和N2分别是N3的整数倍。优选地,N0等于3。
上述实施例中,通过初轨确定、轨道改进、轨道匹配运算,能够快速有效对雷达观测到的空间目标进行识别,并成功和目标数据库中对应目标匹配,有效缩短了空间目标从发现到识别所需要的反应时间,有效增强了预警雷达对空间目标的预警监视能力。
基于相同的技术构思,本申请还提供了一种计算机设备,如图2所示,该计算机设备包括输入输出单元31、处理器32和存储器33,所述存储器33中存储有计算机可读指令,所述计算机可读指令被所述处理器32执行时,使得所述处理器执行上述各实施方式中的所述的空间目标雷达定轨实时识别方法的步骤。
基于相同的技术构思,本申请还提供了一种存储有计算机可读指令的存储介质,所述计算机可读指令被一个或多个处理器执行时,使得一个或多个处理器执行上述各实施方式中的所述的空间目标雷达定轨实时识别方法的步骤。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到上述实施例方法可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件,但很多情况下前者是更佳的实施方式。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质(如ROM/RAM)中,包括若干指令用以使得一台终端(可以是手机,计算机,服务器或者网络设备等)执行本申请各个实施例所述的方法。
上面结合附图对本申请的实施例进行了描述,但是本申请并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,而不是限制性的,本领域的普通技术人员在本申请的启示下,在不脱离本申请宗旨和权利要求所保护的范围情况下,还可做出很多形式,凡是利用本申请说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,这些均属于本申请的保护之内。
Claims (10)
1.一种空间目标雷达定轨实时识别方法,其特征在于,依次包括以下步骤:
S1、获取目标点迹数据;
S2、对目标点迹数据进行坐标转换;
S3、判断目标点迹数据的数量是否大于或等于N0;若为是,则执行S4;若为否,则执行S1;
S4、判断目标点迹数据中是否存在轨道根数;若为否,则执行S5;若为是,则执行S6;
S5、判断目标点迹数据的数量是否大于或等于N1;若为是,则执行S13;若为否,则执行S7;
S6、判断目标点迹数据的数量是否大于或等于N2;若为是,则执行S8;若为否,则执行S7;
S7、根据目标点迹数据执行初轨确定,得到初始轨道根数,执行S1;
S8、判断目标点迹数据的数量是否是N3的倍数;若为是,则执行S9;若为否,则执行S1;
S9、根据目标点迹数据执行轨道改进,得到精密轨道根数;
S10、判断目标点迹数据的数量是否大于或等于N1;若为是,则执行S11;若为否,则执行S1;
S11、根据目标点迹数据执行轨道匹配;
S12、判断轨道匹配是否成功;若为是,则输出空间目标信息;若为否,则执行S13;
S13、提示识别失败;
其中,N0、N1、N2、N3分别为正整数,且N1>N2>N3>N0,且N1和N2分别是N3的整数倍。
2.根据权利要求1所述的空间目标雷达定轨实时识别方法,其特征在于,
所述轨道根数包含半长轴、偏心率、倾角、升交点赤经、近地点幅角和升交角距。
3.根据权利要求1所述的空间目标雷达定轨实时识别方法,其特征在于,
所述根据目标点迹数据执行初轨确定,得到初始轨道根数,包括:根据目标点迹数据执行初轨确定,得到空间目标的运行速度。
4.根据权利要求3所述的空间目标雷达定轨实时识别方法,其特征在于,
所述根据目标点迹数据执行初轨确定,得到空间目标的运行速度,包括以下步骤:
获取3个连续时刻t1、t2和t3的目标点迹数据;其中,t1<t2<t3;
从3个连续时刻的目标点迹数据中确定3个时刻的地心位置矢量r1、r2和r3;
根据3个时刻的地心位置矢量r1、r2和r3,确定空间目标在t2时刻的所述运行速度。
6.根据权利要求1所述的空间目标雷达定轨实时识别方法,其特征在于,
所述根据目标点迹数据执行轨道改进,得到精密轨道根数,包括:
获取航迹数据;
以初始轨道根数为初始值,进行有摄星历计算,得到与实际观测时刻对应的预报数据;
根据所述航迹数据和所述预报数据,得到当前的观测残差;
根据当前的观测残差得到修正的轨道根数,并且计算最近两次观测残差的均方根;
根据最近两次观测残差的均方根,得到迭代参数;若所述迭代参数小于预设阈值,则将修正的轨道根数设定为所述精密轨道根数。
7.根据权利要求6所述的空间目标雷达定轨实时识别方法,其特征在于,
所述观测残差的表达式为:
ΔYj=Y-Yj
其中,j表示迭代次数,ΔYj为第j轮迭代对应的观测残差,Y为航迹数据,Yj为第j轮迭代中计算得到的预报数据。
8.根据权利要求6所述的空间目标雷达定轨实时识别方法,其特征在于,
所述根据所述观测残差得到修正的轨道根数,包括:
根据所述观测残差,得到轨道根数修正值;
所述轨道根数修正值的表达式为:
ΔX=(BTWB)-1BTWΔY
其中,B为观测残差表达式的系数矩阵,W为权矩阵,ΔY为观测残差;
根据所述轨道根数修正值,得到修正的轨道根数。
9.一种计算机设备,其特征在于,包括输入输出单元、存储器和处理器,所述存储器中存储有计算机可读指令,所述计算机可读指令被所述处理器执行时,使得所述处理器执行如权利要求1至8中的任一所述的空间目标雷达定轨实时识别方法中的步骤。
10.一种存储有计算机可读指令的存储介质,其特征在于,所述计算机可读指令被一个或多个处理器执行时,使得一个或多个处理器执行如权利要求1至8中的任一所述的空间目标雷达定轨实时识别方法中的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011390432.XA CN112540367B (zh) | 2020-12-02 | 2020-12-02 | 空间目标雷达定轨实时识别方法、设备和存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011390432.XA CN112540367B (zh) | 2020-12-02 | 2020-12-02 | 空间目标雷达定轨实时识别方法、设备和存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112540367A true CN112540367A (zh) | 2021-03-23 |
CN112540367B CN112540367B (zh) | 2023-11-14 |
Family
ID=75015395
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011390432.XA Active CN112540367B (zh) | 2020-12-02 | 2020-12-02 | 空间目标雷达定轨实时识别方法、设备和存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112540367B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115061106A (zh) * | 2022-08-16 | 2022-09-16 | 中国科学院空天信息创新研究院 | 方锥空域内空间目标可见性确定方法、装置、设备及介质 |
CN115096319A (zh) * | 2022-08-24 | 2022-09-23 | 航天宏图信息技术股份有限公司 | 一种基于光学测角数据的星链卫星初轨确定方法和装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009019984A (ja) * | 2007-07-11 | 2009-01-29 | Nec Corp | 目標観測レーダ装置及び目標追尾方法 |
JP2014013180A (ja) * | 2012-07-04 | 2014-01-23 | Mitsubishi Electric Corp | レーダ処理装置 |
CN105224737A (zh) * | 2015-09-22 | 2016-01-06 | 中国人民解放军63921部队 | 一种空间目标轨道改进初值修正方法 |
-
2020
- 2020-12-02 CN CN202011390432.XA patent/CN112540367B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009019984A (ja) * | 2007-07-11 | 2009-01-29 | Nec Corp | 目標観測レーダ装置及び目標追尾方法 |
JP2014013180A (ja) * | 2012-07-04 | 2014-01-23 | Mitsubishi Electric Corp | レーダ処理装置 |
CN105224737A (zh) * | 2015-09-22 | 2016-01-06 | 中国人民解放军63921部队 | 一种空间目标轨道改进初值修正方法 |
Non-Patent Citations (3)
Title |
---|
JOACHIM ENDER 等: "Radar techniques for space situational awareness", 2011 12TH INTERNATIONAL RADAR SYMPOSIUM (IRS), pages 21 - 26 * |
陶勇 等: "基于轨道匹配和改进的空间目标识别方法", 《雷达科学与技术》, vol. 4, no. 03, 30 June 2006 (2006-06-30), pages 134 - 139 * |
黄晓斌 等: "空间目标的雷达定轨坐标转换问题", 《导航定位学报》, vol. 8, no. 4, 31 August 2020 (2020-08-31), pages 39 - 43 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115061106A (zh) * | 2022-08-16 | 2022-09-16 | 中国科学院空天信息创新研究院 | 方锥空域内空间目标可见性确定方法、装置、设备及介质 |
CN115061106B (zh) * | 2022-08-16 | 2022-11-11 | 中国科学院空天信息创新研究院 | 方锥空域内空间目标可见性确定方法、装置、设备及介质 |
CN115096319A (zh) * | 2022-08-24 | 2022-09-23 | 航天宏图信息技术股份有限公司 | 一种基于光学测角数据的星链卫星初轨确定方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN112540367B (zh) | 2023-11-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Tisdale et al. | Autonomous UAV path planning and estimation | |
CN108279576B (zh) | 一种复合轴目标跟踪仿真测试*** | |
Pasqualetto Cassinis et al. | Cnn-based pose estimation system for close-proximity operations around uncooperative spacecraft | |
CN112540367A (zh) | 空间目标雷达定轨实时识别方法、设备和存储介质 | |
CN111551183A (zh) | 基于天基光学观测数据的geo目标多点择优短弧定轨方法 | |
CN109612438B (zh) | 一种虚拟共面条件约束下的空间目标初轨确定方法 | |
Xie et al. | Dynamic IBVS of a rotary wing UAV using line features | |
CN110187337B (zh) | 一种基于ls和neu-ecef时空配准的高机动目标跟踪方法及*** | |
CN112924997B (zh) | 目标终端定位方法、装置和存储介质及电子设备 | |
CN115790575B (zh) | 一种基于多星协同无源探测的巨型星座目标跟踪方法 | |
Opromolla et al. | Uncooperative spacecraft relative navigation with LIDAR-based unscented Kalman filter | |
CN107300700B (zh) | 敏捷合成孔径雷达卫星聚束模式姿态机动需求计算方法 | |
CN115343744A (zh) | 空中运动目标的光学单双星联合星上定位方法及*** | |
CN111121787A (zh) | 一种基于遥感图像的自主初轨确定方法 | |
CN114137589A (zh) | 无人机飞行定位方法、装置、设备及计算机可读存储介质 | |
CN113932796A (zh) | 高精地图车道线生成方法、装置和电子设备 | |
CN107883925B (zh) | 一种导航星座星间观测目标卫星图像模拟方法 | |
Rohani et al. | Vehicular cooperative map matching | |
CN115683170B (zh) | 基于雷达点云数据融合误差的校准方法 | |
CN114993341B (zh) | 一种基于天基测量的运载火箭弹道估计方法及装置 | |
Zhang et al. | Centimeter-level positioning by instantaneous lidar-aided GNSS ambiguity resolution | |
Li et al. | Sequential unbiased converted measurement non‐linear filter with range rate in ECEF coordinates | |
CN115273070A (zh) | 无人车初始位姿的获取方法、装置、电子设备及存储介质 | |
Waggoner et al. | Analysis of tdoa/fdoa state estimation accuracy of cislunar objects for space situational awareness | |
Dai et al. | Unmanned ground vehicle‐unmanned aerial vehicle relative navigation robust adaptive localization algorithm |
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 |