61
УДК 621.391, 621.396, 621.369
В. А. Пахотин, С. В. Молостова, В. А. Бессонов
ФИЛЬТРАЦИЯ РАДИОИМПУЛЬСОВ МЕТОДОМ МАКСИМАЛЬНОГО ПРАВДОПОДОБИЯ
Представлено решение задачи разрешения двух радиоимпульсов методом максимального правдоподобия. Получены выражения для оптимальной обработки совокупности двух радиоимпульсов в области их ортогональности и в области, когда радиоимпульсы неортогональны, а также выражения для дисперсии параметров радиоимпульсов (дисперсии Рао-Крамера). Приведены результаты модельных расчетов, иллюстрирующих возможности разрешения двух радиоимпульсов методом максимального правдоподобия. Дана оценка области решения задачи разрешения двух радиоимпульсов методом максимального правдоподобия.
This article presents a solution to the problem of resolution of two radio impulses through the maximum likelihood method. The authors have obtained expressions for the optimum processing of a set of two radio impulses in the domains of their orthogonality and non-orthogonality, as well as expressions for the dispersion of radio impulse parameters (Rao-Cramer dispersion). The article presents the results of model calculations showing the possibilities of resolving two radio impulses using the maximum likelihood method. The authors assess the area of solution of the problem of resolving two radio impulses using the maximum likelihood method.
Ключевые слова: теория оптимального приема, метод максимального правдоподобия, корреляционная функция, область ортогональности, неортогональности сигналов, дисперсия Рао-Крамера.
Key words: optimum reception theory, maximum likelihood method, correlation function, domain of orthogonality, non-orthogonality of signals, Rao-Cramer dispersion.
© Пахотин В. А., Молостова С. В., Бессонов В. А., 2015
Вестник Балтийского федерального университета им. И. Канта. 2015. Вып. 4. С. 61 — 67.
62
Введение
В области обработки сигналов одной из важных задач становится разделение (фильтрация) радиоимпульсов, содержащихся в принятом сообщении, с оценкой времени приема и амплитуды. Для решения этой задачи используется корреляционная обработка. Однако она имеет релеевское ограничение, позволяющее разделять лишь ортогональные радиоимпульсы, то есть такие радиоимпульсы, различие по времени приема которых более длительности радиоимпульса. В системах локации, если отраженные от двух близких по дальности целей радиоимпульсы являются неортогонльными, то принимается решение о приеме сигнала от одной цели. В системах связи с временным разделением каналов, в системах телевещания для приема лишь ортогональных радиоимпульсов создают между ними защитные интервалы, в результате чего пропускная способность каналов связи, каналов телевещания уменьшается.
В настоящей работе задача разделения (фильтрации) радиоимпульсов решается на основе положений теории оптимального приема [1; 2] методом максимального правдоподобия. В этом случае возможно разделение неортогональных радиоимпульсов с оценкой амплитуд, фаз, времен приема. В работе приводятся основы теории, результаты модельных расчетов при решении задачи разделения двух неортогональных радиоимпульсов, содержащихся в принятом сообщении. На основе выражения для дисперсии Рао-Крамера оценивается область, в которой возможно решение задачи неортогональных радиоимпульсов.
Основы теории
Запишем принятое сообщение у(I) в комплексном вице:
у(1) = -к) + Ъгё^-ч) + иш (I), (1)
где и\,и 2 — комплексные амплитуды двух радиоимпульсов; С — круговая частота; 11 ,12 — время приема первого и второго радиоимпульсов; и ш (I) — аддитивный шум с нормальным распределением, средним значением, равным нулю, дисперсией о и интервалом корреляции -к .
На основании (1) запишем логарифм функции правдоподобия [1; 2]:
2
1п( Ь(Х') = -0-1
2о -к о
С(.1 -11), Т] „ С( I-12)
у(I) -Т ёС('+ Т2 е
Ж, (2)
где X' — вектор оцениваемых параметров двух радиоимпульсов; штрихами обозначены оцениваемые параметры радиоимпульсов.
Логарифм функции правдоподобия представляет собой плотность распределения, максимум которой дает наиболее правдоподобную оценку параметров сигнала. В связи с этим, дифференцируя (2) по амплитудам 1!\ , и2 и приравнивая дифференциалы к нулю, можно получить систему двух уравнений. Решая эту систему относительно оценочных амплитуд, получим
и 1 =
1 _
и
и ' у(')е-'а('-'2) _ и у(')е-ш( 1<)
и 2 =-
(3)
1 - и
где черта сверху означает интегрирование по '; и = е
_ 'Ю{'2-';)
норми-
рованный коэффициент корреляции двух радиоимпульсов; и — комплексно-сопряженный коэффициент корреляции.
Оценочные значения амплитуд и1 , и2 в (3) зависят от неизвестных оценочных времен приема '[, '2 . В связи с этим подставим выражения из (3) в функционал правдоподобия
аС;, '2) = {
у(') - и1 ;еш(+и2 2е
гю(-'2 )
й'..
(4)
Этой подстановкой исключается зависимость функционала от амплитуд и1 ,и2, и функционал (4) будет представлять собой поверхность в двухмерном пространстве '[, '2. Максимум этой поверхности
определяет решение — оценочные значения амплитуд V; ,и2 и времен приема '[, '2. Значение функционала (4) в минимуме определяет
дисперсию шума с2.
Дадим оценку дисперсий амплитуд радиоимпульсов в минимуме
поверхности функционала при и1 =и1, и2 =и2, '[ ='2 = '2. Для этого получим элементы информационной матрицы Фишера согласно выражению [1; 2]
Jy =-М
( й 2(1п( 1(Л))) ^
йЛйЛ.
' 1 У
' = 1,2; 1 = 1,2,
63
где М — оператор математического ожидания; Л' = и1 , ЛЛ = и2 .
2
64
Выполнив дифференцирование, получим информационную матрицу Фишера:
Т =
т
(
о2-,.
1 Я
*
Я 1
\
(6)
Диагональные элементы матрицы О, обратной матрице Фишера, определяют дисперсии амплитуд Ои1, Ои 2 :
О , = О 2 =-
и1 и 2
о
N (1 -
Я
(7)
где N =
т
количество некоррелированных отсчетов шума.
Выражение для дисперсий (7) при Я = 0 определяет дисперсию
амплитуд ортогональных радиоимпульсов ( О о ); следовательно, относительная дисперсия амплитуд будет
О, 1
д..
1 -
Я
(8)
При изменении модуля коэффициента корреляции в пределах 0-0,9 относительная дисперсия увеличивается незначительно (0-7 дБ).
В области, где
Я
> 0,9, относительная дисперсия резко увеличивается.
Таким образом, областью полученного решения можно считать область, в которой модуль коэффициента корреляции меняется в пределах 0-0,9.
Результаты модельных расчетов
В модельных расчетах приняты следующие значения параметров двух радиоимпульсов: амплитуды и1 = 1, и2 = 0,8, начальные фазы <р1 = 30°, р2 = 230°, время приема Ь =80 мс, Ь = 90 мс длительность радиоимпульсов Т = 80 мс, частота / = 100 КГц, интервал между отсчетами -к = 1 мс.
На рисунке 1 показан вид суперпозиции двух радиоимпульсов при коэффициенте корреляции |.к| = 0,875 . При использовании классического корреляционного анализа корреляционная функция имеет сложный вид (рис. 2) и не позволяет дать оценку времен приема и амплитуд двух радиоимпульсов.
т
к
2
Рис. 1. Суперпозиция Рис. 2. Огибающая функции
двух неортогональных радиоимпульсов. корреляции при R = 0,875
R = 0,875
65
На рисунке 3 показаны три корреляционные функции. Корреляционная функция (ряд 3) получена классическим корреляционным анализом. Корреляционные функции (ряд 1 и 2) получены с использованием метода максимального правдоподобия. Максимумы этих корреляционных функций дают оценки амплитуд и времен приема двух радиоимпульсов. Данный результат показывает, что на основе метода максимального правдоподобия возможно создание адаптивного корреляционного фильтра, который фильтрует (разделяет) радиоимпульсы с оценкой времен приема и амплитуд.
Рис. 3. Огибающие корреляционных функций при R = 0,875
На рисунке 4 приведена зависимость оценок времени приема , /2 при изменении времени приема второго радиоимпульса /2. Данный рисунок позволяет оценить разрешающую способность метода максимального правдоподобия. Классический корреляционный анализ разрешает радиоимпульсы, если |/2 - > 80 мс (в этом случае /2 > 160 мс).
Метод максимального правдоподобия позволяет разрешать радиоимпульсы, если |/2 - /¡| > 2 мс . Следовательно, разрешающая способность может быть увеличена в ~ 40 раз.
66
Значения Т2 в ис.
Рис. 4. Зависимости времени приема от Т2. ис/иш = 28 дБ
Заключение
В настоящей работе представлено решение задачи разрешения двух радиоимпульсов методом максимального правдоподобия. Изложены элементы теории. Приведены результаты модельных расчетов, иллюстрирующие возможности метода. Основной результат работы заключается в следующем. Метод максимального правдоподобия позволяет разрешать неортогональные радиоимпульсы при коэффициенте взаимной корреляции, достигающем значения 0,9 при отношении сигнал/шум более 10 дБ. Это позволяет в системах связи с временным разделением радиоимпульсов увеличить пропускную способность каналов связи.
Список литературы
1. Перов А. И. Статистическая теория радиотехнических систем : учебное пособие для вузов. М., 2003.
2. Пахотин В. А., Бессонов В. А., Молостова С. В., Власова К. В. Курс лекций для радиофизических специальностей. Калининград, 2008.
3. Пахотин В. А., Пахотина К. В., Жукова Н. В. Метод обработки данных, полученных при приеме ионосферных сигналов // Геомагнетизм и аэрономия. 2004. Т. 44, № 4. С. 511-517.
4. Пахотин В. А., Бессонов В. А., Иванова С. В. и др. Частотное разделение сигналов в области высокой корреляции несущих частот // Известия вузов России. Радиоэлектроника. 2005. Вып. 4. С. 117 — 123.
5. Власова К. В., Пахотин В. А., Власов А. А., Брух Я. Р. Решение задачи разрешения по дальности в локации методом максимального правдоподобия // XIV Международная научно-практическая конференция «Радиолокация, навигация, связь». Воронеж, 2007. С. 344 — 353.
Об авторах
Валерий Анатольевич Пахотин — д-р физ.-мат. наук, проф., Балтийский федеральный университет им. И. Канта, Калининград.
E-mail: [email protected]
Светлана Валерьевна Молостова — канд. физ.-мат. наук, доц., Балтийский федеральный университет им. И. Канта, Калининград.
E-mail: [email protected]
Владимир Александрович Бессонов — канд. физ.-мат. наук, доц., Балтийский федеральный университет им. И. Канта, Калининград.
E-mail: [email protected]
About the authors
Prof. Valery Pakhotin, I. Kant Baltic Federal University, Kaliningrad.
E-mail: [email protected]
Dr Svetlana Molostova, Associate Professor, I. Kant Baltic Federal University, Kaliningrad.
E-mail: [email protected]
Dr Vladimir Bessonov, Associate Professor, I. Kant Baltic Federal University, Kaliningrad.
E-mail: [email protected]
67