CN106643298A - 一种基于预置命中点的大气层内反导拦截器中制导方法 - Google Patents
一种基于预置命中点的大气层内反导拦截器中制导方法 Download PDFInfo
- Publication number
- CN106643298A CN106643298A CN201611074394.0A CN201611074394A CN106643298A CN 106643298 A CN106643298 A CN 106643298A CN 201611074394 A CN201611074394 A CN 201611074394A CN 106643298 A CN106643298 A CN 106643298A
- Authority
- CN
- China
- Prior art keywords
- missile
- impact
- pip
- interception device
- preset point
- 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.)
- Granted
Links
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F41—WEAPONS
- F41G—WEAPON SIGHTS; AIMING
- F41G5/00—Elevating or traversing control systems for guns
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F41—WEAPONS
- F41H—ARMOUR; ARMOURED TURRETS; ARMOURED OR ARMED VEHICLES; MEANS OF ATTACK OR DEFENCE, e.g. CAMOUFLAGE, IN GENERAL
- F41H11/00—Defence installations; Defence devices
- F41H11/02—Anti-aircraft or anti-guided missile or anti-torpedo defence installations or systems
Landscapes
- Engineering & Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)
Abstract
一种基于预置命中点的大气层内反导拦截器中制导方法,射前根据反导拦截器拦截能力与性能指标确定拦截范围,对进攻弹弹道进行实时预报,获得初始预置命中点和发射俯仰角,进行射中实时解算,得到每一时刻的预置命中点,直到中制导结束。本发明方法适用于大气层内反导拦截器中制导,计算效率高,能够在有限的能量下提高拦截成功率。
Description
技术领域
本发明涉及一种基于预置命中点的大气层内反导拦截器中制导方法,属于精确制导与控制领域。
背景技术
基于预置命中点的反导拦截器中制导算法是反导拦截制导与指控领域的热点,为提高命中精度,当前的主要解决思路是:在中段制导的过程中确定来袭目标经过的路径点之一,引导拦截器在与目标相同的时间到达路径点,通过将跟踪动目标转换为跟踪静目标,可减少拦截器能量消耗与过载需求。
上述方法计算量较低,但仍有以下不足之处:
1、仅适用于大气层外对惯性弹道进行拦截的中段反导系统,如“标准”-3型反导拦截器,不适用于大气层内反导拦截器中制导。
2、使用黄金比例分割法进行迭代收敛,虽然计算量较小,但收敛速度较慢,计算效率低,未充分利用现代计算机的性能以提高收敛速度。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提供一种基于预置命中点的大气层内反导拦截器中制导方法,提高了计算效率,最大限度地发挥拦截器的性能。
本发明的技术解决方案是:一种基于预置命中点的大气层内反导拦截器中制导方法,包括以下步骤:
步骤一:对反导拦截器进行射前预处理,具体步骤如下:
(1.1)根据反导拦截器拦截能力与性能指标获得拦截高度的上界Hmax、下界Hmin及最大拦截射程Rmax;
(1.2)对进攻弹弹道进行实时预报;
(1.3)在预报得到的弹道上选取处于反导拦截器拦截范围边界上的点作为初始预置命中点PIP0,根据初始预置命中点PIP0和反导拦截器部署点确定发射俯仰角,并计算反导拦截器飞行至初始预置命中点PIP0所需时间MISSILE_tgo;
(1.4)利用公式TARGET_tgo-MISSILE_tgo=tremain计算反导拦截器的待发射时间tremain,若tremain≤0,则反导拦截器立刻发射且转入步骤二,将初始预置命中点PIP0赋值给预置命中点PIP;若tremain>0,则进入(1.5),其中TARGET_tgo为进攻弹到达初始预置命中点PIP0的时间;
(1.5)重复(1.2)-(1.4),更新初始预置命中点PIP0,直到tremain≤0;
步骤二:进行射中实时解算,得到每一时刻的预置命中点,具体步骤如下:
(2.1)对进攻弹弹道进行实时预报,获得t时刻进攻弹到达预置命中点PIP的时间进入(2.2);
(2.2)计算t时刻反导拦截器飞行至预置命中点PIP所需时间进入(2.3);
(2.3)当时,为当前时刻反导拦截器飞行至PIP所需时间,预置命中点不变,t的值加1,判断中制导是否结束,如果结束,则弹中实时解算结束,如果中制导未结束,则返回(2.1);如果Δtgo>σ,则进入(2.4);其中σ为发射阈值;
(2.4)如果则令否则K为迭代速率系数,根据和进攻弹当前速度确定t时刻新的预置命中点PIP,返回(2.2)。
所述步骤(1.2)和(2.1)中,利用如下公式预报t时刻的进攻弹弹道:
其中(x,y,z)为t时刻雷达测量得到的进攻弹位置坐标,(vx0,vy0,vz0)为t时刻雷达测量得到的进攻弹速度,(ax,ay,az)为t时刻雷达解算得到的进攻弹加速度。
所述步骤(1.3)中,计算反导拦截器飞行至PIP0所需时间MISSILE_tgo的方法如下:
(3.1)通过多次打靶试验,获得包含反导拦截器飞行时间、飞行距离、风场数据以及发射俯仰角的射表,利用基于高斯核函数的最小二乘法对射表数据进行拟合,得到飞行平均速度相对于发射俯仰角及风场数据的函数
(3.2)利用发射俯仰角及当地风场数据获得飞行平均速度;
(3.3)利用公式确定反导拦截器飞行至PIP0所需时间MISSILE_tgo,其中S0为反导拦截器部署点到初始预置命中点PIP0的距离,为反导拦截器飞行平均速度。
所述步骤(2.3)中,迭代速率系数K的确定方法如下:
(4.1)设计关于K的模糊查询表,确定模糊规则,根据模糊规则确定模糊查询表中各个元素的值;
模糊查询表
所述模糊规则为:横轴档位和纵轴档位均与K的值正相关;
(4.2)设定每一档对应的Δtgo范围和Δtgo变化速率绝对值的范围,模糊查询表中横轴档位对应Δtgo,纵轴档位对应Δtgo的变化速率绝对值;
(4.3)按照当前时刻计算得到的Δtgo和Δtgo变化速率绝对值,从模糊查询表中选取K值。
与现有技术相比,本发明具有如下有益效果:
(1)由于大气层内拦截器拦截范围有限,传统在中段制导的过程中确定来袭目标经过的路径点的方法可能导致能量提前耗尽,本发明方法通过射前预处理确定反导拦截器的发射时刻和发射俯仰角,适用于大气层内反导拦截器中制导,能够在有限的能量下提高拦截成功率。
(2)利用模糊优化算法计算进攻弹到达预置命中点的时间,相较于黄金比例分割法,提高了计算效率和收敛速度。
(3)本发明方法利用基于高斯核函数的最小二乘法对射表数据进行拟合,获得飞行平均速度相对于发射俯仰角及风场数据的函数,相较于传统插值方法,计算效率高,同时考虑了风场干扰,求解更加精确。
附图说明
图1为末段拦截初始预置命中点设置示意图;
图2为射前预处理流程图;
图3为射中实时解算流程图。
具体实施方式
当前反导拦截制导算法的应用领域为大气层外中段拦截的惯性弹道的远程导弹,且由于使用黄金比例分割法改变进攻弹与拦截器的剩余飞行时间,计算量虽然较小,但迭代速度较慢,同时无法根据拦截状态进行自适应改变迭代步长。基于以上问题,针对大气层内反导拦截器中制导需求,本发明提出了一种基于预置命中点的大气层内反导拦截器中制导方法,包括射前预处理和射中实时解算两部分,在射中实时解算时使用模糊优化算法改变迭代步长,增加收敛速度。
本发明的步骤如下:
步骤一:射前预处理:
(1.1)不同于对弹道导弹进行中段拦截的制导算法,末段拦截器拦截时间有限。为了最大限度增加末段拦截器飞行时间,第一个预置命中点理想状态应位于反导拦截器最大拦截范围的边界上。因此,本发明首先根据反导拦截器拦截能力与性能指标获得拦截高度的上界Hmax、下界Hmin及最大拦截射程Rmax,如图1所示。
(1.2)对进攻弹弹道进行实时预报。此中制导算法的拦截目标范围为中近程弹道导弹,在末段时基本处于无机动或短时间机动阶段,可认为仅受重力影响,因此其t时刻的弹道预报公式如下:
其中(x,y,z)为t时刻雷达测量得到的进攻弹位置坐标,(vx0,vy0,vz0)为t时刻雷达测量得到的进攻弹速度,(ax,ay,az)为t时刻雷达解算得到的进攻弹加速度。
(1.3)在预报得到的弹道上选取处于反导拦截器拦截范围边界上的点作为初始预置命中点PIP0,根据初始预置命中点PIP0和反导拦截器部署点确定发射俯仰角,并计算反导拦截器飞行至PIP0所需时间MISSILE_tgo。
不同于大气层外中段拦截器在中制导段发动机推力曲线与推力时间基本固定,从而可以计算出其平均速度Vaver;末段大气层内拦截器在初始发动机推力段过后,其不再具备轴向推力,速度随着飞行时间的增加而受空气阻力影响下降,但空气阻力的大小又与其飞行状态,即攻角、侧滑角等因素有关。为充分考虑其在各种飞行状态下的平均速度,使用最小二乘法进行拟合。
使用最小二乘法计算反导拦截器飞行至PIP0所需时间MISSILE_tgo的方法如下:
(a)通过多次打靶试验,获得包含反导拦截器飞行时间、飞行距离、风场数据以及发射俯仰角的射表,利用基于高斯核函数的最小二乘法对射表数据进行拟合,得到飞行平均速度相对于发射俯仰角及风场数据的函数
(b)利用发射俯仰角及当地风场数据获得飞行平均速度;
(c)利用公式确定反导拦截器飞行至PIP0所需时间MISSILE_tgo,其中S0为反导拦截器部署点到初始预置命中点PIP0的距离,为反导拦截器飞行平均速度。
使用最小二乘法进行拟合时,基函数可使用线性模型,也可以使用高斯核函数;通过拟合可得到平均误差最小状态下的函数,这样即增强了算法鲁棒性。
(1.4)利用公式TARGET_tgo-MISSILE_tgo=tremain计算反导拦截器的待发射时间tremain,若tremain≤0,则反导拦截器立刻发射且转入步骤二,将初始预置命中点PIP0赋值给预置命中点PIP;若tremain>0,则进入(1.5),其中TARGET_tgo为进攻弹到达初始预置命中点PIP0的时间;
(1.5)重复(1.2)-(1.4),更新初始预置命中点PIP0,直到tremain≤0tremain<0;
射前预处理如图2所示。
步骤二:进行射中实时解算,得到每一时刻的预置命中点,具体步骤如下:
(2.1)对进攻弹弹道进行实时预报,获得t时刻进攻弹到达预置命中点PIP的时间(即进攻弹剩余飞行时间),进入(2.2);
(2.2)计算t时刻反导拦截器飞行至预置命中点PIP所需时间计算法方法与(1.3)一样,只是S0改为反导拦截器部署点到预置命中点PIP的距离,进入(2.3);
(2.3)令当Δtgo≤σ时,为当前时刻反导拦截器飞行至PIP所需时间(即拦截器剩余飞行时间),预置命中点不变,t的值加1,判断中制导是否结束,如果结束,则射中实时解算结束,如果未结束,则进入(2.1);如果Δtgo>σ,否则进入(2.4);其中σ为正实数,取0.5s~1s;
(2.4)如果则令否则K为迭代速率系数,根据和进攻弹当前速度确定t时刻新的预置命中点PIP,返回(2.2)。
为了较好的平衡计算时间与迭代优化时间,引入模糊优化算法获取迭代速率系数K,模糊优化算法具体内容如下:
(a)设计关于K的模糊查询表,根据模糊规则确定的模糊查询表,模糊查询表的一种形式如下:
模糊查询表
(b)设定每一档对应的Δtgo和Δtgo变化速率绝对值的范围,模糊查询表中横轴档位对应Δtgo,纵轴档位对应Δtgo的变化速率绝对值。
一种设定方式为:横轴档位0对应Δtgo的范围为(0,1),横轴档位1对应Δtgo的范围为(1,10),横轴档位2对应Δtgo的范围为(10,30),横轴档位3对应Δtgo的范围为(30,200)。
纵轴档位0对应Δtgo的变化速率绝对值范围为(0,0.5),纵轴档位1对应Δtgo的变化速率绝对值范围为(0.5,1),纵轴档位2对应Δtgo的变化速率绝对值范围为(1,3),纵轴档位3对应Δtgo的变化速率绝对值范围为(3,10)。
(c)按照当前时刻计算得到的Δtgo和Δtgo变化速率绝对值,从模糊查询表中选取K值。
射中实时解算流程如图3所示。
上述模糊优化算法方法以模糊数学为基础,可保证在初期进行大幅收敛,加快收敛速度。
仿真试验表明,本发明方法良好地应用于大气层内拦截,提高了反导系统拦截概率,并提高了计算速度。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (4)
1.一种基于预置命中点的大气层内反导拦截器中制导方法,其特征在于包括以下步骤:
步骤一:对反导拦截器进行射前预处理,具体步骤如下:
(1.1)根据反导拦截器拦截能力与性能指标获得拦截高度的上界Hmax、下界Hmin及最大拦截射程Rmax;
(1.2)对进攻弹弹道进行实时预报;
(1.3)在预报得到的弹道上选取处于反导拦截器拦截范围边界上的点作为初始预置命中点PIP0,根据初始预置命中点PIP0和反导拦截器部署点确定发射俯仰角,并计算反导拦截器飞行至初始预置命中点PIP0所需时间MISSILE_tgo;
(1.4)利用公式TARGET_tgo-MISSILE_tgo=tremain计算反导拦截器的待发射时间tremain,若tremain≤0,则反导拦截器立刻发射且转入步骤二,将初始预置命中点PIP0赋值给预置命中点PIP;若tremain>0,则进入(1.5),其中TARGET_tgo为进攻弹到达初始预置命中点PIP0的时间;
(1.5)重复(1.2)-(1.4),更新初始预置命中点PIP0,直到tremain≤0;
步骤二:进行射中实时解算,得到每一时刻的预置命中点,具体步骤如下:
(2.1)对进攻弹弹道进行实时预报,获得t时刻进攻弹到达预置命中点PIP的时间进入(2.2);
(2.2)计算t时刻反导拦截器飞行至预置命中点PIP所需时间进入(2.3);
(2.3)当时,为当前时刻反导拦截器飞行至PIP所需时间,预置命中点不变,t的值加1,判断中制导是否结束,如果结束,则弹中实时解算结束,如果中制导未结束,则返回(2.1);如果Δtgo>σ,则进入(2.4);其中σ为发射阈值;
(2.4)如果则令否则K为迭代速率系数,根据和进攻弹当前速度确定t时刻新的预置命中点PIP,返回(2.2)。
2.根据权利要求1所述的一种基于预置命中点的大气层内反导拦截器中制导方法,其特征在于:所述步骤(1.2)和(2.1)中,利用如下公式预报t时刻的进攻弹弹道:
其中(x,y,z)为t时刻雷达测量得到的进攻弹位置坐标,(vx0,vy0,vz0)为t时刻雷达测量得到的进攻弹速度,(ax,ay,az)为t时刻雷达解算得到的进攻弹加速度。
3.根据权利要求1所述的一种基于预置命中点的大气层内反导拦截器中制导方法,其特征在于:所述步骤(1.3)中,计算反导拦截器飞行至PIP0所需时间MISSILE_tgo的方法如下:
(3.1)通过多次打靶试验,获得包含反导拦截器飞行时间、飞行距离、风场数据以及发射俯仰角的射表,利用基于高斯核函数的最小二乘法对射表数据进行拟合,得到飞行平均速度相对于发射俯仰角及风场数据的函数
(3.2)利用发射俯仰角及当地风场数据获得飞行平均速度;
(3.3)利用公式确定反导拦截器飞行至PIP0所需时间MISSILE_tgo,其中S0为反导拦截器部署点到初始预置命中点PIP0的距离,为反导拦截器飞行平均速度。
4.根据权利要求1所述的一种基于预置命中点的大气层内反导拦截器中制导方法,其特征在于:所述步骤(2.3)中,迭代速率系数K的确定方法如下:
(4.1)设计关于K的模糊查询表,确定模糊规则,根据模糊规则确定模糊查询表中各个元素的值;
模糊查询表
所述模糊规则为:横轴档位和纵轴档位均与K的值正相关;
(4.2)设定每一档对应的Δtgo范围和Δtgo变化速率绝对值的范围,模糊查询表中横轴档位对应Δtgo,纵轴档位对应Δtgo的变化速率绝对值;
(4.3)按照当前时刻计算得到的Δtgo和Δtgo变化速率绝对值,从模糊查询表中选取K值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611074394.0A CN106643298B (zh) | 2016-11-29 | 2016-11-29 | 一种基于预置命中点的大气层内反导拦截器中制导方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611074394.0A CN106643298B (zh) | 2016-11-29 | 2016-11-29 | 一种基于预置命中点的大气层内反导拦截器中制导方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106643298A true CN106643298A (zh) | 2017-05-10 |
CN106643298B CN106643298B (zh) | 2018-04-27 |
Family
ID=58813168
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611074394.0A Active CN106643298B (zh) | 2016-11-29 | 2016-11-29 | 一种基于预置命中点的大气层内反导拦截器中制导方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106643298B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109857145A (zh) * | 2018-11-27 | 2019-06-07 | 北京航空航天大学 | 一种基于迭代预测命中点的增程型拦截弹预测制导方法 |
CN114357807A (zh) * | 2022-03-11 | 2022-04-15 | 北京航空航天大学 | 一种用于大气层外拦截的最优中制导方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0655599B1 (fr) * | 1993-11-25 | 1998-07-08 | AEROSPATIALE Société Nationale Industrielle | Système de défense antiaérien et missile de défense pour un tel système |
US7446291B1 (en) * | 2005-10-03 | 2008-11-04 | Lockheed Martin Corporation | Augmented proportional navigation guidance law using angular acceleration measurements |
CN101614507A (zh) * | 2008-06-29 | 2009-12-30 | 范海涛 | 一种导弹反拦截方法 |
CN104266546A (zh) * | 2014-09-22 | 2015-01-07 | 哈尔滨工业大学 | 一种基于视线的有限时间收敛主动防御制导控制方法 |
-
2016
- 2016-11-29 CN CN201611074394.0A patent/CN106643298B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0655599B1 (fr) * | 1993-11-25 | 1998-07-08 | AEROSPATIALE Société Nationale Industrielle | Système de défense antiaérien et missile de défense pour un tel système |
US7446291B1 (en) * | 2005-10-03 | 2008-11-04 | Lockheed Martin Corporation | Augmented proportional navigation guidance law using angular acceleration measurements |
CN101614507A (zh) * | 2008-06-29 | 2009-12-30 | 范海涛 | 一种导弹反拦截方法 |
CN104266546A (zh) * | 2014-09-22 | 2015-01-07 | 哈尔滨工业大学 | 一种基于视线的有限时间收敛主动防御制导控制方法 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109857145A (zh) * | 2018-11-27 | 2019-06-07 | 北京航空航天大学 | 一种基于迭代预测命中点的增程型拦截弹预测制导方法 |
CN114357807A (zh) * | 2022-03-11 | 2022-04-15 | 北京航空航天大学 | 一种用于大气层外拦截的最优中制导方法及装置 |
CN114357807B (zh) * | 2022-03-11 | 2022-05-27 | 北京航空航天大学 | 一种用于大气层外拦截的最优中制导方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN106643298B (zh) | 2018-04-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108168381B (zh) | 一种多枚导弹协同作战的控制方法 | |
CN107818219B (zh) | 一种面向突防的多导弹协同弹道规划方法 | |
CN104215935A (zh) | 一种基于决策加权融合的雷达炮弹目标识别方法 | |
CN106643298B (zh) | 一种基于预置命中点的大气层内反导拦截器中制导方法 | |
CN108763836A (zh) | 柱形弹丸撞击下空间碎片防护结构弹道极限方程获取方法 | |
Chusilp et al. | Investigating an iterative method to compute firing angles for artillery projectiles | |
CN106886625B (zh) | 一种基于固定翼鸭舵的双旋稳定弹的气动外形设计方法 | |
CN116360489A (zh) | 一种基于成型制导律的协同中制导律设计方法 | |
CN114819055A (zh) | 一种基于lstm网络的导弹拦截点预测方法 | |
CN112286059B (zh) | 一种具有攻击角和视场角约束的二阶前置角重塑制导方法 | |
Irwin et al. | Investigations into the aerodynamic properties of a battle damaged wing | |
CN112464451B (zh) | 基于作战仿真系统的防空导弹武器命中概率修正方法 | |
CN114153143A (zh) | 一种导弹非奇异固定时间滑模制导律的设计方法 | |
CN105987652A (zh) | 姿态角速率估算系统及应用其的弹药 | |
Chen et al. | Impact envelope simulation for kinetic kill missile intercepting short range high speed target | |
Li et al. | Achieving Interactive Shooting Training: Mathematical Modeling of Gun-launched Missile Ballistic Trajectories in LVC Architecture | |
CN114777577B (zh) | 一种微型战术导弹抗二级发动机推力扰动的控制方法 | |
CN107145761B (zh) | 带动力滑翔导弹多样式弹道规划方法 | |
Hong-He et al. | Modeling and Simulation of External Ballistics of New Air-Defence Guided Projectile | |
KR102136267B1 (ko) | 지대지 유도탄의 궤적성형을 위한 가표적 선정 방법 및 이를 수행하는 장치 | |
Yang et al. | Two-stage calculation model of single pulse vehicle tracking target and flight track | |
Guodong | The Study of the Modeling simulation for the Rocket-Assisted Cartridge | |
Nguyen et al. | EVALUATING EFFECTIVENESS OF EXTENDING ENGAGEMENT ZONE OF SURFACE-TO-AIR MISSILE USING THE COMBINATIONAL GUIDANCE LAW CV-2DGOC | |
Zhai et al. | Exterior Trajectory Computation and Modeling of Kinetic Energy Ammunition Based on Constant Terminal Effect | |
Li et al. | A Study of Aerodynamic Characteristics of an Anti-tank Missile |
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 |