一种基于混沌模型与dnn的飞行器智能轨迹规划方法
技术领域
1.本发明涉及一种基于混沌模型与dnn的飞行器智能轨迹规划方法，属于飞行器轨迹规划技术领域。
背景技术：2.随着陆海空天探测技术的全面发展，现代飞行器越来越难以实现隐身飞行，发展在被探测条件下的新型突防技术是未来的必然趋势，而机动飞行是最重要也最有效的一种突防手段。
3.一方面，人工智能在飞行器中的应用日趋成熟。因此，可借助网络学习带来的精度和成本优势对飞行状态进行预测与分析，为突防决策和机动飞行轨迹规划提供必要的参考信息。
4.另一方面，传统的轨迹规划方法所产生的飞行弹道具有较强的规律性，现有的模糊预测与轨迹预报技术已经能够达到较高的精度。
5.因此，亟需一种能够突破现有轨迹规划模式的高混沌性与高敏捷性轨迹规划方法。
技术实现要素：6.为解决背景技术中存在的问题，本发明提供一种基于混沌模型与dnn的飞行器智能轨迹规划方法。
7.实现上述目的，本发明采取下述技术方案：一种基于混沌模型与dnn的飞行器智能轨迹规划方法，所述方法包括如下步骤：
8.s1：获得飞行状态量与飞行控制量之间的网络映射关系；
9.s2：将过程约束转化为飞行高度和倾侧角约束；
10.s3：通过chebyshev映射模型和飞行器模型测量偏差生成导引指令中的最终混沌信号δ；
11.s4：计算满足任务约束的倾侧角剖面中值，并计算出实时需要的倾侧角大小σc；
12.s5：将最终混沌信号δ和最大倾侧角变化量σ
max
(v)-σc叠加到实时计算得到的需要倾侧角指令上，得到混沌化处理后的倾侧角大小|σ|＝|σc+δ
·
(σ
max
(v)-σc)|；
13.s6：采用视线角误差走廊确定倾侧角方向；
14.s7：对倾侧角幅值和倾侧角速率幅值进行约束。
15.与现有技术相比，本发明的有益效果是：
16.本发明克服了传统飞行轨迹设计时的模型局限性，通过引入混沌模型系统和飞行器自身模型的不确定性，增加了导引指令和飞行轨迹的不可预测性；同时又通过智能神经网络对飞行状态的实时判断，限制了不确定性对任务可达性的影响。因此，本发明在保证预定任务可行性的前提下增加了飞行轨迹的混沌性，并在导引指令计算中满足过程约束，具有一定的工程实用价值。
附图说明
17.图1是本发明dnn网络训练结构图。
具体实施方式
18.下面将结合本发明实施例中的附图，对本发明中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是发明的一部分实施例，而不是全部的实施例，基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
19.为了更为科学准确的描述，给出混沌与轨迹规划的概念：
20.其中：
21.混沌：在非线性动态系统中出现的确定性的类似随机的过程。
22.轨迹规划：为满足飞行器空间位置等状态量约束条件下的飞行程序设计过程。
23.一种基于混沌模型与dnn的飞行器智能轨迹规划方法，所述方法包括如下步骤：
24.s1：获得飞行状态量与飞行控制量之间的网络映射关系；
25.s101：在初始偏差和环境扰动条件下利用飞行器动力学模型采用轨迹积分的形式生成飞行轨迹样本；
26.所述初始偏差和环境扰动包括：初始高度偏差(
±
1000m以内正态分布)、初始速度偏差(
±
100m/s以内正态分布)、初始飞行路径角偏差(
±1°
以内正态分布)、气动阻力系数偏差(
±
10％以内正态分布)、气动升力系数偏差(
±
10％以内正态分布)以及大气密度偏差(
±
10％以内正态分布)。
27.所述飞行器动力学模型为：
[0028][0029]
式(2)中：
[0030]
r为地心距；
[0031]
为飞行器飞行纬度；
[0032]
v为飞行器飞行速度；
[0033]
γ为飞行路径角；
[0034]
ψ为飞行器飞行航向角；
[0035]
θ为飞行器飞行经度；
[0036]
d为飞行器飞行气动阻力；
[0037]
ωe为地球自转角速率；
[0038]
l为飞行器飞行气动升力；
[0039]
σ为倾侧角。
[0040]
s102：以地心距r、速度v、飞行路径角γ、攻角剖面α(v)和剩余航程s等飞行状态量作为输入，以倾侧角剖面中值σ
mid
作为输出，采用dnn进行网络模型学习，获得飞行状态量与
飞行控制量之间的网络映射关系如下：
[0041]
σ
mid
＝net(r,v,γ,α(v),s)
ꢀꢀꢀ
(1)
[0042]
式(1)中：
[0043]
σ
mid
为倾侧角剖面中值；
[0044]
net(
·
)为学习完成的神经网络模型；
[0045]
r为地心距；
[0046]
v为飞行器飞行速度；
[0047]
γ为飞行路径角；
[0048]
α(v)为攻角剖面；
[0049]
s为剩余航程。
[0050]
所述dnn为双隐含层的深度神经网络，激活函数选择为relu函数，损失函数选择为均方差形式。
[0051]
s2：将过载、动压和驻点热流等过程约束转化为飞行高度和倾侧角约束；
[0052]
设定过程约束包括：最大动压约束q
max
、最大过载约束n
max
以及最大驻点热流约束q
max
，将过程约束转化为飞行高度和倾侧角约束的过程如下：
[0053][0054]
式(3)中：
[0055]hqmax
(v)为由最大动压约束转化而来的飞行高度约束；
[0056]hnmax
(v)为由最大过载约束转化而来的飞行高度约束；
[0057]hqmax
(v)为由最大驻点热流约束转化而来的飞行高度约束；
[0058]hs
为基准高度；
[0059]
ln(
·
)为以自然对数为底的对数函数；
[0060]
ρ0为海平面大气密度；
[0061]
v为飞行器飞行速度；
[0062]
l为飞行器飞行气动升力；
[0063]
d为飞行器飞行气动阻力；
[0064]cq
为热流密度系数；
[0065]
r0为地球平均半径；
[0066]
σ
max
(v)为倾侧角约束；
[0067]
σ
imax
(v)(i＝q,n,q)为对应转化后的倾侧角大小约束；
[0068]
q为动压约束；
[0069]
n为过载约束；
[0070]
q为驻点热流约束
[0071]himax
(v)为h
qmax
(v)、h
nmax
(v)和h
qmax
(v)的泛指。
[0072]
s3：通过chebyshev映射模型和飞行器模型测量偏差生成导引指令中的最终混沌信号δ；
[0073]
s301：所述chebyshev映射模型采用k阶的chebyshev映射模型如下：
[0074]
x
n+1
＝cos(k
·
arccosxn)
ꢀꢀꢀ
(4)
[0075]
式(4)中：
[0076]
xn为第n次计算得到的混沌信号；
[0077]
x
n+1
为第n+1次计算得到的混沌信号，定义区间为(-1,1)，n可按照仿真步长进行次数提取；
[0078]
s302：将上述混沌信号δ叠加至导引指令中，由映射模型与飞行器模型测量偏差组合得到：
[0079]
δ＝x
n+1
·
(1+w)
ꢀꢀꢀ
(5)
[0080]
式(5)中：
[0081]
δ为最终混沌信号；
[0082]
w为飞行器模型测量偏差，且
[0083][0084]
式(6)中：
[0085]
f(r,v,α)为根据飞行状态测量得到的飞行器前缘翼面压力；
[0086]fref
(r,v,α)为提前装订的标称条件下的飞行器前缘翼面压力；
[0087]fmax
为f
ref
(r,v,α)的最大值，仿真试验时将f(r,v,α)等效为气动阻力。
[0088]
s4：计算满足任务约束的倾侧角剖面中值，并计算出实时需要的倾侧角大小σc；
[0089]
s401：设定控制指令为攻角和倾侧角，其中，设定攻角与速度相关的函数如下：
[0090][0091]
式(7)中：
[0092]
α为飞行攻角；
[0093]
α
max
为控制热流约束的飞行允许最大攻角；
[0094]
v为飞行器飞行速度；
[0095]v1
和v2为分段处对应速度大小；
[0096]
α
l
为最大升阻比对应的攻角；
[0097]
将倾侧角设定为两线性段组合起来的与速度相关的函数如下：
[0098][0099]
式(8)中：
[0100]
σm为倾侧角；
[0101]
σ0为初始状态下按照准平衡滑翔条件算出来的倾侧角；
[0102]
σf为终端状态下按照准平衡滑翔条件算出来的倾侧角；
[0103]
σ
mid
为根据航程进行调整的倾侧角剖面中值；
[0104]v0
为指定的初始速度；
[0105]vf
为指定的终端速度；
[0106]vmid
为v0和vf的中值；
[0107]
通过上述即可求得倾侧角剖面中值σ
mid
；
[0108]
s402：由实时速度插值倾侧角剖面即可得到实时需要的倾侧角大小σc＝σ。
[0109]
s5：将最终混沌信号和最大倾侧角变化量σ
max
(v)-σc叠加到实时计算得到的需要倾侧角指令上，得到混沌化处理后的倾侧角大小|σ|＝|σc+δ
·
(σ
max
(v)-σc)|；
[0110]
s6：采用视线角误差走廊确定倾侧角方向；
[0111]
s6所述视线角误差走廊采用分段线性形式如下：
[0112][0113]
式(9)中：
[0114]
ψ
los
为飞行器与目标连线的视线角；
[0115]
θf为地面目标点的经度；
[0116]
θ为飞行器飞行经度；
[0117]
为飞行器飞行纬度；
[0118]
φf为地面目标点的纬度；
[0119]
δψ为视线角偏差；
[0120]
ψ为飞行器飞行航向角；
[0121]
δψ
t
为由视线角偏差的上下限确定的视线角偏差边界；
[0122]
δψh为设定的视线角偏差的上限；
[0123]
δψ
l
为设定的视线角偏差的下限；
[0124]
v为飞行器飞行速度；
[0125]v3
和v4为按照速度分段形式的视线角偏差走廊的分段速度点；
[0126]
σ
n+1
为当前周期的倾侧角；
[0127]
σn为上一周期的倾侧角；
[0128]
sgn(
·
)为符号函数。
[0129]
s7：对倾侧角幅值和倾侧角速率幅值进行约束。
[0130]
s7所述倾侧角幅值和倾侧角速率幅值采用如下形式：
[0131][0132]
式(10)中：
[0133]
σ
n+1
为当前周期的倾侧角；
[0134]
σn为上一周期的倾侧角；
[0135]
τs为制导周期；
[0136]
δσ
max
为倾侧角速率幅值约束值；
[0137]
sgn(
·
)为符号函数；
[0138]
σ
max
为倾侧角幅值约束值。
[0139]
本发明涉及一种面向突防任务需求的飞行器智能轨迹规划方法，通过混沌模型与深度神经网络使飞行器在满足任务可达性的同时，具备一定的轨迹混沌性。
[0140]
本发明采用混沌模型与神经网络生成飞行轨迹，提供一种既能满足终端任务约束又具有一定不可预测性的轨迹规划方法。
[0141]
对于本领域技术人员而言，显然本发明不限于上述示范性实施例的细节，而且在不背离本发明的精神或基本特征的情况下，能够以其他的形式实现本发明。因此，无论从哪一点来看，均应将实施例看作是示范性的，而且是非限制性的，本发明的范围由所附权利要求而不是上述说明限定，因此旨在将落在权利要求的等同条件的含义和范围内的所有变化囊括在本发明内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。
[0142]
此外，应当理解，虽然本说明书按照实施方式加以描述，但并非每个实施方式仅包含一个独立的技术方案，说明书的这种叙述方式仅仅是为清楚起见，本领域技术人员应当将说明书作为一个整体，各实施例中的技术方案也可以经适当组合，形成本领域技术人员可以理解的其他实施方式。