V sistemu diferencialnih enačb, ki opisuje gibanje satelitov GLONASS (rus. Globalnaya Navigazionnaya Sputnikovaya Sistema), se pospešek, ki ga povzroča vpliv Lune in Sonca pogosto upošteva kot konstanta v obdobju integracije. V pričujoči raziskavi predpostavljamo, da pospešek zaradi vpliva Lune in Sonca ni konstanten, temveč se v obdobju integracije linearno spreminja. Po tej predpostavki izračunamo linearne funkcije v treh oseh pospeška zaradi vpliva Lune in Sonca za interval 30 minut in jih nato implementiramo v diferencialne enačbe. V dokumentu GLONASS-ICD (angl. Interface Control Document) je za določitev tirnice priporočena uporaba numerične integracije Runge-Kutta četrtega reda za reševanje sistema diferencialnih enačb. Izračun položaja in hitrosti satelita GLONASS je v tej študiji izveden z metodo Runge-Kutta četrtega reda pri integraciji naprej in nazaj, pri čemer so začetni pogoji navedeni v datoteki s satelita oddanih efemerid
|