1.本技术涉及火箭制导技术领域，尤其是涉及一种基于攻角变化的姿态控制系统设计方法。
背景技术：2.随着商业发射需求的增加，由于可回收液体运载火箭具有提高利用效率、降低飞行成本的优势，国内外研究机构及商业航天公司均开展了相应的研制工作。
3.可回收火箭在返回过程中发动机朝下，与传统的低阻力流线体飞行器气动特性差异较大，存在气动线性度差的问题，即不同攻角状态下气动力和力矩系数对攻角的导数变化较大。在进行姿态控制系统设计时，这一气动特性将带来较大的参数误差。
4.采用传统设计方法对误差的适应能力有限，在飞行环境不确定性较大的情况下，控制性能和姿态控制精度相对较差，故而有待改进。
技术实现要素：5.为了改善传统设计方法的控制性能和姿态控制精度相对较差等问题，本技术提供一种基于攻角变化的姿态控制系统设计方法。
6.本技术提供一种基于攻角变化的姿态控制系统设计方法，采用如下的技术方案：
7.一种基于攻角变化的姿态控制系统设计方法，包括以下步骤：
8.步骤一，计算飞行攻角包络；
9.步骤二，设计各时刻点下各攻角状态的校正网络；
10.步骤三，对每一特征时刻点重复步骤二，将各状态下的校正网络组合为参数矩阵装订于飞控软件；
11.步骤四，火箭飞行过程中实时辨识飞行攻角；
12.步骤五，利用当前飞行时刻及实时攻角插值校正网络矩阵用于姿态控制。
13.进一步地，步骤一的具体操作如下：
14.根据总体提供的风场数据使用要求，计算火箭飞行过程中的攻角包络，对于任一特征时刻点ti，存在攻角变化范围[0,α
imax
]，其中i取值为1到n，其中n表示特征时刻点个数，n的取值为正整数，特征时刻点ti的选取特征依据包括但不限于：火箭的飞行弹道和气动特性。
[0015]
进一步地，步骤二的具体操作如下：
[0016]
对特征时刻点ti，针对攻角变化范围在[0,α
imax
]内的每一个状态各设计一组校正网络参数g
ij
，其中j取值为1到m，m的取值为正整数，g
ij
的下角标i表示第i个特征时刻点，g
ij
的下角标j表示第j个攻角状态点，对每一个特征时刻点，攻角变化范围不同，选取的攻角状态点个数应该相同，即对每个时间点m为定值。
[0017]
进一步地，步骤三的具体操作如下：
[0018]
对每一个特征时刻点重复步骤二直至完成该飞行时段内所有特征点校正网络设
计，将每一个特征时刻点的每一个攻角状态下的校正网络参数组合为一个mxkxn的矩阵e
jz
，将矩阵e
jz
装订至飞控软件中，其中k表示任一状态下的校正网络参数的长度，且对每一个状态k值保持不变。
[0019]
进一步地，步骤四的具体操作如下：
[0020]
飞行过程中，从火箭启控开始实时辨识飞行攻角，时刻t
ri
的攻角辨识方法如下：通过火箭惯组反馈的箭体过载ny进行计算，具体计算公式如下：
[0021][0022]
式中，m为箭体质量，cn为法向力系数，q为动压，sm为参考面积，p为摆动发动机推力，为发动机俯仰通道摆角，为俯仰角加速度，xz为质心位置，xa为火箭惯组安装位置；
[0023]
在箭体过载ny已知的前提下，根据上述式可以求出时刻t
ri
的法向力系数cn，由于cn为随攻角变化的参数，以cn插值法向力系数矩阵，即可得到时刻t
ri
的飞行攻角α
ri
。
[0024]
进一步地，用于攻角辨识的箭体过载ny选取为h拍内的平均值，h的取值可以通过仿真试验进行确定。
[0025]
进一步地，步骤五的具体操作如下：
[0026]
以当前飞行时刻t
ri
和飞行攻角α
ri
插值校正网络矩阵e
jz
，得到当前状态下的校正网络参数g
jz_ri
用于火箭的姿态控制。
[0027]
综上所述，本技术包括以下至少一种有益技术效果：
[0028]
在飞行环境不确定性较大的情况下，能够实时辨识火箭飞行过程中的攻角变化，根据攻角大小进行校正网络调整，提高姿态控制系统的鲁棒性以及控制性能，姿态控制精度也相对较好，对误差的适应能力较好。
附图说明
[0029]
为了更清楚地说明本技术实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0030]
图1是本技术实施例的基于攻角变化的姿态控制系统设计方法的流程示意图。
具体实施方式
[0031]
为了使本技术所要解决的技术问题、技术方案及有益效果更加清楚明白，以下结合附图及实施例，对本技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本技术，并不用于限定本技术。
[0032]
需要说明的是，当元件被称为“固定于”或“设置于”另一个元件，它可以直接在另一个元件上或者间接在该另一个元件上。当一个元件被称为是“连接于”另一个元件，它可以是直接连接到另一个元件或间接连接至该另一个元件上。
[0033]
需要理解的是，术语“长度”、“宽度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本技术和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本技术的限制。
[0034]
此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。在本技术的描述中，“多个”的含义是两个或两个以上，除非另有明确具体的限定。
[0035]
以下结合附图1对本技术作进一步详细说明。
[0036]
本技术实施例公开一种基于攻角变化的姿态控制系统设计方法。参照图1，基于攻角变化的姿态控制系统设计方法包括以下步骤：
[0037]
步骤一，计算飞行攻角包络，具体操作如下：根据总体提供的风场数据使用要求，计算火箭飞行过程中的攻角包络，对于任一特征时刻点ti，存在攻角变化范围[0,α
imax
]，其中i取值为1到n，其中n表示特征时刻点个数，n的取值为正整数，特征时刻点ti的选取特征依据包括但不限于：火箭的飞行弹道和气动特性。
[0038]
步骤二，设计各时刻点下各攻角状态的校正网络，具体操作如下：对特征时刻点ti，针对攻角变化范围在[0,α
imax
]内的每一个状态各设计一组校正网络参数g
ij
，其中j取值为1到m，m的取值为正整数，g
ij
的下角标i表示第i个特征时刻点，g
ij
的下角标j表示第j个攻角状态点，对每一个特征时刻点，攻角变化范围不同，选取的攻角状态点个数应该相同，即对每个时间点m为定值。
[0039]
步骤三，对每一特征时刻点重复步骤二，将各状态下的校正网络组合为参数矩阵装订于飞控软件，具体操作如下：对每一个特征时刻点重复步骤二直至完成该飞行时段内所有特征点校正网络设计，将每一个特征时刻点的每一个攻角状态下的校正网络参数组合为一个mxkxn的矩阵e
jz
，将矩阵e
jz
装订至飞控软件中，其中k表示任一状态下的校正网络参数的长度，且对每一个状态k值保持不变，飞控软件即火箭的飞行控制软件。
[0040]
步骤四，火箭飞行过程中实时辨识飞行攻角，具体操作如下：飞行过程中，从火箭启控开始实时辨识飞行攻角，时刻t
ri
的攻角辨识方法如下：通过火箭惯组反馈的箭体过载ny进行计算，具体计算公式如下：
[0041][0042]
式中，m为箭体质量，cn为法向力系数，q为动压，sm为参考面积，p为摆动发动机推力，为发动机俯仰通道摆角，为俯仰角加速度，xz为质心位置，xa为火箭惯组安装位置；
[0043]
在箭体过载ny已知的前提下，根据上述式可以求出时刻t
ri
的法向力系数cn，由于cn为随攻角变化的参数，以cn插值法向力系数矩阵，即可得到时刻t
ri
的飞行攻角α
ri
。
[0044]
为避免箭体过载ny变化过快造成飞行攻角α
ri
突变，用于攻角辨识的箭体过载ny选取为h拍内的平均值，h的取值可以通过仿真试验进行确定。
[0045]
步骤五，利用当前飞行时刻及实时攻角插值校正网络矩阵用于姿态控制，具体操作如下：以当前飞行时刻t
ri
和飞行攻角α
ri
插值校正网络矩阵e
jz
，得到当前状态下的校正网络参数g
jz_ri
用于火箭的姿态控制。
[0046]
本技术实施例一种基于攻角变化的姿态控制系统设计方法的实施原理为：在飞行环境不确定性较大的情况下，能够实时辨识火箭飞行过程中的攻角变化，根据攻角大小进行校正网络调整，提高姿态控制系统的鲁棒性以及控制性能，姿态控制精度也相对较好，对误差的适应能力较好。
[0047]
以上所述仅为本技术的较佳实施例而已，并不用以限制本技术，凡在申请的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本技术的保护范围之内。