带扩张观测器的新型滑模导引律
2015-04-22王华吉简金蕾雷虎民马卫东
王华吉,简金蕾,雷虎民,李 解,马卫东
(空军工程大学 防空反导学院,西安 710051)
带扩张观测器的新型滑模导引律
王华吉,简金蕾,雷虎民,李 解,马卫东
(空军工程大学 防空反导学院,西安 710051)
在精确制导问题中,为克服目标机动和弹体动态特性对制导精度的影响,建立了平面内的弹目相对运动模型,在此基础上建立考虑导弹动态特性的制导模型;为提高末制导精度,设计了考虑导弹动态特性和目标机动的自适应滑模导引律;为了实现该导引律,利用带有滤波器的扩张观测器估计视线角速率、视线角加速度、目标机动加速度及其变化率等制导信息。仿真结果表明,扩张观测器收敛速度快、估计精度高,且具有较强的抗干扰能力;在不同机动条件下,所设计的考虑导弹动态特性的含扩张观测器的改进滑模控制律相比于比例导引律、增广比例导引律和滑模导引律具有较好的导引性能。
滑模导引律;扩张观测器;自动驾驶仪;导引精度
0 引言
在未来战争中,高速机动目标将成为现有导弹防御体系的主要威胁[1],为迎接挑战,越来越多的国家重视对具有拦截高速机动目标能力的新一代寻的导弹研制,而提高末制导精度成为关键问题。变结构控制因具有抗干扰和抗参数摄动的特性,使得它在导引律设计中得到广泛应用。而目标的逃逸机动和弹体动态特性是影响制导性能的重要因素[2],一般的变结构导引律在设计中没有考虑弹体动态特性[3];同时,也将目标机动视为外界干扰来处理,从而因切换项参数选择不当造成抖振现象。
为提高制导性能,本文设计了考虑弹体动态特性和目标机动的变结构导引律。由于导引头无法测量目标机动,所以需要对目标机动加速度进行估计。而观测器在制导信息的估计中得到很好地应用。李雅静等[4]在导弹视线角和视线角速率可测情况下,提出了多状态相关系数矩阵加权组合的方案,对CB观测器进行了改进,用于估计弹目相对距离、相对速度和目标机动加速度。姚郁等[5]将目标机动加速度当作不确定性扩张成新的一阶状态,设计了扩张状态观测器来观测系统状态和估计目标加速度。马克茂等[6]针对末制导过程中制导信息的获取问题,设计了高增益观测器对视线角速率和目标机动加速度进行了估计,取得很好的估计效果。扩张观测器(extended state observer,ESO)是一种非线性状态观测器,是自抗扰控制的重要组成部分。通过把系统中的内外扰动扩张成系统新的一阶状态,再利用特定的非光滑非线性误差反馈,适当选择观测器参数,得到系统所有状态的观测值[5]。
本文采用雷达导引头,假设弹目距离、弹目距离变化率、视线角和视线角速率可测量得到。将系统中的机动目标加速度当作不确定性,并扩张成新的一阶状态,设计了二阶扩张状态观测器来观测系统状态;采用带有滤波器的扩张状态观测器来估计制导信息,然后对所设计的新型滑模制导律进行了实现。仿真结果显示,扩张观测器估计精度高、收敛速度快、具有较强的抗干扰能力,并将估计值应用到制导律的实现,取得了不错的效果。
1 导弹-目标相对运动关系
为研究方便,假设导弹在飞行过程中不发生滚转,将导弹在三维空间中的运动解耦成纵向平面运动和侧向平面运动。由于侧向平面的运动与纵向平面的运动类似,本文仅针对纵向平面拦截情况进行分析。纵向平面的弹目相对运动关系如图1所示。图1中,Vm、Vt分别为导弹与目标的速度;am、at分别为导弹和目标的加速度大小;θm、θt分别为导弹和目标的弹道倾角;R为弹目相对距离;q为弹目视线角;矢量iL、jL为惯性坐标系上的单位方向矢量。
由图1可得,弹目相对运动学模型:
(1)
(2)
(3)
(4)
其中
(5)
(6)
2 考虑导弹动态特性的滑模制导律设计
导弹自动驾驶仪十分复杂,假如将其直接引入到导弹制导律的设计中,将使制导律的设计变得过于复杂而无法完成。为了兼容制导律的设计难度和真实情况的逼近度,本文将导弹自动驾驶仪简化为一阶惯性环节。
2.1 制导模型推导
若将导弹自动驾驶仪考虑成一阶惯性环节,则导弹制导指令和实际过载关系为
(7)
将上式变换到时域有
(8)
又由于
(9)
为便于分析,综合式(6)、式(8)和式(9)得到
(10)
对式(10)第2式求导得
(11)
由式(8)和式(11)得
(12)
(13)
在视线坐标系下,导弹和目标的相对加速度为
a=(atR-amR)iL+(atq-amq)jL
(14)
对式(14)求导得
(15)
(16)
(17)
(18)
由式(13)得
(19)
将式(19)写成矩阵形式为
(20)
其中
2.2 新型滑模导引律的设计
(21)
(22)
选取滑模面为
S=CX,C=[c1c2]
(23)
为保证系统状态能到达变结构,且到达的过程中有优良的动态特性,本文采用趋近律方式推导导引律。由于式(22)所示系统为时变系统,所以构造了对时变参数具有自适应能力的变结构趋近律,即
对式(23)求导得
(24)
则
(25)
将上面a1,a2,f和c2=1代入式(25)得
(26)
(27)
本文采用高增益连续函数S/[|S|+δ]代替符号函数sgn(S),实现准滑动模态控制[11],可有效地消弱抖动。
3 带滤波器的扩张观测器设计
3.1 扩张观测器
(28)
其中,Vq为系统状态量;amq为控制量;atq为系统中的未知干扰;y为系统的可测输出,若将atq作为系统的扩张状态,则新系统为[7]
(29)