Прогнозирование помехоустойчивости спутниковой связи по результатам мониторинга индекса мерцаний ионосферы

Известно [1-2], что помехоустойчивость систем спутниковой связи (ССС) может снижаться на несколько порядков при возникновении мерцаний (сцинтилляций, замираний) принимаемых сигналов. Они вызваны возмущениями ионосферы (обычно в районах полярных и экваториальных широт), которые проявляются в образовании интенсивных мелкомасштабных неоднородностей электронной концентрации. Эти мелкомасштабные неоднородности обусловливают рассеяние радиоволн и их многолучевое распространение в точку приема. Мерцания принимаемых в ССС сигналов обычно подчиняются m -распределению Накагами и могут быть очень глубокими. По экспериментальным данным [2], в области экватора индекс мерцаний S 4 может возрастать с обычных (для средних широт) значений S₄ = 0..Д2 до величины S₄ =1.

Поэтому m -параметр Накагами, связанный с ^4 как m = 5) , может достигать на экваторе значения m = 1, соответствующего релеевским замираниям. Однако иногда достаточно глубокие мерцания принимаемых в ССС сигналов (S₄ >0,3) могут наблюдаться и в средних широтах. В связи с этим возникает необходимость организации мониторинга возмущений ионосферы, вызванных ростом флуктуаций электронной концентрации в ее мелкомасштабных неоднородностях и индекса мерцаний принимаемых сигналов, в целях прогнозирования помехоустойчивости ССС.

Постановка задачи

Согласно [3] помехоустойчивость приема сигналов при передаче дискретных сообщений и распределении Накагами принимаемых сигналов определяется зависимостью Р_ош =и(^²)’ m) вероятности ошибочного приема информационных символов от среднего энергетического отношения «сигнал/шум» на входе приемника (^²) и m -параметра Накагами.

Индекс мерцаний S 4 в трансионосферном (спутниковом) канале связи (определяющий параметр Накагами m = s;h проще всего можно определить путем непосредственных измерений с помощью двухчастотного приемника NovAtel GPStanion – 6-спутниковой радионавигационной системы (СРНС) GPS / ГЛОНАСС [4]. Однако временные ряды измерений индекса мерцаний на выходе двухчастотного приемника GPStanion – 6 формируются с недостаточной точностью.

Это обусловлено тем, что интервал дискретизации рядов ^ 4 (^ ) составляет 1 мин. [4], тогда как для измерения параметров мелкомасштабных (с размерами порядка 100…1000 м) неоднородностей ионосферы необходима частота отсчетов не ниже 50 Гц (то есть через 0,02 с) [5]. Кроме того, в ряды измерений во времени индекса мерцаний принимаемых сигналов S'4 (0, помимо многолучевости из-за рассеяния на неоднородностях ионосферы, вносит свой вклад и многолучевость из-за отражения на неровностях рельефа местности.

Поэтому для повышения точности определения индекса мерцаний S₄ (/) принимаемого сигнала целесообразно использовать возможности [6] двухчастотного приемника GPStanion-6 по проведению измерений временных рядов полного электронного содержания (ПЭС) ионосферы n? ^ с последующим выделением из них рядов вариаций ПЭС, обусловленных мелкомасштабными неоднородностями ^N_TM^ . Это позволит оценить влияние статистических характеристик (например, среднеквадратического отклонения) мелкомасштабных вариаций ПЭС ионосферы ^ATM (0 ) на изменение во времени индекса мерцаний принимаемого в ССС сигнала ^4(0 = ^(О"д™(0), параметра Накагами трансионосферного канала связи m^ = s;4n и на помехоустойчивость приема сигналов P^tVv^h-^m^Y

Целью статьи является разработка структуры построения комплекса прогнозирования помехоустойчиво сти системы спутниковой связи по результатам измерений статистических характеристик мелкомасштабных вариаций полного электронного содержания ионосферы на выходе двухчастотного приемника спутниковой радионавигационной системы и индекса мерцаний принимаемого сигнала.

Принцип построения и работы комплекса прогнозирования помехоустойчивости ССС

Схема на рис.1 иллюстрирует принцип построения и работы комплекса прогнозирования помехоустойчивости ССС по данным мониторинга ионосферы с помощью СРНС. В состав комплекса входит двухчастотный приемник СРНС типа GPStation – 6, блок определения среднеквадратического отклонения (СКО) мелкомасштабных вариаций ПЭС ионосферы °" AW и блок прогнозирования помехоустойчивости "^ОШ .

Суть работы двухчастотной СРНС состоит в следующем [5]. При прохождении через ионо- сферу радиосигналов, излучаемых на двух несущих частотах: f1 ≈ 1,6 ГГц и f2 ≈ 1,2 ГГц с космических аппаратов (КА) СРНС, они испытывают разные задержки Ад и Ад, а также изменения фаз А^ и Ад2 . Это позволяет непрерывно опре- делять по навигационным измерениям двухчастотного приемника СРНС значение ПЭС ионосферы NT ~ (А^2 - А^) вдоль радиотрассы «КА СРНС – приемник СРНС» в любой момент t на выходе приемника.

КА СРНС

КА ССС КА СРНС

Дг/^. ^₂(t), ^ф_l(t), &ф₂(1)

Приемник ССС

^^NTM(t)

РошО)

Блок определения СКО мел комасипабных вариаций ПЭС ионосферы прогнозирования помехоустойчивости ССС

N,(0

Двухчастотный приемник СРНС

Рис. 1. Принцип построения и работы комплекса прогнозирования помехоустойчивости ССС по данным мониторинга ионосферы

Следует отметить, что временные задержки А г_1;2 и изменения фаз А^_;2 навигационных сигналов измеряются двухчастотным приемником NovAtel GPStation-6 непрерывно с интервалом дискретизации А/_д = 0,02 с, а величина ПЭС ( N_T ) – с интервалом 1 с [4]. Однако представляется возможным использовать протокол RANGE в составе платформы NovAtel GPStation-6 так, чтобы и величина ПЭС ионосферы N_T(0 определялась непрерывно во времени с минимальным интервалом дискретизации А?_д = 0,02 с, достаточным для измерения вариаций мелкомасштабных неоднородностей А^См О .

На рис. 1 показаны положения перемещающегося КА СРНС в моменты времени ^и ?! , соответствующие пересечению радиотрассой области возмущений ионосферы с интенсивными мелкомасштабными неоднородностями электронной концентрации А^ТМ в слое F. Исходными данными для блока определения СКО мелкомасштабных вариаций ПЭС ионосферы ^n™ (Р является временной ряд ПЭС ионосферы NT (t) , формируемый на выходе двухчастотного при- емника СРНС. Изменения О" Д/У составляют основу для вычисления изменений во времени индекса мерцаний 54(0 = ^(^д™(0) и параметра Накагами в радиоканале спутниковой навигации m

От приемника ССС к блоку прогнозирования помехоустойчивости поступают сведения о среднем энергетическом отношении «сигнал/шум» на входе приемника (^2) , несущей частоте /о и угле прихода сигналов в текущий момент времени ^с (0 . Потребность в этих сведениях будет обоснована ниже.

Определение СКО мелкомасштабных вариаций ПЭС в ионосфере

В общем случае временной ряд ПЭС на выходе двухчастотного приемника GPStation-6 будет представлять собой аддитивную смесь несколь- ких составляющих [5-6]. Известно [5], что в общем случае формируемый на выходе двухчастотного приемника СРНС временной ряд вариаций ПЭС представляет собой аддитивную смесь нескольких составляющих:

NT(t) = NMO + OTKp(t) + АЖср(О + (1) WO + WO, где NT0(t) – ПЭС однородной ионосферы (фона); ∆NТкр(t), ∆NТср(t) и ∆NТм(t) – вариации ПЭС, обусловленные ионосферными неоднородностями крупного (100…1000 км), среднего (10…100 км) и мелкого (100…1000 м) масштабов соответственно; ∆NТш(t) – вариации, обусловленные шумовой погрешностью измерений двухчастотного приемника.

В качестве примера на рис. 2 приведен отрезок временного ряда вариаций ПЭС N_Т(t) ионосферы (выраженных в единицах измерения TECU = 10¹⁶ м^-2) за время наблюдения 1 ч (с 11:30 до 12:30), полученный 25.07.2015 г. на выходе двухчастотного приемника NovAtel GPStation-6. Этот приемник размещался на базе Северо-Кавказского федерального университета (г. Ставрополь) и принимал навигационные сигналы от КА № 25 СРНС GPS.

Анализ рис. 2 показывает, что в интервале с 11:50 до 12:05 наблюдаются значительные колебания ПЭС (единицей измерения ПЭС является TECU = 1016 м-2) с характерным периодом до 2030 с. Такие колебания обусловлены мелкомасштабными неоднородностями электронной концентрации [5]. Подионосферная точка КА № 25 GPS в период времени 11:50 … 12:05 перемещалась практически вдоль широты 47° по траектории г. Каховка – г. Бердянск (положения этого КА СРНС в моменты времени 11:50 и 12:05 показаны на рис.1).

Согласно методике [6] на первом этапе обработки временного ряда вариаций ПЭС ионосферы (1) на трассе «КА СРНС – двухчастотный приемник СРНС» (полученного на выходе этого приемника) осуществляется сглаживание N_Т(t) скользящим временным окном с шириной равной t_сгл1 = 60 с. Полученный в результате сглаживания тренд <∆N_Т> = N_Т0(t) + ∆N_Ткр(t) + ∆N_Тср(t) вычитается из исходного ряда (1). В результате ряд ПЭС принимает вид комбинации мелкомасштабных вариаций и шумов: N_Т(t) = ∆N_Тм(t) + ∆N_Тш(t). На втором этапе обработки полученный ряд вариаций ПЭС N_Т(t) = ∆N_Тм(t) + ∆N_Тш(t) сглаживается простой скользящей средней с окном t_сгл2 = 0,1 с, что обеспечивает удаление шумовой составляющей ∆N_Тш(t).

Рис. 2. Вариации полного электронного содержания ионосферы

В итоге после двух этапов обработки из ряда (1) формируется временной ряд вариаций ПЭС N_Т(t) = ∆N_Тм(t), обусловленных только мелкомасштабными неоднородностями ионосферы, представленный на рис. 3.

На третьем этапе обработки из этого ряда (см. рис.3) формируются результаты расчета СКО вариаций ПЭС в мелкомасштабных неоднородностях ионосферы ^JO’ представленные на рис. 4. Из рис. 4 видно, что в период с 11:50 по

12:05 значение СКО мелкомасштабных вариаций ПЭС изменяется с 0,03 TECU до 0,18 TECU.

Прогнозированиепомехоустойчивости ССС

В блоке прогнозирования помехоустойчивости ССС на основе полученных значений СКО мелкомасштабных вариаций ПЭС ионосферы ^лг™ (О определяются изменения во времени вероятности ошибки ^ОШ (О при приеме инфор-

Рис. 3. Временной ряд мелкомасштабных вариаций

Рис. 4. Изменение среднеквадратического отклонения мелкомасштабных вариаций полного электронного содержания ионосферы

мационных символов навигационного сообщения (см. рис. 1). Заметим, что согласно рис.1 на трассах спутниковой навигации под углом ^н (О и спутниковой связи под углом «с (О длина трансионосферного участка распространения радиоволн (РРВ) будет различной. Зная СКО мелкомасштабных вариаций ПЭС ионосферы на трассе спутниковой навигации ^лг™ (0 > можно определить СКО мелкомасштабных вариаций ПЭС ионосферы, в радианах, на трассах спутниковой связи Сдг (0 , проходящей через ту же область возмущенной ионосферы, как

(А = с.т (AJcosec aAOcosec а„(А.

т^Цо) ^{x 7} N V^х/ n ^x /

При одинаковых углах трансионосферного РРВ на трассах спутниковой навигации и связи «c(0 = «нО) обеспечивается равенство и СКО мелкомасштабных вариаций ПЭС на этих трассах °"а7-,_М|с| (О ^ДТ™ (0 "

Известно [7-8], что для описания процесса распространения радиоволн через неоднородную ионосферу в ССС используется метод фазово- го экрана, устанавливающий зависимость СКО флуктуаций фазового фронта волны на выходе неоднородной ионосферы ^Ф от СКО мелкомасштабных вариаций ПЭС ионосферы :

(7_&Nt = 80,8л-о-_дл,_г /с/₀, (2)

где c – скорость света [м/с], /о – несущая частота [Гц], 80,8 – размерный коэффициент м³/с². На основе полученного временного ряда СКО мелкомасштабных вариаций ПЭС ^ (О можно получить временной ряд СКО флуктуаций фазового фронта волны на выходе неоднородной ионосферы (в радианах) как

^"^ (О 80,8^ZY7"д/v^^ (0 ^ ^/о * (3)

На рис. 5 представлены результаты расчета согласно (3) изменений во времени СКО флуктуаций фазового фронта волны с несущей частотой /о ® 1 ГГц (широко используемой в ССС) на основе данных об СКО мелкомасштабных вариаций ПЭС ионосферы, приведенных на рис. 4 (то есть при ^,«,.,(0 = ^(0).

Рис. 5. Изменение во времени среднеквадратического отклонения флуктуаций фазового фронта волны на выходе неоднородной ионосферы

Согласно рис. 5 в период с 11:50 по 12:05 значение СКО мелкомасштабных вариаций ПЭС изменяется с 0,03 рад до 0,52 рад. Известно [7-8], что СКО флуктуаций фазового фронта волны на выходе неоднородной ионосферы ОД определяет глубину замираний принимаемого сигнала в трансионосферном радиоканале. Эти замирания характеризуются индексом мерцаний. Последний при выполнении условия дальней зоны связан с О_ф зависимостью [9]

^^[l-expC^^r5.   (4)

В соответствии с выражениями (4) и (3) можно получить зависимость от ^GM . (О изменения во времени индекса мерцаний в трансионосферном радиоканале:

54=[1-ехр(-2^)Г =

• — ^ СХр [ 161 6 7Г ^"ддг    (^) / ^/о ] 1 ^’^ (5)

На рис. 6 приведен временной ряд значений индекса мерцаний, рассчитанных согласно (5) для несущей частоты /о ~ 1 ГГц по данным измерений СКО флуктуаций ПЭС (см. рис. 4).

Согласно рис. 5 в период с 11:50 по 12:05 значение СКО индекса мерцаний принимаемых сигналов в трансионосферном канале связи может достигать S₄ ~ 0,6; что свидетельствует о возникновении в спутниковом радиоканале достаточно глубоких замираний (когда т = S₄² « 3, то есть близких к релееевским).

В соответствии с выражением (5) можно получить зависимость от ^n™ (О изменения параметра m во времени в трансионосферном радиоканале:

m(t) = S;2(0 = Vl-exp(-2(yX -

• ~ 1^ — охр [—161,6 л* О"дуу )(О^^/о]} 1* (6)

Величина m-параметра распределения Накага-ми, зависящая от индекса мерцаний *^4 в трансио-

Рис. 6. Изменение во времени индекса мерцаний сигналов в трансионосферном канале связи

носферном канале связи, сильно влияет на помехоустойчивость ССС. Последняя при некогерентном оптимальном приеме ортогональных (в усиленном смысле) сигналов описывается выражением вида [3]

где (h²} = (E)/N₀ = h² = E / N_o – отношение средней энергии замирающего сигнала (^) к спектральной плотности мощности шума N₀на входе приемника, равное соответствующему энергетическому отношению «сигнал/шум» h² = E/N_n при отсутствии замираний.

В соответствии с выражениями (6)-(7) можно определить изменение во времени помехоустойчи- вости спутниковых систем связи и навигации как рош (0 = | exp

^4^

В соответствии с изменением во времени индекса мерцаний S д в спутниковом канале связи (см. рис. 6) и выражением (8) на рис. 7 приведены результаты прогнозирования помехоустойчивости ССС, использующих несущую частоту /о ~ 1 ГГц, при достаточно большом отношении «сигнал/шум» (ОСШ) на входе приемника (/г²) =35 дБ.

Анализ рис. 7 показывает, что в период возмущений ионосферы с 11:50 до 12:05, сопровождаемых ростом мелкомасштабных неоднородностей электронной концентрации и вариаций ПЭС (см. рис. 3-4), вероятность ошибочного приема информационных символов в ССС с несущей частотой /о * 1 ГГц может возрастать до величины

Р_ош «2-Ю-⁴, более чем на порядок превышающей допустимое значение Рош доп ^ ^ *

На рис. 8, по аналогии с рис. 6 и в соответствии с выражениями (6)-(8), представлены результаты прогнозирования помехоустойчивости ССС, использующих более высокую несущую частоту f₀«2 ГГц при обычном (с учетом энергетического запаса спутниковой радиолинии) отношении «сигнал/шум» на входе приемника ССС Ql"^= 20 дБ.

Анализ рис . 8 показывает, что на интервале возмущений ионосферы с 11:50 до 12:05 вероятность ошибочного приема сигналов в ССС с несущей частотой К «2 ГГц и отношении «сигнал/шум» на входе приемника (^²) = 20 дБ может возрастать до величины Р_ош « 5-10 и превышать допустимое значение Р_ошдоп = Ю"⁵ в течение 13 мин.

На рис. 9 представлены результаты прогнозирования помехоустойчивости ССС, использующих несущую частоту /о « 4 ГГц при отношении «сигнал/шум» на входе приемника ССС (й2) =20 дБ.

Анализ рис. 9 показывает, что в период возмущений ионосферы с 11:50 до 12:05 вероятность ошибочного приема сигналов в ССС с несущей частотой /о « 2 ГГц при (/г²) = 20 дБ может возрастать до величины ^ОШ ~ ^ ^ > превышающей допустимое значение ^ОШ ДОП ^О на два порядка.

Анализ приведенных на рис. 7-9 изменений во времени вероятности ошибочного приема сигналов в ССС при возмущениях ионосферы показывает, что "^ош СО слабо зависит от реализуемого значения ^0 на входе приемника и определяется индексом мерцаний S₄² (0 = т(О, изменения которого зависят согласно (6) от отношения

Рис. 8. Изменение во времени вероятности ошибки в ССС при несущей частоте J_o ~ 2 ГГц и //?² \ = 20 дБ

Рис. 9. Изменение во времени вероятности ошибки в ССС при несущей частоте J_o х; 4 ГГц и и?² \ = 20 дБ

СКО мелкомасштабных вариаций ПЭС ионосферы О"л„ (0 к выбранной несущей частоте передаваемого сигнала ^ .

Структура комплекса прогнозирования помехоустойчивости ССС

Алгоритм прогнозирования помехоустойчивости ССС по результатам измерения СКО мелкомасштабных вариаций ПЭС О' Ал заключается в последовательности следующих операций.

• 1.    Получение временного ряда вариаций ПЭС N_r (7) – см. рис. 2, с помощью двухчастотного приемника СРНС.

• 2.    Сглаживание ряда вариаций ПЭС N_T (0 методом простой скользящей средней с окном 4гл1 = ⁶⁰ С и удаление из исходного ряда (1) полученного в результате сглаживания тренда (N_T (7)^ = ZW _т0 (?) + АЛ^_/кр (7) + AN_TCP (7).

• 3.    Получение ряда мелкомасштабных вариаций ПЭС ^^тм (0 – см. рис. 3, путем сглаживания, полученного на предыдущем шаге ряда N_T (0 — АА_ш (/) + ZSJVjpy (t) методом простой скользящей средней с окном ^сгл2 “0,1 С.

• 4.    Вычисление элементов временного ряда СКО мелкомасштабных вариаций ПЭС ° АЛ™ ^ согласно рис. 4, на интервале ^СКО = 60 с.

• 5.    Вычисление элементов временного ряда значений индекса мерцаний ^ 4 (см. рис.6);

• 6.    Вычисление элементов временного ряда значений вероятности ошибки ^_ош« при приеме информационных символов (см. рис. 7-9).

Структурная схема комплекса, реализующего описанный алгоритм с учетом описанных выше операций, представлена на рис. 10.

Выводы

Таким образом, разработана структурная схема комплекса, позволяющего прогнозировать помехоустойчивость систем спутниковой связи в реальном масштабе времени по результатам мониторинга среднеквадратического отклонения мелкомасштабных вариаций полного электронного содержания ионосферы на выходе двухчастотного приемника СРНС. Практическая ценность полученных результатов состоит в том, что они базируются на экспериментальных данных мониторинга ионосферы с помощью двухчастотного приемника NovAtel GPStation-6 и указывают на возможность образования в средних широтах областей с повышенным значением СКО мелкомасштабных вариаций электронной концентрации, которые могут вызывать глубокие мерцания принимаемых сигналов и недопустимое снижение помехоустойчивости систем спутниковой связи в течение 13 мин. Эти результаты могут быть использованы для адаптации частотно-временных параметров спутниковых систем связи к изменению условий трансионосферного распространения радиоволн с целью повышения помехоустойчивости в условиях ионосферных возмущений типа образования мелкомасштабных неоднородностей электронной концентрации.

Рис. 10. Структурная схема комплекса прогнозирования помехоустойчивости ССС по данным трансионосферного зондирования