卫星六根数参数预测卫星轨迹

以下是使用C++编写的基于卫星轨道六根数预测卫星轨迹的代码实现：

```cpp
#include
#include

using namespace std;

// 地球引力常数 (m³/s²)
const double MU = 3.986004418e14;

// 三维向量结构体
struct Vector3D {
    double x, y, z;
    Vector3D(double x_ = 0, double y_ = 0, double z_ = 0) : x(x_), y(y_), z(z_) {}
};

// 轨道六根数结构体
struct KeplerianElements {
    double a;        // 半长轴 (m)
    double e;        // 偏心率
    double i;        // 轨道倾角 (rad)
    double Omega;    // 升交点赤经 (rad)
    double omega;    // 近地点幅角 (rad)
    double M0;       // 历元平近点角 (rad)
    double t0;       // 历元时间 (s)
};

// 绕Z轴旋转
Vector3D rotate_z(const Vector3D& v, double theta) {
    double cos_t = cos(theta);
    double sin_t = sin(theta);
    return Vector3D(
        v.x * cos_t - v