UTC, на яке виробляємо обчислення координат навігаційного супутника. Для отримання координат навігаційного супутника необхідно щомиті на кожному новому кроці додавати до t одну секунду).
У рівності розраховується середня аномалія навігаційного супутника для даного моменту часу:
M k=M 0 + nt k,
де M 0 - середня аномалія на початку відліку, яку видобувають із навігаційного повідомлення; n -Середня кутова швидкість рухомій точки, яка називається середнім рухом.
Середня аномалія для тіла, що рухається по необуреним орбіті, визначається множенням його середнього руху на інтервал часу після проходження перицентра.
Таким чином, середня аномалія - ??кутова відстань від перицентра гіпотетичного тіла, що рухається з постійною кутовою швидкістю, рівної середньому руху.
Середня кутова швидкість рухається точки n розраховується за формулою:
n=n 0 +? n,
де п ° - розрахункова середня рух, розраховується за формулою (6):
n 0 =? (?/A 3)
де А - велика піввісь орбіти навігаційного супутника, яку видобувають із навігаційного повідомлення, А=(? А) 2;
? =3.986005 € 10 14 м 3/с 2 (гравітаційна стала, помножена на масу землі);
? n - зміна середнього руху від початку відліку часу, передається в кадрі навігаційного повідомлення.
Методом ітерацій вирішується рівняння Кеплера для розрахунку ексцентричної аномалії E k за формулою
де е - ексцентриситет, який отримують із навігаційного повідомлення, початкове наближення E 0 визначається за формулою
Процес продовжується до виконання умови | Е до + 1 - E k | lt; 10 - 7
Для визначення зв'язку між правдивими і ексцентричними аномаліями використовуються формули
де ексцентрична аномалія, обчислена на кроці 3 при вирішенні рівняння Кеплера, е - ексцентриситет, який отримують із навігаційного повідомлення.
Далі, з урахуванням знаків sinv і cosv, визначається номер чверті кута v. І з урахуванням номера чверті для визначення v обчислюється або arcsin (v), або arccos (v)
U '= v + w,
де v - істинна аномалія, обчислена на кроці 5;
?- Аргумент перицентра - відстань перицентра від вузла, що витягають із навігаційного повідомлення.
За формулою розраховується аргумент широти, виправлений за рахунок гравітаційного обурення орбіти навігаційного супутника:
U=U + C us sin 2U + С uс. cos 2U ',
де C uc і C us витягуються з навігаційного повідомлення.
За рівняння обчислюється радіус-вектор навігаційного супутника для даної епохи (рівняння орбіти, задане за допомогою ексцентриситету):
s=A (1-ecosE) + C rs sin2U + C rc cos2U,
де C rc, C ks витягуються з навігаційного повідомлення.
За формулою обчислюється нахил обуреної площини орбіти до площини екватора
i=i 0 + C is sin 2U + C ic cos 2U + IDOT? < b align="justify"> t k ,
де C ic, C is, i 0, IDOT витягуються з навігаційного повідомлення.
За допомогою системи рівнянь обчислюється положення навігаційного супутника в орбітальній площині (напрям х - по висхідному вузлу, а y - перпендикулярне до висхідного вузла напрямок):
х =ScosU y =SsinU
За формулою обчислюється довгота висхідного вузла орбіти навігаційного супутника:
? =? 0 + (? -? e ) t k - ? e t ue
За допомогою системи рівнянь обчислюються прямокутні координати навігаційного супутника на момент часу t k з урахуванням
часу поширення сигналу від навігаційного супутника до споживача:
X=x cos (? +? e?) - y cos i sin (? +? e?),
у=x sin (? +? e?) + y cos i cos (? +? e?),
z=y 'sin i
Розглянемо рівняння псевдодальності наступно...
