Способ определения ортогональных составляющих вектора скорости и способ определения координат космического аппарата с использованием земных станций
Владельцы патента RU 2750228:
Федеральное государственное казенное военное образовательное учреждение высшего образования "Военная академия связи имени Маршала Советского Союза С.М. Буденного" Министерства обороны Российской Федерации (RU)
Заявленная группа изобретений относится к области космонавтики, а именно к технике выполнения траекторных измерений, определения координат и ортогональных составляющих вектора скорости космического аппарата (КА), и могут быть использованы на наземных и бортовых комплексах управления полетом КА для точного определения текущих параметров движения КА. Техническим результатом является повышение точности определения координат и ортогональных составляющих вектора скорости КА. Способ определения вектора скорости основного КА включает: измерение в наземной радиотехнической станции (НРТС) K значений номиналов частот принятых радиосигналов n-х земных станций (ЗС) In после их ретрансляции основным S1 и смежным S2 КА соответственно, расчет координат основного КА х1, у1, z1, вычисление ортогональных составляющих вектора скорости основного КА
Способ определения координат основного КА х1, у1, z1 включает: измерение в НРТС K для каждой n-й ЗС In значений временных задержек Δtn между принятыми радиосигналами после их ретрансляции основным S1 и смежным КА S2 соответственно, расчет координат основного КА x1, y1, z1. 2 н. и 2 з.п. ф-лы, 2 прил., 7 ил.
Заявленные объекты объединены единым изобретательским замыслом и заявленный способ определения координат космического аппарата (КА) с использованием земных станций (ЗС) предназначен для реализации способа определения ортогональных составляющих вектора скорости КА с использованием ЗС.
Изобретения относятся к области космонавтики, а именно к технике выполнения траекторных измерений и определения параметров орбиты КА, и может быть использовано на наземных и бортовых комплексах управления полетом КА для точного определения текущих параметров движения КА.
Известен способ определения ортогональных составляющих вектора скорости КА [1]. Для реализации данного способа последовательно выполняют следующие этапы:
размещают на позициях с известными координатами НРТС и не мене трех ПОРС;
принимают и записывают радиосигналы КА совместно с метками времени с помощью НРТС и указанных ПОРС в момент времени t0;
передают из каждой ПОРС в НРТС записанные радиосигналы КА совместно с метками времени;
с помощью корреляционной обработки измеряют значения частотных сдвигов между радиосигналами, записанными НРТС и каждой из ПОРС;
рассчитывают разности радиальных скоростей КА относительно НРТС и каждой из ПОРС;
вычисляют составляющие вектора скорости КА используя указанные разности радиальных скоростей.
К недостаткам способа [1] относят:
длительное время определения ортогональных составляющих вектора скорости КА, обусловленное необходимостью размещать на позициях с известными координатами не менее трех ПОРС;
относительно-невысокую точность определения ортогональных составляющих вектора скорости КА, обусловленную необходимостью синхронной записи радиосигналов КА с помощью НРТС и не менее трех ПОРС.
Известен способ определения ортогональных составляющих вектора скорости КА [2]. Для реализации данного способа последовательно выполняют следующие этапы:
размещают на позициях с известными координатами НРТС и не мене двух ППОРС;
в момент времени t0 синхронно излучают тестовые радиосигналы с помощью НРТС и указанных ППОРС;
принимают и записывают тестовые радиосигналы после их ретрансляции КА с помощью НРТС и указанных ППОРС;
с помощью корреляционной обработки измеряют значения частотных сдвигов между переданными и принятыми радиосигналами для НРТС и каждой из ППОРС;
рассчитывают радиальные скорости КА относительно НРТС и каждой из ППОРС;
передают из каждой ППОРС в НРТС рассчитанное значение радиальной скорости;
вычисляют ортогональные составляющие вектора скорости КА используя рассчитанные радиальные скорости.
К недостаткам способа [2] относят:
длительное время определения ортогональных составляющих вектора скорости КА, обусловленное необходимостью размещать на позициях с известными координатами не менее двух ППОРС;
относительно-невысокую точность определения ортогональных составляющих вектора скорости КА, обусловленную необходимостью синхронной излучения и записи тестовых радиосигналов КА с помощью НРТС и не менее двух ППОРС.
Из известных способов наиболее близким аналогом (прототипом) предлагаемого способа по технической сущности является способ определения ортогональных составляющих вектора скорости КА [3]. Для реализации данного способа последовательно выполняют следующие этапы:
размещают на позициях с известными координатами наземную радиотехническую станцию (НРТС) и не мене двух излучающих опорных реперных станций (ИОРС);
в момент времени t0 синхронно излучают тестовые радиосигналы с помощью НРТС и указанных ИОРС;
принимают и записывают тестовые радиосигналы после их ретрансляции КА с помощью НРТС;
с помощью корреляционной обработки измеряют значения частотных сдвигов между переданными и принятыми радиосигналами для НРТС и каждой из ИОРС;
рассчитывают суммы радиальных скоростей КА относительно НРТС и каждой из ИОРС;
вычисляют составляющие вектора скорости КА используя указанные суммы радиальных скоростей.
Недостатками способа прототипа [3] являются:
длительное время определения ортогональных составляющих вектора скорости КА, обусловленное необходимостью размещать на позициях с известными координатами не менее двух ИОРС;
относительно-невысокая точность определения ортогональных составляющих вектора скорости КА обусловленная необходимостью синхронной излучения тестовых радиосигналов с помощью НРТС и не менее двух ИОРС.
Известен способ определения координат КА [1]. Для реализации данного способа последовательно выполняют следующие этапы:
размещают на позициях с известными координатами наземную радиотехническую станцию (НРТС) и не менее трех приемных опорных реперных станций (ПОРС);
принимают и записывают радиосигналы КА совместно с метками времени с помощью НРТС и указанных ПОРС в момент времени t0;
передают из каждой ПОРС в НРТС записанные радиосигналы КА совместно с метками времени;
с помощью корреляционной обработки измеряют значения взаимных временных задержек между радиосигналами, записанными НРТС и каждой из ПОРС;
вычисляют координаты КА используя указанные разности дальностей.
К недостаткам способа [1] относят:
длительное время определения координат КА, обусловленное необходимостью размещать на позициях с известными координатами не менее трех ПОРС;
относительно-невысокую точность определения координат КА, обусловленную необходимостью синхронной записи радиосигналов КА с помощью НРТС и не менее трех ПОРС.
Известен способ определения координат КА [2]. Для реализации данного способа последовательно выполняют следующие этапы:
размещают на позициях с известными координатами НРТС и не мене двух приемо-передающих опорных реперных станций (ППОРС);
в момент времени t0 синхронно излучают тестовые радиосигналы с помощью НРТС и указанных ППОРС;
принимают и записывают тестовые радиосигналы после их ретрансляции КА с помощью НРТС и указанных ППОРС;
с помощью корреляционной обработки измеряют значения взаимных временных задержек между переданными и принятыми радиосигналами для НРТС и каждой из ППОРС;
рассчитывают дальности от КА до НРТС и каждой из ППОРС;
передают из каждой ППОРС в НРТС рассчитанное значение дальности;
вычисляют координаты КА используя рассчитанные дальности.
К недостаткам способа [2] относят:
длительное время определения координат КА, обусловленное необходимостью размещать на позициях с известными координатами не менее двух ППОРС;
относительно-невысокую точность определения координат КА, обусловленную необходимостью синхронной излучения и записи тестовых радиосигналов КА с помощью НРТС и не менее двух ППОРС.
Из известных способов наиболее близким аналогом (прототипом) предлагаемого способа по технической сущности является способ определения координат КА [3]. Для реализации данного способа последовательно выполняют следующие этапы:
размещают на позициях с известными координатами НРТС и не мене двух излучающих опорных реперных станций (ИОРС);
в момент времени t0 синхронно излучают тестовые радиосигналы с помощью НРТС и указанных ИОРС;
принимают и записывают тестовые радиосигналы после их ретрансляции КА с помощью НРТС;
с помощью корреляционной обработки измеряют значения взаимных временных задержек между переданными и принятыми радиосигналами для НРТС и каждой из ИОРС;
рассчитывают суммы дальностей от КА до НРТС и каждой из ИОРС;
вычисляют координаты КА используя указанные суммы дальностей.
Недостатками способа прототипа [3] являются:
длительное время определения координат КА, обусловленное необходимостью размещать на позициях с известными координатами не менее двух ИОРС;
относительно-невысокая точность определения координат КА обусловленная необходимостью синхронного излучения тестовых радиосигналов с помощью НРТС и не менее двух ИОРС.
Целью заявляемых технических решений является сокращение временных затрат на определение координат и ортогональных составляющих вектора скорости КА и повышение точности определения координат и ортогональных составляющих вектора скорости КА.
Поставленная цель в заявленном способе определения ортогональных составляющих вектора скорости КА с использованием ЗС достигается тем, что в известном способе определения ортогональных составляющих вектора скорости КА (по патенту №2652603) включающим: размещение НРТС K на позиции с известными координатами xK, yK, zK, прием в момент времени t0 с помощью НРТС K радиосигналов, переданных опорными реперными станциями (ОРС) и ретранслированных основным КА S1, вычисление координат х1, у1, z1 основного КА S1 в момент времени t0, расчет ортогональных составляющих вектора скорости основного КА в момент времени t0 на основе частотных сдвигов радиосигналов, известных координат НРСТ и ОРС, предварительно заданной частоте сдвига рабочей частоты основного 
а также вычисленных координат основного КА x1, y1, z1, дополнительно выбирают смежный КА S2 с известными значениями координат х2, y2, z2, и ортогональных составляющих вектора скорости 
в момент времени t0 и заданной частотой сдвига 
рабочей частоты, а в качестве ОРС используют выбранные ЗС In, размещенные на земной поверхности на позициях с известными координатами 
где n=1…N-номер 3C, N≥3.
Вычисляют координаты основного КА x1, y1, z1.
Для каждой n-й ЗС In измеряют в НРТС K за счет корреляционной обработки радиосигналов значения номиналов частот 
и 
принятых радиосигналов после их ретрансляции основным КА S1 и смежным КА S2 соответственно. На основе известных координат НРТС xK, yK, zK, координат x2, y2, z2 и ортогональных составляющих вектора скорости 
смежного КА S2 в момент времени t0, координат не менее трех ЗС 
рассчитывают значения радиальных скоростей 
смежного КА S2 относительно каждой из n-й ЗС In и НРТС K соответственно.
Используя измеренные значения номиналов частот 
принятых радиосигналов n-х ЗС после их ретрансляции смежным КА S2, рассчитанные доплеровские сдвиги частот на входе 
и на выходе 
смежного КА S2 за счет его сближения или удаления с или от n-й ЗС In и НРТС K, заданную частоту сдвига рабочей частоты смежного КА 
вычисляют значения номиналов частот излучаемых каждой n-й ЗС fn.
На основе известных координат НРТС xK, yK, zK, координат ЗС 
и вычисленных координат основного КА х1, y1, z1 в момент времени t0 рассчитывают расстояния 
от и n-х ЗС In и НРТС K до основного КА S1.
Вычисляют ортогональные составляющие вектора скорости основного КА 
по известным координатам НРТС xK, yK, zK, координатам ЗС 
и координатам смежного КА х2, у2, z2, известным ортогональным составляющим вектора скорости смежного КА 
заданной частоте сдвига рабочей частоты основного КА 
рассчитанным значениям номиналов частот излучаемых каждой n-й ЗС fn, расстояниям 
от и n-х ЗС In и НРТС K до основного КА.
Смежный КА S2 выбирают так, чтобы его диапазон частот на линии "вверх" 
имел одинаковые участки с диапазоном частот на линии "вверх" основного КА 
а зона покрытия смежного КА Ω2 пересекалась с зоной покрытия основного КА Ω1.
ЗС In выбирают такие, значения номиналов излучаемых частот fn которых входят в диапазоны частот на линии "вверх" основного КА 
и смежного КА 
Каждая ЗС In находится в зонах покрытия как основного КА Ω1, так и смежного КА Ω2, при этом взаимные расстояния 
между n-й и m-й ЗС, где m=1…N, m≠n максимальны.
Благодаря перечисленной новой совокупности существенных признаков, за счет использования НРТС K на позиции с известными координатами xK, yK, zK, смежного КА S2 с известными значениями координат х2, y2, z2, и ортогональных составляющих вектора скорости 
в момент времени t0, а также ЗС In, размещенных на земной поверхности на позициях с известными координатами 
достигается цель изобретения: снижение времени определения ортогональных составляющих вектора скорости КА, а также повышение точности определения ортогональных составляющих вектора скорости КА.
Поставленная цель в заявленном способе определения координат КА с использованием ЗС достигается тем, что в известном способе определения координат КА (по патенту №2652603) включающим: размещение НРТС K на позиции с известными координатами xK, yK, zK, выбор начальных значений параметров орбиты основного КА S1, прием в момент времени t0 с помощью НРТС K радиосигналов, переданных опорными реперными станциями (ОРС) и ретранслированных основным КА S1, вычисление на основе временных задержек радиосигналов координат основного КА x1, y1, z1 в момент времени t0, дополнительно выбирают смежный КА S2 с известными координатами х2, у2, z2 в момент времени t0, а в качестве ОРС используют выбранные ЗС In,
размещенные на земной поверхности на позициях с известными координатами 
где n=1…N-номер ЗС, N≥3.
Для каждой n-й ЗС In измеряют в НРТС за счет корреляционной обработки радиосигналов значения временной задержки Δtn между принятыми радиосигналами после их ретрансляции основным и смежным КА соответственно.
Рассчитывают N значений длин 
траекторий 
на основе известных координат ЗС 
смежного КА х2, y2, z2 и НРТС xK, yK, zK.
Используя длины 
траекторий 
а также измеренные временные задержки Δtn рассчитывают длины 
траекторий 
Вычисляют координаты основного КА х1, у1, z1 на основе рассчитанных длин 
траекторий InS1K, известных координат НРТС xK, yK, zK, ЗС 
и смежного КА х2, y2, z2.
Благодаря перечисленной новой совокупности существенных признаков, за счет использования НРТС K на позиции с известными координатами xK, yK, zK, смежного КА S2 с известными координатами х2, у2, z2, в момент времени t0, а также ЗС In, размещенных на земной поверхности на позициях с известными координатами где n-1…N - номер ЗС, N≥3 достигается цель изобретения: снижение времени определения координат КА, а также повышение точности определения координат КА.
Заявленные изобретения поясняются чертежами, на которых показаны: на фиг. 1 - структурная схема подсистемы определения ортогональных составляющих вектора скорости основного КА с использованием ЗС;
на фиг. 2 - схема алгоритма расчета ортогональных составляющих вектора скорости основного КА с использованием трех ЗС;
на фиг. 3 - структурная схема подсистемы определения координат основного КА с использованием ЗС;
на фиг. 4 - схема алгоритма расчета координат основного КА с использованием трех ЗС;
на фиг. 5 - схема выбора ЗС с учетом диапазонов частот на линии "вверх" основного КА и смежного КА, а также номиналов излучаемых частот ЗС;
на фиг. 6 - схема выбора ЗС с учетом зон покрытия основного и смежного КА;
на фиг. 7 - пример схемы топологии размещения трех ЗС.
Теория полета КА, или, как ее еще называют, астродинамика, небесная механика, космическая баллистика, основана на законах И. Кеплера и законе всемирного тяготения И. Ньютона.
В первом приближении движение КА представляется как невозмущенное - такое движение, которое происходило бы только под влиянием силы притяжения Земли по закону Ньютона, т.е. точно соответствует задаче двух тел (Земля - КА) в небесной механике. Это движение называется движением по Кеплеровой орбите, так как подчиняется трем законам Кеплера [4].
Достоинством Кеплеровой орбиты является простота вычисления координат и вектора скорости КА в прогнозируемый момент времени. Это предопределило широкое использование элементов Кеплеровой орбиты. В настоящем изобретении элементы Кеплеровой орбиты КА выступают в качестве априорных данных об КА. С помощью этих элементов ориентируют приемные антенны НРСТ на КА. Кроме того элементы Кеплеровой орбиты КА служат для устранения двузначности определения координат КА.
Недостатком Кеплеровой орбиты является относительно низкая точность определения координат и вектора скорости КА, которая недостаточна для выполнения ряда прикладных задач, например, для корректировки движения КА.
Более точно движение КА описывается с помощью возмущенной орбиты [4, 5], которая задается:
каноническими параметрами КА, включающими координаты КА х1, у1, z1 и ортогональные составляющие вектора его скорости 
в начальный момент времени t0,
факторами, приводящими к отклонениям КА от идеальной (Кеплеровой) орбиты и называемыми возмущающими факторами.
Например, для геостационарного КА достаточно учитывать только три фактора, приводящих к отклонениям КА от идеальной (Кеплеровой) орбиты - влияние Солнца, Луны и нецентральности гравитационного поля Земли.
При прогнозе координат КА размещенных на других (не геостационарных) орбитах дополнительно учитывают ряд факторов, приводящих к отклонениям КА от идеальной (Кеплеровой) орбиты. В качестве таких факторов, например, для КА на низких орбитах выступают: влияние сопротивления атмосферы Земли, светового давления, притяжения планет и др.
Таким образом, определение координат КА х1, у1, z1 и ортогональных составляющих вектора его скорости 
в начальный момент времени t0 с высокой точностью является важной задачей, которая решена в заявленных технических решениях.
Подсистема определения ортогональных составляющих вектора скорости основного КА 
в момент времени t0 содержит НРТС K, смежный КА S2 с известными координатами х2, у2, z2, и ортогональными составляющими вектора скорости 
в момент времени t0, а также N≥3 выбранных ЗС In, n=1…N (см. фиг. 1) на позициях с известными координатами 
излучающие радиосигналы в направлении основного и смежного КА.
На фиг. 1 введены следующие обозначения: 
- расстояние между основным КА S1 и первой ЗС I1; 
- расстояние между основным КА S1 и НРСТ K; 
- расстояние между смежным КА S2 и первой ЗС I1; 
- расстояние между смежным КА S2 и НРСТ K, 
- радиальная скорость основного КА S1 в направлении первой ЗС I1, 
- радиальная скорость основного КА S1 в направлении НРСТ K, 
- радиальная скорость смежного КА S2 в направлении первой ЗС I1, 
- радиальная скорость смежного КА S2 в направлении НРСТ K, 
- векторы скорости основного и смежного КА соответственно, 
- углы между направлением на НРТС и векторами 
соответственно; 
- углы между направлениями на первую ЗС и векторами 
соответственно.
Для каждой n-й ЗС возможно ввести обозначения: 
- радиальная скорость основного КА S1 в направлении n-ю ЗС In, 
- радиальная скорость смежного КА S2 в направлении n-ю ЗС In, 
- углы между направлениями на n-ю ЗС и векторами 
соответственно.
На фиг. 1 представлено три ЗС, как минимально необходимое количество ЗС для однозначного одномоментного определения ортогональных составляющих вектора скорости основного КА 
по предлагаемому способу.
Основополагающей предпосылкой настоящего изобретения является наличие помимо основного КА S1, через который организуется канал связи между земными станциями, смежного КА S2, который способен ретранслировать те же самые радиоизлучения что и основной, но с большим ослаблением и другой частотой переноса. Таким образом, возможно получение за счет корреляционной обработки радиосигналов значения временных задержек Δtn между радиосигналами принятыми от основного S1 и смежного КА S2, а также значения номиналов частот 
принятых радиосигналов после их ретрансляции основным КА S1 и смежным КА S2 соответственно для каждой из выбранных ЗС In.
Для определения ортогональных составляющих вектора скорости основного КА 
в момент времени t0 используют значения номиналов частот 
принятых радиосигналов после их ретрансляции основным КА S1 и смежным КА S2 соответственно, обусловленные различными радиальными скоростями основного и смежного КА относительно каждой из п-й ЗС In и НРСТ K [6].
На основе известных координат НРТС xK, yK, zK, координат х2, у2, z2 и ортогональных составляющих вектора скорости 
смежного КА S2 в момент времени t0, координат ЗС 
рассчитывают значения радиальных скоростей 
смежного КА S2 относительно каждой из n-й ЗС In и НРТС K соответственно.
Используя значения радиальных скоростей 
смежного КА S2 относительно каждой из n-й ЗС In и НРТС K, значения номиналов частот 
принятых радиосигналов после их ретрансляции смежным КА S2, значение заданной частоты сдвига рабочей частоты смежного КА 
вычисляют значения номиналов частот излучаемых каждой n-й ЗС fn.
Вычисляют ортогональные составляющие вектора скорости основного КА 
по известным координатам НРТС xK, yK, zK, координатам ЗС 
и координатам смежного КА х2, у2, z2, известным ортогональным составляющим вектора скорости смежного КА 
заданной частоте сдвига рабочей частоты основного КА 
рассчитанным значениям номиналов частот излучаемых каждой n-й ЗС fn, расстояниям 
от и n-х ЗС In и НРТС K до основного КА.
Для одномоментного и однозначного определения ортогональных составляющих вектора скорости основного КА необходимо использование трех ЗС. Дальнейшее увеличение количества ЗС будет приводить к повышению точности определения ортогональных составляющих вектора скорости основного КА.
В качестве примера в приложении А представлена алгоритм определения ортогональных составляющих вектора скорости основного КА с использованием трех ЗС. Выходными результатами представленного алгоритма выступают составляющие вектора скорости основного КА 
в момент времени t0. Подсистема определения координат основного КА х1 у1, z1 в момент времени t0 содержит НРТС K, смежный КА S2 с известными координатами х2, у2, z2 в момент времени t0 и N≥3 выбранных ЗС In, n=1…N (см. фиг. 3) на позициях с известными координатами 
излучающие радиосигналы в направлении основного S1 и смежного КА S2.
На фиг. 3 введены следующие обозначения: 
- расстояние между основным КА S1 и первой ЗС I1; 
- расстояние между основным КА S1 и НРСТ K; 
- расстояние между смежным КА S2 и первой ЗС I1; 
- расстояние между смежным КА S2 и НРСТ K.
Для каждой n-й ЗС возможно ввести обозначения: 
- расстояние между основным КА S1 и n-й ЗС In; 
- расстояние между смежным КА S2 и n-й ЗС In.
На фиг. 3 представлено три ЗС, как минимально необходимое количество ЗС для однозначного одномоментного определения координат основного КА х1, у1, z1 по предлагаемому способу.
Основополагающей предпосылкой настоящего изобретения является наличие помимо основного КА S1, через который организуется канал связи между земными станциями, смежного КА S2, который способен ретранслировать те же самые радиоизлучения что и основной, но с большим ослаблением и другой частотой переноса. Таким образом, возможно получение за счет корреляционной обработки радиосигналов значения временных задержек Δtn между радиосигналами принятыми от основного S1 и смежного КА S2 для каждой из выбранных ЗС In.
Для определения координат основного КА x1, y1, z1 в момент времени t0 используют временные задержки Δtn, обусловленные разностью длин 
и 
траекторий 
[6] для каждой n-й ЗС.
На основе измеренных временных задержек Δtn и вычисленных длин 
траекторий InS2K рассчитывают длины 
траекторий InS1K.
Каждой из длин 
траекторий InS1K ставят в соответствие поверхность положения (ПП) которая является поверхностью второго порядка - эллипсоидом вращения с фокусами, совпадающими с НРТС K и ЗС In соответственно.
Координаты точки пересечения не менее трех ПП соответствуют искомым координатам КА x1, у1, z1 в момент времени t0. Таким образом, для однозначного одномоментного определения координат основного КА х1, у1, z1 необходимо наличие не менее трех ЗС In. Дальнейшее увеличение количества ЗС будет приводить к повышению точности определения координат КА x1, у1, z1.
В качестве примера в приложении Б представлен алгоритм определения координат основного КА х1, у1 и z1 с использованием трех ЗС F1…F11. Выходными результатами представленного алгоритма выступают координаты основного КА х1, у1, z1 в момент времени t0.
На фиг. 5 в качестве примера представлена схема выбора ЗС с учетом диапазонов частот на линии "вверх" основного КА 
и смежного КА 
а также диапазонов частот ЗС F5…F12, расположенных в районах зон покрытия основного КА 
и смежного КА 
По оси абсцисс на фиг. 5 отложены частоты f, по оси ординат - амплитуды А.
Анализ схемы (фиг. 5) показывает, что в диапазон частот на линии "вверх" основного КА 
входят диапазоны частот ЗС F1…F11, а в диапазон частот на линии "вверх" смежного КА 
входят диапазоны частот F5…F12 ЗС. Таким образом, одновременно в диапазоны частот на линии "вверх" основного КА 
и смежного КА 
входят диапазоны частот ЗС F5…F11. На схеме (фиг. 5) введены обозначения указанных ЗС - I1…I7, а также значений номиналов их излучаемых частот - f1…f7.
На фиг. 6 в качестве примера представлена схема выбора ЗС с учетом зон покрытия основного Ω1 и смежного КА Ω2.
Анализ схемы (фиг. 6) показывает, что в зоне покрытия основного КА Ω1 находятся ЗС I1, I2, I3, I4 и I6, а в зоне покрытия смежного КА Q.2 находятся ЗС I1, I2, I3, I4 и I7. Таким образом, одновременно в зонах покрытия основного Ω1 и смежного КА Ω2 находятся ЗС I1, I2, I3 и I4.
При выборе ЗС учитывают топологию размещения всех ЗС In, от которой зависит точность определения координат х1, у1, z1 и ортогональных составляющих вектора скорости КА 
Под топологией размещения ЗС In понимают их взаимное расположение на поверхности Земли. При этом в качестве важных показателей при выборе ЗС выступают взаимные расстояния 
между n-й и m-й ЗС, где m=I…N, m≠n, которые должны быть максимальными.
Сопоставляя фиг. 5 и фиг. 6 для определения координат х1, у1, z1 и ортогональных составляющих вектора скорости КА 
выбирают ЗС I1, I2 и I3.
На фиг. 7 в качестве примера представлена схемы топологии размещения трех ЗС I1, I2 и I3. Указанные ЗС выбирались так, что бы взаимные расстояния 
были максимальными.
Имитационное моделирование на основе программ для ЭВМ [7, 8] заявленных способов показало возможность снижения времени определения координат х1, у1, z1 и ортогональных составляющих вектора скорости КА 
в 5…10 раз, а также повышения точности определения координат х1, у1, z1 и ортогональных составляющих вектора скорости КА 
по сравнению со способом прототипом на 10…20%.
Источники информации
1. Агиевич С.Н., Беспалов В.Л., Дедовская Э.Г., Матюхин А.С., Подъячев П.А., Севидов В.В. Способ определения параметров орбиты искусственного спутника Земли с использованием приемных опорных реперных станций. Патент №2702098 МПК G01S 5/00 (2006.01). Бюл. №28 от 04.10.19. Заявка №2018127491 от 25.07.18.
2. Агиевич С.Н., Ватутин В.М., Матюхин А.С., Модин М.И., Севидов В.В. Способ определения параметров орбиты искусственного спутника Земли с использованием приемо-передающих опорных реперных станций. Патент №2708883. МПК G01S 5/00 (2006.01). Бюл. №35 от 12.12.19. Заявка №2018134855 от 01.10.18.
3. Балабанов В.В., Беспалов В.Л., Кельян А.Х., Пономарев А.А., Севидов В.В., Чемаров А.О. Способ определения параметров орбиты искусственного спутника Земли. Патент на изобретение №2652603, опубл. 27.04.2018 Бюл. №12.
4. Машбиц Л.М. Компьютерная картография и зоны спутниковой связи, - 2-е изд., перераб и доп. - М.: Горячая линия - Телеком, 2009. - 236 с.
5. Волков Р.В., Малышев С.Р., Симонов А.Н., Севидов В.В. Определение канонических параметров спутников-ретрансляторов по радиосигналам опорных реперных станций // Труды Военно-космической академии им. А.Ф. Можайского. 2016. Вып. 655. С. 88-92.
6. Волков Р.В., Саяпин В.Н., Севидов В.В. Модель измерения временной задержки и частотного сдвига радиосигнала, принятого от спутника-ретранслятора при определении местоположения земной станции // T-Comm: Телекоммуникации и транспорт. 2016. Том 10. №9. С. 14-18.
7. Волков Р.В., Саяпин В.Н., Севидов В.В. Модель движения искусственного спутника Земли // Программы для ЭВМ. Базы данных. Топологии интегральных микросхем. 2016. №2. С. 112.
8. Севидов В.В. Определение координат и параметров движения источника радиоизлучения на основе разностно-временных и разностно-доплеровских измерений // Программы для ЭВМ. Базы данных. Топологии интегральных микросхем. 2015. №11. С. 2.
Приложение А
Алгоритм определения ортогональных составляющих вектора скорости основного КА с использованием трех ЗС
Для определения ортогональных составляющих вектора скорости основного КА 
возможно использовать НРТС K, смежный КА S2 с известными координатами х2, у2, z2, и ортогональными составляющими вектора скорости 
в момент времени t0, а также не менее трех выбранных ЗС In, (см. фиг. 1), размещенных на позициях с известными координатами 
излучающие радиосигналы в направлении основного и смежного КА.
В качестве примера в настоящем приложении представлен вариант с тремя ЗС (n=1…3) как минимально необходимого количества ЗС для однозначного одномоментного определения вектора скорости основного КА 
Предполагается, что координаты основного КА x1, y1, z1 - рассчитаны в соответствии с алгоритмом, представленном в приложении А.
С помощью корреляционной обработки радиосигналов в НРТС K измеряют значения номиналов частот 
принятых радиосигналов n-х ЗС In после их ретрансляции основным и смежным КА.
Для номиналов частот 
справедливы аналитические выражения:
где fn - значения номиналов частот излучаемых каждой n-й ЗС In;
- доплеровские сдвиги частот на входе основного и смежного КА соответственно за счет его сближения (удаления) с (от) n-й ЗС In;
- заданные частоты сдвига рабочих частот основного и смежного КА соответственно; 
- доплеровские сдвиги частот на выходе основного и смежного КА соответственно за счет его сближения (удаления) с (от) НРТС K.
Предполагают, что нестабильность генератора частот КА известна и компенсируется. Влияние других эффектов на изменение частоты, например, гравитационный и релятивистский эффекты в рамках рассматриваемой задачи, пренебрежимо мало и поэтому не учитывают.
Для расчета ортогональных составляющих вектора скорости основного КА 
с использованием трех ЗС разработан алгоритм, схема которого представлена на фиг. 2.
На этапе 1 производят ввод исходных данных, в качестве которых выступают: момент времени измерения t0, координаты НРТС xK, yK, zK; координаты основного КА х1, у1 и z1; координаты х2, y2, z2 и ортогональными составляющими вектора скорости 
смежного КА S2 в момент времени t0, координаты трех ЗС 
значения частот сдвига рабочих частот 
основного и смежного КА соответственно; значения номиналов частот 
принятых радиосигналов n-х ЗС In после их ретрансляции основным и смежным КА соответственно.
На этапе 2 рассчитывают значения радиальных скоростей 
смежного КА S2 относительно каждой из n-й ЗС In и НРТС К соответственно. Для значений радиальных скоростей 
смежного КА S2 возможно записать формулы:
Согласно теореме о скалярном произведении векторов [6] справедливы равенства:
Модуль вектора скорости смежного КА 
равен:
а расстояния от и n-х ЗС и НРТС до смежного КА рассчитывают, как
Уравнения (Б.3) и (Б.4) с учетом (Б.5)…(Б.9) преобразуют к виду:
На этапе 3 рассчитывают значения номиналов частот fn излучаемых каждой n-й ЗС In.
Доплеровские сдвиги частот на входе 
и на выходе 
смежного КА S2 за счет его сближения (удаления) с (от) n-й ЗС In и НРТС K возможно представить в виде:
Для расчета значений номиналов частот излучаемых каждой из ЗС fn выражения (А.2) с учетом уравнений (А.10) и (А.11) преобразуют к виду:
На этапе 4 рассчитывают расстояния 
от и n-х ЗС In и НРТС K до основного КА S1 по формулам:
На этапе 5 вычисляют ортогональные составляющие вектора скорости основного КА 
Доплеровские сдвиги частот на входе
и на выходе 
основного КА S1 за счет его сближения (удаления) с (от) n-й ЗС In и НРТС K возможно представить в виде:
Для значений радиальных скоростей 
основного КА S1 возможно записать формулы:
Согласно теореме о скалярном произведении векторов [6] справедливы равенства:
Модуль вектора скорости основного КА 
равен:
Уравнения (А.14) и (А.15) с учетом (А.12), (А.13), (А.16)…(А.20) преобразуют к виду:
Выражения (А.1) с учетом уравнений (А.12)…(А.22) для частного случая, когда n=1…3, преобразуют в систему линейных уравнений:
где коэффициенты при переменных и свободные члены равны:
Систему из трех линейных уравнений с тремя неизвестными (А.23) решают одним из известных методов, например, методом Крамера. Результатом решения системы уравнений (А.23) выступают ортогональные составляющие вектора скорости основного КА 
На этапе 6 осуществляют вывод результатов, в качестве которых выступают ортогональные составляющие вектора скорости основного КА 
в момент времени t0.
В общем случае, когда количество ЗС N>3, алгоритм определения ортогональных составляющих вектора скорости КА 
остается прежним, с той лишь разницей, что система уравнений (А.23) будет содержать более трех уравнений. Тогда такую систему уравнений решают, например, методом наименьших квадратов.
Приложение Б
Алгоритм определения координат основного КА с использованием трех ЗС
Алгоритм определения координат основного КА х1, у1, z1 с использованием ЗС In, где n-1…N - номер ЗС, N≥3 размещенных на позициях с известными координатами 
основан на том, что каждой из временных задержек Δtn соответствует разности длин 
траекторий 
В качестве примера в настоящем приложении представлен вариант с тремя ЗС (n=1…3) как минимально необходимого количества для однозначного одномоментного определения координат основного КА х1, у1, z1
Длины 
траекторий InS2K вычисляют по формулам:
Используя длины 
траекторий InS2K, а также измеренные временные задержки Δtn рассчитывают длины 
траекторий InS1K:
где с=3×108 м/с - скорость света в вакууме.
Каждой из длин 
траекторий InS1K ставят в соответствие поверхность положения (ПП) основного КА S1, которая является поверхностью второго порядка - эллипсоидом вращения с фокусами, совпадающими с НРТС K и ЗС In соответственно.
Координаты точки пересечения трех ПП соответствуют искомым координатам основного КА x1, y1, z1 в момент времени t0. Таким образом, для однозначного одномоментного определения координат основного КА x1, y1, z1 необходимо наличие не менее трех ЗС.
Для расчета координат КА x1, y1, z1 с использованием трех ЗС разработан алгоритм, схема которого представлена на фиг. 4.
На этапе 1 производят ввод исходных данных, в качестве которых выступают: момент времени измерения t0; координаты НРТС xK, yK, zK координаты смежного КА х2, у2, z2; координаты трех ЗС 
и 
временные задержки между радиосигналами принятыми от основного и смежного КА для каждой из выбранных ЗС 
порог точности δ0 расчета координат основного КА.
На этапе 2 рассчитывают длины 
траекторий InS1K по формулам (Б.1, Б.2).
На этапе 3 выбирают, на основе элементов Кеплеровой орбиты основного КА, координаты опорной точки 
как первое приближение к координатам основного КА.
На этапе 4 рассчитывают длины 
траекторий 
при условии равенства координат основного КА координатам опорной точки 
по формулам:
где 
- расстояния от опорной точки 
до n-й ЗС In, 
- расстояние от опорной точки 
до НРТС K,
Расстояния 
в свою очередь рассчитываются по формулам:
На этапе 5 рассчитывают невязки k1, k2 и k3 как разницы между определенными на этапе 4 длинами 
траекторий 
и длинами 
траекторий InS1K, определенными на этапе 2 соответственно:
На этапе 6 определяют поправки к координатам основного КА Δх1, Δу1, Δz1.
Для определения поправок к координатам основного КА Δx1, Δу1, Δz1 предварительно формируют систему линейных уравнений при разложении в ряд Тейлора функций 
с точностью до первых членов:
где частные производные рассчитываются согласно выражениям
Решая систему линейных уравнений (Б.3) одним из известных методов, например, методом Крамера, получают поправки к координатам основного КА Δх1, Δу1, Δz1.
На этапе 7 рассчитывают координаты новой опорной точки 
Этапы 4-7 в совокупности составляют первую итерацию. Далее итерации повторяют, используя каждый раз новую опорную точку, полученную на предыдущей итерации. Количество необходимых итераций зависит требуемой точности определения координат КА. С точностью определения координат КА напрямую связан шаг итерации dш.
На этапе 8 определяют шаг итерации dш как расстояние между текущей и предыдущей опорными точками:
На этапе 8 сравнивают dш с порогом δ0, задаваемом на этапе 1.
Необходимое число итераций, как правило, составляет 2…4. В качестве координат основного КА х1, у1 и z1 выбирают значения координат опорной точки на последней итерации, вывод которых осуществляют на этапе 10.
В общем случае, когда количество ЗС N>3, алгоритм определения координат основного КА остается прежним, с той лишь разницей, что система уравнений (Б.3) будет содержать более трех уравнений. Тогда такую систему уравнений решают, например, методом наименьших квадратов.
1. Способ определения ортогональных составляющих вектора скорости космического аппарата (КА) с использованием земных станций (ЗС), заключающийся в том, что размещают наземную радиотехническую станцию (НРТС) K на позиции с известными координатами xK, yK, zK, принимают в момент времени t0 с помощью НРТС K радиосигналы, переданные опорными реперными станциями (ОРС) и ретранслированные основным КА S1, вычисляют координаты х1, у1, z1 основного КА S1 в момент времени t0, на основе частотных сдвигов радиосигналов, известных координат НРСТ и ОРС, предварительно заданной частоте сдвига рабочей частоты основного КА 
рассчитывают ортогональные составляющие вектора скорости основного КА 
в момент времени t0, отличающийся тем, что дополнительно выбирают смежный КА S2 с известными значениями координат х2, у2, z2 и ортогональных составляющих вектора скорости 
в момент времени t0 и заданной частотой сдвига 
рабочей частоты, а в качестве ОРС используют выбранные ЗС In, размещенные на земной поверхности на позициях с известными координатами 
, где n=1…N - номер ЗС, N≥3, вычисляют координаты основного КА х1, у1, z1, для каждой n-й ЗС In измеряют в НРТС K за счет корреляционной обработки радиосигналов значения номиналов частот 
принятых радиосигналов после их ретрансляции основным КА S1 и смежным КА S2 соответственно, на основе известных координат НРТС xK, yK, zK, координат x2, y2, z2 и ортогональных составляющих вектора скорости 
смежного КА S2 в момент времени t0, координат не менее трех ЗС рассчитывают значения радиальных скоростей 
смежного КА S2 относительно каждой из n-й ЗС In и НРТС K соответственно, используя измеренные значения номиналов частот 
принятых радиосигналов n-х ЗС после их ретрансляции смежным КА S2, рассчитанные доплеровские сдвиги частот на входе 
и на выходе 
смежного КА S2 за счет его сближения или удаления с или от n-й ЗС In и НРТС K, заданную частоту сдвига рабочей частоты смежного КА вычисляют значения номиналов частот, излучаемых каждой n-й ЗС fn, на основе известных координат НРТС xK, yK, zK, координат ЗС 
и вычисленных координат основного КА х1, y1, z1 в момент времени t0 рассчитывают расстояния 
от и n-х ЗС In и НРТС K до основного КА S1, причем ортогональные составляющие вектора скорости основного КА 
вычисляют по известным координатам НРТС xK, yK, zK, координатам ЗС и координатам смежного КА х2, у2, z2, известным ортогональным составляющим вектора скорости смежного КА 
заданной частоте сдвига рабочей частоты основного КА 
рассчитанным значениям номиналов частот излучаемых каждой n-й ЗС fn, расстояниям 
от и n-х ЗС In и НРТС К до основного КА.
2. Способ по п. 1, отличающийся тем, что смежный КА S2 выбирают такой, чтобы его диапазон частот на линии "вверх" 
имел одинаковые участки с диапазоном частот на линии "вверх" основного КА 
а зона покрытия смежного КА Ω2 пересекается с зоной покрытия основного КА Ω1.
3. Способ по п. 1, отличающийся тем, что выбирают такие ЗС In, значения номиналов излучаемых частот fn которых входят в диапазоны частот на линии "вверх" основного КА 
и смежного КА 
а каждая ЗС In должна находится в зонах покрытия как основного КА Ω1, так и смежного КА Ω2, при этом взаимные расстояния 
между n-й и m-й ЗС, где m=1…N, m≠n максимальны.
4. Способ определения координат космического аппарата (КА) с использованием земных станций (ЗС), заключающийся в том, что размещают наземную радиотехническую станцию (НРТС) K на позиции с известными координатами xK, yK, zK, выбирают начальные значения параметров орбиты основного КА S1, принимают в момент времени t0 с помощью НРТС K радиосигналы, переданные опорными реперными станциями (ОРС) и ретранслированные основным КА S1, на основе временных задержек радиосигналов системы вычисляют координаты основного КА х1, у1, z1 в момент времени t0, отличающийся тем, что дополнительно выбирают смежный КА S2 с известными координатами х2, y2, z2 в момент времени t0, а в качестве ОРС используют выбранные земные станции (ЗС) In, размещенные на земной поверхности на позициях с известными координатами 
где n=1…N - номер ЗС, N≥3, для каждой n-й ЗС In измеряют в НРТС за счет корреляционной обработки радиосигналов значения временной задержки Δtn между принятыми радиосигналами после их ретрансляции основным и смежным КА соответственно, рассчитывают N значений длин 
траекторий InS2K на основе известных координат ЗС смежного КА x2, y2, z2 и НРТС xK, yK, zK, используя длины , траекторий InS2K, а также измеренные временные задержки Δtn, рассчитывают длины 
траекторий InS1K, вычисляют координаты основного КА x1, y1, z1 на основе рассчитанных длин траекторий InS1K, известных координат НРТС xK, yK, zK, ЗС 
и смежного КА х2, у2, z2.
