RU2622511C1 - Method of clock synchronization and device for its implementation - Google Patents
Method of clock synchronization and device for its implementation Download PDFInfo
- Publication number
- RU2622511C1 RU2622511C1 RU2016104061A RU2016104061A RU2622511C1 RU 2622511 C1 RU2622511 C1 RU 2622511C1 RU 2016104061 A RU2016104061 A RU 2016104061A RU 2016104061 A RU2016104061 A RU 2016104061A RU 2622511 C1 RU2622511 C1 RU 2622511C1
- Authority
- RU
- Russia
- Prior art keywords
- output
- frequency
- input
- voltage
- signal
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G04—HOROLOGY
- G04R—RADIO-CONTROLLED TIME-PIECES
- G04R20/00—Setting the time according to the time information carried or implied by the radio signal
- G04R20/02—Setting the time according to the time information carried or implied by the radio signal the radio signal being sent by a satellite, e.g. GPS
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
Description
Предлагаемые способ и устройство относится к технике связи и радиолокации и может быть использованы в радиоинтерферометрии со сверхдлинными базами (РСДБ), а также в службе единого времени и частоты.The proposed method and device relates to communication and radar technology and can be used in radio interferometry with extra-long bases (VLBI), as well as in the service of a single time and frequency.
Известны способы и устройства синхронизации часов (авт. свид. СССР №№591.799, 614.416, 970.300, 1.180.835, 1.244.632, 1.278.800; патенты РФ №№2.003.157, 2.040.035, 2.177.167, 2.182.341, 2.248.669, 2.292.574, 2.301.437, 2.439.643, 2.535.653; патент США №5.519.759; патент Германии №3.278.943; патент ЕР №0.564.220; патент WO №99/97.826 и другие.Known methods and devices for clock synchronization (ed. Certificate of the USSR No. 591.799, 614.416, 970.300, 1.180.835, 1.244.632, 1.278.800; RF patents No. 2.003.157, 2.040.035, 2.177.167, 2.182. 341, 2.248.669, 2.292.574, 2.301.437, 2.439.643, 2.535.653; US patent No. 5.519.759; German patent No. 3.278.943; patent EP No. 0.564.220; patent WO No. 99/97.826 and others.
Из известных способов и устройств наиболее близкими к предлагаемым являются "Способ синхронизации часов и устройство для его реализации" (патент РФ №2.439.643, G04С 11/02, 2014), которые и выбраны в качестве базовых объектов.Of the known methods and devices closest to the proposed are the "Method of clock synchronization and device for its implementation" (RF patent No. 2.439.643, G04C 11/02, 2014), which are selected as the base objects.
Указанные способ и устройство обеспечивают сличение шкал времени, разнесенных на большие расстояния, и основаны на использовании дуплексного метода связи через геостационарный ИСЗ-ретранслятор.The aforementioned method and device provide a comparison of time scales spaced over long distances, and are based on the use of the duplex method of communication through a geostationary satellite repeater.
Однако потенциальные возможности известных технических решений используются не в полной мере.However, the potential capabilities of the known technical solutions are not fully utilized.
Используя радиопросвечивание атмосферы (тропосферы и части ионосферы) с помощью четырех частот можно оценить скорость распространения и направление прихода ионосферного возмущения. Ионосферное возмущение может быть вызвано солнечной, геомагнитной и сейсмической активностью, предвестником землетрясения, извержения вулканов, цунами, процессами грозовой активности, динамикой мощных штормовых циклонов, ядерными и иными крупными взрывами и пожарами, большими аварийными выбросами на атомных электростанциях, запусками космических аппаратов и ракет, излучениями мощных радиопередающих комплексов радиолокационного и связного назначения.Using radio illumination of the atmosphere (the troposphere and parts of the ionosphere) with the help of four frequencies, we can estimate the propagation velocity and direction of arrival of the ionospheric disturbance. The ionospheric disturbance can be caused by solar, geomagnetic and seismic activity, a harbinger of an earthquake, volcanic eruptions, tsunamis, thunderstorm processes, the dynamics of powerful storm cyclones, nuclear and other large explosions and fires, large accidental releases from nuclear power plants, spacecraft and rocket launches, radiations of powerful radio transmitting complexes of radar and connected purpose.
Ионосфера - это ионизированный атмосферный слой в диапазоне высот 5-500 км, который содержит свободные электроны. Наличие этих электронов вызывают задержку распространения ретранслированного сигнала, которая прямо пропорциональна концентрации электронов и обратно пропорционально квадрату частоты ретранслированного сигнала.The ionosphere is an ionized atmospheric layer in the altitude range of 5-500 km, which contains free electrons. The presence of these electrons causes a delay in the propagation of the relay signal, which is directly proportional to the electron concentration and inversely proportional to the square of the frequency of the relay signal.
Для вычисления ионосферной поправки используется измерение псевдодальности на двух частотах:To calculate the ionospheric correction, the measurement of pseudorange at two frequencies is used:
где ,Where ,
f1, f2 - частоты ретранслированных сигналов (f1=fc, f2=f3, fк1, fк2);f 1 , f 2 - frequencies of relayed signals (f 1 = f c , f 2 = f 3 , f k1 , f k2 );
Д1, Д2 - измеренные псевдодальности на частотах fi и f2 соответственно.D 1 , D 2 - measured pseudorange at frequencies fi and f 2, respectively.
Ионосферная поправка псевдодальности устраняет систематическую ошибку порядка 5 метров в определении вектора положения покоящегося наблюдателя.The ionospheric correction of the pseudorange eliminates a systematic error of the order of 5 meters in determining the position vector of the observer at rest.
Тропосфера - самый нижний слой атмосферы (до высоты 8-13 км). Она также обуславливает задержку распространения сигнала от ИСЗ-ретранслятора. Задержка сигнала в тропосфере также вызвана эффектами рефракции. В отличие от ионосферной задержки тропосферная задержка не зависит от частоты сигнала, она зависит от метеопараметров (давления, температуры, влажности), также от высоты ИСЗ-ретранслятора над горизонтом. Для вычисления тропосферной поправки измерения псевдодальности используют измерения температуры, давления воздуха и парциального давления водяного пара. Эти измерения доступны в сети Интернет для каждого наземного пункта А, В, С, Д и Е, фиг. 1.The troposphere is the lowest layer of the atmosphere (up to an altitude of 8-13 km). It also causes a delay in the propagation of the signal from the satellite repeater. Signal delay in the troposphere is also caused by refraction effects. Unlike the ionospheric delay, the tropospheric delay does not depend on the signal frequency, it depends on meteorological parameters (pressure, temperature, humidity), and also on the height of the satellite repeater above the horizon. To calculate the tropospheric correction, measurements of pseudorange use measurements of temperature, air pressure, and partial pressure of water vapor. These measurements are available on the Internet for each point A, B, C, D and E, FIG. one.
Соотношение для вычисления тропосферной поправки псевдодальности наземного наблюдателя имеет вид:The ratio for calculating the tropospheric correction of the pseudorange of the ground observer has the form:
где Т - температура в К;where T is the temperature in K;
Р - давление воздуха [мб];P is the air pressure [mb];
В - парциальное давление водяного пара [мб];In - partial pressure of water vapor [mb];
Θ - зенитный угол направления на ИСЗ-ретранслятор.Θ - zenith angle of direction to the satellite repeater.
Тропосферные задержки вызывают ошибки измерения псевдодальности в 1 м. Определение значения полного электронного содержания (ПЭС) ионосферы осуществляется по двухчастотным измерениям дальности между ИСЗ-ретранслятором наземным приемникомTropospheric delays cause measurement errors of 1 meter pseudorange. The determination of the total electronic content (TEC) of the ionosphere is carried out by two-frequency range measurements between the satellite repeater by the ground receiver
где f1, f2, λ1, λ2 - частоты и длины волн ретранслированных сигналов;where f 1, f 2 , λ 1 , λ 2 - frequency and wavelength of the relayed signals;
L1λ1. L2λ2 - фазовый путь трансионосферных сигналов (L1, L2 - число полных оборотов фазы);L 1 λ 1 . L 2 λ 2 is the phase path of transionospheric signals (L 1 , L 2 is the number of full rotations of the phase);
Θ - зенитный угол луча приемник - ИСЗ-ретранслятор.Θ - zenith angle of the beam; receiver - satellite repeater.
Совокупность лучей приемник - ИСЗ-ретранслятор в заданном регионе образует приемную решетку, каждый i-й элемент которой в момент времени t характеризуется изменением значения ПЭС Yi(t) и положением соответствующей ионосферной точки Xi(t), Yi(t) и Zi(t).The set of receiver-AES repeater beams in a given region forms a receiving array, each i-th element of which at time t is characterized by a change in the TEC value Y i (t) and the position of the corresponding ionospheric point X i (t), Y i (t) and Z i (t).
Временные ряды ПЭС отражают как регулярные изменения ПЭС в точке регистрации, так и вариации ПЭС, вызванные ионосферными возмущениями различного характера.The time series of TECs reflect both regular changes in TECs at the registration point and variations in TECs caused by ionospheric disturbances of a different nature.
Технической задачей изобретения является расширение функциональных возможностей известных технических решений путем определения скорости распространения и направления прихода ионосферного возмущения, вызванного различными дестабилизирующими факторами, за счет восстановления пространственного распределения полного электронного содержания ионосферы по данным радиопросвечивания атмосферы ретранслированными сигналами.An object of the invention is to expand the functionality of known technical solutions by determining the propagation velocity and direction of arrival of the ionospheric disturbance caused by various destabilizing factors, by restoring the spatial distribution of the total electronic content of the ionosphere according to atmospheric radio transmission signals relayed signals.
Поставленная задача решается тем, что способ синхронизации часов, основанный, в соответствии с ближайшим аналогом, на одновременном приеме разнесенными наземными пунктами шумоподобных СВЧ-сигналов с борта искусственного спутника Земли, когерентном преобразовании их к видеочастоте, цифровой регистрации принятых сигналов и определении временной задержки прихода одного и того же сигнала в пункты синхронизации методом корреляционной обработки зарегистрированных сигналов, по величине которой производят сличение шкал времени, при этом в начальный момент времени t1 по часам первого пункта с помощью кодовой последовательности формируют шумоподобный СВЧ-сигнал, регистрируют его на этом же пункте, сформированный сигнал преобразуют на частоту ω1 усиливают его по мощности, излучают усиленный сигнал в направлении на искусственный спутник Земли-ретранслятор, в тот же момент времени t1 по часам второго пункта с помощью такой же кодовой последовательности формируют такой же шумоподобный СВЧ-сигнал, регистрируют его на втором пункте, принимают бортовой аппаратурой искусственного спутника Земли-ретранслятора сигнала на частоте ω1, переизлучают его на первый и второй пункты на частоте ω2 с сохранением фазовых соотношений, в произвольный момент времени t3 по часам второго пункта аналогично формируют и регистрируют шумоподобный СВЧ-сигнал, сформированный сигнал преобразуют на частоту ω1, усиливают его по мощности, излучают усиленный сигнал в направлении того же искусственного спутника Земли-ретранслятора, в тот же момент времени t3 по часам первого пункта с помощью той же кодовой последовательности формируют такой же шумоподобный СВЧ-сигнал, регистрируют его на первом пункте, принимают шумоподобный сигнал по основному каналу на частоте ω2, преобразуют по частоте с использованием напряжения второго гетеродина, сдвинутого по фазе на +90°, выделяют первое и второе напряжения второй промежуточной частоты, второе напряжение второй промежуточной частоты сдвигают по фазе на -90°, суммируют с первым напряжением второй промежуточной частоты, перемножают полученное первое суммарное напряжение с принимаемым шумоподобным сигналом, выделяют гармоническое напряжение на частоте ωГ2 второго гетеродина, детектируют его и используют для разрешения дальнейшей обработки принимаемого сигнала, отличается от ближайшего аналога тем, что принимаемый шумоподобный сигнал по зеркальному каналу на частоте ωЗ1 преобразуют по частоте с использованием напряжения второго гетеродина и напряжения второго гетеродина, сдвинутого по фазе на +90°, выделяют третье и четвертое напряжения второй промежуточной частоты соответственно, сдвигают четвертое напряжение второй промежуточной частоты на +90°, суммируют с третьим напряжением второй промежуточной частоты, перемножают полученное второе суммарное напряжение с принимаемым шумоподобным сигналом, выделяют второе гармоническое напряжение на частоте ωГ2 второго гетеродина, детектируют его и используют для разрешения цифровой регистрации второго суммарного напряжения, определяют временную задержку прихода одного и того же сигнала на частоте ωЗ в первый и третий пункты синхронизации методом корреляционной обработки зарегистрированных сигналов, по величине которой производят сличение шкал времени, принимают шумоподобный сигнал по первому комбинационному каналу на частоте ωК1, преобразуют по частоте с использованием напряжения второго гетеродина и напряжение второго гетеродина, сдвинутого по фазе на +90°, выделяют пятое и шестое напряжения второй промежуточной частоты соответственно, сдвигают шестое напряжение второй промежуточной частоты на -90°, суммируют с пятым напряжением второй промежуточной частоты, перемножают полученное третье суммарное напряжение с принимаемым шумоподобным сигналом, выделяют третье гармоническое напряжение на частоте 2ωГ2 второго гетеродина, детектируют его и используют для разрешения цифровой регистрации третьего суммарного напряжения, определяют временную задержку прихода одного и того же сигнала на частоте ωК1 в первый и четвертый пункты синхронизации методом корреляционной обработки зарегистрированных сигналов, по величине которой производят сличение шкал времени, принимают шумоподобный сигнал по второму комбинационному каналу на частоте ωК2, преобразуют по частоте с использованием напряжения второго гетеродина и напряжения второго гетеродина, сдвинутого по фазе на +90°, выделяют седьмое и восьмое напряжение второй промежуточной частоты соответственно, сдвигают восьмое напряжение второй промежуточной частоты на +90°, суммируют с седьмым напряжением второй промежуточной частоты, перемножают полученное четвертое суммарное напряжение с принимаемым шумоподобным сигналом, выделяют четвертое гармоническое напряжение на частоте 2ωГ2 второго гетеродина, детектируют его и используют для разрешения цифровой регистрации четвертого суммарного напряжения, определяют временную задержку прихода одного и того же сигнала на частоте ωК2 в первый и пятый пункты синхронизации методом корреляционной обработки зарегистрированных сигналов, по величине которой производят сличение шкал времени, отличается от ближайшего аналога тем, что анализируют данные о полном электронном содержании в ионосфере Земли, которые получают в результате обработки сигналов, принимаемых двухчастотными приемниками наземных пунктов, с последующим формированием временных рядов полного электронного содержания и их фильтрацией в диапазоне периодов колебаний, соответствующих отклику ионосферы на воздействие источника ионосферного возмущения, при этом используют протяженную приемную антенную систему и последовательно проверяют гипотезу о значениях направления прихода и скорости распространения плоского фронта ионосферного возмущения путем формирования диаграммы направленности приемной системы и ее сканирования в заданном секторе, обзора пространства направление прихода - скорость распространения ионосферного возмущения за счет синтеза выходного сигнала приемной системы при синфазном суммировании рядов вариаций полного электронного содержания отдельных антенн системы с временными сдвигами, рассчитанными исходя из проверяемых значений направления ионосферного возмущения и расстояний, пройденных фронтом ионосферного возмущения между антенными системы в проверяемом направлении внутри сферического слоя ионосферы Земли, решение о правильности проверяемой гипотезы и обнаружении ионосферного возмущения принимается при превышении суммарным сигналом заданного порогового уровня, соответствующие значения направления прихода и фазовой скорости распространения ионосферного возмущения считаются оценочными значениями.The problem is solved in that the method of clock synchronization, based, in accordance with the closest analogue, on the simultaneous reception of spaced ground-based points of noise-like microwave signals from the artificial Earth satellite, coherently converting them to a video frequency, digital recording of received signals and determining the time delay of arrival of one of the same signal to synchronization points by the method of correlation processing of registered signals, the magnitude of which compares time scales, while in the initial time t 1 according to the clock of the first point using a code sequence form a noise-like microwave signal, register it at the same point, the generated signal is converted to a frequency ω 1 amplify it in power, emit an amplified signal in the direction of the artificial Earth repeater, at the same time point t 1, the clock of the second item with the same code sequence form the same noise-like microwave signal, it is recorded in the second paragraph, taking onboard equipment artificially of Earth relay signal of the satellite at a frequency ω 1, re-emit it at first and second points at the frequency ω 2 while preserving phase relationships, at an arbitrary time t 3 by the clock of the second paragraph similarly formed and recorded noise-like microwave signal generated signal is converted to frequency ω 1, increase its power emit the amplified signal in the direction of the same artificial satellite repeater earth, in the same moment of time t 3 by the clock of the first item by using the same code sequence that is formed oh same noise-like microwave signal is recorded it at the first location, receiving a noise-like signal in the main channel at frequency ω 2, is converted in frequency using the voltage of the second local oscillator shifted in phase by + 90 °, isolated first and second voltage of the second intermediate frequency, the second voltage of the second intermediate frequency is phase shifted by -90 °, summed with the first voltage of the second intermediate frequency, multiply the obtained first total voltage with the received noise-like signal, emit harmonic the voltage at a frequency ω G2 of the second local oscillator, it is detected and used to resolve further processing of the received signal, differs from the closest analogue in that the received noise-like signal through a mirror channel at a frequency of ω З1 is converted in frequency using the voltage of the second local oscillator and the voltage of the second local oscillator shifted in phase by + 90 °, the third and fourth voltages of the second intermediate frequency are isolated, respectively, the fourth voltage of the second intermediate frequency is shifted by + 90 °, summed with retim voltage of the second intermediate frequency, multiplies the resultant second sum voltage with the received spread spectrum signal, allocate the second harmonic voltage at frequency ω r2 of the second local oscillator is detected it and is used to permit digital recording of the second sum voltage, determine the time delay of arrival of the same at the frequency of signal ω З in the first and third synchronization points by the method of correlation processing of registered signals, the magnitude of which is used to compare the scales in belts, receive a noise-like signal through the first combination channel at a frequency ω K1 , convert the frequency using the voltage of the second local oscillator and the voltage of the second local oscillator, phase shifted by + 90 °, emit the fifth and sixth voltages of the second intermediate frequency, respectively, shift the sixth voltage of the second intermediate frequency by -90 °, summarized with the fifth voltage of the second intermediate frequency, multiplied the obtained third total voltage with the received noise-like signal, emit the third harmonic voltage at a frequency of 2ω G2 of the second local oscillator, detect it and use it to digitally record the third total voltage, determine the time delay of the arrival of the same signal at a frequency of ω K1 at the first and fourth synchronization points by the method of correlation processing of recorded signals, the magnitude of which is produced comparison of time scales, a noise-like signal is received over the second combination channel at a frequency ω K2 , converted in frequency using the voltage of the second local oscillator, and the voltage of the second local oscillator phase-shifted by + 90 °, the seventh and eighth voltage of the second intermediate frequency are allocated, respectively, the eighth voltage of the second intermediate frequency is shifted by + 90 °, summed with the seventh voltage of the second intermediate frequency, the resulting fourth total voltage is multiplied with the received noise-like signal , allocate the fourth harmonic voltage at frequency 2ω T2 of the second local oscillator is detected and it is used to permit registration of the fourth digital sum voltage I, determine the time delay of the arrival of the same signal at a frequency ω K2 at the first and fifth synchronization points by the method of correlation processing of recorded signals, the magnitude of which compares time scales, differs from the closest analogue in that they analyze data on the total electronic content in the ionosphere Earth, which is obtained as a result of processing signals received by dual-frequency receivers of ground points, followed by the formation of time series of the full electronic content and their ph by leaching in the range of oscillation periods corresponding to the response of the ionosphere to the action of the ionospheric disturbance source, an extended receiving antenna system is used and the hypothesis about the values of the arrival direction and propagation velocity of the plane front of the ionospheric disturbance is tested by forming the radiation pattern of the receiving system and scanning it in a given sector, space survey direction of arrival - ionospheric disturbance propagation velocity due to synthesis output signal of the receiving system during in-phase summation of the series of variations in the total electronic content of individual antennas of the system with time shifts calculated based on the verified values of the direction of the ionospheric disturbance and the distances traveled by the front of the ionospheric disturbance between the antenna of the system in the checked direction inside the spherical layer of the Earth’s ionosphere, the decision on the correctness of the checked hypotheses and the detection of ionospheric disturbances is accepted when the total signal exceeds a predetermined threshold match, the corresponding values of direction of arrival and phase velocity of ionospheric disturbance estimated values are considered.
Поставленная задача решается тем, что устройство синхронизации часов, содержащее, в соответствии с ближайшим аналогом, геостационарный ИСЗ-ретранслятор, первый и второй наземные пункты, каждый из которых содержит последовательно включенные стандарт частоты и времени, первый гетеродин, первый смеситель, второй вход которого через переключатель соединен с первым выходом генератора псевдошумового сигнала, усилитель первой промежуточной частоты, первый усилитель мощности, дуплексер, вход-выход которого связан с приемопередающей антенной, второй усилитель мощности, второй смеситель, второй вход которого через второй гетеродин соединен с первым выходом стандарта частоты и времени, первый усилитель второй промежуточной частоты, первый сумматор, первый перемножитель, второй вход которого соединен с выходом второго усилителя мощности, первый узкополосный фильтр, первый амплитудный детектор, первый ключ, второй вход которого соединен с выходом первого сумматора, второй клиппер, второй вход которого соединен с третьим выходом стандарта частоты и времени, второе буферное запоминающее устройство и первый измеритель задержек и их производных, выход которого является первым выходом наземного пункта, последовательно подключенные к второму выходу генератора псевдослучайного сигнала, первый клиппер, второй вход которого соединен с вторым выходом стандарта частоты и времени, и первое буферное запоминающее устройство, выход которого соединен с вторым входом первого измерителя задержек и их производных, последовательно подключенные к выходу второго гетеродина первый фазовращатель на +90°, третий смеситель, второй вход которого соединен с выходом второго усилителя мощности, второй усилитель второй промежуточной частоты и второй фазовращатель на -90°, выход которого соединен с вторым входом первого сумматора, последовательно подключенные к выходу второго усилителя второй промежуточной частоты третий фазовращатель на +90°, второй сумматор, второй вход которого соединен с выходом первого усилителя второй промежуточной частоты, второй перемножитель, второй вход которого соединен с выходом второго усилителя мощности, второй узкополосный фильтр, второй амплитудный детектор, второй ключ, второй вход которого соединен с выходом второго сумматора, третий клиппер, второй вход которого соединен с третьим выходом стандарта частоты и времени, третье буферное запоминающее устройство и третий измеритель задержек и их производных, второй вход которого соединен с выходом первого буферного запоминающего устройства, а выход является вторым выходом наземного пункта, к выходу первого перемножителя последовательно подключены третий узкополосный фильтр, третий амплитудный детектор, третий ключ, второй вход которого соединен с выходом первого сумматора, четвертый клиппер, второй вход которого соединен с третьим выходом стандарта частоты и времени, четвертое буферное запоминающее устройство и третий измеритель задержек и их производных, второй вход которого соединен с выходом первого буферного запоминающего устройства, а выход является третьим выходом наземного пункта, к выходу второго перемножителя последовательно подключены четвертый узкополосный фильтр, четвертый амплитудный детектор, четвертый ключ, второй вход которого соединен с выходом второго сумматора, пятый клиппер, второй вход которого соединен с третьим выходом стандарта частоты и времени, пятое буферное запоминающее устройство и четвертый измеритель задержек и их производных, второй вход которого соединен с выходом первого буферного запоминающего устройства, а выход является четвертым выходом наземного пункта, отличается от ближайшего аналога тем, что оно снабжено тремя блоками обработки сигналов, тремя блоками формирования временных рядов полного электронного содержания, третьим сумматором, пороговым блоком и компьютером, причем к выходу первого ключа последовательно подключены первый блок обработки сигналов, второй вход которого соединен с выходом четвертого ключа, первый блок формирования временных рядов полного электронного содержания, третий сумматор, пороговый блок и коммутатор, второй вход которого соединен с блоком управления антенной, к выходу первого ключа последовательно подключены второй блок обработки сигналов, второй вход которого соединен с выходом второго ключа, и второй блок формирования временных рядов полного электронного содержания, выход которого соединен с вторым входом третьего сумматора к выходу первого ключа последовательно подключены третий блок обработки сигналов, второй вход которого соединен с выходом третьего ключа, и третий блок формирования временных рядов полного электронного содержания, выход которого соединен с третьим входом третьего сумматора.The problem is solved in that the clock synchronization device, containing, in accordance with the closest analogue, a geostationary satellite repeater, first and second ground stations, each of which contains a frequency and time standard, a first local oscillator, a first mixer, and a second input through the switch is connected to the first output of the pseudo-noise signal generator, an amplifier of the first intermediate frequency, a first power amplifier, a duplexer, the input-output of which is connected to the transceiver antenna , a second power amplifier, a second mixer, the second input of which is connected through the second local oscillator to the first output of the frequency and time standard, the first amplifier of the second intermediate frequency, the first adder, the first multiplier, the second input of which is connected to the output of the second power amplifier, the first narrow-band filter, the first an amplitude detector, the first key, the second input of which is connected to the output of the first adder, the second clipper, the second input of which is connected to the third output of the frequency and time standard, the second buffer a flicker device and a first delay meter and their derivatives, the output of which is the first output of a ground station, connected in series to the second output of the pseudo-random signal generator, a first clipper, the second input of which is connected to the second output of the frequency and time standard, and a first buffer storage device, the output of which connected to the second input of the first delay meter and their derivatives, serially connected to the output of the second local oscillator, the first phase shifter + 90 °, the third mixer, the second the first input of which is connected to the output of the second power amplifier, the second amplifier of the second intermediate frequency and the second phase shifter -90 °, the output of which is connected to the second input of the first adder, sequentially connected to the output of the second amplifier of the second intermediate frequency third phase shifter + 90 °, the second adder the second input of which is connected to the output of the first amplifier of the second intermediate frequency, the second multiplier, the second input of which is connected to the output of the second power amplifier, the second narrow-band filter, a swarm amplitude detector, a second key, the second input of which is connected to the output of the second adder, a third clipper, the second input of which is connected to the third output of the frequency and time standard, a third buffer memory and a third delay meter and their derivatives, the second input of which is connected to the output of the first buffer storage device, and the output is the second output of a ground station, the third narrow-band filter, the third amplitude detector, the third key are connected in series to the output of the first multiplier the second input of which is connected to the output of the first adder, the fourth clipper, the second input of which is connected to the third output of the frequency and time standard, the fourth buffer memory and the third delay meter and their derivatives, the second input of which is connected to the output of the first buffer memory, and the output is the third output of a ground station, the fourth narrow-band filter, the fourth amplitude detector, the fourth key, the second input of which are connected to the output of the second multiplier connected to the output of the second adder, the fifth clipper, the second input of which is connected to the third output of the frequency and time standard, the fifth buffer memory and the fourth delay meter and their derivatives, the second input of which is connected to the output of the first buffer memory, and the output is the fourth output of the ground paragraph, differs from the closest analogue in that it is equipped with three signal processing units, three units for generating time series of full electronic content, the third adder ohm, a threshold unit and a computer, and the first signal processing unit, the second input of which is connected to the output of the fourth key, the first unit for generating time series of full electronic content, the third adder, the threshold unit and the switch, the second input of which is connected to the antenna control unit, the second signal processing unit, the second input of which is connected to the output of the second key, and the second block for generating time series are sequentially connected to the output of the first key of the complete electronic content, the output of which is connected to the second input of the third adder to the output of the first key, a third signal processing unit is connected in series, the second input of which is connected to the output of the third key, and the third block of time series formation of the complete electronic content, the output of which is connected to the third input of the third adder.
Геометрическая схема расположения наземных пунктов А, В, С, D, Е и ИСЗ-ретранслятора S изображена на фиг. 1, где введены следующие обозначения: О - центр масс Земли; d1, d2, d3, d4 - базы интерферометра, r - радиус-вектор ИСЗ-ретранслятора. Временная диаграмма дуплексного метода сличения часов первого А и второго В пунктов представлена на фиг. 2, где введены следующие обозначения: S, А, В - шкалы времени ИСЗ-ретранслятора и пунктов А и В соответственно. Структурная схема аппаратуры одного из пунктов (А), реализующей предлагаемый способ синхронизации часов, представлена на фиг. 3. Частотная диаграмма, иллюстрирующая преобразование сигналов по частоте, показана на фиг. 4. Геометрия определения координат удаленного точечного источника ионосферного возмущения изображена на фиг. 5.The geometric arrangement of ground posts A, B, C, D, E and the satellite repeater S is shown in FIG. 1, where the following notation is introduced: O is the center of mass of the Earth; d 1 , d 2 , d 3 , d 4 are the bases of the interferometer, r is the radius vector of the satellite repeater. The timing diagram of the duplex method of comparing the clocks of the first A and second B points is shown in FIG. 2, where the following notation is introduced: S, A, B are the time scales of the satellite repeater and points A and B, respectively. The structural diagram of the equipment of one of items (A) that implements the proposed method for clock synchronization is shown in FIG. 3. A frequency diagram illustrating frequency conversion of signals is shown in FIG. 4. The geometry of determining the coordinates of a remote point source of ionospheric disturbance is shown in FIG. 5.
Аппаратура наземного пункта А содержит последовательно включенные стандарт 1 частоты и времени, первый гетеродин 2.1, первый смеситель 5, второй вход которого через переключатель 4 соединен с первым выходом генератора 3 псевдослучайного сигнала, усилитель 6 первой промежуточной частоты, первый усилитель 7 мощности, дуплексер 8, вход-выход которого связан с приемопередающей антенной 9, второй усилитель 12 мощности, второй смеситель 13, второй вход которого через второй гетеродин 2.2 соединен с первым выходом стандарта частоты и времени, первый усилитель 14 второй промежуточной частоты, первый сумматор 22, первый перемножитель 23, второй вход которого соединен с выходом второго усилителя 12 мощности, первый узкополосный фильтр 24, первый амплитудный детектор 25, первый ключ 26, второй вход которого соединен с выходом первого сумматора 22, второй клиппер 15, второй вход которого соединен с третьим выходом стандарта 1 частоты и времени, второе буферное запоминающее устройство 16 и первый измеритель 17 задержек и их производных, выход которого является первым выходом 1 наземного пункта, к второму выходу генератора 3 псевдослучайного сигнала последовательно подключены первый клипер 10, второй вход которого соединен с вторым входом стандарта частоты и времени, и первое буферное запоминающее устройство 11, выход которого соединен с вторым входом первого измерителя 17 задержек и их производных. К выходу второго гетеродина 22 последовательно подключены первый фазовращатель 18 на +90°, третий смеситель 19, второй вход которого соединен с выходом второго усилителя 12 мощности, второй усилитель 20 второй промежуточной частоты и второй фазовращатель 21 на -90°, выход которого соединен с вторым входом первого сумматора 22. К выходу второго усилителя 20 второй промежуточной частоты последовательно подключены третий фазовращатель 27 на +90°, второй сумматор 28, второй вход которого соединен с выходом первого усилителя 14 второй промежуточной частоты, второй перемножитель 29, второй вход которого соединен с выходом второго усилителя 12 мощности, второй узкополосный фильтр 30, второй амплитудный детектор 31, второй ключ 32, второй вход которого соединен с выходом второго сумматора 28, третий клиппер 33, второй вход которого соединен с третьим выходом стандарта 1 частоты и времени, третье буферное запоминающее устройство 34 и второй измеритель 35 задержек и их производных, второй вход которого соединен с выходом первого буферного запоминающего устройства 11, а выход является вторым выходом II наземного пункта. К выходу первого перемножителя 23 последовательно подключены третий узкополосный фильтр 36, третий амплитудный детектор 37, третий ключ 38, второй вход которого соединен с выходом первого сумматора 22, четвертый клиппер 30, второй вход которого соединен с третьим выходом стандарта 1 частоты и времени, четвертое буферное запоминающее устройство 40 и третий измеритель 41 задержек и их производных, второй вход которого соединен с выходом первого буферного запоминающего устройства 11, а выход является третьим Ш выходом наземного пункта. К выходу второго перемножителя 29 последовательно подключены четвертый узкополосный фильтр 42, четвертый амплитудный детектор 43, четвертый ключ 44, второй вход которого соединен с выходом второго сумматора 28, пятый клиппер 45, второй вход которого соединен с третьим выходом стандарта 1 частоты и времени, пятое буферное запоминающее устройство 46 и четвертый измеритель 47 задержек и их производных, второй вход которого соединен с выходом первого буферного запоминающего устройства 11, а выход является четвертым IV выходом наземного пункта.The equipment of ground station A contains serially connected frequency and
К выходу первого ключа 26 последовательно подключены первый блок 48.1 обработки сигналов, второй вход которого соединен с выходом четвертого ключа 44, второй блок 49.1 формирования временных рядов полного электронного держания, третий сумматор 50, пороговый блок 51 и компьютер 52, второй вход которого соединен с блоком управления антенной 9. К выходу первого ключа 26 последовательно подключены второй блок 48.2 обработки сигналов, второй вход которого соединен с выходом второго ключа 32, и второй блок 49.2 формирования временных рядов полного электронного содержания, выход которого соединен с вторым входом третьего сумматора 50. К выходу первого ключа 26 последовательно подключены третий блок 48.3 обработки сигналов, второй вход которого соединен с выходом третьего ключа 38, и третий блок 49.3 формирования временных рядов полного электронного содержания, выход которого соединен с третьим входом третьего сумматора 50.The first signal processing unit 48.1 is connected in series to the output of the
Синхронизацию часов по предлагаемому способу осуществляют следующим образом.Clock synchronization by the proposed method is as follows.
В момент времени t1 A по часам первого пункта А с помощью кодовой последовательности формируют шумоподобный СВЧ-сигнал (сигнал α1) (фиг. 2):At time t 1 A, according to the clock of the first point A, a noise-like microwave signal is generated using a code sequence (signal α 1 ) (Fig. 2):
uс(t)=Uсcos[ωсt+ϕk(t)+ϕc], 0≤t≤Тc,u s (t) = U s cos [ω s t + ϕ k (t) + ϕ c ], 0≤t≤T c ,
где Uс, ωс, ϕс, Тс - амплитуда, несущая частота, начальная фаза и длительность сигнала;where U s , ω s , ϕ s , T s - amplitude, carrier frequency, initial phase and signal duration;
ϕk(t)={0,π} - манипулируемая составляющая фазы, отображающая закон фазовой манипуляции в соответствии с кодовой последовательностью M(t), причем ϕk(t)=const при kτЭ<t<(k+1)τЭ и может изменяться скачком при t=kτЭ, т.е. на границах между элементарными посылками (K=1, 2, …N-1);ϕk (t) = {0, π} is the manipulated component of the phase that displays the law of phase manipulation in accordance with the code sequence M (t), and ϕ k (t) = const for kτ Э <t <(k + 1) τ Э and can change abruptly at t = kτ Oe , i.e. at the borders between elementary premises (K = 1, 2, ... N-1);
tЭ, N - длительность и количество элементарных посылок, из которых составлен сигнал длительностью Tс(Тс=NτЭ), в генераторе 3 с помощью стандарта 1 частоты и времени.t E , N is the duration and number of chips that make up a signal of duration T s (T s = Nτ E ) in
Указанный сигнал поступает на вход клиппера 10, а затем регистрируется в буферном запоминающем устройстве 11. Регистрация синхронизируется стандартом 1 частоты и времени.The specified signal is input to the
Сформированный сигнал uc(t) поступает на первый вход первого смесителя 5, на второй вход которого подается напряжение первого гетеродина 2.1The generated signal u c (t) is fed to the first input of the
uГ1(t)=UГ1cos(ωГ1t+ϕГ1).u Г1 (t) = U Г1 cos (ω Г1 t + ϕ Г1 ).
На выходе смесителя 5 образуются напряжения комбинационных частот. Усилителем 6 выделяется напряжение первой промежуточной (суммарной) частотыAt the output of the
uпр1(t)=Uпр1cos[ϕпр1t+ϕk(t)+ϕпр1], 0≤t≤Тс,u CR1 (t) = U CR1 cos [ϕ CR1 t + ϕ k (t) + ϕ CR1 ], 0≤t≤T s ,
где ;Where ;
К1 - коэффициент передачи смесителя;To 1 - gear ratio of the mixer;
ωпр1=ωс+ωГ1 - первая промежуточная (суммарная) частота;ω CR1 = ω s + ω G1 - the first intermediate (total) frequency;
ϕпр1=ϕс+ϕГ1,ϕ pr1 = ϕ with + ϕ G1 ,
которое после усиления в усилителе 7 мощности через дуплексер 8 и приемопередающую антенну 9 излучается в направлении ИСЗ-ретранслятора на частоте ω1=ωпр1.which, after amplification in the
В тот же момент времени t1 A=t1 B по часам второго пункта В с помощью такой же кодовой последовательности M(t) формируют такой же шумоподобный СВЧ-сигнал (сигнал β1). Регистрируют его на втором пункте В (сигнал β1, который однако не отправляют на ретрансляцию). Принимают бортовой аппаратурой ИСЗ-ретранслятора на частоте ω1 (сигнал α1), переизлучают его на пункты А и В на частоте ω2 с сохранением фазовых соотношений на интервале tс.At the same time t 1 A = t 1 B according to the clock of the second point B using the same code sequence M (t) form the same noise-like microwave signal (signal β 1 ). Register it at the second point (signal β 1 , which, however, is not sent for relay). Accept on-board equipment of the satellite repeater at a frequency of ω 1 (signal α 1 ), re-emit it to points A and B at a frequency of ω 2 while maintaining phase relationships in the interval t s .
Ретранслированный сигнал (сигнал α2) на частоте ω2 Relay signal (signal α 2 ) at a frequency of ω 2
u2(t)=U2cos[ω2t+ϕk(t)+ϕ2], 0≤t≤Тс,u 2 (t) = U 2 cos [ω 2 t + ϕ k (t) + ϕ 2 ], 0≤t≤T s ,
принимается приемопередающей антенной 9 и через дуплексер 8 и усилитель 12 мощности поступает на первые входы второго 13 и третьего 19 смесителей и перемножителя 23. На вторые входы смесителей 13 и 19 подаются напряжения второго гетеродина 2.2:is received by the
uГ2(t)=UГ2(ωГ2t+ϕГ2),u Г2 (t) = U Г2 (ω Г2 t + ϕ Г2 ),
uГ3(t)=UГ2cos(ωГ2t+ϕГ2+90°).u Г3 (t) = U Г2 cos (ω Г2 t + ϕ Г2 + 90 °).
Причем частоты ωГ1 и ωГ2 первого 2.1 и второго 2.2 гетеродинов разнесены на вторую промежуточную частотуMoreover, the frequencies ω G1 and ω G2 of the first 2.1 and second 2.2 local oscillators are spaced at the second intermediate frequency
ωГ2-ωГ1=ωпр2.w G1 G2 -ω = ω WP2.
На выходах смесителей 13 и 19 образуются напряжения комбинационных частот. Усилителями 14 и 20 выделяются напряжения второй промежуточной (разностной) частоты:At the outputs of the
uпр2(t)=Uпр2cos[ωпр2(t)-ϕK1(t)+ϕпр2],u CR2 (t) = U CR2 cos [ω CR2 (t) -ϕ K1 (t) + ϕ CR2 ],
uпр3(t)=Uпр2cos[ωпр2(t)-ϕК1(t)+ϕпр2+90°], 0≤t≤Тс,u CR3 (t) = U CR2 cos [ω CR2 (t) -ϕ K1 (t) + ϕ CR2 + 90 °], 0≤t≤T s ,
где ;Where ;
ωпр2=ωГ2-ω2 - вторая промежуточная (разностная) частота;ω CR2 = ω G2 -ω 2 - the second intermediate (difference) frequency;
ϕпр2=ϕГ2-ϕ2.ϕ pr2 = ϕ G2 -ϕ 2 .
Напряжение uпр3(t) с входа усилителя 20 второй промежуточной частоты поступает на вход фазовращателя 21 на -90°, на выходе которого образуется напряжениеThe voltage u pr3 (t) from the input of the
uпр4(t)=Uпр2cos[ωпр2t-ϕk(t)+ϕпр2+90°-90°]=Uпр2cos[ωпр2t-ϕk(t)+ϕпр2], 0≤t≤Тс.u CR4 (t) = U CR2 cos [ω CR2 t-ϕ k (t) + ϕ CR2 + 90 ° -90 °] = U CR2 cos [ω CR2 t-ϕ k (t) + ϕ CR2 ], 0≤ t≤T s .
Напряжения uпр2(t) и uпр4(t) с выхода усилителя 14 и фазовращателя 21 на -90° поступают на два входа первого сумматора 22, на выходе которого образуется первое суммарное напряжениеVoltages u pr2 (t) and u pr 4 (t) from the output of
uΣ1(t)=UΣ1cos[ωпр2t-ϕk1(t)+ϕпр2], 0≤t≤Тс,u Σ1 (t) = U Σ1 cos [ω CR2 t-ϕ k1 (t) + ϕ CR2 ], 0≤t≤T s ,
где UΣ1=2Uпр2,where U Σ1 = 2U pr2 ,
которое поступает на второй вход перемножителя 23. На выходе последнего образуется гармоническое напряжениеwhich is fed to the second input of the
u1(t)=U1cos(ωГ2t+ϕГ2), 0≤t≤Tc,u 1 (t) = U 1 cos (ω Г2 t + ϕ Г2 ), 0≤t≤T c ,
где ;Where ;
K2 - коэффициент передачи перемножителя,K 2 - transfer coefficient of the multiplier,
которое выделяется узкополосным фильтром 24 (частота настройки ωн которого выбирается равной частоте второго гетеродина 2.2 ωн=ωГ2), детектируется амплитудным детектором 25 и поступает на управляющий вход ключа 26, открывая его. В исходном состоянии ключ 26 всегда закрыт.which is allocated by a narrow-band filter 24 (tuning frequency ω n which is chosen equal to the frequency of the second local oscillator 2.2 ω n = ω G2 ), is detected by the
Напряжение uΣ(t) с выхода сумматора 22 через открытый ключ 26 поступает на вход клиппера 15, где оно клиппируется и записывается в буферное запоминающее устройство 16. Регистрация синхронизируется стандартом 1 частоты и времени.The voltage u Σ (t) from the output of the
На втором шаге (при передаче сигнала из пункта В) переключатель 4 должен быть разомкнут и сигнал α3 из генератора 3 через клиппер 10 поступает в то же запоминающее устройство 11. Ретранслированный сигнал α4 записывается, как и α2, в запоминающее устройство 16.In the second step (when transmitting the signal from point B), the
В произвольный момент времени t3 В=t2 В+Θ по часам второго пункта аналогично формируют и регистрируют шумоподобный СВЧ-сигнал (сигнал (β3). Сформированный сигнал преобразуют на частоту ω1 усиливают его по мощности, излучают усиленный сигнал в направлении того же ИСЗ-ретранслятора.At an arbitrary point in time t 3 V = t 2 V + Θ, a noise-like microwave signal (signal (β 3 ) is generated and recorded in a similar manner by the clock of the second point. The generated signal is converted to a frequency ω 1, amplifies it by power, emits an amplified signal in the direction the same satellite repeater.
В тот же момент времени t3 B=t3 A по часам первого пункта А с помощью той же кодовой последовательности формируют такой же шумоподобный СВЧ-сигнал (сигнал α3). Регистрируют его на первом пункте А. Принимают бортовой аппаратурой ИСЗ-ретранслятора сигнал на частоте ω1 (сигнал α3), переизлучают его на пункты А и В на частоте ω2 с сохранением фазовых соотношений, принимают ретранслированный сигнал на обоих пунктах, преобразуют его на видеочастоту, регистрируют в моменты времени t4 A и t4 B соответственно (сигнал α4, β4).At the same time t 3 B = t 3 A , the same noise-like microwave signal (signal α 3 ) is formed using the same code sequence using the same code sequence. Register it at the first point A. Accept the on-board equipment of the satellite relay signal at a frequency of ω 1 (signal α 3 ), re-emit it to points A and B at a frequency of ω 2 while maintaining phase relationships, receive a relay signal at both points, convert it to video frequency, recorded at time t 4 A and t 4 B, respectively (signal α 4 , β 4 ).
Корреляционной обработкой двух пар зарегистрированных сигналов в измерителе 17 определяют на каждом пункте следующие временные задержки:The correlation processing of two pairs of registered signals in the
τt=β1⊗β2=t2 B-t1 B=a1+b2+(ΔВ И+ΔВ П+ΔS)+Δt,τ t = β 1 ⊗ β 2 = t 2 B -t 1 B = a 1 + b 2 + (Δ В И + Δ В П + ΔS) + Δt,
τ2=α3⊗α4=t4 A-t3 A=a3+b2+(ΔВ И+ΔА П+ΔS)-Δt,τ 2 = α 3 ⊗ α 4 = t 4 A -t 3 A = a 3 + b 2 + (Δ B AND + Δ A P + ΔS) -Δt,
τ3=α1⊗α2=t2 A-t2 A=a1+a2+(ΔФ И+ΔА П+ΔS),τ 3 = α 1 ⊗ α 2 = t 2 A -t 2 A = a 1 + a 2 + (Δ Ф И + Δ А П + ΔS),
τ4=β3⊗β4=t4 B-t3 B=b2+b3+(ΔВ И+ΔВ П+ΔS),τ 4 = β 3 ⊗ β 4 = t 4 B -t 3 B = b 2 + b 3 + (Δ B AND + Δ B P + ΔS),
и соответствующие им частоты интерференции Fi (i=1, 2, 3, 4), которые определяют производные этих задержек:and the corresponding interference frequencies F i (i = 1, 2, 3, 4), which determine the derivatives of these delays:
, ,
где ,Where ,
aj, bj (j=1, 2, 3) - время распространения сигнала между ИСЗ и пунктами А и В соответственно (фиг. 1);a j , b j (j = 1, 2, 3) is the propagation time of the signal between the satellite and points A and B, respectively (Fig. 1);
ΔА И, ΔВ И - задержки сигналов в излучающей аппаратуре обоих пунктов;Δ A AND , Δ B AND - signal delays in the radiating equipment of both points;
ΔA П, ΔВ П - задержки сигналов в приемно-регистрирующей аппаратуре;Δ A P , Δ V P - signal delay in the receiving and recording equipment;
ΔS - задержка сигналов в бортовом ИСЗ-ретрансляторе;ΔS - signal delay in the onboard satellite repeater;
Δt=t1 B-t1 A - искомая разность показаний часов в один и тот же физический момент.Δt = t 1 B -t 1 A is the desired difference in the clock readings at the same physical moment.
Полагая aj и bj линейными функциями с производными , , получаем:Assuming a j and b j linear functions with derivatives , we get:
, ,
гдеWhere
, ,
, ,
, ,
, ,
ΔА,В', ΔА,В'' _ задержки сигнала в атмосфере на частотах ω1 и ω2 соответственно; Δ A, B ', Δ A, B' '_ signal delay in the atmosphere at frequencies ω 1 and ω 2, respectively;
ν - релятивистская поправка (эффект Саньяка);ν - relativistic correction (Sagnac effect);
ω - угловая скорость вращения Земли;ω is the angular velocity of the Earth;
с - скорость света;c is the speed of light;
D - площадь четырехугольника OA'S'B', образуемого в экваториальной плоскости центром масс Земли, проекциями пунктов А, В и ИСЗ-ретранслятора S.D is the area of the quadrangle OA'S'B ', formed in the equatorial plane by the center of mass of the Earth, the projections of points A, B and the satellite S.
Поправки γ на подвижность ИСЗ-ретранслятора во время единичного измерения проще всего свести к нулю соответствующим выбором свободного параметра Θ:Corrections γ on the mobility of the satellite repeater during a single measurement is most easily reduced to zero by the corresponding choice of the free parameter Θ:
, ,
который следует в начале измерений рассчитывать по приближенным эфемеридным данным, а затем уточнить по результатам текущих измерений.which should be calculated at the beginning of measurements by approximate ephemeris data, and then clarified by the results of current measurements.
Что касается поправки δ на аппаратурные задержки, то ее можно найти путем калибровки по методу «нулевой базы».As for the correction δ for hardware delays, it can be found by calibration using the “zero base” method.
Атмосферная поправка ε также учитывается.The atmospheric correction ε is also taken into account.
На пункте В аппаратура работает аналогично, только порядок шагов там обратный. Для вычисления разности показаний часов Δt теперь достаточно обменяться между пунктами полученными цифровыми данными, что можно делать по обычным телефонным или телеграфным каналам связи.At point B, the equipment works similarly, only the order of steps there is the opposite. To calculate the difference between the clock readings Δt, it is now sufficient to exchange the received digital data between the points, which can be done via ordinary telephone or telegraph communication channels.
Описанная выше работа устройства, реализующего предлагаемый способ, соответствует приему полезных сигналов по основному каналу на частоте ω2 (фиг. 4).The above operation of the device that implements the proposed method, corresponds to the reception of useful signals on the main channel at a frequency of ω 2 (Fig. 4).
Если шумоподобный сигнал принимается по зеркальному каналу на частоте ωЗ If the noise-like signal is received by the image channel at frequency ω H
uЗ(t)=UЗcos[ωЗt+ϕK2(t)+ϕЗ], 0≤t≤TЗ,u З (t) = U З cos [ω З t + ϕ K2 (t) + ϕ З ], 0≤t≤T З ,
то усилителями 14 и 20 второй промежуточной частоты выделяются следующие напряжения:the
uПР5(t)=UПР5cos[ωПР2t+ϕK2(t)+ϕПР5],u PR5 (t) = U PR5 cos [ω PR2 t + ϕ K2 (t) + ϕ PR5 ],
uПР6(t)=UПР5cos[ωПР2t+ϕK2(t)+ϕПP5-90°], 0≤t≤ТЗ,u PR6 (t) = U PR5 cos [ω PR2 t + ϕ K2 (t) + ϕ ПP5 -90 °], 0≤t≤T З ,
где ;Where ;
ωПР2=ωЗ-ωГ2 - вторая промежуточная (разностная) частота;ω PR2 = ω З -ω G2 - the second intermediate (difference) frequency;
ϕПР5=ϕЗ _ϕГ2.ϕ PR5 = ϕ З _ ϕ Г2 .
Напряжение uПР6(t) с выхода усилителя 20 второй промежуточной частоты поступает на входы фазовращателей 21 на -90° и 27 на +90°, на выходе которых образуются следующие напряжения:The voltage u PR6 (t) from the output of the
uПР7(t)=UПР5cos[ωПР2t+ϕK2(t)+ϕПР5-90°-90°]=-UПР5cos[ωПР2t+ϕK2(t)+ϕПР5],u PR7 (t) = U PR5 cos [ω PR2 t + ϕ K2 (t) + ϕ PR5 -90 ° -90 °] = - U PR5 cos [ω PR2 t + ϕ K2 (t) + ϕ PR5 ],
uПР8(t)=UПР5cos[ωПР2+ϕК2(t)+ϕПР5-90°+90°]=UПР5cos[ωПР2t+ϕК2(t)+ϕK2(t)+ϕПР5], 0≤t≤TЗ.u PR8 (t) = U PR5 cos [ω PR2 + ϕ K2 (t) + ϕ PR5 -90 ° + 90 °] = U PR5 cos [ω PR2 t + ϕ K2 (t) + ϕ K2 (t) + ϕ PR5 ], 0≤t≤T Z.
Напряжения uПР5(t) и uПР7(t), поступающие на два входа сумматора 22, на его выходе компенсируются.The voltage u PR5 (t) and u PR7 (t) supplied to the two inputs of the
Напряжения uПР5(t) и uПР8(t) поступают на два входа сумматора 28, на выходе которого образуется второе суммарное напряжениеVoltages u PR5 (t) and u PR8 (t) are supplied to two inputs of the
uΣ2(t)=UΣcos[ωПР2t+ϕK2(t)+ϕПР5], 0≤t≤ТЗ,u Σ2 (t) = U Σ cos [ω PR2 t + ϕ K2 (t) + ϕ PR5 ], 0≤t≤T З ,
где UΣ2=2UПР5,where U Σ2 = 2U PR5 ,
которое поступает на второй вход второго перемножителя 29, на первый вход которого подается принимаемый шумоподобный сигнал uЗ(t) с выхода второго усилителя 12 мощности. На выходе перемножителя 29 образуется гармоническое напряжениеwhich is fed to the second input of the
u2(t)=U2cos(ωГ2t+ϕГ2), 0≤t≤TЗ,u 2 (t) = U 2 cos (ω Г2 t + ϕ Г2 ), 0≤t≤T З ,
где ,Where ,
которое выделяется вторым узкополосным фильтром 30, детектируется вторым амплитудным детектором 31 и поступает на управляющий вход второго ключа 32, открывая его. В исходном состоянии ключи 32, 38 и 44 всегда закрыты.which is allocated by the second narrow-
Следует отметить, что частота ωН1 настройки узкополосного фильтра 30 выбирается равной частоте второго гетеродина 22 ωН1=ωГ2, а частота ωН2 настройки узкополосных фильтров 36 и 42 выбирается равной второй гармонике частоты второго гетеродина 22 ωН2=2ωГ2.It should be noted that the tuning frequency ω Н1 of the narrow-
Второе суммарное напряжение uΣ2(t) с выхода второго сумматора 28 через открытый второй ключ 32 поступает на вход второго клипера 33, где оно клиппируется и записывается во второе буферное запоминающее устройство 34. Регистрация синхронизируется стандартом 1 частоты и времени. Затем определяют временную задержку прихода одного и того же шумоподобного сигнала на частоте ωЗ в первый А и третий С пункты синхронизации методом корреляционной обработки зарегистрированных сигналов, по величине которой производят сличение шкал времени.The second total voltage u Σ2 (t) from the output of the
Если шумоподобный сигналIf a noise-like signal
uK1(t)=UK1cos[ωK1t+ϕK3(t)+ϕK1], 0≤t≤ТK1,u K1 (t) = U K1 cos [ω K1 t + ϕ K3 (t) + ϕ K1 ], 0≤t≤T K1 ,
принимается по первому комбинационному каналу на частоте ωК1, то усилителями 14 и 20 второй промежуточной частоты выделяются следующие напряжения:is received on the first combinational channel at a frequency ω K1 , then the
uПР9(t)=UПР9cos[ωПР2t-ϕK3(t)+ϕПР9],u PR9 (t) = U PR9 cos [ω PR2 t-ϕ K3 (t) + ϕ PR9 ],
uПР10(t)=UПР9cos[ωПР2t-ϕK3(t)+ϕПР10+90°], 0≤t≤ТK1,u PR10 (t) = U PR9 cos [ω PR2 t-ϕ K3 (t) + ϕ PR10 + 90 °], 0≤t≤T K1 ,
где ;Where ;
ωПР2=2ωГ2-ωК1 - вторая промежуточная (разностная) частота;ω PR2 = 2ω Г2 -ω К1 - the second intermediate (difference) frequency;
ϕПР9=ϕГ2-ϕК1.ϕ PR9 = ϕ G2 -ϕ K1 .
Напряжение uПР10(t) с выхода усилителя 20 второй промежуточной частоты поступает на входы фазовращателей 21 на -90° и 27 на +90°, на выходе которых образуются следующие напряжения:The voltage u PR10 (t) from the output of the
uПР11(t)=UПР9cos[ωПР2t-ϕК3(e)+ϕПР9+90°-90°]=UПР9cos[ωПР2t-ϕК3(t)+ϕПР9],u PR11 (t) = U PR9 cos [ω PR2 t-ϕ K3 (e) + ϕ PR9 + 90 ° -90 °] = U PR9 cos [ω PR2 t-ϕ K3 (t) + ϕ PR9 ],
UПР12(t)=UПР9cos[ωПР2t-ϕK3(t)+ϕПР9-90°+90°]=-UПР9cos[ωПР2t-ϕK3(t)+ϕПР9], 0≤t≤TK1,U PR12 (t) = U PR9 cos [ω PR2 t-ϕ K3 (t) + ϕ PR9 -90 ° + 90 °] = - U PR9 cos [ω PR2 t-ϕ K3 (t) + ϕ PR9 ], 0 ≤t≤T K1 ,
Напряжения uПР9(t) и uПР12(t), поступающие на два входа сумматора 28, на его выходе компенсируются.The voltage u PR9 (t) and u PR12 (t) supplied to the two inputs of the
Напряжения uПР9(t) и uПР11(t) поступают на два входа сумматора 27, на выходе которого образуется третье суммарное напряжениеVoltages u PR9 (t) and u PR11 (t) are supplied to two inputs of the
uΣ3=UΣ3cos[ωПР2t-ϕK3(t)+ϕПР9], 0≤t≤ТК1,u Σ3 = U Σ3 cos [ω PR2 t-ϕ K3 (t) + ϕ PR9 ], 0≤t≤T K1 ,
где UΣ3=2UПР9,where U Σ3 = 2U PR9 ,
которое поступает на второй вход первого перемножителя 23, на второй вход которого подается принимаемый шумоподобный сигнал uK1(t) с выхода второго усилителя 12 мощности. На выходе перемножителя 23 образуется гармоническое напряжениеwhich is fed to the second input of the
uЗ(t)=UЗcos(2ωГ2t+ϕГ2), 0≤t≤TK1,u З (t) = U З cos (2ω Г2 t + ϕ Г2 ), 0≤t≤T K1 ,
где ;Where ;
которое выделяется третьим узкополосным фильтром 36, детектируется третьим амплитудным детектором 31 и поступает на управляющий вход третьего ключа 38, открывая его. При этом суммарное напряжение uΣ3(t) с выхода сумматора 22 через открытый ключ 38 поступает на вход клипера 39, где оно клиппируется и записывается в третье буферное запоминающее устройство 40. Регистрация синхронизируется стандартом 1 частоты и времени. Затем определяют временные задержки прихода одного и того же шумоподобного сигнала на частоте ωК1 в первый А и четвертый D пункты синхронизации методом корреляционной обработки зарегистрированных сигналов, по величине которой производят сличение шкал времени. Если шумоподобный сигналwhich is allocated by the third narrow-
uK2(t)=UK2cos[ωПР2t+ϕК4(t)+ϕК2], 0≤t≤ТК2,u K2 (t) = U K2 cos [ω PR2 t + ϕ K4 (t) + ϕ K2 ], 0≤t≤T K2 ,
принимается по второму комбинационному каналу на частоте ωК2 (фиг. 4), то усилителями 14 и 20 второй промежуточной частоты выделяются следующие напряжения:is taken along the second combinational channel at a frequency ω K2 (Fig. 4), then the following voltages are allocated by
uПР13(t)=UПР13cos[ωПР2t+ϕK4(t)+ϕПР13],u PR13 (t) = U PR13 cos [ω PR2 t + ϕ K4 (t) + ϕ PR13 ],
uПР14(t)=UПР13cos[ϕПР2t+ϕК4(t)+ϕПР13-90°], 0≤t≤ТК2,u PR14 (t) = U PR13 cos [ϕ PR2 t + ϕ K4 (t) + ϕ PR13 -90 °], 0≤t≤T K2 ,
где ;Where ;
ωПР2=ωK2-2ωГ2 - вторая промежуточная (разностная) частота;ω PR2 = ω K2 -2ω G2 - the second intermediate (difference) frequency;
ϕПР13=ϕK2-ϕГ2.ϕ PR13 = ϕ K2 -ϕ G2 .
Напряжение uПР14(t) с выхода усилителя 20 второй промежуточной частоты поступает на входы фазовращателей 21 на -90° и 27 на +90°, на выходе которых образуются следующие напряжения:The voltage u PR14 (t) from the output of the
UПР15(t)=UПР13cos[ωПР2t+ϕK4(t)+ϕПР13-90°-90°]=-UПР13cos[ωПР2t+ϕК4(t)+ϕПР13],U PR15 (t) = U PR13 cos [ω PR2 t + ϕ K4 (t) + ϕ PR13 -90 ° -90 °] = - U PR13 cos [ω PR2 t + ϕ K4 (t) + ϕ PR13 ],
uПР16(t)=uПР13cos[ωПР2t+ϕK4(t)+ϕПР13-90°+90°]=UПР13cos[ωПР2t+ϕK4(t)+ϕПР13].u PR16 (t) = u PR13 cos [ω PR2 t + ϕ K4 (t) + ϕ PR13 -90 ° + 90 °] = U PR13 cos [ω PR2 t + ϕ K4 (t) + ϕ PR13 ].
Напряжения uПР13(t) и uПР15(t), поступающие на два входа сумматора 22, на выходе компенсируются.The voltage u PR13 (t) and u PR15 (t) supplied to the two inputs of the
Напряжения uПР13(t) и uПР16(t) поступают на два входа сумматора 28, на выходе которого образуется четвертое суммарное напряжениеVoltages u PR13 (t) and u PR16 (t) are supplied to two inputs of the
uΣ4(t)=UΣ4cos[ωПР2t+ϕK4(t)+ϕПР13], 0≤t≤ТК2,u Σ4 (t) = U Σ4 cos [ω PR2 t + ϕ K4 (t) + ϕ PR13 ], 0≤t≤T K2 ,
где UΣ4=2UПР13,where U Σ4 = 2U PR13 ,
которое поступает на второй вход второго перемножителя 29, на первый вход которого подается принимаемый шумоподобный сигнал uK2(t) с выхода второго усилителя 12 мощности. На выходе перемножителя 29 образуется гармоническое напряжениеwhich is fed to the second input of the
u4(t)=U4cos(2ωГ2t+ϕГ2), 0≤t≤TK2,u 4 (t) = U 4 cos (2ω Г2 t + ϕ Г2 ), 0≤t≤T K2 ,
где ,Where ,
которое выделяется четвертым узкополосным фильтром 42, детектируется четвертым амплитудным детектором 43 и поступает на управляющий вход ключа 44, открывая его. При этом суммарное напряжение uΣ4(t) с выхода сумматора 28 через открытый ключ 44 поступает на вход клипера 45, где оно клиппируется и записывается в четвертое буферное запоминающее устройство 46. Регистрация синхронизируется стандартом 1 частоты и времени. Затем определяют временную задержку прихода одного и того же шумоподобного сигнала на частоте ωК2 в первый А и пятый Е пункты синхронизации методом корреляционной обработки зарегистрированных сигналов, по величине которой производят сличение шкал времени.which is allocated by the fourth narrow-
Для выделения характерных ионосферных возмущений ряда ПЭС подвергают процедуре специальной фильтрации в диапазоне периодов, соответствующих масштабу возмущения.To distinguish characteristic ionospheric disturbances, a number of TECs are subjected to special filtering in a range of periods corresponding to the scale of the disturbance.
Для этого напряжения с выхода ключей 26 и 44, 26 и 32, 26 и 38 поступают на входы блоков 48.1, 48.2 и 48.3 обработки сигналов, где анализируются данные о полном электронном содержании в ионосфере Земли. В блоках 49.1, 49.2 и 49.3 формируются временные ряды полного электронного содержания, соответствующие отклики ионосферы на воздействие источника ионосферного возмущения.For this, the voltage from the output of the
Для каждой пары проверяемых значений (α, v) формируется диаграмма направленности приемной системы и соответствующим образом ориентируется в фазовом пространстве [α, v], за счет синфазного суммирования в сумматоре 50 отдельных рядов ΔYc(t) приемной системыFor each pair of checked values (α, v), the radiation pattern of the receiving system is formed and accordingly oriented in the phase space [α, v], due to the common-mode summation in the adder of 50 separate rows ΔY c (t) of the receiving system
где p - количество антенн приемной системы,where p is the number of antennas of the receiving system,
временной сдвиг τi определяется как разность времени tj, j-го отсчета i-го суммарного ряда ПЭС и времени t0 регистрации ионосферного возмущения центральной антенной приемной системы τi=tj-t6, и выбирается исходя из минимизации выражения, описывающего динамику распространения возмущенияthe time shift τ i is defined as the difference between the time t j , the j-th count of the i-th total TEC series and the time t 0 of recording the ionospheric disturbance of the central antenna of the receiving system τ i = t j -t 6 , and is selected based on the minimization of the expression describing the dynamics outrage
где Ар, - расстояние, пройденное фронтом волны между i-м и центральным элементом приемной системы.where Ap, is the distance traveled by the wave front between the ith and the central element of the receiving system.
Для протяженных приемных систем расстояние Δрi рассчитывается с учетом кривизны Земли. С этой целью в заданном направлении α прихода волны ионосферного возмущения на высоте hmax задается удаленный точечный источник (обозначен точкой Е на фиг. 5), который будет являться полюсом ортодромической системы координат, экватор которой (сильная жирная линия на фиг. 5) проходит через центральный элемент приемной системы (точка А на фиг. 5). Тогда фронт волны, распространяющейся от удаленного точечного источника и проходящий через i-й элемент приемной системы (точка В на фиг. 5), будет представлять собой широтный круг (жирная прерывистая линия), параллельный экватору полученной ортодромической системы. Такая модель соответствует плоской волне ионосферного возмущения, распространяющегося на сфере Земли.For extended receiving systems, the distance Δp i is calculated taking into account the curvature of the Earth. To this end, in a given direction α of the arrival of the ionospheric disturbance wave at a height h max, a remote point source (indicated by point E in Fig. 5) is set, which will be the pole of the orthodromic coordinate system, the equator of which (the strong bold line in Fig. 5) passes through the central element of the receiving system (point A in Fig. 5). Then the front of the wave propagating from a remote point source and passing through the ith element of the receiving system (point B in Fig. 5) will be a latitudinal circle (a thick broken line) parallel to the equator of the resulting orthodromic system. Such a model corresponds to a plane wave of ionospheric disturbance propagating over the Earth's sphere.
Геометрические координаты (Хс, Yc, Zc) удаленного источника ионосферного возмущения определяется с использованием правил сферической тригонометрии. При этом рассматривается сферический треугольник, вершиной А которого является центральный элемент приемной системы с известными координатами (Хo, Yo, Zo). Вершиной С этого треугольника является северный полюс геоцентрической системы координат (О, О, R+hmax), где R - радиус Земли. Необходимо определить координаты третьей вершины Е, которая и будет являться удаленным источником. Чтобы удаленный источник Е являлся полюсом ортодромической системы координат, угловой размер стороны АЕ сферического треугольника задается равным π/2. В полученном сферическом треугольнике известны две стороны АС и АЕ, а также угол между ними <LCAE=α, что является типовой задачей сферического треугольника. С использованием теоремы косинусов сторон сферического треугольника определяется третья сторона и координаты (Хс, Yc, Zc) удаленного источника Е.The geometrical coordinates (X c , Y c , Z c ) of the remote source of the ionospheric disturbance are determined using the rules of spherical trigonometry. In this case, a spherical triangle is considered, the vertex A of which is the central element of the receiving system with known coordinates (X o , Y o , Z o ). The vertex C of this triangle is the north pole of the geocentric coordinate system (O, O, R + h max ), where R is the radius of the Earth. It is necessary to determine the coordinates of the third vertex E, which will be the remote source. So that the remote source E is the pole of the orthodromic coordinate system, the angular size of the side AE of the spherical triangle is set to π / 2. In the obtained spherical triangle, two sides AC and AE are known, as well as the angle between them <LCAE = α, which is a typical task of a spherical triangle. Using the cosine theorem of the sides of a spherical triangle, the third side and the coordinates (X c , Y c , Z c ) of the remote source E.
Решение о правильности проверяемой системы принимается при превышении суммарным сигналом заданного порогового уровня в пороговом блоке 51. При этом считается, что обнаружено ионосферное возмущение, а составляющим значения α и v, определяемые в компьютере 52, для которых суммарный сигнал приемной системы превысил пороговое значение, считаются оценками направления прихода и фазовой скорости распространения обнаруженного ионосферного возмущения.The decision on the correctness of the system under test is made when the total signal exceeds a predetermined threshold level in the
Способ синхронизации часов позволяет:The clock synchronization method allows you to:
- достичь предельной точности измерений (около ±0,1 нс) с помощью РСДБ техники и техники ретрансляции, которая уже широко используется на практике;- achieve extreme measurement accuracy (about ± 0.1 ns) using VLBI technology and relay technology, which is already widely used in practice;
- формировать необходимые для проведения измерения СВЧ-сигналы на наземных пунктах, что дает возможность постепенно наращивать точность измерений за счет оптимизации структуры сигнала и усовершенствования наземной техники регистрации без вмешательства в бортовую аппаратуру ИСЗ;- to form the microwave signals necessary for the measurement at ground points, which makes it possible to gradually increase the accuracy of measurements by optimizing the signal structure and improving the ground-based recording technique without interfering with the satellite onboard equipment;
- повысить оперативность измерений, т.е. довести интервал времени от начала измерений до получения результатов вплоть до нескольких десятков секунд (практически до времени корреляционной обработки сигналов);- increase the efficiency of measurements, i.e. bring the time interval from the beginning of measurements to obtain results up to several tens of seconds (almost to the time of correlation signal processing);
- избежать установки на борту ИСЗ высокостабильных хранителей времени и измерителей временных интервалов, ограничить бортовую аппаратуру только системой фазостабильной ретрансляции СВЧ-сигналов.- to avoid the installation on board of a satellite of highly stable time-keepers and time interval meters, to limit the on-board equipment to only a phase-stable microwave signal relay system.
Таким образом, предлагаемые способ и устройство по сравнению с прототипами и другими техническими решениями аналогичного назначения обеспечивают определение скорости распространения и направление прихода ионосферного возмущения, вызванного различными дестабилизирующими факторами, за счет восстановления пространственного распределения полного электронного содержания по данным радиопросвечивания атмосферы ретранслированными сигналами.Thus, the proposed method and device, in comparison with prototypes and other technical solutions of a similar purpose, provides the determination of the propagation velocity and direction of arrival of the ionospheric disturbance caused by various destabilizing factors by restoring the spatial distribution of the total electronic content according to atmospheric radio-transmission data by relayed signals.
Тем самым функциональные возможности известных технических решений расширены.Thus, the functionality of known technical solutions is expanded.
Claims (2)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2016104061A RU2622511C1 (en) | 2016-02-08 | 2016-02-08 | Method of clock synchronization and device for its implementation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2016104061A RU2622511C1 (en) | 2016-02-08 | 2016-02-08 | Method of clock synchronization and device for its implementation |
Publications (1)
Publication Number | Publication Date |
---|---|
RU2622511C1 true RU2622511C1 (en) | 2017-06-16 |
Family
ID=59068434
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2016104061A RU2622511C1 (en) | 2016-02-08 | 2016-02-08 | Method of clock synchronization and device for its implementation |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2622511C1 (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109613815A (en) * | 2018-12-24 | 2019-04-12 | 北京无线电计量测试研究所 | A kind of time interval measurement device based on time-stretching |
CN114461012A (en) * | 2022-01-19 | 2022-05-10 | 许昌许继软件技术有限公司 | Method and device for acquiring running time stamps of different clock domains of embedded system |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1840365A1 (en) * | 1985-04-01 | 2006-09-27 | Российский институт радионавигации и времени | System for synchronizing clocks |
RU2310221C1 (en) * | 2006-03-30 | 2007-11-10 | Институт Прикладной Астрономии Российской Академии Наук | Device for synchronizing clock |
RU2386159C2 (en) * | 2007-09-10 | 2010-04-10 | Российская Федерация, от имени которой выступает Министерство обороны Российской Федерации | Clock synchronisation system |
RU2439643C1 (en) * | 2010-07-07 | 2012-01-10 | Учреждение Российской академии наук Институт прикладной астрономии РАН | Method of clock synchronisation and device for its realisation |
-
2016
- 2016-02-08 RU RU2016104061A patent/RU2622511C1/en active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1840365A1 (en) * | 1985-04-01 | 2006-09-27 | Российский институт радионавигации и времени | System for synchronizing clocks |
RU2310221C1 (en) * | 2006-03-30 | 2007-11-10 | Институт Прикладной Астрономии Российской Академии Наук | Device for synchronizing clock |
RU2386159C2 (en) * | 2007-09-10 | 2010-04-10 | Российская Федерация, от имени которой выступает Министерство обороны Российской Федерации | Clock synchronisation system |
RU2439643C1 (en) * | 2010-07-07 | 2012-01-10 | Учреждение Российской академии наук Институт прикладной астрономии РАН | Method of clock synchronisation and device for its realisation |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109613815A (en) * | 2018-12-24 | 2019-04-12 | 北京无线电计量测试研究所 | A kind of time interval measurement device based on time-stretching |
CN114461012A (en) * | 2022-01-19 | 2022-05-10 | 许昌许继软件技术有限公司 | Method and device for acquiring running time stamps of different clock domains of embedded system |
CN114461012B (en) * | 2022-01-19 | 2024-05-10 | 许昌许继软件技术有限公司 | Method and device for acquiring running time stamps of different clock domains of embedded system |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US4463357A (en) | Method and apparatus for calibrating the ionosphere and application to surveillance of geophysical events | |
Ruiz-Rodon et al. | Nearly zero inclination geosynchronous SAR mission analysis with long integration time for earth observation | |
Alizadeh et al. | Ionospheric effects on microwave signals | |
Coleman | A ray tracing formulation and its application to some problems in over-the-horizon radar | |
Park et al. | A generic level 1 simulator for spaceborne GNSS-R missions and application to GEROS-ISS ocean reflectometry | |
Liou | Radio occultation method for remote sensing of the atmosphere and ionosphere | |
Gladyshev et al. | A hardware-software complex for modelling and research of near navigation based on pseudolites | |
Spits | Total Electron Content reconstruction using triple frequency GNSS signals | |
Galushko et al. | Frequency-and-angular HF sounding and ISR diagnostics of TIDs | |
Mitchell et al. | Ionospheric data assimilation applied to HF geolocation in the presence of traveling ionospheric disturbances | |
RU2560094C2 (en) | Method of determining propagation speed and direction of arrival of ionospheric perturbation | |
RU2622511C1 (en) | Method of clock synchronization and device for its implementation | |
Solomentsev et al. | Three-dimensional assimilation model of the ionosphere for the European region | |
Sweeney et al. | Enabling Mars radio occultation by smallsats | |
RU2535653C1 (en) | Clock synchronisation method and device therefor | |
Bernhardt et al. | Ionospheric applications of the scintillation and tomography receiver in space (CITRIS) mission when used with the DORIS radio beacon network | |
RU2656617C1 (en) | Method of sensing the plasma layer of the geomagnetic tail and ionosphere of the earth | |
Kunitsyn et al. | Earthquake prediction research using radio tomography of the ionosphere | |
Bernhardt et al. | Atmospheric studies with the Tri-Band Beacon instrument on the COSMIC constellation | |
Geeraert | Multi-satellite orbit determination using interferometric observables with RF localization applications | |
RU2624634C1 (en) | Method of determining speed of distribution and direction of ionospheric perturbation | |
Kikuchi et al. | Recent status of SELENE-2/VLBI instrument | |
Garrison et al. | SNOOPI: Demonstrating Earth remote sensing using P-band signals of opportunity (SoOp) on a CubeSat | |
RU2655164C2 (en) | System for determining propagation speed and direction of arrival of ionospheric disturbance | |
Lay et al. | New lightning‐derived vertical total electron content data provide unique global ionospheric measurements |