二、空间碎片聚类-轨道计算与J2000坐标系实现

1. 整体思路

在空间碎片监测、卫星对地观测等任务中，需要精确知道卫星和空间目标在某一时刻的位置。通常我们使用开普勒轨道六要素（半长轴、偏心率、倾角、升交点赤经、近地点幅角、真近点角）来描述轨道，并通过轨道动力学外推得到任意时刻的位置。本文实现了一套基于J2000地心惯性系的轨道计算工具，核心功能包括：

• 根据初始真近点角计算平近点角，并利用平均角速度外推任意时刻的平近点角。

• 求解开普勒方程得到偏近点角，进而得到真近点角。

• 在轨道平面内计算位置坐标，再通过三次旋转（3-1-3欧拉角）变换到J2000坐标系。

• 封装为函数GetJ2000Position，输入当前时间、轨道要素、起始时间，输出位置结构体。

该模块是整个聚类分析系统的基础，后续的可见性判断、聚类计算都依赖准确的位置信息。

2. 代码解析

2.1 全局参数

在运行任何计算前，需要定义以下全局常数（通常放在主脚本或单独初始化）：

matlab

global mu angle2rad rad2angle earthRadius atmosphereHeight au; mu = 398600.4418; % 地心引力常数 km^3/s^2 angle2rad = pi / 180; % 角度转弧度 rad2angle = 180 / pi; % 弧度转角度 earthRadius = 6371.