CN102927851B - 一种基于轨迹在线规划的末制导方法 - Google Patents
一种基于轨迹在线规划的末制导方法 Download PDFInfo
- Publication number
- CN102927851B CN102927851B CN201210472131.0A CN201210472131A CN102927851B CN 102927851 B CN102927851 B CN 102927851B CN 201210472131 A CN201210472131 A CN 201210472131A CN 102927851 B CN102927851 B CN 102927851B
- Authority
- CN
- China
- Prior art keywords
- gamma
- track
- tan
- angle
- sum
- 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.)
- Expired - Fee Related
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种基于轨迹在线规划的末制导方法,属于制导技术领域。本发明首先建立飞行器的运动学和动力学模型,然后采用两段相切圆弧作为侧向参考轨迹,进行初始时刻侧向轨迹设计,再通过旋转平移坐标变换,将任意时刻侧向轨迹问题转化为初始时刻侧向轨迹问题;由侧向参考轨迹得到的地面剩余航程及俯冲段弹道特性,设计纵向参考轨迹为一段椭圆曲线;根据得到的当前实时俯冲段纵向和侧向轨迹计算末制导的倾侧角μ和攻角α,并将其输入飞行器模型,对飞行器轨迹进行实时规划调整,使其满足期望的终端条件,从而实现末制导。本发明计算量小,优化效率高,计算时间短;在线规划的俯冲段轨迹满足终端约束并且可以控制末速。
Description
技术领域
本发明涉及一种基于轨迹在线规划的末制导方法,属于制导技术领域。
背景技术
再入飞行器的俯冲段为整个再入过程的末段,该阶段通常执行精确打击任务。制导***的主要目的是产生合适的指令,使得打击点的脱靶量为零。随着制导技术的快速发展,对制导***的要求也越来越高。
现代战争中,很多任务不仅要求导弹能够准确命中目标,还期望命中目标时具有特定的姿态。例如:要求反导导弹头对头对来袭导弹进行拦截;希望反坦克导弹能够垂直命中薄弱的前装甲;反舰导弹期望能从侧面对舰船进行攻击。针对带落角约束的制导律,国内外学者进行了大量研究。研究主要分为两类:第一类主要是具有偏置项的比例导引律,Kim等人对该方法进行了相应的研究,提出了形式非常简单的制导律,不过这种制导律局限于目标是固定的情况;第二类主要是基于最优控制理论,通过最小化性能指标求出最优问题解得形式,Cho等人对此类方法进行了研究,特别针对固定或运动速度较慢的目标提出了能够同时控制落角和攻击时间的制导律。
在某些特定情况下,作战任务对飞行末速也有着较为严格的要求。我国学者赵汉元对该问题做了一定研究,首先设计理想速度曲线,然后把实际速度控制到理想的速度曲线上去。然而他在推导过程中运用了相应假设并采用了经验公式,实际应用误差较大。因此,需要提出简单易行并可靠性较高的制导方法来解决该领域的问题。
发明内容
本发明为解决多约束条件末制导过程中无法控制末速的问题,提出了一种基于轨迹在线规划的末制导方法。该方法根据俯冲段的终端条件设计参考飞行轨迹,并在飞行过程中根据即时状态不断更新参考轨迹并得到相应的制导控制量。
本发明的技术方案具体如下:
步骤1,建立飞行器的运动学和动力学模型:
其中,x,y,z是地面坐标系下的位置坐标,V是飞行速度,γ,χ分别为弹道倾角和弹道偏角,μ为倾侧角,m是飞行器质量,g是重力加速度,L D分别为升力和阻力,其中, ρ为大气密度,Cx,Cy分别为阻力系数和升力系数,是关于攻角和马赫的函数,Sref为飞行器的参考面积。
由于地面航程s满足:
则得到以下模型:
步骤2,对俯冲段飞行轨迹进行侧向平面和纵向平面设计。
步骤2.1,侧向轨迹设计
本发明采用两段相切圆弧作为侧向参考轨迹,通过转弯达到减速目的。
①初始时刻侧向轨迹设计
初始时刻侧向参考轨迹由两段相切圆弧得到,且第一段圆弧与X轴相切。
两段圆弧的圆心坐标分别为(0,Z1)、(X2,Z2),俯冲段侧向轨迹初、末点坐标分别为(X0,Z0)、(Xf,Zf)。两段圆弧的半径分别为R1,R2,圆心角分别为ψ1,ψ2。ψc为末端的弹道偏角,令其中X0,Z0,Xf,Zf,ψc为设计需满足的终端约束。
由几何关系可得到:
R1=Z1
整理得到:
Z2=Zf+tanψ3(X2-Xf)
从而得出,Z1,Z2,R1,R2,ψ1,ψ2均由第二段圆弧圆心的坐标X2确定,即俯冲段侧向轨迹由X2确定。以俯冲段末速误差最小为性能指标,对X2进行寻优。将得到的最优X2代入上述方程组,得到R1,R2,ψ1,ψ2,(0,Z1)、(X2,Z2),从而确定了初始时刻的侧向轨迹。
②任意时刻侧向轨迹设计
只要两段圆弧的圆心坐标、第一段圆弧的起点坐标和第二段圆弧的终点坐标确定,那么侧向轨迹也就唯一确定。通过旋转平移进行坐标变换,将任意时刻侧向轨迹问题转化为初始时刻侧向轨迹问题。
俯冲段任意时刻t1的弹道偏角为ψc0。将该时刻的侧向弹道曲线整体平移,使t1时刻的俯冲段初点(X0,Z0)与原点重合,然后再将曲线顺时针旋转ψc0,则将t1时刻的侧向轨迹转化为初始侧向剖面。
平移后末点坐标为 旋转矩阵为 平移后的弹道偏角为进而得到然后按照初始侧向轨迹规划方法进行规划。最后将得到的侧向轨迹进行反向旋转和反向平移即可得到真实坐标系下的侧向弹道,假设x,z分别为规划坐标系中的侧向轨迹坐标,X,Z分别为真实坐标系中的侧向轨迹坐标,得到 将弹道偏角通过转化,得到实际坐标系下的弹道偏角:
步骤2.2,地面航程计算
地面剩余航程为两段圆弧的弧长之和。由两段圆弧的半径和圆心角,分别求出两段圆弧的弧长,两段圆弧的航程分别为:
s1=R1ψ1 s2=R2ψ2
则剩余航程为 ssum=s1+s2
步骤2.3,纵向参考轨迹设计
根据俯冲段弹道特性,设计纵向参考轨迹为一段椭圆曲线,横轴为地面航程,纵轴为高度,得到椭圆曲线为:其中,s是实时地面航程;h为高度;a是椭圆长半轴长度;hm为俯冲末点距椭圆长半轴的垂直距离;b是椭圆短半轴长度。
①当期望的弹道倾角末值γf不为-90度时,hm≠0,高度轨迹为:
将高度对地面航程求导,得到弹道倾角参考剖面:
俯冲段飞行过程中,将当前的实时高度hn、弹道倾角γn以及剩余航程ssum和期望的末点弹道倾角γf、末点高度hf=0分别带入高度轨迹和弹道倾角参考剖面,得到:
其中sn,sf分别为当前的航程值和末点航程值,令得到:
将hm,a,b带入高度轨迹方程,得到高度轨迹。
②当期望的弹道倾角末值γf为-90度时,高度轨迹为:
将高度对地面航程求导,得到弹道倾角参考剖面:
俯冲段飞行过程中,将当前的实时高度hn、弹道倾角γn以及剩余航程ssum分别带入高度轨迹和弹道倾角参考剖面,解出a,b:
将a,b带入高度轨迹方程,得到高度轨迹。
步骤3,得到当前实时俯冲段纵向和侧向轨迹之后,计算末制导的控制变量。
具体方法为:将高度轨迹方程对航程s求二阶导数,得到:
对tanγ求时间的一阶导数:
同时,
令
得到
由于侧向运动加速度为向心加速度,
第一段圆弧运动加速度:
第二段圆弧运动加速度:
其中VH=Vcosγ。
令侧向运动加速度为aze,则有其中R=R1或者R2;
从而得到末制导的倾侧角表达式:
升力表达式:
得到升力系数为:
根据升力系数对飞行器气动数据插值,得到末制导的攻角α。
步骤4,将步骤3得到的倾侧角μ和攻角α输入步骤1建立的飞行器模型,对飞行器轨迹进行实时规划调整,使其满足期望的终端条件,从而实现末制导。
有益效果
本发明只需要对一个变量(第二段圆弧圆心坐标)进行寻优,即可得到满足条件的俯冲段飞行轨迹并且得到解析的制导指令,大大减小了计算量,提高了优化效率,缩减了计算时间。在线规划的俯冲段轨迹满足终端约束并且可以控制末速。
附图说明
图1为本发明的俯冲段纵向参考剖面;
图2为本发明的俯冲段初始侧向参考剖面;
图3为本发明的俯冲段任意时刻侧向参考剖面;
图4为本发明的俯冲段轨迹规划流程图;
图5为具体实施方式中一组回归射面的轨迹规划实例;
图6为具体实施方式中规划得到的纵向轨迹;
图7为具体实施方式中规划得到的侧向轨迹;
图8为具体实施方式中轨迹弹道偏角变化曲线;
图9为具体实施方式中轨迹弹道倾角变化曲线;
图10为具体实施方式中攻角变化曲线;
图11为具体实施方式中倾侧角变化曲线;
图12为具体实施方式中速度变化曲线。
具体实施方式
为了更好的说明本发明的目的和优点,下面结合附图和实例进一步详细说明。
选用某型号飞行器,以垂直打击(即弹道倾角末值为-90度)地面目标为例进行介绍。起始位置坐标为(0,30,0)km,起始弹道倾角和弹道偏角均为0度,初始速度为1700m/s。末端位置坐标为(100,0,30)km,末端弹道倾角为-90度,弹道偏角为30度,末速要求为900m/s。考虑到以上终端条件,俯冲段飞行轨迹可以分为在纵向平面和侧向平面分别进行设计。
1.纵向剖面
俯冲段纵向参考剖面如图1所示,其中横轴为地面航程,纵轴为高度,图中的椭圆曲线即为纵向参考轨迹。
在仿真过程中,根据即时状态和期望的终端条件。每个仿真步长更新一次纵向轨迹。
2.侧向剖面
①初始时刻侧向轨迹
俯冲段初始侧向参考剖面如图2所示,是由两段相切的圆弧得到,并且第一段圆弧也与X轴相切。两段圆弧的半径分别为R1,R2,圆心角分别为ψ1,ψ2。只要两段圆弧的圆心坐标、第一段圆弧的起点坐标和第二段圆弧的终点坐标确定,那么侧向轨迹也就唯一确定。
通过几何关系得到,只要第二段圆弧圆心的x坐标X2确定,其他未知量均可确定。因此,俯冲段侧向轨迹可由X2一个变量确定。以俯冲段末速误差最小为性能指标,对X2进行寻优(可以采用matlab自带的fmincon指令),即可得到初始时刻的侧向轨迹。
②任意时刻侧向轨迹求取
在任意的飞行时刻,侧向轨迹规划初始点不在原点且初始弹道偏角不为0。通过旋转平移进行坐标变换,可以将任意时刻轨迹规划问题转化为初始时刻轨迹规划问题。
图3为俯冲段任意时刻(比如t1时刻)的侧向参考剖面,该时刻弹道偏角为ψc0。只需将侧向弹道曲线整体平移,使(X0,Z0)点与原点重合,然后再将曲线顺时针旋转ψc0,即可转化为初始侧向剖面问题。然后按照初始侧向轨迹规划方法进行规划。最后将得到的侧向轨迹进行反向旋转和反向平移即可得到真实坐标系下的侧向弹道。
图4为本发明的流程图。由终端状态分别进行纵向、侧向轨迹规划,使得规划轨迹满足末点和落角的精确,并达到要求的末速,进而得到需要的制导指令。
图5是一组回归射面的俯冲段轨迹规划实例。其中,三维图中的X,Y,Z坐标轴分别为射程轴、高度轴和侧向轴。初点坐标为(0e3,30e3,0e3)m,末点坐标为(100e3,0e3,0e3)m,弹道倾角末值为-90deg,弹道偏角末值的取值从20deg到160deg选取,X2定为90e3m。图中实线为三维飞行轨迹,虚线为飞行轨迹在纵向平面和侧向平面内的投影。从图中可以看出,对于给定的终端条件,本发明提出的规划轨迹算法可以快速准确地得到期望的俯冲段轨迹。由于X2固定,因此该实例没有考虑飞行末速这一指标要求。
图6-图12给出了考虑末速的俯冲段轨迹规划实例,初点坐标为(0e3,30e3,0e3)m,末点坐标为(100e3,0e3,30e3)m,弹道倾角末值为-90deg,弹道偏角末值为30deg,要求末速为900m/s,通过对X2值进行寻优得到满足末速要求的轨迹。图6和图7分别给出了规划轨迹的高度和侧向轨迹,高度轨迹是一条椭圆曲线,侧向是两段圆弧相切轨迹。图8和图9分别给出了弹道偏角和弹道倾角的变化曲线,其中图8中的拐点对应侧向轨迹中两段圆弧相切的点,弹道偏角末值为30度,满足了指标要求;弹道倾角末值为-90度,满足了垂直打击的要求。图10和图11分别给出了攻角和倾侧角的变化曲线,可以看出控制量在两段圆弧相切的位置均有一次“跳变”,从适应侧向轨迹的变化。图12给出了速度的变化曲线,最后末速达到了要求的900m/s。
综上所述,该发明只需确定一个变量,即可准确的规划出一条满足要求的俯冲段轨迹,并得到解析的制导指令,具有很高的工程应用价值。
Claims (2)
1.一种基于轨迹在线规划的末制导方法,其特征在于:包括以下步骤:
步骤1,建立飞行器的运动学和动力学模型:
其中,x,y,z是地面坐标系下的位置坐标,s为地面航程,V是飞行速度,γ,χ分别为弹道倾角和弹道偏角,μ为倾侧角,m是飞行器质量,g是重力加速度,L、D分别为升力和阻力,其中, ρ为大气密度,Cx,Cy分别为阻力系数和升力系数,Sref为飞行器的参考面积;
步骤2,对俯冲段飞行轨迹进行侧向平面和纵向平面设计;
步骤2.1,设计两段相切圆弧作为侧向参考轨迹;
①初始时刻侧向轨迹设计
初始时刻侧向参考轨迹的第一段圆弧与X轴相切;两段圆弧的圆心坐标分别为(0,Z1)、(X2,Z2),俯冲段侧向轨迹初、末点坐标分别为(X0,Z0)、(Xf,Zf);两段圆弧的半径分别为R1,R2,圆心角分别为ψ1,ψ2;ψc为末端的弹道偏角,令其中X0,Z0,Xf,Zf,ψc为设计需满足的终端约束;
由几何关系:
R1=Z1
得到:
Z2=Zf+tanψ3(X2-Xf)
以俯冲段末速误差最小为性能指标,对X2进行寻优;将得到的最优X2代入上述方程组,得到R1,R2,ψ1,ψ2,(0,Z1)、(X2,Z2),从而确定了初始时刻的侧向轨迹;
②俯冲段任意时刻t1的弹道偏角为ψc0;将该时刻的侧向弹道曲线整体平移,使t1时刻的俯冲段初点(X0,Z0)与原点重合,然后再将曲线顺时针旋转ψc0,则将t1时刻的侧向轨迹转化为初始侧向剖面;
平移后末点坐标为 旋转矩阵为 平移后的弹道偏角为得到然后按照初始侧向轨迹规划方法进行规划;最后将得到的侧向轨迹进行反向旋转和反向平移得到真实坐标系下的侧向弹道;规划坐标系中的侧向轨迹坐标为x,z,X,Z分别为真实坐标系中的侧向轨迹坐标,得到 以及实际坐标系下的弹道偏角:
步骤2.2,地面剩余航程为两段圆弧的弧长之和;两段圆弧的航程分别为:
s1=R1ψ1 s2=R2ψ2
则剩余航程为ssum=s1+s2
步骤2.3,纵向参考轨迹设计
根据俯冲段弹道特性,设计纵向参考轨迹为一段椭圆曲线,横轴为地面航程,纵轴为高度,得到椭圆曲线为:其中,s是实时地面航程;h为高度;a是椭圆长半轴长度;hm为俯冲末点距椭圆长半轴的垂直距离;b是椭圆短半轴长度;
①当期望的弹道倾角末值γf不为-90度时,hm≠0,高度轨迹为:
将高度对地面航程求导,得到弹道倾角参考剖面:
俯冲段飞行过程中,将当前的实时高度hn、弹道倾角γn以及剩余航程ssum和期望的末点弹道倾角γf、末点高度hf=0分别带入高度轨迹和弹道倾角参考剖面,得到:
其中sn,sf分别为当前的航程值和末点航程值,令得到:
将hm,a,b带入高度轨迹方程,得到高度轨迹;
②当期望的弹道倾角末值γf为-90度时,高度轨迹为:
将高度对地面航程求导,得到弹道倾角参考剖面:
俯冲段飞行过程中,将当前的实时高度hn、弹道倾角γn以及剩余航程ssum分别带入高度轨迹和弹道倾角参考剖面,解出a,b:
将a,b带入高度轨迹方程,得到高度轨迹;
步骤3,得到当前实时俯冲段纵向和侧向轨迹之后,计算末制导的控制变量;
末制导的倾侧角为:
其中,R=R1或者R2,VH=Vcosγ,
升力系数
根据升力系数对飞行器气动数据插值,得到末制导的攻角α;
步骤4,将步骤3得到的倾侧角μ和攻角α输入步骤1建立的飞行器模型,对飞行器轨迹进行实时规划调整,使其满足期望的终端条件,从而实现末制导。
2.根据权利要求1所述的一种基于轨迹在线规划的末制导方法,其特征在于:所述步骤3中求取倾侧角和升力系数的具体方法为:
将高度轨迹方程对航程s求二阶导数,得到:
对tanγ求时间的一阶导数:
同时,
令
得到
由于侧向运动加速度为向心加速度,
第一段圆弧运动加速度:
第二段圆弧运动加速度:
其中VH=Vcosγ;
令侧向运动加速度为aze,则有其中R=R1或者R2;从而得到末制导的倾侧角表达式:
升力表达式:
得到升力系数为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210472131.0A CN102927851B (zh) | 2012-11-20 | 2012-11-20 | 一种基于轨迹在线规划的末制导方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210472131.0A CN102927851B (zh) | 2012-11-20 | 2012-11-20 | 一种基于轨迹在线规划的末制导方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102927851A CN102927851A (zh) | 2013-02-13 |
CN102927851B true CN102927851B (zh) | 2014-08-20 |
Family
ID=47642712
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210472131.0A Expired - Fee Related CN102927851B (zh) | 2012-11-20 | 2012-11-20 | 一种基于轨迹在线规划的末制导方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102927851B (zh) |
Families Citing this family (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103512426B (zh) * | 2013-09-06 | 2015-05-06 | 北京理工大学 | 一种次优的带末角约束制导方法 |
CN103587723B (zh) * | 2013-11-07 | 2015-12-02 | 北京临近空间飞行器***工程研究所 | 一种再入初始段解析式纵向在线轨迹设计及跟踪方法 |
CN104156595B (zh) * | 2014-08-11 | 2015-06-24 | 北京航天自动控制研究所 | 飞行器的飞行轨迹指令的确定方法和装置 |
CN104298246B (zh) * | 2014-09-05 | 2017-02-08 | 湖北航天技术研究院总体设计所 | 一种盘旋飞行轨迹的导引律设计方法及*** |
CN104729504B (zh) * | 2015-03-09 | 2015-12-02 | 北京航天自动控制研究所 | 一种无动力高超声速飞行器的倾斜角的确定方法 |
CN104965519B (zh) * | 2015-06-10 | 2018-08-10 | 北京理工大学 | 一种基于贝塞尔曲线的带落角约束的末制导方法 |
CN105759830B (zh) * | 2016-04-29 | 2017-03-22 | 北京航天自动控制研究所 | 一种升力式飞行器高动态下压段制导方法 |
CN106294280B (zh) * | 2016-08-22 | 2018-11-20 | 麻毅威 | 一种弹道规划方法 |
CN107238389A (zh) * | 2017-06-27 | 2017-10-10 | 中国航空工业集团公司沈阳飞机设计研究所 | 飞机航线规划方法 |
CN107621198B (zh) * | 2017-08-28 | 2019-04-12 | 北京航空航天大学 | 一种绕多禁飞区的多步决策弹道规划方法 |
CN109253730B (zh) * | 2018-11-07 | 2021-03-23 | 广州大学 | 可重复使用运载器末端能量管理段的三维轨迹在线规划方法及*** |
CN110457647B (zh) * | 2019-07-29 | 2020-09-18 | 上海机电工程研究所 | 一种旋转导弹弹目遭遇时间估算方法 |
CN111366044B (zh) * | 2019-12-29 | 2022-10-04 | 湖北航天飞行器研究所 | 一种平飞过渡段制导控制方法 |
CN111348223B (zh) * | 2020-05-25 | 2020-08-21 | 北京星际荣耀空间科技有限公司 | 一种控制弹道顶点高度的闭路制导方法、装置及设备 |
CN112525003B (zh) * | 2020-11-30 | 2022-12-27 | 北京宇航***工程研究所 | 一种带落角约束的扩展比例导引方法 |
CN112904888B (zh) * | 2021-01-11 | 2024-04-09 | 北京临近空间飞行器***工程研究所 | 多目标参数联合制导的方法 |
CN113051743B (zh) * | 2021-03-18 | 2023-05-26 | 中国人民解放***箭军工程大学 | 一种基于轨迹在线规划的末制导*** |
CN113064448B (zh) * | 2021-03-26 | 2022-04-05 | 北京理工大学 | 一种飞行器的制导方法及*** |
CN113176563B (zh) * | 2021-04-12 | 2023-02-28 | 中国人民解放军96901部队26分队 | 一种雷达景象匹配末制导导弹飞行弹道规划方法 |
CN113759966B (zh) * | 2021-08-26 | 2023-10-20 | 北京理工大学 | 一种三维空间内终端速度可控的末制导方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1361408A (zh) * | 2000-12-23 | 2002-07-31 | 林清芳 | 实时惯性测量装置模拟器 |
CN102645933A (zh) * | 2012-05-02 | 2012-08-22 | 中国人民解放军海军航空工程学院 | 飞行器实现大空域机动的组合过载控制方法 |
CN102706217A (zh) * | 2012-04-17 | 2012-10-03 | 北京理工大学 | 一种控制多枚导弹攻击角度和攻击时间的方法 |
-
2012
- 2012-11-20 CN CN201210472131.0A patent/CN102927851B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1361408A (zh) * | 2000-12-23 | 2002-07-31 | 林清芳 | 实时惯性测量装置模拟器 |
CN102706217A (zh) * | 2012-04-17 | 2012-10-03 | 北京理工大学 | 一种控制多枚导弹攻击角度和攻击时间的方法 |
CN102645933A (zh) * | 2012-05-02 | 2012-08-22 | 中国人民解放军海军航空工程学院 | 飞行器实现大空域机动的组合过载控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102927851A (zh) | 2013-02-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102927851B (zh) | 一种基于轨迹在线规划的末制导方法 | |
CN108168381B (zh) | 一种多枚导弹协同作战的控制方法 | |
US11286065B2 (en) | Method for designing reentry trajectory based on flight path angle planning | |
CN102980449B (zh) | 一种多枚导弹协同作战的控制方法 | |
CN103245257B (zh) | 基于Bezier曲线的多约束飞行器导引方法 | |
CN106371312B (zh) | 基于模糊控制器的升力式再入预测-校正制导方法 | |
CN103995540A (zh) | 一种高超声速飞行器的有限时间轨迹快速生成方法 | |
CN103557867A (zh) | 一种基于稀疏a*搜索的三维多uav协同航迹规划方法 | |
CN105652879A (zh) | 一种无副翼无人机自主飞行控制方法 | |
CN111580547B (zh) | 一种高超声速飞行器编队控制方法 | |
CN104019701B (zh) | 一种利用直接力气动力复合控制的前向拦截制导方法 | |
CN104317305A (zh) | 一种面向复杂战场威胁的射前航迹确定方法 | |
CN105425812A (zh) | 一种基于双模型下的无人机自动着舰轨迹控制方法 | |
CN111027143A (zh) | 一种基于深度强化学习的舰载机进近引导方法 | |
CN103983143A (zh) | 包含速度过程约束和多终端约束的空地导弹投放下滑段制导方法 | |
CN104503471A (zh) | 一种机动飞行器多终端约束反演滑模末制导方法 | |
CN111290278A (zh) | 一种基于预测滑模的高超声速飞行器鲁棒姿态控制方法 | |
CN115857538A (zh) | 三维空间下满足落角约束的多飞行器协同制导方法 | |
CN104656659B (zh) | 舰载机滑跃起飞自动飞行控制方法 | |
CN113741509B (zh) | 一种高超声速滑翔飞行器下压段能量管理方法 | |
You et al. | Distributed synergetic guidance law for multiple missiles with angle-of-attack constraint | |
CN113111433B (zh) | 一种双线程嵌入式实时轨迹优化与制导方法 | |
CN102700727B (zh) | 一种基于速度控制的对空拦截飞行器制导方法 | |
CN106292700B (zh) | 一种大落地倾角条件下应用的侧向导引方法 | |
CN115454145A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140820 Termination date: 20191120 |
|
CF01 | Termination of patent right due to non-payment of annual fee |