CN115629618A - 一种基于落点选择和伪谱法的分离体最优弹道规划的方法 - Google Patents

一种基于落点选择和伪谱法的分离体最优弹道规划的方法 Download PDF

Info

Publication number
CN115629618A
CN115629618A CN202211379122.7A CN202211379122A CN115629618A CN 115629618 A CN115629618 A CN 115629618A CN 202211379122 A CN202211379122 A CN 202211379122A CN 115629618 A CN115629618 A CN 115629618A
Authority
CN
China
Prior art keywords
trajectory
feasible
planning
separation body
ballistic
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
Application number
CN202211379122.7A
Other languages
English (en)
Inventor
刘李雷
黄晓平
唐梦莹
岳小飞
王志军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
CASIC Rocket Technology Co
Original Assignee
CASIC Rocket Technology Co
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by CASIC Rocket Technology Co filed Critical CASIC Rocket Technology Co
Priority to CN202211379122.7A priority Critical patent/CN115629618A/zh
Publication of CN115629618A publication Critical patent/CN115629618A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/10Simultaneous control of position or course in three dimensions
    • G05D1/107Simultaneous control of position or course in three dimensions specially adapted for missiles

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本发明涉及一种基于落点选择和伪谱法的分离体最优弹道规划的方法,包括:选取经验可行的分离体落点信息;设计分离体弹道规划约束条件;分离体弹道规划,根据分离体弹道规划约束条件进行求解优化;将分离体落点信息输入分离体弹道进行求解优化,若得到一个可行弹道作为分离体飞行的标准弹道;若不能得到可行弹道,则重新获取可行的分离体落点信息直至得到可行弹道。本发明提出一种基于落点选择和伪谱法的分离体最优弹道规划的方法,在射前离线确定分离体的弹道,使得分离体的理论落点位于选定的安全区域中。

Description

一种基于落点选择和伪谱法的分离体最优弹道规划的方法
技术领域
本发明属于运载火箭弹道设计领域的延伸和拓展具体是分离体落区弹道规划算法技术领域,具体涉及一种基于落点选择和伪谱法的分离体最优弹道规划的方法。
背景技术
随着航天技术的发展,特别是卫星设计与生产的集成化和小型化、发射需求的密集化,市场对于低成本、运载能力500kg以下的轻型运载火箭的需求有了大幅提升。为了追求低成本与灵活性,具备快速响应能力的轻型运载火箭逐步开始采用简易发射装置或发射车作为发射平台,更有甚者,考虑到发射任务的特殊性或者为了确保发射任务的隐蔽性,发射点可能有全地域和全地形的需求。但是发射点和发射平台的灵活多变给落区的选择带来的较大的挑战,比如某一发射点可能无法满足某类任务的落点要求,因此本质上限制了发射任务的灵活性。为了解决灵活发射情况下的落区控制问题,在不影响火箭最优弹道的前提下,最为可行的方案是分离体本身应当具备主动落区选择和控制能力。
在任务规划之初,首先应确保规划出的火箭标准弹道应完全发挥火箭的可用运载能力,在此基础上通过落区规划算法规划出各级分离体满足落区要求的分离体标准弹道,最后在飞行试验中分离体按飞行程序分离后,按照分离体标准弹道,通过在线计算姿控和制导指令受控在指定的落点附近着陆。为此,需要在射前确定分离体的弹道,使得分离体的理论落点位于安全区域中,保证落点附近人员和建筑的安全。
发明内容
针对现有技术的不足,本发明提出一种基于落点选择和伪谱法的分离体最优弹道规划的方法,在射前离线确定分离体的弹道,使得分离体的理论落点位于选定的安全区域中。
为实现上述目的,本发明提供的技术方案如下:
一种基于落点选择和伪谱法的分离体最优弹道规划的方法,包括如下步骤:
选取经验可行的分离体落点信息;
设计分离体弹道规划约束条件;
分离体弹道规划,根据分离体弹道规划约束条件采用Legendre伪谱法对分离体弹道进行求解优化;
将所述分离体落点信息输入分离体弹道进行求解优化,若得到一个可行弹道则将所述可行弹道输出,作为分离体飞行的标准弹道;若不能得到可行弹道,则重新获取可行的分离体落点信息直至得到可行弹道。
进一步地,所述选取经验可行的分离体落点信息包括
根据发射点、分离体分离点以及人物活动区域的位置关系,要求分离体落点散布在人、物频繁活动区域周边,选取经验可行的分离体落点n个,所述n个分离体落点对应n条分离体弹道,所述n为自然数1,2......i。
进一步地,所述每一条分离体弹道规划中设计分离体约束条件包括
获取分离体总体参数,包括质量和质心参数;
气动和载荷约束,根据分离体通过气动力做弱机动进而改变落点的特点,在规划弹道时利用满足要求的攻角和侧滑角作为控制变量得到在分离体上产生的横法向过载值,所述横、法向过载要求其随时间的积分最小,为此选取如下二次型目标泛函:
Figure BDA0003927570890000031
由于分离体在飞行过程中所受到的横法向过载均由气动力产生,且与攻角、侧滑角的大小直接相关,因此式可等效为如下形式:
Figure BDA0003927570890000032
式中:
k1,k2,kα,kβ——人为确定的系数;
ny——横向过载;
nz——法向过载;
α——攻角(不考虑风);
β——侧滑角(不考虑风);
所述满足要求的攻角和侧滑角跟俯仰舵偏角和偏航舵偏角的关系即分离体所受约束条件如下:
Figure BDA0003927570890000033
式中
q=0.5ρV2——弹道动压
s——参考面积
l——参考长度
Cm——俯仰力矩系数
Cn——偏航力矩系数
Ma——马赫数
δp——俯仰舵偏角
δr——偏航舵偏角
δmax——最大允许舵偏角
k——舵效限制(一般小于1,需要留取一定舵偏余量用于姿态调整);
分离体落点要求满足约束如下:
Figure BDA0003927570890000041
式中
tf——着陆时间(规划前不确定)
(x,y,z)——落点在发射系中的坐标
i——落点编号,Xi,Yi,Zi是可行落点的坐标,在发射系中的坐标
动力学方程构建。一般分离体质心运动方程可表示如下:
Figure BDA0003927570890000051
式中
V——分离体速度,初值为分离时刻标准弹道速度,向量;
X——分离体在发射系中的坐标;
m——分离体质量。
进一步地,所述根据分离体约束条件对分离体弹道进行求解优化包括
为简单起见,记控制变量(向量)为
Figure BDA0003927570890000054
记状态变量(向量)为
Figure BDA0003927570890000055
将静态约束公式(1.3)的内容记为向量g(s,u)≤0,0作为向量,根据公式(1.3)能表达成小于等于0的标量形式;将终端约束公式(1.4)的内容记为向量hi(s)=0,0作为向量;对应第i个落点的目标泛函公式(1.2)记为Ji(u),将动力学与运动学方程公式(1.5)记为
Figure BDA0003927570890000052
那么用于规划对应第i个落点的优化问题可按下列提法给出:
Figure BDA0003927570890000053
记分离时刻为t0=0,着陆时间为tf,做如下微分变换:
Figure BDA0003927570890000061
利用Legendre伪谱法选择离散点τ01,…,τN,其中τ0=-1,τN=1,那么状态量可由插值点近似表示为
Figure BDA0003927570890000062
上式中s对τ求导可得
Figure BDA0003927570890000063
并将目标泛函转化为离散形式后,可以给出弹道规划问题的提法如下:
Figure BDA0003927570890000071
由上,弹道规划优化问题已经转变为关于控制变量序列u(τ0),u(τ1),…,u(τN)的最优化问题,通过一般的优化算法方便求解,终止条件为优化算法收敛或者运行时间超过一定值。
进一步地,所述若得到一个可行弹道则将所述可行弹道输出,作为分离体飞行的标准弹道,包括
所述n个分离体落点对应n条分离体弹道分别采用Legendre伪谱法对分离体弹道进行求解优化,若得到任何一条可行弹道则终止计算,将所述可行弹道输出,作为分离体飞行的标准弹道。
进一步地,若不能得到可行弹道,则重新获取可行的分离体落点信息直至得到可行弹道,包括
所述n个分离落点在现有约束条件下均不能得到可行弹道(不可达(算法不收敛,超时终止)),则为规划失败落点;根据优化后的控制变量序列u(τ0),u(τ1),…,u(τN)的取值,输出期望的分离体落点相对规划失败落点的可行方向,根据输出的可行方向重新获取可行的分离体落点信息,重新进行弹道规划,直至得到满足要求的分离体弹道为止。
进一步地,所述根据输出的可行方向重新选取可行的分离体落点信息,包括
在优化后的控制变量序列u(τ0),u(τ1),,u(τN)中如果攻角连续K个取值均达到kδmax,则落点相对分离体初始位置沿射向过远,则选取沿速度方向较规划失败落点近的落点;同样地,若攻角连续K个均达到-kδmax,则落点相对分离体初始位置沿射向过近,则选取沿速度方向较规划失败落点远的落点;同样,如果侧滑角连续K个取值均达到kδmax或-kδmax,则落点相对分离体初始位置沿速度垂向过远,则选取沿速度垂向较规划失败落点近的落点。
另一方面,本发明还提供一种计算机可读的存储介质,所述计算机可读的存储介质包括存储的程序,其中,所述程序运行时执行前述的基于落点选择和伪谱法的分离体最优弹道规划的方法。
再者,本发明还提供一种电子装置,包括存储器和处理器,所述存储器中存储有计算机程序,所述处理器被设置为通过所述计算机程序执行前述的基于落点选择和伪谱法的分离体最优弹道规划的方法。
本发明所构思的以上技术方案与现有技术相比,能够取得下列有益效果:
本发明设计了一种基于落点筛选和伪谱法的分离体落点最优弹道生成算法和装置。主要内容按实施的先后顺序包括:首先根据人口分布确定可行落点,之后按照气动约束、载荷约束按照弹道不飞跃人口稠密区优先,其他弹道次之的顺序依次尝试利用伪谱法生成最优弹道,直到最优弹道有解收敛。若算法无法给出最优弹道,则给出落点的大致可行方向,为重新选择落点提供参考。本发明作为火箭无动力带气动舵分离体的落点约束最优弹道生成手段,算法的输入为标准弹道对应的分离体质心位置和速度,输出为分离体最优弹道或者落点调整远近的方向指示建议,根据输出的可行方向指示建议重新获取可行的分离体落点信息,重新进行弹道规划,直至得到满足要求的分离体弹道为止,做到了在射前确定分离体的弹道,使得分离体的理论落点位于安全区域中。
附图说明
图1为分离体可行落点与不可行区域的相对位置关系示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,均属于本发明保护的范围。
本发明的实施例提供一种基于落点选择和伪谱法的分离体最优弹道规划的方法,包括如下步骤:
步骤S1,选取经验可行的分离体落点信息;根据发射点、分离体分离点以及人物活动区域的位置关系,要求分离体落点散布在人、物频繁活动区域周边,选取经验可行的分离体落点n个,所述n个分离体落点对应n条分离体弹道,所述n为自然数1,2......i......n,,一般人为主观根据发射点位置,附近人物活动区域选取,理论上不一定可达,可行落点实际上不止一个,但应当都在射面附近。
步骤S2,设计分离体弹道规划约束条件,获取分离体总体参数,包括质量和质心参数;考虑气动和载荷约束,根据分离体通过气动力做弱机动进而改变落点的特点,在规划弹道时利用满足要求的攻角和侧滑角作为控制变量得到在分离体上产生的横法向过载值,考虑到分离体本质上为一薄壳体,其所能够承受的横法向过载有限,在设计规划算法的时候应当尽量降低理论最优弹道上分离体壳体所受到的横法向过载,因此所述横、法向过载要求其随时间的积分最小,为此选取如下二次型目标泛函:
Figure BDA0003927570890000101
由于分离体在飞行过程中所受到的横法向过载均由气动力产生,且与攻角、侧滑角的大小直接相关,因此式可等效为如下形式:
Figure BDA0003927570890000102
(1.1)式和(1.2)式中:
k1,k2,kα,kβ——人为确定的系数;首先取一个值,然后根据结果调整,k1大于k2,取值范围[0.5,5],根据火箭的具体运载能力的不同进行改变,但必须是正值,不等于0,一般从1开始选取;
kα与k1成正相关,kβ与k2成正相关,具体比值大小与气动参数有关,不同的火箭取值不同。实际优化时可直接选取kα和kβ,不一定要从k1和k2的取值开始,在实际操作时也是直接选取kα和kβ开始试探求解,在此写出公式(1.1)中的系数k1和k2,是为了说明如何将过载最小这一优化目标转化为公式(1.2),其实公式(1.1)并不参与最后的优化求解。
ny——横向过载;
nz——法向过载;
α——攻角(不考虑风)
β——侧滑角(不考虑风);
在构造约束条件时,不考虑环境热、结构强度的情况下,从控制专业考虑,分离体能否按照设计的弹道程序姿态稳定飞行,很大程度上取决于在规划出来的弹道上,分离体所受合力矩是否能够达到平衡;由于分离体是一旋转对称体,在选取约束条件时只需要考虑俯仰和偏航两个通道;若以质心作为参考点,所述满足要求的攻角和侧滑角跟俯仰舵偏角和偏航舵偏角的关系即分离体所受约束条件如下:
Figure BDA0003927570890000111
式中所有符号都是标量
q=0.5ρV2——弹道动压,ρ为当地大气密度,V为来流速度(没有风时是箭体速度),是标量,其他公式例如(1.5)中的速度V是箭体速度,是向量;
s——参考面积
l——参考长度
Cm——俯仰力矩系数,与(α,β,Ma,δp,δr)相关的俯仰力矩系数
Cn——偏航力矩系数
Ma——马赫数
δp——俯仰舵偏角
δr——偏航舵偏角
δmax——最大允许舵偏角
k——舵效限制(一般小于1,需要留取一定舵偏余量用于姿态调整);
分离体落点要求满足约束如下:
Figure BDA0003927570890000121
式中
tf——着陆时间(规划前不确定)
(x,y,z)——落点在发射系中的坐标
i——落点编号,Xi,Yi,Zi是可行落点的坐标,在发射系中的坐标动力学方程构建。一般分离体质心运动方程可表示如下:
Figure BDA0003927570890000131
式中
V——分离体速度,初值为分离时刻标准弹道速度,向量;头上“.”代表求导;
X——分离体在发射系中的坐标,三维位置向量,代表发射系中xyz三个标量;头上“.”代表求导;
m——分离体质量;
F——是作用在分离体上的合力。
步骤S3,分离体弹道规划,根据分离体弹道规划约束条件采用Legendre伪谱法对分离体弹道进行求解优化;
优化问题提法如下:
为简单起见,记控制变量(向量)为
Figure BDA0003927570890000133
记状态变量(向量)为
Figure BDA0003927570890000134
将静态约束公式(1.3)的内容记为向量g(s,u)≤0,0作为向量,根据公式(1.3)能表达成小于等于0的标量形式;将终端约束公式(1.4)的内容记为向量hi(s)=0,0作为向量;对应第i个落点的目标泛函公式(1.2)记为Ji(u),将动力学与运动学方程公式(1.5)记为
Figure BDA0003927570890000132
那么用于规划对应第i个落点的优化问题可按下列提法给出:
Figure BDA0003927570890000141
上式中的minuJi(u)式物理意义为找到一个控制变量u使得目标泛函最小;s.t.代表约束条件(简写);
公式(1.6)的解法为:记分离时刻为t0=0,着陆时间为tf,做如下微分变换:
Figure BDA0003927570890000142
Figure BDA0003927570890000143
的物理意义为如果t=0,τ就等于-1,如果t=tf,τ就等于1;x(1)、y(1)和z(1)代表落点坐标,τ等于1时候的落点坐标值。
利用Legendre伪谱法选择离散点τ01,…,τN,其中τ0=-1,τN=1,τ0=-1,是分离时刻,τN=1是落点时刻,那么状态量可由插值点近似表示为
Figure BDA0003927570890000151
s(τ)是状态向量,u(τ)是控制变量,在对状态向量和控制变量进行插值时,需要选择离散点,依据离散点拟合曲线,在拟合曲线时,对不同离散点的权重是不一样的,ωi(τ)是权重;
上式中s对τ求导可得
Figure BDA0003927570890000152
Di(τ)即是ωi(τ)求导再做累加;
并将目标泛函转化为离散形式后,可以给出弹道规划问题的提法如下:
Figure BDA0003927570890000153
由上,弹道规划问题已经转变为关于控制变量序列u(τ0),u(τ1),…,u(τN)的最优化问题,可通过一般的优化算法方便求解,终止条件为优化算法收敛或者运行时间超过一定值,具体值需结合具体平台考虑,一般性能的计算机,能正常收敛的情况下,半小时左右属正常。发散的情况则要强制停止解算。即是把公式(1.6)转换成公式(1.10)的离散形式,便于计算机解算。实际是通过优化该问题最后输出沿弹道的攻角和侧滑角,进而计算出对应的俯仰角和偏航角。
步骤S4,将所述分离体落点信息输入分离体弹道进行求解优化,若得到一个可行弹道则将所述可行弹道输出,作为分离体飞行的标准弹道;若不能得到可行弹道,则重新获取可行的分离体落点信息直至得到可行弹道。
具体地,若得到一个可行弹道则将所述可行弹道输出,作为分离体飞行的标准弹道,本实施例中四个分离体落点对应四条分离体弹道分别采用Legendre伪谱法对分离体弹道进行求解优化,若得到任何一条可行弹道则终止计算,将所述可行弹道输出,作为分离体飞行的标准弹道。
假设可行落点与不可行区域的相对位置关系如图1所示,图中射点为O,分离体与箭体(继续飞行的部分)分离点为F,散布在人、物频繁活动区域G周边的四个可行落区分别以A、B、C和D为中心。箭体首先在射面内沿弹道OF飞行至分离点F,之后分离体从F点开始以可行弹道飞行至A、B、C和D四个落点(假设四个落点均可达),这四条弹道都可以利用本发明的方法步骤S1-S3给出,若无法给出可行弹道,即是运行时间超过一定值,强制停止解算的情况,则对应落点为不可达落点,从落点方案中剔除。在实际生成弹道过程中,考虑到计算时间限制,一般只要得到一个可行弹道之后便可终止计算,并将该可行弹道输出,作为分离体飞行的标准弹道。
若不能得到可行弹道,则重新获取可行的分离体落点信息直至得到可行弹道,本实施例中四个分离落点在现有约束条件下均不能得到可行弹道(即不可达,算法不收敛,超时终止),则为规划失败落点;根据优化后的控制变量序列u(τ0),u(τ1),…,u(τN)的取值(即是没有收敛或强制终止之后得到的控制变量序列),输出期望的分离体落点相对规划失败落点的可行方向,根据输出的可行方向重新获取可行的分离体落点信息,具体给出可行弹道方向提示,再通过调整期望落点信息输入后,重新进行分离体弹道规划,直至得到满足要求的分离体弹道为止。虽然算法不收敛或强制终止,但每一个规划失败落点还是会有一个最后的结果u(τ0),u(τ1),…,u(τN)输出,只是结果不能够满足终端或者静态约束条件,据此结果的具体情况来调整下一次优化的方向。
根据输出的可行方向重新选取可行的分离体落点信息,具体是:在优化后的控制变量序列u(τ0),u(τ1),…,u(τN)中如果攻角连续K个取值均达到kδmax(公式(1.3)中),则落点相对分离体初始位置沿分离体射向过远,即抬头飞-飞得远,则选取沿速度方向较规划失败落点近的落点,具体,K的数量根据值插值点的密度,如果密度大,K值选取多一些,如果一直处于最大值,则落点相对分离体初始位置沿射向过远,往近调;同样地,若攻角连续K个均达到-kδmax(公式(1.3)中),则落点相对分离体初始位置沿射向过近,则选取沿速度方向较规划失败落点远的落点;同样,如果侧滑角连续K个取值均达到kδmax或-kδmax,则落点相对分离体初始位置沿分离体速度垂向过远,则选取沿速度垂向较规划失败落点近的落点。分离体射向就是分离体的速度在水平面上投影指向,分离体速度垂向就是分离体的速度在水平面上投影指向的垂直方向。
另一方面,本发明还提供一种计算机可读的存储介质,所述计算机可读的存储介质包括存储的程序,其中,所述程序运行时执行前述的基于落点选择和伪谱法的分离体最优弹道规划的方法。
再者,本发明还提供一种电子装置,包括存储器和处理器,所述存储器中存储有计算机程序,所述处理器被设置为通过所述计算机程序执行前述的基于落点选择和伪谱法的分离体最优弹道规划的方法。

Claims (9)

1.一种基于落点选择和伪谱法的分离体最优弹道规划的方法,其特征在于包括如下步骤:
选取经验可行的分离体落点信息;
设计分离体弹道规划约束条件;
分离体弹道规划,根据分离体弹道规划约束条件采用Legendre伪谱法对分离体弹道进行求解优化;
将所述分离体落点信息输入分离体弹道进行求解优化,若得到一个可行弹道则将所述可行弹道输出,作为分离体飞行的标准弹道;若不能得到可行弹道,则重新获取可行的分离体落点信息直至得到可行弹道。
2.根据权利要求1所述的分离体最优弹道规划的方法,其特征在于所述选取经验可行的分离体落点信息包括
根据发射点、分离体分离点以及人物活动区域的位置关系,要求分离体落点散布在人、物频繁活动区域周边,选取经验可行的分离体落点n个,所述n个分离体落点对应n条分离体弹道,所述n为自然数1,2......i。
3.根据权利要求2所述的分离体最优弹道规划的方法,其特征在于所述每一条分离体弹道规划中设计分离体约束条件包括
获取分离体总体参数,包括质量和质心参数;
气动和载荷约束,根据分离体通过气动力做弱机动进而改变落点的特点,在规划弹道时利用满足要求的攻角和侧滑角作为控制变量得到在分离体上产生的横法向过载值,所述横、法向过载要求其随时间的积分最小,为此选取如下二次型目标泛函:
Figure FDA0003927570880000021
由于分离体在飞行过程中所受到的横法向过载均由气动力产生,且与攻角、侧滑角的大小直接相关,因此式等效为如下形式:
Figure FDA0003927570880000022
式中:
k1,k2,kα,kβ——人为确定的系数;
ny——横向过载;
nz——法向过载;
α——攻角(不考虑风)
β——侧滑角(不考虑风);
所述满足要求的攻角和侧滑角跟俯仰舵偏角和偏航舵偏角的关系即分离体所受约束条件如下:
Figure FDA0003927570880000023
式中
q=0.5ρV2——弹道动压
s——参考面积
l——参考长度
Cm——俯仰力矩系数
Cn——偏航力矩系数
Ma——马赫数
δp——俯仰舵偏角
δr——偏航舵偏角
δmax——最大允许舵偏角
k——舵效限制(一般小于1,需要留取一定舵偏余量用于姿态调整);
分离体落点要求满足约束如下:
Figure FDA0003927570880000031
式中
tf——着陆时间(规划前不确定)
(x,y,z)——落点在发射系中的坐标
i——落点编号,Xi,Yi,Zi是可行落点的坐标,在发射系中的坐标动力学方程构建。一般分离体质心运动方程可表示如下:
Figure FDA0003927570880000032
式中
V——分离体速度,初值为分离时刻标准弹道速度,向量;
X——分离体在发射系中的坐标;
m——分离体质量。
4.根据权利要求3所述的分离体最优弹道规划的方法,其特征在于所述根据分离体约束条件对分离体弹道进行求解优化包括
记控制变量(向量)为
Figure FDA0003927570880000044
记状态变量(向量)为
Figure FDA0003927570880000045
将静态约束公式(1.3)的内容记为向量g(s,u)≤0;将终端约束公式(1.4)的内容记为向量hi(s)=0;对应第i个落点的目标泛函公式(1.2)记为Ji(u),将动力学与运动学方程公式(1.5)记为
Figure FDA0003927570880000041
那么用于规划对应第i个落点的优化问题可按下列提法给出
Figure FDA0003927570880000042
记分离时刻为t0=0,着陆时间为tf,做如下微分变换:
Figure FDA0003927570880000043
利用Legendre伪谱法选择离散点τ01,…,τN,其中τ0=-1,τN=1,那么状态量可由插值点近似表示为
Figure FDA0003927570880000051
上式中s对τ求导可得
Figure FDA0003927570880000052
并将目标泛函转化为离散形式后,可以给出弹道规划问题的提法如下:
Figure FDA0003927570880000053
由上,弹道规划优化问题已经转变为关于控制变量序列u(τ0),u(τ1),…,u(τN)的最优化问题,通过一般的优化算法方便求解,终止条件为优化算法收敛或者运行时间超过一定值。
5.根据权利要求4所述的分离体最优弹道规划的方法,其特征在于所述若得到一个可行弹道则将所述可行弹道输出,作为分离体飞行的标准弹道,包括所述n个分离体落点对应n条分离体弹道分别采用Legendre伪谱法对分离体弹道进行求解优化,若得到任何一条可行弹道则终止计算,将所述可行弹道输出,作为分离体飞行的标准弹道。
6.根据权利要求4所述的分离体最优弹道规划的方法,其特征在于所述若不能得到可行弹道,则重新获取可行的分离体落点信息直至得到可行弹道,包括
所述n个分离落点在现有约束条件下均不能得到可行弹道,则为规划失败落点;根据优化后的控制变量序列u(τ0),u(τ1),…,u(τN)的取值输出拟可行的分离体落点相对规划失败落点的可行方向,根据输出的可行方向重新获取可行的分离体落点信息,重新进行弹道规划,直至得到满足要求的分离体弹道为止。
7.根据权利要求6所述的分离体最优弹道规划的方法,其特征在于所述根据输出的可行方向重新选取可行的分离体落点信息,包括
在优化后的控制变量序列u(τ0),u(τ1),…,u(τN)中如果攻角连续K个取值均达到kδmax,则落点相对分离体初始位置沿速度指向过远,则选取沿速度方向较规划失败落点近的落点;同样地,若攻角连续K个均达到-kδmax,则落点相对分离体初始位置沿速度指向过近,则选取沿速度方向较规划失败落点远的落点;同样,如果侧滑角连续K个取值均达到kδmax或-kδmax,则落点相对分离体初始位置沿速度垂向过远,则选取沿速度垂向较规划失败落点近的落点。
8.一种计算机可读的存储介质,其特征在于,所述计算机可读的存储介质包括存储的程序,其中,所述程序运行时执行上述权利要求1至7任一项中所述的方法。
9.一种电子装置,包括存储器和处理器,其特征在于,所述存储器中存储有计算机程序,所述处理器被设置为通过所述计算机程序执行所述权利要求1至7任一项中所述的方法。
CN202211379122.7A 2022-11-04 2022-11-04 一种基于落点选择和伪谱法的分离体最优弹道规划的方法 Pending CN115629618A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211379122.7A CN115629618A (zh) 2022-11-04 2022-11-04 一种基于落点选择和伪谱法的分离体最优弹道规划的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211379122.7A CN115629618A (zh) 2022-11-04 2022-11-04 一种基于落点选择和伪谱法的分离体最优弹道规划的方法

Publications (1)

Publication Number Publication Date
CN115629618A true CN115629618A (zh) 2023-01-20

Family

ID=84909233

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211379122.7A Pending CN115629618A (zh) 2022-11-04 2022-11-04 一种基于落点选择和伪谱法的分离体最优弹道规划的方法

Country Status (1)

Country Link
CN (1) CN115629618A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116362163A (zh) * 2023-06-01 2023-06-30 西安现代控制技术研究所 一种非奇异多约束弹道快速优化方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2612029C1 (ru) * 2016-02-09 2017-03-02 Акционерное общество "Федеральный научно-производственный центр "Нижегородский научно-исследовательский институт радиотехники" Способ распознавания баллистических целей и определения координат точек их пуска и падения для обзорных радиолокационных станций
CN110822994A (zh) * 2019-10-25 2020-02-21 北京航空航天大学 一种带落角约束的线性伪谱散布控制制导方法
CN113154955A (zh) * 2020-12-28 2021-07-23 航天科工火箭技术有限公司 一种自旋稳定的火箭分离体残骸落区精准控制系统及方法
CN114036727A (zh) * 2021-10-22 2022-02-11 湖北航天技术研究院总体设计所 一种基于末制导的滑翔弹道快速规划方法
CN114491811A (zh) * 2022-01-20 2022-05-13 航天科工火箭技术有限公司 一种运载火箭分离体的弹道设计方法
CN114624999A (zh) * 2022-01-21 2022-06-14 航天科工火箭技术有限公司 一种固体火箭一级分离体落区控制系统及方法
CN116817676A (zh) * 2022-12-09 2023-09-29 航天科工火箭技术有限公司 一种运载火箭分离体落区自动规划方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2612029C1 (ru) * 2016-02-09 2017-03-02 Акционерное общество "Федеральный научно-производственный центр "Нижегородский научно-исследовательский институт радиотехники" Способ распознавания баллистических целей и определения координат точек их пуска и падения для обзорных радиолокационных станций
CN110822994A (zh) * 2019-10-25 2020-02-21 北京航空航天大学 一种带落角约束的线性伪谱散布控制制导方法
CN113154955A (zh) * 2020-12-28 2021-07-23 航天科工火箭技术有限公司 一种自旋稳定的火箭分离体残骸落区精准控制系统及方法
CN114036727A (zh) * 2021-10-22 2022-02-11 湖北航天技术研究院总体设计所 一种基于末制导的滑翔弹道快速规划方法
CN114491811A (zh) * 2022-01-20 2022-05-13 航天科工火箭技术有限公司 一种运载火箭分离体的弹道设计方法
CN114624999A (zh) * 2022-01-21 2022-06-14 航天科工火箭技术有限公司 一种固体火箭一级分离体落区控制系统及方法
CN116817676A (zh) * 2022-12-09 2023-09-29 航天科工火箭技术有限公司 一种运载火箭分离体落区自动规划方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
LINGXIA MU 等: "Landing Trajectory Generation using Gauss Pseudo-spectral Method", IEEE, 31 December 2021 (2021-12-31) *
李惠峰;葛亚杰;李昭莹;: "高超声速飞行器自适应BTT末制导律", 北京航空航天大学学报, no. 05, 7 May 2013 (2013-05-07) *
梁栋;刘良栋;何英姿;: "月球精确软着陆最优标称轨迹在轨制导方法", 中国空间科学技术, no. 06, 25 December 2011 (2011-12-25) *
黄诘 等: "基于Radau伪谱方法的轨迹优化", 海军航空工程学院学报, vol. 30, no. 4, 31 December 2015 (2015-12-31) *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116362163A (zh) * 2023-06-01 2023-06-30 西安现代控制技术研究所 一种非奇异多约束弹道快速优化方法
CN116362163B (zh) * 2023-06-01 2023-09-01 西安现代控制技术研究所 一种多约束弹道快速优化方法

Similar Documents

Publication Publication Date Title
CN111306998B (zh) 一种参数摄动自适应的制导火箭弹垂直攻击制导方法
CN109508030B (zh) 一种考虑多禁飞区约束的协同解析再入制导方法
Azinheira et al. A backstepping controller for path‐tracking of an underactuated autonomous airship
CN111580547B (zh) 一种高超声速飞行器编队控制方法
CN112550770B (zh) 一种基于凸优化的火箭软着陆轨迹规划方法
CN111336871B (zh) 一种基于迂回式飞行的垂直攻击制导方法
JPH09506962A (ja) 発射前、発射後動作制御装置を有するへリコプタ一体型発射及び飛行制御
CN110425943B (zh) 面向变质心飞行器的工程化再入制导方法
CN110032206B (zh) 远程制导飞行器大落角攻顶控制方法及控制系统
CN115629618A (zh) 一种基于落点选择和伪谱法的分离体最优弹道规划的方法
CN112507461B (zh) 一种运载火箭动力软着陆段发动机开机方法
CN107977009A (zh) 一种考虑耦合的吸气式飞行器姿态控制律设计方法
CN112182772A (zh) 火箭推进控制方法、设备及存储介质
CN116991170B (zh) 短距起飞和垂直降落飞机的降落阶段自适应控制设计方法
CN114200827A (zh) 一种超音速大机动靶标的多约束双通道控制方法
CN113504723B (zh) 一种基于逆强化学习的运载火箭减载控制方法
Yang Analysis of optimal midcourse guidance law
Lu et al. Abort guidance during powered descent for crewed lunar missions
CN111897219B (zh) 基于在线逼近器的倾转四旋翼无人机过渡飞行模式最优鲁棒控制方法
CN112629339A (zh) 一种基于直接法的火箭软着陆轨迹规划方法
CN112660426A (zh) 一种火箭软着陆制导方法
CN113359861B (zh) 一种无人飞艇编队飞行控制方法及系统
CN115657458A (zh) 一种基于能量匹配的飞行器爬升轨迹控制方法
CN112327905B (zh) 一种基于直接升力的空中加油对接飞行控制方法
CN114237299A (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