(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202211146730.3 (22)申请日 2022.09.21 (71)申请人 北京宇航系统工程研究所 地址 100076 北京市丰台区南大红门路1号 (72)发明人 赵永志　宋征宇　肖耘　吴义田　 徐珊姝　陈晓飞　袁赫　刘银　 李静琳　张普卓　杜昊昱　余光学　 朱海洋　王紫扬　邓舞燕　袁晗　 张荣升　 (51)Int.Cl. F42B 15/00(2006.01) F42B 15/01(2006.01) G06F 17/16(2006.01) G06F 30/15(2020.01) G06F 30/28(2020.01) G06F 113/08(2020.01)G06F 119/14(2020.01) (54)发明名称 非轴对称运载火箭减载控制方法及存储设 备 (57)摘要 本发明提供的非轴对称运载火箭减载控制 方法及存储设备， 包括获取标准弹道设计数据、 发射点地理信息数据、 发射窗口预报高空风数 据； 计算时间 时刻的发射数据、 高空风数据以 及发射系到箭体系的转换矩阵， 地心系到发射系 的转换矩阵， 当地水平系到地心系的转换矩阵； 计算 时刻的单位风矢量在发射系的投影； 计算 时间序列对应的滚动程序角理论值序列； 对滚动 程序角理论值序列进行数值拟合得到工程化的 滚动程序角装订值序列； 计算滚动程序角装订值 序列 的一阶差分得到滚动程序角速率序列； 将 滚动程序角插值表、 滚动程序角速率插值表作为 诸元装订至运载火箭上； 具有实现降低飞行摆角 需求和降低飞行载荷的目的， 适用于运载火箭领 域。 权利要求书3页 说明书11页 附图11页 CN 115265292 A 2022.11.01 CN 115265292 A 1.非轴对称运载火箭减载控制方法， 其特 征在于： 包括以下步骤： 获取标准弹道设计数据、 发射 点地理信息数据、 发射窗口预报高空风数据； 初始化火箭 一级飞行 段特征的时间序列， 记为： ； 基于标准 弹道设计数据、 发射窗口预报高空风数据， 计算时间 时刻的发射数据和高空 风数据； ； 计算时间 时刻的发射系 到箭体系 的转换矩阵 ， 地心系 到发射系 的转换 矩阵 ， 当地水平系 到地心系 的转换矩阵 ； 计算 时刻的单位 风矢量在发射系的投影； 通过滚动程序角理论 值 ， 计算时间序列 对应的滚动程序角理论 值序列 ； 对滚动程序角理论值序列 进行数值拟合， 并对程序角曲线切入和切出部分采用二 次曲线进行过渡， 得到 工程化的滚动程序角装订 值序列 ； 计算滚动程序角装订值序列 一阶差分,并补0作为第一个元素， 得到滚动程序角速率 序列 ； 将时间序列 分别与滚动程序角装订值序列 、 滚动程序角速率序列 一一对应， 组 合得到滚动程序角插值表、 滚动程序角速率插值表； 将滚动程序角插值表、 滚动程序角速率插值表作为诸元装订至运载火箭上。 2.根据权利要求1所述的非轴 对称运载火箭减载控制方法， 其特征在于： 所述标准弹道 设计数据包括： 随时间变化的发射系 俯仰姿态角插值表 ， 发射系偏航姿态角插值 表 ， 发射系滚动姿态角插值表 ,箭下点大地经度插值表 ， 箭下点 地心纬度插值表 ， 飞行高度插值表 ， 射向 ； 所述发射 点地理信息数据包括： 发射 点大地经度 ， 发射点地理纬度 ； 所述发射 窗口预报高空风数据包括： 火箭发射窗口实测的随飞行高度变化的风速插值 表 和风向插值表 。 3.根据权利要求2所述的非轴对称运载火箭减载控制方法， 其特征在于： 所述 时刻的 发射数据包括： 发射系俯仰姿态角 ,发射系偏航姿态角 ,发射系滚动姿态角 , 箭下点大地经度 ,箭下点地心纬度 ,飞行高度 ； 所述 时刻的高空风数据包括： 高度处的风速 和风向。权　利　要　求　书 1/3 页 2 CN 115265292 A 24.根据权利要求3所述的非轴对称运载火箭减载控制方法， 其特征在于： 所述 时刻的 发射系 到箭体系 的转换矩阵的计算表达式为：  式 （1） 所述 时刻的地心系 到发射系 的转换矩阵的计算表达式为：  式 （2） 所述 时刻的当地水平系 到地心系 的转换矩阵的计算表达式为：  式 （3） 。 5.根据权利要求4所述的非轴对称运载火箭减载控制方法， 其特征在于： 所述 时刻的 单位风矢量在发射系的投影的计算表达式为：  式 （4） 式 （4） 中， 表示 时刻的单位 风矢量在发射系的投影。 6.根据权利要求5所述的非轴 对称运载火箭减载控制方法， 其特征在于： 所述滚动程序 角理论值 的计算表达式为： 式 （5） 式 （5） 中，权　利　要　求　书 2/3 页 3 CN 115265292 A 3