文章快速检索  
  高级检索
过虚拟交班点的能量最优制导律
李晨迪1,2, 王江1,2, 李斌1,2, 何绍溟1,2, 张彤3     
1. 北京理工大学 宇航学院, 北京 100081;
2. 北京理工大学 无人机自主控制技术北京市重点实验室, 北京 100081;
3. 北方华安工业集团, 齐齐哈尔 161006
摘要: 针对固定目标的导弹过虚拟交班点制导问题,在希尔伯特空间下,基于最优化理论设计了有无终端落角约束两种情况下的全局能量最优制导律。通过对模型进行线性化,将提出的最优制导模型转化为线性二次型最优控制问题,在此基础上利用零控脱靶量(ZEM)概念对系统模型进行降阶,并推导出解析解。设计的制导律可以使导弹准确经过虚拟交班点,并实现期望终端落角。仿真结果表明,与经典制导律对比,该制导律可以显著减少全局控制能量消耗。
关键词: 最优制导     虚拟交班点     终端落角约束     轨迹跟踪     比例导引     弹道成型制导律    
Energy-optimal guidance law with virtual hand-over point
LI Chendi1,2, WANG Jiang1,2, LI Bin1,2, HE Shaoming1,2, ZHANG Tong3     
1. School of Aerospace Engineering, Beijing Institute of Technology, Beijing 100081, China;
2. Beijing Key Laboratory of UAV Autonomous Control, Beijing Institute of Technology, Beijing 100081, China;
3. Hua An Industry Group Company Limited, Qiqihar 161006, China
Abstract: Aiming at the problem of fixed target missile guidance with virtual hand-over point, a global energy optimal guidance law with/without terminal angle constraint is designed based on the optimization theory in Hilbert space. By linearizing the model, the proposed optimal guidance model is transformed into a linear quadratic opti-mal control problem. Then the concept of Zero-Effort-Miss (ZEM) is used to reduce the order of the system, deriving the analytical solution. The proposed guidance law can ensure the missile pass through the virtual hand-over point accurately and reached the desired terminal angle. The simulation results show that, compared with the classical guidance law, the proposed guidance law can significantly reduce the energy consumption of global control.
Keywords: optimal guidance     virtual hand-over point     terminal angle constraint     path following     proportional navigation guidance     trajectory shaping guidance    

随着现代军事技术的发展,对精确制导武器的需求与日俱增,超视距导弹武器系统一般采用复合制导体制[1-3]。复合制导中,导弹由中制导转为末制导时,交接班问题是实现精确制导的关键。为了使导引头开机时顺利截获目标,实现中末制导交班,往往对此时导弹的位置等信息有一定要求[1],为了满足对导弹中末制导交班的弹道约束,本文研究过虚拟交班点的制导律,为导弹提供可靠的中末制导交班条件。

近年来,基于轨迹跟踪的制导律在无人飞行器领域被广泛研究[4-5]。根据跟踪方式,轨迹跟踪法可分为两类:连续轨迹跟踪和路径点(弹道点)跟踪。与连续轨迹跟踪相比,路径点跟踪只需要选取有限个弹道特征点,通过使导弹经过选定的弹道点而得到期望弹道,计算量较小。连续轨迹跟踪往往需要通过复杂计算实时得到待飞轨迹长度,或者基于虚拟目标进行轨迹跟踪[6]。文献[7]提出了一种两点间线性二次最优轨迹跟踪制导律,用于跟踪两个弹道点之间的直线,其路径改变点同样基于需用过载最小得到,但该制导律不是全局最优。文献[8]提出了一种数值优化方法,使得轨迹跟踪制导律的能量消耗最小,计算复杂度较高。文献[9]以反舰导弹为应用背景,通过离线参数优化得到每个弹道点的最优边界条件,例如角度和过载,从而在两个弹道点之间采用最优攻击角度制导律,得到了全局能量最优的轨迹跟踪制导律。

有时,需要导弹以特定的终端落角对目标进行攻击,以增强毁伤效果[10],近年来考虑终端角度约束的制导律受到广泛关注。文献[11-12]通过求解含剩余飞行时间的线性二次最优控制方程得到含终端角度约束的最优制导律的一般形式。文献[13]应用Schwarz不等式研究了带落角约束的任意加权最优制导律,得到了制导律的一般表达式。文献[14]针对高速/低速目标拦截问题,采用时变偏置角速率设计了一种三维联合偏置比例导引律。文献[15]研究了一种无剩余飞行时间的偏置比例导引律,在比例导引的基础上加上角度约束偏置项,通过改变导弹过载实现终端角度约束。文献[16]针对空地导弹具有终端角度约束的制导问题,提出了一种在有限时间内稳定的新型二阶滑模变结构制导律。

路径点跟踪制导问题在每两个相邻路径点之间可以看作一般的制导问题,因此近几十年研究得到的一些经典导弹制导律同样可以应用于路径点跟踪问题。对于不考虑制导动力学滞后的制导问题,无终端落角约束时,导航比为3的比例导引为最优制导律[17],当考虑终端落角约束时,弹道成型制导律则成为最优制导律[18-20]。但是在多个弹道点之间连续采用上述制导律的全局能量消耗是否为最优尚未得到研究验证[21]。基于此,本文将过虚拟交班点的精确制导问题视为路径点跟踪的特殊情况,即包括虚拟交班点与目标点两个弹道点,以气动力控制的导弹为研究对象,分别研究了无落角约束与考虑终端落角约束情况下的全局最优制导问题,推导出了解析解。最后,对非线性模型进行仿真,结果表明,与经典最优制导律相比,无落角约束与考虑终端落角约束的情况下,该制导律可以将控制能量消耗分别减小49%和22%。

1 运动模型建立 1.1 非线性运动模型

为方便描述,将虚拟交班点与目标落点视为两个导弹需要经过的路径点,记作Pi(i∈{1, 2}),导弹与目标路径点Pi之间的相对几何关系如图 1所示。图中:XIOYI为惯性坐标系;Va分别为导弹速度与法向过载;θ为弹道倾角;riqiσi分别为导弹与第i个目标路径点之间的相对距离、弹目视线角和速度方向误差角。

图 1 弹目相对运动关系 Fig. 1 Missile-target engagement relationship

导弹与目标路径点之间的相对运动方程为

$ \left\{ {\begin{array}{*{20}{l}} {{{\dot r}_i} = - V\cos {\sigma _i}}\\ {{{\dot q}_i} = - \frac{{V\sin {\sigma _i}}}{{{r_i}}}\;\;\;\;i \in \left\{ {1,2} \right\}}\\ {\dot \theta = \frac{a}{V}} \end{array}} \right. $ (1)
 

由几何关系有

$ {\sigma _i} = \theta - {q_i} $ (2)
 
1.2 模型线性化

本文将基于线性模型推导最优制导律的解析解。引入参考坐标系,对建立的非线性运动学模型进行线性化。如图 1所示,将坐标系OXIYI绕原点旋转qR角,得到参考坐标系OXRYR。令yi=yPiyM为导弹与Pi之间垂直于XR轴方向的相对位移,由几何知识,yi=0 $ \Leftrightarrow $ ri=0。

在参考坐标系XROYR下,弹目相对运动方程可表示为

$ \left\{ {\begin{array}{*{20}{l}} {{{\dot y}_i} = {v_i}}\\ {{{\dot v}_i} = - {a_ \bot }\cos \left( {{q_i} - {q_{\rm{R}}}} \right)} \end{array}\quad i \in \left\{ {1,2} \right\}} \right. $ (3)
 

式中:viai分别为导弹与Pi之间垂直于XR轴方向的相对速度与过载,由几何关系可得

$ {a_ \bot } = a\cos {\sigma _i} $ (4)
 

在工程实际当中,飞行弹道由无控段、中制导段与末制导段3部分组成,对于中制导段与末制导段,可在每段飞行过程中分别视为弹目视线角变化较小,因此通过选取合适的qR可以使qiqR为小量。同时,在实际情况中,导弹速度方向误差角σi也为小量。基于以上合理假设,式(3)可改写为

$ \left\{ {\begin{array}{*{20}{l}} {{{\dot y}_i} = {v_i}}\\ {{{\dot v}_i} = - a} \end{array}\;\;\;\;i \in \left\{ {1,2} \right\}} \right. $ (5)
 

记系统状态向量为x=[x1, x2]T=[y1, v1, y2, v2]T,输出向量为y=[y1, y2]T,线性化的系统状态方程可以写为

$ \left\{ {\begin{array}{*{20}{l}} {\mathit{\boldsymbol{\dot x}} = \mathit{\boldsymbol{Ax}} + \mathit{\boldsymbol{B}}a}\\ {\mathit{\boldsymbol{y}} = \mathit{\boldsymbol{Cx}}} \end{array}} \right. $ (6)
 

式中:

$ \left\{ {\begin{array}{*{20}{l}} {\mathit{\boldsymbol{A}} = \left[ {\begin{array}{*{20}{c}} {{\mathit{\boldsymbol{A}}_1}}&{\bf{0}}\\ {\bf{0}}&{{\mathit{\boldsymbol{A}}_2}} \end{array}} \right],{\mathit{\boldsymbol{A}}_i} = \left[ {\begin{array}{*{20}{l}} 0&1\\ 0&0 \end{array}} \right]}\\ {\mathit{\boldsymbol{B}} = {{\left[ {0, - 1,0, - 1} \right]}^{\rm{T}}}}\\ {\mathit{\boldsymbol{C}} = \left[ {1,0,1,0} \right]} \end{array}} \right. $ (7)
 
2 最优控制问题描述

为得到控制量的全局最优解,提出以下形式的性能指标:

$ J = \left\{ {\begin{array}{*{20}{l}} {\int_t^{{t_{{\rm{f}},1}}} {{a^2}\left( \tau \right){\rm{d}}\tau } + \int_{{t_{{\rm{f}},1}}}^{{t_{{\rm{f}},2}}} {{a^2}\left( \tau \right){\rm{d}}\tau } \quad t \le {t_{{\rm{f,1}}}},}\\ {\int_t^{{t_{{\rm{f}},2}}} {{a^2}\left( \tau \right){\rm{d}}\tau } \quad {t_{{\rm{f}},1}} < t \le {t_{{\rm{f}},2}}} \end{array}} \right. $ (8)
 

式中:tf, i为导弹到达第i个路径点的时刻。

本文研究的过虚拟交班点的固定目标最优控制问题可以表述如下。

问题1   过定点的最优制导律问题

对线性化的系统(6),解得使性能指标(8)取极小值的制导指令a,并满足约束:

$ {y_i}\left( {{t_{{\rm{f}},i}}} \right) = 0\;\;\;\;i \in \left\{ {1,2} \right\} $ (9)
 

问题2  带终端落角约束的过定点最优制导律问题

对于带有终端落角约束的最优控制问题,在上述无落角约束情况下,需额外满足约束

$ \theta \left( {{t_{{\rm{f}},2}}} \right) = {\theta _{\rm{d}}} $ (10)
 

式中:θd为期望落角。

3 过定点的最优制导律设计

本节根据取控制能量最小的性能指标,选取零控脱靶量(ZEM)[19]作为唯一状态变量,对系统进行变换,从而实现系统降阶,对上述两个最优控制问题进行推导,以得到最优制导律的解析解,并对其进行分析。

3.1 系统降阶

为了方便对最优解的解析形式进行推导,对状态向量进行坐标变换,使系统降阶。以第i个点为例,进行变换

$ {Z_i} = \left\{ {\begin{array}{*{20}{l}} {{\mathit{\boldsymbol{C}}_i}{\mathit{\boldsymbol{ \boldsymbol{\varPhi} }}_i}\left( {{t_{{{\rm{f}}},i}},t} \right)\mathit{\boldsymbol{x}}_i^{\rm{T}}}&{t \le {t_{{\rm{f}},i}}}\\ {{Z_i}\left( {{t_{{\rm{f}},i}}} \right)}&{t > {t_{{\rm{f}},i}}} \end{array}\;\;\;\;i \in \left\{ {1,2} \right\}} \right. $ (11)
 

式中:变换矩阵$ {\mathit{\boldsymbol{ \boldsymbol{\varPhi} }}_i}\left( {{t_{{\rm{f}}, i}}, t} \right)$

$ {\mathit{\boldsymbol{ \boldsymbol{\varPhi} }}_i}\left( {{t_{{{\rm{f}}},i}},t} \right) = \left[ {\begin{array}{*{20}{c}} 1&{{t_{{\rm{f}},i}} - t}\\ 0&1 \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} 1&{{t_{{\rm{go}},i}}}\\ 0&1 \end{array}} \right] $ (12)
 

其中:tgo, i为导弹至第i个点的剩余飞行时间(time-to-go),将式(12)代入式(11)可得

$ {Z_i} = \left\{ {\begin{array}{*{20}{l}} {{y_i} + {t_{{\rm{go}},i}}{v_i}}&{t \le {t_{{\rm{f}},i}}}\\ {{Z_i}\left( {{t_{{\rm{f}},i}}} \right)}&{t > {t_{{\rm{f}},i}}} \end{array}\;\;\;\;i \in \left\{ {1,2} \right\}} \right. $ (13)
 

得到Zi即零控脱靶量,它表示从当前时刻起,导弹不再执行任何控制指令的情况下到达第i个点的脱靶量。

对式(13)求导可得

$ {{\dot Z}_i} = \left\{ {\begin{array}{*{20}{l}} { - {t_{{\rm{go,}}i}}a}&{t \le {t_{{\rm{f}},i}}}\\ 0&{t > {t_{{\rm{f}},i}}} \end{array}\quad i \in \left\{ {1,2} \right\}} \right. $ (14)
 

经过以上状态变换,系统(6)由4阶降为2阶。

3.2 最优制导律设计

对于变换后的系统,最优控制问题1转化为在相同性能指标下,对系统(14)求最优解,即找到适当的a使Zi尽快减小到0,同时控制能量最小。

根据ZEM的物理意义,新的系统状态的终端约束为

$ {Z_i}\left( {{t_{{\rm{f}},i}}} \right) = 0\;\;\;\;i \in \left\{ {1,2} \right\} $ (15)
 

根据线性系统理论,式(14)可改写为

$ {Z_i}\left( {{t_{{\rm{f}},i}}} \right) = {Z_i}\left( t \right) = \int_t^{{t_{{\rm{f}},i}}} - \left( {{t_{{\rm{f,}}i}} - \tau } \right)a\left( \tau \right){\rm{d}}\tau \;\;\;t \le {t_{{\rm{f,}}i}} $ (16)
 

将式(15)代入式(16),可得

$ {Z_i}\left( t \right) = \int_t^{{t_{{\rm{f}},i}}} {\left( {{t_{{\rm{f,}}i}} - \tau } \right)a\left( \tau \right){\rm{d}}\tau } \;\;\;t \le {t_{{\rm{f,}}i}} $ (17)
 

在希尔伯特空间$ H=L^{2}\left[t, t_{\mathrm{f}, 2}\right]$中考虑此控制问题,定义内积${\rm{(}}f{\rm{ }}|{\rm{ }}g{\rm{) = }}{Z_i}(t){\rm{ = }}\int_t^{{t_{{\rm{f}}, i}}} f (\tau )g(\tau ){\rm{d}}\tau $。则此问题转化为在满足约束方程的所有元素中,求解最小范数元素,应用定理1[22]求解上述问题的解析解。

定理1   H为希尔伯特空间,{α1, α2, …, αn}为H中的一个线性无关集,c1, c2, …, cn为一组已知常数,在H中满足约束条件(x|α1)=c1, (x|α2)=c2,…,(x|αn)=cn的所有解构成一个集合V,此时V是一个余维数为n的线性流型,若xminV,并且具有最小范数,则有

$ {x_{\min }} = \sum\limits_{i = 1}^n {{b_i}{\alpha _i}} $

bi满足:

$ \left\{ {\begin{array}{*{20}{l}} {\left( {{\alpha _1}|{\alpha _1}} \right){b_1} + \left( {{\alpha _2}|{\alpha _1}} \right){b_2} + \cdots + \left( {{\alpha _n}|{\alpha _1}} \right){b_n} = {c_1}}\\ {\left( {{\alpha _1}|{\alpha _2}} \right){b_1} + \left( {{\alpha _2}|{\alpha _2}} \right){b_2} + \cdots + \left( {{\alpha _n}|{\alpha _2}} \right){b_n} = {c_2}}\\ {\;\; \vdots }\\ {\left( {{\alpha _1}|{\alpha _n}} \right){b_1} + \left( {{\alpha _2}|{\alpha _n}} \right){b_2} + \cdots + \left( {{\alpha _n}|{\alpha _n}} \right){b_n} = {c_n}} \end{array}} \right. $

对于最优控制问题1,式(17)中a(τ)与tf, iτ分别对应定理1中的xαi,由定理1可得,含有拉格朗日乘子λi的最优解形式为

$ a = \left\{ {\begin{array}{*{20}{l}} {{\lambda _1}\left( {{t_{{\rm{f}},1}} - t} \right) + {\lambda _2}\left( {{t_{{\rm{f}},2}} - t} \right)}&{t \le {t_{{\rm{f}},1}}}\\ {{\lambda _2}\left( {{t_{{\rm{f}},2}} - t} \right)}&{{t_{{\rm{f}},1}} < t \le {t_{{\rm{f}},2}}} \end{array}} \right. $ (18)
 

式中:拉格朗日乘子λi即定理1中的bi

ttf, 1为例,求解拉格朗日乘子的表达式。将式(18)代入式(17),可得

$ \begin{array}{l} {Z_1}\left( t \right) = {\lambda _1}\int_t^{{t_{1,1}}} {\left( {{t_{{\rm{f}},1}} - \tau } \right)} \left( {{t_{{\rm{f}},1}} - \tau } \right){\rm{d}}\tau + \\ \;\;\;\;\;\;\;\;\;\;\;{\lambda _2}\int_t^{{t_{{\rm{f}},1}}} {\left( {{t_{{\rm{f}},i}} - \tau } \right)} \left( {{t_{{\rm{f}},2}} - \tau } \right){\rm{d}}\tau = \\ \;\;\;\;\;\;\;\;\;\;\;{\lambda _1}\int_t^{{t_{{\rm{f}},1}}} {{{\left( {{t_{{\rm{f}},1}} - \tau } \right)}^2}} {\rm{d}}\tau + \\ \;\;\;\;\;\;\;\;\;\;\;{\lambda _2}\int_t^{{t_{{\rm{f}},1}}} {\left[ {\left( {{t_{{\rm{f}},2}} - {t_{{\rm{f}},1}}} \right)\left( {{t_{{\rm{f}},1}} - \tau } \right) + } \right.} \\ \;\;\;\;\;\;\;\;\;\;\;\left. {{{\left( {{t_{{\rm{f}},1}} - \tau } \right)}^2}} \right]{\rm{d}}\tau = {\lambda _1}\frac{1}{3}t_{{\rm{go}},1}^3 + {\lambda _2}\left[ {\frac{1}{2}\left( {{t_{{\rm{f}},2}} - } \right.} \right.\\ \;\;\;\;\;\;\;\;\;\;\;\left. {\left. {{t_{{\rm{f}},1}}} \right)t_{{\rm{go}},1}^2 + \frac{1}{3}t_{{\rm{go}},1}^3} \right] = {\lambda _1}\frac{1}{3}t_{{\rm{go}},1}^3 + \\ \;\;\;\;\;\;\;\;\;\;\;{\lambda _2}\left[ {\frac{1}{2}\left( {{t_{{\rm{go}},2}} - {t_{{\rm{go}},1}}} \right)t_{{\rm{go}},1}^2 + \frac{1}{3}t_{{\rm{go}},1}^3} \right] = \\ \;\;\;\;\;\;\;\;\;\;\;{\lambda _1}\frac{1}{3}t_{{\rm{go}},1}^3 + {\lambda _2}\left( {\frac{1}{2}{t_{{\rm{go}},2}}t_{{\rm{go}},1}^2 - \frac{1}{6}t_{{\rm{go}},1}^3} \right) \end{array} $ (19)
 

同理可以求得

$ {Z_2}\left( t \right) = {\lambda _2}\frac{1}{3}t_{{\rm{go}},2}^3 $ (20)
 

定义$ \mathit{\boldsymbol{\lambda }} = {\left[ {{\lambda _1}, {\lambda _2}} \right]^{\rm{T}}}, \mathit{\boldsymbol{Z}} = {\left[ {{Z_1}, {Z_2}} \right]^{\rm{T}}}$, 由式(19)和式(20)可得

$ \mathit{\boldsymbol{Z}} = \mathit{\boldsymbol{G\lambda }} $ (21)
 

式中:

$ \mathit{\boldsymbol{G}} = \left[ {\begin{array}{*{20}{c}} {\frac{1}{3}t_{{\rm{go}},1}^3}&{\frac{1}{2}{t_{{\rm{go}},2}}t_{{\rm{go}},1}^2 - \frac{1}{6}t_{{\rm{go}},1}^3}\\ {\frac{1}{2}{t_{{\rm{go}},2}}t_{{\rm{go}},1}^2 - \frac{1}{6}t_{{\rm{go}},1}^3}&{\frac{1}{3}t_{{\rm{go}},2}^3} \end{array}} \right] $ (22)
 

由式(21)可得

$ \mathit{\boldsymbol{\lambda }} = {\mathit{\boldsymbol{G}}^{ - 1}}\mathit{\boldsymbol{Z}} $ (23)
 

$t_{\mathrm{f}, 1} <t \leqslant t_{\mathrm{f}, 2} $时,Z1=0,ZEM向量降为一维,即Z=[Z2],此时$\boldsymbol{G}=\left[\frac{1}{3} t_{\mathrm{go}, 2}^{3}\right] $

将式(23)代入式(18),可最终得到

$ a = \left\{ {\begin{array}{*{20}{l}} {{{\left( {{\mathit{\boldsymbol{G}}^{ - 1}}\mathit{\boldsymbol{Z}}} \right)}^{\rm{T}}}{{\left[ {{t_{{\rm{go}},1}},{t_{{\rm{go}},2}}} \right]}^{\rm{T}}}}&{t \le {t_{{\rm{f}}.{\rm{i}}}}}\\ {\frac{{3{Z_2}}}{{t_{{\rm{go}},2}^2}}}&{{t_{{\rm{i}},1}} < t \le {t_{{\rm{f}},2}}} \end{array}} \right. $ (24)
 
3.3 讨论与分析

为了更好地理解上述推导得到的制导律,进行如下分析与讨论。

讨论1   用tf, i=t+tgo, i计算到达时间,当导弹位于路径点附近时,可以视为得到精确的到达时间。

讨论2   考虑到实际应用问题,对式(24)中的Z作如下处理:由于qiqR为小量,所以有

$ {q_i} - {q_{\rm{R}}} = \frac{{{y_i}}}{{{r_i}}} $ (25)
 

求导得

$ {{\dot q}_i} = \frac{{{v_i}{r_i} - {y_i}{{\dot r}_i}}}{{r_i^2}} = \frac{{{v_i}{t_{{\rm{go}},i}} + {y_i}}}{{r_i^2/V}} = \frac{{{Z_i}}}{{Vt_{{\rm{go}},i}^2}} $ (26)
 

由式(26)可得

$ {Z_i} = Vt_{{\rm{go}},i}^2{{\dot q}_i} $ (27)
 

在实际工程应用中,由弹目视线角速率、速度和弹目距离等可测物理量,可计算得到制导指令。

讨论3   当tf, 1 < ttf, 2时,a=3Z2/tgo, 22,与式(27)联立,可得

$ a = 3V{{\dot q}_2}\quad {t_{{\rm{f}},1}} < t \le {t_{{\rm{f}},2}} $ (28)
 

这与经典的比例导引制导律(Proportional Navigation Guidance, PNG)形式一致。根据文献[18],在只有一个目标点时,导航比为3的比例导引为最优制导律,与式(28)一致。进一步可以得出结论,仅在相邻两点之间连续采用比例导引并不是能量最优的,因为比例导引只保证当前段弹道能量最优,没有终端落角约束,因此总的控制能量消耗与本文提出的最优制导律不同。

讨论4   将ttf, 1时刻的过载指令表达式展开得

$ \begin{array}{l} a = \frac{{6\left( {2t_{{\rm{go}},2}^2 - {t_{{\rm{go}},1}}{t_{{\rm{go}},2}}} \right)}}{{t_{{\rm{go}},1}^2\left( {{t_{{\rm{go}},2}} - {t_{{\rm{go}},1}}} \right)\left( {4{t_{{\rm{go}},2}} - {t_{{\rm{go}},1}}} \right)}}{Z_1} - \\ \;\;\;\;\;\frac{{ - 6}}{{\left( {{t_{{\rm{go}},2}} - {t_{{\rm{go}},1}}} \right)\left( {4{t_{{\rm{go}},2}} - {t_{{\rm{go}},1}}} \right)}}{Z_2} \end{array} $ (29)
 

将式(27)代入式(29),可得

$ a = {N_1}V{{\dot q}_1} - {N_2}V{{\dot q}_2} $ (30)
 

式中:

$ \left\{ \begin{array}{l} {N_1} = \frac{{6{t_{{\rm{go}},2}}\left( {2{t_{{\rm{go}},2}} - {t_{{\rm{go}},1}}} \right)}}{{\left( {{t_{{\rm{go}},2}} - {t_{{\rm{go}},1}}} \right)\left( {4{t_{{\rm{go}},2}} - {t_{{\rm{go}},1}}} \right)}}\\ {N_2} = \frac{{6t_{{\rm{go}},2}^2}}{{\left( {{t_{{\rm{go}},2}} - {t_{{\rm{go}},1}}} \right)\left( {4{t_{{\rm{go}},2}} - {t_{{\rm{go}},1}}} \right)}} \end{array} \right. $ (31)
 

式(30)与含有偏置项的比例导引(BPN)形式一致,将其看作导航比随时间变化的BPN,根据文献[14]有

$ {{\dot q}_1} = {{\dot q}_0}{\left( {\frac{{{t_{{\rm{go}},1}}}}{{{t_{{\rm{f}},1}}}}} \right)^{{N_1} - 2}} + \frac{{{N_2}{{\dot q}_2}}}{{{N_1} - 2}}\left[ {1 - {{\left( {\frac{{{t_{{\rm{go}},1}}}}{{{t_{{\rm{f}},1}}}}} \right)}^{{N_1} - 2}}} \right] $ (32)
 

式中:q0为与初始条件有关的常数。

由式(32)可知$\mathop {\lim }\limits_{{t_{{\rm{go}}.1}} \to {0^ + }} {N_1} = 3 > 2 $,因此,$\lim \limits_{t_{\mathrm{go}, 1} \rightarrow 0^{+}}\left(\frac{t_{\mathrm{go}, 1}}{t_{\mathrm{f}, 1}}\right)^{N_{1}-2}=0 $,当导弹接近虚拟交班点,即路径点1时,

$ \mathop {\lim }\limits_{{t_{{\rm{go}},1}} \to {0^ + }} {{\dot q}_1} = \frac{{{N_2}}}{{{N_1} - 2}}{{\dot q}_2} $ (33)
 

此时,制导指令为

$ \begin{array}{l} \mathop {\lim }\limits_{{t_{{\rm{go}},1}} \to {0^ + }} a = \frac{{{N_1}{N_2}}}{{{N_1} - 2}}V{{\dot q}_2} - {N_2}V{{\dot q}_2} = \frac{{2{N_2}}}{{{N_1} - 2}}V{{\dot q}_2} = \\ \;\;\;\;\;\;\;\;\;\;\;\frac{{6t_{{\rm{go}},2}^2}}{{2t_{{\rm{go}},2}^2 - t_{{\rm{go}},1}^2 + 2{t_{{\rm{go}},1}}{t_{{\rm{go}},2}}}}V{{\dot q}_2} \end{array} $ (34)
 

由于

$ \mathop {\lim }\limits_{{t_{{\rm{go}},1}} \to {0^ + }} \frac{{6t_{{\rm{go}},2}^2}}{{2t_{{\rm{go}},2}^2 - t_{{\rm{go}},1}^2 + 2{t_{{\rm{go}},1}}{t_{{\rm{go}},2}}}} = 3 $ (35)
 

可得

$ \mathop {\lim }\limits_{{t_{{\rm{go}},1}} \to {0^ + }} a = 3V{{\dot q}_2} $ (36)
 

与式(28)一致,因此,在导弹经过虚拟交班点时,式(30)中的偏置项$ N_{2} V \dot{q}_{2}$具有抑制制导指令跳变的作用,使得式(24)得到的制导指令连续。

4 考虑终端落角约束的过定点最优制导律设计

在一些情况下,需要导弹以一定的终端落角攻击目标,以达到更好的毁伤效果。本节在第3节的基础上加入终端落角约束,对有落角约束的过定点最优制导律进行设计。

4.1 最优制导律设计

经过系统降阶后,第2节中的最优控制问题2转化为,对系统(14)求最优解,使目标函数达到极小值,并在满足约束(15)的同时,还满足以下终端落角约束:

$ \theta \left( {{t_{{\rm{f}},2}}} \right) = {\theta _d} $ (37)
 

由线性系统理论可得

$ \left\{ {\begin{array}{*{20}{l}} {{Z_i}\left( t \right) = \int_t^{{t_{{\rm{f}},i}}} {\left( {{t_{{\rm{f}},i}} - \tau } \right)a\left( \tau \right){\rm{d}}\tau } }&{t \le {t_{{\rm{f}},i}}}\\ {{\theta _{\rm{d}}} - \theta \left( t \right) = \int_t^{{t_{{\rm{f}},2}}} {\frac{{a\left( \tau \right)}}{V}{\rm{d}}\tau } }&{t \le {t_{{\rm{f}},2}}} \end{array}} \right. $ (38)
 

过载指令由两部分组成:

$ a = {a_Z} + {a_\theta } $ (39)
 

式中:aZaθ分别为与ZEM和弹道倾角相关的调节项。根据定理1,有

$ \left\{ \begin{array}{l} {a_Z} = \left\{ {\begin{array}{*{20}{l}} {{\lambda _1}\left( {{t_{{\rm{f}},1}} - t} \right) + {\lambda _2}\left( {{t_{{\rm{f}},2}} - t} \right)}&{t \le {t_{{\rm{f}},1}}}\\ {{\lambda _2}\left( {{t_{{\rm{f}},2}} - t} \right)}&{{t_{{\rm{f}},1}} < t \le {t_{\rm{f}}}} \end{array}} \right.\\ {a_\theta } = {\lambda _3}/V\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;t \le {t_{{\rm{f}},2}} \end{array} \right. $ (40)
 

同样地,以ttf, 1为例求解拉格朗日乘子,将式(39)代入式(38)可得

$ \left\{ \begin{array}{l} {Z_1}(t) = {\lambda _1}\int_t^{{t_{{\rm{f}},1}}} {\left( {{t_{{\rm{f}},1}} - \tau } \right)} \left( {{t_{{\rm{f}},1}} - \tau } \right){\rm{d}}\tau + \\ \;\;\;\;\;\;\;\;\;\;{\lambda _2}\int_t^{{t_{{\rm{f}},1}}} {\left( {{t_{{\rm{f}},1}} - \tau } \right)} \left( {{t_{{\rm{f}},2}} - \tau } \right){\rm{d}}\tau + \\ \;\;\;\;\;\;\;\;\;\frac{{{\lambda _3}}}{V}\int_t^{{t_{{\rm{f}},1}}} {\left( {{t_{{\rm{f}},1}} - \tau } \right){\rm{d}}\tau } \\ {Z_2}(t) = {\lambda _1}\int_t^{{t_{{\rm{f}},1}}} {\left( {{t_{{\rm{f}},2}} - \tau } \right)} \left( {{t_{{\rm{f}},1}} - \tau } \right){\rm{d}}\tau + \\ \;\;\;\;\;\;\;\;\;\;\;{\lambda _2}\int_t^{{t_{{\rm{f,2}}}}} {\left( {{t_{{\rm{f,2}}}} - \tau } \right)} \left( {{t_{{\rm{f,2}}}} - \tau } \right){\rm{d}}\tau + \\ \;\;\;\;\;\;\;\;\;\;\frac{{{\lambda _3}}}{V}\int_t^{{t_{{\rm{f}},2}}} {\left( {{t_{{\rm{f}},2}} - \tau } \right){\rm{d}}\tau } \\ {\theta _{\rm{d}}} - \theta (t) = \frac{{{\lambda _1}}}{V}\int_t^{{t_{{\rm{f}},1}}} {\left( {{t_{{\rm{f}},1}} - \tau } \right){\rm{d}}\tau } + \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\frac{{{\lambda _2}}}{V}\int_t^{{t_{{\rm{f}},2}}} {\left( {{t_{{\rm{f}},2}} - \tau } \right){\rm{d}}\tau } + \frac{{{\lambda _3}}}{{{V^2}}}\int_t^{{t_{{\rm{f}},2}}} {\rm{d}} \tau \end{array} \right. $ (41)
 

由式(41)可得

$ \left[ {\begin{array}{*{20}{c}} \mathit{\boldsymbol{Z}}\\ {{\theta _{\rm{d}}} - \theta } \end{array}} \right] = \mathit{\boldsymbol{G\lambda }},\mathit{\boldsymbol{G}} = \left[ {\begin{array}{*{20}{c}} {{\mathit{\boldsymbol{G}}_1}}&{{\mathit{\boldsymbol{G}}_{12}}}\\ {{\mathit{\boldsymbol{G}}_{21}}}&{{G_2}} \end{array}} \right] $ (42)
 

式中:$ \boldsymbol{G}_{1} \in {\bf R}^{2 \times 2}, \boldsymbol{G}_{12} \in {\bf R}^{2 \times 1}, \boldsymbol{G}_{21} \in {\bf R}^{1 \times 2}$,由3.2节可知

$ {\mathit{\boldsymbol{G}}_1} = \left[ {\begin{array}{*{20}{c}} {\frac{1}{3}t_{{\rm{go}},1}^3}&{\frac{1}{2}{t_{{\rm{go}},2}}t_{{\rm{go}},1}^2 - \frac{1}{6}t_{{\rm{go,1}}}^3}\\ {\frac{1}{2}{t_{{\rm{go}},2}}t_{{\rm{go}},1}^2 - \frac{1}{6}t_{{\rm{go}},1}^3}&{\frac{1}{3}t_{{\rm{go}},2}^3} \end{array}} \right] $ (43)
 

下面求解G12G21G2:

$ \left\{ {\begin{array}{*{20}{l}} {\int_t^{{t_{{\rm{f}},i}}} {\left( {{t_{{\rm{f}},i}} - \tau } \right)} {\rm{d}}\tau = - \left. {\frac{1}{2}{{\left( {{t_{{\rm{f}},i}} - \tau } \right)}^2}} \right|_t^{{t_{{\rm{f}},i}}} = \frac{{t_{{\rm{go}},i}^2}}{2}}\\ {\int_t^{{t_{{\rm{f}},2}}} {\rm{d}} \tau = {t_{{\rm{go}},2}}} \end{array}} \right. $ (44)
 

可得

$ \left\{ {\begin{array}{*{20}{l}} {{\mathit{\boldsymbol{G}}_{12}} = \mathit{\boldsymbol{G}}_{21}^{\rm{T}} = {{\left[ {\frac{{t_{{\rm{go}},1}^2}}{{2V}},\frac{{t_{{\rm{go}},2}^2}}{{2V}}} \right]}^{\rm{T}}}}\\ {{G_2} = \frac{{{t_{{\rm{go}}.2}}}}{{{V^2}}}} \end{array}} \right. $ (45)
 

由式(42)~式(45)可以解得拉格朗日乘子为

$ \mathit{\boldsymbol{\lambda }} = {\mathit{\boldsymbol{G}}^{ - 1}}\left[ {\begin{array}{*{20}{c}} \mathit{\boldsymbol{Z}}\\ {{\theta _{\rm{d}}} - \theta } \end{array}} \right] $ (46)
 

可得过载指令为

$ \begin{array}{l} a = {\mathit{\boldsymbol{\lambda }}^{\rm{T}}}{\left[ {{t_{{\rm{go}},1}},{t_{{\rm{go}},2}},\frac{1}{V}} \right]^{\rm{T}}} = \\ \;\;\;\;\;{\left( {{\mathit{\boldsymbol{G}}^{ - 1}}\left[ {\begin{array}{*{20}{c}} \mathit{\boldsymbol{Z}}\\ {{\theta _{\rm{d}}} - \theta } \end{array}} \right]} \right)^{\rm{T}}}{\left[ {{t_{{\rm{go}},1}},{t_{{\rm{go,}}2}},\frac{1}{V}} \right]^{\rm{T}}} \end{array} $ (47)
 

tf, 1 < ttf, 2时,Z1=0,ZEM向量降为一维,即Z=[Z2],此时

$ \mathit{\boldsymbol{G}} = \left[ {\begin{array}{*{20}{c}} {\frac{{t_{{\rm{go}},2}^3}}{3}}&{\frac{{t_{{\rm{go}},2}^2}}{{2V}}}\\ {\frac{{t_{{\rm{go}},2}^2}}{{2V}}}&{\frac{{{t_{{\rm{go}},2}}}}{{{V^2}}}} \end{array}} \right] $ (48)
 

过载指令为

$ \begin{array}{l} a = {\mathit{\boldsymbol{\lambda }}^{\rm{T}}}{\left[ {{t_{{\rm{go}},2}},\frac{1}{V}} \right]^{\rm{T}}} = \\ \;\;\;\;\;\;{\left( {{\mathit{\boldsymbol{G}}^{ - 1}}\left[ {\begin{array}{*{20}{c}} {{Z_2}}\\ {{\theta _{\rm{d}}} - \theta } \end{array}} \right]} \right)^{\rm{T}}}{\left[ {{t_{{\rm{go}},2}},\frac{1}{V}} \right]^{\rm{T}}} \end{array} $ (49)
 
4.2 讨论与分析

注意到当tf, 1 < ttf, 2时,模型转化为针对一个固定目标,考虑终端落角约束的最优制导问题,由式(49)可得

$ a = {\lambda _2}{t_{{\rm{go,}}2}} + \frac{{{\lambda _3}}}{V} = \frac{{6{Z_2}}}{{t_{{\rm{go,}}2}^2}} - \frac{{2V}}{{{t_{{\rm{go,}}2}}}}\left( {{\theta _{\rm{d}}} - \theta } \right) $ (50)
 

对式(1)的第3项积分可得

$ {\theta _{\rm{d}}} - \theta = \frac{1}{V}\int_t^{{t_{\rm{f}}}} {a{\rm{d}}\tau } $ (51)
 

将式(5)代入式(51),可得

$ \begin{array}{l} {\theta _{\rm{d}}} - \theta = \frac{1}{V}\int_t^{{t_{\rm{f}}}} {a{\rm{d}}\tau } = \frac{1}{V}\int_t^{{t_{\rm{f}}}} - {{\dot v}_2}{\rm{d}}\tau = \\ \;\;\;\;\;\;\;\;\frac{1}{V}\left( {{v_2} - {v_{2{\rm{f}}}}} \right) = \frac{1}{V}\left( {{{\dot y}_2} - {{\dot y}_{2{\rm{f}}}}} \right) \end{array} $ (52)
 

又由式(13)可得

$ {Z_2} = {y_2} + {t_{{\rm{go,}}2}}{v_2} = {y_2} + {{\dot y}_2}{t_{{\rm{go}},2}} $ (53)
 

将式(52)、式(53)代入式(50),可得

$ \begin{array}{l} a = \frac{6}{{t_{{\rm{go}},2}^2}}\left( {{y_2} + {{\dot y}_2}{t_{{\rm{go}},2}}} \right) - \frac{{2V}}{{{t_{{\rm{go}},2}}}} \cdot \frac{1}{V}\left( {{{\dot y}_2} - {{\dot y}_{2{\rm{f}}}}} \right) = \\ \;\;\;\;\frac{{6{y_2} + 4{{\dot y}_2}{t_{{\rm{go}},2}} + 2{{\dot y}_{{\rm{2f}}}}{t_{{\rm{go}},2}}}}{{t_{{\rm{go}},2}^2}} \end{array} $ (54)
 

与经典的弹道成型制导律(Trajectory Shaping Guidance, TSG)形式一致。在只有一个目标点,考虑终端落角约束的制导问题中,TSG是能量最优的[18]

5 仿真验证

本节以非线性系统为研究对象,对提出的制导律进行仿真验证,并与经典最优制导律进行对比。飞行弹道包括无控段、中制导段与末制导段,导弹在40 s启控后进入中制导段,经过虚拟交班点进入末制导段,仿真参数如表 1所示。

表 1 仿真参数 Table 1 Simulation Index
参数 数值
导弹初始坐标/m (0,0)
启控时间/s 40
虚拟交班点/m (22000,5000)
目标点/m (25000,0)
速度/(m·s-1) 200
初始弹道倾角/(°) 40
5.1 无落角约束的最优制导律仿真

本节针对过虚拟交班点的无落角约束最优制导律(24)进行仿真验证,并在中制导与末制导段分别采用最优比例导引,即导航比为3的比例导引律,与本文提出的最优制导律进行对比。仿真结果如图 2所示。

图 2 无落角约束情况下的仿真结果对比 Fig. 2 Comparison of simulation results without angle constraint

图 2分别给出了2种制导律下的弹道、弹道倾角、过载指令和控制能量曲线。图 2(a)中五角星表示需要经过的路径点1(虚拟交班点位置)与路径点2(目标),可看出两种制导律都可以准确经过目标路径点,但通过对比2条弹道曲线末制导段发现,本文提出的最优制导律比最优比例导引制导律弹道更为平直;对比图 2(c),该最优制导律在末制导阶段的过载指令远小于最优比例导引。且最优比例导引的过载指令在虚拟交班点处出现跳变,本文提出的制导律下过载指令为连续的,最大需用过载不超过分段比例导引的36%。图 2(d)为2种制导律下的控制能量消耗,可以看出,本文提出的最优制导律能量消耗大大减少,比最优比例导引所需的能量消耗减少49%。

5.2 考虑终端落角约束的最优制导律仿真

本节对考虑终端落角约束的最优制导律进行仿真验证。导弹期望落角选取为-90°。采用PNG与TSG切换的制导律进行对比,即导弹在启控后采用PNG,经过虚拟交班点进入末制导阶段后,由PNG切换为TSG以实现期望落角。仿真结果对比如图 3所示。

图 3 考虑终端落角约束情况下的仿真结果对比 Fig. 3 Comparison of simulation results with angle constraint

对比发现,图 3(a)图 2(a)第一段弹道相同,而最优制导律在加入终端落角约束后,末制导段弹道轨迹与弹道倾角变化较为明显,因为本文提出的最优制导律在开始时就将落角约束考虑进去,从而得到全局最优的制导律。图 3(c)给出了两种制导律下的过载指令,可以看到在经过虚拟交班点时,PNG+TSG下过载指令与无落角约束时一样出现了不连续性,而本文提出的最优制导律为连续的。图 3(d)为2种制导律所消耗的控制能量,对比发现本文提出的制导律所需控制能量比无落角约束时明显增大,但仍比PNG+TSG下能量消耗减少22%。

6 结论

1) 本文在希尔伯特空间基于最优化理论,提出了一种针对固定目标,过虚拟交班点的全局最优制导律。在无落角约束下,能准确经过虚拟交班点,到达固定目标点,且全局所需控制能量比最优比例导引制导律减少了49%。

2) 在最优制导律的基础上考虑终端落角约束,推导出含终端落角约束的全局最优制导律,该制导律可以实现期望落角,且全局所需控制能量比PNG+TSG下减少了22%。

3) 该制导律形式简单,计算制导指令所需的信息较少,具有良好的工程应用价值。

参考文献
[1] 刘海军, 王丽娜. 复合制导防空导弹中末制导交班问题研究[J]. 现代防御技术, 2006, 34(2): 29-33.
LIU H J, WANG L N. Study on handover problem of compound guidance missile weapon[J]. Modern Defense Technology, 2006, 34(2): 29-33. (in Chinese)
Cited By in Cnki (8) | Click to display the text
[2] 姚郁, 郑天宇, 贺风华, 等. 飞行器末制导中的几个热点问题与挑战[J]. 航空学报, 2015, 36(8): 2696-2716.
YAO Y, ZHENG T Y, HE F H. Several hot issues and challenges in terminal guidance of flight vehicles[J]. Acta Aeronautica et Astronautica Sinica, 2015, 36(8): 2696-2716. (in Chinese)
Cited By in Cnki | Click to display the text
[3] 唐明南, 郑学合, 王天舒, 等. 多平台信息下机动目标拦截中末制导交班研究[J]. 系统工程与电子技术, 2017, 39(12): 2757-2764.
TANG M N, ZHENG X H, WANG T S, et al. Handover of guidance from midcourse to terminal in the maneuvering target interception based on multi-platform information[J]. Systems Engineering and Electronics, 2017, 39(12): 2757-2764. (in Chinese)
Cited By in Cnki | Click to display the text
[4] SHEHAB S, RODRIGUES L. Preliminary results on UAV path following using piecewise-affine control[C]//IEEE Conference on Control Applications. Piscataway, NJ: IEEE Press, 2005.
[5] YU M J, ATKINS E M, KOLMANOVSKY I, et al. Trim-commanded adaptive control for waypoint-defined trajectory following[C]//AIAA Guidance, Navigation, and Control (GNC) Conference. Reston, VA: AIAA, 2013.
[6] 张友安, 梁勇, 刘京茂, 等. 基于轨迹成型的攻击角度与时间控制[J]. 航空学报, 2018, 39(9): 322009.
ZHANG Y A, LIANG Y, LIU J M, et al. Trajectory reshaping based impact angle and time control[J]. Acta Aeronautica et Astronautica Sinica, 2018, 39(9): 322009. (in Chinese)
Cited By in Cnki | Click to display the text
[7] WHANG I H, HWANG T W. Horizontal waypoint guid ance design using optimal control[J]. IEEE Transactions on Aerospace and Electronic Systems, 2002, 38(3): 1116-1120.
Click to display the text
[8] MOON G, KIM Y. Optimum flight path design passing through waypoints for autonomous flight control system[J]. Engineering Optimization, 2013, 37(7): 755-774.
Click to display the text
[9] RYOO C K, SHIN H S, TAHK M J. Energy optimal waypoint guidance synthesis for antiship missiles[J]. IEEE Transactions on Aerospace & Electronic Systems, 2010, 46(1): 80-95.
[10] 孙胜, 张华明, 周荻. 考虑自动驾驶仪动特性的终端角度约束滑模导引律[J]. 宇航学报, 2013, 34(1): 69-78.
SUN S, ZHANG H M, ZHOU D. Sliding mode guidance law with autopilot lag for terminal angle constrained trajectories[J]. Journal of Astronautics, 2013, 34(1): 69-78. (in Chinese)
Cited By in Cnki (30) | Click to display the text
[11] RYOO C K, CHO H, TAHK M J. Optimal guidance laws with terminal impact angle constraint[J]. Journal of Guidance, Control, and Dynamics, 2005, 28(4): 724-732.
Click to display the text
[12] RYOO C K, CHO H, TAHK M J. Time-to-go weighted optimal guidance with impact angle constraints[J]. IEEE Transactions on Control Systems Technology, 2006, 14(3): 483-492.
Click to display the text
[13] 张友安, 黄诘, 孙阳平. 带有落角约束的一般加权最优制导律[J]. 航空学报, 2014, 35(3): 848-856.
ZHANG Y A, HUANG J, SUN Y P. Generalized weighted optimal guidance laws with impact angle constraints[J]. Acta Aeronautica et Astronautica Sinica, 2014, 35(3): 848-856. (in Chinese)
Cited By in Cnki | Click to display the text
[14] 闫梁, 赵继广, 沈怀荣, 等. 带末端碰撞角约束的三维联合偏置比例制导律设计[J]. 航空学报, 2014, 35(7): 1999-2010.
YAN L, ZHAO J G, SHEN H R, et al. Three-dimensional united biased proportional navigation guidance law for interception of targets with angular constraints[J]. Acta Aeronautica et Astronautica Sinica, 2014, 35(7): 1999-2010. (in Chinese)
Cited By in Cnki | Click to display the text
[15] ERER K S, MERTTOPCUOGLU O. Indirect impact-angle-control against stationary targets using biased pure proportional navigation[J]. Journal of Guidance, Control, and Dynamics, 2012, 35(2): 700-704.
Click to display the text
[16] 郭建国, 韩拓, 周军, 等. 基于终端角度约束的二阶滑模制导律设计[J]. 航空学报, 2017, 38(2): 210-219.
GUO J G, HAN T, ZHOU J, et al. Second-order sliding-mode guidance law with impact angle constraint[J]. Acta Aeronautica et Astronautica Sinica, 2017, 38(2): 210-219. (in Chinese)
Cited By in Cnki | Click to display the text
[17] JEON I S, LEE J I. Optimality of proportional navigation based on nonlinear formulation[J]. IEEE Transactions on Aerospace and Electronic Systems, 2010, 46(4): 2051-2055.
Click to display the text
[18] ZARCHAN P. Tactical and strategic missile guidance[M]. Reston, VA: AIAA, 2012: 569-577.
[19] SHAFERMAN V, SHIMA T. Linear quadratic guidance laws for imposing a terminal intercept angle[J]. Journal of Guidance, Control, and Dynamics, 2008, 31(5): 1400-1412.
Click to display the text
[20] RATNOO A, HAYOUN S Y, GRANOT A, et al. Path following using trajectory shaping guidance[C]//AIAA Guidance, Navigation, and Control (GNC) Conference. Reston, VA: AIAA, 2013.
[21] HE S, LEE C-H, SHIN H-S. Minimum-effort waypoint-following guidance[J/OL]. Journal of Guidance, Control, and Dynamics, (2019-01-30)[2019-06-27]. http://dspace.lib.cranfield.ac.uk/handle/1826/13904.
[22] 王日爽. 泛函分析与最优化理论[M]. 北京: 北京航空航天大学出版社, 2003: 98-105.
WANG R S. Functional analysis and optimization theory[M]. Beijing: Beihang University Press, 2003: 98-105. (in Chinese)
http://dx.doi.org/10.7527/S1000-6893.2019.23249
中国航空学会和北京航空航天大学主办。
0

文章信息

李晨迪, 王江, 李斌, 何绍溟, 张彤
LI Chendi, WANG Jiang, LI Bin, HE Shaoming, ZHANG Tong
过虚拟交班点的能量最优制导律
Energy-optimal guidance law with virtual hand-over point
航空学报, 2019, 40(12): 323249.
Acta Aeronautica et Astronautica Sinica, 2019, 40(12): 323249.
http://dx.doi.org/10.7527/S1000-6893.2019.23249

文章历史

收稿日期: 2019-06-27
退修日期: 2019-08-01
录用日期: 2019-09-06
网络出版时间: 2019-09-17 10:02

相关文章

工作空间