升力式飞行器的再入标准轨迹指令的确定方法
【技术领域】
[0001] 本发明涉及升力式飞行器技术，特别涉及指一种升力式飞行器的再入标准轨迹指 令的确定方法。
【背景技术】
[0002] 对于具有高升阻比特性的升力式滑翔飞行器，现有技术中一般采用基于标准轨迹 跟踪制导方案，标准轨迹剖面可选取阻力加速度-速度剖面形式，在跟踪过程中需要实时 计算阻力加速度指令、阻力加速度指令对时间的一阶导数、阻力加速度指令对时间的二阶 导数等参数。在实时计算阻力加速度指令、阻力加速度指令对时间的一阶导数、阻力加速度 指令对时间的二阶导数等参数的过程中经常采用差分求导的方式。然而，采用差分求导时 可能存在指令突变或非可导等问题出现，导致跟踪制导时参数设计困难，且产生的倾侧角 指令抖动的问题。
[0003] 因此，现有技术中存在对基于标准剖面的跟踪制导方案中能够产生跟踪所需的轨 迹指令及轨迹指令的导数项从而避免跟踪制导参数设计困难的需要。
【发明内容】
[0004] 本发明是一种适用于升力式飞行器的再入标准轨迹指令计算方法，主要应用于标 准剖面跟踪制导方案中，用于产生跟踪所需的轨迹指令及轨迹指令的导数项，通过本发明 所提方法，能够实时生成跟踪制导所需的轨迹指令，减小在线指令生成复杂度，降低对跟踪 控制律参数设计的敏感度。
[0005] 为解决现有技术中存在的技术问题，本发明的实施例提供一种升力式飞行器的再 入标准轨迹指令的确定方法，包括：根据确定的再入标准轨迹来确定高度一速度、阻力加速 度一速度的映射关系；利用数学拟合工具，通过数学拟合来确定高度一速度、阻力加速度一 速度的解析形式的映射表达式，并根据所确定的映射表达式来确定阻力加速度指令DcJP 高度指令Ηα，
[0006] Dcx= b ! · Ve+b2 · Ve2+b3 · Ve3+b4 · Ve4+b5 公式（I)
[0007] Hcx= c ! · Ve+c2 · Ve2+c3 · Ve3+c4 · Ve4+c5 公式（2)
[0008] 其中，Dcx为阻力加速度指令，H α为高度指令，V e= V d/V。，Vd为飞行器对地速度，V c 为归一化常数，其值为第一宇宙速度7910m/s，Id1, b2, b3, b4, b5, C1, c2, c3, c4, C5为待定系数；
[0009] 利用数学拟合来确定高度一速度、阻力加速度一速度的解析形式的映射表达式中 的待定系数的拟合值，得到确定的阻力加速度指令Dc3、高度指令1与速度的拟合映射表达 式；
[0010] 根据确定的阻力加速度指令Dc3、高度指令Hcx与速度的拟合映射表达式，确定再入 标准轨迹的阻力加速度导数?。、、二阶导数?εχ和当地弹道倾角指令Θεχ。 toon] 其中，利用数学拟合来确定确定高度一速度、阻力加速度一速度的解析形式的映 射表达式中的待定系数的拟合值包括：
[0012] 将预先获取的一组标准轨迹参数映射表代入公式（1)和（2)，并利用数学拟合工 具，在给定精度下求解出所述待定系数h b2, b3, b4, b5, C1, c2, c3, c4, C5所分别对应的拟合值
[0013] 将确定的拟合值分别代入公式（1)和（2)替换掉待定系数Id1, b2, b3, b4, b5, C1, c2, c3，c4, c5，得到确定的阻力加速度指令Dc3、高度指令Hcx与速度的拟合映射表达式。
[0014] 其中，根据确定的阻力加速度指令Dcx、高度指令Hcx与速度的拟合映射表达式，确 定再入标准轨迹的阻力加速度导数?εχ、二阶导数De\和当地弹道倾角指令?^的步骤包 括：
[0015] 基于公式（1)，分别确定阻力加速度导数]^和二阶导数￡^，以及当地弹道倾角 指令I。
[0016]
[0017] 公式 (4)
[0018] 基于公式⑵和下述公式（5)，确定当地弹道倾角指令
[0019]
[0020] 其中，基于平衡滑翔原理，升力式飞行器在大气层内滑翔飞行时dVd/dt ~_DCX，
[0021] 本发明的上述技术方案的有益效果如下：
[0022] 上述方案中，针对升力式飞行器再入滑翔过程中的再入标准指令生成问题，先通 过数学拟合确定再入标准轨迹的阻力加速度一速度、高度一速度映射关系，进而基于平衡 滑翔原理，通过简化速度、倾角与阻力加速度的关系，直接采用解析形式得到了阻力加速度 导数、二阶导数和当地弹道倾角指令的计算表达式，避免了多次分段高阶拟合上述导数项， 有效降低了计算复杂度，提高了设计的通用性，可直接用于滑翔飞行过程中的制导方案设 计中。
【附图说明】
[0023] 图1为本发明实施例的升力式飞行器的再入标准轨迹指令的确定方法的流程图。
【具体实施方式】
[0024] 为使本发明要解决的技术问题、技术方案和优点更加清楚，下面将结合附图及具 体实施例进行详细描述。
[0025] 针对升力式飞行器滑翔飞行过程中的再入标准指令生成问题，本发明提供了适用 于升力式飞行器的再入轨迹指令的确定方法，避免了多次分段高阶拟合轨迹指令导数项， 有效降低计算复杂度，提高设计通用性，为滑翔飞行过程中的制导系统设计提供指令。
[0026] 如图1所示，示出了根据本发明的升力式飞行器的再入标准轨迹指令的确定方 法。在步骤100,确定升力式飞行器的再入标准轨迹，并根据确定的标准轨迹确定高度一速 度、阻力加速度一速度的映射关系。具体地，根据升力式飞行器的再入的过程和终端约束条 件和飞行器总体参数、气动特性，计算标准轨迹，并根据计算的标准轨迹确定高度一速度、 阻力加速度一速度映射关系。根据本发明的实施例，在升力式飞行器的滑翔段飞行过程中， 需满足严格的热流、过载、动压等过程约束，以及终端速度、终端高度、航程等约束条件限 制，并基于平衡滑翔原理，采用多种方式快速生成一条满足过程约束和终端约束的再入标 准轨迹，从而得到高度-速度剖面和阻力加速度-速度剖面标准轨迹。升力式飞行器的再 入标准轨迹的设计方法属于现有技术，在此不作详细介绍。
[0027] 在步骤102,利用数学拟合工具，通过数学拟合来确定高度一速度、阻力加速度一 速度的解析形式的映射表达式，并根据所确定的映射表达式来确定阻力加速度指令DcJP 高度指令Hcx。本发明的实施例中，利用数学拟合工具，对高度一速度、阻力加速度一速度映 射关系进行分析，可找出合适的映射关系表达式，如：
[0028] Dcx= b ! · Ve+b2 · Ve2+b3 · Ve3+b4 · Ve4+b5 公式（I)
[0029] Hcx= c ! · Ve+c2 · Ve2+c3 · Ve3+c4 · Ve4+c5 公式（2)
[0030] 其中，Dcx为阻力加速度指令，H α为高度指令，V e= V d/V。，Vd为飞行器对地速度，V c 为归一化常数，较佳地可取为第一宇宙速度7910m/s，Id1, b2, b3, b4, b5, C1, c2, c3, c4,