航天器姿态容错控制：执行器故障与饱和的解决方案

1. 航天器姿态控制中的执行器容错挑战

航天器姿态控制系统是确保卫星、空间站等航天器在轨稳定运行的核心子系统。在实际太空环境中，执行机构（如反作用飞轮、控制力矩陀螺或推力器）常面临两类典型问题：物理饱和与突发故障。执行器饱和源于力矩输出能力的物理限制，当控制指令超出硬件极限时，系统动态特性将严重恶化；而执行器故障则包括效率下降（Loss of Effectiveness）和附加偏差（Additive Bias）两种形式，可能由机械磨损、电子元件老化或空间辐射等因素引发。

传统PD控制器在理想工况下表现良好，但面对上述非理想条件时存在明显局限：首先，固定增益难以适应执行器效能下降带来的控制力度不足；其次，积分项累积会导致饱和情况下的"wind-up"现象；更重要的是，缺乏故障检测机制使得系统难以及时响应突发异常。线性虚拟输入的反步控制（Backstepping Control）虽然通过非线性设计提升了鲁棒性，但对故障参数的时变特性仍缺乏自适应能力。

2. 系统建模与问题表述

2.1 四元数姿态动力学框架

采用四元数q=[q0,qv]描述姿态（q0为标量部，qv=[q1,q2,q3]为矢量部），其运动学方程为：

matlab复制function dq = kinematics(q,w)
    Omega = [0    -w(1) -w(2) -w(3);
             w(1)  0     w(3) -w(2);
             w(2) -w(3)  0     w(1);
             w(3)  w(2) -w(1)  0];
    dq = 0.5*Omega*q;
end

动力学方程考虑惯量矩阵J和外部扰动τd：

matlab复制function dw = dynamics(J,w,tau,tau_d)
    dw = inv(J)*(-cross(w,J*w) + tau + tau_d);
end

2.2 执行器故障建模

故障模型包含效率矩阵E(t)∈[0,1]^3和偏差向量b(t)：

matlab复制function tau_f = faulty_actuator(tau,E,b)
    tau_f = diag(E)*tau + b;