import numpy as npimport matplotlibpyplot as pltfrom mpl_toolkitsmplot3d import Axes3Dfrom astroqueryjplhorizons import Horizons# 从JPL Horizons数据库中获取木星的轨道数据obj = Horizonsid=599 location=sun epochs=sta

如果无法连接JPL Horizons数据库，可以考虑使用其他的天体轨道数据源，比如NASA的HORIZONS系统。以下是使用NASA HORIZONS系统获取木星和其四个卫星的轨道数据并进行可视化的代码示例：

import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
from astroquery.jplhorizons import Horizons

# 从NASA HORIZONS系统中获取木星的轨道数据
obj = Horizons(id='5', location='@sun', epochs={'start':'2022-01-01', 'stop':'2022-12-31', 'step':'1d'}, id_type='majorbody')
eph = obj.ephemerides()

# 提取木星的位置和速度数据
jupiter_pos = np.array([eph['x'], eph['y'], eph['z']])
jupiter_vel = np.array([eph['vx'], eph['vy'], eph['vz']])

# 获取木卫一至木卫四的轨道数据
moons = ['Io', 'Europa', 'Ganymede', 'Callisto']
moon_pos = []
moon_vel = []

for moon in moons:
    obj = Horizons(id=moon, location='@5', epochs={'start':'2022-01-01', 'stop':'2022-12-31', 'step':'1d'}, id_type='majorbody')
    eph = obj.ephemerides()
    moon_pos.append([eph['x'], eph['y'], eph['z']])
    moon_vel.append([eph['vx'], eph['vy'], eph['vz']])

# 创建交互式可视化
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
ax.set_title('Jupiter and its Moons')

# 绘制木星轨道
ax.plot(jupiter_pos[0], jupiter_pos[1], jupiter_pos[2], label='Jupiter')

# 绘制木卫一至木卫四的轨道
for i in range(len(moons)):
    ax.plot(moon_pos[i][0], moon_pos[i][1], moon_pos[i][2], label=moons[i])

# 设置坐标轴标签
ax.set_xlabel('X [AU]')
ax.set_ylabel('Y [AU]')
ax.set_zlabel('Z [AU]')

# 添加图例
ax.legend()

# 显示可视化
plt.show()

请注意，这里使用的是NASA HORIZONS系统，所以id_type参数设置为majorbody，并且木星的ID为5。使用此代码示例，您应该能够获取木星和其四个卫星的轨道数据并进行可视化