УДК 621.391, 621.396, 621.369
С. В. Молостова, К. В. Власова В. А. Бессонов, И. В. Либерман
МОДЕЛЬ ИОНОСФЕРНОГО СИГНАЛА С МУЛЬТИПЛИКАТИВНОЙ ПОМЕХОЙ
Показано, что нестационарность параметров составляющих ионосферного сигнала приводит к мультипликативной помехе, оказывающей значительное влияние на получение достоверных данных при приеме ионосферных сигналов. Дана оценка дисперсии амплитуд при приеме ионосферных сигналов с мультипликативной помехой. Представлены результаты исследования эффективности методов доплеровской фильтрации при приеме ионосферных сигналов.
It is shown that unsteadiness parameters of components of the ionospheric signal leads to multiplicative noise. It largely affects the reliable data when receiving ionospheric signals. The estimation of the variance of the amplitudes for the reception of ionospheric signals with multiplicative noise. The results of studying the effectiveness of the methods of Doppler filtering when receiving ionospheric signals.
Ключевые слова: ионосфера, распространение радиоволн, теория оптимального приема, мультипликативная помеха, доплеровская фильтрация.
Key words: ionosphere, radio wave propagation, the theory of optimal reception, the multiplicative noise, Doppler filtering.
При ионосферном распространении радиоволн в точку приема приходит несколько лучей, отраженных от разных областей (слоев) ионосферы. Однако каждый отраженный луч не является стационарным во времени: его параметры меняются в зависимости от времени и состояния ионосферы. Амплитуда отдельного луча меняется в зависимости от поглощения, внутримодовой интерференции, факторов фокусировки, дефокусировки, от рассеяния и поляризации. Фазовая длина траектории также варьирует в зависимости от геометрического изменения длины траектории при пересечении траектории луча отдельными неоднородностями ионосферы. Это приводит к изменению начальной фазы отдельного луча и к доплеровскому сдвигу частоты отдельного луча. Таким образом, отдельные лучи ионосферного сигнала имеют как амплитудные, так и частотные изменения. Квазипериод их изменений практически одинаков. Это требует создания новых методов обработки ионосферных сигналов.
Метод максимального правдоподобия, а также классический метод доплеровской фильтрации позволяют решать задачу частотного разделения лучевой структуры сигнала при условии стационарности его параметров. Однако для приема ионосферных сигналов данный метод оказывается неприемлемым. Вследствие нестационарности возникает сильная мультипликативная помеха, которая приводит к существенно-
85
Вестник Балтийского федерального университета им. И. Канта. 2012. Вып. 4. С. 85—89.
му уменьшению максимума функционала правдоподобия. Результаты обработки при этом становятся недостоверными. В настоящее время методов обработки сигналов с мультипликативной помехой практически не существует.
В данной статье представлены анализ влияния мультипликативной помехи на обработку ионосферных сигналов и экспериментальные данные, подтверждающие существенное снижение эффективности обработки ионосферных сигналов.
Примем в качестве модели сигнала при ионосферном распространении радиоволн следующую:
£ (Л Ґ) =^ит (ґ)ехр(7®и (ґ)ґ) , (1)
т=1
где ит (ґ) — амплитуда m — составляющей сигнала; сот (ґ) — круговая частота m-составляющей ионосферного сигнала, также выступает функцией времени.
Временные зависимости амплитуды и частоты отдельных лучей ионосферного сигнала определяют модуляционные изменения, которые в общем случае являются неизвестными. Тогда логарифм функции правдоподобия может быть записан в следующем виде:
2
У(Ґ) -^0'т (Ґ) ЄХР (т (Ґ )Ґ )
(2)
При такой общей постановке решение основных задач при помощи теории оптимального приема оказывается сложным. В этих случаях рекомендуется сначала провести статистическое усреднение логарифма функции правдоподобия [1], а затем — решение задач при помощи теории оптимального приема. Однако при таком подходе оптимальность решения существенно теряется. В качестве нулевого приближения можно рассматривать стационарную модель сигнала
т _ м т
^(Л, 1) = X У™ еХР(г®т?) . (3)
т = 1
Здесь амплитуда и частота сигнала не меняются на интервале обработки. Но при этом не будет соответствия левой и правой частей аргумента логарифма функции правдоподобия. Принятое сообщение тогда будет иметь следующий вид:
м м
У(!) =%0т (I )е'-т'-” + иш (I) = ХО„те'-т' Пт 0) + С) • (4)
где и 0
ляющей сигнала; цт (I) = ^ т(|) е‘(Шщ (1) ш°т )1 — мультипликативная помеха.
^От
При ионосферном распространении радиоволн мультипликативная помеха становится основной. Она определяет дополнительные линии
т=1
то — комплексная стационарная амплитуда и частота m-состав
спектра сигнала и значение максимума функционала. В теоретическом плане методы обработки сигналов при наличии мультипликативной помехи практически в литературе не описаны. Логарифм функции правдоподобия при нулевом приближении можно записать в виде:
1п (ИХ,,) ) = -—- ]
1К 0
У (0 'о т ЄХР И®'о »0
(5)
В этом случае сигнал определен выражением (3). Для увеличения степени соответствия левой и правой частей функционала можно использовать модель с амплитудой, меняющейся по известному закону, но с постоянной частотой. В этом случае логарифм функции правдоподобия будет иметь вид:
2
1п ИХ I) )=-А
2--К о
У(^) - Е и 'о» (і) ЄХР И 'от 1)
(6)
где /т(г) — известная функция изменения амплитуды с различными параметрами.
Однако усложнение вида функций сигнала в правой части приводит к существенному усложнению решения. Покажем это на примере однолучевого и двулучевого сигнала. При однолучевом сигнале логарифм функции правдоподобия запишется в виде:
1пД(Х,г)) = - 1 Цу (г) -и'0 /(г)ехр(гаг) Л. (7)
2<С ТК 0
Принятое сообщение определяется выражением у ( О = иоМ( 0 ехР (' г) + иш,
где ц( г) — функция, определяющая эффект модуляции сигнала за счет изменения ионосферных параметров (мультипликативная помеха).
Дифференцируя (7) по амплитуде и '0 и приравнивая дифференциал к нулю, получим
т
I" у (г) / ( г) ехр (га' г) Л
и' = -0----------------------------------------. (8)
| / 2№
0
Подставляя выражение(8) в (7) получим зависимость функционала правдоподобия от частоты а' и функции /(?):
А = |у( 0| - и '0 у*(г)/(г)ехР ('). (9)
Второе слагаемое определяет функцию корреляции принятого сообщения у ( г) и формы сигнала /( г)ехр(га'() . Если они совпадают, то функционал (9) будет иметь минимальное значение, определяемое дисперсией аддитивного шума с2. Если различия между принятым сообще-
87
2
т=1
т=1
88
С. В. Молостова, К.В. Власова, В.А. Бессонов, И.В. Либерман нием и формой сигнала существенные, тогда коэффициент корреляции
близок к нулевому значению и функционал Д практически равен | у(? Удобнее работать с нормированным обратным функционалом
Д1 =
* г2
(10)
В этом случае функционал Д1 меняется в пределах от 1 до да. При Д1 = 1 принятое сообщение и форма сигнала имеют коэффициент корреляции, равный нулю. Если коэффициент корреляции равен 1, тогда значение функционала будет определяться отношением энергии сигнала к дисперсии аддитивного шума:
2 '
ДЬ
I2Ц)
-2— при МО = I(0 . (11)
а
Таким образом, в зависимости от того, насколько правильно отображает функция I (¿) изменения амплитуды ионосферного сигнала, значение функционала Д1 будут меняться. Следовательно, подбирая функциональную зависимость /(¿) в виде простых функций можно оценивать степень их соответствия реальной модулирующей функции по значению функционала Д1.
Рассмотрим случай приема двух лучей ионосферного сигнала. Принятое сообщение запишем в виде
у(0 = и1^10)ещ( + и2$2( ¿)еШ2‘. (12)
Логарифм функции правдоподобия представим как
- 1 Т 2
1пЦ(Л,0) = _—Цу(0 _и\ тещ _и'2 Шещ\ Л. (13)
2а —к о
Дифференцируя (13) по амплитудам Ц ', и Ц '2 и приравнивая дифференциалы к нулю, получим систему уравнений:
'ь г = и 'ь /2(о+и '2 тте , (14)
у а) гг($ )е_щ ‘ = и 'ь /ко /2(1)е-'Ц'2_щ >+и \ /(Т),
где черта сверху означает интегрирование в пределах 0 Т.
Решением этой систем уравнений будут следующие соотношения:
и \ =
У(о те_,щ‘ /2\1) _ я у (о/2(0е_
1,(0 /22(0 _ Я
и, = У(0/2(0е_щ'2'/ь2(0_|Я У(0/,(0е_
и 2 = --------------------- •
/12(0 /ДО _ Я
2
где R = У1( t)f2(t)e '(°2 0l)t — функция корреляции двух составляющих ионосферного сигнала.
Подставляя решения (15) в функционал правдоподобия (13) получим:
Д(©>'2 ) = |Я0|2 - и \ y*(t) /(t)e° - U \ y*(t) f2(t)e°2t. (16)
Данный функционал зависит от функций f1(t) и f2(t) и от частот 0\, О '2. Он представляет собой поверхность в двумерном пространстве о\, О '2, минимум которой зависит от функций fj(t), f2 (t), от их близости к функциям p(t) и <p2(t).
Таким образом, значение функционала в зависимости от степени соответствия функций p., (t) и f (t), а также p2 (t) и f2 (t) будет меняться в пределах от единицы до отношения сигнал/шум по мощности.
Список литературы
1. Перов А. И. Статистическая теория радиотехнических систем : учебное пособие для вузов. М., 2003.
2. Пахотин В. А., Бессонов В. А., Молостова С. В., Власова К. В. Теоретические основы оптимальной обработки сигналов : курс лекций для радиофизических специальностей. Калининград, 2008.
Об авторах
Светлана Валерьевна Молостова — канд. физ.-мат. наук, доц., Балтийский федеральный университет им. И. Канта, Калининград.
E-mail: p_ksenia@mail.ru
Ксения Валерьевна Власова — канд. физ.-мат. наук, доц., Балтийская государственная академия рыбопромыслового флота, Калининград.
E-mail: p_ksenia@mail.ru
Владимир Александрович Бессонов — канд. физ.-мат. наук, доц., Балтийский федеральный университет им. И. Канта, Калининград.
E-mail: VBessonov@kantiana.ru
Ирина Владимировна Либерман — канд. физ.-мат. наук, ст. преп., Балтийский федеральный университет им. И. Канта, Калининград.
E-mail: VPakhotin@kantiana.ru
About authors
Svetlana Molostova — PhD, associate professor, I. Kant Baltic Federal University, Kaliningrad.
E-mail: p_ksenia@mail.ru
Ksenia Vlasova — PhD, associate professor, Baltic State Academy RF.
E-mail: p_ksenia@mail.ru
Vladimir Bessonov — PhD, associate professor, I. Kant Baltic Federal University, Kaliningrad.
E-mail: VBessonov@kantiana.ru
Irina Liberman — PhD, read., I. Kant Baltic Federal University, Kaliningrad.
E-mail: VPakhotin@kantiana.ru