CN115755955B - 基于截获概率的空空导弹协同探测方法 - Google Patents
基于截获概率的空空导弹协同探测方法 Download PDFInfo
- Publication number
- CN115755955B CN115755955B CN202211362631.9A CN202211362631A CN115755955B CN 115755955 B CN115755955 B CN 115755955B CN 202211362631 A CN202211362631 A CN 202211362631A CN 115755955 B CN115755955 B CN 115755955B
- Authority
- CN
- China
- Prior art keywords
- target
- coordinate system
- probability
- projectile
- representing
- 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.)
- Active
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 61
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 50
- 238000005457 optimization Methods 0.000 claims abstract description 33
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 11
- 239000002245 particle Substances 0.000 claims abstract description 11
- 238000010586 diagram Methods 0.000 claims description 39
- 238000000034 method Methods 0.000 claims description 28
- 239000011159 matrix material Substances 0.000 claims description 20
- 230000001133 acceleration Effects 0.000 claims description 10
- 238000005259 measurement Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 238000005315 distribution function Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 claims description 4
- 230000005484 gravity Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000005312 nonlinear dynamic Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)
Abstract
基于截获概率的空空导弹协同探测方法,具体涉及一种基于目标截获概率的空空导弹协同探测的方法,本发明为了解决多拦截器探测时未给出截获概率,导致协同探测准确率低的问题,获取考虑目标机动时目标捕获时刻目标运动区域及其概率分布;选取一枚导弹作为领弹,四枚导弹作为从弹,根据目标运动区域计算领弹和从弹编队的队形参数;根据目标运动区域概率分布和队形参数计算二维平面内的协同截获概率;以队形参数为优化参数,截获概率为优化目标,利用粒子群优化算法优化协同截获概率,得到截获概率最大值,实现最佳的多弹协同探测。本发明属于多弹协同探测领域。
Description
技术领域
本发明涉及一种协同探测方法,具体涉及一种基于目标截获概率的空空导弹协同探测的方法,属于多弹协同探测领域。
背景技术
现代战争中随着复杂的战场环境及目标强机动特性等因素的影响,空空导弹不仅要求自身的机动性能和精确制导性能,还需完成多平台间的协同作战,以实现整体作战效能的最大化。在多平台协同作战中,可以有效扩大侦查范围、提高探测精度、增强能力,并为协同末制导提供良好的中末交接班条件。现有多拦截器(多拦截器包括空空导弹,除了空空导弹之外还包括多无人机、多战斗机等)中,选取编队队形和拦截器视场角与目标区域的包含关系,再通过对队形参数进行调整,在中末制导时刻实现对目标的捕获,但该方法中并未给出截获概率,导致协同探测准确率低,因此,存在一定的局限性。
发明内容
本发明为了解决多拦截器探测时未给出截获概率,导致协同探测准确率低的问题,进而提出了一种基于截获概率的空空导弹协同探测方法。
本发明采取的技术方案是:
它包括以下步骤:
S1、获取考虑目标机动时,目标捕获时刻目标运动区域及其概率分布;
S2、选取一枚导弹作为领弹,四枚导弹作为从弹,根据目标运动区域计算领弹和从弹编队的队形参数;
S3、根据目标运动区域概率分布和队形参数计算二维平面内的协同截获概率;
S4、以队形参数为优化参数,截获概率为优化目标,利用粒子群优化算法优化协同截获概率,得到截获概率最大值。
进一步地,S1中获取考虑目标机动时,目标捕获时刻目标运动区域及其概率分布,具体过程为:
S11、根据惯性坐标系下的目标运动学方程和弹道坐标系下的目标动力学方程计算得到目标离散运动方程;
S12、在中制导段时,通过卫星测量得到目标在惯性坐标系下的位置,根据位置得到目标的控制输入;
S13、计算控制输入的均值和方差,根据均值和方差利用雅可比矩阵获得目标运动状态方差离散方程;
S14、根据目标离散运动方程和目标运动状态方差离散方程计算得到惯性坐标系下捕获时刻的目标运动区域,目标运动区域服从正态分布;
S15、根据目标运动区域计算惯性坐标系下的目标运动区域概率分布。
进一步地,S11中目标离散运动方程为:
X(k+1)=X(k)+ΔT×f(X(k),u(k),k) (1)
其中,X表示目标在惯性坐标系下的状态量,X=(x,y,z,V,θ,ψ)T,(x,y,z)表示目标在惯性坐标系下的位置坐标;V表示目标在惯性坐标系下的速度;θ表示目标在惯性坐标系下的弹道倾角;ψ表示目标在惯性坐标系下的弹道偏角;X(k+1)表示k+1时刻的目标信息;X(k)表示k时刻的目标信息;ΔT表示离散时间步长;f(·)表示非线性函数;u(k)表示k时刻目标的控制输入,u=(ax,ay,az)T,ax,ay,az表示目标在弹道坐标系下的x,y,z三个方向的加速度。
进一步地,S13中目标运动状态方差离散方程为:
Bk表示非线性函数f(·)对控制量uT的偏导数在X=X(k)时刻的雅可比矩阵;
进一步地,S14中目标运动区域为:
进一步地,S15中目标运动区域概率分布为:
进一步地,S2中选取一枚导弹作为领弹,四枚导弹作为从弹,根据目标运动区域计算领弹和从弹编队的队形参数,具体过程为:
S21、选取一枚导弹作为领弹M1,将S14的目标运动区域通过转换矩阵转换为领弹弹道坐标系下的目标运动区域,得到对应的均值和方差,且领弹弹道坐标系的x轴始终指向目标运动区域的均值,根据均值和方差构建目标三维运动区域,根据目标三维运动区域对目标可能出现区域的正态分布概率取±3σ,得到目标可能出现的近似估计区域;
S22、假设导弹的导引头为圆形视场,将目标三维运动区域的三维平面图分别沿着领弹弹道坐标系的不同坐标轴进行投影,得到三张二维平面投影示意图,分别是沿领弹弹道坐标系的x轴投影得到的二维平面投影示意图、沿领弹弹道坐标系的y轴投影得到的二维平面投影示意图、沿领弹弹道坐标系的z轴投影得到的二维平面投影示意图;
S23、选取四枚导弹作为从弹,从弹分别为M2、M3、M4、M5,从弹关于领弹M1空间对称分布,并与领弹M1保持人字形编队,并将领弹与四枚从弹进行S22中的投影,将投影结果标记在S22中三张二维平面投影示意图内,得到包含近似估计区域和领弹与从弹位置的三张二维平面投影示意图,设置领弹和从弹编队的队形参数r表示领弹和从弹之间的相对距离,ξ表示领弹和从弹的相对距离矢量与领弹弹道坐标系平面之间的夹角;/>表示领弹和从弹的相对距离矢量在领弹弹道坐标系平面的投影与M1Y轴之间的夹角;
S24、在包含近似估计区域和领弹与从弹位置的沿领弹弹道坐标系的x轴投影得到的二维平面投影示意图上计算领弹与从弹的探测区域,利用约束限制相邻导弹的探测区域间存在交集,并使领弹与四枚从弹的交集探测区域小于领弹的探测区域;
S25、根据探测区域,在包含近似估计区域和领弹与从弹位置的沿领弹弹道坐标系的x轴投影得到的二维平面投影示意图上建立极坐标系,定义(ρ′0,θ′0,r′0)作为极坐标系下的变量,ρ′0,θ′0分别表示从弹M3在领弹弹道坐标系空间的位置在领弹弹道坐标系M1YZ平面的投影时在极坐标下的极轴和极角;r′0表示领弹空间圆锥形探测区域在M1YZ平面投影的半径,通过(ρ′0,θ′0,r′0)计算得到队形参数
进一步地,S3中根据目标运动区域概率分布和队形参数计算二维平面内的协同截获概率,具体过程为:
以S23得到的包含近似估计区域和领弹与从弹位置的沿领弹弹道坐标系的x轴投影得到的二维平面投影示意图为基础,将所述示意图的目标运动区域概率分布作为积分函数,将S24计算的存在交集的探测区域作为积分区间,利用双重积分计算得到所述示意图的截获概率:
其中,P表示二维平面投影示意图的截获概率;Pq表示二维平面投影示意图某个象限的截获概率;σy表示目标捕获时刻运动状态方差转化至领弹弹道系中的y轴方向的方差;σz表示目标捕获时刻运动状态方差转化至领弹弹道系中的z轴方向的方差;d表示从弹开启导引头时与投影平面的距离;ρ′为极坐标的极轴;θ′为极坐标的极角;某个象限中象限圆的圆心坐标为(ρ′0,θ′0);RM2表示象限圆的半径。
进一步地,S3中根据目标运动区域概率分布和队形参数计算二维平面内的协同截获概率时,对领弹和从弹协同探测过程有如下要求:
1)、从弹的速度方向与领弹的速度方向保持一致;
2)、采用空间类人字形编队飞行,增大协同探测区域;
3)、从弹与领弹之间相互通信;
4)、从弹空间位置关于领弹对称。
进一步地,S4中以队形参数为优化参数,截获概率为优化目标,利用粒子群优化算法优化协同截获概率,得到截获概率最大值,具体过程为:
以导引头视场角为变化量,以队形参数为优化参数,截获概率为优化目标,利用粒子群优化算法获得不同导引头视场角下的队形参数及对应的截获概率,当导引头视场角越大时,协同截获概率越大,取协同截获概率最大值,此时空空导弹弹的协同探测效果最佳。
有益效果:
本发明在考虑目标机动时,假设目标初始运动状态以及目标机动均满足高斯分布,根据惯性坐标系下的目标运动学方程和弹道坐标系下的目标动力学方程,通过对目标运动方程离散化得到了目标捕获时刻的目标运动区域及其概率分布;从空空导弹中选取一枚导弹作为领弹,四枚导弹作为从弹,多弹以空间人字形构型,根据目标运动区域计算领弹和从弹编队的队形参数;以目标运动区域概率密度分布为积分函数,以多弹探测区域交集为积分区间,根据目标运动区域概率分布和队形参数计算二维平面内的协同截获概率;以队形参数为优化参数,截获概率为优化目标,利用粒子群优化算法优化得到中末制导交接班时刻协同截获概率。本发明提出基于截获概率的空空导弹协同探测方法,利用截获概率计算中末制导交接班时刻协同截获概率,提高了协同探测的准确率。
附图说明
图1是本发明的流程图;
图2是沿领弹弹道坐标系z轴投影时探测区域与目标运动区域的示意图;
图3是沿领弹弹道坐标系y轴投影探测区域与目标运动区域的示意图;
图4是沿领弹弹道坐标系x轴投影探测区域与目标运动区域的示意图;
图6是极坐标系下未满足约束的示意图;
图7是目标运动近似估计区域的二维平面示意图;
图8是极坐标系下协同截获概率求解示意图;
图9是视场角对截获概率的影响曲线图;
具体实施方式
具体实施方式一:结合图1-图9说明本实施方式,本实施方式所述一种基于截获概率的空空导弹协同探测方法,它包括以下步骤:
S1、获取考虑目标机动时,目标捕获时刻目标运动区域及其概率分布,具体过程为:
本发明通过考虑协同探测的目标运动信息以及中制导卫星测量误差,利用目标运动学方程和动力学方程得到协同探测的目标(简称目标)捕获时刻运动区域的概率分布。
惯性坐标系下的目标运动学方程:
其中,表示位置坐标(x,y,z)三个方向的导数,(x,y,z)表示目标在惯性坐标系下的位置坐标;V表示目标在惯性坐标系下的速度;θ表示目标在惯性坐标系下的弹道倾角;ψ表示目标在惯性坐标系下的弹道偏角。
弹道坐标系下的目标动力学方程:
其中,m表示目标的质量;表示速度的导数;P表示发动机的推力;α表示目标的攻角;β表示目标的侧滑角;G表示目标的阻力;/>表示弹道倾角的导数;γ表示目标的速度倾斜角;L表示目标的升力;g表示重力加速度;/>表示弹道偏角的导数;Z表示目标的侧向力。/>
基于过载的定义,描述目标动力学方程可用过载矢量在弹道坐标系各轴上的分量nx,ny,nz表示为:
过载为作用在目标上的除了重力以外所有外力的合力与目标重力的比值。由于加速度等于过载和重力加速度的乘积,即目标在弹道坐标系下的x,y,z三个方向的加速度为ax=nxg,ay=nyg,az=nzg,从而能够得到目标非线性运动学方程:
目标非线性动力学方程:
定义状态量X=(x,y,z,V,θ,ψ)T,u=(ax,ay,az)T,则上述方程(4)和(5)简化为
由于目标运动学方程和目标动力学方程是非线性的,因此,采用数值运算方法对公式(6)离散化,定义离散时间步长为ΔT,则目标离散运动方程为:
X(k+1)=X(k)+ΔT×f(X(k),u(k),k) (7)
其中,X(k+1)表示k+1时刻的目标信息;X(k)表示k时刻的目标信息;u(k)表示k时刻目标的控制输入。
在中制导段时,卫星测量到的直接信息是目标在惯性坐标系下的位置坐标(x,y,z),依次经过差分和滤波运算得到目标在惯性坐标系下的速度V和加速度a估计信息,根据目标在惯性系下的速度V信息求出目标的弹道倾角和弹道偏角ψ=arctan(vy/vx),vx,vy,vz分别为目标的速度V在惯性坐标系的三个分量。根据计算出的弹道倾角θ和弹道偏角ψ将惯性坐标系下估计出的目标的加速度信息a利用公式(8)转至弹道坐标系下,得到公式(9),进而由公式(9)得到控制输入u=(ax,ay,az)T。
惯性坐标系和弹道坐标系的转换矩阵:
其中,aIx,aIy,aIz表示目标在惯性坐标系下的x,y,z三个方向的加速度。
综上,根据卫星在k时刻测得的目标位置信息计算出k时刻目标的其他信息,以及控制输入u的均值和方差。假设控制输入u的均值和方差分别为uT,以及初始时刻控制输入u的均值和方差为uT0,/>则利用雅可比矩阵能够获得惯性坐标系下的目标运动状态方差离散方程:
其中,QX(k+1)表示k+1时刻目标运动状态的方差;QX(k)表示k时刻目标运动状态的方差;Ak表示非线性函数f(·)对X的偏导数在X=X(k)时刻的雅可比矩阵;
Bk表示非线性函数f(·)对控制量uT的偏导数在X=X(k)时刻的雅可比矩阵;
假设k时刻惯性坐标系下的目标运动状态的均值和方差分别为X(k),QX(k),通过公式(7)和公式(10)迭代计算出捕获时刻tb目标运动状态的均值和方差/>即考虑目标机动预测出的捕获时刻的目标运动分布区域/>为:
其中,N(·)表示正态分布函数;
由考虑目标机动得到的目标运动区域服从正态分布,得到目标运动区域中X概率密度为:
S2、选取一枚导弹作为领弹,四枚导弹作为从弹,根据目标运动区域计算领弹和从弹编队的队形参数,具体过程为:
同构或异构空空导弹协同探测时,需要确定导弹数量、弹间距以及弹间相对信息,从而在每枚导弹的探测视场角以及最大探测距离确定的情况下使得截获概率最大。选取一枚导弹作为领弹M1,将上述计算得到的惯性坐标系下的目标运动区域(公式(13))通过领弹的弹道坐标系和惯性坐标系的转换矩阵(公式(15))转换至领弹的弹道坐标系下,且领弹弹道坐标系的x轴始终指向目标预测区域的均值。
转换矩阵:
其中,θv为领弹的弹道倾角,ψv为领弹的弹道偏角。
经过转换得到:
若导弹存在姿态测量误差,由惯性坐标系转换至领弹弹道坐标系的转换矩阵就会存在误差,则公式(15)变为:
其中,Δθv为由姿态测量误差产生的领弹弹道倾角偏量,Δψv为由姿态测量误差产生的领弹弹道偏量。根据公式(18)更新存在姿态测量误差时领弹弹道坐标系下目标运动区域的均值和方差。
本发明采用捕获时刻考虑目标机动特性产生的目标在领弹弹道坐标系下运动状态的均值和方差(公式16、17)构建目标三维运动区域。为了简化问题的求解,根据目标三维运动区域对目标可能出现区域的正态分布概率取±3σ,得到目标可能出现的近似估计区域。
假设导弹的导引头为圆形视场,则将三维截获概率求解问题转化为二维平面截获概率求解问题,即将三维平面图分别沿着领弹弹道坐标系的不同坐标轴进行投影,得到三张二维平面投影示意图,分别是沿领弹弹道坐标系的x轴投影得到的二维平面投影示意图、沿领弹弹道坐标系的y轴投影得到的二维平面投影示意图、沿领弹弹道坐标系的z轴投影得到的二维平面投影示意图,则沿领弹弹道坐标系的x轴投影得到的二维平面投影示意图中包括近似估计区域。
考虑到目标运动区域关于均值点对称,为了保证对目标的截获,选择四枚导弹作为从弹(M2、M3、M4、M5)关于领弹空间对称,并组成人字形编队飞行。为了确定队形参数,将领弹与四枚从弹进行上述投影,将投影结果标记在每张二维平面投影示意图内,得到包含近似估计区域和领弹与从弹位置的三张二维平面投影示意图,如图2-图4所示。在领弹的弹道坐标系下重新设定领弹和从弹编队的三个队形参数如图5所示,r表示领弹和从弹之间的相对距离;ξ表示领弹和从弹在领弹弹道坐标系下的相对距离矢量与领弹弹道坐标系的M1YZ平面之间的夹角;/>表示领弹和从弹在领弹弹道坐标系下的相对距离矢量在领弹弹道坐标系的M1YZ平面的投影与M1Y轴之间的夹角。在队形保持控制中将上述三个参数转为相应的坐标,从而实现队形控制律的设计。
为了计算包含近似估计区域和领弹与从弹位置的沿领弹弹道坐标系的x轴投影得到的二维平面投影示意图(图4)中领弹与四枚从弹在投影平面上探测区域圆的半径大小,假设每个从弹开启导引头时距离投影平面为d=dmax-6σx,其中,dmax为从弹导引头最大截获距离,σx为目标捕获时刻运动状态方差转化至领弹弹道系中x轴方向的方差,θm为从弹导引头视线角,则领弹在投影平面投影的圆半径为RM1=(d-rsinξ)tanθm,从弹在投影平面投影的圆半径为RM2=dtanθm。此外,为了保证相邻导弹的视场投影区域有交集防止出现漏扫,且中间圆没有被完全覆盖,防止交叠区域过大,需要对相邻导弹的中心距进行约束,如图6所示。
根据探测区域,在图4上建立极坐标系,采用极坐标对协同探测截获概率进行求解,定义变量(ρ′0,θ′0,r′0)作为极坐标系下的求解变量,ρ′0,θ′0分别为从弹M3在领弹弹道坐标系空间的位置在领弹弹道坐标系M1YZ平面的投影时在极坐标下的极轴和极角;r′0为领弹空间圆锥形探测区域在M1YZ平面投影的半径大小。再通过变量(ρ′0,θ′0,r′0)反解出导弹的队形参数得到/>
S3、计算二维平面内的协同截获概率,具体过程为:
根据图4求解截获概率时,对多空空导弹协同探测过程有如下要求:
1)、从弹的速度方向与领弹的速度方向保持一致;
2)、采用空间类人字形编队飞行,增大协同探测区域;
3)、从弹与领弹之间相互通信;
4)、从弹空间位置关于领弹对称。
根据图4能够得到图7,如图7所示,以图4为基础,图4中的M1yz平面的概率分布(公式(14))作为积分函数,领弹和从弹投影区域的并集作为积分区间,从而通过双重积分求出二维平面下的截获概率,具体过程为:
由于图8中四个象限的区域面积和概率分布函数均相等,因此只选择第一象限进行截获概率的求解。在极坐标下定义第一象限圆的圆心坐标为(ρ′0,θ′0),半径为RM2,则第一象限圆的极坐标方程表示为:
θ′为极坐标的极角;
极坐标下目标运动状态概率分布密度表示为:
σy表示目标捕获时刻运动状态方差转化至领弹弹道系中的y轴方向的方差;σz表示目标捕获时刻运动状态方差转化至领弹弹道系中的z轴方向的方差
则截获概率计算公式为:
P=4P1 (23)
则存在约束(下述的等式中的s.t中的表达式)时的协同截获概率优化问题的数学描述如下:
其中,P表示二维平面投影示意图的截获概率;P1表示二维平面投影示意图某个象限的截获概率。
S4、以队形参数为优化参数,截获概率为优化目标,利用粒子群优化算法优化得到中末制导交接班时刻协同截获概率最大值。
上述约束条件中,对圆的半径r′0约束的下限是保证在象限中没有漏扫区域;约束上限是保证中间圆没有被完全覆盖,防止交叠区域过大。
经过上述的优化,通过粒子群优化算法求解协同截获概率及队形参数的值。导引头视场角θm、目标捕获时刻运动状态方差转化至领弹弹道系中的x轴方向方差σx、y轴方向方差σy、z轴方向方差σz以及导引头最大探测距离dmax的变化均会对最大截获概率产生影响,这里以导引头视场角为变化量,以队形参数为优化参数,截获概率为优化目标,利用粒子群优化算法获得不同导引头视场角下的多弹编队队形参数及对应的截获概率,求解结果如表1和图9所示:
表1不同导引头视场角下的仿真结果
从图9中能够看出,导引头视场角越大,截获概率越大,增大导引头的视场角可以显著提高截获概率,从表1可以看出,在上述参数设置下,当导引头视场角取为6°时,截获概率接近于1,因此导引头视场角的大小对截获概率有很大的影响,所以当协同截获概率最大时,多个空空导弹的协同探测效果最佳。
Claims (8)
1.基于截获概率的空空导弹协同探测方法,其特征在于:它包括以下步骤:
S1、获取考虑目标机动时,目标捕获时刻的目标运动区域及其概率分布;
S2、选取一枚导弹作为领弹,四枚导弹作为从弹,根据目标运动区域计算领弹和从弹编队的队形参数,具体过程为:
S21、选取一枚导弹作为领弹M1,将S1的目标运动区域通过转换矩阵转换为领弹弹道坐标系下的目标运动区域,得到对应的均值和方差,且领弹弹道坐标系的x轴始终指向目标运动区域的均值,根据均值和方差构建目标三维运动区域,根据目标三维运动区域对目标可能出现区域的正态分布概率取±3σ,得到目标可能出现的近似估计区域;
S22、假设导弹的导引头为圆形视场,将目标三维运动区域的三维平面图分别沿着领弹弹道坐标系的不同坐标轴进行投影,得到三张二维平面投影示意图,分别是沿领弹弹道坐标系的x轴投影得到的二维平面投影示意图、沿领弹弹道坐标系的y轴投影得到的二维平面投影示意图、沿领弹弹道坐标系的z轴投影得到的二维平面投影示意图;
S23、选取四枚导弹作为从弹,从弹分别为M2、M3、M4、M5,从弹关于领弹M1空间对称分布,并与领弹M1保持人字形编队,并将领弹与四枚从弹进行S22中的投影,将投影结果标记在S22中三张二维平面投影示意图内,得到包含近似估计区域和领弹与从弹位置的三张二维平面投影示意图,设置领弹和从弹编队的队形参数r表示领弹和从弹之间的相对距离,ξ表示领弹和从弹的相对距离矢量与领弹弹道坐标系平面之间的夹角;/>表示领弹和从弹的相对距离矢量在领弹弹道坐标系平面的投影与M1Y轴之间的夹角;
S24、在包含近似估计区域和领弹与从弹位置的沿领弹弹道坐标系的x轴投影得到的二维平面投影示意图上计算领弹与从弹的探测区域,利用约束限制相邻导弹的探测区域间存在交集,并使领弹与四枚从弹的交集探测区域小于领弹的探测区域;
S25、根据探测区域,在包含近似估计区域和领弹与从弹位置的沿领弹弹道坐标系的x轴投影得到的二维平面投影示意图上建立极坐标系,定义(ρ′0,θ′0,r′0)作为极坐标系下的变量,ρ′0,θ′0分别表示从弹M3在领弹弹道坐标系空间的位置在领弹弹道坐标系M1YZ平面的投影时在极坐标下的极轴和极角;r′0表示领弹空间圆锥形探测区域在M1YZ平面投影的半径,通过(ρ′0,θ′0,r′0)计算得到队形参数
S3、根据目标运动区域概率分布和队形参数计算二维平面内的协同截获概率,具体过程为:
以S23得到的包含近似估计区域和领弹与从弹位置的沿领弹弹道坐标系的x轴投影得到的二维平面投影示意图为基础,将所述示意图的目标运动区域概率分布作为积分函数,将S24计算的存在交集的探测区域作为积分区间,利用双重积分计算得到所述示意图的截获概率:
s.t 0≤θ′0≤π/2
其中,P表示二维平面投影示意图的截获概率;Pq表示二维平面投影示意图某个象限的截获概率;σy表示目标捕获时刻运动状态方差转化至领弹弹道系中的y轴方向的方差;σz表示目标捕获时刻运动状态方差转化至领弹弹道系中的z轴方向的方差;d表示从弹开启导引头时与投影平面的距离;ρ′为极坐标的极轴;θ′为极坐标的极角;某个象限中象限圆的圆心坐标为(ρ′0,θ′0);RM2表示象限圆的半径;
S4、以队形参数为优化参数,截获概率为优化目标,利用粒子群优化算法优化协同截获概率,得到截获概率最大值。
2.根据权利要求1中所述的基于截获概率的空空导弹协同探测方法,其特征在于:S1中获取考虑目标机动时,目标捕获时刻的目标运动区域及其概率分布,具体过程为:
S11、根据惯性坐标系下的目标运动学方程和弹道坐标系下的目标动力学方程计算得到目标离散运动方程;
S12、在中制导段时,通过卫星测量得到目标在惯性坐标系下的位置,根据位置得到目标的控制输入;
S13、计算控制输入的均值和方差,根据均值和方差利用雅可比矩阵获得目标运动状态方差离散方程;
S14、根据目标离散运动方程和目标运动状态方差离散方程计算得到惯性坐标系下捕获时刻的目标运动区域,目标运动区域服从正态分布;
S15、根据目标运动区域计算惯性坐标系下的目标运动区域概率分布。
3.根据权利要求2中所述的基于截获概率的空空导弹协同探测方法,其特征在于:S11中目标离散运动方程为:
X(k+1)=X(k)+ΔT×f(X(k),u(k),k) (1)
其中,X表示目标在惯性坐标系下的状态量,X=(x,y,z,V,θ,ψ)T,(x,y,z)表示目标在惯性坐标系下的位置坐标;V表示目标在惯性坐标系下的速度;θ表示目标在惯性坐标系下的弹道倾角;ψ表示目标在惯性坐标系下的弹道偏角;X(k+1)表示k+1时刻的目标信息;X(k)表示k时刻的目标信息;ΔT表示离散时间步长;f(·)表示非线性函数;u(k)表示k时刻目标的控制输入,u=(ax,ay,az)T,ax,ay,az表示目标在弹道坐标系下的x,y,z三个方向的加速度。
7.根据权利要求6中所述的基于截获概率的空空导弹协同探测方法,其特征在于:S3中根据目标运动区域概率分布和队形参数计算二维平面内的协同截获概率时,对领弹和从弹协同探测过程有如下要求:
1)、从弹的速度方向与领弹的速度方向保持一致;
2)、采用空间类人字形编队飞行,增大协同探测区域;
3)、从弹与领弹之间相互通信;
4)、从弹空间位置关于领弹对称。
8.根据权利要求7中所述的基于截获概率的空空导弹协同探测方法,其特征在于:S4中以队形参数为优化参数,截获概率为优化目标,利用粒子群优化算法优化协同截获概率,得到截获概率最大值,具体过程为:
以导引头视场角为变化量,以队形参数为优化参数,截获概率为优化目标,利用粒子群优化算法获得不同导引头视场角下的队形参数及对应的截获概率,当导引头视场角越大时,协同截获概率越大,取协同截获概率最大值,此时空空导弹弹的协同探测效果最佳。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211362631.9A CN115755955B (zh) | 2022-11-02 | 2022-11-02 | 基于截获概率的空空导弹协同探测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211362631.9A CN115755955B (zh) | 2022-11-02 | 2022-11-02 | 基于截获概率的空空导弹协同探测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115755955A CN115755955A (zh) | 2023-03-07 |
CN115755955B true CN115755955B (zh) | 2023-06-02 |
Family
ID=85355436
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211362631.9A Active CN115755955B (zh) | 2022-11-02 | 2022-11-02 | 基于截获概率的空空导弹协同探测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115755955B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101286071A (zh) * | 2008-04-24 | 2008-10-15 | 北京航空航天大学 | 基于微粒群优化和遗传算法的多无人机三维编队重构方法 |
WO2012048287A2 (en) * | 2010-10-08 | 2012-04-12 | Loctronix Corporation | Doppler aided inertial navigation |
CN112558634A (zh) * | 2020-12-08 | 2021-03-26 | 北京航空航天大学 | 一种多导弹协同制导方法及系统 |
CN114371737A (zh) * | 2022-01-17 | 2022-04-19 | 哈尔滨工业大学 | 一种角度约束的智能弹药最优对策协同制导方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111221352B (zh) * | 2020-03-03 | 2021-01-29 | 中国科学院自动化研究所 | 基于多无人机协同博弈对抗的控制系统 |
-
2022
- 2022-11-02 CN CN202211362631.9A patent/CN115755955B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101286071A (zh) * | 2008-04-24 | 2008-10-15 | 北京航空航天大学 | 基于微粒群优化和遗传算法的多无人机三维编队重构方法 |
WO2012048287A2 (en) * | 2010-10-08 | 2012-04-12 | Loctronix Corporation | Doppler aided inertial navigation |
CN112558634A (zh) * | 2020-12-08 | 2021-03-26 | 北京航空航天大学 | 一种多导弹协同制导方法及系统 |
CN114371737A (zh) * | 2022-01-17 | 2022-04-19 | 哈尔滨工业大学 | 一种角度约束的智能弹药最优对策协同制导方法 |
Non-Patent Citations (9)
Title |
---|
Three-dimensional consensus algorithm for nonsingular distributed cooperative guidance strategy;xiaoyan yang等;《Aerospace Science and Technology》;第1-18页 * |
中远程制导火箭弹协同攻击动力学分析与控制;杨小艳;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》(第01(2021)期);第C032-56页 * |
众多导弹协同作战制导控制的研究;吴胜亮;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》(第06(2014)期);第C032-58页 * |
分布式MIMO雷达目标检测与跟踪及实验研究;余万玉;《中国优秀硕士学位论文全文数据库信息科技辑》(第07(2020)期);第I136-728页 * |
协同制导导弹截获概率仿真分析;王向坡等;《第八届中国指挥控制大会论文集》;第622-627页 * |
多拦截器协同制导控制问题研究;王家炜;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》(第02(2015)期);第C032-62页 * |
多阶段高斯伪谱法在编队最优控制中的应用;王芳;林涛;张克;崔乃刚;;宇航学报;第36卷(第11期);第1262-1269页 * |
导弹协同制导中末交接班目标截获概率研究;付昭旺;于雷;寇英信;徐安;;弹道学报;第24卷(第04期);第31-36页 * |
编队协同空空攻击导弹截获概率仿真分析;张明明等;《航空兵器》;第28卷(第01期);第72-76页 * |
Also Published As
Publication number | Publication date |
---|---|
CN115755955A (zh) | 2023-03-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yanushevsky | Modern missile guidance | |
CN109597423B (zh) | 一种基于参考视线角信号的多约束末制导律的设计方法 | |
CN110412874B (zh) | 针对机动目标和时延通信的多导弹协同制导律设计方法 | |
CN108534614B (zh) | 一种三维全向实时预测制导方法 | |
CN111649624B (zh) | 一种空间微型精确制导武器控制方法 | |
CN110345814B (zh) | 一种不依赖自身导引头测量信息的末制导算法 | |
CN112033234B (zh) | 一种多导弹协同制导律设计方法 | |
CN108073742B (zh) | 基于改进粒子滤波算法的拦截导弹末段飞行状态估计方法 | |
CN112648886A (zh) | 一种组合制导目标拦截方法及系统 | |
CN115857538A (zh) | 三维空间下满足落角约束的多飞行器协同制导方法 | |
CN114020021A (zh) | 一种多导弹分布式协同制导律的设计方法及系统 | |
CN114035616A (zh) | 一种飞行器对移动目标打击控制方法及系统 | |
CN115755955B (zh) | 基于截获概率的空空导弹协同探测方法 | |
CN109780933B (zh) | 一种单兵制导火箭动态目标预测导引方法 | |
CN110471283B (zh) | 一种带碰撞角约束的三维鲁棒制导律构建方法 | |
CN111930142A (zh) | 速度不可控条件下的多弹编队协同控制方法 | |
CN116795130A (zh) | 一种针对机动目标的避障制导方法 | |
CN114153143A (zh) | 一种导弹非奇异固定时间滑模制导律的设计方法 | |
Zuoe et al. | Study on vertical attack penetration probability of anti-ship missile | |
CN112379599A (zh) | 一种针对空中机动目标的三维协同中制导方法 | |
Sultan et al. | The Reachable Set of a Drone: Exploring the Position Isochrones for a Quadcopter | |
Wang et al. | Optimization of tactical aircraft weapon delivery using Tactics Templates | |
Tang et al. | Research on Target State Estimation and Terminal Guidance Algorithm in the Process of Multi-UAV Cooperative Attack | |
Li et al. | Anti-jamming Trajectory Planning of Infrared Imaging Air-to-air Missile | |
CN117891271B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |