CN103090728A - 一种基于滑模控制的带末角约束制导方法 - Google Patents
一种基于滑模控制的带末角约束制导方法 Download PDFInfo
- Publication number
- CN103090728A CN103090728A CN201310004137XA CN201310004137A CN103090728A CN 103090728 A CN103090728 A CN 103090728A CN 201310004137X A CN201310004137X A CN 201310004137XA CN 201310004137 A CN201310004137 A CN 201310004137A CN 103090728 A CN103090728 A CN 103090728A
- Authority
- CN
- China
- Prior art keywords
- gamma
- prime
- guidance
- angle
- aircraft
- 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
Images
Landscapes
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明涉及一种基于滑模控制的带末角约束的制导方法,属于制导技术领域。首先建立一种新的飞行器运动学和动力学模型,然后以制导末时刻,飞行器位置坐标与目标位置坐标(xf,yf)距离最小,并且为期望的末端弹道倾角γf为设计目标,根据反步法设计虚拟控制量,使得滑模函数及其导数在飞行末时刻同时到0;根据Lyapunov方法求解得到辅助控制量弹道倾角变化率γ',再将其转化为攻角α,输入最初建立的飞行器新模型,对飞行器轨迹进行实时调整,使其满足期望的终端条件,从而实现末制导。本发明方法考虑了飞行器的气动特性对制导过程的影响,更接近实际情况,且需要信息量少,可获得的弹道倾角末值范围广,得到的控制量变化平滑,易于姿态控制***进行跟踪。
Description
技术领域
本发明涉及一种基于滑模控制的带末角约束的制导方法,属于制导技术领域。
背景技术
在现代战争中,作战环境日趋多样,对制导律的要求也越来越高。某些特殊任务要求制导律不仅能够实现脱靶量最小,还期望实现以特定的末角对目标进行打击,例如:希望反坦克导弹能够垂直命中前装甲,从而使得毁伤力度最大;反导导弹能够头对头直接碰撞来袭导弹,从而尽最大可能摧毁作战弹头;某些反舰导弹希望从侧面对舰船进行攻击。传统的制导律(如比例制导)由于其信息需求量小,从而被广泛应用于实际作战***中。然而,此类制导律通常无法满足末角约束。因此,近年来,针对带末角约束的制导方法,各国学者展开了深入研究。
其中研究较为广泛的一类方法是以比例制导为基础的。Kim等人在传统的比例制导基础上加了一个时变偏置项,用来满足末角约束的要求;Lu等人提出了一种自适应的带末角约束制导律,该方法通过实时在线调整比例系数从而实现期望的末角约束;Ratnoo等人提出了另一种时变比例系数的带末角约束的制导律,该制导律分为两段,第一段为指向制导,该段比例系数是时变的,第二段为传统的比例制导,该段比例系数为常值。另一类广泛研究的方法是基于最优控制的。Ratnoo等人提出了一种基于状态相关Riccati方程(SDRE)的带末角约束的制导律,将制导问题转化成了非线性规划问题;Ryoo等人针对恒速导弹打击固定目标问题,给出了具有末角约束的最优制导律;Oza等人基于模型预测静态规划(MPSP)方法,针对地面目标提出了一种次优的带末角约束的制导律。
由于滑模控制方法对外部扰动具有鲁棒性,该方法也被广泛应用于制导律设计。Hou等人基于自适应滑模算法提出了一种打击地面固定目标的带末角约束的制导控制方法,然而该方法可获得的末角范围较窄;Harl等人基于滑模与反步法设计了一种同时对末角和飞行时间进行控制的末制导律;Kumar等人基于终端滑模设计了一种全方位带落角约束的制导控制方法,该方法可以用来打
击静止、常值速度以及机动目标。然而以上两种方法都进行了相应假设(如导
弹速度恒定),在真实情况下会对其结果造成较大影响。因此,需要提出一种简
单易行且符合实际情况的制导方法来解决该领域的问题。
发明内容
本发明为解决带末角约束的末制导问题,提出了一种基于滑模控制的带末角约束的制导方法。该方法根据末制导段的终端约束设计滑模函数,并结合反步法和Lyapunov方法得到制导控制量。
本发明的技术方案具体如下:
步骤1,建立二维平面飞行器的运动学和动力学模型:
其中,x,y是地面坐标系下的位置坐标(即射程和高度),V是飞行速度,γ为弹道倾角,m是飞行器质量,g是重力加速度,LD分别为升力和阻力,其中, ρ为大气密度,Cx,Cy分别为阻力系数和升力系数,是关于攻角和马赫的函数,Sref为飞行器的参考面积。
设计新的独立变量Y=y0-y (5)其中,y0是飞行器的初始高度。末制导段飞行器高度y单调递减,Y单调递增。
以Y作为独立变量,得到新模型如下:
步骤2,设计带末角约束的制导律
设计的目标为:在制导末时刻,飞行器位置坐标与目标位置坐标(xf,yf)距离最小,并且飞行器的弹道倾角为期望的末端弹道倾角γf。其中下标f表示变量末值。
步骤2.1,设计滑模函数
根据终端约束,设计滑模函数如下:
S1=x-xf-x'f(Y-Yf) (11)
将S1对Y求导,得到
S1'=x'-x'f (12)
为达到设计目标,设计控制律,使得滑模函数S1及其导数S1'在飞行末时刻同时收敛到0。
步骤2.2,求解辅助控制量
首先,根据反步法设计虚拟控制量,使得S1,S1'在飞行末时刻同时到0;然后根据Lyapunov方法求解得到辅助控制量。
将S1'作为虚拟控制量,为使得Y到达Yf时,S1,S1'同时到0,设计S1'为如下形式:
选取弹道倾角变化率γ'为辅助控制量,利用Lyapunov方法求解,得到辅助控制量使从某一时刻开始直到制导结束所设计的S1'形式一直成立,即式(13)成立。求解得到的γ'为:
γ'=-Msin2γ-ksin2γsgn(S2);
其中, k>0, γ为即时的弹道倾角。
步骤3,将辅助控制量转化为实际控制量
将步骤2中得到的作为辅助控制量的弹道倾角变化率γ'转化为攻角α。
将步骤2得到的γ'和即时状态代入新模型,即代入式(9),得到升力L,然后计算出升力系数Cy。攻角α和升力系数Cy存在一一对应的关系,由升力系数对飞行器气动数据插值,得到末制导段需要的攻角。
步骤4,将步骤3得到的攻角α输入步骤1建立的飞行器新模型,对飞行器轨迹进行实时调整,使其满足期望的终端条件,从而实现末制导。
有益效果
本发明有四方面优点:1.考虑了飞行器的气动特性(气动阻力和重力)对制导过程的影响,更接近实际情况。2.只需知道飞行的初始条件、末端要求和飞行器的即时状态信息便可实现末制导,需要信息量少。3.可获得的弹道倾角末值范围广。4.得到的控制量变化平滑,易于姿态控制***进行跟踪。
附图说明
图1为本发明的末制导段二维弹目几何关系;
图2为本发明方法的流程图;
图3为具体实施方式中不同末角约束条件下的飞行轨迹曲线;
图4为具体实施方式中不同末角约束条件下的弹道倾角变化曲线;
图5为具体实施方式中不同末角约束条件下的攻角变化曲线;
图6为具体实施方式中大气密度拉偏情况下的飞行轨迹曲线;
图7为具体实施方式中大气密度拉偏情况下的弹道倾角变化曲线;
图8为具体实施方式中大气密度拉偏情况下的攻角变化曲线;
图9为具体实施方式中大气密度拉偏情况下的速度变化曲线;
图10为具体实施方式中存在***滞后情况下的飞行轨迹曲线;
图11为具体实施方式中存在***滞后情况下的弹道倾角变化曲线;
图12为具体实施方式中存在***滞后情况下的攻角变化曲线。
具体实施方式
为了更好的说明本发明的目的和优点,下面结合附图和实例对技术方案做进一步详细说明。
1.带末角约束的制导律设计
选用某升力式再入飞行器,以垂直打击(即弹道倾角末值为-90度)地面目标为例进行介绍。飞行器起始位置坐标(x0,y0)为(0,20)km,起始弹道倾角为-3deg,初始速度为1700m/s。末端位置坐标(xf,yf)为(100,0)km,末端弹道倾角为-90deg。其他需要用到的参数选择如下:
n=4,∈=1e-3,p=0.3
末制导段二维弹目几何关系如图1所示,其中横轴为射程,纵轴为高度,三角形代表目标。本发明设计过程如图2所示。首先建立飞行器模型,从模型中获得飞行的即时状态,考虑终端约束条件,设计滑模函数S1;然后利用反步法设计滑模函数S2;利用Lyapunov方法求解得到辅助控制量(弹道倾角变化率),进而得到实际控制量(攻角)。在仿真过程中,将得到的控制指令输入给飞行器,从而控制其飞行轨迹满足期望的要求。具体步骤为:
步骤1,建立二维平面飞行器的运动学和动力学模型:
其中,x,y是地面坐标系下的位置坐标(即射程和高度),V是飞行速度,γ为弹道倾角,m是飞行器质量,g是重力加速度,LD分别为升力和阻力,其中, ρ为大气密度,Cx,Cy分别为阻力系数和升力系数,是关于攻角和马赫的函数,Sref为飞行器的参考面积。
以上模型是以时间作为独立变量,而在末制导研究中,飞行时间往往不是大家关心的主要问题。因此,设计新的独立变量为:
Y=y0-y (5)其中,y0是飞行器的初始高度。末制导段飞行器高度y单调递减,Y单调递增。
以Y作为独立变量,得到新模型如下:
步骤2,设计带末角约束的制导律
设计的目标为:在制导末时刻,飞行器位置坐标与目标位置坐标(xf,yf)距离最小,并且飞行器的弹道倾角为期望的末端弹道倾角γf。其中下标f表示变量末值。
步骤2.1,设计滑模函数
根据终端约束,设计滑模函数如下:
S1=x-xf-x'f(Y-Yf) (11)
将S1对Y求导,得到
S1'=x'-x'f (12)
在制导过程中,飞行器高度y逐渐接近目标高度yf,并在末时刻与其相等。因此,由式(5)可知,在飞行末点Y=y0-yf=Yf。从式(11)可以看出,如果在制导末时刻S1收敛到0,则可得到x=xf,即飞行器末点位置坐标与目标位置坐标距离最小。由式(6)可知,x'仅与弹道倾角γ有关,因此,根据式(12),如果在制导末时刻S1'也同时收敛到0,则可得到x'=x'f,即末角约束得到满足。
为达到设计目标,设计控制律,使得滑模函数S1及其导数S1'在飞行末时刻同时收敛到0。
步骤2.2,求解控制律
首先,根据反步法设计虚拟控制量,使得S1,S1’在飞行末时刻同时到0;然后根据Lyapunov方法求解得到辅助控制量。
以S1'为虚拟控制量,使得在Y到达Yf时,S1,S1'同时到0,其中S1'设计为如下形式:
将式(14)的积分初值分别设为S1b,Yb,可以得到如下积分等式
求解式(15)可以得到
将式(16)等式两边取指数即可得到S1的表达式如下
将式(17)对Y求导,即可得到下式
其中Yb=pYf,0<p<1,S1b是S1在Y=Yb时刻的值。
若在飞行过程中从Y=Yb时刻开始到制导末时刻,式(13)一直成立,则根据式(17)和式(18)可知:S1,S1'最后同时为0。
由于初始时刻式(13)不成立,因此选取弹道倾角变化率γ'为辅助控制量,利用Lyapunov方法求解,得到辅助控制量使得从Y=Yb时刻开始直到制导过程结束式(13)一直成立。具体求解过程为:
设计新的滑模函数如下
将S2对Y求导得到
令 得到
定义半正定的Lyapunov函数
将V1对Y求导得到
辅助控制量γ'=γeq'+γdisc' (24)
其中,γeq'是等效控制,使得S2'=0;γdisc'是非连续切换控制,使得滑模面S2=0具有吸引性。
等效控制γeq′=-Msin2γ (25)
将式(25)带入式(24),再将结果带入式(23),得到
非连续切换控制γdisc'=-ksin2γsgn(S2)k>0 (27)
将式(27)带入式(26)得到V1'=-k|S2| (28)
根据Lyapunov稳定性理论,得出***是渐进稳定的。
将式(25)和式(27)带入式(24),再将结果带入式(21)得到:
S2'=-ksgn(S2) (29)
从式(29)得到,切换增益k表示滑模函数向滑模面的趋近速度。由于设计期望为在Y=Yb时刻,滑模面S2=0。因此,切换增益计算如下
其中S2(0)是滑模函数S2的初值。从而滑模函数S2以常值速度k趋近于滑模面S2=0,并在Y=Yb时到达滑模面。在此之后,***状态一直保持在滑模面S2=0上,即式(13)一直成立,从而保证在制导末时刻,S1,S1'同时到0。
由于制导过程在Y=Yf时刻结束,上述方法所能获得的弹道倾角末值范围为:(-180,0)deg,0deg和-180deg的反正切值不存在因而排除在外。
为了减弱在沿着滑模面滑动过程中,由于使用切换函数sgn(·)引起的抖振现象,用饱和函数取代切换函数,其中∈是边界层。该设计能使滑模函数收敛到边界层内,若边界层厚度选取足够小,则滑模函数近似收敛到0。
步骤3,将辅助控制量转化为实际控制量
将步骤2中得到的作为辅助控制量的弹道倾角变化率γ'转化为攻角α。
将步骤2得到的γ'和即时状态带入式(9),得到升力L,然后计算出升力系数Cy。攻角α和升力系数Cy存在一一对应的关系,由升力系数对飞行器气动数据插值,即可得到末制导段需要的攻角。
步骤4,将步骤3得到的攻角α输入步骤1建立的飞行器模型,对飞行器轨迹进行实时调整,使其满足期望的终端条件,从而实现末制导。
2.验证本发明提出的制导律的有效性
针对不同情况对该发明的有效性进行验证。首先,验证该发明提出的制导律可满足不同的弹道倾角末值需求;然后,验证该发明提出的制导律能够应对大气密度拉偏的情况;最后,验证在***存在高阶滞后的情况下该发明提出的制导律仍然有着很高的精度。以下仿真结果,最终脱靶量均在1e-3m以内,末角误差均在1e-4deg以内。
①不同的末端弹道倾角的情况
在本实施例中,目标横坐标为60km,要求的末端弹道倾角分别为-10deg,-30deg,-90deg,-150deg,-170deg。图3为不同末角约束条件下的飞行轨迹曲线,图4为不同末角约束条件下的弹道倾角变化曲线,图5为不同末角约束条件下的攻角变化曲线。由图3和图4中可以看出,飞行器在击中目标的同时,所有的末角约束情况都得到了满足。在实际应用中,通常会对控制量进行限幅,在本发明中,攻角变化范围均限定在(-30,+30)deg之间。图5中与末角要求为-150deg和-170deg相对应的攻角曲线均存在控制量饱和现象,这是由于为了实现该末角约束,飞行轨迹需要很大程度的弯曲,导致控制量需求很大。
②大气密度拉偏的情况
在本实施例中,分别将大气密度正向、负向拉偏20%进行仿真,并与标称情况进行对比分析。图6为大气密度拉偏情况下的飞行轨迹曲线,图7为大气密度拉偏情况下的弹道倾角变化曲线,图8为大气密度拉偏情况下的攻角变化曲线,图9为大气密度拉偏情况下的速度变化曲线。由图6可以看出,大气密度拉偏对实际飞行轨迹几乎没有产生任何影响。结合图7,可以看出在这三种情况下本发明提出的制导律均能准确满足末角约束。图8和图9说明大气密度拉偏对控制量和飞行速度有较大影响,大气密度越大,飞行速度下降越快,相应的控制量峰值越小。这是因为飞行器所受的阻力随着大气密度的增大而增大。
③存在***滞后的情况
在实际情况中,制导***给出的指令不能被立即实现,而是需要经过一定的时间滞后。因此,在本实例中引进了***滞后来验证本发明所提出的制导律的鲁棒性。由于一阶滞后往往不能表现实际情况,本例引入五阶滞后如下式所示:
其中,τ为滞后时间常数,αc为制导***给出的指令,αr为实现的指令。
图10为在***滞后情况下的飞行轨迹曲线,图11为在***滞后情况下的弹道倾角变化曲线,图12为在***滞后情况下的攻角变化曲线。从图12可以看出,虽然***存在5阶滞后,然而实现的控制量曲线和期望的控制量曲线几乎完全重合,这是由于本发明提出的制导律产生的控制指令变化非常平缓。图10和图11说明在这种情况下制导律对给定末端指标要求仍然可以精确达到。因此,该制导律在实际应用中也可以达到非常好的效果。
综上所述,该发明提出的制导律形式简单,鲁棒性强,不但能够使得脱靶量最小,而且可以精确地满足末角约束,具有很高的工程应用价值。
Claims (5)
1.一种基于滑模控制的带末角约束制导方法,其特征在于:具体实现步骤如下:
步骤1,建立二维平面飞行器的运动学和动力学模型:
其中,x,y是地面坐标系下的位置坐标,V是飞行速度,γ为弹道倾角,m是飞行器质量,g是重力加速度,LD分别为升力和阻力, ρ为大气密度,Cx,Cy分别为阻力系数和升力系数,Sref为飞行器的参考面积;
设计新的独立变量Y=y0-y (5)其中,y0是飞行器的初始高度;
以Y作为独立变量,得到新模型如下:
步骤2,设计带末角约束的制导律
设计的目标为:在制导末时刻,飞行器位置坐标与目标位置坐标(xf,yf)距离最小,并且飞行器的弹道倾角为期望的末端弹道倾角γf;其中下标f表示变量末值;
步骤2.1,设计滑模函数
根据终端约束,设计滑模函数如下:
S1=x-xf-x'f(Y-Yf) (11)
将S1对Y求导,得到
S1'=x'-x'f (12)
步骤2.2,求解辅助控制量
为使得S1,S1'在飞行末时刻同时到0,根据反步法设计虚拟控制量;
将S1'作为虚拟控制量,为使得Y到达Yf时,'S1S1'同时到0,设计S1'为如下形式:
再选取弹道倾角变化率γ'为辅助控制量,利用Lyapunov方法求解,得到辅助控制量使从某一时刻开始直到制导结束所设计的S1’形式一直成立;
求解得到的γ'为:
γ'=-Msin2γ-ksin2γsgn(S2);
其中, k>0, γ为即时的弹道倾角;
步骤3,将步骤2中得到的作为辅助控制量γ'转化为攻角α;
步骤4,将步骤3得到的攻角α输入步骤1建立的飞行器新模型,对飞行器轨迹进行实时调整,使其满足期望的终端条件,实现末制导。
2.根据权利要求1所述的一种基于滑模控制的带末角约束制导方法,其特征在于:γ'的具体求解过程为:
设计新的滑模函数如下
将S2对Y求导得到
令 得到
定义半正定的Lyapunov函数
将V1对Y求导得到
辅助控制量γ'=γeq'+γdisc'(24)
其中,γeq'是等效控制,使得S2'=0;γdisc'是非连续切换控制,使得滑模面S2=0具有吸引性;
等效控制γeq'=-Msin2γ (25)
从而得到
非连续切换控制γdisc'=-ksin2γSgn(S2)k>0 (27)
V1'=-k|S2| (28)
根据Lyapunov稳定性理论,***是渐进稳定的;
S2'=-ksgn(S2) (29)
切换增益为
其中S2(0)是滑模函数S2的初值。
3.根据权利要求1所述的一种基于滑模控制的带末角约束制导方法,其特征在于:用饱和函数取代切换函数,其中∈是边界层。
4.根据权利要求1所述的一种基于滑模控制的带末角约束制导方法,其特征在于:末制导段飞行器高度y单调递减,Y单调递增。
5.根据权利要求1所述的一种基于滑模控制的带末角约束制导方法,其特征在于:步骤3所述辅助控制量γ'转化为攻角α的具体方法为:将步骤2得到的γ'和即时状态代入新模型,得到升力L,然后计算出升力系数Cy;由升力系数对飞行器气动数据插值,得到末制导段需要的攻角。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310004137.XA CN103090728B (zh) | 2013-01-07 | 2013-01-07 | 一种基于滑模控制的带末角约束制导方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310004137.XA CN103090728B (zh) | 2013-01-07 | 2013-01-07 | 一种基于滑模控制的带末角约束制导方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103090728A true CN103090728A (zh) | 2013-05-08 |
CN103090728B CN103090728B (zh) | 2014-07-23 |
Family
ID=48203607
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310004137.XA Expired - Fee Related CN103090728B (zh) | 2013-01-07 | 2013-01-07 | 一种基于滑模控制的带末角约束制导方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103090728B (zh) |
Cited By (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103512426A (zh) * | 2013-09-06 | 2014-01-15 | 北京理工大学 | 一种次优的带末角约束制导方法 |
CN103728976A (zh) * | 2013-12-30 | 2014-04-16 | 北京航空航天大学 | 一种基于广义标控脱靶量概念的多过程约束和多终端约束末制导律 |
CN103838914A (zh) * | 2013-12-30 | 2014-06-04 | 北京航空航天大学 | 一种高超声速飞行器滑翔段弹道解析求解方法 |
CN104155990A (zh) * | 2014-08-15 | 2014-11-19 | 哈尔滨工业大学 | 考虑攻角约束的高超声速飞行器俯仰通道姿态控制方法 |
CN104156595A (zh) * | 2014-08-11 | 2014-11-19 | 北京航天自动控制研究所 | 飞行器的飞行轨迹指令的确定方法和装置 |
CN104634182A (zh) * | 2014-12-16 | 2015-05-20 | 北京控制工程研究所 | 一种跳跃式再入标准弹道在线修正的跟踪制导方法 |
CN104792232A (zh) * | 2015-04-28 | 2015-07-22 | 北京理工大学 | 一种带有落角约束的最小过载末端导引方法 |
CN104950898A (zh) * | 2015-06-10 | 2015-09-30 | 北京理工大学 | 一种再入飞行器全阶非奇异终端滑模姿态控制方法 |
CN104965519A (zh) * | 2015-06-10 | 2015-10-07 | 北京理工大学 | 一种基于贝塞尔曲线的带落角约束的末制导方法 |
CN105159308A (zh) * | 2015-08-28 | 2015-12-16 | 北京控制工程研究所 | 一种可重复使用飞行器着陆段制导与控制律一体化耦合设计方法 |
CN105785762A (zh) * | 2016-03-17 | 2016-07-20 | 北京航空航天大学 | 一种基于自适应反步滑模的两轴惯性稳定平台高精度控制方法 |
CN105867403A (zh) * | 2016-05-18 | 2016-08-17 | 北京航天自动控制研究所 | 一种使飞行器满足强姿态约束条件的方法 |
CN106091816A (zh) * | 2016-05-27 | 2016-11-09 | 北京航空航天大学 | 一种基于滑模变结构理论的半捷联空空导弹制导方法 |
CN106370059A (zh) * | 2016-08-26 | 2017-02-01 | 方洋旺 | 一种随机快速光滑二阶滑模末制导方法 |
CN106647276A (zh) * | 2016-12-30 | 2017-05-10 | 中国人民解放军国防科学技术大学 | 一种无动力飞行器弹道末端位置与姿态平滑控制方法 |
CN108050888A (zh) * | 2017-11-27 | 2018-05-18 | 北京航空航天大学 | 一种具有禁飞区约束的显式制导方法 |
CN108628332A (zh) * | 2018-05-28 | 2018-10-09 | 浙江工业大学 | 基于反正切增强型指数趋近律和快速终端滑模面的四旋翼飞行器有限时间控制方法 |
CN109283932A (zh) * | 2018-09-18 | 2019-01-29 | 浙江工业大学 | 一种基于积分反步滑模的四旋翼飞行器姿态控制方法 |
CN109471454A (zh) * | 2018-12-07 | 2019-03-15 | 湖北航天飞行器研究所 | 一种指定攻击倾角的微型作业飞行器的末端制导段进入方法 |
CN110618608A (zh) * | 2018-06-20 | 2019-12-27 | 河南科技大学 | 一种复合制导跟踪控制方法与装置 |
CN111306998A (zh) * | 2020-03-24 | 2020-06-19 | 北京理工大学 | 一种参数摄动自适应的制导火箭弹垂直攻击制导方法 |
CN111649624A (zh) * | 2020-07-30 | 2020-09-11 | 哈尔滨工业大学 | 一种空间微型精确制导武器控制方法 |
CN113051743A (zh) * | 2021-03-18 | 2021-06-29 | 中国人民解放***箭军工程大学 | 一种基于轨迹在线规划的末制导*** |
CN113656887A (zh) * | 2021-07-28 | 2021-11-16 | 上海机电工程研究所 | 基于坐标偏置优化的弹目运动仿真能力拓展方法和*** |
CN113759966A (zh) * | 2021-08-26 | 2021-12-07 | 北京理工大学 | 一种三维空间内终端速度可控的末制导方法 |
CN114020019A (zh) * | 2021-11-03 | 2022-02-08 | 北京理工大学 | 飞行器的制导方法与装置 |
CN114216368A (zh) * | 2021-12-07 | 2022-03-22 | 南京理工大学 | 一种变速条件下弹群协同制导方法及*** |
CN114237299A (zh) * | 2021-12-22 | 2022-03-25 | 北京航空航天大学 | 一种飞行器制导方法、装置、电子设备及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5263662A (en) * | 1992-05-19 | 1993-11-23 | United Technologies Corporation | Helicopter integrated fire and flight control system having turn coordination control |
CN102410785A (zh) * | 2011-12-13 | 2012-04-11 | 林德福 | 激光末制导炮弹半实物仿真用的转台 |
-
2013
- 2013-01-07 CN CN201310004137.XA patent/CN103090728B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5263662A (en) * | 1992-05-19 | 1993-11-23 | United Technologies Corporation | Helicopter integrated fire and flight control system having turn coordination control |
CN102410785A (zh) * | 2011-12-13 | 2012-04-11 | 林德福 | 激光末制导炮弹半实物仿真用的转台 |
Non-Patent Citations (1)
Title |
---|
张亚松等: "带落角约束的滑模变结构制导律研究", 《电光与控制》 * |
Cited By (46)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103512426B (zh) * | 2013-09-06 | 2015-05-06 | 北京理工大学 | 一种次优的带末角约束制导方法 |
CN103512426A (zh) * | 2013-09-06 | 2014-01-15 | 北京理工大学 | 一种次优的带末角约束制导方法 |
CN103728976B (zh) * | 2013-12-30 | 2016-04-27 | 北京航空航天大学 | 一种基于广义标控脱靶量概念的多过程约束和多终端约束末制导律 |
CN103728976A (zh) * | 2013-12-30 | 2014-04-16 | 北京航空航天大学 | 一种基于广义标控脱靶量概念的多过程约束和多终端约束末制导律 |
CN103838914A (zh) * | 2013-12-30 | 2014-06-04 | 北京航空航天大学 | 一种高超声速飞行器滑翔段弹道解析求解方法 |
CN103838914B (zh) * | 2013-12-30 | 2017-01-18 | 北京航空航天大学 | 一种高超声速飞行器滑翔段弹道解析求解方法 |
CN104156595A (zh) * | 2014-08-11 | 2014-11-19 | 北京航天自动控制研究所 | 飞行器的飞行轨迹指令的确定方法和装置 |
CN104156595B (zh) * | 2014-08-11 | 2015-06-24 | 北京航天自动控制研究所 | 飞行器的飞行轨迹指令的确定方法和装置 |
CN104155990A (zh) * | 2014-08-15 | 2014-11-19 | 哈尔滨工业大学 | 考虑攻角约束的高超声速飞行器俯仰通道姿态控制方法 |
CN104155990B (zh) * | 2014-08-15 | 2016-09-14 | 哈尔滨工业大学 | 考虑攻角约束的高超声速飞行器俯仰通道姿态控制方法 |
CN104634182B (zh) * | 2014-12-16 | 2016-02-10 | 北京控制工程研究所 | 一种跳跃式再入标准弹道在线修正的跟踪制导方法 |
CN104634182A (zh) * | 2014-12-16 | 2015-05-20 | 北京控制工程研究所 | 一种跳跃式再入标准弹道在线修正的跟踪制导方法 |
CN104792232A (zh) * | 2015-04-28 | 2015-07-22 | 北京理工大学 | 一种带有落角约束的最小过载末端导引方法 |
CN104792232B (zh) * | 2015-04-28 | 2016-04-20 | 北京理工大学 | 一种带有落角约束的最小过载末端导引方法 |
CN104950898A (zh) * | 2015-06-10 | 2015-09-30 | 北京理工大学 | 一种再入飞行器全阶非奇异终端滑模姿态控制方法 |
CN104965519A (zh) * | 2015-06-10 | 2015-10-07 | 北京理工大学 | 一种基于贝塞尔曲线的带落角约束的末制导方法 |
CN104965519B (zh) * | 2015-06-10 | 2018-08-10 | 北京理工大学 | 一种基于贝塞尔曲线的带落角约束的末制导方法 |
CN105159308A (zh) * | 2015-08-28 | 2015-12-16 | 北京控制工程研究所 | 一种可重复使用飞行器着陆段制导与控制律一体化耦合设计方法 |
CN105785762A (zh) * | 2016-03-17 | 2016-07-20 | 北京航空航天大学 | 一种基于自适应反步滑模的两轴惯性稳定平台高精度控制方法 |
CN105867403A (zh) * | 2016-05-18 | 2016-08-17 | 北京航天自动控制研究所 | 一种使飞行器满足强姿态约束条件的方法 |
CN106091816A (zh) * | 2016-05-27 | 2016-11-09 | 北京航空航天大学 | 一种基于滑模变结构理论的半捷联空空导弹制导方法 |
CN106370059A (zh) * | 2016-08-26 | 2017-02-01 | 方洋旺 | 一种随机快速光滑二阶滑模末制导方法 |
CN106647276A (zh) * | 2016-12-30 | 2017-05-10 | 中国人民解放军国防科学技术大学 | 一种无动力飞行器弹道末端位置与姿态平滑控制方法 |
CN106647276B (zh) * | 2016-12-30 | 2019-09-17 | 中国人民解放军国防科学技术大学 | 一种无动力飞行器弹道末端位置与姿态平滑控制方法 |
CN108050888A (zh) * | 2017-11-27 | 2018-05-18 | 北京航空航天大学 | 一种具有禁飞区约束的显式制导方法 |
CN108050888B (zh) * | 2017-11-27 | 2019-10-29 | 北京航空航天大学 | 一种具有禁飞区约束的显式制导方法 |
CN108628332B (zh) * | 2018-05-28 | 2021-08-03 | 浙江工业大学 | 基于反正切增强型指数趋近律和快速终端滑模面的四旋翼飞行器有限时间控制方法 |
CN108628332A (zh) * | 2018-05-28 | 2018-10-09 | 浙江工业大学 | 基于反正切增强型指数趋近律和快速终端滑模面的四旋翼飞行器有限时间控制方法 |
CN110618608A (zh) * | 2018-06-20 | 2019-12-27 | 河南科技大学 | 一种复合制导跟踪控制方法与装置 |
CN110618608B (zh) * | 2018-06-20 | 2022-04-15 | 河南科技大学 | 一种复合制导跟踪控制方法与装置 |
CN109283932A (zh) * | 2018-09-18 | 2019-01-29 | 浙江工业大学 | 一种基于积分反步滑模的四旋翼飞行器姿态控制方法 |
CN109471454A (zh) * | 2018-12-07 | 2019-03-15 | 湖北航天飞行器研究所 | 一种指定攻击倾角的微型作业飞行器的末端制导段进入方法 |
CN109471454B (zh) * | 2018-12-07 | 2021-11-26 | 湖北航天飞行器研究所 | 一种指定攻击倾角的微型作业飞行器的末端制导段进入方法 |
CN111306998A (zh) * | 2020-03-24 | 2020-06-19 | 北京理工大学 | 一种参数摄动自适应的制导火箭弹垂直攻击制导方法 |
CN111306998B (zh) * | 2020-03-24 | 2021-01-26 | 北京理工大学 | 一种参数摄动自适应的制导火箭弹垂直攻击制导方法 |
CN111649624A (zh) * | 2020-07-30 | 2020-09-11 | 哈尔滨工业大学 | 一种空间微型精确制导武器控制方法 |
CN113051743A (zh) * | 2021-03-18 | 2021-06-29 | 中国人民解放***箭军工程大学 | 一种基于轨迹在线规划的末制导*** |
CN113051743B (zh) * | 2021-03-18 | 2023-05-26 | 中国人民解放***箭军工程大学 | 一种基于轨迹在线规划的末制导*** |
CN113656887A (zh) * | 2021-07-28 | 2021-11-16 | 上海机电工程研究所 | 基于坐标偏置优化的弹目运动仿真能力拓展方法和*** |
CN113656887B (zh) * | 2021-07-28 | 2024-04-05 | 上海机电工程研究所 | 基于坐标偏置优化的弹目运动仿真能力拓展方法和*** |
CN113759966A (zh) * | 2021-08-26 | 2021-12-07 | 北京理工大学 | 一种三维空间内终端速度可控的末制导方法 |
CN113759966B (zh) * | 2021-08-26 | 2023-10-20 | 北京理工大学 | 一种三维空间内终端速度可控的末制导方法 |
CN114020019A (zh) * | 2021-11-03 | 2022-02-08 | 北京理工大学 | 飞行器的制导方法与装置 |
CN114020019B (zh) * | 2021-11-03 | 2023-12-01 | 北京理工大学 | 飞行器的制导方法与装置 |
CN114216368A (zh) * | 2021-12-07 | 2022-03-22 | 南京理工大学 | 一种变速条件下弹群协同制导方法及*** |
CN114237299A (zh) * | 2021-12-22 | 2022-03-25 | 北京航空航天大学 | 一种飞行器制导方法、装置、电子设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN103090728B (zh) | 2014-07-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103090728B (zh) | 一种基于滑模控制的带末角约束制导方法 | |
CN104516356B (zh) | 基于rrt的动态障碍规避算法 | |
CN106681348A (zh) | 考虑全捷联导引头视场约束的制导控制一体化设计方法 | |
CN104019701B (zh) | 一种利用直接力气动力复合控制的前向拦截制导方法 | |
CN103983143B (zh) | 包含速度过程约束和多终端约束的空地导弹投放下滑段制导方法 | |
CN103245257B (zh) | 基于Bezier曲线的多约束飞行器导引方法 | |
CN111591470B (zh) | 一种适应推力可调模式的飞行器精确软着陆闭环制导方法 | |
CN102927851A (zh) | 一种基于轨迹在线规划的末制导方法 | |
CN111336871B (zh) | 一种基于迂回式飞行的垂直攻击制导方法 | |
CN114020019B (zh) | 飞行器的制导方法与装置 | |
CN107621198B (zh) | 一种绕多禁飞区的多步决策弹道规划方法 | |
CN105629734A (zh) | 一种近空间飞行器的轨迹跟踪控制方法 | |
CN104199303A (zh) | 一种基于矢量场制导的平流层卫星平面路径跟踪控制方法 | |
CN115406312B (zh) | 考虑视场角和舵机延时约束的导弹制导控制一体化方法 | |
Zhao et al. | New three-dimensional second-order sliding mode guidance law with impact-angle constraints | |
CN111897223B (zh) | 一种考虑自动驾驶仪动态特性的速度追踪制导方法 | |
CN107766967A (zh) | 一种基于多项式拟合法的拦截弹预测制导方法 | |
CN117171877A (zh) | 基于时机博弈的高超声速飞行器机动突防策略设计方法 | |
CN113758380A (zh) | 考虑舵机动态性能和攻击角约束的制导控制一体化方法 | |
CN109827478A (zh) | 一种带落角约束与过载约束的制导方法 | |
CN108050888B (zh) | 一种具有禁飞区约束的显式制导方法 | |
CN110471283A (zh) | 一种带碰撞角约束的三维鲁棒制导律构建方法 | |
CN107422744B (zh) | 一种基于径向速度控制的交会时间控制方法 | |
Viswanath et al. | Disturbance observer based sliding mode control for proportional navigation guidance | |
CN114815878A (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: 20140723 Termination date: 20160107 |
|
CF01 | Termination of patent right due to non-payment of annual fee |