CN110703793B - 采用姿态角测量的飞行器积分比例导引攻击机动目标的方法 - Google Patents

采用姿态角测量的飞行器积分比例导引攻击机动目标的方法 Download PDF

Info

Publication number
CN110703793B
CN110703793B CN201911107385.0A CN201911107385A CN110703793B CN 110703793 B CN110703793 B CN 110703793B CN 201911107385 A CN201911107385 A CN 201911107385A CN 110703793 B CN110703793 B CN 110703793B
Authority
CN
China
Prior art keywords
aircraft
line
sight
angle
maneuvering target
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
Application number
CN201911107385.0A
Other languages
English (en)
Other versions
CN110703793A (zh
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.)
Naval Aeronautical University
Original Assignee
Naval Aeronautical University
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 Naval Aeronautical University filed Critical Naval Aeronautical University
Priority to CN201911107385.0A priority Critical patent/CN110703793B/zh
Publication of CN110703793A publication Critical patent/CN110703793A/zh
Application granted granted Critical
Publication of CN110703793B publication Critical patent/CN110703793B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course or altitude of land, water, air, or space vehicles, e.g. automatic pilot
    • G05D1/08Control of attitude, i.e. control of roll, pitch, or yaw
    • G05D1/0808Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft

Abstract

本发明是关于一种采用姿态角测量的飞行器积分比例导引攻击机动目标的方法,属于飞行器制导技术领域,该方法包括以下步骤:对飞行器的姿态角以及飞行器相对于机动目标运动的视线角速率进行测量;根据所述视线角速率构建视线角信号,并根据所述视线角信号构建积分型导引律;将所述积分型导引律的输出信号作为所述飞行器的期望航向角,使得姿态角能够稳定的对期望航向角进行跟踪,完成对所述机动目标的攻击。该方法解决了现有技术中由于比例导引规律末段指令偏大,进而导致脱靶量过大的问题,使得飞行器可以准确的攻击机动目标。

Description

采用姿态角测量的飞行器积分比例导引攻击机动目标的方法
技术领域
本发明涉及飞行器制导技术领域,具体而言,涉及一种采用姿态角测量的飞行器积分比例导引攻击机动目标的方法。
背景技术
飞行器末段精确制导律的研究具有较高的军事与民用价值。传统的比例导引由于末段脱靶量较大的原因,反而没有直接被采用。但基于比例导引思想的改进型导引方法研究,则引起工程研究者的广泛兴趣。
在很多传统设计中,依然以姿态稳定控制设计为主,因此如何在不引入过载控制回路的情况下,依然实现比例导引,也是一个具有很高工程价值的问题。
在目前已有的比例导引方法中,一般是直接根据导引头提供的视线角速率信号,然后转换为过载信号,提供给飞行器过载稳定回路进行过载跟踪,实现准确攻击目标的目的;但是采用该比例导引方法使得比例导引规律末段指令偏大,进而导致脱靶量过大。
需要说明的是,在上述背景技术部分发明的信息仅用于加强对本发明的背景的理解,因此可以包括不构成对本领域普通技术人员已知的现有技术的信息。
发明内容
本发明的目的在于提供一种采用姿态角测量的飞行器积分比例导引攻击机动目标的方法,进而至少在一定程度上克服由于相关技术的限制和缺陷而导致的比例导引规律末段指令偏大,进而导致脱靶量过大的问题。
根据本发明的一个方面,提供一种采用姿态角测量的飞行器积分比例导引攻击机动目标的方法,包括以下步骤:
步骤S10,对飞行器的姿态角以及飞行器相对于机动目标运动的视线角速率进行测量;
步骤S20,根据所述视线角速率构建视线角信号,并根据所述视线角信号构建积分型导引律;
步骤S30,将所述积分型导引律的输出信号作为所述飞行器的期望航向角,使得姿态角能够稳定的对期望航向角进行跟踪,完成对所述机动目标的攻击。
在本发明的一个示例实施例中,对飞行器相对于机动目标运动的视线角速率进行测量包括:
Figure BDA0002271729540000021
Figure BDA0002271729540000022
Figure BDA0002271729540000023
其中,
Figure BDA0002271729540000024
为视线角速率,
Figure BDA0002271729540000025
为qε的导数;x、y、z为所述飞行器在空间三维坐标系上的位置坐标;xT、yT、zT为所述机动目标在三维坐标系上的位置坐标,R为飞行器与所述机动目标之间的距离。
在本发明的一个示例实施例中,根据所述视线角速率构建视线角信号包括:
Figure BDA0002271729540000026
其中,u1为视线角信号;
Figure BDA0002271729540000027
为视线角速率;∫dt为积分符号。
在本发明的一个示例实施例中,根据所述视线角信号构建积分型导引律包括:
Figure BDA0002271729540000031
其中,u2=∫u1dt;
k1、k2、k3、k4、ε为常参数。
本发明一种采用姿态角测量的飞行器积分比例导引攻击机动目标的方法,一方面,根据飞行器导引头测量得到的视线角速率,并根据视线角速率信号构建积分型比例导引规律,然后将导引信号转换为姿态信号,提供给飞行器姿态稳定回路进行跟踪,解决了现有技术中由于比例导引规律末段指令偏大,进而导致脱靶量过大的问题,使得飞行器可以准确的攻击机动目标;另一方面,可以实现准确攻击机动目标,能够解决采用一套传统的姿态控制体系实现飞行器全部飞行控制与制导的任务,因此具有很高的工程实用价值。
应当理解的是,以上的一般描述和后文的细节描述仅是示例性和解释性的,并不能限制本发明。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分,示出了符合本发明的实施例,并与说明书一起用于解释本发明的原理。显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得与本发明的附图类似的其他的附图。
图1是本发明提供的一种采用姿态角测量的飞行器积分比例导引攻击机动目标方法的流程图。
图2是本发明实施例所提供方法的飞行器与目标在航向平面相对运动曲线。
图3是本发明实施例所提供方法的脱靶量曲线。
图4是本发明实施例所提供方法的脱靶量放大曲线。
图5是本发明实施例所提供方法的实际姿态角与期望姿态角的对比曲线。
图6是本发明实施例所提供方法的积分型比例导引律的输出。
具体实施方式
现在将参考附图更全面地描述示例实施方式。然而,示例实施方式能够以多种形式实施,且不应被理解为限于在此阐述的范例;相反,提供这些实施方式使得本发明将更加全面和完整,并将示例实施方式的构思全面地传达给本领域的技术人员。所描述的特征、结构或特性可以以任何合适的方式结合在一个或更多实施方式中。在下面的描述中,提供许多具体细节从而给出对本发明的实施方式的充分理解。在其它情况下,不详细示出或描述公知技术方案以避免喧宾夺主而使得本发明的各方面变得模糊。
本发明提供了一种采用姿态角测量的飞行器积分比例导引攻击机动目标的方法,根据飞行器导引头测量得到的视线角速率,考虑末端角速率信号的盲区问题,根据视线角速率信号构建积分型比例导引规律,然后将导引信号转换为姿态信号,提供给飞行器姿态稳定回路进行跟踪,实现飞行器准确攻击机动目标。
在目前已有的比例导引方法中,一般是直接根据导引头提供的视线角速率信号,然后转换为过载信号,提供给飞行器过载稳定回路进行过载跟踪,实现准确攻击目标的目的。但有的飞行器系统采用传统的姿态控制体制,因此没有设计过载稳定回路,为了末段导引单独引入过载稳定回路又代价太大,如存在过载精确测量器件昂贵与过载稳定回路设计复杂的问题。同时传统的比例导引规律末段指令偏大导致脱靶量过大,本发明提供了积分型比例导引,能够提高攻击目标的精确度。因此本发明提供了采用积分比例导引与姿态角测量稳定跟踪的方法,实现准确攻击机动目标,能够解决采用一套传统的姿态控制体系实现飞行器全部飞行控制与制导的任务,因此具有很高的工程实用价值。
下面,将结合附图对本发明的一种采用姿态角测量的飞行器积分比例导引攻击机动目标的方法进行进一步的解释以及说明。参考图1所示,该采用姿态角测量的飞行器积分比例导引攻击机动目标的方法可以包括以下步骤:
步骤S10,对飞行器的姿态角、姿态角速率以及相对目标运动的视线角速率进行测量。
具体的,首先,可以采用惯性导航组合测量飞行器的姿态角信号,下面以航向通道与航向平面导引为研究对象。假设采用惯性导航组合元器件测量得到飞行器航向通道的姿态角,记作ψc。同时测量其姿态角速率,记作为ωy;也可在飞行器上安装姿态角度测量陀螺仪与速率陀螺仪,分别测量姿态角ψc与姿态角速率ωy
然后,采用导引头测量飞行器相对目标运动的视线角速率,记作
Figure BDA0002271729540000054
其数学上表示视线角qε的导数,其中qε的定义为:
Figure BDA0002271729540000051
Figure BDA0002271729540000052
Figure BDA0002271729540000053
其中,x、y、z为飞行器在空间三维坐标系上的位置坐标;xT、yT、zT为目标在三维坐标系上的位置坐标。
步骤S20,根据相对目标运动的视线角速率构建视线角信号,并根据视线角信号构建积分型导引律。
具体的,首先,根据视线角速率
Figure BDA0002271729540000061
构造视线角信号u1
其中,
Figure BDA0002271729540000062
∫dt为积分符号。
此处需要补充说明的是,此时的u1信号由于是由测量信号积分而成,而测量信号不可避免地带有延迟与噪声,因此u1与理想的视线角信号qε不同,主要是包含延迟与噪声。
其次,在视线角信号u1的基础上,再积分,构造u2信号如下:
u2=∫u1dt;
进一步的,再根据u2设计积分型导引律u:
Figure BDA0002271729540000063
其中,k1、k2、k3、k4以及ε为常参数,具体的选取详见后文案例实施。
考虑到导引头测量元器件的盲区,将积分型导引律改进处理为:
Figure BDA0002271729540000064
其中,ua0为R=20时刻ua的值;R为目标与飞行器之间的距离。
步骤S30,将所述积分型导引律的输出信号作为所述飞行器的期望航向角,使得姿态角能够稳定的对期望航向角进行跟踪,完成对所述机动目标的攻击。
具体的,此处主要涉及的是比例导引信号与期望姿态信号的转换方法,在此直接将上述积分型比例导引信号的输出u作为飞行器姿态稳定回路的输入信号,也就是飞行器的期望航向角
Figure BDA0002271729540000065
使得飞行器的实际航向姿态角ψc能够稳定跟踪期望航向角
Figure BDA0002271729540000066
也就是稳定跟踪积分型比例导引律的输出u,即设置
Figure BDA0002271729540000071
进一步的,下文以常规的姿态稳定回路为例说明姿态稳定回路跟踪律的设计过程与飞行器姿态稳定跟踪飞行的过程。
dyc=ka1e+ka2∫edt+ka3wy
其中,dyc为飞行器偏航通道的控制量,即偏航舵偏指令信号。e为误差信号,其定义为:
Figure BDA0002271729540000072
edt为误差信号的积分,ψc为步骤一所测量的姿态角信号,wy为步骤一所测量的姿态角速率信号。参数ka1与ka2、ka3的选取和飞行器的空气动力学特性有关,在此不是本专利要涉及与保护的内容,其选取与设计方法不作详细介绍。dyc即为最终飞行器航向通道的控制输出,用于控制飞行器航向舵,使得飞行器运动轨迹发生变化,从而命中目标。
最后,是关于本发明的目标模拟与参数调整。
具体的,可以选取需要攻击的目标运动参数,模拟不同速度、不同初始位置情况的机动目标与飞行器末段相遇的场景,进行上述积分导引与姿态匹配的制导规律参数设计,通过多种目标运动模拟,判断不同制导律参数情况下的脱靶量大小,从而选定最终的导引律参数。
其中脱靶量的定义为如下的简单形式:
当Δx=xT-x<0时刻,飞行器与目标的距离R即为脱靶量。
最终,依靠多次仿真模拟,挑选满意的参数,即组成最终的积分型导引律方法。
案例实施与计算机仿真模拟结果分析
首先,选取k1=3、k2=0.2、k3=0.3、k4=0.2、ε=0.2。
然后,选取ka1=2.2与ka2=0.7、ka3=0.8。
以某型陆地机动目标为例,说明案例实施的过程。假设初始的目标位置为xT(0)=5000、yT(0)=1以及zT(0)=300,此时目标匀速运动,其速度大小为20m/s,其方向为与x轴正方向成20度夹角。其它初始位置与目标速度大小、速度方向的情况,参数挑选原理相同,在此不一一说明。
从图2可以看出,在航向平面内,飞行器与目标能够相互接近;由图3可以看出,飞行器与目标间的相对距离是逐渐减小的;由图4可以看出,最终的脱靶量为3.1米,这对于尺寸大于5米的目标,如船舶等是完全满足要求的。由图5可以看出,实际姿态角和期望姿态角是基本一致的,而且角度小于60度,因此可以看出该类积分型比例导引律设计是合理了,飞行器的姿态稳定跟踪回路能够完全实现上述期望姿态角,从而使得飞行器最终命中目标。由图6可以看出,导引律的输出值初始段比较小,而末段较大,主要是由于离目标较近导致的,也是合理的。
在上述基础上,考虑目标运动速度改变、或者目标初始位置的变化,对上述导引律参数进行微调,最终确定整套积分型导引律的参数,从而完成积分型比例导引与姿态跟踪的匹配设计。
本领域技术人员在考虑说明书及实践本发明后,将容易想到本发明的其他实施例。本申请旨在涵盖本发明的任何变型、用途或者适应性变化,这些变型、用途或者适应性变化遵循本发明的一般性原理并包括本发明的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的,本发明的真正范围和精神由权利要求指出。

Claims (1)

1.一种采用姿态角测量的飞行器积分比例导引攻击机动目标的方法,其特征在于,包括以下步骤:
步骤S10,对飞行器的姿态角以及飞行器相对于机动目标运动的视线角速率进行测量;
步骤S20,根据所述视线角速率构建视线角信号,并根据所述视线角信号构建积分型导引律;
步骤S30,将所述积分型导引律的输出信号作为所述飞行器的期望航向角,使得姿态角能够稳定的对期望航向角进行跟踪,完成对所述机动目标的攻击;
其中,对飞行器相对于机动目标运动的视线角速率进行测量包括:
Figure FDA0003530295510000011
Figure FDA0003530295510000012
Figure FDA0003530295510000013
其中,
Figure FDA0003530295510000014
为视线角速率,
Figure FDA0003530295510000015
为qε的导数;x、y、z为所述飞行器在空间三维坐标系上的位置坐标;xT、yT、zT为所述机动目标在三维坐标系上的位置坐标,R为飞行器与所述机动目标之间的距离;
根据所述视线角速率构建视线角信号包括:
Figure FDA0003530295510000016
其中,u1为视线角信号;
Figure FDA0003530295510000017
为视线角速率;∫dt为积分符号;
根据所述视线角信号构建积分型导引律包括:
Figure FDA0003530295510000018
其中,u2=∫u1dt;
k1、k2、k3、k4、ε为常参数。
CN201911107385.0A 2019-11-13 2019-11-13 采用姿态角测量的飞行器积分比例导引攻击机动目标的方法 Active CN110703793B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911107385.0A CN110703793B (zh) 2019-11-13 2019-11-13 采用姿态角测量的飞行器积分比例导引攻击机动目标的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911107385.0A CN110703793B (zh) 2019-11-13 2019-11-13 采用姿态角测量的飞行器积分比例导引攻击机动目标的方法

Publications (2)

Publication Number Publication Date
CN110703793A CN110703793A (zh) 2020-01-17
CN110703793B true CN110703793B (zh) 2022-06-24

Family

ID=69205941

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911107385.0A Active CN110703793B (zh) 2019-11-13 2019-11-13 采用姿态角测量的飞行器积分比例导引攻击机动目标的方法

Country Status (1)

Country Link
CN (1) CN110703793B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111399529B (zh) * 2020-04-02 2021-05-14 上海交通大学 一种基于非线性滑模与前置的飞行器复合导引方法
CN112162569B (zh) * 2020-09-09 2022-02-18 北京航空航天大学 一种飞行器绕多禁飞区路径规划与决策方法
CN115493589B (zh) * 2022-11-04 2023-03-28 开拓导航控制技术股份有限公司 视线角解算方法和系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103994699A (zh) * 2014-05-23 2014-08-20 中国人民解放军海军航空工程学院 一种超声速导弹爬升弹道过载指令自动生成的方法
CN107132765A (zh) * 2017-06-01 2017-09-05 烟台南山学院 一种基于轨迹规划的攻击角度与攻击时间控制方法
CN107726921A (zh) * 2017-08-30 2018-02-23 湖北航天技术研究院总体设计所 一种适用于目标机动条件下的主动雷达角制导方法
CN109597423A (zh) * 2019-01-08 2019-04-09 北京航空航天大学 一种基于参考视线角信号的多约束末制导律的设计方法
CN109814095A (zh) * 2019-01-23 2019-05-28 北京仿真中心 一种多目标红外仿真系统的动态空间位置模拟方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6653972B1 (en) * 2002-05-09 2003-11-25 Raytheon Company All weather precision guidance of distributed projectiles

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103994699A (zh) * 2014-05-23 2014-08-20 中国人民解放军海军航空工程学院 一种超声速导弹爬升弹道过载指令自动生成的方法
CN107132765A (zh) * 2017-06-01 2017-09-05 烟台南山学院 一种基于轨迹规划的攻击角度与攻击时间控制方法
CN107726921A (zh) * 2017-08-30 2018-02-23 湖北航天技术研究院总体设计所 一种适用于目标机动条件下的主动雷达角制导方法
CN109597423A (zh) * 2019-01-08 2019-04-09 北京航空航天大学 一种基于参考视线角信号的多约束末制导律的设计方法
CN109814095A (zh) * 2019-01-23 2019-05-28 北京仿真中心 一种多目标红外仿真系统的动态空间位置模拟方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Cooperative Control of Multi-UAV with Time Constraint in The Threat Environment;Pei-bei, M.et al.;《Proceedings of 2014 IEEE Chinese Guidance, Navigation and Control Conference》;20140831;第2424-2428页 *
基于变结构控制理论的导弹平滑导引律研究;李静等;《海军航空工程学院学报》;20070531;第22卷(第3期);第313-316页 *
多导弹攻击时间和攻击角度协同制导研究综述;马培蓓等;《飞航导弹》;20180630;第59-63页 *
带落角约束的超声速飞行器制导控制一体化设计;王建华等;《电机与控制学报》;20170228;第21卷(第2期);第76-86页 *
带落角限制的虚拟目标比例导引律设计;顾文锦等;《飞行力学》;20060630;第24卷(第2期);第43-46页 *

Also Published As

Publication number Publication date
CN110703793A (zh) 2020-01-17

Similar Documents

Publication Publication Date Title
CN110703793B (zh) 采用姿态角测量的飞行器积分比例导引攻击机动目标的方法
CN110764534B (zh) 基于非线性转换的前置导引与姿态稳定匹配制导的方法
CN110764523B (zh) 基于抗饱和平滑变换的比例积分前置导引攻击目标的方法
CN109740198B (zh) 一种基于解析预测的滑翔飞行器三维再入制导方法
CN110926278B (zh) 一种采用多前置角叠加与末端修正的前置导引方法
CN110687931B (zh) 一种方位角姿态与前置导引切换的一体化机动导引方法
CN110044321B (zh) 利用地磁信息和角速率陀螺解算飞行器姿态的方法
CN109445449B (zh) 一种高亚音速无人机超低空飞行控制系统及方法
CN110032206B (zh) 远程制导飞行器大落角攻顶控制方法及控制系统
CN110425943B (zh) 面向变质心飞行器的工程化再入制导方法
CN111324149A (zh) 一种基于视线角信息与前置导引信息的复合制导方法
CN107102547B (zh) 一种基于滑模控制理论的rlv着陆段制导律获取方法
CN110017808B (zh) 利用地磁信息和加速计解算飞行器姿态的方法
CN113359856B (zh) 一种无人机指定航向目标点导引方法及系统
CN111102886B (zh) 一种小微飞行器无陀螺精确制导方法
CN114489101A (zh) 一种无人机的末端制导控制方法及系统
CN110823016B (zh) 一种转捩研究用高精度三维空间制导方法
CN110879604A (zh) 一种带落角控制的飞行器航向导引方法
CN114815888B (zh) 一种仿射形式的制导控制一体化控制方法
CN110972541B (zh) 一种提高脱靶平面内弹目相对高度的导引方法
CN110017809B (zh) 利用地磁信息和光流传感器解算飞行器姿态的方法
CN111913491B (zh) 一种基于视线角非线性抗饱和与不确定性补偿的制导方法
CN111474948B (zh) 一种带时间控制的前置导引与姿态控制制导的方法
CN114020019A (zh) 飞行器的制导方法与装置
CN113587740A (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