CN111442697A - 一种基于伪谱法修正的过重补制导方法和弹道整形制导方法 - Google Patents
一种基于伪谱法修正的过重补制导方法和弹道整形制导方法 Download PDFInfo
- Publication number
- CN111442697A CN111442697A CN202010082200.1A CN202010082200A CN111442697A CN 111442697 A CN111442697 A CN 111442697A CN 202010082200 A CN202010082200 A CN 202010082200A CN 111442697 A CN111442697 A CN 111442697A
- Authority
- CN
- China
- Prior art keywords
- guidance
- angle
- missile
- pseudo
- law
- 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.)
- Pending
Links
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F42—AMMUNITION; BLASTING
- F42B—EXPLOSIVE CHARGES, e.g. FOR BLASTING, FIREWORKS, AMMUNITION
- F42B15/00—Self-propelled projectiles or missiles, e.g. rockets; Guided missiles
- F42B15/01—Arrangements thereon for guidance or control
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F41—WEAPONS
- F41G—WEAPON SIGHTS; AIMING
- F41G3/00—Aiming or laying means
-
- 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
- G06F17/13—Differential equations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Software Systems (AREA)
- Operations Research (AREA)
- Chemical & Material Sciences (AREA)
- Aviation & Aerospace Engineering (AREA)
- Combustion & Propulsion (AREA)
- Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)
Abstract
本发明提供了一种基于伪谱弹道优化结果修正的过重补制导方法和弹道整形制导方法,属于航天技术、武器技术、制导控制领域。具体包括:一、建立交战模型,确定导弹气动参数和大气模型参数;二、建立导弹运动方程组并建立最优控制问题;三、使用伪谱法对最优控制问题进行求解;四、通过多项式拟合等方式拟合伪谱最优结果;五、根据伪谱最优结果修正典型弹道整形制导律和典型重力补偿制导律。优点在于:提升导弹的末速,大大增强导弹的攻击效能;通过对伪谱最优结果的拟合,将制导律分段,前一段提升导弹飞行高度,后一段调整制导参数使其满足落角约束;解决了典型重力补偿制导律的“黑洞”现象,对实际应用有重要意义。
Description
技术领域
本发明提供了一种基于伪谱弹道优化结果修正的过重补制导方法和弹道整形制导方法,属于航天技术、武器技术、制导控制领域。
背景技术
双脉冲中远程空空弹的制导问题需要解决的是能量管理问题,需要合理控制导弹的速度曲线,同时提高导弹飞行高度,使导弹飞行的大部分时间都处于空气稀薄、密度较小的高空区域。为了实现对高价值高风险目标的打击,现代精确制导飞行器需要进一步提升作战能力,因此工程中对其提出了更高的要求,比如限制导弹落角。因此,有必要研究可提高导弹弹道高度同时通过一定手段可调节导弹落角的制导律。
现今,常见的可提升导弹弹道高度的典型制导律包括:弹道整形制导律和过重力补偿制导律等,均属于增广比例导引法。由于对比离线伪谱最优仿真结果以及典型制导律仿真结果发现,在满足了一定的落角约束后,使用典型制导律进行制导的导弹无法达到较高的高度,且重力补偿制导律在一定射程范围内会出现无法击中目标的现象。因此,需要对典型制导律进行改进,改进典型制导律可提升弹道高度,同时通过调整制导参数控制导弹落角。
发明内容
针对中远程空空导弹的上述问题,本发明基于现有典型制导律,根据伪谱法离线解算的最优结果,提出了一种基于伪谱法修正的过重补制导方法和弹道整形制导方法,该制导方法在导弹发射后的短时间内拟合最优结果,从而提升弹道高度,且导弹落角可以通过调整制导参数进行调整,大幅提升典型制导律的制导能力。
具体步骤如下:
步骤一:建立空空弹运动学模型
在空空导弹的中制导段,可以将其运动分解在侧向和纵向两个平面内。对于侧向平面,由于同一高度的大气参数相同,且不存在重力在水平方向的分量,横向运动只取决于导弹目标的终端位置,通常采用比例导引法,它的飞行弹道是平直的,运动幅度不大。主要考虑纵向平面内的高度变化,因此将导弹的运动限制在铅锤面内。系统的运动学方程如式(1)所示,
式中,r为水平位置,h为高度,v为速度,γ为弹道倾角,T为推力,α为攻角,L和D为导弹受到的升力和阻力。对于双脉冲空空导弹来说,导弹的质量流量和推力可写成随时间变化的分段函数,
步骤二:建立最优控制问题并使用伪谱法进行解算
使用伪谱法可获得二次点火问题的一系列最优弹道,此处的弹道优化问题由于推力的多脉冲形式,是一个多段的最优控制问题。选择状态量写为式(4)。
X=[x1 x2 x3 x4 x5 x6]T=[r h V γ m α]T (4)
根据导弹的动力学建模,写出微分方程为方程(5)。
由于该优化问题是一个多段最优控制问题,需要约束状态量在不同段之间的连接条件,并且需对控制变量加以约束。为保证飞行品质和可操作性,还需增加路径约束为舵偏角约束。
目标函数的选择对弹道起着决定性的影响,通常考虑最短飞行时间或是最大末速度,这两者互相矛盾需要取舍,故设计一个加权系数k来观察不同的目标函数的效果,
minJ=k*tf-v(tf) (6)
k的取值决定了目标函数对两个因素的侧重:k=0,对应最大末速度指标;0<k<1,更侧重末速度指标,而且由于末速度的变化范围更大,所以此时基本等同于最大末速度情况;k>>1,更侧重飞行时间指标;k→∞,对应最短飞行时间指标。
步骤三:拟合最优制导规律
改变弹目初始状态,使用伪谱法解算最优制导律。观察攻角随时间变化的曲线图,发现在发射后的10s之内,攻角呈现线性变化规律,初始攻角随射程的变化而变化,根据此规律对典型制导律进行改进。
将典型制导律进行分段,前10s内按照攻角线性变化的规律推导此时控制量过载的值,得到改进弹道整形制导律的公式(7)如下,
改进重力补偿制导律的公式(8)如下,
其中,CL为引力系数,α0为初始攻角,k为攻角随时间变化律,t为时间,m为质量,g为重力加速度,T为当前推力,Vc为弹目相对速度,γ为弹道倾角,λ为视线角,λf为预设末端视线角,tgo为剩余飞行时间,N'为比例导引系数,c为重力补偿系数。
对改进弹道整形制导律来说,通过调整预设末端视线角λf控制导弹落角大小;对改进重力补偿制导律来说,通过调整重力补偿系数c控制导弹落角大小。
步骤四:选择典型制导律的制导参数
改进弹道整形制导律和改进重力补偿制导律中初始攻角α0和攻角随时间变化律k可根据不同弹目初始状态下的伪谱最优结果进行调整,使改进典型制导律的仿真结果尽可能接近伪谱最优结果,如式(9)。
在对落角约束较宽松的条件下,随着初始条件的变化,仿真结果显示最优落角并不一定是最大落角,同理可根据初始条件的变化对落角的变化进行拟合,如式(10)所示,其中γf,min为根据实际状况要求的落角约束,如约束为-50°,γf,op为无落角约束时的伪谱最优结果,γf则为最终选取的落角,可对制导参数进行调整使得导弹落角达到γf。
γf=max(γf,op,γf,min)
在实际应用中,经过离线的仿真拟合,可计算出制导参数γf和c与初始条件和当前状态的关系式,从而实现制导律的在线应用。
至此,本发明提出的控制规律已经可以仿真应用。
下面通过一个仿真的伪代码说明控制律的仿真方法:
1.获取包括弹目速度、位置和姿态等在内的初始状态信息。
2.根据弹目初始状态计算初始攻角α0和攻角随时间变化律k。
3.根据当前状态量推算出制导参数γf和c的值,并判断当前飞行时间是否大于10s,根据飞时间按照不同规律计算控制量。
4.用四阶龙格库塔法对给定的控制量进行动力学仿真,判断是否击中或错过目标,若击中或错过,程序运行结束,未击中或错过,继续。
5.回到步骤3,直到击中或错过目标。
综上所述,通过上述四个步骤,推导出了本发明方法,即一种基于伪谱法修正的过重补制导方法和弹道整形制导方法,算例仿真结果表明本发明方法能够通过拟合伪谱最优结果提高弹道高度,调整制导参数则可约束导弹落角,大大提升了导弹的射程,且消除了典型重力补偿制导律的“黑洞”现象,综合性能优异。
本发明的优点在于:
(1)提出了一种基于伪谱法修正的过重补制导方法和弹道整形制导方法,能够提升导弹的末速,大大增强导弹的攻击效能。
(2)通过对伪谱最优结果的拟合,将制导律分段,前一段提升导弹飞行高度,后一段调整制导参数使其满足落角约束。
(3)解决了典型重力补偿制导律的“黑洞”现象,对实际应用有重要意义。
附图说明
图1是摘要附图。
图2是空空弹在纵平面内的受力分析示意图。
图3是射程不同时的伪谱最优弹道曲线。
图4是射程不同时的伪谱最优速度曲线。
图5是射程不同时的伪谱最优弹道倾角曲线。
图6是射程不同时的伪谱最优攻角曲线。
图7是射程不同时的伪谱最优动压曲线。
图8是铅锤面内的弹目交战示意图。
图9是不同制导律导弹末速对比。
图10是不同制导律飞行时间对比。
图11是不同制导律最大高度对比。
图12是不同制导律弹道曲线对比。
图13是不同制导律速度曲线对比。
图14是不同制导律弹道倾角曲线对比。
图15是不同制导律攻角曲线对比。
图16是不同制导律动压曲线对比。
图17是黑洞现象弹道曲线图。
图18是黑洞现象速度变化图
图19是基于伪谱最优结果可提高弹道高度的改进典型制导方法仿真流程图。
上述图中,涉及到的符号、代号说明如下:
图2中,L、D、G和T分别为导弹所受升力、阻力、重力和推力,V为导弹速度向量,α、γ和ψ分别为攻角、弹道倾角和俯仰角,下同。图9中,λ为视线角。图18中,CL为引力系数,α0为初始攻角,k为攻角随时间变化律,t为时间,m为质量,g为重力加速度,Vc为弹目相对速度,λf为预设末端视线角,tgo为剩余飞行时间,N'为比例导引系数,c为重力补偿系数。
具体实施方式
下面将结合附图和实施案例对本发明作进一步的详细说明。
针对中远程空空导弹的上述问题,本发明基于现有典型制导律,根据伪谱法离线解算的最优结果,提出了一种基于伪谱法修正的过重补制导方法和弹道整形制导方法,该制导方法在导弹发射后的短时间内拟合最优结果,从而提升弹道高度,且导弹落角可以通过调整制导参数进行调整,大幅提升典型制导律的制导能力。
具体步骤如下:
步骤一:建立空空弹运动学模型
在空空导弹的中制导段,可以将其运动分解在侧向和纵向两个平面内。对于侧向平面,由于同一高度的大气参数相同,且不存在重力在水平方向的分量,横向运动只取决于导弹目标的终端位置,通常采用比例导引法,它的飞行弹道是平直的,运动幅度不大。主要考虑纵向平面内的高度变化,因此将导弹的运动限制在铅锤面内,铅锤面内导弹受力分析示意图如图2所示。系统的运动学方程如式(1)所示,
式中,r为水平位置,h为高度,v为速度,γ为弹道倾角,T为推力,α为攻角,L和D为导弹受到的升力和阻力。对于双脉冲空空导弹来说,导弹的质量流量和推力可写成随时间变化的分段函数,
步骤二:建立最优控制问题并使用伪谱法进行解算
使用伪谱法可获得二次点火问题的一系列最优弹道,此处的弹道优化问题由于推力的多脉冲形式,是一个多段的最优控制问题。选择状态量写为式(4)。
X=[x1 x2 x3 x4 x5 x6]T=[r h V γ m α]T (4)
根据导弹的动力学建模,写出微分方程为方程(5)。
由于该优化问题是一个多段最优控制问题,需要约束状态量在不同段之间的连接条件,并且需对控制变量加以约束。为保证飞行品质和可操作性,还需增加路径约束为舵偏角约束。
目标函数的选择对弹道起着决定性的影响,通常考虑最短飞行时间或是最大末速度,这两者互相矛盾需要取舍,故设计一个加权系数k来观察不同的目标函数的效果,
minJ=k*tf-v(tf) (6)
k的取值决定了目标函数对两个因素的侧重:k=0,对应最大末速度指标;0<k<1,更侧重末速度指标,而且由于末速度的变化范围更大,所以此时基本等同于最大末速度情况;k>>1,更侧重飞行时间指标;k→∞,对应最短飞行时间指标。
步骤三:拟合最优制导规律
使用伪谱法对最优制导问题进行解算,改变目标初始位置,使用伪谱法解算最优制导律,仿真弹道曲线、速度曲线、弹道倾角曲线、攻角曲线以及动压曲线如图3-图7所示。观察攻角随时间变化的曲线图,发现在发射后的10s之内,攻角近似呈现线性变化规律,初始攻角随射程的变化而变化,根据此规律对典型制导律进行改进。
将典型制导律进行分段,前10s内按照攻角线性变化的规律推导此时控制量过载的值,得到改进弹道整形制导律的公式(7)如下,
改进重力补偿制导律的公式(8)如下,
其中,CL为引力系数,α0为初始攻角,k为攻角随时间变化律,t为时间,m为质量,g为重力加速度,T为当前推力,Vc为弹目相对速度,γ为弹道倾角,λ为视线角,如图8所示,λf为预设末端视线角,tgo为剩余飞行时间,N'为比例导引系数,c为重力补偿系数。
对改进弹道整形制导律来说,通过调整预设末端视线角λf控制导弹落角大小;对改进重力补偿制导律来说,通过调整重力补偿系数c控制导弹落角大小。
步骤四:选择典型制导律的制导参数
改进弹道整形制导律和改进重力补偿制导律中初始攻角α0和攻角随时间变化律k可根据不同弹目初始状态下的伪谱最优结果进行调整,使改进典型制导律的仿真结果尽可能接近伪谱最优结果,如式(9)。
在对落角约束较宽松的条件下,随着初始条件的变化,仿真结果显示最优落角并不一定是最大落角,同理可根据初始条件的变化对落角的变化进行拟合,如式(10)所示,其中γf,min为根据实际状况要求的落角约束,如约束为-50°,γf,op为无落角约束时的伪谱最优结果,γf则为最终选取的落角,可对制导参数进行调整使得导弹落角达到γf。
γf=max(γf,op,γf,min)
本算例中,由于初始条件仅有目标的初始位置rt0有变化,故α0和γf可拟合为随rt0变化的多项式,k取-1.8,如式(13),式中rt0单位为千米[km],α0单位为度[°],
仿真对比伪谱最优结果、改进弹道整形制导律、改进重力补偿制导律、弹道整形制导律和重力补偿制导律的制导效果。其中制导参数对改进弹道整形制导律和典型弹道整形制导律来说为预设末端视线角λf,对改进重力补偿制导律和典型重力补偿制导律来说为重力补偿系数c。
对比不同制导律的仿真结果,如图9-图11分别为导弹末速、飞行时间和最大高度的对比图,其中脱靶量过大的点未绘制,不具有比较意义。攻击距离为150km的弹道对比曲线、速度对比曲线、弹道倾角对比曲线、攻角对比曲线和动压对比曲线图如图12-图16所示,改进典型制导律的制导效果相比典型制导律有很大提升,弹道更接近伪谱最优结果。改进典型制导律的最大高度接近伪谱最优结果,飞行时间略长,导弹末速稍小。
对重力补偿制导律来说,改进重力补偿制导律消除了典型重力补偿制导律的“黑洞”现象,“黑洞”是只导弹在射程较小和较大的情况下均可击中目标,中间射程则由于落速过小无法击中目标的现象。如图17-图18所示。随着射程的增加,为使导弹击中目标,典型重力补偿制导律需要通过增大重补系数来提高弹道高度,弹道高度越高,由于空气密度减小故空气阻力减小,但同时航程增加,故能量的损失与弹道高度之间不是线性关系。而改进重力补偿制导律则由于发射后拟合了伪谱最优结果,可使导弹在较优的空域飞行,从而消除了“黑洞”现象。
在实际应用中,经过离线的仿真拟合,可计算出制导参数γf和c与初始条件和当前状态的关系式,从而实现制导律的在线应用。
至此,本发明提出的控制规律已经可以仿真应用。下面通过一个仿真的伪代码说明控制律的仿真方法,流程图如图19所示:
1.获取包括弹目速度、位置和姿态等在内的初始状态信息。
2.根据弹目初始状态计算初始攻角α0和攻角随时间变化律k。
3.根据当前状态量推算出制导参数γf和c的值,并判断当前飞行时间是否大于10s,根据飞时间按照不同规律计算控制量。
4.用四阶龙格库塔法对给定的控制量进行动力学仿真,判断是否击中或错过目标,若击中或错过,程序运行结束,未击中或错过,继续。
5.回到步骤3,直到击中或错过目标。
综上所述,通过上述四个步骤,推导出了本发明方法,即一种基于伪谱法修正的过重补制导方法和弹道整形制导方法,算例仿真结果表明本发明方法能够通过拟合伪谱最优结果提高弹道高度,调整制导参数则可约束导弹落角,大大提升了导弹的射程,且消除了典型重力补偿制导律的“黑洞”现象,综合性能优异。
Claims (1)
1.一种基于伪谱法修正的过重补制导方法和弹道整形制导方法,其特征在于:该制导方法在导弹发射后的短时间内拟合最优结果,从而提升弹道高度,且导弹落角可以通过调整制导参数进行调整,大幅提升典型制导律的制导能力。具体步骤如下:
步骤一:建立空空弹运动学模型
在空空导弹的中制导段,可以将其运动分解在侧向和纵向两个平面内。对于侧向平面,由于同一高度的大气参数相同,且不存在重力在水平方向的分量,横向运动只取决于导弹目标的终端位置,通常采用比例导引法,它的飞行弹道是平直的,运动幅度不大。主要考虑纵向平面内的高度变化,因此将导弹的运动限制在铅锤面内。系统的运动学方程如式(1)所示,
式中,r为水平位置,h为高度,v为速度,γ为弹道倾角,T为推力,α为攻角,L和D为导弹受到的升力和阻力。对于双脉冲空空导弹来说,导弹的质量流量和推力可写成随时间变化的分段函数,
步骤二:建立最优控制问题并使用伪谱法进行解算
使用伪谱法可获得二次点火问题的一系列最优弹道,此处的弹道优化问题由于推力的多脉冲形式,是一个多段的最优控制问题。选择状态量写为式(4)。
X=[x1 x2 x3 x4 x5 x6]T=[r h V γ m α]T (4)
根据导弹的动力学建模,写出微分方程为方程(5)。
由于该优化问题是一个多段最优控制问题,需要约束状态量在不同段之间的连接条件,并且需对控制变量加以约束。为保证飞行品质和可操作性,还需增加路径约束为舵偏角约束。
目标函数的选择对弹道起着决定性的影响,通常考虑最短飞行时间或是最大末速度,这两者互相矛盾需要取舍,故设计一个加权系数k来观察不同的目标函数的效果,
min J=k*tf-v(tf) (6)
k的取值决定了目标函数对两个因素的侧重:k=0,对应最大末速度指标;0<k<1,更侧重末速度指标,而且由于末速度的变化范围更大,所以此时基本等同于最大末速度情况;k>>1,更侧重飞行时间指标;k→∞,对应最短飞行时间指标。
步骤三:拟合最优制导规律
改变弹目初始状态,使用伪谱法解算最优制导律。观察攻角随时间变化的曲线图,发现在发射后的10s之内,攻角呈现线性变化规律,初始攻角随射程的变化而变化,根据此规律对典型制导律进行改进。
将典型制导律进行分段,前10s内按照攻角线性变化的规律推导此时控制量过载的值,得到改进弹道整形制导律的公式(7)如下,
改进重力补偿制导律的公式(8)如下,
其中,CL为引力系数,α0为初始攻角,k为攻角随时间变化律,t为时间,m为质量,g为重力加速度,T为当前推力,Vc为弹目相对速度,γ为弹道倾角,λ为视线角,如图8所示,λf为预设末端视线角,tgo为剩余飞行时间,N'为比例导引系数,c为重力补偿系数。
对改进弹道整形制导律来说,通过调整预设末端视线角λf控制导弹落角大小;对改进重力补偿制导律来说,通过调整重力补偿系数c控制导弹落角大小。
步骤四:选择典型制导律的制导参数
改进弹道整形制导律和改进重力补偿制导律中初始攻角α0和攻角随时间变化律k可根据不同弹目初始状态下的伪谱最优结果进行调整,使改进典型制导律的仿真结果尽可能接近伪谱最优结果,如式(9)。
在对落角约束较宽松的条件下,随着初始条件的变化,仿真结果显示最优落角并不一定是最大落角,同理可根据初始条件的变化对落角的变化进行拟合,如式(10)所示,其中γf,min为根据实际状况要求的落角约束,如约束为-50°,γf,op为无落角约束时的伪谱最优结果,γf则为最终选取的落角,可对制导参数进行调整使得导弹落角达到γf。
γf=max(γf,op,γf,min)
在实际应用中,经过离线的仿真拟合,可计算出制导参数γf和c与初始条件和当前状态的关系式,从而实现制导律的在线应用。
至此,本发明提出的控制规律已经可以仿真应用。
下面通过一个仿真的伪代码说明控制律的仿真方法:
a)获取包括弹目速度、位置和姿态等在内的初始状态信息。
b)根据弹目初始状态计算初始攻角α0和攻角随时间变化律k。
c)根据当前状态量推算出制导参数γf和c的值,并判断当前飞行时间是否大于10s,根据飞时间按照不同规律计算控制量。
d)用四阶龙格库塔法对给定的控制量进行动力学仿真,判断是否击中或错过目标,若击中或错过,程序运行结束,未击中或错过,继续。
e)回到步骤c),直到击中或错过目标。
综上所述,通过上述四个步骤,推导出了本发明方法,即一种基于伪谱法修正的过重补制导方法和弹道整形制导方法,算例仿真结果表明本发明方法能够通过拟合伪谱最优结果提高弹道高度,调整制导参数则可约束导弹落角,大大提升了导弹的射程,且消除了典型重力补偿制导律的“黑洞”现象,综合性能优异。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010082200.1A CN111442697A (zh) | 2020-02-07 | 2020-02-07 | 一种基于伪谱法修正的过重补制导方法和弹道整形制导方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010082200.1A CN111442697A (zh) | 2020-02-07 | 2020-02-07 | 一种基于伪谱法修正的过重补制导方法和弹道整形制导方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111442697A true CN111442697A (zh) | 2020-07-24 |
Family
ID=71627104
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010082200.1A Pending CN111442697A (zh) | 2020-02-07 | 2020-02-07 | 一种基于伪谱法修正的过重补制导方法和弹道整形制导方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111442697A (zh) |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113311871A (zh) * | 2021-06-02 | 2021-08-27 | 上海机电工程研究所 | 跳跃-滑翔导弹的制导律优化方法及系统 |
CN113587740A (zh) * | 2021-07-22 | 2021-11-02 | 北京航天飞腾装备技术有限责任公司 | 一种基于弹目视线角的被动反辐射导引方法及系统 |
CN113739635A (zh) * | 2021-08-26 | 2021-12-03 | 北京航天飞腾装备技术有限责任公司 | 一种实现导弹大扇面角发射的制导方法 |
CN113742643A (zh) * | 2021-07-29 | 2021-12-03 | 北京航空航天大学 | 一种基于系统状态约束的低能耗弹道优化方法 |
CN114020018A (zh) * | 2021-11-03 | 2022-02-08 | 北京航空航天大学 | 导弹控制策略的确定方法、装置、存储介质及电子设备 |
CN114115325A (zh) * | 2020-08-31 | 2022-03-01 | 南京理工大学 | 一种基于hp-RPM算法的滑翔弹在线闭环制导方法 |
CN114440707A (zh) * | 2022-02-09 | 2022-05-06 | 北京信息科技大学 | 顶部与侧面协同拦截的飞行器制导方法、装置及系统 |
CN114646238A (zh) * | 2022-03-30 | 2022-06-21 | 南京理工大学 | 飞行体状态感知自适应方案弹道跟踪方法 |
CN114719689A (zh) * | 2021-05-21 | 2022-07-08 | 北京航空航天大学 | 基于弹道解析解的运载火箭助推段制导方法及装置 |
CN114935277A (zh) * | 2022-03-05 | 2022-08-23 | 南京理工大学 | 一种滑翔增程制导炮弹理想弹道的在线规划方法 |
CN115047767A (zh) * | 2022-06-14 | 2022-09-13 | 北京中科飞鸿科技股份有限公司 | 一种考虑高度表使用约束的俯冲弹道优化方法 |
CN115079565A (zh) * | 2022-08-23 | 2022-09-20 | 中国人民解放军国防科技大学 | 变系数的带落角约束制导方法、装置和飞行器 |
CN115309059A (zh) * | 2022-10-10 | 2022-11-08 | 北京航天自动控制研究所 | 一种考虑引力补偿的直接制导方法 |
CN115857325A (zh) * | 2022-10-09 | 2023-03-28 | 北京航天自动控制研究所 | 一种基于多项式弹道的能量管理制导方法、存储介质 |
CN116379847A (zh) * | 2023-06-05 | 2023-07-04 | 西安现代控制技术研究所 | 一种导弹在线能量管理制导律构建方法 |
CN117073473A (zh) * | 2023-10-17 | 2023-11-17 | 中国空气动力研究与发展中心空天技术研究所 | 一种基于时间约束的导弹视角规划制导方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1767893A1 (en) * | 2005-09-23 | 2007-03-28 | Saab Ab | Missile guidance system |
CN104217041A (zh) * | 2013-05-29 | 2014-12-17 | 北京航空航天大学 | 一种多约束在线高斯伪谱末制导方法 |
CN105841550A (zh) * | 2016-04-15 | 2016-08-10 | 哈尔滨工业大学 | 一种具有高度约束的高置修正比例导引律方法 |
CN108107899A (zh) * | 2018-01-19 | 2018-06-01 | 北京航空航天大学 | 线性伪谱外太空飞行器能量最优大角度姿态机动控制方法 |
CN109270960A (zh) * | 2018-12-05 | 2019-01-25 | 中南大学 | 基于Radau伪谱法的在线最优反馈再入制导方法 |
-
2020
- 2020-02-07 CN CN202010082200.1A patent/CN111442697A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1767893A1 (en) * | 2005-09-23 | 2007-03-28 | Saab Ab | Missile guidance system |
CN104217041A (zh) * | 2013-05-29 | 2014-12-17 | 北京航空航天大学 | 一种多约束在线高斯伪谱末制导方法 |
CN105841550A (zh) * | 2016-04-15 | 2016-08-10 | 哈尔滨工业大学 | 一种具有高度约束的高置修正比例导引律方法 |
CN108107899A (zh) * | 2018-01-19 | 2018-06-01 | 北京航空航天大学 | 线性伪谱外太空飞行器能量最优大角度姿态机动控制方法 |
CN109270960A (zh) * | 2018-12-05 | 2019-01-25 | 中南大学 | 基于Radau伪谱法的在线最优反馈再入制导方法 |
Cited By (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114115325A (zh) * | 2020-08-31 | 2022-03-01 | 南京理工大学 | 一种基于hp-RPM算法的滑翔弹在线闭环制导方法 |
CN114115325B (zh) * | 2020-08-31 | 2023-02-28 | 南京理工大学 | 一种基于hp-RPM算法的滑翔弹在线闭环制导方法 |
CN114719689B (zh) * | 2021-05-21 | 2023-07-28 | 北京航空航天大学 | 基于弹道解析解的运载火箭助推段制导方法及装置 |
CN114719689A (zh) * | 2021-05-21 | 2022-07-08 | 北京航空航天大学 | 基于弹道解析解的运载火箭助推段制导方法及装置 |
CN113311871B (zh) * | 2021-06-02 | 2022-12-13 | 上海机电工程研究所 | 跳跃-滑翔导弹的制导律优化方法及系统 |
CN113311871A (zh) * | 2021-06-02 | 2021-08-27 | 上海机电工程研究所 | 跳跃-滑翔导弹的制导律优化方法及系统 |
CN113587740A (zh) * | 2021-07-22 | 2021-11-02 | 北京航天飞腾装备技术有限责任公司 | 一种基于弹目视线角的被动反辐射导引方法及系统 |
CN113742643A (zh) * | 2021-07-29 | 2021-12-03 | 北京航空航天大学 | 一种基于系统状态约束的低能耗弹道优化方法 |
CN113742643B (zh) * | 2021-07-29 | 2023-06-23 | 北京航空航天大学 | 一种基于系统状态约束的低能耗弹道优化方法 |
CN113739635B (zh) * | 2021-08-26 | 2023-01-24 | 北京航天飞腾装备技术有限责任公司 | 一种实现导弹大扇面角发射的制导方法 |
CN113739635A (zh) * | 2021-08-26 | 2021-12-03 | 北京航天飞腾装备技术有限责任公司 | 一种实现导弹大扇面角发射的制导方法 |
CN114020018A (zh) * | 2021-11-03 | 2022-02-08 | 北京航空航天大学 | 导弹控制策略的确定方法、装置、存储介质及电子设备 |
CN114020018B (zh) * | 2021-11-03 | 2024-02-27 | 北京航空航天大学 | 导弹控制策略的确定方法、装置、存储介质及电子设备 |
CN114440707B (zh) * | 2022-02-09 | 2023-10-20 | 北京信息科技大学 | 顶部与侧面协同拦截的飞行器制导方法、装置及系统 |
CN114440707A (zh) * | 2022-02-09 | 2022-05-06 | 北京信息科技大学 | 顶部与侧面协同拦截的飞行器制导方法、装置及系统 |
CN114935277A (zh) * | 2022-03-05 | 2022-08-23 | 南京理工大学 | 一种滑翔增程制导炮弹理想弹道的在线规划方法 |
CN114935277B (zh) * | 2022-03-05 | 2023-08-04 | 南京理工大学 | 一种滑翔增程制导炮弹理想弹道的在线规划方法 |
CN114646238A (zh) * | 2022-03-30 | 2022-06-21 | 南京理工大学 | 飞行体状态感知自适应方案弹道跟踪方法 |
CN114646238B (zh) * | 2022-03-30 | 2023-08-04 | 南京理工大学 | 飞行体状态感知自适应方案弹道跟踪方法 |
CN115047767A (zh) * | 2022-06-14 | 2022-09-13 | 北京中科飞鸿科技股份有限公司 | 一种考虑高度表使用约束的俯冲弹道优化方法 |
CN115079565B (zh) * | 2022-08-23 | 2022-10-25 | 中国人民解放军国防科技大学 | 变系数的带落角约束制导方法、装置和飞行器 |
CN115079565A (zh) * | 2022-08-23 | 2022-09-20 | 中国人民解放军国防科技大学 | 变系数的带落角约束制导方法、装置和飞行器 |
CN115857325A (zh) * | 2022-10-09 | 2023-03-28 | 北京航天自动控制研究所 | 一种基于多项式弹道的能量管理制导方法、存储介质 |
CN115309059B (zh) * | 2022-10-10 | 2023-02-03 | 北京航天自动控制研究所 | 一种考虑引力补偿的直接制导方法 |
CN115309059A (zh) * | 2022-10-10 | 2022-11-08 | 北京航天自动控制研究所 | 一种考虑引力补偿的直接制导方法 |
CN116379847A (zh) * | 2023-06-05 | 2023-07-04 | 西安现代控制技术研究所 | 一种导弹在线能量管理制导律构建方法 |
CN117073473A (zh) * | 2023-10-17 | 2023-11-17 | 中国空气动力研究与发展中心空天技术研究所 | 一种基于时间约束的导弹视角规划制导方法及系统 |
CN117073473B (zh) * | 2023-10-17 | 2024-01-02 | 中国空气动力研究与发展中心空天技术研究所 | 一种基于时间约束的导弹视角规划制导方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111442697A (zh) | 一种基于伪谱法修正的过重补制导方法和弹道整形制导方法 | |
CN108168381B (zh) | 一种多枚导弹协同作战的控制方法 | |
CN103090728B (zh) | 一种基于滑模控制的带末角约束制导方法 | |
CN112558631B (zh) | 基于测量信息的带大落角约束变参数制导方法 | |
CN111306998B (zh) | 一种参数摄动自适应的制导火箭弹垂直攻击制导方法 | |
CN111473696B (zh) | 一种基于落点估计的制导火箭垂直攻击制导方法 | |
CN111336871B (zh) | 一种基于迂回式飞行的垂直攻击制导方法 | |
CN107766967B (zh) | 一种基于多项式拟合法的拦截弹预测制导方法 | |
CN109737812B (zh) | 空对地制导武器侧向攻击方法和装置 | |
CN110717245B (zh) | 基于落角落速约束的拟滑翔弹道的设计方法 | |
CN111176315B (zh) | 一种基于l1自适应控制的变体巡飞弹控制方法 | |
CN111813137B (zh) | 一种靶机人在环控制方法 | |
CN113758380A (zh) | 考虑舵机动态性能和攻击角约束的制导控制一体化方法 | |
CN113359819B (zh) | 一种带有碰撞角约束和加速度限制的最优制导律 | |
CN114675673A (zh) | 一种空中动目标追踪方法及系统 | |
CN112836383B (zh) | 一种基于靶式数据修正的武器效能建模方法 | |
CN116301028B (zh) | 基于吸气式高超声速平台的多约束在线飞行轨迹规划中段导引方法 | |
CN116360489A (zh) | 一种基于成型制导律的协同中制导律设计方法 | |
CN112050693B (zh) | 考虑攻击角和视场约束的半捷联制导控制一体化设计方法 | |
RU2310152C1 (ru) | Способ стрельбы боевой машины по цели и система для его реализации | |
CN114815878A (zh) | 基于实时优化和深度学习的高超声速飞行器协同制导方法 | |
CN116956571A (zh) | 一种基于伪谱法优化的制导方法 | |
CN116793151A (zh) | 一种基于神经网络和伪谱法的制导方法 | |
CN114580138B (zh) | 一种基于贝塞尔牛顿的导弹多约束末制导律设计方法 | |
CN117932793B (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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20200724 |
|
WD01 | Invention patent application deemed withdrawn after publication |