[发明专利]一种航天器惯性积在轨估计方法

说明书

一种航天器惯性积在轨估计方法，基于在轨扰动力矩估计方法，通过对重力梯度力矩的估计，准确获得了卫星惯性积的大小及符号，利用在轨估计方法获得的惯性积可有效消除惯性积不准确带来的附加补偿力矩，通过在轨扰动力矩实际估计获得卫星惯性积，相比地面仿真计算准确性更高，提高了卫星姿态控制性能。

技术领域

本发明涉及一种航天器惯性积在轨估计方法，属于卫星姿态控制技术领域。

背景技术

对于一类大惯量卫星，其转动惯量非对角线元素(惯性积)的数值也较大。因此，大型卫星惯性积不论是在整星稳态控制时的重力梯度补偿，还是姿态机动时的前馈力矩补偿，都有很重要的影响。然而，惯性积通常很难通过地面试验准确获得。

发明内容

本发明解决的技术问题是：针对目前现有技术中，传统大惯量卫星的惯性积难以通过地面试验准确获得的问题，提出了一种航天器惯性积在轨估计方法。

本发明解决上述技术问题是通过如下技术方案予以实现的：

一种航天器惯性积在轨估计方法，步骤如下：

(1)设计在轨扰动力矩估计器，根据卫星单轴角动量、卫星单轴姿态角速度、卫星单轴控制力矩计算观测器状态变量导数，同时根据卫星单轴角动量、观测器状态变量、观测器带宽计算卫星扰动力矩；

(2)设计针对陀螺信号噪声的结构滤波器，滤除陀螺高频噪声影响以提高估计结果平稳性；

(3)根据步骤(1)所得卫星扰动力矩，计算卫星惯性积Jyz、Jzx。

所述步骤(1)中，在轨扰动力矩估计器计算观测器状态变量导数与卫星扰动力矩的具体方法如下：

式中，J为卫星单轴角动量，ω为卫星单轴姿态角速度，l为观测器带宽，w为观测器状态变量，τ_c为卫星单轴控制力矩，为卫星扰动力矩，w.为观测器状态变量导数。

所述步骤(2)中，结构滤波器滤除陀螺高频噪声影响的具体公式为：

式中，ωz＝ωp，阻尼比ζpζz≥0，s为拉氏量。

所述步骤(3)中，卫星惯性积Jyz、Jzx的计算方法如下：

式中，为卫星X轴估计的扰动力矩，为卫星Z轴估计的扰动力矩，ω_o为卫星轨道角速度。

本发明与现有技术相比的优点在于：

本发明提供的一种航天器惯性积在轨估计方法，通过在轨扰动力矩估计方法获得重力梯度力矩的估计，不需要对卫星进行复杂的姿态机动即可获得yz项与xz项的惯性积，同时通过在轨扰动力矩实际估计获得卫星惯性积，相比地面仿真计算准确性更高，利用在轨估计方法获得的惯性积可有效消除惯性积不准确带来的附加补偿力矩，从而提高系统控制性能。

附图说明

图1为发明提供的错误惯性积计算所得X轴重力梯度力矩计算值和估计值示意图；

图2为发明提供的正确惯性积计算所得X轴重力梯度力矩计算值和估计值示意图；

具体实施方式

一种航天器惯性积在轨估计方法，通过在轨扰动力矩估计方法获得重力梯度力矩的估计，准确获得了惯性积的大小及符号，根据估计结果通过在轨注入的方式消除了惯性积不准确带来的附加补偿力矩，从而验证了本发明所述方法的准确性和有效性。

在轨估计的具体步骤如下：