如何使用Skyfield生成模拟轨道？

pip install skyfield

from skyfield.api import load, Topos
from skyfield.data import mpc
from skyfield.framelib import ecliptic_frame

# 加载星历数据
data = load('de421.bsp')

# 加载天体对象
earth = data['earth']
sun = data['sun']
moon = data['moon']

# 创建观测地点对象
location = earth + Topos(latitude_degrees=37.7749, longitude_degrees=-122.4194)

# 计算轨道位置
ts = load.timescale()
t = ts.utc(2022, 1, range(1, 31))  # 设置时间范围

# 计算天体在观测地点的位置
astrometric = location.at(t).observe(sun)

# 获取天体的赤道坐标
ra, dec, distance = astrometric.radec()

# 获取天体的黄道坐标
lon, lat, distance = astrometric.frame_latlon(ecliptic_frame)

# 获取天体的视角直径
angular_diameter = sun.angular_diameter().degrees

# 打印轨道位置信息
for i in range(len(t)):
    print(f"时间：{t[i].utc_jpl()}，赤道坐标：{ra[i].hours}h {dec[i].degrees}°，黄道坐标：{lon[i].degrees}° {lat[i].degrees}°，视角直径：{angular_diameter[i]}°")