Алгоритм расчета координат потребителя в системе GPS

Алгоритм расчета координат приемника потребителя поизмеренным псевдозадержкам и эфемеридной информацииПсевдозадержки и поправочные коэффициентыPDj – псевдозадержка дальномерного кода j-го спутникаaf0, af1, af2 – коэффициенты полинома поправок времениTGD – групповая задержка сигнала в аппаратуре спутникаЭфемеридная информация√A – корень квадратный из большой полуоси эллипса орбитые – эксцентриситет орбитыW – прямое восхождение восходящего узла орбиты спутникаW' – скорость изменения прямого восхождения восходящегоузла орбиты спутникаi – угол наклона плоскости орбиты к плоскости экватораw – аргумент перигеяi' – скорость изменения угла наклонаM0 – средняя аномалия на референcный моментDn – отклонение значения среднего движения отпредвычисленногоСuс и Cus – амплитуды косинусоидального и синусоидальногочленов в формуле для поправки в аргумент широтыСrс и Crs – амплитуды косинусоидального и синусоидальногочленов в формуле для поправки в радиус орбитыCic и Cis – амплитуды косинусоидального и синусоидальногочленов в формуле для поправки в угол наклона орбитыtOE – референсное время вычисления эфемеридtOC – референсное время вычисления поправок к часам спутника  T  t   PD  t T  t   T  t   DTDT  a  a  T  t   t   a  T  t   t   T  DTTjjt precsysrjprecjsvf2mjf0sysjmjprecf1jsvjprecsysOC2jprecOCGDRDTR  Ce0 A sin E k  Ce0 A sin Mk , где1/2C  2 m/cjt K  Tsys t precn  n 0  Dn 2  4.442807633 10t10с/мOEm A3 Dn14где m  3.986005 103м с2M K  M0  n  t KE K, n 1  E K, n E K, n  e0  sin E K, n   M K1  e0 cos E K, n1/2Алгоритм расчета координат и скорости движенияприемника потребителя (продолжение)K  arctg21  e0sin E K cos E K   e0 K  K  U K  CUC cos 2 K  CUS sin 2 KrK  A 1  e0 cos E K   rKrK  CRC cos 2K  CRS sin 2KVUK m e0 sin KA21  e0m 1  e0 cos KAгде 3  7.2921151467 105рад/сек.XSVK  rK  cos U K cos K  sin U K sin K cosi K U K   K  U KVRK K  0     3  t K  3 t OE21  e0i K  i0  i K   i  t Ki K  CIC cos 2 K  CIS sin 2 KYSVK  rK  cos U K sin K  sin U K cos K cos i K ZSVK  rK sin U K sin i KXSVK  VRK  cos UK cos K  sin U K sin K cosi K   VUK  sin UK cos K  cos UK sin K cosi K   3YSVK  VRK  cos U K sin  K  sin U K cos  K cosi K  YSVK VUK  sin U K sin  K  cos U K cos  K cosi K   3XSVKZSVK  VRK sin UK sini K  VUK cos UK sini KАлгоритм пересчета геоцентрических координат, X, Y, Zв геодезические: B – широту, L – долготу, h – высотуZK tg(B) 22X K  YKYK, h tg(L) XK2 N 1  Nh2XK2 YKcos(B)N-122N  a 1   sin B1 22;   2  2.