CN106021679A - 一种末端机动小过载制导指令的奇异处理方法 - Google Patents

一种末端机动小过载制导指令的奇异处理方法 Download PDF

Info

Publication number
CN106021679A
CN106021679A CN201610319674.7A CN201610319674A CN106021679A CN 106021679 A CN106021679 A CN 106021679A CN 201610319674 A CN201610319674 A CN 201610319674A CN 106021679 A CN106021679 A CN 106021679A
Authority
CN
China
Prior art keywords
guidance
power
overload
ycx
force
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
Application number
CN201610319674.7A
Other languages
English (en)
Other versions
CN106021679B (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 Aerospace Automatic Control Research Institute
Original Assignee
Beijing Aerospace Automatic Control Research Institute
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 Aerospace Automatic Control Research Institute filed Critical Beijing Aerospace Automatic Control Research Institute
Priority to CN201610319674.7A priority Critical patent/CN106021679B/zh
Publication of CN106021679A publication Critical patent/CN106021679A/zh
Application granted granted Critical
Publication of CN106021679B publication Critical patent/CN106021679B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)

Abstract

本发明公开了一种末端机动小过载制导指令的奇异处理方法,该方法包括:根据预设的末端机动段制导律和方案,计算得到纵向制导力Fycx和横向制导力Fzcx;对制导力进行奇异处理:如果Fycx<Kc·m·g0,则Fyn=sign(Kc·m·g0)·|Kc·m·g0|,否则,Fyn=Fycx;根据奇异处理后的制导力计算制导所需的倾侧角;根据纵向制导力Fycx和横向制导力Fzcx计算制导需用的法向力。通过使用本发明所提供的方法,可以解决高升阻比飞行器在末端机动飞行段由于小过载制导指令导致的奇异问题,降低不必要的大幅弹体滚转,提高制导精度。

Description

一种末端机动小过载制导指令的奇异处理方法
技术领域
本发明涉及高超声速飞行器制导技术领域,特别涉及一种末端机动小过载制导指令的奇异处理方法。
背景技术
升力式高超声速飞行器由于升力大,当飞行器接近目标时,该飞行器在末端机动段飞行时极易出现小过载制导指令的需求,此时制导律给出的纵向和横向制导力需求均比较小,可能会出现纵向制导力|Fycx|≤0.1g0,横向制导力|Fzcx|≤0.1g0的情况,此时如果还是按照γcx=tg-1(Fzcx/Fycx)计算倾侧角指令γcx,将会导致倾侧角大幅滚转,甚至到90°的情况,从而对姿态控制和制导尤为不利。
小过载制导指令引起的倾侧实际上是不必要的机动,其产生的侧向加速度并不大,虽然倾侧大,但弹体法向过载小,其对机动飞行的贡献基本可以忽略,不必要的滚转反而对控制系统造成不利,甚至影响制导精度。
因此,制导系统设计需要考虑消除由于小过载制导指令导致的奇异问题,提高制导的有效性。
发明内容
有鉴于此,本发明提供一种末端机动小过载制导指令的奇异处理方法,从而可以解决高升阻比飞行器在末端机动飞行段由于小过载制导指令导致的奇异问题,降低不必要的大幅弹体滚转,提高制导精度。
本发明的技术方案具体是这样实现的:
一种末端机动小过载制导指令的奇异处理方法,该方法包括:
根据预设的末端机动段制导律和方案,计算得到纵向制导力Fycx和横向制导力Fzcx
对制导力进行奇异处理:如果Fycx<Kc·m·g0,则Fyn=sign(Kc·m·g0)·|Kc·m·g0|,否则,Fyn=Fycx;其中,Fyn为限幅后的纵向制导力;Kc为设计参数;m为飞行质量;g0为地球引力加速度;
根据奇异处理后的制导力计算制导所需的倾侧角;
根据纵向制导力Fycx和横向制导力Fzcx计算制导需用的法向力。
较佳的,所述Kc的取值大于1.0。
较佳的,使用如下的公式计算得到制导所需的倾侧角
&gamma; c x d = tg - 1 ( F z c x F y n ) , &gamma; c x d &Element; ( - &pi; , &pi; &rsqb; .
较佳的,使用如下的公式计算得到制导需用的法向力Fn
F n = s i g n ( F y c x ) &CenterDot; F y c x 2 + F z c x 2 .
如上可见,在本发明中的末端机动小过载制导指令的奇异处理方法中,由于先根据预设的末端机动段制导律和方案,计算得到纵向制导力和横向制导力;然后对制导力进行奇异处理,并根据奇异处理后的制导力计算制导所需的倾侧角,最后可以根据纵向制导力和横向制导力计算制导需用的法向力,因此可以根据计算得到的倾侧角和法向力实现飞行器在滑翔过程中的倾侧翻转,从而可以从制导力需求和倾侧机动的效率入手,当弹体法向为小过载的条件下时(即出现横法向小过载需求时),消除倾侧角指令的奇异性,尽量不进行倾侧机动,降低不必要的大幅弹体滚转,因而可以避免无意义的姿态机动,提高飞行器倾侧机动飞行的有效性,完全可避免高升阻比飞行器在末端机动飞行段由于小过载制导指令导致的奇异问题,降低不必要的大幅弹体滚转,提高制导精度。与现有技术中的采用倾侧角进行直接限幅的方式相比,尤其当飞行器真正需要大倾侧转弯时,采用本发明中的末端机动小过载制导指令的奇异处理方法,将不会因为倾侧角别限制而影响制导精度,也不会因为过载变化而引起指令跳跃。
附图说明
图1为本发明实施例中的末端机动小过载制导指令的奇异处理方法的流程示意图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白,以下参照附图并举实施例,对本发明进一步详细说明。
本实施例提供了一种末端机动小过载制导指令的奇异处理方法,该方法适用于以气动力飞行为主的升力式飞行器。
图1为本发明实施例中的末端机动小过载制导指令的奇异处理方法的流程示意图。如图1所示,本发明实施例中的末端机动小过载制导指令的奇异处理方法包括:
步骤101,根据预设的末端机动段制导律和方案,计算得到纵向制导力Fycx和横向制导力Fzcx
在本发明的技术方案中,在本步骤中,将首先按照预设的末端机动段制导律和方案,进行纵向制导力Fycx和横向制导力Fzcx需求的计算。
另外,在本发明的技术方案中,末端机动制导段可以采用最优比例导引进行制导。因此,在本发明的一个较佳实施例中,可以按照最优比例导引进行制导计算,得到纵向和横向所需要的制导力,即纵向制导力Fycx和横向制导力Fzcx。其中,可以使用本领域中的常用方法按照最优比例导引进行制导计算,因此在此不再赘述。
步骤102,对制导力进行奇异处理:如果Fycx<Kc·m·g0,则Fyn=sign(Kc·m·g0)·|Kc·m·g0|,否则,Fyn=Fycx
其中,Fyn为限幅后的纵向制导力;Kc为设计参数;m为飞行质量,是一个常数;g0为地球引力加速度,较佳的,在本发明的一个具体实施例中,所述g0取值为9.8055。
在本步骤中,将对制导力进行奇异处理,即当弹体法向为小过载的条件下时(即出现横法向小过载需求时,例如,Fycx<Kc·m·g0),则需要消除倾侧角指令的奇异性,尽量不进行倾侧机动,降低不必要的大幅弹体滚转,从而避免无意义的姿态机动,提高飞行器倾侧机动飞行的有效性。
较佳的,在本发明的一个具体实施例中,所述Kc的取值大于1.0。Kc值的选取原则是能够消除倾侧角指令的奇异性,但又不至于影响制导精度。Kc值过大,虽然能够消除倾侧角指令的奇异性,但无法满足飞行器真正需要的倾侧指令,最终导致制导精度差。
步骤103,根据奇异处理后的制导力计算制导所需的倾侧角。
较佳的,在本发明的一个具体实施例中,所述倾侧角可以使用如下所述的公式计算得到:
&gamma; c x d = tg - 1 ( F z c x F y n ) , &gamma; c x d &Element; ( - &pi; , &pi; &rsqb; - - - ( 1 )
其中,为制导所需的倾侧角。
步骤104,根据纵向制导力Fycx和横向制导力Fzcx计算制导需用的法向力。
在本发明的技术方案中,由于飞行器运动所需的力由弹体提供,因此,较佳的,在本发明的一个具体实施例中,所述制导需用的法向力可以使用如下所述的公式计算得到:
F n = s i g n ( F y c x ) &CenterDot; F y c x 2 + F z c x 2 - - - ( 2 )
通过上述的步骤101~104,即可得到制导所需的倾侧角和法向力,从而可以根据计算得到的倾侧角和法向力实现飞行器在滑翔过程中的倾侧翻转。
综上可知,在本发明中的末端机动小过载制导指令的奇异处理方法中,由于先根据预设的末端机动段制导律和方案,计算得到纵向制导力和横向制导力;然后对制导力进行奇异处理,并根据奇异处理后的制导力计算制导所需的倾侧角,最后可以根据纵向制导力和横向制导力计算制导需用的法向力,因此可以根据计算得到的倾侧角和法向力实现飞行器在滑翔过程中的倾侧翻转,从而可以从制导力需求和倾侧机动的效率入手,当弹体法向为小过载的条件下时(即出现横法向小过载需求时),消除倾侧角指令的奇异性,尽量不进行倾侧机动,降低不必要的大幅弹体滚转,因而可以避免无意义的姿态机动,提高飞行器倾侧机动飞行的有效性,完全可避免高升阻比飞行器在末端机动飞行段由于小过载制导指令导致的奇异问题,降低不必要的大幅弹体滚转,提高制导精度。与现有技术中的采用倾侧角进行直接限幅的方式相比,尤其当飞行器真正需要大倾侧转弯时,采用本发明中的末端机动小过载制导指令的奇异处理方法,将不会因为倾侧角别限制而影响制导精度,也不会因为过载变化而引起指令跳跃。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明保护的范围之内。

Claims (4)

1.一种末端机动小过载制导指令的奇异处理方法,其特征在于,该方法包括:
根据预设的末端机动段制导律和方案,计算得到纵向制导力Fycx和横向制导力Fzcx
对制导力进行奇异处理:如果Fycx<Kc·m·g0,则Fyn=sign(Kc·m·g0)·|Kc·m·g0|,否则,Fyn=Fycx;其中,Fyn为限幅后的纵向制导力;Kc为设计参数;m为飞行质量;g0为地球引力加速度;
根据奇异处理后的制导力计算制导所需的倾侧角;
根据纵向制导力Fycx和横向制导力Fzcx计算制导需用的法向力。
2.根据权利要求1所述的方法,其特征在于:
所述Kc的取值大于1.0。
3.根据权利要求2所述的方法,其特征在于,使用如下的公式计算得到制导所需的倾侧角
&gamma; c x d = tg - 1 ( F z c x F y n ) , &gamma; c x d &Element; ( - &pi; , &pi; &rsqb; .
4.根据权利要求3所述的方法,其特征在于,使用如下的公式计算得到制导需用的法向力Fn
F n = s i g n ( F y c x ) &CenterDot; F y c x 2 + F z c x 2 .
CN201610319674.7A 2016-05-13 2016-05-13 一种末端机动小过载制导指令的奇异处理方法 Expired - Fee Related CN106021679B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610319674.7A CN106021679B (zh) 2016-05-13 2016-05-13 一种末端机动小过载制导指令的奇异处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610319674.7A CN106021679B (zh) 2016-05-13 2016-05-13 一种末端机动小过载制导指令的奇异处理方法

Publications (2)

Publication Number Publication Date
CN106021679A true CN106021679A (zh) 2016-10-12
CN106021679B CN106021679B (zh) 2017-04-19

Family

ID=57100806

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610319674.7A Expired - Fee Related CN106021679B (zh) 2016-05-13 2016-05-13 一种末端机动小过载制导指令的奇异处理方法

Country Status (1)

Country Link
CN (1) CN106021679B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108801081A (zh) * 2018-06-15 2018-11-13 上海航天控制技术研究所 一种针对btt导弹的非奇异滚转指令生成算法
CN110209197A (zh) * 2019-06-25 2019-09-06 湖北航天技术研究院总体设计所 一种飞行器控制系统设计方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060022091A1 (en) * 2004-07-30 2006-02-02 Peck Mason A Method and system for optimizing torque in a cmg array
CN104266546A (zh) * 2014-09-22 2015-01-07 哈尔滨工业大学 一种基于视线的有限时间收敛主动防御制导控制方法
CN105022403A (zh) * 2015-05-11 2015-11-04 北京航天自动控制研究所 滑翔飞行器的纵向轨迹控制增益的确定方法
CN105466660A (zh) * 2015-12-29 2016-04-06 清华大学 一种采用加速度观测器的飞机气动参数分区分步辨识方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060022091A1 (en) * 2004-07-30 2006-02-02 Peck Mason A Method and system for optimizing torque in a cmg array
CN104266546A (zh) * 2014-09-22 2015-01-07 哈尔滨工业大学 一种基于视线的有限时间收敛主动防御制导控制方法
CN105022403A (zh) * 2015-05-11 2015-11-04 北京航天自动控制研究所 滑翔飞行器的纵向轨迹控制增益的确定方法
CN105466660A (zh) * 2015-12-29 2016-04-06 清华大学 一种采用加速度观测器的飞机气动参数分区分步辨识方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
探月飞船预测-校正再入制导率设计;李慧峰等;《空间控制技术与应用》;20090228;第35卷(第1期);第19-24页 *
李慧峰等: "探月飞船预测-校正再入制导率设计", 《空间控制技术与应用》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108801081A (zh) * 2018-06-15 2018-11-13 上海航天控制技术研究所 一种针对btt导弹的非奇异滚转指令生成算法
CN110209197A (zh) * 2019-06-25 2019-09-06 湖北航天技术研究院总体设计所 一种飞行器控制系统设计方法

Also Published As

Publication number Publication date
CN106021679B (zh) 2017-04-19

Similar Documents

Publication Publication Date Title
CN104615144B (zh) 基于目标规划的高超声速飞行器再入轨迹在线优化方法
US11286065B2 (en) Method for designing reentry trajectory based on flight path angle planning
US10156846B2 (en) Lane changing apparatus and method of autonomous vehicle
US9207681B2 (en) Automatic recovery method for an unmanned aerial vehicle
CN106371312B (zh) 基于模糊控制器的升力式再入预测-校正制导方法
CN104648695B (zh) 一种基于倾侧角可用性的再入走廊最优规划方法
CN105159308B (zh) 一种可重复使用飞行器着陆段制导与控制律一体化耦合设计方法
CN104326093B (zh) 光学成像小卫星姿态控制系统及工作模式在轨切换方法
CN106643341B (zh) 基于准平衡滑翔原理的力热控制耦合设计方法
CN103625637B (zh) 一种大型飞机侧向阵风缓和方法
CN109018440A (zh) 垂直起降运载火箭子级着陆段精确垂直软着陆四次多项式制导方法
CN103587681A (zh) 抑制侧滑角信号常值偏差影响的高超声速飞行器控制方法
CN106483974A (zh) 一种固定翼无人机近距离几何避障方法
CN110733670B (zh) 一种短航程低过载的再入轨迹设计方法
CN106021679B (zh) 一种末端机动小过载制导指令的奇异处理方法
CN105947239B (zh) 一种基于矢量投影的航天器光学敏感器视场分析方法
CN104656659B (zh) 舰载机滑跃起飞自动飞行控制方法
CN104085539A (zh) 成像定标的姿态控制方法
CN105138808B (zh) 基于摄动理论的滑翔弹道误差传播分析方法
CN111077903B (zh) 一种基于窗口决策的高容错性的天钩回收方法
CN109101725A (zh) 航天器受控再入落区预示方法
US20100019078A1 (en) Missile guidance system
CN105759830B (zh) 一种升力式飞行器高动态下压段制导方法
CN104656661A (zh) 一种公务机下降着陆控制方法
CN102880049A (zh) 一种基于帆板挠性形变测量的自适应振动控制方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170419

Termination date: 20200513

CF01 Termination of patent right due to non-payment of annual fee