Рассчитать положение и скорость спутника
Я использовал этот код, но есть некоторая ошибка в конечном результате положения и скорости, потому что я использую тот же TLE и эпоху STK для проверки моего результата.
Результат STK:https://drive.google.com/file/d/1fMTmDLBMc2n7DRlevAUPFmzhAoVIjoq5/view?usp=sharing
from sgp4.api import Satrec,jday
s = '1 44792U 98067QX 21309.33315742 .00311477 00000-0 12417-2 0 9999'
t = '2 44792 51.6322 309.1372 0007423 71.1960 288.9849 15.85183673111730'
satellite = Satrec.twoline2rv(s, t)
#5 Nov 2021 07:59:44.801 GMT
jd, fr = jday(2021,11,5,7,59,44.801)
e,r, v = satellite.sgp4(jd,fr)
print(r) # True Equator Mean Equinox position (km)
print(v) # True Equator Mean Equinox velocity (km/s)
вывод кода: для r -> (4224.540458337468, -5191.398656066208, 0.0054742096980043954) для v -> (3.7080450717209392, 3.0331606529487622, 6.0556337625858685)