Frågan Tracking SpaceIL Beresheet Moon lander länkar till vissa datafiler som kommer från webbplatsen: http://live.spaceil.com (även länkad här):
rådata:
http://live.spaceil.com/data/data_s1.txt ( 22/2 - 12/4)
http://live.spaceil.com/data/data_s2.txt (4/4 - 12/4)
http://live.spaceil.com/data/data_m.txt (22/2 - 16/12)
Här är plot data i data_s1.txt
(rymdskepp) och data_m.txt
(Moon). Det finns korrekta tillståndsvektorer för rymdfarkosten $ [x, y, z, v_x, v_y, v_z] $ men endast positionsvektorer för månen. $ [x, y, z] $ .
Det finns dock ett verkligt problem; de stämmer inte korrekt! Rymdfarkosten kretsar inte om månen. Inte ens nära!
Jag tror att det kan finnas något som är knepigt med data. felet är mycket större än 1/81 så det är mer än bara barycenter kontra jordcentrets koordinater. Något är fel.
Fråga: Varför står inte SpaceILs rymdfarkoster och månbanor i linje?
Python-skript:
importera numpy som npimport matplotlib.pyplot som pltfrån mpl_toolkits.mplot3d import Axes3Dwith open ('data_s1.txt', 'r') som infile: slines = infile.readlines () slines = [[float (x) för x i line.strip (). split (',') [1:]] för rad i linjer [1: -1]] sstates = np.array (zip (* slines)) med öppen ('data_m.txt', 'r') som infil: mlines = infile.readlines () mlines = [[float (x) för x i linje.strip (). delad (',') [1:]] för rad i mlines [1: -1]] mstates = np.array (zip (* mlines)) om True: fig = plt.figure (figsize = [10, 8]) # [12, 10] ax = fig.add_subplot (1, 1, 1, projection = '3d') x, y, z = sstates [: 3] ax. plot (x, y, z, '-r') x, y, z = mstatus [: 3] ax.plot (x, y, z, '-b') ax.set_xlim (-405000, 405000) ax. set_ylim (-405000, 405000) ax.set_zlim (-405000, 405000) plt.show ()