CN112540618B - 高动态旋转飞行器复合制导控制系统及方法 - Google Patents

高动态旋转飞行器复合制导控制系统及方法 Download PDF

Info

Publication number
CN112540618B
CN112540618B CN201910894011.1A CN201910894011A CN112540618B CN 112540618 B CN112540618 B CN 112540618B CN 201910894011 A CN201910894011 A CN 201910894011A CN 112540618 B CN112540618 B CN 112540618B
Authority
CN
China
Prior art keywords
aircraft
target
guidance
module
line
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
CN201910894011.1A
Other languages
English (en)
Other versions
CN112540618A (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.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201910894011.1A priority Critical patent/CN112540618B/zh
Publication of CN112540618A publication Critical patent/CN112540618A/zh
Application granted granted Critical
Publication of CN112540618B publication Critical patent/CN112540618B/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, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • 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
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • 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
    • G05D1/0816Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/10Simultaneous control of position or course in three dimensions
    • G05D1/101Simultaneous control of position or course in three dimensions specially adapted for aircraft

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本发明公开了一种高动态旋转飞行器复合制导控制系统及方法,所述系统包括制导模块(1)、飞行器姿态测量模块(2)、舵机指令解算模块(3)和执行模块(4),其中,所述制导模块采用雷达/红外的双模复合导引模式,在飞行器飞行的过程中根据飞行器与目标的距离自动切换制导模式,能够实现全天候、多环境的使用,提高了制导精度;通过设置捷联导引头和视线角速度解算模块,提高了舱体的空间利用率,降低了干扰性;另外,通过对控制舱进行了一体化设计,增加了战斗部空间,提高了飞行品质。

Description

高动态旋转飞行器复合制导控制系统及方法
技术领域
本发明涉及旋转飞行器制导控制领域,具体涉及高动态旋转飞行器复合制导控制系统及方法。
背景技术
传统的高动态旋转飞行器多为单一制导体制,可在常规环境条件下稳定地完成制导任务,而无法满足全天候、大跨域的战场环境要求,且单一制导体制的飞行器容易受到外部干扰,导致制导精度降低,甚至造成脱靶等不利结果。
此外,现有高动态飞行器的制导控制舱部件较多,舱内结构复杂,舱体体积大,影响战斗部的体积,降低了飞行器的威力。
因此,有必要提供一种高动态旋转飞行器复合制导控制系统及方法,其能够实现全天候、多环境的使用,制导精度高,且舱体结构简单。
发明内容
为了克服上述问题,本发明人进行了锐意研究,设计出一种高动态旋转飞行器复合制导控制系统及方法,该系统采用雷达/红外的双模复合导引模式,在飞行器飞行的过程中根据飞行器与目标的距离自动切换制导模式,能够实现全天候、多环境的使用,提高了制导精度;通过设置捷联导引头和视线角速度解算模块,提高了舱体的空间利用率,降低了干扰性;另外,通过对控制舱进行了一体化设计,增加了战斗部空间,提高了飞行品质,从而完成了本发明。
具体来说,本发明的目的在于提供以下方面:
第一方面,提供一种高动态旋转飞行器复合制导控制系统,
所述系统包括制导模块1、飞行器姿态测量模块2、舵机指令解算模块3和执行模块4;其中,
所述制导模块1用于实时获取飞行器与目标的位置信息,并传递至舵机指令解算模块3,
所述飞行器姿态测量模块2用于实时获取飞行器的飞行参数信息,并传递至舵机指令解算模块3,
所述舵机指令解算模块3用于实时根据获取的飞行器与目标的位置信息及飞行器姿态信息获得打舵指令,
所述执行模块4用于实时根据接收到的舵机指令调整飞行器姿态。
第二方面,提供一种高动态旋转飞行器复合制导控制方法,优选采用第一方面所述的系统进行,其中,
步骤1,通过制导模块1中的雷达导引系统11实时获得飞行器与目标之间的距离,由制导模式选择模块5选择制导模式;
步骤2,通过导引头6实时获得飞行器与目标的视线角,并由视线角速度解算模块7实时解算得到飞行器与目标的视线角速度,并传递至舵机指令解算模块3;
步骤3,通过飞行器姿态测量模块2实时获得飞行器的飞行参数信息,并传递至舵机指令解算模块3;
步骤4,通过舵机指令解算模块3实时获得打舵指令,并传递至执行模块4;
步骤5,通过执行模块4实时调整飞行器的姿态。
本发明所具有的有益效果包括:
(1)本发明提供的高动态旋转飞行器复合制导控制系统,采用雷达导引系统和红外导引系统的双模复合制导模式,可以根据飞行器与目标的距离进行自动选择,实现了全天候、多环境的使用,提高了制导精度;
(2)本发明提供的高动态旋转飞行器复合制导控制系统,采用捷联导引头和视线角速度解算模块结合的方式实时获取飞行器与目标的视线角速度,能够节省控制舱的空间,增加战斗部的空间,提升飞行器的空间利用率,抗干扰性强,提高飞行品质,降低脱靶量;
(3)本发明提供的高动态旋转飞行器复合制导控制系统,获得的滚转角精确性高,显著提高了制导精度;
(4)本发明提供的高动态旋转飞行器复合制导控制方法,可以满足全天候、大跨域的战场要求,制导过程不易受到内外部干扰的影响,稳定性高。
附图说明
图1示出本发明一种优选实施方式的高动态旋转飞行器复合制导控制系统的整体结构示意图;
图2示出本发明实验例1中不考虑外部扰动影响时,计算机实时给出的飞行器与目标的视线角速度的真实值与解算值之间的比较结果;
图3示出本发明实验例1中考虑外部扰动影响时,计算机实时给出的飞行器与目标的视线角速度的真实值与解算值之间的比较结果;
图4示出本发明实验例2中不考虑外部扰动影响时,计算机实时给出的飞行器与目标的视线角速度的真实值与解算值之间的比较结果;
图5示出本发明实验例2中考虑外部扰动影响时,计算机实时给出的飞行器与目标的视线角速度的真实值与解算值之间的比较结果;
图6示出本发明实验例3中不考虑外部扰动影响时,计算机实时给出的飞行器与目标的视线角速度的真实值与解算值之间的比较结果;
图7示出本发明实验例3中考虑外部扰动影响时,计算机实时给出的飞行器与目标的视线角速度的真实值与解算值之间的比较结果。
附图标号说明:
1-制导模块;
11-雷达导引系统;
12-红外导引系统;
2-飞行器姿态测量模块;
21-地磁传感器;
22-惯性陀螺;
3-舵机指令解算模块;
4-执行模块;
5-制导模式选择模块;
6-导引头;
7-视线角速度解算模块;
8-电源供电模块。
具体实施方式
下面通过附图和实施方式对本发明进一步详细说明。通过这些说明,本发明的特点和优点将变得更为清楚明确。其中,尽管在附图中示出了实施方式的各种方面,但是除非特别指出,不必按比例绘制附图。
本发明提供了一种高动态旋转飞行器复合制导控制系统,如图1所示,所述系统包括制导模块1、飞行器姿态测量模块2、舵机指令解算模块3和执行模块4。
其中,所述高动态是指飞行器的转速高,一般转速在10r/s以上。
根据本发明一种优选的实施方式,所述制导模块1用于实时获取飞行器与目标的位置信息,并传递至舵机指令解算模块3,
所述制导模块1采用复合制导模式,以提高指导精度。
其中,所述飞行器与目标的位置信息包括飞行器与目标的视线角信息以及飞行器与目标的距离。
在进一步优选的实施方式中,所述制导模块1包括雷达导引系统11和红外导引系统12。
其中,所述雷达导引系统由天线、接收机、驱动装置、基准电压发生器、相位比较器组成。其中,天线接收目标反射的雷达回波,当天线轴正好对准目标时输出为零。如果天线未对正目标,则有一个与失调角成比例的信号输出,此信号经接收机接收放大后进入相位比较器。同时,由基准电压发生器发射的与天线同频率的基准信号,经过放大后同时进入相位比较器。相位比较器根据基准信号确定出偏差分量,雷达导引系统根据偏差分量确定飞行器与目标的视线角信息,根据雷达回波时间确定飞行器与目标间的距离。
所述红外导引系统主要由红外探测系统和电子线路两部分组成,是一个使光学系统跟踪目标的机电装置,其将红外信号转换为光信号,光信号通过导引头中的光电转换器转化为电信号,通过电子线路将电信号传递至舵机指令解算模块。
在本发明中,对雷达导引系统和红外导引系统的结构不作特别限定,可以为现有技术中常见的具有上述功能的雷达探测器和红外探测器。
本发明人经过研究发现,雷达导引系统虽然可以全天时、全天候工作,但容易受到侦查干扰且测量精度低;红外导引系统虽测量精度高,抗干扰能力强,但是受气候影响较大,云雾天气探测距离急剧减小。
为解决上述问题,根据本发明一种优选的实施方式,所述高动态旋转飞行器的复合制导系统还包括制导模式选择模块5,其根据飞行器与目标之间的距离选择制导模式。
其中,飞行器在飞向目标的过程中,其上的敏感元件及舵机等期间开始上电工作时称为启控,在启控后的飞行阶段称之为制导段,所述制导段包括中制导段和末制导段,一般来说,带有激光导引头的制导飞行器在距离目标3km,并且激光导引头开始捕获目标时,飞行器进入末制导段。
在进一步优选的实施方式中,所述制导模式选择模块5按照以下规则选择制导模式:
当飞行器与目标之间的距离小于等于10Km且大于2.99Km时,选择雷达导引系统进行制导;
当飞行器与目标之间的距离小于等于2.99Km且大于0.99Km时,选择雷达导引系统和红外导引系统进行复合制导;
当飞行器与目标之间的距离小于等于0.99km时,选择雷达导引系统和红外导引系统进行复合制导,其中,飞行器与目标的视线角信息由红外导引系统获得,飞行器与目标之间的距离由雷达导引系统获得。
其中,当飞行器与目标之间的距离小于等于2.99Km且大于0.99Km时,选择雷达导引系统和红外导引系统进行复合制导,此时的飞行器与目标的视线角由两种制导模式下得到的视线角加权得到,飞行器与目标之间的距离由雷达导引系统获得。
在更进一步优选的实施方式中,当飞行器与目标之间的距离小于等于2.99Km且大于0.99Km时,复合制导模式下的视线角通过下式(一)获得:
Figure BDA0002209658140000071
其中,
Figure BDA0002209658140000072
表示复合制导模式下的飞行器与目标的视线角,q1表示雷达导引系统进行制导模式下的飞行器与目标的视线角,q2表示红外导引系统进行制导模式下的飞行器与目标的视线角,a为复合制导模式下雷达导引系统所占权重,b为复合制导模式下红外导引系统所占权重。
优选地,所述a和b分别由下式(二)和式(三)获得:
a=0.5x-0.5(二);
b=1-a(三)
其中,x表示飞行器与目标之间的距离。
根据本发明一种优选的实施方式,在所述制导模式选择模块5与舵机指令解算模块3之间还设置有导引头6,所述导引头6为捷联导引头。
其中,与传统的平台导引头相比,捷联导引头可直接固连于飞行器上,具有体积小、抗干扰性强、稳定性好等优点,可适用于各类高动态飞行器中。
本发明人研究发现,捷联导引头实现视场稳定的难度较大,只能得到飞行器与目标的视线角,而很难得到飞行器与目标的视线角速度,由于飞行器与目标的视线角速度作为制导指令中重要的一环,如果计算精度低将导致后续控制精度降低,从而影响飞行品质,甚至提高脱靶量。
因此,在进一步优选的实施方式,在所述导引头6与舵机指令解算模块3之间还设置有视线角速度解算模块7,以根据接收到的视线角信息实时解算得到飞行器与目标的视线角速度。
在更进一步优选的实施方式中,所述视线角速度解算模块7通过下式(四)、(五)和(六)实时获得飞行器与目标的视线角速度信息:
Figure BDA0002209658140000081
Figure BDA0002209658140000082
Figure BDA0002209658140000083
其中,qg表示由导引头获得的飞行器与目标的视线角;q0表示飞行器与目标的视线角的估测值,即在解算过程中通过上述式(四)、(五)、(六)估算出来的飞行器与目标的视线角的估测值;q01表示飞行器与目标的视线角速度估测值,即在解算过程中通过上述式(四)、(五)、(六)估算出来的飞行器与目标的视线角速度的估测值;
Figure BDA0002209658140000084
表示x2的导数;
Figure BDA0002209658140000085
表示x1的导数;
Figure BDA0002209658140000086
表示x0的导数,上一时刻得到的数值作为下一时刻迭代的初始值;
初始时刻时,x0=0,x1=0,x2=0,每隔0.001s作为积分步长,进行迭代,并且得到x0、x1和x2在下一时刻的数值;
具体来说,第一个迭代时刻中,将初始时刻x0=0,x1=0,x2=0和接收到的qg值代入到式(四)、(五)、(六)中,从而解算出
Figure BDA0002209658140000087
进而得到下一时刻的初始值x0、x1和x2;再将得到的x0、x1、x2和接收到的qg值代入到式(四)、(五)、(六)中,从而得到再下一时刻对应的
Figure BDA0002209658140000088
如此持续循环迭代即可持续得到每次积分得到的对应x0、x1、x2
其中,
Figure BDA0002209658140000089
表示飞行器与目标的视线角速度,实时输出该飞行器与目标的视线角速度给舵机指令解算模块3即可用于解算需用过载。
其中,所述a0、a1、a2、δ、k1和k2均为设计参数,优选地,a0=1~1.5、a1=7~10、a2=10~15、δ=1~2、k1=0.1~0.4和k2=0.2~0.4;
更优选地,所述a0=1.1、a1=8.5、a2=11.5、δ=1.5、k1=0.3、k2=0.3。
在本发明中,当存在外部扰动时,捷联导引头对飞行器与目标的视线角的测量便存在一定的误差,而所有的干扰均可等效为飞行器与目标的视线角的增量。
在本发明中,所述的需用过载指的是用以控制飞行器上舵机工作量的指标数据,飞行器上的舵机根据计算得到的需用过载打舵工作。
根据本发明一种优选的实施方式,所述飞行器姿态测量模块2用于实时获取飞行器的飞行参数信息,并传递至舵机指令解算模块3。
在本发明中,所述飞行参数信息包括飞行器的滚转角信息、俯仰角信息和偏航角信息。
在进一步优选的实施方式中,所述飞行器姿态测量模块2包括地磁传感器21和惯性陀螺22,其中,
所述地磁传感器21用于实时获得飞行器的滚转角信息,
所述惯性陀螺22用于实时获得飞行器的俯仰角信息和偏航角信息。
其中,地磁传感器为磁阻传感器,基于磁铁材料具有磁阻效应的原理工作。主要由镍铁导磁合金组成,通常组成惠斯通电桥来感应外界磁场。具有灵敏度高、响应时间短等特点,可为高速运动的载体提供测量信息。
本发明人研究发现,地磁传感器输出地磁信号在载体坐标系下的三轴分量,结合载体坐标系与导航坐标系的转换关系、惯性陀螺测得的飞行器俯仰角及偏航角信息,可以获得飞行器精确的滚转角信息,即输出滚转角信息,进而传递至舵机指令解算模块3。
在更进一步优选的实施方式中,所述输出滚转角信息通过下式(七)获得:
Figure BDA0002209658140000101
其中,Hx*表示地磁信号在载体坐标系的x轴分量;Hy*表示地磁信号在载体坐标系的y轴分量;Hz*表示地磁信号在载体坐标系的z轴分量;Hx表示地磁信号在导航坐标系中的x轴分量;Hy表示地磁信号在导航坐标系中的y轴分量;Hz表示地磁信号在导航坐标系中的z轴分量;其中,Hx*、Hy*、Hz*、Hx、Hy、Hz由地磁传感器直接测得;ψ表示偏航角,θ表示俯仰角,均由惯性陀螺直接测得;γ表示滚转角。
其中,载体坐标系与导航坐标系的转换关系如下式(八)所示:
Figure BDA0002209658140000102
其中,Hx*表示地磁信号在载体坐标系的x轴分量;Hy*表示地磁信号在载体坐标系的y轴分量;Hz*表示地磁信号在载体坐标系的z轴分量;Hx表示地磁信号在导航坐标系中的x轴分量;Hy表示地磁信号在导航坐标系中的y轴分量;Hz表示地磁信号在导航坐标系中的z轴分量;ψ表示偏航角;θ表示俯仰角;γ表示滚转角。
根据本发明一种优选的实施方式,所述舵机指令解算模块3用于实时根据获取的飞行器与目标的位置信息及飞行器姿态信息获得打舵指令。
其中,所述舵机指令解算模块包括实时获得飞行器的需用过载,并结合飞行器的实时姿态信息获得打舵指令。
优选地,所述舵机指令解算模块为微型计算机。
在进一步优选的实施方式中,所述舵机指令解算模块3采用比例导引制导律进行制导控制,即
Figure BDA0002209658140000111
其中,导航比N取值为4,aM为飞行器的需用过载,V为飞行器的速度,
Figure BDA0002209658140000112
为飞行器与目标的视线角速度,即为视线角速度解算模块7实时输出的飞行器与目标的视线角速度。
其中,飞行器的速度通过现有技术中常用方法测得。
根据本发明一种优选的实施方式,所述执行模块4用于实时根据接收到的舵机指令调整飞行器姿态。
在进一步优选的实施方式中,所述执行模块4为舵机,优选为气动舵机,更优选为冷气式舵机。
其中,所述冷气式舵机由四个舵片、两个气瓶及相应的操纵机构组成。舵机经过小型化设计,减小舵片面积、提前舵片位置、缩短气瓶长度、增加气瓶数量,在充分利用空间的同时保证充足的气源。
在本发明中,优选采用型号为AF09的冷气式舵机。
舵机接收指令信号,根据信号宽度进行不同时长的偏转,控制飞行器以稳定的姿态沿弹道飞行并命中目标。
具体地,舵机分为两对(即具有四个舵机),其中一对根据舵机指令工作,以调节姿态;另一对舵机根据地磁传感器获得的飞行器滚转角工作,提供重力补偿。
提供重力补偿的方法如下:
当滚转角在0°~180°之间时,以地磁纵轴为标准,当地磁纵轴转过
Figure BDA0002209658140000113
时,舵机控制舵片偏转+5°,给出向上的控制力,平均控制力为:
Figure BDA0002209658140000121
方向沿重力向上;
其中,
Figure BDA0002209658140000122
表示飞行器滚转角度,由地磁传感器测得,F表示舵机提供的控制力;
Figure BDA0002209658140000123
表示滚转角在0°~180°时的平均控制力。
当滚转角在180°~360°之间时,以地磁纵轴为标准,当地磁纵轴转过
Figure BDA0002209658140000124
时,舵机控制舵片偏转-5°,给出向上的控制力,平均控制力为:
Figure BDA0002209658140000125
其中,
Figure BDA0002209658140000126
表示飞行器的滚转角度,由地磁传感器测得,
Figure BDA0002209658140000127
表示滚转角在180°~360°时的平均控制力。
Figure BDA0002209658140000128
即可确定γ和γ*的值。
其中,
Figure BDA0002209658140000129
表示滚转角在0°~180°时的平均控制力,
Figure BDA00022096581400001210
表示滚转角在180°~360°时的平均控制力,m表示飞行器质量。
根据本发明一种优选的实施方式,所述高动态旋转飞行器复合制导控制系统还包括电源供电模块8,以向各模块提供所需的电压,保证正常工作。
其中,电源供电模块接入飞行器上装载的热电源,并对整个电路的输入输出进行整合,防止出现短路造成的烧坏系统的问题。电源供电模块可向各模块提供其所需的额定电压,保证元件正常工作;还可以向有特定需求的模块提供复位电压信号。
本发明还提供了一种高动态旋转飞行器的复合制导控制方法,优选采用上述的复合制导控制系统进行,所述方法包括以下步骤:
步骤1,通过制导模块1中的雷达导引系统11实时获得飞行器与目标之间的距离,由制导模式选择模块5选择制导模式。
其中,所述制导模块1包括雷达导引系统11和红外导引系统12。
根据本发明一种优选的实施方式,所述制导模式选择模块5按照以下规则选择制导模式:
当飞行器与目标之间的距离小于等于10Km且大于2.99Km时,选择雷达导引系统进行制导;
当飞行器与目标之间的小于等于2.99Km且大于0.99Km时,选择雷达导引系统和红外导引系统进行复合制导;
当飞行器与目标之间的距离小于等于0.99km时,选择雷达导引系统和红外导引系统进行复合制导,其中,飞行器与目标的视线角信息由红外导引系统获得,飞行器与目标之间的距离由雷达导引系统获得。
其中,当飞行器与目标之间的距离小于等于2.99Km且大于0.99Km时,选择雷达导引系统和红外导引系统进行复合制导,此时的飞行器与目标的视线角由两种制导模式下得到的视线角加权得到,飞行器与目标之间的距离由雷达导引系统获得。
在更进一步优选的实施方式中,当飞行器与目标之间的距离小于等于2.99Km且大于0.99Km时,复合制导模式下的视线角通过下式(一)获得:
Figure BDA0002209658140000131
其中,
Figure BDA0002209658140000132
表示复合制导模式下的飞行器与目标的视线角,q1表示雷达导引系统进行制导模式下的飞行器与目标的视线角,q2表示红外导引系统进行制导模式下的飞行器与目标的视线角,a为复合制导模式下雷达导引系统所占权重,b为复合制导模式下红外导引系统所占权重。
优选地,所述a和b分别由下式(二)和式(三)获得:
a=0.5x-0.5(二);
b=1-a(三)
其中,x表示飞行器与目标之间的距离。
步骤2,通过导引头6实时获得飞行器与目标的视线角,并由视线角速度解算模块7实时解算得到飞行器与目标的视线角速度,并传递至舵机指令解算模块3。
其中,所述导引头6为捷联导引头。
根据本发明一种优选的实施方式,所述视线角速度解算模块7通过下式(四)、(五)和(六)实时获得飞行器与目标的视线角速度信息:
Figure BDA0002209658140000141
Figure BDA0002209658140000142
Figure BDA0002209658140000143
其中,qg表示由导引头获得的飞行器与目标的视线角;q0表示飞行器与目标的视线角的估测值,即在解算过程中通过上述式(四)、(五)、(六)估算出来的飞行器与目标的视线角的估测值;q01表示飞行器与目标的视线角速度估测值,即在解算过程中通过上述式(四)、(五)、(六)估算出来的飞行器与目标的视线角速度的估测值;
Figure BDA0002209658140000144
表示x2的导数;
Figure BDA0002209658140000145
表示x1的导数;
Figure BDA0002209658140000146
表示x0的导数,上一时刻得到的数值作为下一时刻迭代的初始值;
初始时刻时,x0=0,x1=0,x2=0,每隔0.001s作为积分步长,进行迭代,并且得到x0、x1和x2在下一时刻的数值;
具体来说,第一个迭代时刻中,将初始时刻x0=0,x1=0,x2=0和接收到的qg值代入到式(四)、(五)、(六)中,从而解算出
Figure BDA0002209658140000147
进而得到下一时刻的初始值x0、x1和x2;再将得到的x0、x1、x2和接收到的qg值代入到式(四)、(五)、(六)中,从而得到再下一时刻对应的
Figure BDA0002209658140000151
如此持续循环迭代即可持续得到每次积分得到的对应x0、x1、x2
其中,
Figure BDA0002209658140000152
表示飞行器与目标的视线角速度,实时输出该飞行器与目标的视线角速度给舵机指令解算模块3即可用于解算需用过载。
其中,所述a0、a1、a2、δ、k1和k2均为设计参数,本申请中优选地,a0=1~1.5、a1=7~10、a2=10~15、δ=1~2、k1=0.1~0.4和k2=0.2~0.4;
更优选地,所述a0=1.1、a1=8.5、a2=11.5、δ=1.5、k1=0.3、k2=0.3。
步骤3,通过飞行器姿态测量模块2实时获得飞行器的飞行参数信息,并传递至舵机指令解算模块3。
其中,所述飞行器姿态测量模块2包括地磁传感器21和惯性陀螺22,所述地磁传感器21用于实时获得飞行器的滚转角信息,
所述惯性陀螺22用于实时获得飞行器的俯仰角信息和偏航角信息。
本发明人研究发现,地磁传感器输出地磁信号在载体坐标系下的三轴分量,结合载体坐标系与导航坐标系的转换关系、惯性陀螺测得的飞行器俯仰角及偏航角信息,可以获得飞行器精确的滚转角信息,即输出滚转角信息。
根据本发明一种优选的实施方式,所述输出滚转角信息通过下式(七)获得:
Figure BDA0002209658140000153
其中,Hx*表示地磁信号在载体坐标系的x轴分量;Hy*表示地磁信号在载体坐标系的y轴分量;Hz*表示地磁信号在载体坐标系的z轴分量;Hx表示地磁信号在导航坐标系中的x轴分量;Hy表示地磁信号在导航坐标系中的y轴分量;Hz表示地磁信号在导航坐标系中的z轴分量;其中,Hx*、Hy*、Hz*、Hx、Hy、Hz由地磁传感器直接测得;ψ表示偏航角,θ表示俯仰角,均由惯性陀螺直接测得;γ表示滚转角。
其中,载体坐标系与导航坐标系的转换关系如下式(八)所示:
Figure BDA0002209658140000161
其中,Hx*表示地磁信号在载体坐标系的x轴分量;Hy*表示地磁信号在载体坐标系的y轴分量;Hz*表示地磁信号在载体坐标系的z轴分量;Hx表示地磁信号在导航坐标系中的x轴分量;Hy表示地磁信号在导航坐标系中的y轴分量;Hz表示地磁信号在导航坐标系中的z轴分量;ψ表示偏航角;θ表示俯仰角;γ表示滚转角。
步骤4,通过舵机指令解算模块3实时获得打舵指令,并传递至执行模块4。
其中,所述舵机指令解算模块3采用比例导引制导律进行制导控制,即aM=NVq;其中,导航比N取值为4,aM为飞行器的需用过载,V为飞行器的速度,
Figure BDA0002209658140000162
为飞行器与目标的视线角速度,即为视线角速度解算模块7实时输出的飞行器与目标的视线角速度。
其中,飞行器的速度通过现有技术中常用方法测得。
步骤5,通过执行模块4实时调整飞行器的姿态。
其中,所述执行模块4为舵机,优选为气动舵机,更优选为冷气式舵机。
舵机接收指令信号,根据信号宽度进行不同时长的偏转,控制飞行器以稳定的姿态沿弹道飞行并命中目标。
在发明中所述的高动态旋转飞行器的复合制导控制方法,可以实现全天候、多环境下的精确制导,不易受到内外部干扰。
实施例
实施例1
通过计算机进行飞行器的模拟仿真实验,高动态旋转飞行器的模拟条件为:旋转飞行器的飞行速度为200m/s,转速为4πrad/s,飞行器与目标的距离为10Km;
具体仿真的制导控制方法为:
(1)根据飞行器与目标的距离,选择单模制导模式(雷达导引系统);
(2)通过计算机能够直接模拟该旋转飞行器,并且可以实时给出该旋转飞行器对应的飞行器与目标的视线角、飞行器与目标的视线角速度、飞行器的俯仰角、偏航角,其中,给出的飞行器与目标的视线角速度为该飞行器的导引头获得的真实值,将其作为真实值实时输送至该旋转飞行器;
(3)通过下式(四)、(五)和(六)实时解算得到飞行器与目标的视线角速度信息的解算值:
Figure BDA0002209658140000171
Figure BDA0002209658140000172
Figure BDA0002209658140000173
其中,qg表示由导引头获得的飞行器与目标的视线角;q0表示飞行器与目标的视线角的估测值,即在解算过程中通过上述式(四)、(五)、(六)估算出来的飞行器与目标的视线角的估测值;q01表示飞行器与目标的视线角速度估测值,即在解算过程中通过上述式(四)、(五)、(六)估算出来的飞行器与目标的视线角速度的估测值;
Figure BDA0002209658140000181
表示x2的导数;
Figure BDA0002209658140000182
表示x1的导数;
Figure BDA0002209658140000183
表示x0的导数,上一时刻得到的数值作为下一时刻迭代的初始值;
初始时刻时,x0=0,x1=0,x2=0,每隔0.001s作为积分步长,进行迭代,并且得到x0、x1和x2在下一时刻的数值;
具体来说,第一个迭代时刻中,将初始时刻x0=0,x1=0,x2=0和接收到的qg值代入到式(四)、(五)、(六)中,从而解算出
Figure BDA0002209658140000184
进而得到下一时刻的初始值x0、x1和x2;再将得到的x0、x1、x2和接收到的qg值代入到式(四)、(五)、(六)中,从而得到再下一时刻对应的
Figure BDA0002209658140000185
如此持续循环迭代即可持续得到每次积分得到的对应x0、x1、x2
其中,
Figure BDA0002209658140000186
表示飞行器与目标的视线角速度,即为所述飞行器与目标的视线角速度的解算值;
其中,所述a0=1.1、a1=8.5、a2=11.5、δ=1.5、k1=0.3、k2=0.3;
(4)通过下式(七)获得滚转角的解算值:
Figure BDA0002209658140000187
其中,Hx*表示地磁信号在载体坐标系的x轴分量;Hy*表示地磁信号在载体坐标系的y轴分量;Hz*表示地磁信号在载体坐标系的z轴分量;Hx表示地磁信号在导航坐标系中的x轴分量;Hy表示地磁信号在导航坐标系中的y轴分量;Hz表示地磁信号在导航坐标系中的z轴分量;ψ表示偏航角,θ表示俯仰角;γ表示滚转角;
其中,载体坐标系与导航坐标系的转换关系如下式(八)所示:
Figure BDA0002209658140000191
其中,Hx*表示地磁信号在载体坐标系的x轴分量;Hy*表示地磁信号在载体坐标系的y轴分量;Hz*表示地磁信号在载体坐标系的z轴分量;Hx表示地磁信号在导航坐标系中的x轴分量;Hy表示地磁信号在导航坐标系中的y轴分量;Hz表示地磁信号在导航坐标系中的z轴分量;ψ表示偏航角;θ表示俯仰角;γ表示滚转角。
在不考虑内外部扰动的影响下,计算机实时给出的飞行器与目标的视线角速度的真实值与解算值(估测值)之间的比较结果,如图2所示。
由图2可知,本发明中飞行器与目标的视线角速度的解算值(估测值)能够快速逼近真实视线角速度,所需时间不到0.4秒,且在追踪上真实视线角速度后能够随真实视线角速度变化而变化,表明本申请中的视线角速度解算模块的精确性较高。
在考虑内外部扰动的情况下,扰动为幅值为0.3度视线角正弦扰动,频率为0.3Hz,即扰动项为0.3sin(2t),计算机实时给出的飞行器与目标的视线角速度的真实值与解算值(估测值)之间的比较结果,如图3所示。
由图3可知,本发明中飞行器与目标的视线角速度的解算值(估测值)依然能够快速逼近真实视线角速度,所需时间不到0.4秒,且在追踪上真实视线角速度后能够随真实视线角速度变化而变化,表明本申请中的视线角速度解算模块的精确性较高、鲁棒性较强。
实施例2
本实施例与实施例1相似,区别在于:
飞行器与目标的距离为2.99Km;
具体仿真的制导控制方法的步骤(1)中:
根据飞行器与目标的距离,选择双模制导模式(雷达导引系统和红外导引系统进行复合制导),所述复合制导的视线角通过下式(一)获得:
Figure BDA0002209658140000201
其中,
Figure BDA0002209658140000202
表示复合制导模式下的飞行器与目标的视线角,q1表示雷达导引系统进行制导模式下的飞行器与目标的视线角,q2表示红外导引系统进行制导模式下的飞行器与目标的视线角,a为复合制导模式下雷达导引系统所占权重,b为复合制导模式下红外导引系统所占权重;
其中,所述a和b分别由下式(二)和式(三)获得:
a=0.5x-0.5(二);
b=1-a(三)
其中,x表示飞行器与目标之间的距离;
在不考虑内外部扰动的影响下,计算机实时给出的飞行器与目标的视线角速度的真实值与解算值(估测值)之间的比较结果,如图4所示。
由图4可知,本发明中飞行器与目标的视线角速度的解算值(估测值)能够快速逼近真实视线角速度,所需时间不到0.5秒,且在追踪上真实视线角速度后能够随真实视线角速度变化而变化,表明本申请中的视线角速度解算模块的精确性较高。
在考虑内外部扰动的情况下,扰动为幅值为0.5度视线角正弦扰动,频率为0.8Hz,即扰动项为0.5sin(5t),计算机实时给出的飞行器与目标的视线角速度的真实值与解算值之间的比较结果,如图5所示。
由图5可知,本发明中飞行器与目标的视线角速度的解算值依然能够快速逼近真实视线角速度,所需时间不到0.6秒,且在追踪上真实视线角速度后能够随真实视线角速度变化而变化,表明本申请中的视线角速度解算模块的精确性较高、鲁棒性较强。
实施例3
本实施例与实施例1相似,区别在于:
飞行器与目标的距离为0.99Km;
具体仿真的制导控制方法的步骤(1)中:
根据飞行器与目标的距离,选择双模制导模式(雷达导引系统和红外导引系统进行复合制导),其中,飞行器与目标的视线角信息由红外导引系统获得,飞行器与目标之间的距离由雷达导引系统获得。
在不考虑内外部扰动的影响下,计算机实时给出的飞行器与目标的视线角速度的真实值与解算值(估测值)之间的比较结果,如图6所示。
由图6可知,本发明中飞行器与目标的视线角速度的解算值能够快速逼近真实视线角速度,所需时间不到0.4秒,且在追踪上真实视线角速度后能够随真实视线角速度变化而变化,表明本申请中的视线角速度解算模块的精确性较高。
在考虑内外部扰动的情况下,扰动为幅值为1度视线角正弦扰动,频率为1.3Hz,即扰动项为sin(8t),计算机实时给出的飞行器与目标的视线角速度的真实值(估测值)与解算值之间的比较结果,如图7所示。
由图7可知,本发明中飞行器与目标的视线角速度的解算值依然能够快速逼近真实视线角速度,所需时间不到0.4秒,且在追踪上真实视线角速度后能够随真实视线角速度变化而变化,表明本申请中的视线角速度解算模块的精确性较高、鲁棒性较强。
以上结合了优选的实施方式对本发明进行了说明,不过这些实施方式仅是范例性的,仅起到说明性的作用。在此基础上,可以对本发明进行多种替换和改进,这些均落入本发明的保护范围内。

Claims (8)

1.一种高动态旋转飞行器复合制导控制系统,其特征在于,所述系统包括制导模块(1)、飞行器姿态测量模块(2)、舵机指令解算模块(3)和执行模块(4);其中,
所述制导模块(1)用于实时获取飞行器与目标的位置信息,并传递至舵机指令解算模块(3),
所述飞行器姿态测量模块(2)用于实时获取飞行器的飞行参数信息,并传递至舵机指令解算模块(3),
所述舵机指令解算模块(3)用于实时根据获取的飞行器与目标的位置信息及飞行器姿态信息获得打舵指令,
所述执行模块(4)用于实时根据接收到的舵机指令调整飞行器姿态;
所述制导模块(1)采用复合制导模式,包括雷达导引系统(11)和红外导引系统(12);
所述系统还包括制导模式选择模块(5),其根据飞行器与目标之间的距离选择制导模式;
在所述制导模式选择模块(5)与舵机指令解算模块(3)之间还设置有导引头(6),所述导引头(6)为捷联导引头;
在所述导引头(6)与舵机指令解算模块(3)之间还设置有视线角速度解算模块(7),以根据接收到的视线角信息实时解算得到飞行器与目标的视线角速度;
所述视线角速度解算模块(7)通过下式(四)、(五)和(六)实时获得飞行器与目标的视线角速度信息:
Figure FDA0003395172070000011
Figure FDA0003395172070000012
Figure FDA0003395172070000013
其中,qg表示由导引头获得的飞行器与目标的视线角;q0表示飞行器与目标的视线角的估测值,即在解算过程中通过上述式(四)、(五)、(六)估算出来的飞行器与目标的视线角的估测值;q01表示飞行器与目标的视线角速度估测值,即在解算过程中通过上述式(四)、(五)、(六)估算出来的飞行器与目标的视线角速度的估测值;
Figure FDA0003395172070000023
表示x2的导数;
Figure FDA0003395172070000024
表示x1的导数;
Figure FDA0003395172070000025
表示x0的导数,上一时刻得到的x0、x1和x2数值作为下一时刻迭代的初始值;
所述a0、a1、a2、δ、k1和k2均为设计参数,a0=1~1.5、a1=7~10、a2=10~15、δ=1~2、k1=0.1~0.4和k2=0.2~0.4。
2.根据权利要求1所述的系统,其特征在于,所述制导模式选择模块(5)按照以下规则选择制导模式:
当飞行器与目标之间的距离小于等于10Km且大于2.99Km时,选择雷达导引系统进行制导;
当飞行器与目标之间的距离小于等于2.99km且大于0.99Km时,选择雷达导引系统和红外导引系统进行复合制导;
当飞行器与目标之间的距离小于等于0.99km时,选择雷达导引系统和红外导引系统进行复合制导,其中,飞行器与目标的视线角信息由红外导引系统获得,飞行器与目标之间的距离由雷达导引系统获得。
3.根据权利要求2所述的系统,其特征在于,当飞行器与目标之间的距离小于等于2.99Km且大于0.99Km时,复合制导模式下的视线角通过下式(一)获得:
Figure FDA0003395172070000021
其中,
Figure FDA0003395172070000022
表示复合制导模式下的飞行器与目标的视线角,q1表示雷达导引系统进行制导模式下的飞行器与目标的视线角,q2表示红外导引系统进行制导模式下的飞行器与目标的视线角,a为复合制导模式下雷达导引系统所占权重,b为复合制导模式下红外导引系统所占权重。
4.根据权利要求3所述的系统,其特征在于,式(一)中的a和b分别由下式(二)和式(三)获得:
a=0.5x-0.5 (二);
b=1-a (三)
其中,x表示飞行器与目标之间的距离。
5.根据权利要求1所述的系统,其特征在于,所述飞行器姿态测量模块(2)包括地磁传感器(21)和惯性陀螺(22),其中,
所述地磁传感器(21)用于实时获得飞行器的滚转角信息,
所述惯性陀螺(22)用于实时获得飞行器的俯仰角信息和偏航角信息。
6.根据权利要求1所述的系统,其特征在于,所述高动态旋转飞行器复合制导控制系统还包括电源供电模块(8),以向各模块提供所需的电压,保证正常工作。
7.一种高动态旋转飞行器复合制导控制方法,采用权利要求1至6之一所述的系统进行,其特征在于,所述方法包括以下步骤:
步骤1,通过制导模块(1)中的雷达导引系统(11)实时获得飞行器与目标之间的距离,由制导模式选择模块(5)选择制导模式;
步骤2,通过导引头(6)实时获得飞行器与目标的视线角,并由视线角速度解算模块(7)实时解算得到飞行器与目标的视线角速度,并传递至舵机指令解算模块(3);
步骤3,通过飞行器姿态测量模块(2)实时获得飞行器的飞行参数信息,并传递至舵机指令解算模块(3);
步骤4,通过舵机指令解算模块(3)实时获得打舵指令,并传递至执行模块(4);
步骤5,通过执行模块(4)实时调整飞行器的姿态。
8.根据权利要求7所述的方法,其特征在于,步骤2中,所述飞行器与目标的视线角速度信息通过下式(四)、(五)和(六)实时获得:
Figure FDA0003395172070000041
Figure FDA0003395172070000042
Figure FDA0003395172070000043
其中,qg表示由导引头获得的飞行器与目标的视线角;q0表示飞行器与目标的视线角的估测值,即在解算过程中通过上述式(四)、(五)、(六)估算出来的飞行器与目标的视线角的估测值;q01表示飞行器与目标的视线角速度估测值,即在解算过程中通过上述式(四)、(五)、(六)估算出来的飞行器与目标的视线角速度的估测值;
Figure FDA0003395172070000044
表示x2的导数;
Figure FDA0003395172070000045
表示x1的导数;
Figure FDA0003395172070000046
表示x0的导数,上一时刻得到的x0、x1和x2数值作为下一时刻迭代的初始值;
所述a0、a1、a2、δ、k1和k2均为设计参数,a0=1~1.5、a1=7~10、a2=10~15、δ=1~2、k1=0.1~0.4和k2=0.2~0.4。
CN201910894011.1A 2019-09-20 2019-09-20 高动态旋转飞行器复合制导控制系统及方法 Active CN112540618B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910894011.1A CN112540618B (zh) 2019-09-20 2019-09-20 高动态旋转飞行器复合制导控制系统及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910894011.1A CN112540618B (zh) 2019-09-20 2019-09-20 高动态旋转飞行器复合制导控制系统及方法

Publications (2)

Publication Number Publication Date
CN112540618A CN112540618A (zh) 2021-03-23
CN112540618B true CN112540618B (zh) 2022-02-18

Family

ID=75012569

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910894011.1A Active CN112540618B (zh) 2019-09-20 2019-09-20 高动态旋转飞行器复合制导控制系统及方法

Country Status (1)

Country Link
CN (1) CN112540618B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115493589B (zh) * 2022-11-04 2023-03-28 开拓导航控制技术股份有限公司 视线角解算方法和系统
CN116661495B (zh) * 2023-05-30 2024-02-20 北京理工大学 一种飞行器近射程减速控制方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103245256A (zh) * 2013-04-25 2013-08-14 北京理工大学 一种多导弹协同作战制导律设计方法
CN105259908A (zh) * 2015-11-05 2016-01-20 南京航空航天大学 一种雷达引导无人机自动着舰制导与控制系统及其控制方法
CN105539835A (zh) * 2016-01-18 2016-05-04 成都纵横自动化技术有限公司 一种复合翼垂直起降飞行器
CN106406344A (zh) * 2016-09-28 2017-02-15 北京理工大学 一种旋转制导飞行器的增程系统及方法
CN108762286A (zh) * 2018-05-31 2018-11-06 智飞智能装备科技东台有限公司 一种可对多架无人机飞行控制的地面控制系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8380473B2 (en) * 2009-06-13 2013-02-19 Eric T. Falangas Method of modeling dynamic characteristics of a flight vehicle

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103245256A (zh) * 2013-04-25 2013-08-14 北京理工大学 一种多导弹协同作战制导律设计方法
CN105259908A (zh) * 2015-11-05 2016-01-20 南京航空航天大学 一种雷达引导无人机自动着舰制导与控制系统及其控制方法
CN105539835A (zh) * 2016-01-18 2016-05-04 成都纵横自动化技术有限公司 一种复合翼垂直起降飞行器
CN106406344A (zh) * 2016-09-28 2017-02-15 北京理工大学 一种旋转制导飞行器的增程系统及方法
CN108762286A (zh) * 2018-05-31 2018-11-06 智飞智能装备科技东台有限公司 一种可对多架无人机飞行控制的地面控制系统

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Error analysis and correction for hardware-in-the-loop simulation system of laser guidance weapon;Fan Shipeng等;《Infrared and Laser Engineering 》;20130425;904-908 *
Frame angle design of image guided aircraft;Lin Defu等;《Science & Technology Review 》;20150413;55-60 *
临近空间飞行器再入段复合控制律设计;潘彦鹏等;《系统工程与电子技术》;20130822(第11期);全文 *
基于测角信号注入的复合制导半实物仿真方法;赖鹏等;《系统仿真学报》;20150108(第01期);全文 *
小型无人机导航与制导关键技术;任涵;《科技创新与应用》;20180228(第06期);全文 *

Also Published As

Publication number Publication date
CN112540618A (zh) 2021-03-23

Similar Documents

Publication Publication Date Title
CN102506892B (zh) 一种光流多传感器和惯导器件信息融合配置方法
CN108803639A (zh) 一种基于反步法的四旋翼飞行器飞行控制方法
CN103587708A (zh) 超小型无人旋翼飞行器野外定点零盲区自主软着陆方法
Jann Advanced features for autonomous parafoil guidance, navigation and control
CN112540618B (zh) 高动态旋转飞行器复合制导控制系统及方法
CN110470304A (zh) 一种基于无人机光电平台的高精度目标定位测速方法
US7616130B2 (en) Method and a device for processing and displaying aircraft piloting information
CN110068336A (zh) 一种基于磁悬浮控制敏感陀螺平行构型的角运动测量方法
CN111504256A (zh) 一种基于最小二乘法的滚转角实时估计方法
JP2022520118A (ja) 無人航空機の高度から風向および風速度の測定値を決定するためのシステムおよび方法
CN113359856B (zh) 一种无人机指定航向目标点导引方法及系统
Ramos et al. Autonomous flight experiment with a robotic unmanned airship
CN111412793A (zh) 应用于远程制导飞行器上的防侧偏的全射程覆盖控制系统
US4541591A (en) Guidance law to improve the accuracy of tactical missiles
CN109445283B (zh) 一种用于欠驱动浮空器在平面上定点跟踪的控制方法
Du et al. A low-cost attitude estimation system for UAV application
CN115727842B (zh) 一种无人机快速对准方法、系统、计算机设备及存储介质
Yu et al. Design of multi-sensor fusion positioning and navigation method based on federated Kalman filter for low-speed rotorcraft of atmospheric detection
Oshman et al. Mini-UAV altitude estimation using an inertially stabilized payload
Song et al. A high precision autonomous navigation algorithm of UAV based on MEMS sensor
CN110017809B (zh) 利用地磁信息和光流传感器解算飞行器姿态的方法
Peng et al. Real-time attitude estimation for high-speed UAV in high-frequency environmental dithering based on AMCF
CA2135362A1 (en) Method for determining the line-of-sight rates of turn with a rigid seeker head
CN110017831B (zh) 利用地磁信息和声呐传感器解算飞行器姿态的方法
CN111273682B (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