CN114967749B - 一种面向低成本高度表的机动弹道设计方法 - Google Patents
一种面向低成本高度表的机动弹道设计方法 Download PDFInfo
- Publication number
- CN114967749B CN114967749B CN202210672684.4A CN202210672684A CN114967749B CN 114967749 B CN114967749 B CN 114967749B CN 202210672684 A CN202210672684 A CN 202210672684A CN 114967749 B CN114967749 B CN 114967749B
- Authority
- CN
- China
- Prior art keywords
- height profile
- height
- design
- module
- altitude
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 22
- 238000005259 measurement Methods 0.000 claims abstract description 21
- 238000004364 calculation method Methods 0.000 claims abstract description 16
- 230000001133 acceleration Effects 0.000 claims abstract description 15
- 230000036461 convulsion Effects 0.000 claims abstract description 8
- 238000004519 manufacturing process Methods 0.000 claims abstract description 6
- 238000010586 diagram Methods 0.000 description 6
- 230000007547 defect Effects 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000011217 control strategy Methods 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/10—Simultaneous control of position or course in three dimensions
- G05D1/101—Simultaneous control of position or course in three dimensions specially adapted for aircraft
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Software Systems (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- General Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Operations Research (AREA)
- Automation & Control Theory (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Computing Systems (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明公开了一种面向低成本高度表的最优机动弹道设计方法,主要包括高度剖面设计模块、高度剖面计算模块、高度剖面应用模块。其中,高度剖面设计模块与高度剖面计算模块相连,用于构建关于高度‑剩余航程的函数关系,提取关键设计变量。高度剖面计算模块与高度剖面设计模块、高度剖面应用模块、低成本气压式高度表相连,用于根据给定的初始条件、终端条件以及设计参数,通过求解方程组得到高度剖面中涉及的多项式系数,进而确定高度剖面的最终形式。高度剖面应用模块与高度剖面计算模块相连,用于根据当前运动状态和高度剖面设计结果生成实时法向过载。与现有“弹道设计‑高度测量‑弹道修正”的常规思路不同,本发明能够为高度表提供高度方向速度、加速度、加加速度及更高阶导数信息,从弹道规划层面对高度测量设备进行补偿,能够实现飞行器控制精度和生产使用成本间的最佳平衡。
Description
技术领域
本发明属于飞行器总体设计领域,涉及一种适应低成本气压式高度表的在线弹道规划方法。
背景技术
本发明以低成本多空域作战导弹武器为背景,针对其飞行空域大、高度变化快、响应速度高等特点,解决飞行器高控制精度和低成本化间的关键矛盾问题。
对于高空大跨域飞行,通常采用气压式高度表实现实时高度测量,相应的测量信息将作为弹道控制***的反馈,进而实现在线制导和任务规划。气压式高度表能够根据大气环境实时给出飞行器的飞行高度,是弹上在线测量高度的一种最直接有效的手段,对飞行器的弹道、制导、导航和控制***极为重要,在航空航天、军工、民用等领域有着广泛的应用。气压式高度表能够适应较大的飞行空域和高度变化速率,但其在测量精度上的固有不足,需要通过技术手段予以补偿。
目前针对低成本气压式高度表的精度补偿方案主要有以下几种:
方案一:借助高精度测量仪、高精度惯组、GPS/北斗导航仪等设备来弥补。对于低成本飞行器而言,上述补偿手段通常是不可取的。
方案二:使用高度-气压非线性求解器,但受高度和大气环境间非定常、非线性关系的影响,无法满足飞行器的高精度控制需求。
现有基于低成本气压式高度表的弹道控制方法存在以下缺点:
缺点一:高度信息测量精度和实时性差。低成本气压式高度表在大气环境测量感知精度上稍显不足,并且不能适应飞行器在高度上较为快速、大幅的变化,无法为强机动高速飞行弹道模式提供准确的信息输入。
缺点二:提高高度测量信息精度的成本高。低成本高度表在测量精度上的固有不足,通常需要借助高精度测量仪、高精度惯组、GPS/北斗导航仪等设备来弥补;而这些补偿手段会大大增加高度测量***的制造使用成本,不符合经济性要求。
缺点三:高度-大气函数关系适应性差。修正海压法是目前较为常见的气压式高度表补偿修正方法,而函数拟合法是根据气压传感器实现快速高度计算并修正测量数据的常用手段。但实际上高度与气压的关系并非一成不变的,既无法解析确定,又不能采用一种通用的拟合方法全面表述,尤其是飞行器高速飞行且高度剧烈变化的情况。
缺点四:设计约束过强、影响弹道设计灵活性。为了使用低成本气压式高度表,通常需要在导弹武器***的设计过程中施加大量约束条件,如高度区间、高度变化率、高度方向加速度等;这些约束条件大大限制了弹道设计的可行域,无法满足高动态作战模式对弹道制导***提出的需求。
传统高度表精度补偿方式存在成本高、结构复杂、增加总重、抗扰性弱等问题,难以满足低成本作战使用方式。本发明提出一种适应低成本高度表的在线弹道规划方法。与现有“弹道设计-高度测量-弹道修正”的常规思路不同,本发明能够为气压式高度表提供高度方向上的速度、加速度、加加速度及更高阶导数信息,从弹道规划层面对高度测量设备进行补偿,能够实现飞行器控制精度和生产使用成本间的最佳平衡。
发明内容
本发明的目的是提供一种适应低成本高度表的在线弹道规划方法,具有较好的抗干扰效果、适用性和效率性,无需依赖额外的高成本测量装置。
本发明的目的是通过以下技术方案实现的:
(1)设计高度-剩余航程飞行剖面
设计飞行高度为剩余航程的函数:
其中,
n为多项式阶次,
a i 为待求系数,
K为修正系数。本发明取
n=5。
(2)生成高度表反馈信息
根据式,计算高度方向上的速度、加速度、加加速度及跟高阶导数。
(3)生成飞行控制信息
根据飞行器动力学方程计算高度对剩余航程的一/二阶导,然后求出控制过载。
(4)考虑终端约束及高度表使用约束的弹道迭代方法
设
F(
R L )为五阶多项式,确定飞行剖面须求出函数
F(
R L )中的系数,因此有6个未知数需要求解。通过构建6个方程来求解这些未知数。
1)根据当前高度和终端高度,可以得到两个边界条件方程。
2)根据当前飞行路径角和终端飞行路径角,可以得到两个边界条件方程。
3)设剩余航程1/3和2/3处(分别记为和)的高度值分别
h 1 和
h 2 ,则有:
由上述六个方程组成方程组,给定
h 1 和
h 2 后,可通过求解方程组得到多项式
F(
R L )中的系数。定义设计变量为:
(5)高度剖面在线修正方法
根据飞行器实时运动状态和剖面设计结果,计算法向过载,带入弹道控制***后递推至下一步运动状态;巡航上述步骤直至剩余航程为零。
由本发明的技术方案可以看出,本发明能够为气压式高度表提供高度方向速度、加速度、加加速度及更高阶导数信息,从弹道规划层面对高度测量设备进行补偿,能够实现飞行器控制精度和生产使用成本间的最佳平衡;同时仅含两个涉及参数(
h 1 ,
h 2 ),只需求解一个方程组即可完成弹道设计,整体计算量较低,利于实现基于高度测量信息的在线弹道重构,在信息来源和效率性方面均优于传统方案。
附图说明
为了清楚地说明本发明的技术方案,下面将介绍实施例附图。对于本领域的普通技术人员,在不付出创造性劳动的前提下,还可以根据这些附图获得其他附图。
图1为本发明实施例提供的基于高度-剩余航程飞行剖面的反馈控制策略示意图。
图2为本发明实施例提供的不同高度-剩余航程飞行剖面设计结果示意图。
图3为本发明实施例提供的高度变化率曲线设计结果示意图。
图4为本发明实施例提供的高度方向加速度曲线设计结果示意图。
图5为本发明实施例提供的法向控制过载曲线设计结果示意图。
图6为本发明实施例提供的弹道倾角变化率曲线设计结果示意图。
其中,图2~图6采用了相同的初始条件和终端条件,但采用了不同的剖面设计参数。对于本领域的普通技术人员,在不付出创造性劳动的前提下,通过改变初始条件、终端条件以及剖面设计参数(式中的
h 1 和
h 2 ),还可以获得其他附图。
具体实施方式
下面结合本发明技术中的附图,对本发明的技术方案进行详细描述;所描述的实施例并非本发明的全部实施例。基于本发明的实施例,本领域普通技术人员在未做出创造性劳动的前提下,所获得的其他实施例都属于本发明的保护范围。
本发明实施例提供一种面向低成本高度表的机动弹道设计方法,能够为低成本高度表提供高度方向速度、加速度、加加速度及更高阶导数信息,从弹道规划层面对高度测量设备进行补偿,能够实现飞行器控制精度和生产使用成本间的最佳平衡。主要包括:高度剖面设计模块、高度剖面计算模块、高度剖面应用模块;其中:
(1)高度剖面设计模块与高度剖面计算模块相连,用于构建关于高度-剩余航程的函数关系,提取关键设计变量;
设飞行器匀速运动,则其动力学方程为:
其中,
V为速度大小,
γ为弹道倾角,
h为飞行高度,
N y 为法向控制过载,
a e 为地球半径,
g为地球引力加速度大小,
R L 为剩余航程(从当前点到目标点的航程)。
设计飞行高度为剩余航程的函数,见式。
(2)高度剖面计算模块与高度剖面设计模块、高度剖面应用模块、低成本气压式高度表相连,用于根据给定的初始条件、终端条件以及设计参数,通过求解方程组得到高度剖面中涉及的多项式系数,进而确定高度剖面的最终形式;
确定飞行剖面须求出函数
F(
R L )中的系数。由于
F(
R L )为五阶多项式,因此有6个未知数需要求解,需构建6个方程。
1)根据当前高度和终端高度,可以得到两个方程:
2)根据当前飞行路径角和终端飞行路径角,可以得到两个方程:
3)展开式可得:
定义,则有,其中:
综上所述,确定飞行剖面仅需给出两个参数,即式中的
h 1 和
h 2 。
(3)高度剖面应用模块与高度剖面计算模块相连,用于根据当前运动状态和高度剖面设计结果生成实时法向过载。
根据飞行器动力学方程,高度对剩余航程的一阶导和二阶导可表示为:
根据动力学方程可得:
由此可得控制过载为:
如图1所示,整个高度剖面更新过程详述如下:
1)根据式计算法向过载,并带入弹道控制***;
2)根据动力学方程,递推得到下一时刻的飞行状态;
3)根据实时运动状态和终端约束,根据式~重新确定系数;
4)返回步骤1,直至剩余航程为零。
根据上述方案,得到的部分实施例如图2~图6所示。
图2给出了以剩余航程为自变量设计高度剖面,在不同设计参数下(
h 1 ,
h 2 ),飞行器的高度随剩余航程的变化历程。可以看出,在不同设计参数下,弹道设计结果均满足初始高度(15km)和终端高度(10km),并且满足航程约束(177.4km)。
图3给出了以剩余航程为自变量设计高度剖面,在不同设计参数下(
h 1 ,
h 2 ),飞行器的高度变化率随剩余航程的变化历程。可以看出,在不同设计参数下,飞行器的高度变化率连续且非常平滑,其变化率也是连续且平滑的,这对于低成本高度表的测量和解算都是很有利的。
图4给出了以剩余航程为自变量设计高度剖面,在不同设计参数下(
h 1 ,
h 2 ),飞行器的高度随时间变化的二阶导数(高度方向的加速度)随剩余航程的变化历程。可以看出,在不同设计参数下,飞行器在高度方向上的加速度连续且非常平滑,其变化率也是连续且平滑的,这对于低成本高度表的测量和解算都是很有利的,同时说明本发明提出的高度剖面设计方法可以提供高度方向上的高阶导数,能够为高度测量***提供更多可靠的高度信息。
图5给出了以剩余航程为自变量设计高度剖面,在不同设计参数下(
h 1 ,
h 2 ),飞行器的法向控制过载随剩余航程的变化历程。可以看出,在不同设计参数下,飞行器的控制过载连续且非常平滑,其变化率也是连续且平滑的,这对控制***实现飞行控制是很有利的,同时降低过载变化的剧烈程度也有利于降低外部载荷对高度测量装置的干扰。
图6给出了以剩余航程为自变量设计高度剖面,在不同设计参数下(
h 1 ,
h 2 ),飞行器的弹道倾角随剩余航程的变化历程。可以看出,在不同设计参数下,弹道设计结果均满足初始高度(0°)和终端高度(-10°),并且满足航程约束(177.4km);同时飞行器的弹道倾角连续且非常平滑,其变化率也是连续且平滑的,这一点也反映在法向过载上。
本发明具有如下优点:
(1)可根据当前剩余航程,快速输出高度方向速度、加速度、加加速度等高阶导数信息,用于补偿低成本高度表在测量信息上的不足;
(2)在线生成满足低精度高度表使用约束的高度剖面,整体高度变化趋势较为平滑,高度的高阶导数满足连续性和可导性;
(3)实现高度通道上的在线弹道修正,仅通过两个设计变量即可完成弹道重构。
上述仅为本发明诸多具体实施方式中的一种,本发明的保护范围不限于此,任何熟悉本技术领域的技术人员,在本发明所披露的技术范围内可轻易想到的变化或替换都涵盖在本发明的保护范围之内,本发明的保护范围应该以权利要求书的保护范围为准。
Claims (1)
1.一种面向低成本高度表的机动弹道设计方法,其特征在于,能够为低成本高度表提供高度方向速度、加速度、加加速度及更高阶导数信息,从弹道规划层面对高度测量设备进行补偿,能够实现飞行器控制精度和生产使用成本间的最佳平衡;
主要包括:高度剖面设计模块、高度剖面计算模块、高度剖面应用模块;
其中:
1)高度剖面设计模块与高度剖面计算模块相连,用于构建关于高度-剩余航程的函数关系,提取关键设计变量;
设飞行器匀速运动,则其动力学方程为:
其中,V为速度大小,γ为弹道倾角,h为飞行高度,Ny为法向控制过载,ae为地球半径,g为地球引力加速度大小,RL为剩余航程;
设计飞行高度为剩余航程的函数:
其中,n为多项式阶次,ai为待求系数,K为修正系数;
根据剩余航程的函数,计算高度方向上的速度、加速度、加加速度及跟高阶导数;
根据飞行器动力学方程计算高度对剩余航程的一/二阶导,然后求出控制过载;
2)高度剖面计算模块与高度剖面设计模块、高度剖面应用模块、低成本气压式高度表相连,用于根据给定的初始条件、终端条件以及设计参数,通过求解方程组得到高度剖面中涉及的多项式系数,进而确定高度剖面的最终形式;
设F(RL)为五阶多项式,确定飞行剖面须求出函数F(RL)中的系数;
根据当前高度和终端高度,可以得到两个方程:
根据当前飞行路径角和终端飞行路径角,可以得到两个方程:
设剩余航程1/3和2/3处的高度值分别h1和h2,
分别记为
RL1=RL0-(RL0-RLf)/3,和
RL2=RLf-(RL0-RLf)/3,则有:
由上述方程组成方程组,给定h1和h2后,可通过求解方程组得到多项式F(RL)中的系数,定义设计变量为:
Uh=[h1,h2]T;
3)高度剖面应用模块与高度剖面计算模块相连,用于根据当前运动状态和高度剖面设计结果生成实时法向过载。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210672684.4A CN114967749B (zh) | 2022-06-14 | 2022-06-14 | 一种面向低成本高度表的机动弹道设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210672684.4A CN114967749B (zh) | 2022-06-14 | 2022-06-14 | 一种面向低成本高度表的机动弹道设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114967749A CN114967749A (zh) | 2022-08-30 |
CN114967749B true CN114967749B (zh) | 2023-04-07 |
Family
ID=82964591
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210672684.4A Active CN114967749B (zh) | 2022-06-14 | 2022-06-14 | 一种面向低成本高度表的机动弹道设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114967749B (zh) |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109191921B (zh) * | 2018-08-31 | 2020-09-18 | 中国电子科技集团公司第二十八研究所 | 一种基于4d轨迹预测的空管监视数据仿真方法 |
CN113534847B (zh) * | 2021-08-23 | 2022-02-15 | 哈尔滨工业大学 | 一种有动力可重复使用飞行器能量管理轨迹设计方法 |
-
2022
- 2022-06-14 CN CN202210672684.4A patent/CN114967749B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN114967749A (zh) | 2022-08-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102645933B (zh) | 飞行器实现大空域机动的组合过载控制方法 | |
CN106681344A (zh) | 一种用于飞行器的高度控制方法及控制*** | |
CN104764467B (zh) | 空天飞行器惯性传感器误差在线自适应标定方法 | |
CN101441478A (zh) | 一种小型无人机自动着陆拉平控制方法及其装置 | |
CN107121929B (zh) | 基于线性协方差模型预测控制的鲁棒再入制导方法 | |
CN112683446B (zh) | 一种飞机实时重心位置估计方法 | |
CN101713654A (zh) | 跨音速飞行阶段大气攻角与惯性攻角的融合方法 | |
CN110553642B (zh) | 一种提高惯性制导精度的方法 | |
CN111290278B (zh) | 一种基于预测滑模的高超声速飞行器鲁棒姿态控制方法 | |
WO2021102669A1 (zh) | 超低轨道卫星轨道自主维持方法 | |
CN107844128A (zh) | 一种基于复合比例导引的高超声速飞行器巡航段制导方法 | |
Larson et al. | Qualitative evaluation of a flush air data system at transonic speeds and high angles of attack | |
CN113900448A (zh) | 一种基于滑模干扰观测器的飞行器预测校正复合制导方法 | |
CN112034879B (zh) | 一种基于高度-射程比的标准轨迹跟踪制导方法 | |
CN114608583A (zh) | 一种高超声速飞行器动态杆臂效应智能补偿方法 | |
CN114967749B (zh) | 一种面向低成本高度表的机动弹道设计方法 | |
Raab et al. | In-flight testing of MEMS pressure sensors for flight loads determination | |
CN111637878A (zh) | 无人机导航滤波器 | |
CN115574666B (zh) | 一种掠地巡航靶标定高方法 | |
CN109084772B (zh) | 一种基于无迹卡尔曼的视线转率提取方法及*** | |
CN110017809B (zh) | 利用地磁信息和光流传感器解算飞行器姿态的方法 | |
CN113720354A (zh) | 一种基于fbg应变花的递推式传递对准方法 | |
CN115047767A (zh) | 一种考虑高度表使用约束的俯冲弹道优化方法 | |
Haering Jr | Airdata calibration techniques for measuring atmospheric wind profiles | |
RU2649052C2 (ru) | Способ проведения метеорологической подготовки стрельбы артиллерии в горных условиях |
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 |