网站建设合同服务响应时间网站底部友情链接
多弹协同末制导律设计
摘要
本文针对多弹协同任务中的导弹末制导问题,提出了一种基于相对运动模型的末制导律设计方法。该方法在传统比例导航制导基础上引入多弹协同修正因子,以确保各导弹在终端拦截过程中既能保持对目标的精准跟踪,又能协调避免相互干扰。文中详细推导了视线角变化率、传统比例导航制导律以及多弹协同修正部分的数学模型。
关键词
多弹协同;末制导律;比例导航;相对运动模型;协同控制
1. 引言
随着作战环境的复杂化和目标机动能力的提升,单枚导弹作战逐渐向多弹协同作战转变。传统的末制导律如比例导航(Proportional Navigation, PNG)在单弹拦截中虽然成熟,但在多弹协同场景下,导弹间可能出现交叉干扰和碰撞风险。因此,如何在保留单兵制导优点的基础上,实现导弹之间的有效协同与避碰,成为一个亟待解决的技术难题。本文基于相对运动模型,提出一种改进的多弹协同末制导律,并针对视线角及其变化率进行了细致的公式推导。
2. 系统模型与基本假设
考虑一枚导弹与目标之间的相对运动,设导弹位置为 P m = [ x m , y m , z m ] T \mathbf{P}_m = [x_m, y_m, z_m]^T Pm=[xm,ym,zm]T,目标位置为 P t = [ x t , y t , z t ] T \mathbf{P}_t = [x_t, y_t, z_t]^T Pt=[xt,yt,zt]T。在末制导阶段,假设目标与导弹主要在平面内运动,可忽略高度( z z z方向)的变化,则相对位置向量可写为:
R = P t − P m = [ R x , R y ] T . \mathbf{R} = \mathbf{P}_t - \mathbf{P}_m = [R_x, R_y]^T. R=Pt−Pm=[Rx,Ry]T.
导弹与目标的速度分别记为 v m = [ v m , x , v m , y ] T \mathbf{v}_m = [v_{m,x}, v_{m,y}]^T vm=[vm,x,vm,y]T 和 v t = [ v t , x , v t , y ] T \mathbf{v}_t = [v_{t,x}, v_{t,y}]^T vt=[vt,x,vt,y]T,它们的相对速度为:
v r e l = v t − v m . \mathbf{v}_{rel} = \mathbf{v}_t - \mathbf{v}_m. vrel=vt−vm.
3. 导弹末制导律设计
3.1 视线角及其变化率的推导
定义导弹与目标之间的视线角 λ \lambda λ:
λ = arctan ( R y R x ) . \lambda = \arctan\left(\frac{R_y}{R_x}\right). λ=arctan(RxRy).
对 λ \lambda λ 关于时间求导,应用链式法则可得视线角的变化率:
λ ˙ = R x R ˙ y − R y R ˙ x R x 2 + R y 2 . \dot{\lambda} = \frac{R_x \dot{R}_y - R_y \dot{R}_x}{R_x^2 + R_y^2}. λ˙=Rx2+Ry2RxR˙y−RyR˙x.
其中, R ˙ x \dot{R}_x R˙x 和 R ˙ y \dot{R}_y R˙y 分别为相对距离分量的时间导数,有
R ˙ x = v t , x − v m , x , R ˙ y = v t , y − v m , y . \dot{R}_x = v_{t,x} - v_{m,x}, \quad \dot{R}_y = v_{t,y} - v_{m,y}. R˙x=vt,x−vm,x,R˙y=vt,y−vm,y.
3.2 传统比例导航制导律(PNG)
传统PNG制导律提出,加速度指令与视线角变化率成正比:
a P N G = N V m λ ˙ , a_{PNG} = N V_m \dot{\lambda}, aPNG=NVmλ˙,
其中:
- N N N 为比例导航增益;
- V m = ∥ v m ∥ V_m = \|\mathbf{v}_m\| Vm=∥vm∥ 为导弹的速度大小。
3.3 多弹协同修正设计
在多弹协同作战中,除了单弹制导要求外,还需要考虑导弹间的协同配合与避碰,防止拦截过程中出现交叉干扰。为此,设计引入协同修正因子。假设系统共有 k k k 枚导弹,对第 i i i 枚导弹,定义协同因子 α i \alpha_i αi 为:
α i = 1 + β ∑ j = 1 j ≠ i k 1 d i j , \alpha_i = 1 + \beta \sum_{\substack{j=1 \\ j \neq i}}^{k} \frac{1}{d_{ij}}, αi=1+βj=1j=i∑kdij1,
其中:
- d i j = ∥ P m i − P m j ∥ d_{ij} = \|\mathbf{P}_{m_i} - \mathbf{P}_{m_j}\| dij=∥Pmi−Pmj∥ 表示导弹 i i i 与导弹 j j j 的距离;
- β \beta β 为用于调节协同影响程度的正系数。
为了引入导弹间的斥力作用,定义修正斥力项:
f r e p = K ∑ j = 1 j ≠ i k P m i − P m j d i j 2 , f_{rep} = K \sum_{\substack{j=1 \\ j \neq i}}^{k} \frac{\mathbf{P}_{m_i} - \mathbf{P}_{m_j}}{d_{ij}^2}, frep=Kj=1j=i∑kdij2Pmi−Pmj,
其中 K K K 为常数,用于确定斥力大小。
因此,导弹 i i i 的最终加速度指令可写为:
a c i = N V m i λ i ˙ + α i ⋅ f r e p . a_{c_i} = N V_{m_i} \dot{\lambda_i} + \alpha_i \cdot f_{rep}. aci=NVmiλi˙+αi⋅frep.
3.4 制导律完整推导过程
综合上述分析,设计流程如下:
-
计算相对距离与视线角:
导弹与目标的相对位置:
R = [ R x , R y ] T = P t − P m , \mathbf{R} = [R_x, R_y]^T = \mathbf{P}_t - \mathbf{P}_m, R=[Rx,Ry]T=Pt−Pm,
并计算视线角:
λ = arctan ( R y R x ) . \lambda = \arctan\left(\frac{R_y}{R_x}\right). λ=arctan(RxRy). -
求解视线角变化率:
利用导弹和目标速度分量:
λ ˙ = R x ( v t , y − v m , y ) − R y ( v t , x − v m , x ) R x 2 + R y 2 . \dot{\lambda} = \frac{R_x (v_{t,y} - v_{m,y}) - R_y (v_{t,x} - v_{m,x})}{R_x^2 + R_y^2}. λ˙=Rx2+Ry2Rx(vt,y−vm,y)−Ry(vt,x−vm,x). -
基础加速度指令(PNG):
按照传统PNG制导律,得基础加速度指令:
a P N G = N V m λ ˙ . a_{PNG} = N V_m \dot{\lambda}. aPNG=NVmλ˙. -
计算协同修正因子:
对于多弹系统,引入协同因子:
α i = 1 + β ∑ j = 1 j ≠ i k 1 ∥ P m i − P m j ∥ . \alpha_i = 1 + \beta \sum_{\substack{j=1 \\ j \neq i}}^{k} \frac{1}{\|\mathbf{P}_{m_i} - \mathbf{P}_{m_j}\|}. αi=1+βj=1j=i∑k∥Pmi−Pmj∥1. -
计算斥力修正项:
定义导弹间的斥力:
f r e p = K ∑ j = 1 j ≠ i k P m i − P m j ∥ P m i − P m j ∥ 2 . f_{rep} = K \sum_{\substack{j=1 \\ j \neq i}}^{k} \frac{\mathbf{P}_{m_i} - \mathbf{P}_{m_j}}{\|\mathbf{P}_{m_i} - \mathbf{P}_{m_j}\|^2}. frep=Kj=1j=i∑k∥Pmi−Pmj∥2Pmi−Pmj. -
综合得最终加速度指令:
最终对第 i i i 枚导弹,修正后的加速度指令为:
a c i = N V m i λ i ˙ + α i ⋅ f r e p . a_{c_i} = N V_{m_i} \dot{\lambda_i} + \alpha_i \cdot f_{rep}. aci=NVmiλi˙+αi⋅frep.
该制导律在保证单导弹末制导精度的基础上,通过协同修正有助于实现多导弹间的避碰与协同拦截。