CN106352738A - 基于输出一致性的多弹协同制导方法 - Google Patents
基于输出一致性的多弹协同制导方法 Download PDFInfo
- Publication number
- CN106352738A CN106352738A CN201610402960.XA CN201610402960A CN106352738A CN 106352738 A CN106352738 A CN 106352738A CN 201610402960 A CN201610402960 A CN 201610402960A CN 106352738 A CN106352738 A CN 106352738A
- Authority
- CN
- China
- Prior art keywords
- guided missile
- piece
- variable
- matrix
- missile
- 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
Classifications
-
- 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
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)
Abstract
本发明公开了一种基于输出一致性的多弹协同制导方法,该方法为首先给定各导弹法向加速度和视场角的约束条件;基于导弹和目标的相对运动关系,建立导弹状态依赖的线性模型,选择所述导弹状态依赖的线性模型的趋同变量、视线角速度和视场角作为输出变量,根据设定的预测时域建立预测控制模型;基于预测控制模型和各导弹打击时间协同的要求,获得性能指标函数;其次将性能指标函数和约束条件转化为二次规划标准型;获得多弹协同制导的优化模型,利用凸优化方法,获得当前时刻各导弹的法向加速度,最后下一时刻采用更新后的各导弹状态直至各导弹命中目标。本发明能够实现在法向过载及视场角受到限制情况下的最优协同制导。
Description
技术领域
本发明属于制导技术领域,具体涉及基于输出一致性的多弹协同制导方法。
背景技术
随着导弹拦截***的不断发展,依靠单弹突防、实现打击的战术逐渐受到限制。采用多枚导弹对同一目标实现饱和攻击,摧毁具有一定拦截能力目标的方法,已经成为了保证导弹战场生存能力、提高打击任务成功率的主要手段。从该思路出发,设计一种驱动多枚导弹同时命中目标的协同制导方法,对于突防技术的进步具有重要意义。
现有的保证同时打击的协同制导方法,存在以下问题有待解决:
1)单纯针对修正打击时间设计协同的制导方法难以保证精确命中目标,导引末端需要切换导引律以保证打击精度,而舍弃一定打击时间控制精度。此外,切换导引律导致***复杂程度提高,出现非连续情况,不利于实施。
2)协同制导要求各导弹同时命中目标,而导弹速度大小难以通过舵面控制,因此对于距离目标较近的导弹,需要绕行机动,延长命中时间。此举易导致弹道过于弯曲,在制导律设计时若不考虑视场角限制,则难以保证导引头视场始终捕获目标。
3)由于导弹舵面偏转有限,导弹的法向加速度受到限制,导致导弹控制***的输入具有饱和非线性。
4)针对单一导弹设计命中时间控制制导律,而后指定命中时间的协同制导方法,未将导弹群作为整体考虑,无法实现最优的整体表现。
综上所述,现有技术主要针对单一导弹设计命中时间控制制导律,而非将协 同制导问题作为集群问题解决,且未考虑各种物理条件限制。
发明内容
有鉴于此,本发明提供了基于输出一致性的多弹协同制导方法;本发明能够实现在法向过载及视场角受到限制情况下的最优协同制导。
实现本发明的具体实施方案如下:
基于输出一致性的多弹协同制导方法,具体步骤如下:
步骤一、给定各导弹法向加速度和视场角的约束条件;
步骤二、基于导弹和目标的相对运动关系,建立导弹状态依赖的线性模型,采用数值方法对所述导弹状态依赖的线性模型进行离散化,基于离散化后的导弹状态依赖的线性模型,根据设定的预测时域建立预测控制模型;
步骤三、基于预测控制模型和各导弹打击时间协同的要求,获得性能指标函数;
步骤四、基于二次规划形式,将性能指标函数和约束条件转化为二次规划标准型;
步骤五、基于当前时刻各导弹状态、二次规划标准型的各导弹法向加速度和视场角的约束条件和性能指标函数,获得多弹协同制导的优化模型,基于多弹协同制导的优化模型,利用凸优化方法,获得当前时刻各导弹的最优控制序列;
步骤六、从最优控制序列中选择各导弹的法向加速度并代入导弹状态依赖的线性模型,计算下一时刻的各导弹状态,按照步骤五的方式重复直至各导弹命中目标。
进一步地,其特征在于,步骤一具体实现过程如下:
定义编队中第i枚导弹的法向加速度约束条件为其中a(i)为第i枚导弹的法向加速度,为第i枚导弹的最大法向加速度;第i枚导弹的视场角的约束条件为σ(i)为第i枚导弹的速度前置角,为第i枚导弹的速度最大前置角。
进一步地,步骤二具体实现过程如下:
选取第i枚导弹的状态变量 为第i枚导弹状态变量的变化率,第i枚导弹的控制变量u(i)=a(i),其中i,j=1,...,Nv,i≠j,根据导弹和目标的相对运动关系,建立导弹状态依赖的线性模型:
其中a(i)表示第i枚导弹的法向加速度;A(i)为第i枚导弹的状态矩阵、B(i)为第i枚导弹的输入矩阵、y(i)、和分别为第i枚导弹的一致性变量、视线角速度和视场角,分别为第i枚导弹的一致性变量、视线角速度和视场角的输出矩阵,分别为r(i),λ(i),σ(i)对应的赋予变量,其表示形式分别如下:
r(i)表示第i枚导弹的弹目距离,表示第i枚导弹的弹目相对速度,λ(i)表示第i枚导弹的视线角,表示第i枚导弹的视线角速度,σ(i)为第i枚导弹的速度前置角,V表示导弹飞行速度,kr和分别为弹目距离的权重和弹目相对速度的权重;
采用数值方法对状态矩阵A(i)、输入矩阵B(i)进行离散化,获得第i枚导弹的导弹状态依赖的线性模型的离散化形式
定义k时刻导弹群整体状态变量X(k),k时刻导弹群整体的一致性变量、视线角速度和视场角分别为Y(k)、Z1(k)、Z2(k),控制变量U(k),X(k)、Y(k)、Z1(k)、Z2(k)和U(k)分别为
其中,x(i)(k)为k时刻第i枚导弹状态变量;y(i)(k)为k时刻第i枚导弹的一致性变量;为k时刻第i枚导弹的视线角速度;为k时刻第i枚导弹的视场角;u(i)(k)为k时刻第i枚导弹的控制变量,i=1,...,Nv;
定义当前时刻导弹群整体状态矩阵当前时刻导弹群整体状态的一致性变量、视线角速度和视场角的输出矩阵分别为控制矩阵 和分别为
其中blkdiag为分块对角矩阵运算;
定义预测控制模型的预测时域为N,获得预测时域内状态变量预测时域内导弹群整体状态的一致性变量、视线角速度和视场角的输出矩阵的输出变量分别为和控制变量状态矩阵控制矩阵预测时域内导弹群整体状态的一致性变量、视线角速度和视场角的输出矩阵的输出矩阵分别为 和和分别为
其中IM为M阶单位矩阵,M=N,为矩阵的克罗内克积运算;X(k+t)为k+t时刻导弹群整体的一致性变量,Y(k+t)为k+t时刻导弹群整体的视线角速度,Z1(k+t)为k+t时刻导弹群整体的视线角速度,Z2(k+t)为k+t时刻导弹群整体的视线视场角,U(k+t-1)为k+t-1时刻导弹群整体的控制变量;t=1,...,N;
选择所述导弹状态依赖的线性模型的一致性变量、视线角速度和视场角作为输出变量,基于预测时域内状态变量输出变量和控制变量状态矩阵控制矩阵输出矩阵和获得预测控制模型为
其中,X(k)为k时刻导弹群整体状态变量。
进一步地,所述数值方法为欧拉方法、改进欧拉方法、龙格-库塔法或线性多步法。
进一步地,所述凸优化方法为二次规划方法。
进一步地,步骤三具体实现过程如下:根据各导弹打击时间协同的要求,即将各导弹打击时间协同目标转化为各导弹弹目距离和弹目相对速度一致性目标、各导弹控制输入的优化和驱动各导弹命中目标,性能指标函数表达式为
J=J1+J2+J3 (3)
其中,J1为协同变量代价函数,J2为控制变量代价函数,J3为视线角速度代价函数,其分别为:
Im为m阶单位矩阵,u(i)(k)为k时刻第i枚导弹的控制变量,y(i)(k+1)为k+1时刻第i枚导弹的协同变量,为k+1时刻第i枚导弹的视线角速度,Nv为导弹总枚数,N为预测时域,R(i)、Q(i)为各导弹控制变量的权重和视线角速度的权重,为表示导弹通讯拓扑的邻接矩阵Aadj中的元素,
定义扩展拉普拉斯矩阵为其中L为表示导弹通讯拓扑的拉普拉斯矩阵,L=D-Aadj,D为入度矩阵,Im为m阶单位矩阵。
定义导弹群整体的权重,和分别为协控制变量和视线角速度的整体权重
基于扩展拉普拉斯矩阵和导弹群整体的权重,则预测时域上的拉普拉斯矩阵协控制变量权重和视线角速度权重分别为
其中,IN为N阶单位矩阵,为矩阵的克罗内克积运算;
基于预测控制模型,将式(4)带入式(3),获得性能指标函数J
其中为加权范数运算,
进一步地,其特征在于,步骤四具体实现过程如下:
基于性能指标函数,作如下变换:
其中H、F、constant分别为二次项系数、一次项系数和常数项,有
和均为预测时域内输出变量、为预测时域内控制变量,状态矩阵控制矩阵和分别为预测时域内输出矩阵;为加权范数运算,
针对各导弹法向加速度约束条件依据二次规划形式,则二次规划标准型为
其中
其中为N×Nv阶单位矩阵;1N为N×1列向量;为第i枚导弹的最大法向加速度;u(i)为第i枚导弹的控制变量;X(k)为k时刻导弹群整体的一致性变量;Umax为预测时域内控制变量;
针对各导弹视场角约束条件依据二次规划形式,则二次规划标准型为
其中
其中,σ(i)为第i枚导弹的速度前置角,为第i枚导弹的最大速度前置角,i=1,...,Nv。
有益效果:
(1)本发明基于预测控制模型和各导弹打击时间协同的要求,将各导弹打击时间协同目标转化为各导弹弹目距离和弹目相对速度一致性目标,即同时命中的协同制导问题转化为各导弹弹目距离与弹目相对速度一致性问题,避免了对剩余时间的估计以及制导律切换等问题,便于实现。
(2)本发明采用模型预测控制对导弹群整体建立模型,并设计控制律,可保证整体最优性。
(3)本发明基于二次规划的求解方式易于在控制指令生成时引入约束条件限制(如法向加速度、视场角),更贴近实际制导环境。
附图说明
图1为基于输出一致性的多弹协同制导方法流程图。
图2为弹目相对关系几何图。
图3为弹道轨迹图。
图4为法向加速度曲线图。
图5为视场角曲线图。
图6为视线角速度曲线图。
图7为剩余命中时间曲线图。
图8为最大剩余命中时间误差曲线图。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
本发明提供了基于输出一致性的多弹协同制导方法。其基本思想在于:将各导弹命中时间一致的要求在连续的控制指令形成过程中转化为各导弹弹目距离和各导弹弹目相对速度趋于一致的要求,从而达到各导弹命中时间相同的协同 制导控制目标。在制导过程中,本发明采用模型预测控制的方式形成控制指令,使得导弹群整体性能表现最优化,并且能够实现对法向过载和视场角的限制。
具体实现步骤如下:
步骤一:给定不同导弹法向加速度和视场角的约束条件;
如法向加速度、视场角等。对于不同类型的导弹,可给出不同的约束条件。定义编队中第i枚导弹的法向加速度约束条件为其中a(i)为第i枚导弹的法向加速度。设计视场角时,假设导弹侧滑角较小且可忽略,将视场角近似为导弹速度前置角σ,则给出第i枚导弹的视场角的约束条件为σ(i)为第i枚导弹的速度前置角。
步骤二:基于导弹和目标的相对运动关系,建立导弹状态依赖的线性模型,根据导弹状态依赖的线性模型,选择所述模型的一致性变量、视线角速度和视场角作为输出变量,采用数值方法对所述导弹状态依赖的线性模型进行离散化,基于离散化后的导弹状态依赖的线性模型,根据设定的预测时域建立预测控制模型。
二维平面内导弹群攻击目标的几何关系如图1所示,其中M(i)、M(j)、T分别表示第i枚导弹、第j枚导弹、攻击目标,上标(i)、(j)分别表示第i枚导弹、第j枚导弹对应的状态,i,j=1,...,Nv,i≠j,Nv为导弹总枚数,a表示导弹法向加速度。
选择所述导弹状态依赖的线性模型的一致性变量、视线角速度和视场角作为输出变量,选取第i枚导弹的状态变量 为第i枚导弹状态变量的变化率,第i枚导弹的控制变量u(i)=a(i),其中i,j=1,...,Nv,i≠j,根据导弹和目标的相对运动关系,建立导弹状态依赖的线性模型:
其中a(i)表示第i枚导弹的法向加速度;A(i)为第i枚导弹的状态矩阵、B(i)为第i枚导弹的输入矩阵、y(i)、和分别为第i枚导弹的一致性变量、视线角速度和视场角,分别为第i枚导弹的一致性变量、视线角速度和视场角的输出矩阵,分别为r(i),λ(i),σ(i)对应的赋予变量,其表示形式分别如下:
r(i)表示第i枚导弹的弹目距离,表示第i枚导弹的弹目相对速度,表示第i枚导弹的视线角速度,σ(i)为第i枚导弹的速度前置角,V表示导弹飞行速度,kr和分别为弹目距离的权重和弹目相对速度的权重;λ(i)表示第i枚导弹的视线角。
采用数值方法对状态矩阵A(i)、输入矩阵B(i)进行离散化,获得第i枚导弹的导弹状态依赖的线性模型的离散化形式
定义k时刻导弹群整体状态变量X(k),k时刻导弹群整体的一致性变量、视线角速度和视场角分别为Y(k)、Z1(k)、Z2(k),控制变量U(k),X(k)、Y(k)、Z1(k)、Z2(k)和U(k)分别为
其中,x(i)(k)为k时刻第i枚导弹状态变量;y(i)(k)为k时刻第i枚导弹的一致性变量;为k时刻第i枚导弹的视线角速度;为k时刻第i枚导弹的视场角;u(i)(k)为k时刻第i枚导弹的控制变量,i=1,...,Nv;
定义当前时刻导弹群整体状态矩阵当前时刻导弹群整体状态的一致性变量、视线角速度和视场角的输出矩阵分别为控制矩阵 和分别为
其中blkdiag为分块对角矩阵运算。
定义预测控制模型的预测时域为N,获得预测时域内状态变量预测时域内导弹群整体状态的一致性变量、视线角速度和视场角的输出矩阵的输出变量分别为和控制变量状态矩阵控制矩阵预测时域内导弹群 整体状态的一致性变量、视线角速度和视场角的输出矩阵的输出矩阵和和分别为
其中IM为M阶单位矩阵,M=N,为矩阵的克罗内克积运算;X(k+t)为k+t时刻导弹群整体的一致性变量,Y(k+t)为k+t时刻导弹群整体的视线角速度,Z1(k+t)为k+t时刻导弹群整体的视线角速度,Z2(k+t)为k+t时刻导弹群整体的视线视场角,U(k+t-1)为k+t-1时刻导弹群整体的控制变量;t=1,...,N;选择所述导弹状态依赖的线性模型的一致性变量、视线角速度和视场角作为输出变量,基于预测时域内状态变量输出变量和控制变量状态矩阵控制矩阵输出矩阵和获得预测控制模型为
其中,X(k)为k时刻导弹群整体状态变量。
步骤三:基于预测控制模型和各导弹打击时间协同的要求,获得性能指标函数。
根据各导弹打击时间协同的要求,即将各导弹打击时间协同目标转化为各导弹弹目距离和弹目相对速度一致性目标、各导弹控制输入的优化和驱动各导弹命中目标,性能指标函数表达式为
J=J1+J2+J3 (3)
其中
J1为协同变量代价函数,用于保证各导弹的弹目距离趋于一致且弹目相对速度趋于一致,
J2为控制变量代价函数,用于优化控制输入,
J3为视线角速度代价函数,用于驱动各导弹视线角速度归零,保证命中目标。其中,Im为m阶单位矩阵,u(i)(k)为k时刻第i枚导弹的控制变量,y(i)(k+1)为k+1时刻第i枚导弹的协同变量,为k+1时刻第i枚导弹的视线角速度,Nv为导弹总枚数,N为预测时域,R(i)、Q(i)为各导弹控制变量的权重和视线角速度的权重,为表示导弹通讯拓扑的邻接矩阵Aadj中的元素,
定义扩展拉普拉斯矩阵为其中L为表示导弹通讯拓扑的拉普拉斯矩阵,L=D-Aadj,D为入度矩阵,
定义导弹群整体的权重,和分别为协控制变量和视线角速度的整体权重
基于扩展拉普拉斯矩阵和导弹群整体的权重,则预测时域上的拉普拉斯矩阵协控制变量权重和视线角速度权重分别为
基于预测控制模型,将式(4)带入式(3),获得性能指标函数
其中为加权范数运算。
步骤四:基于性能指标函数和约束条件,依据二次规划形式转化为二次规划标准型。
基于性能指标函数,作如下变换:
其中H、F、constant分别为二次项系数、一次项系数和常数项,有
和均为预测时域内输出变量、为预测时域内控制变量,状态矩阵控制矩阵和分别为预测时域内输出矩阵;为加权范数运算,
针对各导弹法向加速度约束条件依据二次规划形式,则二次规划标准型为
其中
其中为N×Nv阶单位矩阵;1N为N×1列向量;为第i枚导弹的最大法向加速度;u(i)为第i枚导弹的控制变量;X(k)为k时刻导弹群整体的一致性变量;Umax为预测时域内控制变量;Λu为和组成的矩阵。
针对各导弹视场角约束条件依据二次规划形式,则二次规划标准型为
其中
其中,σ(i)为第i枚导弹的速度前置角(即视场角),为第i枚导弹的最大速度前置角(即最大视场角),i=1,...,Nv;为由和组成的矩阵。
步骤五:基于当前时刻各导弹状态、二次规划标准型的各导弹法向加速度和视场角的约束条件和性能指标函数,获得多弹协同制导的优化模型,利用凸优化方法(二次规划方法),获得当前时刻各导弹的最优控制序列,即为***输入量。
基于二次规划标准型的当前时刻各导弹状态、各导弹法向加速度和视场角的约束条件、性能指标函数,获得多弹协同制导的优化模型:
步骤六、从最优控制序列中选择各导弹的法向加速度并代入导弹状态依赖的线性模型,计算下一时刻的各导弹状态,按照步骤五的方式重复直至各导弹 命中目标。
本发明采用在线二次规划方法寻求最优解,但寻优方法不仅限于此。得到最优控制序列后,将其首向U(k)作为各导弹的输入量,输入制导***,完成闭环,并采样得到下一时刻的各导弹的状态:视线角λ(i)(k+1)、速度前置角σ(i)(k+1)、航向角θ(i)(k+1)、弹目距离r(i)(k+1),计算得到x(i)(k+1),作为初值,再进行下一时刻的优化问题求解,给出控制指令,重复该过程直到各导弹命中目标。
实施例:
实例场景为4枚导弹组成的导弹群攻击一静止目标,目标位置和各导弹位置、初始速度前置角如表1所示
表1导弹、目标初始运动参数
位置(千米) | 初始速度前置角(°) | |
导弹1 | (4.5528,-3.7266) | 14.3239 |
导弹2 | (2.9435,-1.4304) | 11.4592 |
导弹3 | (2.9793,3.8354) | -22.9183 |
导弹4 | (5.2758,5.9533) | -5.7296 |
目标 | (10,0) |
设置各导弹法向加速度限制为50米/秒2,视场角限制为0.75弧度,并假设各导弹速度为300米/秒。选择模型预测控制预测时域N=20,选择权重Q(i)=1R(i)=0.23、kr=0.67、导弹间的通讯拓扑为连通无向图。采用本发明实现4枚导弹同时打击一静止目标的弹道和各运动状态变化曲线如图3-8所示。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.基于输出一致性的多弹协同制导方法,其特征在于,具体步骤如下:
步骤一、给定各导弹法向加速度和视场角的约束条件;
步骤二、基于导弹和目标的相对运动关系,建立导弹状态依赖的线性模型,采用数值方法对所述导弹状态依赖的线性模型进行离散化,基于离散化后的导弹状态依赖的线性模型,根据设定的预测时域建立预测控制模型;
步骤三、基于预测控制模型和各导弹打击时间协同的要求,获得性能指标函数;
步骤四、基于二次规划形式,将性能指标函数和约束条件转化为二次规划标准型;
步骤五、基于当前时刻各导弹状态、二次规划标准型的各导弹法向加速度和视场角的约束条件和性能指标函数,获得多弹协同制导的优化模型,基于多弹协同制导的优化模型,利用凸优化方法,获得当前时刻各导弹的最优控制序列;
步骤六、从最优控制序列中选择各导弹的法向加速度并代入导弹状态依赖的线性模型,计算下一时刻的各导弹状态,按照步骤五的方式重复直至各导弹命中目标。
2.如权利要求1所述基于输出一致性的多弹协同制导方法,其特征在于,步骤一具体实现过程如下:
定义编队中第i枚导弹的法向加速度约束条件为其中a(i)为第i枚导弹的法向加速度,为第i枚导弹的最大法向加速度;第i枚导弹的视场角的约束条件为σ(i)为第i枚导弹的速度前置角,为第i枚导弹的速度最大前置角。
3.如权利要求1所述基于输出一致性的多弹协同制导方法,其特征在于, 步骤二具体实现过程如下:
选取第i枚导弹的状态变量 为第i枚导弹状态变量的变化率,第i枚导弹的控制变量u(i)=a(i),其中i,j=1,...,Nv,i≠j,根据导弹和目标的相对运动关系,建立导弹状态依赖的线性模型:
其中a(i)表示第i枚导弹的法向加速度;A(i)为第i枚导弹的状态矩阵、B(i)为第i枚导弹的输入矩阵、y(i)、和分别为第i枚导弹的一致性变量、视线角速度和视场角,分别为第i枚导弹的一致性变量、视线角速度和视场角的输出矩阵,分别为r(i),λ(i),σ(i)对应的赋予变量,其表示形式分别如下:
r(i)表示第i枚导弹的弹目距离,表示第i枚导弹的弹目相对速度,λ(i)表示第i枚导弹的视线角,表示第i枚导弹的视线角速度,σ(i)为第i枚导弹的速度前置角,V表示导弹飞行速度,kr和分别为弹目距离的权重和弹目相对速度的权重;
采用数值方法对状态矩阵A(i)、输入矩阵B(i)进行离散化,获得第i枚导弹的导弹状态依赖的线性模型的离散化形式
定义k时刻导弹群整体状态变量X(k),k时刻导弹群整体的一致性变量、视线角速度和视场角分别为Y(k)、Z1(k)、Z2(k),控制变量U(k),X(k)、Y(k)、Z1(k)、Z2(k)和U(k)分别为
其中,x(i)(k)为k时刻第i枚导弹状态变量;y(i)(k)为k时刻第i枚导弹的一致性变量;为k时刻第i枚导弹的视线角速度;为k时刻第i枚导弹的视场角;u(i)(k)为k时刻第i枚导弹的控制变量,i=1,...,Nv;
定义当前时刻导弹群整体状态矩阵当前时刻导弹群整体状态的一致性变量、视线角速度和视场角的输出矩阵分别为控制矩阵 和分别为
其中blkdiag为分块对角矩阵运算;
定义预测控制模型的预测时域为N,获得预测时域内状态变量预测时域内导弹群整体状态的一致性变量、视线角速度和视场角的输出矩阵的输出变量分别为和控制变量状态矩阵控制矩阵预测时域内导弹群整体状态的一致性变量、视线角速度和视场角的输出矩阵的输出矩阵分别为 和和分别为
其中IM为M阶单位矩阵,M=N,为矩阵的克罗内克积运算;X(k+t)为k+t时刻导弹群整体的一致性变量,Y(k+t)为k+t时刻导弹群整体的视线角速度,Z1(k+t)为k+t时刻导弹群整体的视线角速度,Z2(k+t)为k+t时刻导弹群整体的视线视场角,U(k+t-1)为k+t-1时刻导弹群整体的控制变量;t=1,...,N;
选择所述导弹状态依赖的线性模型的一致性变量、视线角速度和视场角作为 输出变量,基于预测时域内状态变量输出变量和控制变量状态矩阵控制矩阵输出矩阵和获得预测控制模型为
其中,X(k)为k时刻导弹群整体状态变量。
4.如权利要求1或3所述基于输出一致性的多弹协同制导方法,其特征在于,所述数值方法为欧拉方法、改进欧拉方法、龙格-库塔法或线性多步法。
5.如权利要求1所述基于输出一致性的多弹协同制导方法,其特征在于,所述凸优化方法为二次规划方法。
6.如权利要求1所述基于输出一致性的多弹协同制导方法,其特征在于,步骤三具体实现过程如下:根据各导弹打击时间协同的要求,即将各导弹打击时间协同目标转化为各导弹弹目距离和弹目相对速度一致性目标、各导弹控制输入的优化和驱动各导弹命中目标,性能指标函数表达式为
J=J1+J2+J3 (3)
其中,J1为协同变量代价函数,J2为控制变量代价函数,J3为视线角速度代价函数,其分别为:
Im为m阶单位矩阵,u(i)(k)为k时刻第i枚导弹的控制变量,y(i)(k+1)为k+1时刻第i枚导弹的协同变量,为k+1时刻第i枚导弹的视线角速度,Nv为导 弹总枚数,N为预测时域,R(i)、Q(i)为各导弹控制变量的权重和视线角速度的权重,为表示导弹通讯拓扑的邻接矩阵Aadj中的元素,
定义扩展拉普拉斯矩阵为其中L为表示导弹通讯拓扑的拉普拉斯矩阵,L=D-Aadj,D为入度矩阵,Im为m阶单位矩阵。
定义导弹群整体的权重,和分别为协控制变量和视线角速度的整体权重
基于扩展拉普拉斯矩阵和导弹群整体的权重,则预测时域上的拉普拉斯矩阵协控制变量权重和视线角速度权重分别为
其中,IN为N阶单位矩阵,为矩阵的克罗内克积运算;
基于预测控制模型,将式(4)带入式(3),获得性能指标函数J
其中为加权范数运算,
7.如权利要求1所述基于输出一致性的多弹协同制导方法,其特征在于,步骤四具体实现过程如下:
基于性能指标函数,作如下变换:
其中H、F、constant分别为二次项系数、一次项系数和常数项,有
和均为预测时域内输出变量、为预测时域内控制变量,状态矩阵控制矩阵和分别为预测时域内输出矩阵;为加权范数运算,
针对各导弹法向加速度约束条件依据二次规划形式,则二次规划标准型为
其中
其中为N×Nv阶单位矩阵;1N为N×1列向量;为第i枚导弹的最大法向加速度;u(i)为第i枚导弹的控制变量;X(k)为k时刻导弹群整体的一致性变量;Umax为预测时域内控制变量;
针对各导弹视场角约束条件依据二次规划形式,则二次规划标准型为
其中
其中,σ(i)为第i枚导弹的速度前置角,为第i枚导弹的最大速度前置角,i=1,...,Nv。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610402960.XA CN106352738B (zh) | 2016-06-08 | 2016-06-08 | 基于输出一致性的多弹协同制导方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610402960.XA CN106352738B (zh) | 2016-06-08 | 2016-06-08 | 基于输出一致性的多弹协同制导方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106352738A true CN106352738A (zh) | 2017-01-25 |
CN106352738B CN106352738B (zh) | 2018-04-03 |
Family
ID=57843113
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610402960.XA Expired - Fee Related CN106352738B (zh) | 2016-06-08 | 2016-06-08 | 基于输出一致性的多弹协同制导方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106352738B (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107179021A (zh) * | 2017-06-14 | 2017-09-19 | 北京理工大学 | 一种驾束制导体制下多弹协同高精度制导控制方法 |
CN107818219A (zh) * | 2017-10-31 | 2018-03-20 | 中国人民解放军国防科技大学 | 一种面向突防的多导弹协同弹道规划方法 |
CN108008736A (zh) * | 2017-11-21 | 2018-05-08 | 山东科技大学 | 飞行器协同控制方法、装置、计算机可读存储介质及终端 |
CN108362171A (zh) * | 2018-02-02 | 2018-08-03 | 北京航空航天大学 | 一种具有攻击时间和攻击角度约束的制导策略 |
CN110162084A (zh) * | 2018-10-11 | 2019-08-23 | 北京机电工程研究所 | 基于一致性理论的飞航导弹集群***编队控制方法 |
CN110296634A (zh) * | 2019-07-04 | 2019-10-01 | 哈尔滨工业大学 | 分体式制导的探测模式及精度分析方法 |
CN110375596A (zh) * | 2019-07-19 | 2019-10-25 | 北京机电工程研究所 | 一种面向协同制导仿真***的时间一致性控制方法 |
CN110716585A (zh) * | 2019-10-31 | 2020-01-21 | 西安电子科技大学 | 基于分布式一致性与自组网的自主协同控制方法及*** |
CN111121539A (zh) * | 2019-12-20 | 2020-05-08 | 上海航天控制技术研究所 | 一种基于高速数据链的双目协同探测方法及装置 |
CN111506101A (zh) * | 2019-10-21 | 2020-08-07 | 北京理工大学 | 基于通信网络拓扑结构的飞行器协同制导控制方法及*** |
CN111707148A (zh) * | 2020-06-28 | 2020-09-25 | 北京航空航天大学 | 一种时变时延条件下多导弹协同制导方法及*** |
CN112526872A (zh) * | 2020-12-04 | 2021-03-19 | 北京理工大学 | 带大落角约束中制导和末制导交接与制导信息的处理方法 |
CN112558631A (zh) * | 2020-12-04 | 2021-03-26 | 北京理工大学 | 基于测量信息的带大落角约束变参数制导方法 |
CN112902767A (zh) * | 2021-01-28 | 2021-06-04 | 西安交通大学 | 一种多弹时间协同的导弹制导方法及*** |
CN114216368A (zh) * | 2021-12-07 | 2022-03-22 | 南京理工大学 | 一种变速条件下弹群协同制导方法及*** |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102890506A (zh) * | 2011-07-19 | 2013-01-23 | 北京理工大学 | 一种基于约束规划的小天体接近段制导控制方法 |
CN102945000A (zh) * | 2012-11-29 | 2013-02-27 | 北京理工大学 | 基于可观测性约束的行星着陆轨迹随机优化方法 |
CN103245257A (zh) * | 2013-04-11 | 2013-08-14 | 北京航空航天大学 | 基于Bezier曲线的多约束飞行器导引方法 |
CN103245256A (zh) * | 2013-04-25 | 2013-08-14 | 北京理工大学 | 一种多导弹协同作战制导律设计方法 |
CN103412485A (zh) * | 2013-07-22 | 2013-11-27 | 西北工业大学 | 基于滚动优化策略的刚体航天器姿态机动路径规划方法 |
CN103499974A (zh) * | 2013-09-27 | 2014-01-08 | 中国人民解放军空军工程大学 | 一种双机协同被动雷达探测路径规划方法 |
CN103884237A (zh) * | 2014-04-08 | 2014-06-25 | 哈尔滨工业大学 | 基于目标概率分布信息的多对一协同制导方法 |
CN103995539A (zh) * | 2014-05-15 | 2014-08-20 | 北京航空航天大学 | 一种无人机自主编队评价指标与mpc编队控制方法 |
CN105202972A (zh) * | 2015-09-16 | 2015-12-30 | 北京理工大学 | 一种基于模型预测控制技术的多导弹协同作战制导方法 |
-
2016
- 2016-06-08 CN CN201610402960.XA patent/CN106352738B/zh not_active Expired - Fee Related
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102890506A (zh) * | 2011-07-19 | 2013-01-23 | 北京理工大学 | 一种基于约束规划的小天体接近段制导控制方法 |
CN102945000A (zh) * | 2012-11-29 | 2013-02-27 | 北京理工大学 | 基于可观测性约束的行星着陆轨迹随机优化方法 |
CN103245257A (zh) * | 2013-04-11 | 2013-08-14 | 北京航空航天大学 | 基于Bezier曲线的多约束飞行器导引方法 |
CN103245256A (zh) * | 2013-04-25 | 2013-08-14 | 北京理工大学 | 一种多导弹协同作战制导律设计方法 |
CN103412485A (zh) * | 2013-07-22 | 2013-11-27 | 西北工业大学 | 基于滚动优化策略的刚体航天器姿态机动路径规划方法 |
CN103499974A (zh) * | 2013-09-27 | 2014-01-08 | 中国人民解放军空军工程大学 | 一种双机协同被动雷达探测路径规划方法 |
CN103884237A (zh) * | 2014-04-08 | 2014-06-25 | 哈尔滨工业大学 | 基于目标概率分布信息的多对一协同制导方法 |
CN103995539A (zh) * | 2014-05-15 | 2014-08-20 | 北京航空航天大学 | 一种无人机自主编队评价指标与mpc编队控制方法 |
CN105202972A (zh) * | 2015-09-16 | 2015-12-30 | 北京理工大学 | 一种基于模型预测控制技术的多导弹协同作战制导方法 |
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107179021A (zh) * | 2017-06-14 | 2017-09-19 | 北京理工大学 | 一种驾束制导体制下多弹协同高精度制导控制方法 |
CN107179021B (zh) * | 2017-06-14 | 2018-08-28 | 北京理工大学 | 一种驾束制导体制下多弹协同高精度制导控制方法 |
CN107818219A (zh) * | 2017-10-31 | 2018-03-20 | 中国人民解放军国防科技大学 | 一种面向突防的多导弹协同弹道规划方法 |
CN107818219B (zh) * | 2017-10-31 | 2018-09-25 | 中国人民解放军国防科技大学 | 一种面向突防的多导弹协同弹道规划方法 |
CN108008736A (zh) * | 2017-11-21 | 2018-05-08 | 山东科技大学 | 飞行器协同控制方法、装置、计算机可读存储介质及终端 |
CN108008736B (zh) * | 2017-11-21 | 2020-12-01 | 山东科技大学 | 飞行器协同控制方法、装置、计算机可读存储介质及终端 |
CN108362171A (zh) * | 2018-02-02 | 2018-08-03 | 北京航空航天大学 | 一种具有攻击时间和攻击角度约束的制导策略 |
CN108362171B (zh) * | 2018-02-02 | 2019-03-29 | 北京航空航天大学 | 一种具有攻击时间和攻击角度约束的制导策略 |
CN110162084A (zh) * | 2018-10-11 | 2019-08-23 | 北京机电工程研究所 | 基于一致性理论的飞航导弹集群***编队控制方法 |
CN110296634A (zh) * | 2019-07-04 | 2019-10-01 | 哈尔滨工业大学 | 分体式制导的探测模式及精度分析方法 |
CN110296634B (zh) * | 2019-07-04 | 2021-07-13 | 哈尔滨工业大学 | 分体式制导的探测模式及精度分析方法 |
CN110375596A (zh) * | 2019-07-19 | 2019-10-25 | 北京机电工程研究所 | 一种面向协同制导仿真***的时间一致性控制方法 |
CN111506101A (zh) * | 2019-10-21 | 2020-08-07 | 北京理工大学 | 基于通信网络拓扑结构的飞行器协同制导控制方法及*** |
CN111506101B (zh) * | 2019-10-21 | 2021-04-20 | 北京理工大学 | 基于通信网络拓扑结构的飞行器协同制导控制方法及*** |
CN110716585B (zh) * | 2019-10-31 | 2020-08-21 | 西安电子科技大学 | 基于分布式一致性与自组网的自主协同控制方法及*** |
CN110716585A (zh) * | 2019-10-31 | 2020-01-21 | 西安电子科技大学 | 基于分布式一致性与自组网的自主协同控制方法及*** |
CN111121539A (zh) * | 2019-12-20 | 2020-05-08 | 上海航天控制技术研究所 | 一种基于高速数据链的双目协同探测方法及装置 |
CN111121539B (zh) * | 2019-12-20 | 2022-04-22 | 上海航天控制技术研究所 | 一种基于高速数据链的双目协同探测方法及装置 |
CN111707148A (zh) * | 2020-06-28 | 2020-09-25 | 北京航空航天大学 | 一种时变时延条件下多导弹协同制导方法及*** |
CN112526872A (zh) * | 2020-12-04 | 2021-03-19 | 北京理工大学 | 带大落角约束中制导和末制导交接与制导信息的处理方法 |
CN112558631A (zh) * | 2020-12-04 | 2021-03-26 | 北京理工大学 | 基于测量信息的带大落角约束变参数制导方法 |
CN112526872B (zh) * | 2020-12-04 | 2021-09-24 | 北京理工大学 | 带大落角约束中制导和末制导交接与制导信息的处理方法 |
CN112902767A (zh) * | 2021-01-28 | 2021-06-04 | 西安交通大学 | 一种多弹时间协同的导弹制导方法及*** |
CN112902767B (zh) * | 2021-01-28 | 2022-05-20 | 西安交通大学 | 一种多弹时间协同的导弹制导方法及*** |
CN114216368A (zh) * | 2021-12-07 | 2022-03-22 | 南京理工大学 | 一种变速条件下弹群协同制导方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN106352738B (zh) | 2018-04-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106352738A (zh) | 基于输出一致性的多弹协同制导方法 | |
CN108168381B (zh) | 一种多枚导弹协同作战的控制方法 | |
CN105759612B (zh) | 带落角约束的微分对策反拦截机动突防/精确打击导引方法 | |
CN103090728B (zh) | 一种基于滑模控制的带末角约束制导方法 | |
CN112902767B (zh) | 一种多弹时间协同的导弹制导方法及*** | |
CN110928329A (zh) | 一种基于深度q学习算法的多飞行器航迹规划方法 | |
CN109084641B (zh) | 导弹制导方法及装置 | |
CN109597423A (zh) | 一种基于参考视线角信号的多约束末制导律的设计方法 | |
CN104215935B (zh) | 一种基于决策加权融合的雷达炮弹目标识别方法 | |
CN114063644B (zh) | 基于鸽群反向对抗学习的无人作战飞机空战自主决策方法 | |
CN114003050B (zh) | 一种基于微分博弈的三体对抗策略的主动防御制导方法 | |
CN113282061A (zh) | 一种基于课程学习的无人机空中博弈对抗的解决方法 | |
CN113741525A (zh) | 基于策略集合maddpg多无人机协同攻防对抗方法 | |
CN110187640A (zh) | 针对机动目标和允许通信时滞的多导弹协同作战制导律设计方法 | |
CN116858039A (zh) | 一种高超声速飞行器博弈制导方法、***、设备及介质 | |
Kung et al. | Missile guidance algorithm design using particle swarm optimization | |
CN106091816A (zh) | 一种基于滑模变结构理论的半捷联空空导弹制导方法 | |
CN112379600A (zh) | 一种基于通信时变延迟的分布式协同制导律构建方法 | |
CN112099348A (zh) | 一种基于观测器和全局滑模的碰撞角控制制导方法 | |
CN113359819B (zh) | 一种带有碰撞角约束和加速度限制的最优制导律 | |
Shi et al. | Research on self-adaptive decision-making mechanism for competition strategies in robot soccer | |
CN109543135A (zh) | 一种基于信息最大化的多auv协同同时导引方法 | |
Yang et al. | Ballistic missile maneuver penetration based on reinforcement learning | |
CN113110576A (zh) | 一种可连续切换的自适应固定时间收敛协同制导方法 | |
Saleem et al. | Three dimensional nonlinear guidance law for exact impact time control |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180403 Termination date: 20190608 |
|
CF01 | Termination of patent right due to non-payment of annual fee |