CN111879178A - 一种提高卫星制导弹药制导精度的方法 - Google Patents

一种提高卫星制导弹药制导精度的方法 Download PDF

Info

Publication number
CN111879178A
CN111879178A CN202010782551.3A CN202010782551A CN111879178A CN 111879178 A CN111879178 A CN 111879178A CN 202010782551 A CN202010782551 A CN 202010782551A CN 111879178 A CN111879178 A CN 111879178A
Authority
CN
China
Prior art keywords
guidance
ammunition
tracker
satellite
function
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.)
Pending
Application number
CN202010782551.3A
Other languages
English (en)
Inventor
韩丁丁
郭强
卢朝林
董楠
王刚
李兴国
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xi'an Ruigao Measurement And Control Technology Co Ltd
Original Assignee
Xi'an Ruigao Measurement And Control Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xi'an Ruigao Measurement And Control Technology Co Ltd filed Critical Xi'an Ruigao Measurement And Control Technology Co Ltd
Priority to CN202010782551.3A priority Critical patent/CN111879178A/zh
Publication of CN111879178A publication Critical patent/CN111879178A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F41WEAPONS
    • F41GWEAPON SIGHTS; AIMING
    • F41G7/00Direction control systems for self-propelled missiles
    • F41G7/34Direction control systems for self-propelled missiles based on predetermined target position data
    • F41G7/346Direction control systems for self-propelled missiles based on predetermined target position data using global navigation satellite systems, e.g. GPS, GALILEO, GLONASS

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • General Engineering & Computer Science (AREA)
  • Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种提高卫星制导弹药制导精度方法,本发明的方法估计弹目视线角速率的结果相对真实值虽然也有发散,但发散幅度很小且与全程与真实值接近,这意味着制导指令与真实的需用指令相匹配,对弹体的制导精度有利;有效抑制了弹药的脱靶量,提高了制导精度,增强了弹药的平均毁伤能力。

Description

一种提高卫星制导弹药制导精度的方法
技术领域:
本实用导弹控制领域,涉及一种提高卫星制导弹药制导精度方法。
背景技术
卫星制导系统具有成本低、工作时间长、全天候工作、定位精度不受射程影响等特点,满足了现代战争的作战需求,因此广泛应用于现有的制导武器上。通过卫星制导系统,制导弹药可以实时的获取自身的位置和速度信息,利用这些信息通过制导律计算得到制导指令。目前工程上使用频率最高的速度追踪、比例导引、弹道成型等制导律的关键都是计算弹目视线角q或弹目视线角速率
Figure BDA0002620758330000011
目前工程的计算方法是通过卫星制导系统获得的制导弹药的位置以及预先装定的目标位置计算弹目视线角q,并直接微分得到弹目视线角速率
Figure BDA0002620758330000012
如果卫星制导系统没有噪声,这种计算方法的精度很高,但实际使用中卫星制导系统获得的位置和速度都是包含噪声的,计算得到弹目视线角q一定程度上已经偏离了真实值,而直接微分后将放大噪声的影响,恶化了弹目视线角速率
Figure BDA0002620758330000013
的精度,对制导精度的影响很大。因此,本专利提出一种提高卫星制导体制下弹目视线角及角速率信息精度的方法,通过快速跟踪微分器估计方法,对原卫星信号进行跟随和过滤,以计算得到精度更高的制导信息。
发明内容
本发明的目的在于克服上述现有技术的缺点,提供提高卫星制导弹药制导精度方法。本发明的方法估计弹目视线角速率的结果相对真实值虽然也有发散,但发散幅度很小且与全程与真实值接近,这意味着制导指令与真实的需用指令相匹配,对弹体的制导精度有利;有效抑制了弹药的脱靶量,提高了制导精度,增强了弹药的平均毁伤能力。
首先定义发射坐标系、弹目视线角及弹目视线角速率。
发射坐标系:
原点o——发射点
x轴——水平面内的发射方向
y轴——位于包含水平面内的发射方向的铅锤面内,垂直于x轴向上的方向
z轴——与x、y轴构成右手坐标系
弹目视线角:一般用q表示,表征弹目连线与发射坐标系xoz平面的夹角。弹目视线角是工程常用的速度追踪制导律的重要计算部分。
弹目视线角速率:一般用
Figure BDA0002620758330000021
表示,表征弹目视线角q在空间中旋转的速率。弹目视线角速率是工程常用的比例导引制导律的重要计算部分。
设定在发射坐标系中,制导弹药的真实实时位置为
Figure BDA0002620758330000022
真实实时速度
Figure BDA0002620758330000023
目标位置为(xT,yT,zT)。
计算中间量:
Figure BDA0002620758330000024
真实弹目视线角:
Figure BDA0002620758330000025
真实弹目视线角速率:
Figure BDA0002620758330000026
目前工程上卫星制导体制下的弹目视线角及角速率信息计算方法:
设定在发射坐标系中,卫星导航系统获得的弹药实时位置为
Figure BDA0002620758330000031
实时速度
Figure BDA0002620758330000032
目标位置为(xT,yT,zT)。
计算中间量:
Figure BDA0002620758330000033
计算弹目视线角:
Figure BDA0002620758330000034
计算弹目视线角速率:
Figure BDA0002620758330000035
由此可知,如果
Figure BDA0002620758330000036
包含噪声,qG
Figure BDA0002620758330000037
将偏离真实值qR
Figure BDA0002620758330000038
尤其是
Figure BDA0002620758330000039
由于微分的噪声放大特性,使得
Figure BDA00026207583300000310
在飞行的末段时间内波动剧烈,如果使用比例导引、弹道成型等需要使用弹目视线角速率的制导律时,
Figure BDA00026207583300000311
的剧烈波动经常使弹药在飞行终端脱靶。
现使用最速微分跟踪器对qG
Figure BDA00026207583300000312
的真实值进行估计:
Figure BDA00026207583300000313
Figure BDA00026207583300000314
算法中的符号表征意义如下:
qG(k)——第k个qG的解算值
Figure BDA00026207583300000315
——第k、k+1个qG的估计值
Figure BDA00026207583300000316
——第k、k+1个
Figure BDA00026207583300000317
的估计值
h——估计时间步长,对于卫星导航系统,h=0.1s
r——速度因子,可以取r=1~10
h0——滤波因子,可以取h0=2~20h
在最速微分跟踪器中,
Figure BDA00026207583300000318
函数公式如下:
Figure BDA0002620758330000041
Figure BDA0002620758330000042
a0=hx2
y0=x1+a0
Figure BDA0002620758330000043
Figure BDA0002620758330000044
a=(a0+y0)fsg(y0,d)+a2(1-fsg(y0,d))
Figure BDA0002620758330000045
函数中的符号表征意义如下:
x1、x2、d、y0、a1、a2、a——无特殊含义,是用于计算fhan函数结果所需的中间量
Figure BDA0002620758330000046
函数中的fsg(a,d)函数公式如下:
Figure BDA0002620758330000047
fsg(a,d)中的sgn(x)函数公式如下:
Figure BDA0002620758330000048
弹载计算机在执行最速微分跟踪器算法时,需要先对
Figure BDA0002620758330000049
的初值赋值,初值可以设为0。当卫星导航系统初始化完毕,开始正常输出信号时,最速微分跟踪器开始迭代计算。跟踪器需要一定的收敛时间 (<1s),对于一般的制导弹药,卫星导航系统输出信号到制导弹药开始控制的间隔时间较长(>3s),跟踪器的收敛时间对控制没有影响。弹载计算机可以使用跟踪器输出的信号
Figure BDA00026207583300000410
计算制导指令。
一、概念定义
Figure BDA00026207583300000411
Figure BDA0002620758330000051
二、符号列表
Figure BDA0002620758330000052
Figure BDA0002620758330000061
有益效果:
本发明的方法估计弹目视线角速率的结果相对真实值虽然也有发散,但发散幅度很小且与全程与真实值接近,这意味着制导指令与真实的需用指令相匹配,对弹体的制导精度有利;有效抑制了弹药的脱靶量,提高了制导精度,增强了弹药的平均毁伤能力。
附图说明:
图1为单次仿真的结果对比(飞行末段)图;
图2为图1飞行末段局部放大图;
图3为蒙特卡洛仿真的结果对比(目标点x=35050m,z=0m)图;
具体实施方式:
下面结合附图对本发明做进一步详细描述:
在实际应用中,如果制导过程正常,则真实弹目视线角速率的末端幅值在±3°/s,在工程常用的卫星导航系统的噪声级别下,末端弹目视线角速率的计算值可发散至±100°/s,由于该值发散所引起的脱靶量在总脱靶量中的占比超过60%。使用本专利方案可以将末端弹目视线角速率的发散程度抑制在±10°/s,远低于传统计算方法的发散程度,并将大大降低由于末端弹目视线角速率发散引起的脱靶量。卫星制导弹药的CEP值可以在原基础上降低40-50%。
使用某制导弹药数学模型对本专利的方案进行仿真。设该制导弹药使用卫星导航系统,定位误差符合正态分布N(0,(3.3m)2),测速误差符合正态分布N(0,(0.66m/s)2),制导律选择比例导引制导律。分别从单次仿真的结果和蒙特卡洛仿真的结果分析对比本专利方案的性能。图1表示了单次仿真中,飞行的最后24s的数据对比,具体对比了弹目视线角速率的真实值、使用直接微分算法的结果以及使用本专利方案得到的微分结果;图2对比了在蒙特卡洛实验中,使用直接微分算法和本专利方案得到的落点分布以及CEP值。
以上所述,仅是本发明的较佳实施例而已,并非对本发明作任何形式上的限制,虽然本发明已以较佳实施例揭露如上,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的方法及技术内容作出些许的更动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,仍属于本发明技术方案的范围内。

Claims (1)

1.一种提高卫星制导弹药制导精度方法,其特征在于:使用最速微分跟踪器对qG
Figure FDA0002620758320000011
的真实值进行估计:
Figure FDA0002620758320000012
Figure FDA0002620758320000013
算法中的符号表征意义如下:
qG(k)——第k个qG的解算值
Figure FDA0002620758320000014
——第k、k+1个qG的估计值
Figure FDA0002620758320000015
——第k、k+1个
Figure FDA0002620758320000016
的估计值
h——估计时间步长,对于卫星导航系统,h=0.1s
r——速度因子,可以取r=1~10
h0——滤波因子,可以取h0=2~20h
在最速微分跟踪器中,
Figure FDA0002620758320000017
函数公式如下:
Figure FDA0002620758320000018
Figure FDA0002620758320000019
a0=hx2
y0=x1+a0
Figure FDA00026207583200000110
Figure FDA00026207583200000111
a=(a0+y0)fsg(y0,d)+a2(1-fsg(y0,d))
Figure FDA00026207583200000112
函数中的符号表征意义如下:
x1、x2、d、y0、a1、a2、a——无特殊含义,是用于计算fhan函数结果所需的中间量;
Figure FDA00026207583200000113
函数中的fsg(a,d)函数公式如下:
Figure FDA0002620758320000021
fsg(a,d)中的sgn(x)函数公式如下:
Figure FDA0002620758320000022
弹载计算机在执行最速微分跟踪器算法时,需要先对
Figure FDA0002620758320000023
的初值赋值,初值可以设为0;当卫星导航系统初始化完毕,开始正常输出信号时,最速微分跟踪器开始迭代计算;跟踪器需要一定的收敛时间(<1s),对于一般的制导弹药,卫星导航系统输出信号到制导弹药开始控制的间隔时间较长(>3s),跟踪器的收敛时间对控制没有影响;弹载计算机可以使用跟踪器输出的信号
Figure FDA0002620758320000024
计算制导指令。
CN202010782551.3A 2020-08-06 2020-08-06 一种提高卫星制导弹药制导精度的方法 Pending CN111879178A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010782551.3A CN111879178A (zh) 2020-08-06 2020-08-06 一种提高卫星制导弹药制导精度的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010782551.3A CN111879178A (zh) 2020-08-06 2020-08-06 一种提高卫星制导弹药制导精度的方法

Publications (1)

Publication Number Publication Date
CN111879178A true CN111879178A (zh) 2020-11-03

Family

ID=73211479

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010782551.3A Pending CN111879178A (zh) 2020-08-06 2020-08-06 一种提高卫星制导弹药制导精度的方法

Country Status (1)

Country Link
CN (1) CN111879178A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114625002A (zh) * 2022-02-28 2022-06-14 浙江零跑科技股份有限公司 一种基于模型预测控制的车辆横纵向集成控制方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE2428704A1 (de) * 1973-06-19 1975-01-16 Emi Ltd Leit- und/oder fuehrungssystem
CA2793235A1 (en) * 2011-11-01 2013-05-01 Ge Aviation Systems Llc Methods for adjusting a relative navigation system
CN106406344A (zh) * 2016-09-28 2017-02-15 北京理工大学 一种旋转制导飞行器的增程系统及方法
CN108931155A (zh) * 2018-07-09 2018-12-04 北京航天控制仪器研究所 一种不依赖卫星导航增程制导弹药自主制导系统

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE2428704A1 (de) * 1973-06-19 1975-01-16 Emi Ltd Leit- und/oder fuehrungssystem
CA2793235A1 (en) * 2011-11-01 2013-05-01 Ge Aviation Systems Llc Methods for adjusting a relative navigation system
CN106406344A (zh) * 2016-09-28 2017-02-15 北京理工大学 一种旋转制导飞行器的增程系统及方法
CN108931155A (zh) * 2018-07-09 2018-12-04 北京航天控制仪器研究所 一种不依赖卫星导航增程制导弹药自主制导系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
王佩等: "跟踪-微分器在全捷联制导中的应用分析", 《兵器装备工程学报》 *
王涛等: "考虑阻力加速度的再入预测-校正制导算法", 《宇航学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114625002A (zh) * 2022-02-28 2022-06-14 浙江零跑科技股份有限公司 一种基于模型预测控制的车辆横纵向集成控制方法
CN114625002B (zh) * 2022-02-28 2024-04-23 浙江零跑科技股份有限公司 一种基于模型预测控制的车辆横纵向集成控制方法

Similar Documents

Publication Publication Date Title
CN110686564B (zh) 一种红外半捷联导引头制导方法及系统
WO2016118665A1 (en) Real-time ballistic solutions for calculating an aiming adjustment and for indicating a subsonic threshold
CN108279005A (zh) 一种导引头数据失效模式下的制导信息重构方法
CN111879178A (zh) 一种提高卫星制导弹药制导精度的方法
CN116592708A (zh) 适用于多种时空约束耦合下的飞行器弹道成型制导方法
CN110728010A (zh) 一种高炮火控系统射击诸元求取方法
CN113587740A (zh) 一种基于弹目视线角的被动反辐射导引方法及系统
RU2243482C1 (ru) Способ стрельбы боевой машины по цели и система для его реализации
CN115685778A (zh) 基于射频半实物仿真的雷达寻的飞行器脱靶量估算方法
CN111649734B (zh) 一种基于粒子群算法的捷联导引头目标定位方法
CN109376364A (zh) 基于扩展卡尔曼滤波的高速旋转弹气动参数辨识方法
RU2310152C1 (ru) Способ стрельбы боевой машины по цели и система для его реализации
Lim Predicting the accuracy of unguided artillery projectiles
Kerkhoff et al. The systemic error in the vertical component of handgun bullet trajectory reconstructions
CN112346474A (zh) 一种有限时间收敛的微分对策制导律的设计方法
CN113188507B (zh) 红外成像反坦克导弹弹目距离估计方法
CN116182647A (zh) 适应图像导引头的制导弹药带落角约束制导信息提取方法
CN114415724B (zh) 容许速度衰减的三维多导弹协同比例导引律设计方法
Zha et al. Study On the Dynamic Modeling and the Correction Method of the Self-Propelled Artillery
Abdin et al. Optimized Gain Proportional Navigation law for compensation a Delay in Line-Of-Sight Rate
RU2763897C1 (ru) Способ подготовки к выполнению задачи стрельбы на поражение из минометов
RU2692844C1 (ru) Способ повышения точности стрельбы артиллерийского вооружения боевой машины по цели (варианты) и система для его реализации
CN117889868B (zh) 一种融合红外导引头信息的导弹位置精确估计方法
CN118778660A (zh) 采用目标能观性主动增强制导律的制导方法
Xiaohui et al. Water column correction method of naval gun shooting at sea based on measuring radar

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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20201103