Оптимальная обработка сигналов с затуханием

Тип работы:
Реферат
Предмет:
ТЕХНИЧЕСКИЕ НАУКИ


Узнать стоимость

Детальная информация о работе

Выдержка из работы

УДК 621. 372. 542
ОПТИМАЛЬНАЯ ОБРАБОТКА СИГНАЛОВ С ЗАТУХАНИЕМ Пахотин В. А., Власова К. В., Кочмарский А. В.
Балтийский федеральный университет имени И. Канта, г. Калининград, Россия (236 040, г. Калининград, ул. Ал. Невского, 14-а) p_ksenia@mail. ru
В настоящей работе рассматривается обработка экспоненциально затухающих сигналов на основе теории оптимального приема. Показано, что при применении метода максимального правдоподобия можно анализировать сигналы, имеющие в своей основе две составляющие. Они могут незначительно отличаться по частоте, что делает невозможным их разрешение с помощью спектрального анализа. Сделан вывод о зависимости точности оценки параметров двух сигналов от отношения сигнал/шум. Получены формулы для обнаружения, определения параметров и оценки дисперсии параметров сигналов. Результаты математических расчетов подтверждены модельными исследованиями.
Ключевые слова: экспоненциально затухающие сигналы, метод максимального правдоподобия, дисперсия параметров сигналов.
OPTIMUM PROCESSING OF SIGNALS WITH ATTENUATION Pahotin V.A., Vlasova K.V., Kochmarsky A.V.
The Baltic federal university of a name of I. Kanta, Kaliningrad, Russia (236 040, Kaliningrad, street Al. Nevskogo, 14) p ksenia@mail. ru
In the present work processing exponents fading signals on the basis of the theory of optimum reception is considered. It is shown that at application of a method of the maximum credibility it is possible to analyze the signals having in the basis two components. They can slightly differ on frequency that does impossible their permission by means of the spectral analysis. The signal/noise is made you-waters about dependence of accuracy of an estimation of parameters of two signals on the relation. Formulas for detection, definitions of parameters and an estimation of a dispersion of parameters of signals are received. Results of mathematical calculations are confirmed by modeling researches.
Key words: exponentsfading signals, a method of the maximum credibility, a dispersion of parameters of signals.
Наряду с импульсными и непрерывными сигналами на практике часто встречаются сигналы с затуханием по экспоненциальному закону. Это, например, относится к сигналам свободной индукции, получаемых в спектрометрах ядерного магнитного резонанса (ЯМР), ядерного квадрупольного резонанса (ЯКР) /1, 2/. Для обработки такого рода сигналов, как правило, используется спектральный анализ /3/ и лишь в частных случаях метод линейного предсказания /3/, минимума дисперсии /3/, «MUSIC», /3/, которые характеризуются высоким разрешением. При использовании спектрального анализа спектр экспоненциально затухающего сигнала оказывается комплексным. Мнимая часть этого спектра определяет линию поглощения, а действительная часть спектра определяет сигнал рассеяния.
В настоящей работе рассмотрен вопрос обработки экспоненциально затухающего сигнала на основе положений теории оптимального приема /4, 5/. Она позволяет получить более точную оценку параметров сигнала по сравнению со спектральным анализом и определить возможность обнаружения слабого сигнала с оценкой отношения сигнал/шум.
Рассмотрим решение задачи обнаружения экспоненциально затухающего сигнала на основании положений теории оптимального приема. Пусть принятое сообщение имеет вид:
y (t) = U0e~at+iat + Um (t), (1)
где U0 — комплексная амплитуда сигнала- a — коэффициент затухания сигнала- ш — круговая частота сигнала-
Um (t) — аддитивный гаусовский шум со средним значением, равным нулю, дисперсией а2и интервалом корреляции хк.
Запишем функцию правдоподобия в двух случаях: при наличии сигнала в принятом сообщении!*^!), где Я- вектор параметров сигнала, и при отсутствии сигнала в принятом сообщении L2 (Я)
L1(X) = const exp (-2^/0°°|У (0 — Uoe~at+iat|2dt), (2)
L2(l) = const exp (-2^ f^ytofdt^,
гдеconst — константа, несущественная для процесса обработки сигнала. Она может быть определена из условия нормировки.
Если выполняется условие
Li (A) & gt- L2 (я), (3)
то в принятой реализации сигнал присутствует. Если условие (3) не выполняется, тогда сигнала в принятой реализации нет. Подставляя в (3) функции правдоподобия (2), можно получить неравенство
U* /" f (t)e~at+iMdt & gt- ^ Jo°° e~2atdt = h0, (4)
где U *- комплексно сопряженная амплитуда сигнала.
Правая часть этого выражения определяет половину энергии сигнала# = Щ~. Она
принимается в качестве порогового значения h0. Левая часть неравенства (4) определяет структуру оптимального приемника, на выходе которого формируется функция q. Если функция q будет больше порогового значения, то сигнал в принятом сообщении присутствует. Функция q является случайной и распределена по гауссовскому закону в связи с линейными преобразованиями при обработке (4). Математическое ожидание q равно нулю при отсутствии сигнала в принятом сообщении и равно энергии сигнала при наличии сигнала в принятом сообщении. Дисперсия случайной величины q одинакова для случая отсутствия сигнала в принятой реализации и при наличии сигнала в реализации. Она определяется выражением
D" = M (q2) = (5)
где M — оператор математического ожидания.
Среднее значение и дисперсия случайной величины q позволяет записать выражения для плотности распределения случайной величины q при наличии сигнала в принятой реализации P1(q'-) и при отсутствии сигнала в принятой реализации P2(q)
pi (q)=72kexp (-^) —
p2(q) = 7Щexp (-20& lt-6>-
Вероятность правильного обнаружения сигнала в принятом сообщении определяется выражением
Пр. обн. = р1(д)йц. (7)
Вероятность ошибочного решения определяется выражением
Рош. = .1 РгШЧ- (8)
Интегральные зависимости (7) и (8) можно представить в виде функции ошибок
Ф (ад
Ф№о) =. й, е-х2/1& lt-1х. (9)
В этом случае вероятность правильного обнаружения сигнала и вероятность ошибки могут быть представлены в виде
Рпр. обн. = Ф (- = Ф (- 1 Т $¦
Р& gt-ш. = Ф ЫУ = Ф (Т $) (10)
где N — количество некоррелированных отсчетов шума, содержащихся на интервале То =
гг 1
ле То = -
а
Таким образом, вероятность обнаружения экспоненциально затухающего сигнала определяется отношением сигнал/шум У/а до обработки и количеством некоррелированных отсчетов шума на интервале То.
Рассмотрим оценку параметров экспоненциально затухающего сигнала методом максимального правдоподобия. Для этого на основании принятого сообщения (1) запишем логарифм функции правдоподобия
у (ь) — (]ое-аЬешЬ йг, (11)
2о1тк-
где, А — оценочный вектор параметров сигнала.
Штрихами отмечены оценочные параметры сигнала. Дифференцируя (11) по 0о и приравнивая дифференциал к нулю, получим уравнение правдоподобия, из которого следует
до = 2а /" у (ь)е-(аЬе10}Ь йг. (12)
Дифференцируя (11) по параметру р = -а + ш и приравнивая дифференциал к нулю, получим второе уравнение правдоподобия, из которого следует
°о = (2а)2 у (1)е-аЬешЬСйС. (13)
Данные выражения определяют оптимальный алгоритм обработки экспоненциально затухающего сигнала. В ядерной спектроскопии для обработки сигналов используется преобразование Фурье. Однако теория оптимального приема рекомендует (12) проводить обработку такого рода сигналов методом преобразования Лапласа. Для этого выполним операцию математического ожидания в (12) от °о
м (ро) = 2а/" М (у (г))е-ссЬейг = 2а°о /" е-а+а^е1^-^ = ио^а-0О{--у
(14)
Спектр сигнала при преобразовании Лапласа является двумерным (а, ш) и комплексным. В точке максимума при, а = а и ш = ш математическое ожидание от оценочно-
го значения °о дает точное значение °о. Следовательно, решение выражения (12) несмещенное. Дисперсия оценочного значения °о может быть определена в точке максимума согласно выражению
Оо" = М°о — М (°о)2 = 2атка2. (15)
Следовательно, отношение сигнал/шум в в точке максимума спектральной линии при преобразовании Лапласа будет равно
в = 1? о]_ = Ш"!? 06)
^2аТк& lt-?2 о у 2
Оно совпадает с аргументом функции ошибок в (10).
Вторым уравнением правдоподобия является (13). Оно отличается от преобразования Лапласа весовым множителем 1 Смысл этого множителя заключается в том, что для оценки частоты и коэффициента затухания, первые отсчеты принятого сообщения менее информативны, чем отсчеты в конце сообщения. Данное выражение предпочтительнее при оценке частоты или коэффициента затухания.
Оценим дисперсии параметров экспоненциально затухающих сигналов. Для этого определим элементы информационной матрицы Фишера по выражению
/11 =. (18)
Ь = -М-7т±± Ь (17)
где М — оператор математического ожидания-
А — вектор параметров сигнала.
Элементы информационной матрицы Фишера находятся в точке максимума,
гда, А = А. Дифференцируя (11) два раза по амплитуде (] о, получим
1
2аТко'-2
Обратная величина определяет дисперсию амплитуды экспоненциально затухающего сигнала Ои
Ои = 2а%к ог = 2^-. (19)
Оно совпадает со значением дисперсии (15). Аналогичным образом можно получить дисперсию начальной фазы
т-
'-Ф° ~ I ОI 2 м
Дифференцируя (11) по частоте ш и коэффициенту затухания а, можно получить дисперсии О- и Оа
4т2
О- = Ба = ,|2 2. (21)
— а I 012МТд У '-
Таким образом, получены оценки дисперсии параметров экспоненциально затухающего сигнала, оценки согласно выражению Рао — Крамера.
Рассмотрим результаты модельных экспериментов для экспоненциально затухающего сигнала с двумя составляющими. Введем относительные параметры:
и0
а = - - относительная амплитуда-
7 -0
Ь =-----относительная частота-
а
х (ш) = - - относительная оценочная частота.
а
°Фо =. (20)
а
0,6 0,4 0,2 0,0 -0,2 -0,4 -0,6 -0,8 -I
& amp-

-20 -10 0 10 20 ь
Рис. 1. Экспоненциально затухающий сигнал с двумя составляющими
а «
0,6
0,3 0,0 -0,3 -0,6 -0,9
-20,0 -10,0 0,0 10,0 20,0 Ь
Рис. 2. Экспоненциально затухающий сигнал с двумя составляющими. Частоты сближены
1 а
0,8
0,6
0,4
0,2
0,0
-0,2
-0,4
-0,6
-0,8
-1,0 -I
1----1----г
-20 -10 0 10 20 Ь
Рис. 3. Экспоненциально затухающий сигнал с двумя составляющими
Сдвиг функций по параметру& quot- Ь& quot- Рис. 4. Корреляционная функция
Форма сигнала, состоящего из двух составляющих с параметрами: а1 =
1, а2 = 1,2, Ь = 0, Ь2 = 15, представлена на рис. 1. Коэффициент корреляции при таком разнесении составляющих близок к нулю И = 0,072. Шумовая составляющая равна 0,1. Следовательно, составляющие сигналов ортогональны и могут быть разрешены классическими методами. Оценки относительных частотЬ 1 иЬ 2 достаточно точно определяются по пересечению кривой с горизонтальной осью. Однако амплитуды оцениваются с достаточно большой погрешностью.
При сближении относительных частот Ь1 = 0, Ь2 =3 (рис. 2) становится сложно оценить относительные частоты Ь1 и Ь2 и даже структуру сигнала. Коэффициент корреляции двух составляющих равен в этом случае И = 0,26. Это значение коэффициента корреляции можно считать критическим.
На рис. 3 показан сигнал с составляющими, характеризующимися относительными частотами Ь1 = 0, Ь2 = 2. Коэффициент корреляции равен И = 0,46. При всем желании на этом рисунке выделить две составляющие невозможно.
На рис. 4 показана функция корреляции сигнала, созданного тонкой пленкой. Пунктиром разделены области, где возможно классическое решение (справа от пунктирной линии). Слева от пунктирной
линии классическое решение невозможно. В этом случае две составляющие сливаются вместе.
Рассмотрим возможности обработки одиночного сигнала, созданного одним слоем на основе метода максимального правдоподобия. Введем параметры обработки
а =
Уп
хп+Ь
(хп+Ь)2 + 1
Хп+Ь
(хп+Ь)2 +1
, Д (Ь) =
Щ Уп-^+Т-), Уп (хп+Ь) +1)
где, а — оценка относительной амплитуды- уп — выборка данных, N = 500-
Ь — оценка относительной частоты-
хп — выборка, определяющая зависимость от относительной текущей частоты, хп =-- Д (Ь) — нормированный функционал правдоподобия (чем больше значение функционала в максимуме, тем достовернее оценки параметрова и Ь) — у% - средний квадрат выборки данных.
Задавая произвольное значение Ь из области определения, можно определить а, а затем и точку функционала Д (Ь). Перебирая все значения Ь, получим функцию Д (Ь). Максимум этой функции определяет параметры, а и Ь. ЗначениеДтах связано с дисперсией шума и определяет степень достоверности оценок, а и Ь. На рис. 5 показана функция Д (Ь).
Она получена при следующих параметрах сигнала: а = 1, Ь = 2, иш = 0,01. Значение функционала в максимуме достигает Атах = 1025. В этом случае достоверность оценочных значений параметров высокая: а = 1,982, Ь = 2,05.
Помехоустойчивость метода обработки представлена в таблице 1. В зависимости от отношения сигнал/шум (первый столбец таблицы) даны оценочные значения параметров, а и Ь. В последнем столбце даны значения максимума функционала Атах. Из таблицы ясно, что удовлетворительные оценки можно получить при отношении сигнал/шум более 20 дБ. В этом случае Атах опускается до значения -11,5. Это ниеДтах может служить критерием для отбора достоверных оценок. При этом погрешность амплитуды составляет не более 6,8%, погрешность относительной частоты составляет не более 17,5%.
Таблица 1. Определение относительных значений параметров сигналов при различных отношениях сигнал/шум
Рис. 5. Зависимость функционала от относительной частоты
и'/и. дБ, а Ь. Дтах
60 0,99 998 2,0 1,1*105
40 1,0098 2,05 103
26 1,029 2,15 43,1
20 1,068 2,35 11,5
13 1,13 2,7 3,64
На основании данных последнего столбца (Дтах) можно оценивать отношение сигнал/шум, которое будет реализовываться в эксперименте.
Таким образом, в случае сигнала, созданного в спектрометре от одного тонкого слоя, можно использовать метод максимального правдоподобия.
Список литературы
1. Бородин П. М. (ред.). Ядерный магнитный резонанс. — Л.: Изд-во ЛГУ, 1982.
2. Марпл С. Л. -мл. Цифровой спектральный анализ и его приложения. — М.: Мир,
1990.
3. Перов А. И. Статистическая теория радиотехнических систем: учеб. пособие для ВУЗов. — М.: Радиотехника, 2003. — 400 с.
4. Тихонов В. И. Оптимальный прием сигналов. — М.: Радио и связь, 1983. — 320 с.
5. Пахотин В. А., Бессонов В. А., Молостова С. В., Власова К. В. Теоретические основы оптимальной обработки сигналов: курс лекций для радиотехнических специальностей.
— Калининград: Изд-во РГУ им. И. Канта, 2008. — 189 с.
Рецензенты:
Волхонская Е. В., д.т.н., профессор кафедры теоретических основ радиотехники Балтийской государственной академии рыбопромыслового флота, г. Калининград.
Захаров В. Е., д.ф. -м.н., профессор, зав. кафедрой радиофизики и информационной безопасности Балтийского федерального университета им. И. Канта, г. Калининград.

ПоказатьСвернуть
Заполнить форму текущей работой