УДК 621.391, 621.396, 621.369
К. В. Власова, М. А. Никитин,
А. С. Чугайнов, А. В. Кочмарский
ОЦЕНКА ПАРАМЕТРОВ ИОНОСФЕРНОГО СИГНАЛА
Представлена методика оценки параметров двухлучевого ионосферного сигнала методом максимального правдоподобия. Приведены результаты модельных экспериментов, а также расчеты реальных ионосферных сигналов на трассе Москва — Калининград.
The technique of an estimation of parameters of a dual-beam ionospheric signal is presented by a method of the maximum credibility. Results of modeling calculations and results of calculations of real ionospheric signals on line Moscow — Kaliningrad are resulted.
Ключевые слова: ионосфера, распространение радиоволн, метод максимального правдоподобия, фазовая плоскость.
Key words: ionosphere, distribution of radio-waves, method of the maximum credibility, phase plane.
При длительностях зондирующих сигналов менее секунды в точку приема приходит несколько радиоимпульсов с временными задержками ~ 400 мкс. В этом случае ионосферный сигнал становится сложным. Он состоит из ряда радиоимпульсов, отраженных от разных областей ионосферы, с частичными наложениями во времени. Методика оценки параметров такого ионосферного сигнала дана в статье «Модель ионосферного сигнала с мультипликативной помехой», опубликованной в данном журнале (см. с. 85). В нашей работе продолжен анализ ее возможностей. Представлены результаты модельных экспериментов и расчетов параметров ионосферного сигнала на трассе Москва — Калининград. Полученные при этом данные отображены на фазовой плоскости. Показана возможность фазового уточнения частоты сигнала и начальной фазы радиоимпульсов. В отличие от корреляционного анализа начальная фаза оценивается на интервале радиоимпульса, а не в точке начала радиоимпульса.
В статье «Модель ионосферного сигнала с мультипликативной помехой» представлена методика оценки параметров ионосферного сигнала, состоящего из ряда радиоимпульсов, отраженных от разных областей ионосферы. Продолжим анализ возможностей данного метода на примере двухлучевого ионосферного сигнала. Частота О оценивается выражением
Вестник Балтийского федерального университета им. И. Канта. 2012. Вып. 4. С. 78 — 84.
^]{Уп+ 2 + Уп )УЯ + 1 /у + у \у
С' _ п=1 _ \ Уп+2 Уп) Уп+1 , (1)
N 2
£уп+12 Уп+1
п=1
где С' = ехр(/'шА?).
Для оценки амплитуды и фазы сигнала запишем функционал правдоподобия:
м , „ 2
Л1п =Е| уп+т - и' еХР (Шп+т )| . (2)
т_0
Индекс т выполняет операцию суммирования на интервале обработки информации Л1 (М _ 20, Лt _ 20 мкс). Индекс п определяет
скользящий характер обработки. Дифференцируя (2) по и' и приравнивая дифференциал к нулю, получим оценочное значение комплексной амплитуды и':
иП = Уп+т ехр (-міп+т). (3)
Черта сверху означает суммирование на интервале Аґ = 20 мкс. Выражение (3) справедливо для всех трех областей ионосферного сигнала за исключением переходных, где нет соответствия правой и левой частей функционала. На рисунке 1 показана зависимость оценочной фазы сигнала от времени. На интервале времени, где существует сигнал, дисперсия фазы небольшая. Отмечаются различия фаз первой, второй и третьей областей сигнала, а также регулярный наклон фазовой зависимости, связанный с неточностью оценки частоты ю. Наличие этого наклона позволяет провести коррекцию оценки частоты по фазе. Для этого из соотношения
р( ґ ) = Аюґ + р0, (4)
определяется погрешность частоты Аю и добавляется к первоначальной оценке ю' = ю + Аю .
В результате частота ю' существенно уточняется, наклон линий сводится к нулю и фазы всех трех областей оказываются независимыми от времени (рис. 2). Таким образом, имеется возможность оценить начальные фазы ионосферного сигнала в моменты времени ґ01, ґ02, ґ01 + Т . Фазы трех областей сигнала р1, р2, р3 могут быть представлены следующими выражениями:
Р =Роі,
Р2 = р01 + ю(ґ02 — Ґ01) — р02 , (5) Рз =Р02 +ю(Ґ01 + Т),
79
где р01, р02, t01, t02 — начальные фазы и время приема первого и второго радиоимпульсов. При выполнении операции «коррекция часто-
ты» в выражениях (5) частота полностью исключается и фазы трех областей сигнала р1, р2, р3 однозначно оказываются связанными с начальными фазами радиоимпульсов р , р02.
80
(в
е
200
200
оостіонмтм-ілюь-мтон
ОО^ННгННННгНгНгНгНГчІГ^ГЧ
Рис. 1
м
гв
е
Рис. 2
Амплитуды трех областей ионосферного радиосигнала оцениваются по выражению (3). На рисунке 3 показана зависимость амплитуд от времени.
Первая область дает возможность оценить амплитуду первого и1 радиоимпульса, третья — второго и2, а вторая — амплитуду и12:
и 12 = д/ и,2 + и22 + 2ир2 008 (02 — Р01) . (6)
В результате параметры двухлучевого ионосферного сигнала могут быть оценены.
81
Статистика получаемых решений показана на рисунке 4, изображающем фазовую плоскость, на которой точками представлены квадратурные компоненты U cos (^) и U sin (^) . В результате совокупность
точек в области нуля определяет распределение шума в отсутствии сигнала. При наличии сигнала отображаются три области с соответствующими амплитудами и фазами. Отношение «сигнал/шум» в данном случае равно 27 дБ. Статистика параметров сигнала (амплитуд и фаз) при уменьшении данного отношения до 13 дБ показана на рисунке 5.
Рис. 4
82
Рис. 5
Дисперсия амплитуды и фазы увеличилась, однако вполне возможно выделение областей сигнала и оценка средних значений амплитуд и фаз.
Фазовая плоскость является основой для обнаружения сложного, многолучевого импульсного ионосферного сигнала и анализа его структуры. Она может быть преобразована в двумерную поверхность плотности распределения амплитуды фазовых точек сигнала. Кроме того, на фазовой плоскости совокупность радиоимпульсов в комплексной форме имеет простой вид. Последовательные отсчеты сигнала в первой области можно записать как
U1n = U01 eXP O' (0tn + P01)) , f01 < fn < t02 .
Во второй области
U2n = U0ieXP O' (n + 0t02 +P01)) + U02eXP (' (0tn +P02 )) , f02 < fn < f01 + T •
В третьей —
U 3n = U 02 exP (' (atn +P02)), t01 + T < tn < t02 + T .
Если исключить частоту О , то в комплексном виде в первой, второй и третьей областях сигнала выражения будут иметь следующий вид:
U1 = U 01 exP (р01),
U2 = U01 eXP (P01 ) + U02 eXP (p02 ) , (7)
U3 = U02 eXP (P02 ).
Это основа для анализа сложной многолучевой структуры ионосферного сигнала. При наличии трехлучевой структуры выражения (7) запишутся в виде
и1 = U 01 exp (p),
U2 = U01 exp (p ) + U02 exp (p ) ,
Uз = U01 ЄХР (p01 ) + U02 ЄХР (Р02 ) + U03 ЄХР (p03 ) , (8)
U4 = U02 ЄХР (Р02 ) + U03 ЄХР (Р03 ) ,
U5 = U03 ЄХР (p03 )
Решая данную систему относительно комплексных амплитуд радиоимпульсов, можно оценить структуры ионосферного сигнала.
В таблице представлены оценки параметров двухлучевого сигнала в зависимости от отношения «сигнал/шум». В первой ее строке даны модельные значения параметров двухлучевого ионосферного сигнала. В остальных строках представлены оценочные данные сигнала в зависимости от отношения «сигнал/шум».
Результаты модельных исследований структуры сигнала
83
Сигнал/ шум, дБ Ui, в Р, ° *01, мс U 2 , В Р2 , ° *02 , мс /, кГц
Модельные данные 1 10 1 0,5 90 1,2 215
33 1,0089 9,58 1 0,499 89,96 1,2 215,0015
27 1,012 9,94 1 0,498 89,8 1,2 215,0007
21 1,019 10,7 1 0,494 89,8 1,198 214,9991
18 1,027 10,1 1 0,489 88,9 1,193 215,0013
15 1,035 9,75 1 0,485 87,6 1,191 215,0027
13 1,04 9,65 1 0,48 85,9 1,19 215,0041
12 1,04 -1,009 1,01 0,479 64,5 - 215,0322
33 1,0089 9,58 1 0,499 89,96 1,2 215,0015
Оценим с точки зрения теории оптимального приема погрешности параметров двухлучевого ионосферного сигнала.
Согласно теории [1; 2], дисперсии амплитуды Ои, начальной фазы
частоты ^ и времени приема оцениваются выражением Рао — Крамера
^-2 _2 _2 2
= 2 °--------, Ои =—, = -——, А =——-----------, (9)
|и| ЫТ2 (2^)2 Ы 0 |и| N 0 |и| Ыо2
где —2 — дисперсия шума; N — количество некоррелированных по шуму отсчетов; Т — длительность интервала обработки; О — круговая частота.
Принимая — = 0,1, / = 215 кГц, Т = 20 мкс, получим следующие среднеквадратичные отклонения: СКО^ = 0,006, СКОи^ = 0,006, СКО%1 = 1,14°, СКО^ = 2,29°, СКО/ = 0,44, СКО0 = 5 -10_3. Данные погрешности параметров сигнала, в принципе, соответствует оценкам Рао — Краме-
84
ра. Исключением является время приема первого радиоимпульса 001
которое имеет погрешность ~3 мкс, что не соответствует оценке Рао — Крамера. Однако это близко к оценке Вудворда, определяющей дисперсию времени приема по огибающей радиоимпульса:
D0 = °2Т . (10)
00 ~ 2 , V 2
р| N(2^)2
CKO0 , рассчитанное согласно этой формуле, дает значение 0,2 мкс.
Таким образом, модельные расчеты показали возможность оценки параметров двухлучевого импульсного ионосферного сигнала с высокой точностью. Применение метода оценки параметров ограничено отношением «сигнал/шум» ~13 дБ.
Список литературы
1. Перов А. И. Статистическая теория радиотехнических систем : учебное пособие для вузов. М., 2003.
2. Пахотин В. А., Бессонов В. А., Молостова С. В., Власова К. В. Теоретические основы оптимальной обработки сигналов : курс лекций для радиофизических специальностей. Калининград, 2008.
3. Афраймович Э. Л. Интерференционные методы радиозондирования ионосферы. М., 1982.
Об авторах
Ксения Валерьевна Власова — канд. физ.-мат. наук, доц., Балтийская государственная академия рыбопромыслового флота, Калининград.
E-mail: p_ksenia@mail.ru
Михаил Анатольевич Никитин — д-р физ.-мат. наук, проф., Балтийский федеральный университет им. И. Канта, Калининград.
E-mail: MNikitin@kantiana.ru
Александр Сергеевич Чугайнов — асп., Балтийский федеральный университет им. И. Канта, Калининград.
E-mail: p_ksenia@mail.ru
Алексей Викторович Кочмарский — асп., Балтийский федеральный университет им. И. Канта, Калининград.
E-mail: p_ksenia@mail.ru
About authors
Ksenia Vlasova — PhD, associate professor, Baltic State Academy, Kaliningrad. E-mail: p_ksenia@mail.ru
Mikhail Nikitin — Dr, professor, I. Kant Baltic Federal University, Kaliningrad. E-mail: MNikitin@kantiana.ru
Alexander Chugajnov — PhD student, I. Kant Baltic Federal University, Kaliningrad. E-mail: p_ksenia@mail.ru
Aleksey Kochmarsky — PhD student, I. Kant Baltic Federal University, Kaliningrad. E-mail: p_ksenia@mail.ru