一种基于路径规划的挠性航天器姿态自适应控制方法技术

发明专利技术涉及一种基于路径规划的挠性航天器姿态自适应控制方法，首先基于SMPPa姿态机动路径规划方法，对航天器的期望角度进行柔化操作；然后基于特征模型思想的挠性航天器数学模型，利用梯度下降法进行参数在线辨识，确定三轴模型特征参数；接着根据确定的挠性航天器数学模型及其特征参数，确定控制力矩，根据控制力矩控制挠性航天器姿态，最后重复上述步骤，直至角度达到航天器期望角度。本发明专利技术抑制了挠性附件振动对控制性能的影响，提高了姿态控制的稳态精度和动态特性，适用于具有三轴耦合非线性特点的挠性航天器。

An adaptive attitude control method for flexible spacecraft based on Path Planning

【技术实现步骤摘要】
一种基于路径规划的挠性航天器姿态自适应控制方法
本专利技术涉及挠性航天器姿态控制技术，特别是涉及一种基于路径规划的挠性航天器姿态自适应控制方法。
技术介绍
随着空间技术的高速发展和广泛应用，人们对卫星在数据通信传输能力上的要求进一步提高。为了达到这一要求，现代卫星大多配置有大面积的太阳帆板等种类繁多的挠性附件。同时，为适应空间任务的多样化要求，现代航天器通常需要具备快速姿态跟踪或大角度姿态机动的能力。但目前针对挠性航天器这一多输入多输出且具有不确定性的非线性强耦合系统，研究主要集中在对挠性航天器的动力学建模、控制器的设计以及对挠性振动的抑制上，很少涉及对输入的姿态机动路径进行规划。近年来，基于吴宏鑫院士提出的特征建模及黄金分割理论的卫星姿态机动控制的研究，现已取得了较大的进展。由于基于特征建模思想的控制方法不依赖于精确的航天器动力学模型，因此适用于对挠性结构进行控制。但目前的研究主要集中于单轴挠性航天器，对存在三轴耦合非线性特点的挠性航天器大角度机动问题的研究较少。例如，有学者提出了基于特征模型思想对特殊的SISO高阶线性定常系统进行建模，用于挠性结构的控制中，验证了其控制的有效性，但未考虑MIMO的非线性系统的挠性结构控制。有学者提出了一种将特征模型自适应控制应用到反馈线性化控制中的姿态控制策略，考虑了高速航天器强耦合的模型特点，但未考虑高频未建模动态的鲁棒性。
技术实现思路
本专利技术的目的在于提供一种基于路径规划的挠性航天器姿态自适应控制方法，抑制了挠性附件振动对控制性能的影响，提高了姿态控制的稳态精度和动态特性，适用于具有三轴耦合非线性特点的挠性航天器。实现本专利技术目的的技术解决方案为：一种基于路径规划的挠性航天器姿态自适应控制方法，包括如下步骤：步骤1、基于SMPPa姿态机动路径规划方法，对航天器的期望角度进行柔化操作；步骤2、基于特征模型思想的挠性航天器数学模型，利用梯度下降法进行参数在线辨识，确定三轴模型特征参数；步骤3、根据确定的挠性航天器数学模型，确定控制力矩；步骤4、根据控制力矩控制挠性航天器姿态；步骤5、重复步骤2-4，直至角度达到航天器期望角度。本专利技术与现有技术相比，其显著优点为：1)本专利技术利用基于SMPPa姿态机动路径规划方法进行柔化，将柔化后的指令信号作为控制系统的期望输入信号，解决挠性附件振动耦合的问题，可以减少太阳帆板的振动的激发，缓解角速度和角加速度的阈值问题；2)本专利技术针对挠性航天器MIMO非线性强耦合系统，基于特征建模思想建立了的挠性航天器数学模型，在该模型的基础上提出了黄金分割自适应姿态控制方法，能够加强航天器姿态系统的自适应性和鲁棒性，使系统过渡过程趋于平稳，加快了系统的跟踪速度，减小了系统超调和克服挠性附件振动。附图说明图1为基于路径规划的挠性航天器姿态自适应控制示意图。图2为基于余弦函数的加速度无跳变的对称型机动路径(SMPPa)流程图。图3为基于余弦函数的加速度无跳变的对称型机动路径(SMPPa)规划图。图4为基于路径规划的挠性航天器姿态自适应控制参数在线辨识流程图。图5为基于路径规划的挠性航天器姿态自适应控制三轴角度仿真结果图。图6为基于路径规划的挠性航天器姿态自适应控制三轴角度误差仿真结果图。图7为基于路径规划的挠性航天器姿态自适应控制三轴角速度误差仿真结果图。图8为基于路径规划的挠性航天器姿态自适应控制振动模态曲线仿真结果图。图9为基于路径规划的挠性航天器姿态自适应控制方法总体流程图。具体实施方式下面结合附图和具体实施方式进一步说明本专利技术方案。本专利技术首先分析了MIMO非线性强耦合系统的挠性航天器动力学模型，给出了其姿态控制系统的总体结构；然后，为有效抑制挠性附件振动对控制性能的影响，提出了SMPPa姿态机动路径规划方法对指令信号进行柔化，以减少指令信号中的高频成分；然后，采用基于MIMO特征建模理论，设计了基于黄金分割自适应跟踪控制方法；最后将输入期望值经过SMPPa姿态规划方法柔化后的控制信号输入给设计的黄金分割自适应跟踪控制方法，控制挠性航天器完成航天器姿态控制；而对于挠性航天器特征模型的参数辨识采用了梯度下降法的在线辨识。此处先分析MIMO非线性强耦合系统的挠性航天器动力学模型。采用欧拉角对带有大面积太阳帆板和天线等种类繁多的挠性附件的航天器姿态控制系统进行描述，按照y-x-z的旋转方式，可得其运动学方程：基于挠性航天器动力学模型，建立挠性航天器混合坐标下动力学方程：式中，θ,ψ分别为航天器俯仰轴、横滚轴和偏航轴的角度，J∈R3×3为航天器的转动惯量矩阵；ω＝[ω1ω2ω3]T∈R3为刚体转角速度矢量；C0∈R3×3为太阳帆板振动和中心刚体的耦合矩阵；Tc为作用在刚体上的控制力矩；Td＝[Td1Td2Td3]T为作用在刚体上的干扰力矩；η∈Rn×1为太阳帆板振动模态坐标；ξ＝diag(ξ1,ξ2,...,ξn)为太阳帆板振动模态阻尼比；Λ＝diag(Λ1,Λ2,...,Λn)为太阳帆板振动模态频率矩阵；ω×为矢量ω的斜对称矩阵。根据上述过程建立的模型，给出了图1所示挠性航天器姿态控制方法结构框图，从图1可以看出该路径规划的挠性航天器姿态自适应控制方法主要完成路径规划方法设计、特征模型建立、参数在线辨识、自适应控制方法设计，具体设计如下所述。基于路径规划的挠性航天器姿态自适应控制方法，包括如下步骤：步骤1、基于SMPPa姿态机动路径规划方法，对航天器的期望角度进行柔化操作。对于挠性航天器姿态大角度机动过程中近似阶跃的跳变，本专利技术提出了SMPPa姿态机动路径规划方法对输入的指令信号进行柔化，可以减少太阳帆板振动的激发，而且可以缓解角速度和角加速度的阈值问题。假设余弦函数型加速度的最大值为acos_max、周期为2Tcos1；角速度匀速段的时间为Tcos2、角速度为Vcos_max；起始角位置为(度)、目标角位置为(度)。若要求的姿态机动总角度为(度)，则(度)。由于受航天器实际运行过程中测量机构的限制，角速度的测量值存在一个阈值，故在对航天器进行姿态机动路径规划时，对于角加速度和角速度的值要求有：acos_max∈(0,aallowed_max)，Vcos_max∈(0,Vallowed_max)之内，并且要求总的机动时间tcos_max＝2Tcos1+Tcos2尽可能的达到最小；要求航天器姿态机动过程结束后，其角加速度和角速度值最终趋于零，姿态角位置最终稳定在给定角度下面详细介绍SMPPa的规划步骤。实现SMPPa规划的前提条件：已知总机动角度角速度加速段的余弦函数型加速度的最大值acos_max、姿态机动可以达到的最大角加速度值aallowed_max、可达到的最大角速度Vallowed_max和加速段时间Tcos1。余弦角加速度的路径规划(SMPPa)方法实现步骤如下：步骤1.1、设置余弦函数型加速度的最大值为acos_max、周期2Tcos1、姿态机动可以达到的最大角加速度值aallowed_max、可达到的最大角速度Vallowed_max和要求的姿态机动总角度步骤1.2、计算姿态机动路径所能达到的最大角速度值Vcos_max：步骤1.3、根据Vcos_max允许达到的上限Vallowed_max调整acos_max，若计算得到的Vco本文档来自技高网...

【技术保护点】
一种基于路径规划的挠性航天器姿态自适应控制方法，其特征在于，包括如下步骤：步骤1、基于SMPPa姿态机动路径规划方法，对航天器的期望角度进行柔化操作；步骤2、基于特征模型思想的挠性航天器数学模型，利用梯度下降法进行参数在线辨识，确定三轴模型特征参数；步骤3、根据确定的挠性航天器数学模型及其特征参数，确定控制力矩；步骤4、根据控制力矩控制挠性航天器姿态；步骤5、重复步骤2‑4，直至角度达到航天器期望角度。

【技术特征摘要】
1.一种基于路径规划的挠性航天器姿态自适应控制方法，其特征在于，包括如下步骤：步骤1、基于SMPPa姿态机动路径规划方法，对航天器的期望角度进行柔化操作；步骤2、基于特征模型思想的挠性航天器数学模型，利用梯度下降法进行参数在线辨识，确定三轴模型特征参数；步骤3、根据确定的挠性航天器数学模型及其特征参数，确定控制力矩；步骤4、根据控制力矩控制挠性航天器姿态；步骤5、重复步骤2-4，直至角度达到航天器期望角度。2.根据权利要求1所述的基于路径规划的挠性航天器姿态自适应控制方法，其特征在于，步骤1基于SMPPa姿态机动路径规划方法，确定SMPPa的角位置轨迹具体方法为：步骤1.1、设置余弦函数型加速度的最大值为acos_max、周期2Tcos1、姿态机动可以达到的最大角加速度值aallowed_max、可达到的最大角速度Vallowed_max和要求的姿态机动总角度步骤1.2、计算姿态机动路径所能达到的最大角速度值Vcos_max：步骤1.3、根据Vcos_max允许达到的上限Vallowed_max调整acos_max，若计算得到的Vcos_max＞Vallowed_max，令Vcos_max＝Vallowed_max，根据新的Vallowed_max值计算得到较小的acos_max，保证姿态角速度的最大值不超速，acos_max更新公式为：步骤1.4、计算姿态机动路径所能达到的最小角步骤1.5、根据要求的姿态机动总角度匀速段时间Tcos2，若根据下式计算匀速段时间Tcos2：若令匀速段时间Tcos2＝0；步骤1.6、根据匀速段时间Tcos2更新Tcos1，进而更新acos_max和Vcos_max，更新acos_max的公式为：步骤1.7、根据确定的路径参数acos_max、Tcos1、Tcos2、Vcos_max，以及余弦函数的积分规律，得到SMPPa的角位置轨迹3.根据权利要求1所述的基于路径规划的挠性航天器姿态自适应控制方法，其特征在于，步骤2确定三轴模型特征参数的具体方法为：步骤2.1、基于特征模型思想建立挠性航天器数学模型：式中，表示当前的三轴角度的状态矩阵，表示当前俯仰轴、横滚轴和偏航轴的角度输出值，表示下一时刻的三轴角度的状态矩阵，表示下一时刻俯仰轴、横滚轴和偏航轴的角度输出值，表示上一时刻的三轴角度的状态矩阵，表示上一时刻俯仰轴、横滚...

【专利技术属性】
技术研发人员：吴益飞，芦梦兰，高熠，关妍，陈庆伟，郭健，郭毓，金立伟，刘小壮，周梦兰，范成旺，陈鑫，
申请(专利权)人：南京理工大学，
类型：发明
国别省市：江苏,32