RU2749878C1 - Способ определения ортогональных составляющих векторов скоростей и способ определения координат двух космических аппаратов с использованием земных станций - Google Patents

Способ определения ортогональных составляющих векторов скоростей и способ определения координат двух космических аппаратов с использованием земных станций Download PDF

Info

Publication number
RU2749878C1
RU2749878C1 RU2020136690A RU2020136690A RU2749878C1 RU 2749878 C1 RU2749878 C1 RU 2749878C1 RU 2020136690 A RU2020136690 A RU 2020136690A RU 2020136690 A RU2020136690 A RU 2020136690A RU 2749878 C1 RU2749878 C1 RU 2749878C1
Authority
RU
Russia
Prior art keywords
spacecraft
main
coordinates
adjacent
nrts
Prior art date
Application number
RU2020136690A
Other languages
English (en)
Inventor
Владимир Витальевич Севидов
Original Assignee
федеральное государственное казенное военное образовательное учреждение высшего образования "Военная академия связи имени Маршала Советского Союза С.М. Буденного" Министерства обороны Российской Федерации
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by федеральное государственное казенное военное образовательное учреждение высшего образования "Военная академия связи имени Маршала Советского Союза С.М. Буденного" Министерства обороны Российской Федерации filed Critical федеральное государственное казенное военное образовательное учреждение высшего образования "Военная академия связи имени Маршала Советского Союза С.М. Буденного" Министерства обороны Российской Федерации
Priority to RU2020136690A priority Critical patent/RU2749878C1/ru
Application granted granted Critical
Publication of RU2749878C1 publication Critical patent/RU2749878C1/ru

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G3/00Observing or tracking cosmonautic vehicles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/24Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/50Systems of measurement based on relative movement of target
    • G01S13/58Velocity or trajectory determination systems; Sense-of-movement determination systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/03Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S3/00Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
    • G01S3/02Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
    • G01S3/14Systems for determining direction or deviation from predetermined direction
    • G01S3/46Systems for determining direction or deviation from predetermined direction using antennas spaced apart and measuring phase or time difference between signals therefrom, i.e. path-difference systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S5/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Astronomy & Astrophysics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

Предложенная группа изобретений относится к области космонавтики, а именно к технике выполнения траекторных измерений и определения координат и ортогональных составляющих векторов скоростей КА, и может быть использована на наземных и бортовых комплексах управления полетом КА для точного определения текущих параметров движения КА. Техническим результатом изобретения является повышение точности определения координат и ортогональных составляющих векторов скорости двух КА. Способ определения векторов скорости основного КА
Figure 00000254
и смежного КА
Figure 00000255
включает: измерение в наземной радиотехнической станции (НРТС) K значений номиналов частот
Figure 00000256
и
Figure 00000257
принятых радиосигналов n-х ЗС In после их ретрансляции основным S1 и смежным КА S2 соответственно, расчет координат основного х1, у1, z1 и смежного КА х2, у2, z2, вычисление ортогональных составляющих вектора скорости основного КА
Figure 00000258
Способ определения координат основного х1, у1, z1 и смежного КА х2, у2, z2 включает: измерение в НРТС K для каждой n-й ЗС In значений временных задержек Δtn между принятыми радиосигналами после их ретрансляции основным S1 и смежным КА S2 соответственно, расчет координат основного КА x1, y1, z1. 2 н. и 2 з.п. ф-лы, 7 ил.

Description

Заявленные объекты объединены единым изобретательским замыслом и заявленный способ определения координат двух космических аппаратов (КА) с использованием земных станций (ЗС) предназначен для реализации способа определения ортогональных составляющих векторов скорости двух КА с использованием ЗС.
Изобретения относятся к области космонавтики, а именно к технике выполнения траекторных измерений и определения параметров орбиты КА, и могут быть использованы на наземных и бортовых комплексах управления полетом КА для точного определения текущих параметров движения КА.
Известен способ определения ортогональных составляющих вектора скорости КА [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, на основе частотных сдвигов радиосигналов, известных координат НРСТ, ОРС и предварительно заданной частоте сдвига рабочей частоты основного КА
Figure 00000001
рассчитывают ортогональные составляющие вектора скорости основного КА
Figure 00000002
в момент времени t0, дополнительно выбирают смежный КА S2 с неизвестными значениями координат х2, y2, z2, неизвестными ортогональными составляющими вектора скорости
Figure 00000003
в момент времени t0 и заданной частотой сдвига
Figure 00000004
рабочей частоты, а в качестве ОРС используют выбранные ЗС In, размещенные на земной поверхности на позициях с известными координатами xIn, yIn, zIn, где n=1…N - номер ЗС, N≥6, излучающие радиосигналы на известных частотах ƒn в направлении основного и смежного КА.
Вычисляют координаты основного x1, у1, z1 и смежного КА х2, у2, z2 в момент времени t0.
Для каждой n-й ЗС In измеряют в НРТС K за счет корреляционной обработки радиосигналов значения номиналов частот
Figure 00000005
и
Figure 00000006
принятых радиосигналов после их ретрансляции основным КА S1 и смежным КА S2 соответственно.
На основе известных координат НРТС xK, yK, zK, координат ЗС xIn, yIn, zIn и вычисленных координат основного КА x1, у1, z1 в момент времени t0 рассчитывают расстояния
Figure 00000007
и
Figure 00000008
от и n-х ЗС In и НРТС K до основного КА S1.
На основе известных координат НРТС xK, yK, zK, координат ЗС xIn, yIn, zIn и вычисленных координат смежного КА х2, у2, z2 в момент времени t0 рассчитывают расстояния
Figure 00000009
и
Figure 00000010
от и n-х ЗС In и НРТС K до смежного КА S2.
Используя значения номиналов частот
Figure 00000011
и
Figure 00000012
принятых радиосигналов после их ретрансляции основным КА S1 и смежным КА S2 вычисляют значения разностей частот Δƒn, между принятыми радиосигналами от каждой из n-й ЗС In после их ретрансляции основным и смежным КА соответственно.
Вычисляют ортогональные составляющие векторов скоростей основного
Figure 00000013
и смежного КА
Figure 00000014
по известным координатам НРТС xK, yK, zK, координатам ЗС xIn, yIn, zIn и координатам основного х1, y1, z1 и смежного КА х2, у2, z2, заданным частотам сдвига рабочих частот основного
Figure 00000015
и смежного КА
Figure 00000016
известным значениям номиналов частот излучаемых каждой n-й ЗС ƒn, рассчитанным расстояниям
Figure 00000017
и
Figure 00000018
от и n-х ЗС In и НРТС K до основного КА, рассчитанным расстояниям
Figure 00000019
и
Figure 00000020
от и n-х ЗС In и НРТС K до смежного КА S2.
Смежный КА S2 выбирают так, чтобы его диапазон частот на линии "вверх"
Figure 00000021
имел одинаковые участки с диапазоном частот на линии "вверх" основного КА
Figure 00000022
а зона покрытия смежного КА Ω2 пересекалась с зоной покрытия основного КА Ω1.
ЗС In выбирают такие, значения номиналов излучаемых частот ƒn которых входят в диапазоны частот на линии "вверх" основного КА
Figure 00000023
и смежного КА
Figure 00000024
Каждая ЗС In должна находиться в зонах покрытия как основного КА Ω1, так и смежного КА Ω2. Взаимные расстояния
Figure 00000025
между n-й и m-й ЗС, где m=1…N, m≠n должны быть максимальными.
Благодаря перечисленной новой совокупности существенных признаков, за счет использования НРТС K на позиции с известными координатами xK, yK, zK, смежного КА S2 с неизвестными значениями координат х2, у2, z2, и неизвестными ортогональными составляющими вектора скорости
Figure 00000026
в момент времени t0, предварительно заданных частот сдвига рабочих частот основного
Figure 00000027
и смежного КА
Figure 00000028
а также ЗС In, размещенных на земной поверхности на позициях с известными координатами xIn, yIn, zIn, достигается цель изобретения: снижение времени определения ортогональных составляющих векторов скорости двух КА, а также повышение точности определения ортогональных составляющих векторов скорости двух КА.
Поставленная цель в заявленном способе определения координат двух КА с использованием земных станций достигается тем, что в известном способе определения координат КА с использованием земных станций (по патенту №2652603) включающим: размещение НРТС K на позиции с известными координатами xK, yK, zK, выбор начальных значений параметров орбиты основного КА S1, прием в момент времени t0 с помощью НРТС K радиосигналов, переданных опорными реперными станциями (ОРС) и ретранслированных основным КА S1, вычисление на основе временных задержек радиосигналов координат основного КА х1, у1, z1 в момент времени t0, дополнительно выбирают смежный КА S2 с неизвестными координатами x2, y2, z2 в момент времени t0, а в качестве ОРС используют выбранные земные станции (ЗС) In, размещенные на земной поверхности на позициях с известными координатами xIn, yIn, zIn, где n=1…N - номер ЗС, N≥6.
Для каждой n-й ЗС In измеряют в НРТС за счет корреляционной обработки радиосигналов значения временной задержки Δtn между принятыми радиосигналами после их ретрансляции основным S1 и смежным КА S2 соответственно.
Используя задержки Δtn рассчитывают разности длин ΔRn траекторий InS1K и InS2K.
Вычисляют координаты основного х1, у1, z1 и смежного КА х2, у2, z2 на основе рассчитанных разностей длин ΔRn траекторий InS1K и InS2K, известных координат ЗС xIn, yIn, zIn и НРТС xK, yK, zK.
Благодаря перечисленной новой совокупности существенных признаков, за счет использования НРТС K на позиции с известными координатами xK, yK, zK, смежного КА S2 с неизвестными координатами х2, у2, z2, в момент времени t0, а также ЗС In, размещенных на земной поверхности на позициях с известными координатами xIn, yIn, zIn, где n=1…N - номер ЗС, N≥6 достигается цель изобретения: снижение времени определения координат двух КА, а также повышение точности определения координат двух КА.
Заявленные изобретения поясняются чертежами, на которых показаны:
на фиг. 1 - структурная схема подсистемы определения ортогональных составляющих вектора скорости основного КА с использованием ЗС;
на фиг. 2 - схема алгоритма расчета ортогональных составляющих вектора скорости основного КА с использованием шести ЗС;
на фиг. 3 - структурная схема подсистемы определения координат основного КА с использованием ЗС;
на фиг. 4 - схема алгоритма расчета координат основного КА с использованием шести ЗС;
на фиг. 5 - схема выбора ЗС с учетом диапазонов частот на линии "вверх" основного КА и смежного КА, а также номиналов излучаемых частот ЗС;
на фиг. 6 - схема выбора ЗС с учетом зон покрытия основного и смежного КА;
на фиг. 7 - пример схемы топологии размещения шести ЗС.
Теория полета КА, или, как ее еще называют, астродинамика, небесная механика, космическая баллистика, основана на законах И. Кеплера и законе всемирного тяготения И. Ньютона.
В первом приближении движение КА представляется как невозмущенное - такое движение, которое происходило бы только под влиянием силы притяжения Земли по закону Ньютона, т.е. точно соответствует задаче двух тел (Земля - КА) в небесной механике. Это движение называется движением по Кеплеровой орбите, так как подчиняется трем законам Кеплера [4].
Достоинством Кеплеровой орбиты является простота вычисления координат и вектора скорости КА в прогнозируемый момент времени. Это предопределило широкое использование элементов Кеплеровой орбиты. В настоящем изобретении элементы Кеплеровой орбиты КА выступают в качестве априорных данных об КА. С помощью этих элементов ориентируют приемные антенны НРСТ на КА. Кроме того элементы Кеплеровой орбиты КА служат для устранения двузначности определения координат КА.
Недостатком Кеплеровой орбиты является относительно низкая точность определения координат и вектора скорости КА, которая недостаточна для выполнения ряда прикладных задач, например, для корректировки движения КА.
Более точно движение КА описывается с помощью возмущенной орбиты [4, 5], которая задается:
каноническими параметрами КА, включающими координаты КА х1, у1, z1 и ортогональные составляющие вектора его скорости
Figure 00000029
в начальный момент времени t0,
факторами, приводящими к отклонениям КА от идеальной (Кеплеровой) орбиты и называемыми возмущающими факторами.
Например, для геостационарного КА достаточно учитывать только три фактора, приводящих к отклонениям КА от идеальной (Кеплеровой) орбиты - влияние Солнца, Луны и нецентральности гравитационного поля Земли.
При прогнозе координат КА размещенных на других (не геостационарных) орбитах дополнительно учитывают ряд факторов, приводящих к отклонениям КА от идеальной (Кеплеровой) орбиты. В качестве таких факторов, например, для КА на низких орбитах выступают: влияние сопротивления атмосферы Земли, светового давления, притяжения планет и др.
Таким образом, определение координат КА х1, у1, z1 и ортогональных составляющих вектора его скорости
Figure 00000030
в начальный момент времени t0 с высокой точностью является важной задачей, которая решена в заявленных технических решениях.
Подсистема определения ортогональных составляющих векторов скоростей основного
Figure 00000031
и смежного КА
Figure 00000032
в момент времени t0 содержит НРТС K, а также N≥6 выбранных ЗС In, n=1…N (см. фиг. 1) на позициях с известными координатами xIn, yIn, zIn, излучающие радиосигналы в направлении основного и смежного КА.
На фиг. 1 введены следующие обозначения:
Figure 00000033
- расстояние между основным КА S1 и первой ЗС I1,
Figure 00000034
- расстояние между основным КА S1 и НРСТ K;
Figure 00000035
- расстояние между смежным КА S2 и первой ЗС I1,
Figure 00000036
- расстояние между смежным КА S2 и НРСТ K;
Figure 00000037
- радиальная скорость основного КА S1 в направлении первой ЗС I1,
Figure 00000038
- радиальная скорость основного КА S1 в направлении НРСТ K,
Figure 00000039
- радиальная скорость смежного КА S2 в направлении первой ЗС I1,
Figure 00000040
- радиальная скорость смежного КА S2 в направлении НРСТ K,
Figure 00000041
и
Figure 00000042
- векторы скорости основного и смежного КА соответственно,
Figure 00000043
и
Figure 00000044
- углы между направлением на НРТС и векторами
Figure 00000045
и
Figure 00000046
соответственно;
Figure 00000047
и
Figure 00000048
- углы между направлениями на первую ЗС и векторами
Figure 00000049
и
Figure 00000050
соответственно.
Для каждой n-й ЗС возможно ввести обозначения:
Figure 00000051
- расстояния между основным КА S1 и n-й ЗС In;
Figure 00000052
- расстояния между смежным КА S2 и n-й ЗС In;
Figure 00000053
- радиальные скорости основного КА S1 в направлении n-ю ЗС In,
Figure 00000054
- радиальные скорости смежного КА S2 в направлении n-ю ЗС In,
Figure 00000055
и
Figure 00000056
- углы между направлениями на n-ю ЗС и векторами
Figure 00000057
и
Figure 00000058
соответственно.
На фиг. 1 представлено шесть ЗС, как минимально необходимое количество ЗС для однозначного одномоментного определения ортогональных составляющих векторов скоростей основного
Figure 00000059
и смежного КА
Figure 00000060
в момент времени t0 по предлагаемому способу.
Основополагающей предпосылкой настоящего изобретения является наличие помимо основного КА S1, через который организуется канал связи между земными станциями, смежного КА S2, который способен ретранслировать те же самые радиоизлучения что и основной, но с большим ослаблением и другой частотой переноса. Таким образом, возможно получение за счет корреляционной обработки радиосигналов значения временных задержек Δtn между радиосигналами принятыми от основного S1 и смежного КА S2, а также значения номиналов частот
Figure 00000061
и
Figure 00000062
принятых радиосигналов после их ретрансляции основным КА S1 и смежным КА S2 соответственно для каждой из выбранных ЗС In.
Для определения ортогональных составляющих векторов скоростей основного
Figure 00000063
и смежного КА
Figure 00000064
в момент времени t0 используют значения номиналов частот
Figure 00000065
и
Figure 00000066
принятых радиосигналов после их ретрансляции основным КА S1 и смежным КА S2 соответственно, обусловленные различными радиальными скоростями основного и смежного КА относительно каждой из n-й ЗС In и НРСТ K [6].
На основе известных координат НРТС xK, yK, zK, координат ЗС xIn, yIn, zIn, координат основного КА х1, у1, z1 в момент времени t0 рассчитывают расстояния
Figure 00000067
и
Figure 00000068
от НРТС K и n-х ЗС In до основного КА S1.
На основе известных координат НРТС xK, yK, zK, координат ЗС xIn, yIn, zIn, координат смежного КА х2, у2, z2 в момент времени t0 рассчитывают расстояния
Figure 00000069
и
Figure 00000070
от НРТС K и n-х ЗС In до смежного КА S2.
Вычисляют ортогональные составляющие векторов скоростей основного
Figure 00000071
и смежного КА
Figure 00000072
в момент времени t0 по известным координатам НРТС xK, yK, zK, координатам ЗС xIn, yIn, zIn и координатам основного х1, у1, z1 и смежного КА х2, у2, z2, заданным частотам сдвига рабочих частот основного
Figure 00000073
и смежного КА
Figure 00000074
известным номиналам частот излучаемых каждой n-й ЗС ƒn, рассчитанным расстояниям
Figure 00000075
и
Figure 00000076
от и n-х ЗС In и НРТС K до основного КА, рассчитанным расстояниям
Figure 00000077
и
Figure 00000078
от НРТС K и n-х ЗС In до смежного КА S2.
Для одномоментного и однозначного определения ортогональных составляющих векторов скоростей основного
Figure 00000079
и смежного КА
Figure 00000080
в момент времени t0 необходимо использование шести ЗС. Дальнейшее увеличение количества ЗС будет приводить к повышению точности определения ортогональных составляющих вектора скорости основного КА.
В качестве примера в приложении А представлен алгоритм определения ортогональных составляющих векторов скоростей основного
Figure 00000081
и смежного КА
Figure 00000082
в момент времени t0 с использованием шести ЗС. Выходными результатами представленного алгоритма выступают ортогональные составляющие векторов скоростей основного
Figure 00000083
и смежного КА
Figure 00000084
в момент времени t0.
Подсистема определения координат основного х1, у1, z1 и смежного КА х2, y2, z2 в момент времени t0 содержит НРТС K и N≥6 выбранных ЗС In, n=1…N - номер ЗС, N≥6 (см. фиг. 3) на позициях с известными координатами xIn, yIn, zIn, излучающие радиосигналы в направлении основного S1 и смежного КА S2.
На фиг. 3 введены следующие обозначения:
Figure 00000085
- расстояние между основным КА S1 и первой ЗС I1;
Figure 00000086
- расстояние между основным КА S1 и НРСТ K;
Figure 00000087
- расстояние между смежным КА S2 и первой ЗС I1;
Figure 00000088
- расстояние между смежным КА S2 и НРСТ K.
Для каждой n-й ЗС возможно ввести обозначения:
Figure 00000089
- расстояния между основным КА S1 и n-й ЗС In;
Figure 00000090
расстояния между смежным КА S2 и n-й ЗС In.
На фиг. 3 представлено шесть ЗС, как минимально необходимое количество ЗС для однозначного одномоментного определения координат основного x1, у1, z1 и смежного КА х2, у2, z2 в момент времени t0 по предлагаемому способу.
Основополагающей предпосылкой настоящего изобретения является наличие помимо основного КА S1, через который организуется канал связи между земными станциями, смежного КА S2, который способен ретранслировать те же самые радиоизлучения что и основной, но с большим ослаблением и другой частотой переноса. Таким образом, возможно получение за счет корреляционной обработки радиосигналов значения временных задержек Δtn между радиосигналами принятыми от основного S1 и смежного КА S2 для каждой из выбранных ЗС In.
Для определения координат основного x1, у1, z1 и смежного КА х2, у2, z2 в момент времени t0 в момент времени t0 используют временные задержки Δtn, обусловленные разностью длин
Figure 00000091
и
Figure 00000092
траекторий
Figure 00000093
и
Figure 00000094
[6] для каждой n-й ЗС.
На основе измеренных временных задержек Δtn рассчитывают разности длин ΔRn траекторий
Figure 00000095
и
Figure 00000096
Используя рассчитанные разности длин ΔRn траекторий
Figure 00000097
и
Figure 00000098
а также известные координаты НРТС xK, yK, zK и ЗС
Figure 00000099
рассчитывают координаты основного х1, у1, z1 и смежного КА х2, у2, z2 в момент времени t0.
Для однозначного одномоментного определения координат основного х1, у1, z1 и смежного КА х2, у2, z2 в момент времени t0 необходимо наличие не менее шести ЗС In. Дальнейшее увеличение количества ЗС будет приводить к повышению точности определения координат основного х1, у1, z1 и смежного КА х2, у2, z2 в момент времени t0.
В качестве примера в приложении Б представлен алгоритм определения координат основного х1, у1, z1 и смежного КА х2, у2, z2 в момент времени t0 с использованием шести ЗС In. Выходными результатами представленного алгоритма выступают координаты основного x1, у1, z1 и смежного КА х2, у2, z2 в момент времени t0 в момент времени t0.
На фиг. 5 в качестве примера представлена схема выбора ЗС с учетом диапазонов частот на линии "вверх" основного КА
Figure 00000100
и смежного КА
Figure 00000101
а также диапазонов частот ЗС F1…F12, расположенных в районах зон покрытия основного КА
Figure 00000102
и смежного КА
Figure 00000103
По оси абсцисс на фиг. 5 отложены частоты ƒ, по оси ординат - амплитуды А.
Анализ схемы (фиг. 5) показывает, что в диапазон частот на линии "вверх" основного КА
Figure 00000104
входят диапазоны частот ЗС F1…F11, а в диапазон частот на линии "вверх" смежного КА
Figure 00000105
входят диапазоны частот F3…F12 ЗС. Таким образом, одновременно в диапазоны частот на линии "вверх" основного КА
Figure 00000106
и смежного КА
Figure 00000107
входят диапазоны частот ЗС F3…F11. На схеме (фиг. 5) введены обозначения указанных ЗС - I1…I9, а также значений номиналов их излучаемых частот - ƒ1…ƒ9.
На фиг. 6 в качестве примера представлена схема выбора ЗС с учетом зон покрытия основного Ω1 и смежного КА Ω2.
Анализ схемы (фиг. 6) показывает, что в зоне покрытия основного КА Ω1 находятся ЗС I1, I2, I3, I4, I5, I6, I7 и I8, а в зоне покрытия смежного КА Ω2 находятся ЗС I1, I2, I3, I4, I5, I6,I7 и I9. Таким образом, одновременно в зонах покрытия основного Ω1 и смежного КА Ω2 находятся ЗС I1, I2, I3, I4, I5, I6 и I7.
При выборе ЗС учитывают топологию размещения всех ЗС In, от которой зависит точность определения координат основного х1, у1, z1 и смежного КА х2, у2, z2 и ортогональных составляющих векторов скоростей основного
Figure 00000108
и смежного КА
Figure 00000109
Под топологией размещения ЗС In понимают их взаимное расположение на поверхности Земли. При этом в качестве важных показателей при выборе ЗС выступают взаимные расстояния
Figure 00000110
между n-й и m-й ЗС, где m=1…N, m≠n, которые должны быть максимальными.
Сопоставляя фиг. 5 и фиг. 6 для определения координат основного х1, y2, z1 и смежного КА х2, у2, z2 и ортогональных составляющих векторов скоростей основного
Figure 00000111
и смежного КА
Figure 00000112
выбирают ЗС I1, I2, I3, I4, I5 и I6.
На фиг. 7 в качестве примера представлена схемы топологии размещения шести ЗС I1, I2, I3, I4, I5 и I6. Указанные ЗС выбирались так, что бы взаимные расстояния
Figure 00000113
были максимальными.
Имитационное моделирование на основе программ для ЭВМ [7, 8] заявленных способов показало возможность снижения времени определения координат основного x1, у1, z1 и смежного КА х2, у2, z2 и ортогональных составляющих векторов скоростей основного
Figure 00000114
и смежного КА
Figure 00000115
в 2…5 раз, а также повышение точности определения координат основного х1, y1, z1 и смежного КА х2, у2, z2 и ортогональных составляющих векторов скоростей основного
Figure 00000116
и смежного КА
Figure 00000117
по сравнению со способом прототипом на 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.
Приложение А
Алгоритм определения ортогональных составляющих векторов скоростей основного и смежного КА с использованием шести ЗС
Для определения ортогональных составляющих векторов скоростей основного
Figure 00000118
и смежного КА
Figure 00000119
в момент времени t0 возможно использовать НРТС K, а также не менее шести выбранных ЗС In, (см. фиг. 1), размещенных на позициях с известными координатами
Figure 00000120
излучающие радиосигналы на известных частотах ƒn в направлении основного и смежного КА.
В качестве примера в настоящем приложении представлен вариант с шестью ЗС (n=1…6) как минимально необходимого количества ЗС для однозначного одномоментного определения векторов скоростей основного
Figure 00000121
и смежного КА
Figure 00000122
Предполагается, что координаты основного х1, у1, z1 и смежного КА х2, y2, z2 в момент времени t0 - рассчитаны в соответствии с алгоритмом, представленном в приложении А.
С помощью корреляционной обработки радиосигналов в НРТС K измеряют значения номиналов частот
Figure 00000123
и
Figure 00000124
принятых радиосигналов n-х ЗС In после их ретрансляции основным и смежным КА соответственно.
Для номиналов частот
Figure 00000125
и
Figure 00000126
справедливы аналитические выражения:
Figure 00000127
Figure 00000128
где ƒn - значения номиналов частот излучаемых каждой n-й ЗС In;
Figure 00000129
и
Figure 00000130
- доплеровские сдвиги частот на входе основного и смежного КА соответственно за счет его сближения (удаления) с (от) n-й ЗС In;
Figure 00000131
- заданные частоты сдвига рабочих частот основного и смежного КА соответственно;
Figure 00000132
и
Figure 00000133
- доплеровские сдвиги частот на выходе основного и смежного КА соответственно за счет его сближения (удаления) с (от) НРТС K.
Предполагают, что нестабильность генератора частот КА известна и компенсируется. Влияние других эффектов на изменение частоты, например, гравитационный и релятивистский эффекты в рамках рассматриваемой задачи, пренебрежимо мало и поэтому не учитывают.
Для расчета ортогональных составляющих ортогональных составляющих векторов скоростей основного
Figure 00000134
и смежного КА
Figure 00000135
с использованием шести ЗС разработан алгоритм, схема которого представлена на фиг. 2.
На этапе 1 производят ввод исходных данных, в качестве которых выступают: момент времени измерения t0; координаты НРТС xK, yK, zK; координаты основного x1, y1, z1 и смежного КА х2, у2, z2 в момент времени t0; координаты шести ЗС
Figure 00000136
значения частот сдвига рабочих частот
Figure 00000137
и
Figure 00000138
основного и смежного КА соответственно; значения номиналов частот излучаемых каждой n-й ЗС ƒn; значения номиналов частот
Figure 00000139
и
Figure 00000140
принятых радиосигналов n-х ЗС In после их ретрансляции основным и смежным КА соответственно.
На этапе 2 рассчитывают расстояния
Figure 00000141
и
Figure 00000142
от и n-х ЗС In и НРТС K до основного КА S1 по формулам:
Figure 00000143
Figure 00000144
На этапе 3 рассчитывают расстояния
Figure 00000145
и
Figure 00000146
от и n-х ЗС In и НРТС K до смежного КА S2 по формулам:
Figure 00000147
Figure 00000148
На этапе 4 составляют аналитические выражения для радиальных скоростей
Figure 00000149
и
Figure 00000150
основного КА S1 относительно каждой из n-й ЗС In и НРТС K соответственно.
Для значений радиальных скоростей
Figure 00000151
и
Figure 00000152
основного КА S1 возможно записать формулы:
Figure 00000153
Figure 00000154
Согласно теореме о скалярном произведении векторов [6] справедливы равенства:
Figure 00000155
Figure 00000156
Модуль вектора скорости основного КА
Figure 00000157
равен:
Figure 00000158
Уравнения (А.7) и (А.8) с учетом (А.3), (А.4), (А.9)…(А.11) преобразуют к виду:
Figure 00000159
Figure 00000160
На этапе 5 составляют аналитические выражения для радиальных скоростей
Figure 00000161
и
Figure 00000162
смежного КА S2 относительно каждой из n-й ЗС In и НРТС K соответственно.
Для значений радиальных скоростей
Figure 00000163
и
Figure 00000164
смежного КА S2 возможно записать формулы:
Figure 00000165
Figure 00000166
Согласно теореме о скалярном произведении векторов [6] справедливы равенства:
Figure 00000167
Figure 00000168
Модуль вектора скорости смежного КА
Figure 00000169
равен:
Figure 00000170
Уравнения (А.14) и (А.15) с учетом (А.5), (А.6), (А.16)…(А.18) преобразуют к виду:
Figure 00000171
Figure 00000172
На этапе 6 вычисляют значения разностей частот Δƒn, между принятыми радиосигналами от каждой из n-й ЗС In после их ретрансляции основным и смежным КА соответственно по формулам:
Figure 00000173
На этапе 7 вычисляют ортогональные составляющие векторов скоростей основного
Figure 00000174
и смежного КА
Figure 00000175
Для значений разностей частот Δƒn, между принятыми радиосигналами от каждой из n-й ЗС In после их ретрансляции основным и смежным КА с учетом (A.1), (А.2) и (А.21) возможно составить выражения:
Figure 00000176
Доплеровские сдвиги частот на входе
Figure 00000177
и на выходе
Figure 00000178
основного КА S1 за счет его сближения (удаления) с (от) n-й ЗС In и НРТС K возможно представить в виде:
Figure 00000179
Figure 00000180
Доплеровские сдвиги частот на входе
Figure 00000181
и на выходе
Figure 00000182
смежного КА S2 за счет его сближения (удаления) с (от) n-й ЗС In и НРТС K возможно представить в виде:
Figure 00000183
Figure 00000184
Выражения (А.22) с учетом уравнений (А.23)…(А.26) для частного случая, когда n=1…6, преобразуют в систему линейных уравнений:
Figure 00000185
где коэффициенты при переменных и свободные члены равны:
Figure 00000186
Figure 00000187
Figure 00000188
Figure 00000189
Figure 00000190
Figure 00000191
Figure 00000192
Систему из шести линейных уравнений с шестью неизвестными (А.27) решают одним из известных методов, например, методом Крамера. Результатом решения системы уравнений (А.27) выступают ортогональные составляющие векторов скоростей основного
Figure 00000193
и смежного КА
Figure 00000194
в момент времени t0.
На этапе 8 осуществляют вывод результатов, в качестве которых выступают ортогональные составляющие векторов скоростей основного
Figure 00000195
и смежного КА
Figure 00000196
в момент времени t0.
В общем случае, когда количество ЗС N>6, алгоритм определения ортогональных составляющих основного
Figure 00000197
и смежного КА
Figure 00000198
остается прежним, с той лишь разницей, что система уравнений (А.27) будет содержать более шести уравнений. Тогда такую систему уравнений решают, например, методом наименьших квадратов.
Приложение Б
Алгоритм определения координат основного и смежного КА с использованием шести ЗС
Алгоритм определения координат основного x1, у1, z1 и смежного КА x2, у2, z2 в момент времени t0 с использованием ЗС In, где n=1…N - номер ЗС, N≥6, размещенных на позициях с известными координатами
Figure 00000199
основан на том, что каждой из временных задержек Δtn соответствует разности длин
Figure 00000200
и
Figure 00000201
траекторий InS1K и InS2K далее обозначаемые как ΔRn.
В качестве примера в настоящем приложении представлен вариант с шестью ЗС (n=1…6) как минимально необходимого количества для однозначного одномоментного определения координат основного и смежного КА x1, у1, z1 и х2, у2, z2
Разности длин ΔRn траекторий InS1K и InS2K рассчитывают по формулам:
Figure 00000202
где с=3×108 м/с - скорость света в вакууме.
Для расчета координат основного x1, у1, z1 и смежного КА х2, у2, z2 с использованием шести ЗС разработан алгоритм, схема которого представлена на фиг. 4.
На этапе 1 производят ввод исходных данных, в качестве которых выступают: момент времени измерения t0; координаты НРТС xK, yK, zK; координаты шести ЗС
Figure 00000203
временные задержки между радиосигналами принятыми от основного и смежного КА для каждой из выбранных ЗС Δtn; пороги точности δ1 и δ2 расчета координат основного и смежного КА соответственно.
На этапе 2 рассчитывают разности длин ΔRn траекторий InS1K и InS2K по формулам (Б.1).
На этапе 3 выбирают, на основе элементов Кеплеровой орбиты основного и смежного КА, координаты опорных точек S'1 (x'1, у'1, z'1) и S'2 (х'2, у'2, z'2), как первые приближения к координатам основного и смежного КА.
На этапе 4 рассчитывают разности длин ΔR'n траекторий InS'1K и InS'2K, при условии равенства координат основного и смежного КА координатам опорных точек S'1 (x'1, у'1, z'1) и S'2 (х'2, y'2, z'2) по формулам:
Figure 00000204
где
Figure 00000205
- расстояния от опорной точки S'1 до n-й ЗС In,
Figure 00000206
- расстояние от опорной точки S'1 до НРТС K,
Figure 00000207
- расстояния от опорной точки S'2 до n-й ЗС In,
Figure 00000208
- расстояние от опорной точки S'2 до НРТС K.
Расстояния
Figure 00000209
и
Figure 00000210
в свою очередь рассчитываются по формулам:
Figure 00000211
Figure 00000212
Figure 00000213
Figure 00000214
На этапе 5 рассчитывают невязки kn как разницы между определенными на этапе 4 разностями длин ΔR'n траекторий InS'1K и InS'2K, и разностями длин ΔRn траекторий InS1K и InS2K, определенными на этапе 2 соответственно:
Figure 00000215
На этапе 6 определяют поправки к координатам основного КА Δх1, Δу1, Δz1 и смежного КА Δх2, Δу2, Δz2.
Для определения поправок к координатам основного КА Δх1, Δу1, Δz1 и смежного КА Δх2, Δу2, Δz2 предварительно формируют систему линейных уравнений при разложении в ряд Тейлора функций
Figure 00000216
с точностью до первых членов:
Figure 00000217
где частные производные рассчитываются согласно выражениям
Figure 00000218
Figure 00000219
Figure 00000220
Figure 00000221
Figure 00000222
Figure 00000223
Решая систему линейных уравнений (Б.3) одним из известных методов, например, методом Крамера, получают поправки к координатам основного Δх1, Δу1, Δz1 и смежного КА Δх2, Δу2, Δz2.
На этапе 7 рассчитывают координаты новых опорных точек
Figure 00000224
Figure 00000225
Figure 00000226
Этапы 4-7 в совокупности составляют первую итерацию. Далее итерации повторяют, используя каждый раз новые опорные точки, полученные на предыдущей итерации. Количество необходимых итераций зависит требуемой точности определения координат основного и смежного КА. С точностью определения координат КА напрямую связан величинами шагов итераций для основного и смежного КА d1 и d2.
На этапе 8 определяют шаги итерации d1 и d2 как расстояния между текущими и предыдущими опорными точками:
Figure 00000227
Figure 00000228
На этапе 8 сравнивают d1 и d2 с порогами δ1 и δ2, задаваемым на этапе 1.
Необходимое число итераций, как правило, составляет 2…4. В качестве координат основного и смежного КА х1, у1, z1 и х2, у2, z2 выбирают значения координат опорных точек на последней итерации, вывод которых осуществляют на этапе 10.
В общем случае, когда количество ЗС N>6, алгоритм определения координат основного и смежного КА остается прежним, с той лишь разницей, что система уравнений (Б.3) будет содержать более шести уравнений. Тогда такую систему уравнений решают, например, методом наименьших квадратов.

Claims (4)

1. Способ определения ортогональных составляющих векторов скорости двух космических аппаратов (КА) с использованием земных станций (ЗС), заключающийся в том, что размещают наземную радиотехническую станцию (НРТС) K на позиции с известными координатами xK, yK, zK, принимают в момент времени t0 с помощью НРТС K радиосигналы, переданные опорными реперными станциями (ОРС) и ретранслированные основным КА S1, вычисляют координаты x1, y1, z1 основного КА S1 в момент времени t0, на основе частотных сдвигов радиосигналов, известных координат НРСТ, ОРС и предварительно заданной частоте сдвига рабочей частоты основного КА
Figure 00000229
рассчитывают ортогональные составляющие вектора скорости основного КА
Figure 00000230
в момент времени t0, отличающийся тем, что дополнительно выбирают смежный КА S2 с неизвестными значениями координат х2, у2, z2, неизвестными ортогональными составляющими вектора скорости
Figure 00000231
в момент времени t0 и заданной частотой сдвига
Figure 00000232
рабочей частоты, а в качестве ОРС используют выбранные ЗС In, размещенные на земной поверхности на позициях с известными координатами xIn, yIn, zIn, где n=1…N - номер ЗС, N≥6, излучающие радиосигналы на известных частотах ƒn в направлении основного и смежного КА, вычисляют координаты основного x1, у1, z1 и смежного КА х2, у2, z2 в момент времени t0, для каждой n-й ЗС In измеряют в НРТС K за счет корреляционной обработки радиосигналов значения номиналов частот
Figure 00000233
и
Figure 00000234
принятых радиосигналов после их ретрансляции основным КА S1 и смежным КА S2 соответственно, на основе известных координат НРТС xK, yK, zK, координат ЗС xIn, yIn, zIn и вычисленных координат основного КА x1, y1, z1 в момент времени t0 рассчитывают расстояния
Figure 00000235
и
Figure 00000236
от n-х ЗС In и НРТС К до основного КА S1, на основе известных координат НРТС xK, yK, zK, координат ЗС xIn, yIn, zIn и вычисленных координат смежного КА х2, у2, z2 в момент времени t0 рассчитывают расстояния
Figure 00000237
и
Figure 00000238
от n-х ЗС In и НРТС K до смежного КА S2, используя значения номиналов частот
Figure 00000239
и
Figure 00000240
принятых радиосигналов после их ретрансляции основным КА S1 и смежным КА S2, вычисляют значения разностей частот Δƒn между принятыми радиосигналами от каждой из n-й ЗС In после их ретрансляции основным и смежным КА соответственно, причем ортогональные составляющие векторов скоростей основного
Figure 00000241
и смежного КА
Figure 00000242
вычисляют по известным координатам НРТС xK, yK, zK, координатам ЗС xIn, yIn, zIn и координатам основного х1, у1, z1 и смежного КА х2, y2, z2, заданным частотам сдвига рабочих частот основного
Figure 00000243
и смежного КА
Figure 00000244
известным значениям номиналов частот излучаемых каждой n-й ЗС ƒn, рассчитанным расстояниям
Figure 00000245
и
Figure 00000246
от n-х ЗС In и НРТС K до основного КА, рассчитанным расстояниям
Figure 00000247
и
Figure 00000248
от n-х ЗС In и НРТС K до смежного КА S2.
2. Способ по п. 1, отличающийся тем, что смежный КА S2 выбирают такой, что бы его диапазон частот на линии "вверх"
Figure 00000249
имел одинаковые участки с диапазоном частот на линии "вверх" основного КА
Figure 00000250
а зона покрытия смежного КА Ω2 пересекалась с зоной покрытия основного КА Ω1.
3. Способ по п. 1, отличающийся тем, что выбирают такие ЗС In, значения номиналов излучаемых частот ƒn которых входят в диапазоны частот на линии "вверх" основного КА
Figure 00000251
и смежного КА
Figure 00000252
а каждая ЗС In находится в зонах покрытия как основного КА Ω1, так и смежного КА Ω2, при этом взаимные расстояния
Figure 00000253
между n-й и m-й ЗС, где m=1…N, m≠n максимальны.
4. Способ определения координат двух космических аппаратов (КА) с использованием земных станций (ЗС), заключающийся в том, что размещают наземную радиотехническую станцию (НРТС) K на позиции с известными координатами xK, yK, zK, выбирают начальные значения параметров орбиты основного КА S1, принимают в момент времени t0 с помощью НРТС K радиосигналы, переданные опорными реперными станциями (ОРС) и ретранслированные основным КА S1, на основе временных задержек радиосигналов системы вычисляют координаты основного КА х1, y1, z1 в момент времени t0, отличающийся тем, что дополнительно выбирают смежный КА S2 с неизвестными координатами х2, y2, z2 в момент времени t0, а в качестве ОРС используют выбранные земные станции (ЗС) In, размещенные на земной поверхности на позициях с известными координатами xIn, yIn, zIn, где n=X…N - номер ЗС, N≥6, для каждой n-й ЗС In измеряют в НРТС за счет корреляционной обработки радиосигналов значения временных задержек Δtn между принятыми радиосигналами после их ретрансляции основным и смежным КА соответственно, используя задержки Δtn, рассчитывают разности длин ΔRn траекторий InS1K и InS2K, вычисляют координаты основного x1, y1, z1 и смежного КА х2, у2, z2 на основе рассчитанных разностей длин ΔRn траекторий InS1K и InS2K, известных координат ЗС xIn, yIn, zIn и НРТС xK, yK, zK.
RU2020136690A 2020-11-06 2020-11-06 Способ определения ортогональных составляющих векторов скоростей и способ определения координат двух космических аппаратов с использованием земных станций RU2749878C1 (ru)

Priority Applications (1)

Application Number Priority Date Filing Date Title
RU2020136690A RU2749878C1 (ru) 2020-11-06 2020-11-06 Способ определения ортогональных составляющих векторов скоростей и способ определения координат двух космических аппаратов с использованием земных станций

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
RU2020136690A RU2749878C1 (ru) 2020-11-06 2020-11-06 Способ определения ортогональных составляющих векторов скоростей и способ определения координат двух космических аппаратов с использованием земных станций

Publications (1)

Publication Number Publication Date
RU2749878C1 true RU2749878C1 (ru) 2021-06-18

Family

ID=76377412

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2020136690A RU2749878C1 (ru) 2020-11-06 2020-11-06 Способ определения ортогональных составляющих векторов скоростей и способ определения координат двух космических аппаратов с использованием земных станций

Country Status (1)

Country Link
RU (1) RU2749878C1 (ru)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2787890C1 (ru) * 2022-04-05 2023-01-13 федеральное государственное казенное военное образовательное учреждение высшего образования "Военная орденов Жукова и Ленина Краснознаменная академия связи имени Маршала Советского Союза С.М. Буденного" Министерства обороны Российской Федерации Способ определения координат двух космических аппаратов с использованием земных станций и излучающей опорной реперной станции

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5745072A (en) * 1996-09-20 1998-04-28 The Johns Hopkins University Method and apparatus for precise noncoherent doppler tracking of a spacecraft
US5878034A (en) * 1996-05-29 1999-03-02 Lockheed Martin Corporation Spacecraft TDMA communications system with synchronization by spread spectrum overlay channel
WO2002079800A1 (en) * 2001-03-29 2002-10-10 Ses Astra S.A. Ranging system for determining ranging information of a spacecraft
EP1777159A1 (de) * 2005-10-20 2007-04-25 Astrium GmbH Anordnung und Verfahren zur Positions- und Lagebestimmung eines Flugkörpers, insbesondere eines Raumflugkörpers
RU2401437C2 (ru) * 2008-09-16 2010-10-10 Сорокин Виктор Леонидович Способ определения расстояний между космическим аппаратом и станциями
RU2508558C2 (ru) * 2008-07-24 2014-02-27 Сес Астра С.А. Система и способ оценки положения космического аппарата
RU2525343C1 (ru) * 2013-05-27 2014-08-10 Открытое акционерное общество "Российская корпорация ракетно-космического приборостроения и информационных систем" (ОАО "Российские космические системы") Способ одновременного определения шести параметров движения космического аппарата при проведении траекторных измерений и система для его реализации
RU2526401C1 (ru) * 2013-05-07 2014-08-20 Открытое акционерное общество "Российская корпорация ракетно-космического приборостроения и информационных систем" (ОАО "Российские космические системы") Способ радиотехнических доплеровских угломерных измерений космического аппарата и система для осуществления данного способа
RU2652603C1 (ru) * 2017-06-20 2018-04-27 федеральное государственное казенное военное образовательное учреждение высшего образования "Военная академия связи имени Маршала Советского Союза С.М. Буденного" Министерства обороны Российской Федерации Способ определения параметров орбиты искусственного спутника земли
RU2660676C1 (ru) * 2017-01-17 2018-07-09 Российская Федерация, от имени которой выступает Министерство обороны Российской Федерации Доплеровский измеритель скорости космического аппарата
RU2702098C1 (ru) * 2018-07-25 2019-10-04 федеральное государственное казенное военное образовательное учреждение высшего образования "Военная академия связи имени Маршала Советского Союза С.М. Буденного" Министерства обороны Российской Федерации Способ определения параметров орбиты искусственного спутника земли с использованием приемных опорных реперных станций
RU2708883C1 (ru) * 2018-10-01 2019-12-12 федеральное государственное казенное военное образовательное учреждение высшего образования "Военная академия связи имени Маршала Советского Союза С.М. Буденного" Министерства обороны Российской Федерации Способ определения параметров орбиты искусственного спутника земли с использованием приемо-передающих опорных реперных станций

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5878034A (en) * 1996-05-29 1999-03-02 Lockheed Martin Corporation Spacecraft TDMA communications system with synchronization by spread spectrum overlay channel
US5745072A (en) * 1996-09-20 1998-04-28 The Johns Hopkins University Method and apparatus for precise noncoherent doppler tracking of a spacecraft
WO2002079800A1 (en) * 2001-03-29 2002-10-10 Ses Astra S.A. Ranging system for determining ranging information of a spacecraft
EP1777159A1 (de) * 2005-10-20 2007-04-25 Astrium GmbH Anordnung und Verfahren zur Positions- und Lagebestimmung eines Flugkörpers, insbesondere eines Raumflugkörpers
RU2508558C2 (ru) * 2008-07-24 2014-02-27 Сес Астра С.А. Система и способ оценки положения космического аппарата
RU2401437C2 (ru) * 2008-09-16 2010-10-10 Сорокин Виктор Леонидович Способ определения расстояний между космическим аппаратом и станциями
RU2526401C1 (ru) * 2013-05-07 2014-08-20 Открытое акционерное общество "Российская корпорация ракетно-космического приборостроения и информационных систем" (ОАО "Российские космические системы") Способ радиотехнических доплеровских угломерных измерений космического аппарата и система для осуществления данного способа
RU2525343C1 (ru) * 2013-05-27 2014-08-10 Открытое акционерное общество "Российская корпорация ракетно-космического приборостроения и информационных систем" (ОАО "Российские космические системы") Способ одновременного определения шести параметров движения космического аппарата при проведении траекторных измерений и система для его реализации
RU2660676C1 (ru) * 2017-01-17 2018-07-09 Российская Федерация, от имени которой выступает Министерство обороны Российской Федерации Доплеровский измеритель скорости космического аппарата
RU2652603C1 (ru) * 2017-06-20 2018-04-27 федеральное государственное казенное военное образовательное учреждение высшего образования "Военная академия связи имени Маршала Советского Союза С.М. Буденного" Министерства обороны Российской Федерации Способ определения параметров орбиты искусственного спутника земли
RU2702098C1 (ru) * 2018-07-25 2019-10-04 федеральное государственное казенное военное образовательное учреждение высшего образования "Военная академия связи имени Маршала Советского Союза С.М. Буденного" Министерства обороны Российской Федерации Способ определения параметров орбиты искусственного спутника земли с использованием приемных опорных реперных станций
RU2708883C1 (ru) * 2018-10-01 2019-12-12 федеральное государственное казенное военное образовательное учреждение высшего образования "Военная академия связи имени Маршала Советского Союза С.М. Буденного" Министерства обороны Российской Федерации Способ определения параметров орбиты искусственного спутника земли с использованием приемо-передающих опорных реперных станций

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
S. P. Panko, M. S. Tsimbal Measuring the speed of the spacecraft // Research Naukograda. 2015, N 4 (14), pp. 25-29. *
V. V. Sevidov and A. O. Chemarov Determination of the coordinates of repeater satellites in the differential-range-finding geolocation system // News of higher educational institutions of Russia. Radio electronics. 2015, N 3, pp. 41-47. *
ПАНЬКО С.П., ЦИМБАЛ М.С. Измерение скорости космического аппарата // Исследования Наукограда. 2015 г., N 4(14), сс. 25-29. СЕВИДОВ В.В., ЧЕМАРОВ А.О. Определение координат спутников-ретрансляторов в разностно-дальномерной системе геолокации // Известия высших учебных заведений России. Радиоэлектроника. 2015 г., N 3, сс. 41-47. *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2787890C1 (ru) * 2022-04-05 2023-01-13 федеральное государственное казенное военное образовательное учреждение высшего образования "Военная орденов Жукова и Ленина Краснознаменная академия связи имени Маршала Советского Союза С.М. Буденного" Министерства обороны Российской Федерации Способ определения координат двух космических аппаратов с использованием земных станций и излучающей опорной реперной станции
RU2791153C1 (ru) * 2022-04-05 2023-03-03 федеральное государственное казенное военное образовательное учреждение высшего образования "Военная орденов Жукова и Ленина Краснознаменная академия связи имени Маршала Советского Союза С.М. Буденного" Министерства обороны Российской Федерации Способ определения ортогональных составляющих векторов скорости двух космических аппаратов с использованием земных станций и излучающей опорной реперной станции
RU2822687C1 (ru) * 2023-10-30 2024-07-11 федеральное государственное казенное военное образовательное учреждение высшего образования "Военная орденов Жукова и Ленина Краснознаменная академия связи имени Маршала Советского Союза С.М. Буденного" Министерства обороны Российской Федерации Способ определения координат двух космических аппаратов с использованием земных станций и излучающей опорной реперной станции

Similar Documents

Publication Publication Date Title
RU2708883C1 (ru) Способ определения параметров орбиты искусственного спутника земли с использованием приемо-передающих опорных реперных станций
RU2702098C1 (ru) Способ определения параметров орбиты искусственного спутника земли с использованием приемных опорных реперных станций
RU2593274C2 (ru) Способ и система для определения временных изменений при повторной передаче и распространении сигналов, используемых для измерения расстояний, синхронизирования исполнительных механизмов и выполнения геопривязки
RU2652603C1 (ru) Способ определения параметров орбиты искусственного спутника земли
US11821997B2 (en) Techniques for determining geolocations
CN106643744B (zh) 一种基于四程中继跟踪模式的远月面着陆器精密定位方法
RU2525343C1 (ru) Способ одновременного определения шести параметров движения космического аппарата при проведении траекторных измерений и система для его реализации
RU2749878C1 (ru) Способ определения ортогональных составляющих векторов скоростей и способ определения координат двух космических аппаратов с использованием земных станций
RU2750983C1 (ru) Способ определения ортогональных составляющих векторов скоростей и способ определения координат двух космических аппаратов с использованием земных станций
RU2750228C1 (ru) Способ определения ортогональных составляющих вектора скорости и способ определения координат космического аппарата с использованием земных станций
RU2750753C1 (ru) Способ определения ортогональных составляющих вектора скорости и способ определения координат космического аппарата с использованием земных станций
Sweeney et al. Enabling Mars radio occultation by smallsats
RU2801257C1 (ru) Способ определения ортогональных составляющих вектора скорости космического аппарата с использованием земной станции и излучающей опорной реперной станции
RU2788518C1 (ru) Способ определения ортогональных составляющих вектора скорости космического аппарата с использованием земных станций и излучающей опорной реперной станции
RU2803662C1 (ru) Способ определения координат космического аппарата с использованием земных станций и излучающей опорной реперной станции
RU2787890C1 (ru) Способ определения координат двух космических аппаратов с использованием земных станций и излучающей опорной реперной станции
RU2805667C1 (ru) Способ определения координат космического аппарата с использованием земных станций и излучающей опорной реперной станции
RU2791153C1 (ru) Способ определения ортогональных составляющих векторов скорости двух космических аппаратов с использованием земных станций и излучающей опорной реперной станции
Öz et al. TDOA based tracking measurement for geo satellites orbit determination: evaluation for the satellite operators
RU2822690C1 (ru) Способ определения ортогональных составляющих векторов скорости двух космических аппаратов с использованием земных станций и излучающей опорной реперной станции
RU2822687C1 (ru) Способ определения координат двух космических аппаратов с использованием земных станций и излучающей опорной реперной станции
Rowlands et al. Satellite laser altimetry on-orbit calibration techniques for precise geolocation
Cheung et al. Differencing Methods for 3D Positioning of Spacecraft
RU2679890C1 (ru) Способ определения местоположения абонентского терминала с помощью спутника-ретранслятора на низкой околоземной орбите
RU2684740C1 (ru) Способ определения местоположения абонентского терминала с помощью не менее двух спутников-ретрансляторов на низкой околоземной орбите