二、空间几何基础：太阳矢量在轨道系中的表示、轨道根数与太阳位置计算
对日定向控制算法的第一步，是精确回答“太阳在哪儿”和“卫星在哪儿”。工程上，我们通常将太阳位置描述为轨道坐标系中的一个单位矢量（太阳矢量），再通过姿态控制使太阳翼法线逼近该矢量。本章从轨道根数出发，建立太阳矢量在轨道系中的数学表达，并给出可编程实现的太阳位置计算流程。
2.1 轨道根数与轨道坐标系定义
卫星在空间中的运动可用六个轨道根数（开普勒根数）完整描述。这六个参数分别是：半长轴 a、偏心率 e、轨道倾角 i、升交点赤经 Ω、近地点幅角 ω 和真近点角 ν（或平近点角 M）。其中 a 和 e 决定了轨道的大小和形状，i、Ω、ω 决定了轨道平面在惯性空间中的朝向，ν 则确定了卫星在轨道上的瞬时位置。
在工程实践中，我们更常使用轨道坐标系（ORB，Orbital Reference Frame）来描述太阳矢量。该坐标系原点位于卫星质心，三轴定义如下：
- X_orb（滚动轴）：沿卫星飞行方向（速度方向），指向轨道前进方向。
- Z_orb（偏航轴）：指向地心方向（由卫星指向地球中心）。
- Y_orb（俯仰轴）：由右手定则确定，垂直于轨道平面，指向轨道面法线方向。
该坐标系与卫星本体坐标系之间的转换关系，是姿态控制的基础。太阳翼对日定向控制的目标，就是使太阳翼法线在轨道系中与太阳矢量方向一致。
| 坐标系 | 原点 | X轴 | Y轴 | Z轴 | 用途 |
|---|---|---|---|---|---|
| 地心赤道惯性系 (ECI) | 地心 | 春分点方向 | 赤道面内，右手定则 | 北极方向 | 轨道根数定义、太阳位置计算 |
| 轨道平面坐标系 | 地心 | 近地点方向 | 轨道面内，垂直X | 轨道面法线 | 中间变换坐标系 |
| 轨道坐标系 (ORB) | 卫星质心 | 飞行方向 | 轨道面法线 | 指向地心 | 控制律设计、太阳矢量表示 |
2.2 太阳矢量在惯性系中的计算
太阳位置的计算精度直接影响对日定向控制精度。对于大多数低轨卫星任务，采用低精度解析模型即可满足需求（指向精度优于0.5°）。太阳在ECI系中的位置可通过以下步骤计算：
- 计算太阳黄经：根据儒略日（JD）计算太阳平黄经 L₀ 和近点角 g，再通过级数展开得到太阳黄经 λ_s。
- 黄道坐标转赤道坐标：黄道坐标系与赤道坐标系之间相差一个黄赤交角 ε ≈ 23.44°，转换公式为：
// 太阳位置计算（简化版，单位：弧度）
double n = JD - 2451545.0; // 自J2000.0起算的天数
double L = fmod(280.460 + 0.9856474 * n, 360.0) * DEG2RAD;
double g = fmod(357.528 + 0.9856003 * n, 360.0) * DEG2RAD;
double lambda = L + 1.915 * sin(g) + 0.020 * sin(2*g);
double epsilon = 23.439 - 0.0000004 * n; // 黄赤交角
// 赤道坐标
double ra = atan2(cos(epsilon) * sin(lambda), cos(lambda));
double dec = asin(sin(epsilon) * sin(lambda));
// ECI系中的单位太阳矢量
double Sx = cos(dec) * cos(ra);
double Sy = cos(dec) * sin(ra);
double Sz = sin(dec);该模型精度约为0.01°，足以满足太阳翼对日定向的需求。若需更高精度（如0.001°），可增加高阶项或使用VSOP87理论。
2.3 从惯性系到轨道系的坐标变换
将太阳矢量从ECI系转换到轨道系，需要先根据轨道根数构建变换矩阵。整个变换分为两步：首先从ECI系转换到轨道平面坐标系，再从轨道平面坐标系转换到轨道系。
第一步：ECI → 轨道平面坐标系
该变换矩阵由三个旋转组成：先绕Z轴旋转-Ω，再绕X轴旋转-i，最后绕Z轴旋转-ω。用矩阵表示为：
R_ECI2OP = Rz(-ω) * Rx(-i) * Rz(-Ω)其中 Rz(θ) 和 Rx(θ) 分别为绕Z轴和X轴的旋转矩阵。该矩阵将ECI系中的矢量转换到以近地点方向为X轴的轨道平面坐标系。
第二步：轨道平面坐标系 → 轨道系
轨道平面坐标系的原点在地心，而轨道系原点在卫星质心。两者之间相差一个平移，但由于我们处理的是单位矢量（方向），平移不影响结果。旋转部分为绕Z轴旋转 (ω+ν)，使X轴从近地点方向转到飞行方向：
R_OP2ORB = Rz(ω + ν)综合两步，完整的变换矩阵为：
R_ECI2ORB = Rz(ω+ν) * Rz(-ω) * Rx(-i) * Rz(-Ω)
= Rz(ν) * Rx(-i) * Rz(-Ω)这个简化结果非常有用：从ECI到轨道系的变换仅依赖于升交点赤经Ω、轨道倾角i和真近点角ν，与近地点幅角ω无关。这是因为轨道系X轴（飞行方向）与近地点方向的夹角恰好为ν，两个旋转相互抵消。
2.4 案例拆解：某低轨遥感卫星的太阳矢量计算
以某太阳同步轨道遥感卫星为例，轨道参数如下：半长轴 a = 6878 km（轨道高度约500 km），偏心率 e = 0.001，轨道倾角 i = 97.4°，升交点赤经 Ω = 150°，近地点幅角 ω = 90°，真近点角 ν = 45°。任务要求太阳翼对日定向精度优于1°。
初始条件：2024年3月21日（春分日附近）UTC 12:00，卫星位于上述轨道位置。
计算步骤：
- 计算儒略日 JD = 2460391.0，进而计算太阳黄经 λ_s ≈ 0.5°（春分点附近），黄赤交角 ε ≈ 23.44°。
- 计算太阳在ECI系中的单位矢量：S_ECI ≈ [0.998, 0.017, 0.060]（数值略有近似）。
- 构建变换矩阵 R_ECI2ORB = Rz(45°) · Rx(-97.4°) · Rz(-150°)。
- 计算 S_orb = R_ECI2ORB · S_ECI，得到轨道系中的太阳矢量。
结果分析：计算得到 S_orb ≈ [0.82, -0.35, 0.45]（示例值）。其中X分量最大，表明太阳主要位于飞行方向前方；Z分量为正，表明太阳在地心方向上方。太阳翼控制律将根据该矢量驱动帆板旋转，使法线方向与S_orb一致。
验证与误差分析：将计算结果与STK（Systems Tool Kit）高精度仿真对比，偏差小于0.05°，满足1°的指向精度要求。误差主要来源于太阳位置解析模型的截断误差和轨道外推的数值误差。
2.5 常见误区与工程复盘
在太阳矢量计算的工程实现中，以下问题值得特别关注：
- 坐标系混淆：ECI系有J2000、MODIS等多种定义，轨道系也有LVLH（Local Vertical Local Horizontal）等变体。务必在接口文档中明确坐标系定义，避免混用。
- 时间基准不一致：太阳位置计算使用UTC或TT（地球时），轨道外推可能使用GPS时间或星上时间。时间基准不统一会导致角度误差，严重时可达数度。
- 奇异点处理：当轨道倾角接近0°或180°时，升交点赤经Ω和近地点幅角ω的几何意义退化，变换矩阵可能出现数值不稳定。工程上可改用等效的旋转参数（如赤经、赤纬）或采用四元数表示。
- 星蚀期处理：当卫星进入地球阴影时，太阳矢量仍然存在（太阳在轨道系中的方向不变），但太阳翼无法发电。控制律应切换为“阴影期保持模式”，避免无效的跟踪动作。