СПОСОБ ОПРЕДЕЛЕНИЯ ПОГРЕШНОСТЕЙ ПРИ ТРАЕКТОРНЫХ ИЗМЕРЕНИЯХ МЕЖПЛАНЕТНЫХ КОСМИЧЕСКИХ АППАРАТОВ ЗА СЧЕТ РАСПРОСТРАНЕНИЯ РАДИОСИГНАЛОВ В ИОНОСФЕРЕ ЗЕМЛИ И МЕЖПЛАНЕТНОЙ ПЛАЗМЕ Российский патент 2018 года по МПК B64G3/00 

Предлагаемый способ относится к области космонавтики, а именно к системам траекторных измерений межпланетных космических аппаратов.

Известны способы определения погрешностей при траекторных измерениях межпланетных космических аппаратов (патенты РФ №№2.107.306, 2.578.003 С1, 31663 U1; патент США №7.197.381; заявка США №2.010.012.4895; Яковлев О.И. Космическая радиофизика, М.: 1998, 432 с. и другие).

Из известных наиболее близким к предлагаемому является «Способ определения погрешностей при траекторных измерениях межпланетных космических аппаратов за счет распространения радиосигналов в ионосфере Земли и межпланетной плазме» (патент РФ №2.578.003, B64G 3/00, 2016), который и выбран в качестве прототипа.

Известно, что на практике при проведении траекторных измерений межпланетных космических аппаратов (МКА) используют измерение двух параметров движения МКА - наклонной дальности до МКА R и радиальной составляющей скорости МКА относительно измерительной станции R, которые последовательно измеряются несколькими территориально разнесенными измерительными станциями (ИС) [4].

В состав погрешности измерений R и входят погрешности, вызываемые распространением радиоволн в ионосфере Земли и ионизированной межпланетной плазме.

При совершенствовании космической техники аппаратурные погрешности измерений R и уменьшаются и по мере их уменьшения доля погрешностей за счет распространения радиоволн в суммарных погрешностях траекторных измерений увеличивается. Погрешности распространения становятся сравнимыми по величине с аппаратурными погрешностями и даже превышают их.

Учет погрешностей траекторных измерений за счет распространения радиоволн на практике осуществляется расчетным путем. Точность получения результатов при этом оказывается очень низкой. Поэтому актуальной задачей является реальное измерение погрешностей за счет распространения радиоволн объективными методами с необходимой точностью. Для этих целей используется двухчастотный метод в радиоинтерферометрии со сверхдлинными базами (РСДБ) [3].

В результате корреляционной обработки радиосигналов, принимаемых от МКА и опорного квазара, находящегося от него на минимальном угловом расстоянии с точно известными координатами, которые служит информационной основой для определения относительного углового положения МКА и квазара на небесной сфере. Точные положения используемых при измерениях квазаров даются соответствующими каталогами квазаров [5].

Разница в положении квазара, взятого из каталога, и измеренного радиоинтерферометром со сверхдлинными базами (РСДБ), вызвана погрешностями измерений, за счет распространения радиосигналов через межпланетную плазму и ионосферу Земли.

Технической задачей изобретения является повышение достоверности в измерении погрешностей траекторных измерений МКА, вызванных распространением измерительных радиосигналов через ионизированную среду, путем точного и однозначного определения местоположения ближайшего опорного квазара.

Поставленная задача решается тем, что способ определения погрешностей при траекторных измерениях межпланетных космических аппаратов, возникающих за счет прохождения радиосигнала через ионосферу Земли и межпланетную плазму, заключающийся, в соответствии с ближайшим аналогом, в том, что в процессе проведения траекторных измерений дальности и радиальной составляющей скорости МКА, кроме измерительных сигналов с МКА на частоте f_o1, принимают сигналы с квазара на излучаемой им частоте f₀₁ и частоте f₀₂, значительно отличающейся от частоты f₀₁, при этом обеспечивают минимальный сдвиг по времени между измерениями с МКА и квазара, выбирают угловое положение квазара таким образом, чтобы проекция его положения на небесной сфере была максимально приближена к положению МКА во время траекторных измерений так, чтобы трассы прохождения сигналов от МКА и квазара к измерительной станции практически совпадали по измеренным значениям принятых с квазара излучаемых им частот f_пр1 и f_пр2, определяют двухчастотным методом смещение частоты f_пр1 относительно частоты f₀₁, равной Δf₁, по измеренным значениям принятых с квазара частот f_пр1 и f_пр2 равное погрешности в измерениях скорости МКА из-за практического совпадения трасс прохождения сигналов от МКА и квазара, кроме того, по величине смещения частот Δf₁ определяют интегральную электронную концентрацию трассы квазар - измерительная станция, и по ее величине вычисляют временную задержку Δt прохождения сигнала f₀₁, которая практически равна погрешности измерения дальности Δt С, где С - скорость света, при траекторных измерениях МКА полученные данные передают в баллистический центр совместно с результатами траекторных измерений МКА для расчета траектории МКА, отличается от ближайшего аналога тем, что приемные пункты размещают в виде прямоугольного треугольника, формируют радиоинтерферометр с тремя сверхдлинными базами в азимутальной, угломестной и гипотенузной плоскостях соответственно, принимают сигналы от квазара на частоте f₀₂, измеряют разности фаз между ними:

где d₁, d₂, d₃ - сверхдлинные измерительные базы;

λ - длина волны с частотой f₀₂;

α, β, γ - азимут, угол места и угол ориентации квазара,

формируют фазовые шкалы отсчетов угловых координат α, β, γ: точные, но неоднозначные, перемножают сигнал, принимаемый вторым пунктом с сигналом, принимаемым первым пунктом и пропущенным через первый блок регулируемой задержки, выделяют низкочастотное напряжение, пропорциональное первой корреляционной функции R₁(τ), где τ-несущая временная задержка, изменяют временную задержку до получения максимального значения первой корреляционной функции R₁(τ), что соответствует равенству τ=τ_з1, и определяют азимут квазара

перемножают сигнал принимаемый третьим пунктом с сигналом, принимаемым первым пунктом и пропущенным через второй блок регулируемой задержки, выделяют низкочастотное напряжение, пропорциональное второй корреляционной функции R₂(τ), изменяют временную задержку τ до получения максимального значения второй корреляционной функции R₂(τ), что соответствует равенству τ=τ_З2, и определяют угол места квазара

перемножают сигнал, принимаемый вторым пунктом, с сигналом, принимаемым третьим пунктом и пропущенным через третий блок регулируемой задержки, выделяют низкочастотное напряжение, пропорциональное третьей корреляционной функции R₃(τ), изменяют временную задержку τ до получения максимального значения третьей корреляционной функции R₃(τ), что соответствует равенству τ=τ_з2, и определяют угол ориентации квазара

формируют временные шкалы отсчета угловых координат α, β и γ: грубые, но однозначные, по измеренным значениям угловых координат α, β и γ точно и однозначно определяют местоположение квазара.

Взаимное расположение приемных пунктов 3, 4, 5, квазара 1 и МКА 2 показано на фиг. 1. Взаимное расположение измерительной станции 9, ионизированного пространства 8, квазара 1 и МКА 2 показано на фиг. 2. Структурная схема устройства, реализующего предлагаемый способ, представлена на фиг. 3.

Устройство содержит три приемных пункта, каждый из которых состоит из последовательно включенных приемной антенны 11 (12, 13), полосового фильтра 14 (15, 16) и фазометра 17 (18, 19), выход которого подключен к компьютеру 7. К выходу второго полосового фильтра 15 последовательно подключено первый перемножитель 22.1, второй вход которого через первый блок 21.1 регулируемой задержки соединен с выходом первого полосового фильтра 14, первый фильтр 23.1 нижних частот и первый экстремальный регулятор 24.1, выход которого соединен с вторым входом первого блока 21.1 регулируемой задержки, второй вход которого через индикатор 25 азимута подключен к четвертому входу компьютера 7. К выходу третьего полосового фильтра 16 последовательно подключены второй перемножитель 22.2, второй вход которого через второй блок 21.2 регулируемой задержки соединен с выходом первого полосового фильтра 14, второй фильтр 23.2 нижних частот и второй экстремальный регулятор 24.2, выход которого соединен с вторым входом второго блока 21.2 регулируемой задержки, второй выход через индикатор 26 угла места подключен к пятому входу компьютера 7. К выходу второго полосового фильтра 15 последовательно подключены третий перемножитель 22.3, второй вход которого через третий блок 21.3 регулируемой задержки соединен с выходом третьего полосового фильтра 16, третий фильтр 23.3 нижних частот и третий экстремальный регулятор 24.3, выход которого подключен к второму входу третьего блока 21.3 регулируемой задержки, второй выход которого через индикатор 27 угла ориентации подключен к шестому входу компьютера 7. Первый блок 21.1 регулируемой задержки, первый перемножитель 22.1, первый фильтр 23.1 нижних частот и первый экстремальный регулятор 24.1 образуют первый коррелятор 20.1.

Второй блок 21.1 регулируемой задержки, второй перемножитель 22.2, второй фильтр 23.2 нижних частот и второй экстремальный регулятор 24.2. образуют второй коррелятор 20.2.

Третий блок 21.3 регулируемой задержки, третий перемножитель 22.3, третий фильтр 23.3 нижних частот и третий экстремальный регулятор 24.3 образуют третий коррелятор 20.3.

Предлагаемый способ реализуют следующим образом.

Приемные пункты 3,4 и 5 размещают в виде прямоугольного треугольника, образуя три сверхдлинные измерительные базы d₁, d₂ и d₃ в азимутальной, прямоугольной и гипотенузной плоскостях соответственно (фиг. 1). Принимаемые сигналы

u₁(t)=U₁cos(2πf₀₂t+ϕ₁),

u₂(t)=U₂cos(2πf_o2t+ϕ₂),

u₃(t)=U₃cos(2πf_o2t+ϕ₃),

от квазара 1 на частоте f₀₂ выделяются полосовыми фильтрами 14, 15 и 16 соответственно и подаются на входы фазометров 17, 18 и 19. Последние измеряют следующие разности фаз:

где d₁, d₂ и d₃ - сверхдлинные измерительные базы,

λ - длина волны с частотой f₀₂;

α, β и γ - азимут, угол места и угол ориентации квазара,

которые фиксируются компьютером 7. Так формируются фазовые шкалы отсчета угловых координат α, β и γ квазара: точные, но неоднозначные.

Принимаемые сигналы u₁(t), u₂(t) и u₃(t) одновременно поступают на два входа корреляторов 20.1, 20.2 и 20.3. Получаемые на выходе фильтров 23,1, 23.2 и 23.3 корреляционные функции R₁(τ), R₂(τ) и R₃(τ) имеют максимум при значении введенного регулируемого запаздывания:

τ=τ_з1, τ=τ_з2, τ=τ_з3.

Максимальные значения корреляционных функций R₁(τ), R₂(τ) и R₃(τ) поддерживаются автоматически с помощью экстремальных регуляторов 24.1, 24.3 и 24.3, воздействующих на управляющие входы блоков 21.1, 21.2 и 21.3 регулируемой задержки. Индикаторы 25, 26 м 27 градуируются непосредственно в значениях угловых координат α, β и γ квазара:

Значения угловых координат α, β и γ фиксируются компьютером 7. Так формирующиеся временные шкалы отсчета угловых координат α, β и γ квазара: грубые, но однозначные.

По измеренным значениям угловых координат α, β и γ определяется местоположение квазара. В предлагаемом способе измерения погрешностей траекторных измерений МКА используются космические источники излучения (квазары), расположенные на небесной сфере вблизи МКА за счет прохождения радиосигналов через ионизирующую среду путем измерения двух изменений частот, излучаемых квазаром, при приеме их наземной измерительной станцией (фиг. 2).

Соотношения частот определяются двумя выражениями:

где f₀₁ и f₀₂ - частоты сигналов, излучаемых квазаром.

При этом одна частота f₀₁ совпадает с частотой, на которой проводятся траекторные измерения МКА, а другая частота f₀₂ выбирается со значительным отличием от f₀₁. Для обеспечения когерентности частот принимаемых с квазара сигналов, гетеродины приемников измерительной станции 9 должны быть когерентны.

Δf₁ и Δf₂ - изменения значений излучаемых частот за счет прохождения сигналов через ионизированные области космического пространства;

f_пр1 и f_пр2 - значения частот сигналов, принимаемых наземными измерительными станциями;

- соотношение частот, излучаемых квазаром;

- (2) соотношение величины изменения частот, за счет

прохождения сигналов через ионизированное пространство.

Уравнение (2) может быть преобразовано

Решив уравнения (1) и (2) относительно Δf₁, получим величину изменения частоты f₀₁ через значения принимаемых частот f_пр1 и f_пр2

при значении m>2-3, единицей в знаменателе можно пренебречь. Тогда

Ввиду малого значения угла ψ (фиг. 2) линия прохождения сигналов квазара L₁ и МКА L₂ проходят через ионизированную плазму 8 на близком расстоянии и можно считать, что значение интегральной электронной концентрации ионосферы и межпланетной плазмы вдоль линий L₁ и L₂, близки.

Поэтому измеренное вдоль линии L₁ значение сдвига частоты сигнала f₀₁ будет иметь то же значение, что и при прохождении сигнала вдоль линии L₂, т.е. измеренное значение Δf₁ будет являться значением погрешности при измерении радиальной скорости МКА.

Для определения погрешности измерения дальности МКА необходимо решить обратную задачу, т.е. измеренному значению сдвига частоты Δf₁ определить интегральную электронную концентрацию ионосферы вдоль линии L₂, а по ней определить запаздывание радиосигнала Δt по линии L₂, которое фактически будет являться погрешностью измерения дальности до МКА вдоль линии L₂, Δt c=ΔR.

Такие измерения погрешностей за счет распространения радиосигналов в ионизированном пространстве с использованием квазара 1 необходимо проводить в каждом сеансе траекторных измерений МКА 2, а результаты передавать в баллистический центр для учета этих погрешностей при расчете траектории МКА.

Таким образом, предлагаемый способ по сравнению с прототипом и другими техническими решениями аналогичного назначения обеспечивает повышение достоверности в измерении погрешностей траекторных измерений МКА, вызванных распространением измерительных радиосигналов через ионизированную среду. Это достигается за счет точного и однозначного определения местоположения ближайшего квазара.

Следует отметить, что расположение приемных пунктов в виде прямоугольного треугольника продиктовано новой идеологией фазовой пеленгации источников радиоизлучений в пространстве, которая обеспечивает определение местоположения квазара пассивным методом. Причем для точного и однозначного определения местоположения квазара используются фазовые шкалы отсчета угловых координат α, β и γ: точные, но неоднозначные, и временные шкалы отсчета угловых координат α, β и γ: грубые, но однозначные, полученные за счет корреляционной обработки принимаемых сигналов.

Радиоинтерферометры со сверхдлинными базами, реализующие метод прямоугольного треугольника, отличаются новизной, оригинальностью и могут найти широкое практическое применение для точного и однозначного определения местоположения квазаров и МКА.

Литература

1. Патент РФ №2.578.003, В64 G 3/00, 2016.

2. Яковлев О.И. Космическая радиофизика. М.: «Научная книга», 1998.

3. Радиосистемы межпланетных космических аппаратов. / Под ред. А.С. Винницкого. - М.: Радио и связь, 1993.

4. Молотов Е.П. Наземные радиотехнические системы управления космическими аппаратами. - М.: ФИЗМАТЛИТ, 2004 - 256 с.

5. Каталог квазаров ICRF - Международная небесная система отсчета.

Изобретение относится к слежению за полётом межпланетных космических аппаратов (МКА) (2), куда вносит погрешности прохождение радиосигналов от МКА (на частоте f01) и близкого к нему на небесной сфере квазара (1) (на частотах f01 и f02) через ионизированную среду (8). По смещению Δf1 = f01- fпр1 принятой с квазара частоты fпр1 двухчастотным методом определяют погрешность измерений скорости МКА, а также интегральную электронную концентрацию среды (8) и по ее величине - временную задержку сигнала с частотой f01. Приемные пункты (9) размещают в вершинах прямоугольного треугольника-интерферометра с тремя сверхдлинными базами: в азимутальной, угломестной и гипотенузной плоскостях. Принятые интерферометром сигналы, излученные квазаром на частоте f02, обрабатывают, используя блоки регулируемой задержки и корреляционные функции. В итоге точно и однозначно определяют местоположение квазара (1). Техническим результатом изобретения является повышение достоверности определения погрешностей траекторных измерений МКА. 3 ил.

Способ определения погрешностей при траекторных измерениях межпланетных космических аппаратов (МКА) за счет распространения радиосигналов в ионосфере Земли и межпланетной плазме, заключающийся в том, что в процессе проведения траекторных измерений дальности и радиальной составляющей скорости МКА, кроме измерительных сигналов с МКА на частоте f₀₁, принимают сигналы с квазара, излучаемые им на частоте f₀₁ и частоте f₀₂, значительно отличающейся от частоты f₀₁, при этом обеспечивают минимальный сдвиг по времени между измерениями с МКА и квазара, выбирают угловое положение квазара таким образом, чтобы проекция его положения на небесной сфере была максимально приближена к положению МКА во время траекторных измерений, при этом трассы прохождения сигналов от МКА и квазара к измерительной станции практически совпадают, по измеренным значениям принятых с квазара частот f_пр1 и f_пр2 определяют двухчастотным методом смещение Δf₁ принятой с квазара частоты f_пр1 относительно частоты f₀₁, которое считают равным погрешности в измерениях скорости МКА из-за практического совпадения трасс прохождения сигналов от МКА и квазара, кроме того, по величине смещения частоты Δf₁ определяют интегральную электронную концентрацию вдоль трассы квазар - измерительная стация и по ее величине вычисляют временную задержку Δt прохождения сигнала с частотой f₀₁, которая практически равна погрешности измерения дальности Δt С, где С - скорость света, при траекторных измерениях МКА, полученные данные передают в баллистический центр совместно с результатами траекторных измерений МКА для расчета траектории МКА, отличающийся тем, что приемные пункты размещают в виде прямоугольного треугольника, формируют радиоинтерферометр с тремя сверхдлинными базами в азимутальной, угломестной и гипотенузной плоскостях, принимают сигналы от квазара на частоте f₀₂, измеряют разности фаз между ними

где d₁, d₂ и d₃ - сверхдлинные измерительные базы,

λ - длина волны с частотой f₀₂,

α, β и γ - азимут, угол места и угол ориентации квазара,

формируют фазовые шкалы отсчетов угловых координат α, β и γ: точные, но неоднозначные, перемножают сигнал, принимаемый вторым пунктом, с сигналом, принимаемым первым и пунктом, пропускаемым через первый блок регулируемой задержки, выделяют низкочастотное напряжение, пропорциональное первой корреляционной функции R₁(τ), где τ - текущая временная задержка, изменяют временную задержку до получения максимального значения первой корреляционной функции R₁(τ), что соответствует равенству τ=τ_з1, и определяют азимут квазара

перемножают сигнал, принимаемый третьим пунктом, с сигналом, принимаемым первым пунктом и пропущенным через второй блок регулируемой задержки, выделяют низкочастотное напряжение, пропорциональное второй корреляционной функции R₂(τ), изменяют временную задержку τ до получения максимального значения второй корреляционной функции R₂(τ), соответствующего равенству τ=τ_з2, и определяют угол места квазара

перемножают сигнал, принимаемый вторым пунктом, с сигналом, принимаемым третьим пунктом и пропущенным через третий блок регулируемой задержки, выделяют низкочастотное напряжение, пропорциональное третьей корреляционной функции R₃(τ), изменяют временную задержку τ до получения максимального значения третьей корреляционной функции R₃(τ), соответствующего равенству τ=τ_з3, и определяют угол ориентации квазара

формируют грубые, но однозначные временные шкалы отсчета угловых координат α, β и γ, по измеренным значениям угловых координат α, β и γ точно и однозначно определяют местоположение квазара.