导弹稳定控制系统参数自动化设计算法应用研究

(整期优先)网络出版时间:2022-03-16
/ 2

导弹稳定控制系统参数自动化设计算法应用研究

韦汉林 韦世顺

江南机电设计研究所 贵州贵阳 550000

江南机电设计研究所 550000

摘要:针对目前导弹控制系统参数设计时完全依靠设计师主观经验,对所选取特征点的控制系统参数进行逐个‘试凑’整定的耗时、主观、不规范设计问题,本文研究了将粒子群算法(PSO)应用于导弹稳定控制回路参数自动化设计的方法,且针对标准粒子群算法在导弹稳定控制回路参数寻优中表现出收敛快、成功率低的问题,通过引入遗传算法中选择、交叉思想对标准粒子群算法进行了改进,形成粒子群-遗传算法(PSO-GA),扩大了解空间范围,成功解决了标准粒子群算法存在的缺点。通过对PSO算法及PSO-GA算法的研究,实现了导弹稳定控制回路的规范化、自动化设计,具有一定的实用价值。

主题词:粒子群算法 粒子群-遗传算法 稳定控制回路 自动化设计

  1. 引言

导弹稳定控制系统的设计,一般是在整个飞行空域内取一定数量的特征点,然后根据线性化数学模型,对每个特征点逐一进行控制回路参数设计。目前工程设计中,多采用人工‘试凑’的方法,进行多次循环反复设计和校验,直到凑出满足设计指标要求的控制器参数。这种依赖设计师工程经验,参杂了人为主观因素在内的设计,不仅不规范且耗时,若选取的特征点数量太多,仍然采用传统的人工设计方法,是一件非常繁杂且不可完成的任务。

本文研究了粒子群算法及其改进算法在导弹侧向稳定控制回路参数自动化设计中的应用,分别将PSO算法和PSO-GA算法用于某型导弹侧向稳定控制回路自动化设计,并将设计结果在六自由度中进行仿真验证,并与传统人工设计结果进行比较,证实了粒子群优化算法及其改进算法在导弹侧向稳定控制回路中设计的可行性。

  1. 粒子群算法及其改进

    1. 标准PSO算法

设在623159f8b1cfa_html_2516832871563a24.gif 维目标搜索空间中有623159f8b1cfa_html_de7c3ab3314fa66a.gif 个粒子,粒子623159f8b1cfa_html_e154004f1e0ffa3b.gif 的空间位置为623159f8b1cfa_html_3c9012cf55f97624.gif ,将623159f8b1cfa_html_a2a0385d7e509e80.gif 代入目标函数就可以计算出其适应度值,根据适应度值的大小衡量粒子623159f8b1cfa_html_4a7222e84f628e8d.gif 的优劣。单个粒子所经历过的最优位置记为623159f8b1cfa_html_2c5cb4ee096cb56d.gif ,整个粒子群经历过的最优位置记作623159f8b1cfa_html_8ed866626a2289d8.gif ,粒子根据公式(1)来更新自己的速度和位置:

623159f8b1cfa_html_f6add0dd97f65081.gif (1)

式中623159f8b1cfa_html_5de5b8f0d2effcc.gif 为加权系数也称动力常量,为非负数,代表了前一速度对当前速度的影响,本文采用自适应调整方法进行623159f8b1cfa_html_5de5b8f0d2effcc.gif 设计,如式(2);623159f8b1cfa_html_4b081bca12b1c053.gif623159f8b1cfa_html_a83b29c22c65fd9d.gif 为学习因子见式(3);623159f8b1cfa_html_e3f5abfb887a0812.gif 为第623159f8b1cfa_html_4df45dacced3ac26.gif 个粒子的飞翔速度,623159f8b1cfa_html_56eafe9b2f58e1c9.gif623159f8b1cfa_html_c7d9975dede8818e.gif 为第623159f8b1cfa_html_a20a9da353b8776e.gif 个粒子的位置;623159f8b1cfa_html_1bd9c41c3547853.gif 为相互独立的随机数,服从623159f8b1cfa_html_559d7d33a8e635e2.gif 的均匀分布。

623159f8b1cfa_html_810793f8f0e78a0f.gif (2)

623159f8b1cfa_html_30a630910e42ef2d.gif (3)

式中623159f8b1cfa_html_5da0a2ee835add5e.gif 分别为当前迭代次数和最大迭代次数。

由于导弹侧向稳定控制回路参数设计具有多约束条件,计算结果表明:标准粒子群算法在对这种多约束条件下的复杂系统进行控制参数寻优时,由于其‘早熟’收敛特性,容易陷入局部极值,成功率低。

    1. 粒子群-遗传算法(PSO-GA)

遗传算法中交叉是以二进制码的形式进行,这样的操作需要进行二进制码生成和二进制解码操作,操作精度高,但额外增加了计算量。本文直接采用十进制数进行交叉操作,如式(4)。

623159f8b1cfa_html_fb8fec295eed5fa8.gif (4)

式中,623159f8b1cfa_html_745aa3f50f2925a8.gif 为标准PSO寻优结果,作为两个父代粒子,623159f8b1cfa_html_732811090adf041d.gif623159f8b1cfa_html_7133e6667747ae94.gif 为经过选择的两个父代粒子进行交叉后得到的两个新粒子,623159f8b1cfa_html_ce3f5c644779ec33.gif623159f8b1cfa_html_645ae8e1c973ccc.gif 间的随机数。

  1. 系统模型与目标函数

    1. 侧向稳定控制回路结构

本文采用复合稳定回路的结构形式,如图1:

623159f8b1cfa_html_7564ac778745571.gif

  1. 稳定控制回路结构图

图中623159f8b1cfa_html_98058817e2458155.gif 为控制指令,623159f8b1cfa_html_c1fb6ea132983066.gif 为过载响应,623159f8b1cfa_html_b3bcb15ca9b2cba5.gif 为侧向角速度,623159f8b1cfa_html_3539b9089ce088df.gif623159f8b1cfa_html_992ab56d93e73459.gif623159f8b1cfa_html_f48c2b03f64b0fd9.gif623159f8b1cfa_html_1b1d47919cb66961.gif 为调参参数。

舵偏角623159f8b1cfa_html_a3c8afb814d50ac9.gif 到过载的传递函数为:

623159f8b1cfa_html_31e30a4adcc95633.gif (5)

舵偏角623159f8b1cfa_html_a3c8afb814d50ac9.gif 到角速度的传递函数为:

623159f8b1cfa_html_e9fdea8d3b5c8020.gif (6)

三阶期望特征多项式为:

623159f8b1cfa_html_fce9d44c059e3665.gif (7)

其中,参数623159f8b1cfa_html_5002573e7f41d846.gif623159f8b1cfa_html_8fc8767c7f683315.gif623159f8b1cfa_html_7654ec667dfc88c2.gif 为设计参数,本文利用PSO算法和PSO-GA算法进行自动寻优,从而解出调参系数623159f8b1cfa_html_3539b9089ce088df.gif623159f8b1cfa_html_992ab56d93e73459.gif623159f8b1cfa_html_f48c2b03f64b0fd9.gif623159f8b1cfa_html_1b1d47919cb66961.gif

    1. 目标函数

本文定义包含侧向稳定控制回路单位阶跃过载响应的超调量、上升时间、稳态误差、幅值裕度、相位裕度、侧向角速、舵偏角、舵偏角速度,来评价侧向稳定控制系统的综合性能,得到下述目标函数:

623159f8b1cfa_html_4ebf52aa6c613342.gif (8)

式中,623159f8b1cfa_html_77800102e213733b.gif 是过载响应上升时间,623159f8b1cfa_html_756041bf01ee0f37.gif 是过载响应的超调量,623159f8b1cfa_html_11a0d95b5e79d618.gif 为过载响应稳态误差,623159f8b1cfa_html_e3bbcaf66acf5246.gif 分别为阻尼回路幅裕度、幅裕度期望值、阻尼回路相裕度、相裕度期望值、阻尼回路穿越频率、阻尼回路穿越频率最大允许值、加速度回路幅裕度、幅裕度期望值、加速度回路相裕度、相裕度期望值,623159f8b1cfa_html_e27a2bca7306503e.gif623159f8b1cfa_html_db09e446c8182a0f.gif 分别最大侧向角速度、最大舵偏角、最大舵偏角速度及分别对应的指标值,623159f8b1cfa_html_74160a76f187794d.gif 是加权系数,可进行适当调整。

  1. 仿真验证

以某一特征点为例,分别采用PSO和PSO-GA完成对导弹侧向稳定控制回路控制器参数的自化设计。

    1. 算法初始化

算法迭代次数623159f8b1cfa_html_48b69eebfc2bcb79.gif :120;种群规模:30;粒子位置上下边界值623159f8b1cfa_html_88346cbdb6717148.gif623159f8b1cfa_html_8a1239dbcab7618b.gif ;粒子飞翔速度上界623159f8b1cfa_html_9be6a960d6bf9292.gif 。目标函数权重值:623159f8b1cfa_html_2fdb4852063a8aef.gif623159f8b1cfa_html_dc92144b0c272e2e.gif623159f8b1cfa_html_9648ea06186629a5.gif623159f8b1cfa_html_888826cae8b37ce3.gif623159f8b1cfa_html_6fa32d2c10bd2572.gif623159f8b1cfa_html_99d0d993519e8ded.gif ;裕度要求:623159f8b1cfa_html_7adb43c28fa15606.gif623159f8b1cfa_html_849471bc0393ca94.gif623159f8b1cfa_html_6fafc945c2fa873e.gif623159f8b1cfa_html_627455292561ca09.gif ,遗传交叉概率623159f8b1cfa_html_92239abe37a6a948.gif 。表1为该特征点参数情况:

  1. 特征点参数

序号

623159f8b1cfa_html_8e9b0708138e5976.gif

623159f8b1cfa_html_b4947d97c596ba94.gif

623159f8b1cfa_html_a6327aad2dae546d.gif

623159f8b1cfa_html_60faaf99e032c543.gif

623159f8b1cfa_html_19a0fd653e079ab0.gif

623159f8b1cfa_html_7a8fdc536588105d.gif

623159f8b1cfa_html_8eeb7d1b96f10e35.gif

623159f8b1cfa_html_8bb99b1fbff1b6ea.gif

623159f8b1cfa_html_eeef9c21826f38e8.gif

2.49

258.7446

-203.0137

0.5850

0.0087

612.15

123638.8

1.94

9.7423

    1. 特征点定点设计仿真结果

表2为在PSO算法和PSO-GA算法下该特征点控制参数自动寻优设计结果。

  1. 两种算法的控制参数优化结果

序号

算法

623159f8b1cfa_html_e58a4a0274fc2f8b.gif

623159f8b1cfa_html_4932d2160a80d45c.gif

623159f8b1cfa_html_c6a11377cf21b864.gif

PSO

11.517465

0.73831717

0.026537256

PSO-GA

16.863501

0.69685713

0.024

在对该型导弹所选取特征点的定点设计中,PSO算法在某些特征点下出现最多的寻了4次才成功的情况,而PSO-GA算法均一次性成功。

  1. 两种算法的设计结果比较

序号

算法

幅值裕度(dB)

相裕度(°)

穿越频率(rad/s)

上升时间

(s)

超调量

(%)

侧向角速度最大值(°/s)

舵偏角速度最大值(°/s)

PSO

10.2344

58.6538

6.7085

0.2330

10.0004

97.3482

194.4088

PSO-GA

8.9423

59.8023

7.0947

0.2160

10.0001

107.4536

205.7907

从表3可知,标准PSO与PSO-GA算法参数寻优结果满足设计指标要求。图2为单位阶跃指令下过载响应曲线。

623159f8b1cfa_html_ce5987fbafcc0ab8.gif

  1. 单位阶跃响应曲线

    1. 六自由度仿真验证

为验证PSO算法和PSO-GA算法下自动设计得到的结果在导弹六自由度中同样适用,将设计结果带入到导弹六自由度模型中进行仿真,并与传统设计的仿真结果进行比较。

从图3可以看出,采用PSO算法和PSO-GA算法自动化寻优结果在六自由度仿真中的性能与传统人工设计结果相当,控制效果良好,可直接用于导弹稳定控制回路设计。

623159f8b1cfa_html_cb0ef4783958c4e5.gif

  1. 指令-过载响应曲线

  1. 结语

PSO算法和PSO-GA算法自动设计结果满足指标要求,在导弹六自由度中进行了仿真验证,控制效果良好,与传统设计结果相当,证实了本文算法在导弹侧向稳定控制回路设计中的可行性,能有效实现导弹稳定控制系统参数的自动化设计,具有一定的工程应用价值。

  1. 参考文献

[1] 张民,陈欣,陆宇平.基于改进PSO算法的导弹控制参数优化.南京航空航天大学学报,2009.

[2] 沈永福,李呈良,邓方林.基于改进遗传算法的导弹高度控制回路自动化设计.弹箭与制导学报,2002.

[3] 邱亚男,王建琦,贾晓洪.基于PSO算法的自动驾驶仪控制参数设计.航天控制,2009.