2. 轨道控制对象建模

轨道控制仿真的第一步，不是写控制器，而是把“被控对象”描述清楚。被控对象由两部分组成：一是航天器所处的轨道动力学环境，二是航天器自身的质量、推力、比冲等参数。前者决定轨道如何随时间演化，后者决定控制指令能产生多大的实际加速度。两者共同构成仿真中的状态方程与输入通道，任何一处建模偏差都会直接反映为控制精度下降。

一、轨道动力学模型：从二体到多源摄动

最基础的模型是二体问题。在惯性坐标系下，设地球引力常数为 μ，航天器位置矢量为 r，则加速度为 a = −μr/|r|³。二体模型给出理想的开普勒轨道，能量与角动量守恒，轨道为闭合椭圆。它适合用于初步验证控制律、估算速度增量量级，但真实轨道会因摄动而偏离。

J2 摄动来自地球非球形引力场，是低地球轨道最主要的摄动源。其加速度可写为：

a_J2 = −(3/2)·J2·μ·Re²/r⁵ · [ (1−5(z/r)²)·x,
                              (1−5(z/r)²)·y,
                              (3−5(z/r)²)·z ]
其中 J2 ≈ 1.0826×10⁻³，Re 为地球赤道半径。

J2 的典型效应是升交点赤经漂移、近地点幅角旋转和平近点角变化。对太阳同步轨道，正是利用 J2 使升交点赤经每天进动约 0.9856° 来保持光照条件。仿真中若忽略 J2，太阳同步轨道的设计将完全失效。

大气阻力作用于低轨航天器，加速度为 a_drag = −½·ρ·(C_dA/m)·|v_rel|·v_rel。其中 ρ 为大气密度，C_d 为阻力系数，A 为迎风面积，m 为质量。大气密度随高度、太阳活动、地方时变化显著，工程上常用 NRLMSISE-00 或指数模型近似。阻力导致轨道半长轴衰减，是轨道维持必须补偿的主要项。

三体引力在深空或地月转移中不可忽略。以地月系统为例，月球引力加速度为 a_moon = −μ_m(r−r_m)/|r−r_m|³ + μ_mr_m/|r_m|³，第二项为间接项。在近地轨道，三体项量级约为 10⁻⁶ m/s²，通常可忽略；但在平动点附近，它成为主导项。

二、航天器参数建模：质量、推力与比冲

质量是连接动力学与控制的枢纽。仿真中质量不是常数，而是随推进剂消耗变化的状态量。设初始质量 m₀，推进剂质量 m_p，则 m(t) = m₀ − ∫ṁ dt。推进剂消耗率由推力与比冲决定：ṁ = F/(Isp·g₀)，其中 g₀ = 9.80665 m/s²。例如一台 490 N 发动机、Isp = 310 s，则 ṁ ≈ 490/(310×9.80665) ≈ 0.161 kg/s。若点火 100 s，消耗约 16.1 kg 推进剂。

推力模型需要区分连续推力与脉冲推力。连续推力用于电推进或长弧点火，仿真中直接作为加速度输入；脉冲推力用于化学推进的短时点火，常以速度增量 Δv 描述。实际建模时还需考虑推力上升/下降时间、最小点火时长、推力方向安装偏差等非理想因素。

比冲反映推进效率。化学推进 Isp 约 280～320 s，电推进可达 1500～3000 s。比冲越高，相同 Δv 下推进剂消耗越少，但推力通常越小。仿真中应根据任务阶段选择对应参数，而不是全程使用单一值。

三、模型选择与仿真实现要点

模型选择取决于任务精度与轨道高度。对于 LEO 轨道维持仿真，至少应包含二体 + J2 + 大气阻力；对于太阳同步轨道设计，J2 必须保留；对于地月转移，需加入月球三体引力；对于高轨长期演化，还需考虑太阳光压与第三体摄动。

在数值实现上，常用四阶龙格-库塔法积分状态量 [r, v, m]。积分步长需与轨道周期匹配：LEO 周期约 90 min，步长可取 1～10 s；深空轨道可放宽至 60 s 以上。步长过大将丢失 J2 短周期项，步长过小则增加计算量。

完成上述建模后，被控对象即可用一组微分方程描述。下一步才是控制律设计与仿真验证。建模阶段多花时间核对参数与摄动项，往往比后期反复调参更有效。