CN104731104B - 一种高超声速飞行器滑翔飞行段的纵向制导方法 - Google Patents
一种高超声速飞行器滑翔飞行段的纵向制导方法 Download PDFInfo
- Publication number
- CN104731104B CN104731104B CN201510102948.2A CN201510102948A CN104731104B CN 104731104 B CN104731104 B CN 104731104B CN 201510102948 A CN201510102948 A CN 201510102948A CN 104731104 B CN104731104 B CN 104731104B
- Authority
- CN
- China
- Prior art keywords
- attack
- angle
- instruction
- path instructions
- additional angle
- 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.)
- Expired - Fee Related
Links
Landscapes
- Traffic Control Systems (AREA)
Abstract
本发明公开了一种高超声速飞行器滑翔飞行段的纵向制导方法。该方法包括:在滑翔段飞行过程中,制导系统根据导航参数实时生成标准飞行轨迹指令;制导系统根据标准飞行轨迹指令中的飞行速度或马赫给出预置的攻角acx0;根据轨迹指令和导航系统提供的测量值,计算附加攻角指令Dacx;根据预置的攻角acx0和附加攻角指令Dacx,计算得到当前的实际攻角acx。通过使用本发明所提供的高超声速飞行器滑翔飞行段的纵向制导方法,可以效地改善阻力加速度跟踪控制的动态特性,提高动态条件下阻力加速度的跟踪控制精度,抑制飞行轨迹的波动,增强制导系统对各种偏差、不确定干扰和不确定条件的适应能力。
Description
技术领域
本发明涉及高超声速飞行器导航、制导与控制技术,特别涉及一种高超声速飞行器滑翔飞行段的纵向制导方法。
背景技术
在现有技术中,对于具有高升阻比特性的滑翔飞行器,在其高超声速飞行段大多采用预置攻角飞行,对其飞行攻角不作任何调整。当飞行器进行倾侧角翻转时,由于滚动角速率的限制,角速率一般小于20°/秒,因此,如果进行从负60度倾侧角到正60度倾侧角的翻转操作,一般至少需要6秒钟的时间,在偏差条件下甚至将超过10秒,因此将对阻力加速度跟踪制导产生不利的影响,导致动态条件下的阻力加速度跟踪的误差较大。此外,由于飞行器的攻角不作任何调整,因此当飞行器进行倾侧角翻转时,还会导致滑翔飞行器的飞行高度和弹道倾角发生波动,从而对制导产生不利的影响。再次,当飞行过程中存在各种干扰和不确定条件时,飞行器将面临着需用速度过剩或不足的情况,而如果能适当调整攻角则可补偿飞行器的能力不足或增强减速控制效果。
发明内容
有鉴于此,本发明提供一种高超声速飞行器滑翔飞行段的纵向制导方法,从而可以有效地改善阻力加速度跟踪控制的动态特性,提高动态条件下阻力加速度的跟踪控制精度,抑制飞行轨迹的波动,增强制导系统对各种偏差、不确定干扰和不确定条件的适应能力。
本发明的技术方案具体是这样实现的:
一种高超声速飞行器滑翔飞行段的纵向制导方法,该方法包括:
在滑翔段飞行过程中,制导系统根据导航参数实时生成标准飞行轨迹指令;
制导系统根据标准飞行轨迹指令中的飞行速度或马赫给出预置的攻角αcx0;
根据轨迹指令和导航系统提供的测量值,计算附加攻角指令Dαcx;
根据预置的攻角αcx0和附加攻角指令Dαcx,计算得到当前的实际攻角αcx。
较佳的,所述标准飞行轨迹指令包括:阻力加速度指令Dcx、阻力加速度导数指令标称速度指令Vcx和当地弹道倾角指令Θcx;
所述导航系统提供的测量值包括:Dsb、Vd和Θ。
较佳的,所述根据轨迹指令和导航系统提供的测量值,计算附加攻角指令Δαcx包括:
根据轨迹指令和导航系统提供的测量值,计算第一附加攻角指令分量αD;
根据轨迹指令和导航系统提供的测量值,计算第二附加攻角指令分量αΘ;
根据轨迹指令和导航系统提供的测量值,计算第三附加攻角指令分量αV;
根据αD、αΘ和αV,计算得到附加攻角指令Δαcx。
较佳的,通过如下的公式计算得到αD:
并根据预先设置的第一最小值和第一最大值对计算得到的αD进行限幅,使得所述αD满足条件:
其中,所述fd1和fd2分别为第一控制参数和第二控制参数。
较佳的,通过如下的公式计算得到第二附加攻角指令分量αΘ:
αΘ=-(Θ-Θcx)×fΘ;
并根据预先设置的第二最小值和第二最大值对计算得到的αΘ进行限幅,使得所述αΘ满足条件:其中,所述fΘ为第三控制参数。
较佳的,通过如下的公式计算得到第三附加攻角指令分量αV:
αV=-(Vcx-Vd)×fV;
并根据预先设置的第三最小值和第三最大值对计算得到的αV进行限幅,使得所述αV满足条件:其中,所述fV分别为第四控制参数。
较佳的,通过如下的公式计算得到附加攻角指令Δαcx:
Δαcx=αV+αD+αΘ。
较佳的,通过如下的公式计算得到αcx:
αcx=αcx0+Δαcx。
如上可见,在本发明中的高超声速飞行器滑翔飞行段的纵向制导方法中,由于在制导回路中闭环阻力加速度及其导数、速度和当地弹道倾角信息,对预置的标准攻角剖面进行调整,从而可以有效地改善阻力加速度跟踪控制的动态特性,提高动态条件下阻力加速度的跟踪控制精度,抑制飞行轨迹的波动,增强制导系统对各种偏差、不确定干扰和不确定条件的适应能力,可以应用于高超声速飞行器滑翔飞行段的纵向通道的制导。
附图说明
图1为本发明实施例中的高超声速飞行器滑翔飞行段的纵向制导方法的流程示意图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白,以下参照附图并举实施例,对本发明进一步详细说明。
本实施例提供了一种高超声速飞行器滑翔飞行段的纵向制导方法。
图1为本发明实施例中的高超声速飞行器滑翔飞行段的纵向制导方法的流程示意图。如图1所示,本发明实施例中的高超声速飞行器滑翔飞行段的纵向制导方法主要包括如下所述的步骤:
步骤101,在滑翔段飞行过程中,制导系统根据导航参数实时生成标准飞行轨迹指令。
较佳的,在本发明的具体实施例中,所述导航参数可以由导航系统提供。所述导航参数可以包括:待飞航程和速度信息。
较佳的,在本发明的具体实施例中,所述标准飞行轨迹指令可以包括:阻力加速度指令Dcx、阻力加速度导数指令标称速度指令Vcx和当地弹道倾角指令Θcx。
步骤102,制导系统根据标准飞行轨迹指令中的飞行速度或马赫给出预置的攻角αcx0。
在本发明的技术方案中,高超声速飞行器一般采用预置攻角飞行,且在飞行过程中,攻角不作较大调整,仅仅是在一定范围内调整,如下表示例所示:
速度Vcx | αcx0攻角(°) |
1000.0 | 10.0 |
2000.0 | 15.0 |
3000.0 | 15.0 |
4000.0 | 20.5 |
5000.0 | 25.0 |
6000.0 | 30.0 |
表1 标准攻角剖面
步骤103,根据轨迹指令和导航系统提供的测量值,计算附加攻角指令Δαcx。
较佳的,在本发明的具体实施例中,所述导航系统提供的测量值可以包括:阻力加速度Dsb、阻力加速度变化率相对运动速度Vd和当地弹道倾角Θ。
在本发明的技术方案中,可以通过多种方式来实现上述的步骤103,以下将以其中的一种具体实现方式为例,对本发明的技术方案进行详细的介绍。
例如,较佳的,在本发明的具体实施例中,所述步骤103可进一步包括:
步骤31,根据轨迹指令和导航系统提供的测量值,计算第一附加攻角指令分量αD。
在飞行器的飞行过程中,由于滚动角速率的限制,倾侧反转不可能快速实现,倾侧角反转(从负向变化到正向,或从正向变化的负向)穿越0度时,由于倾侧角小攻角大,因此将对该飞行器实际的阻力加速度产生一个较大的干扰,从而将对阻力加速度跟踪制导产生不利的影响。为改善飞行过程中由于倾侧角反转而带来的阻力加速度的动态误差大的问题,可根据阻力加速度及其导数计算出一个第一附加攻角指令分量αD,并根据该第一附加攻角指令分量αD对攻角进行适当的调整,以改善阻力加速度跟踪控制的动态误差。
较佳的,在本发明的具体实施例中,可以通过如下所述的公式计算得到αD:
并根据预先设置的第一最小值和第一最大值对计算得到的αD进行限幅,使得所述αD满足条件:
其中,所述fd1和fd2分别为第一控制参数和第二控制参数。在本发明的技术方案中,所述fd1和fd2的取值可以根据控制理论和仿真实验的结果预先确定或设置。
步骤32,根据轨迹指令和导航系统提供的测量值,计算第二附加攻角指令分量αΘ。
在飞行过程中,为抑制弹道的波动,还需要对弹道倾角进行适当的控制。因此,可以根据轨迹指令和导航系统提供的测量值,计算第二附加攻角指令分量αΘ。
较佳的,在本发明的具体实施例中,可以通过如下所述的公式计算得到第二附加攻角指令分量αΘ:
αΘ=-(Θ-Θcx)×f Θ (2)
并根据预先设置的第二最小值和第二最大值对计算得到的αΘ进行限幅,使得所述αΘ满足条件:
其中,所述fΘ为第三控制参数。在本发明的技术方案中,所述fΘ的取值可以根据控制理论和仿真实验的结果预先确定或设置。
步骤33,根据轨迹指令和导航系统提供的测量值,计算第三附加攻角指令分量αV。
为实现减速和增速控制效果,需要对攻角进行一定调整。例如,增大攻角有助于减速效果,而减小攻角则可补偿速度的不足。因此,根据轨迹指令和导航系统提供的测量值,计算第三附加攻角指令分量αV。
较佳的,在本发明的具体实施例中,可以通过如下所述的公式计算得到第三附加攻角指令分量αV:
αV=-(Vcx-Vd)×fV (3)
并根据预先设置的第三最小值和第三最大值对计算得到的αV进行限幅,使得所述αV满足条件:
其中,所述fV分别为第四控制参数。在本发明的技术方案中,所述fV的取值可以根据控制理论和仿真实验的结果预先确定或设置。
步骤34,根据αD、αΘ和αV,计算得到附加攻角指令Δαcx。
本发明中,附加攻角指令是指攻角的调节量,即变化量。
例如,较佳的,在本发明的具体实施例中,可以通过如下所述的公式计算得到纵向通道的附加攻角指令Δαcx:
Δαcx=αV+αD+αΘ (4)
由上可知,通过上述的步骤31~34,即可计算得到纵向通道的附加攻角指令Δαcx。
步骤104,根据预置的攻角αcx0和附加攻角指令Δαcx,计算得到当前的实际攻角αcx。
例如,较佳的,在本发明的具体实施例中,可以通过如下所述的公式计算得到αcx:
αcx=αcx0+Δαcx (5)
综上可知,在本发明中的高超声速飞行器滑翔飞行段的纵向制导方法中,由于在制导回路中闭环阻力加速度及其导数、速度和当地弹道倾角信息,对预置的标准攻角剖面进行调整,从而可以有效地改善阻力加速度跟踪控制的动态特性,提高动态条件下阻力加速度的跟踪控制精度,抑制飞行轨迹的波动,增强制导系统对各种偏差、不确定干扰和不确定条件的适应能力,可以应用于高超声速飞行器滑翔飞行段的纵向通道的制导。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明保护的范围之内。
Claims (6)
1.一种高超声速飞行器滑翔飞行段的纵向制导方法,其特征在于,该方法包括:
在滑翔段飞行过程中,制导系统根据导航参数实时生成标准飞行轨迹指令;
制导系统根据标准飞行轨迹指令中的飞行速度或马赫给出预置的攻角αcx0;
根据轨迹指令和导航系统提供的测量值,计算附加攻角指令Δαcx,该步骤包括:
根据轨迹指令和导航系统提供的测量值,计算第一附加攻角指令分量αD;
根据轨迹指令和导航系统提供的测量值,计算第二附加攻角指令分量αΘ;
根据轨迹指令和导航系统提供的测量值,计算第三附加攻角指令分量αV;
根据αD、αΘ和αV,计算得到附加攻角指令Δαcx;
根据预置的攻角αcx0和附加攻角指令Δαcx,计算得到当前的实际攻角αcx;
其中,所述标准飞行轨迹指令包括:阻力加速度指令Dcx、阻力加速度导数指令标称速度指令Vcx和当地弹道倾角指令Θcx;
所述导航系统提供的测量值包括:阻力加速度Dsb、阻力加速度变化率相对运动速度Vd和当地弹道倾角Θ。
2.根据权利要求1所述的方法,其特征在于,通过如下的公式计算得到αD:
并根据预先设置的第一最小值和第一最大值对计算得到的αD进行限幅,使得所述αD满足条件:
其中,所述fd1和fd2分别为第一控制参数和第二控制参数。
3.根据权利要求2所述的方法,其特征在于,通过如下的公式计算得到第二附加攻角指令分量αΘ:
αΘ=-(Θ-Θcx)×fΘ;
并根据预先设置的第二最小值和第二最大值对计算得到的αΘ进行限幅,使得所述αΘ满足条件:其中,所述fΘ为第三控制参数。
4.根据权利要求3所述的方法,其特征在于,通过如下的公式计算得到第三附加攻角指令分量αV:
αV=-(Vcx-Vd)×fV;
并根据预先设置的第三最小值和第三最大值对计算得到的αV进行限幅,使得所述αV满足条件:其中,所述fV分别为第四控制参数。
5.根据权利要求4所述的方法,其特征在于,通过如下的公式计算得到附加攻角指令Δαcx:
Δαcx=αV+αD+αΘ。
6.根据权利要求5所述的方法,其特征在于,通过如下的公式计算得到αcx:
αcx=αcx0+Δαcx。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510102948.2A CN104731104B (zh) | 2015-03-09 | 2015-03-09 | 一种高超声速飞行器滑翔飞行段的纵向制导方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510102948.2A CN104731104B (zh) | 2015-03-09 | 2015-03-09 | 一种高超声速飞行器滑翔飞行段的纵向制导方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104731104A CN104731104A (zh) | 2015-06-24 |
CN104731104B true CN104731104B (zh) | 2015-11-04 |
Family
ID=53455098
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510102948.2A Expired - Fee Related CN104731104B (zh) | 2015-03-09 | 2015-03-09 | 一种高超声速飞行器滑翔飞行段的纵向制导方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104731104B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104973250B (zh) * | 2015-07-02 | 2016-01-13 | 北京航天自动控制研究所 | 一种适用于滑翔飞行器的下压弹道攻角剖面确定方法 |
CN105676638B (zh) * | 2016-01-11 | 2018-06-29 | 北京航空航天大学 | 平稳滑翔/准自然频率跳跃滑翔组合机动突防弹道规划方法 |
CN105843232B (zh) * | 2016-04-08 | 2017-03-22 | 北京航天自动控制研究所 | 一种飞行器滑翔减速控制方法 |
CN106020216B (zh) * | 2016-05-13 | 2017-03-22 | 北京航天自动控制研究所 | 一种攻角约束下的平衡滑翔制导力分配方法 |
CN107941087B (zh) * | 2017-10-18 | 2019-08-06 | 北京航空航天大学 | 一种基于阻力剖面的高升阻比高超平稳滑翔再入制导方法 |
CN108398959A (zh) * | 2018-03-28 | 2018-08-14 | 北京临近空间飞行器系统工程研究所 | 一种高超声速滑翔飞行器快速下压制导控制方法 |
CN108536020B (zh) * | 2018-07-17 | 2019-06-21 | 哈尔滨工业大学 | 一种针对垂直起降重复使用运载器的模型参考自适应滑模控制方法 |
CN109240335B (zh) * | 2018-10-29 | 2021-09-07 | 北京控制工程研究所 | 一种空天飞行器进场着陆制导方法 |
CN112904888B (zh) * | 2021-01-11 | 2024-04-09 | 北京临近空间飞行器系统工程研究所 | 多目标参数联合制导的方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103587723A (zh) * | 2013-11-07 | 2014-02-19 | 北京临近空间飞行器系统工程研究所 | 一种再入初始段解析式纵向在线轨迹设计及跟踪方法 |
CN104176268A (zh) * | 2014-08-06 | 2014-12-03 | 北京航空航天大学 | 一种滑翔飞行弹道阻尼控制方法 |
-
2015
- 2015-03-09 CN CN201510102948.2A patent/CN104731104B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103587723A (zh) * | 2013-11-07 | 2014-02-19 | 北京临近空间飞行器系统工程研究所 | 一种再入初始段解析式纵向在线轨迹设计及跟踪方法 |
CN104176268A (zh) * | 2014-08-06 | 2014-12-03 | 北京航空航天大学 | 一种滑翔飞行弹道阻尼控制方法 |
Non-Patent Citations (1)
Title |
---|
基于平衡滑翔的升力再入攻角设计分析及应用;张冉等;《宇航学报》;20121130;第33卷(第11期);1557-1563 * |
Also Published As
Publication number | Publication date |
---|---|
CN104731104A (zh) | 2015-06-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104731104B (zh) | 一种高超声速飞行器滑翔飞行段的纵向制导方法 | |
CN107941087B (zh) | 一种基于阻力剖面的高升阻比高超平稳滑翔再入制导方法 | |
CN107844123B (zh) | 一种非线性飞行器航迹控制方法 | |
CN105159308A (zh) | 一种可重复使用飞行器着陆段制导与控制律一体化耦合设计方法 | |
CN104656450B (zh) | 一种高超声速飞行器平稳滑翔再入弹道设计方法 | |
CN103587723B (zh) | 一种再入初始段解析式纵向在线轨迹设计及跟踪方法 | |
CN105629734B (zh) | 一种近空间飞行器的轨迹跟踪控制方法 | |
CN103090728A (zh) | 一种基于滑模控制的带末角约束制导方法 | |
CN112327922B (zh) | 一种飞翼无人机自主起降综合控制方法 | |
CN103983143B (zh) | 包含速度过程约束和多终端约束的空地导弹投放下滑段制导方法 | |
CN103728976A (zh) | 一种基于广义标控脱靶量概念的多过程约束和多终端约束末制导律 | |
CN104503457B (zh) | 无人机编队飞行转弯防碰控制方法 | |
CN104199303B (zh) | 一种基于矢量场制导的平流层卫星平面路径跟踪控制方法 | |
CN102707624A (zh) | 基于飞行器常规模型的纵向控制器区域设计方法 | |
CN109752955A (zh) | 基于二维位置引导的飞行器轨迹跟踪与抗扰控制系统及方法 | |
CN104331084A (zh) | 一种基于方向舵控滚转策略的气动舵偏范围计算方法 | |
CN103245257A (zh) | 基于Bezier曲线的多约束飞行器导引方法 | |
CN103512426A (zh) | 一种次优的带末角约束制导方法 | |
CN105550402A (zh) | 一种基于攻角或倾侧角变频的高超平稳机动滑翔弹道设计方法 | |
CN106197173A (zh) | 基于扰动估计和补偿的战术导弹鲁棒姿态控制方法 | |
CN107121929A (zh) | 基于线性协方差模型预测控制的鲁棒再入制导方法 | |
CN104536448B (zh) | 一种基于Backstepping法的无人机姿态系统控制方法 | |
CN105022271A (zh) | 一种无人自主飞行器自适应pid控制方法 | |
CN110488875A (zh) | 基于动态逆的无人机跟踪目标初始段航向误差修正方法 | |
CN110209192A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20151104 Termination date: 20210309 |