Scientific journal
Modern problems of science and education
ISSN 2070-7428
"Перечень" ВАК
ИФ РИНЦ = 1,006

PROCESSING OF SIGNALS IN A FMR THE METHOD OF THE MAXIMUM CREDIBILITY.

Vlasova K.V. 1 Pakhotin V.A. 1 Molostova S.V. 1
1 The Baltic federal university of a name of I.Kanta
In the given work the theory of optimum reception in the appendix to signals of a ferromagnetic resonance is considered. The question on increase of the permission of two signals is solved with close resonant frequencies. It is shown that at consecutive application of the theory of optimum reception resolution will be defined not by criterion of Relej, and the relation a signal/noise and frequency of digitization of readout in sample of the data. By means of modeling calculations comparison resolution for a method of the maximum credibility and for a classical method (transformation of Fure) has been spent. The following is received: at factor of correlation of two signals more than 0,09, classical method processings are given by an error, and at factor of correlation more than 0,5 signals become indiscernible. The method of the maximum credibility authentically defines signal parameters even at factor of correlation of an order 0,9. Also the analysis of a noise stability of a method of the maximum credibility has been carried out. Result of research: at the relation a signal/noise ~6 дБ parameters of two signals are defined authentically. At fall of the relation a signal/noise to-3,5 дБ frequencies also are defined authentically.
resolution.
method of the maximum credibility
Ferromagnetic resonance
Сигналы ФМР имеют вид экспоненциально затухающих сигналов. Преобразование Фурье позволят получить комплексный спектр, действительная часть которого определяет линию Лоренца, а мнимая часть - линию поглощения. В ФМР спектроскопии для обнаружения сигнала, а также для оценки параметров используют линию Лоренца. Она описывается выражением:

,  (1.)

где  - вектор параметров сигнала;

a - амплитудное значение, нормированное на коэффициент поглощения ;

b - резонансная частота, нормированная на коэффициент поглощения ;

x - текущая нормированная частота .

Методом максимального правдоподобия можно получить выражения для оптимальной обработки сигнала (1.) /1/

 

(2.)

где  - выборка данных, содержащая сигнал (1.).

Штрихами отмечены оценочные параметры. Задавая  в пределах области определения, можно найти  и определить точку функционала . Перебирая все значения , можно найти зависимость . Минимум этой зависимости определяет ее решение, т.е. оценочные значения  и . Если в выборке данных будет два или три сигнала ФМР, то зависимость  будет содержать два или три минимума, по которым определяются параметры этих сигналов ФМР. Это классический способ оценки параметров. Разрешающая способность этого метода ограничена областью ортогональности двух сигналов ФМР. Если два сигнала ФМР по частоте являются неортогональными, тогда данный метод разрешить эти сигналы не может.

Для получения высокого разрешения двух сигналов ФМР необходимо использовать другой алгоритм для обработки сигналов /1/. Записывая два сигнала ФМР в виде

 ,(3)

можно получить систему уравнений. Их решение следующее:

, (4.)

,

,

где черта сверху означает суммирование по индексу n.

Задавая в (4.) значения  и  из области их определения, можно получить  и . Подставляя их в выражения для функционала  и вычисляя, можно определить точку поверхности функционала в двумерном пространстве , . Перебирая все значения  и  из области определения, можно получить полную поверхность функционала. Минимум поверхности функционала определяет решение, то есть оценочные параметры a`1,a`2 b`1,b`2  .

Функционал  инвариантен относительно перестановки  и . Следовательно, поверхность функционала будет содержать два зеркально симметричных минимума. На иллюстрациях максимумы смотрятся лучше, чем минимумы. В связи с этим, при модельных исследованиях для нахождения минимумов поверхности функционала будет использоваться обратный нормированный функционал:

.  (5.)

В этом случае значения функционала будут изменяться от 1 до ∞ в зависимости от дисперсии шума.

Рассмотрим с помощью модельных расчетов решение задачи разрешения двух сигналов ФМР методом максимального правдоподобия. В модельных расчетах зададим следующие значения параметров:   - переменная величина, приближающаяся к значению . Количество точек обработки N = 400, интервал между отсчетами - . Разрешающая способность метода обработки определяется выражением:

. (6.)

Коэффициент корреляции определяется выражением:

.  (7.)

При отношении сигнал/шум 20 дБ и коэффициенте корреляции двух сигналов ФМР R = 0,02, сигналы практически ортогональны. В этом случае классический метод обработки (2.) и метод максимального правдоподобия (4.) дают практически одинаковые решения. Если сближать относительные частоты, приближая b2 к b1, тогда коэффициент корреляции начинает увеличиваться, что будет указывать на степень неортогональности сигналов ФМР. При достижении определенной границы классический метод не будет решать задачу разрешения двух сигналов, а метод максимального правдоподобия будет давать информацию о параметрах сигнала.

На рис. 1 показана суперпозиция двух сигналов ФМР с относительными частотами  , амплитудные значения оставлены теми же самыми. В этом случае коэффициент корреляции , а разрешающая способность - . На рисунке два сигнала ЯМР сливаются вместе. Визуально трудно определить наличие второго сигнала. На рис. 2 показан результат обработки классическим методом. Отмечается единственный максимум функционала, который дает следующие оценки параметров: . Они не соответствуют модельным значениям. Задача разрешения двух сигналов ФМР оказывается не решенной.

  

 

Рис. 3.Поверхность функционала

На рис. 3 показана поверхность функционала, полученная при заданных параметрах методом максимального правдоподобия. Один из максимумов дает следующие оценки параметров двух сигналов ФМР:  . Задача разрешения двух сигналов ФМР оказывается решенной, а оценки параметров соответствуют модельным значениям.

При заданном интервале между отсчетами  сигналы ФМР по частоте можно сближать лишь до значения . Приведем данные модельных расчетов для этого случая. Коэффициент корреляции , разрешающая способность . Максимум поверхности функционала позволяет определить следующие параметры:  . При этом начинают возрастать шумовые максимумы, и дальнейшее сближение частот окажется невозможным. Таким образом, предел разрешения двух сигналов ФМР в методе максимального правдоподобия будет определяться не критерием Релея (как в классике), а отношением сигнал/шум. При отношении сигнал/шум ~20 дБ разрешающая способность может быть увеличена дополнительно, если увеличить интервал между отсчетами Δx. Однако в спектрометрах интервал Δx не может быть уменьшен произвольно. Он должен быть равен интервалу корреляции шума на выходе спектрометра ФМР.

Если в принятом сообщении содержится лишь один сигнал ФМР, тогда классический метод по максимуму функционала дает следующие оценки сигнала: , значение функционала в максимуме . Метод максимального правдоподобия в этом случае создает поверхность функционала с одним максимумом. Он дает следующие оценки параметров:  , значение функционала . Признаком наличия в выборке данных одного сигнала ФМР является отсутствие зеркального максимума и малая амплитудавторого сигнала  (на уровне шумовых максимумов). Функционалы правдоподобия в двух методах оценки параметров близки по своему значению. В таблице 1 представлены зависимости параметров сигналов ФМР от разности частот . Во втором столбце даны значения коэффициента корреляции двух сигналов ФМР. В третьем  столбце даны значения разрешающей способности А. Как видно из таблицы, метод максимального правдоподобия вполне работоспособен даже при коэффициенте корреляции . Относительные частоты b1 и b2  определяются с высокой точностью, а относительные амплитуды оцениваются с точностью . Классический метод уже при коэффициенте корреляции  оценивает амплитуды с точностью .

 

Таблица 1. Сравнение разрешающей способности

R

A

Метод максимального правдоподобия

Классический метод

a1

b1

a2

b2

a1

b1

a2

b2

15

0,0047

7,3

10,22

100,0

11,97

115,0

10,28

100,0

12,01

115,0

12

0,014

8,

10,22

100,0

12,26

112,0

10,39

100,0

12,41

112,0

10

0,026

10,5

10,23

100,0

11,77

110,0

10,53

100,0

12,05

110,0

8

0,046

13

10,22

100,0

11,97

108,0

10,78

100,0

12,44

108,0

6

0,088

17

10,19

100,0

12,31

106,0

11,28

100,0

13,22

106,0

4

0,18

25

10,19

100,0

12,16

104,0

12,5

100,0

14,09

104,0

3

0,29

33

10,22

100,0

11,99

103,0

14,13

100,5

15,23

102,5

2

0,49

50,5

10,15

100,0

12,15

102,0

17,84

101,5

-

-

1

0,79

100,5

10,37

100,0

11,81

101,0

20,83

100,5

-

-

Δx=0,25

1

0,79

100,5

9,98

100,0

12,31

101,0

20,88

100,5

-

-

0,5

0,86

133

9,95

100,0

12,32

100,75

21,45

100,5

-

-

0,75

0,93

200,5

9,98

100,0

12,37

100,5

21,89

100,25

-

-

При К = 0,49 классический метод не может решить задачу разрешения даже с большой погрешностью в оценке параметров. В последних трех строках таблицы дана информация о параметрах сигналов ФМР при обработке с уменьшенным до Δx=0,25 интервалом между отсчетами. Даже при разности частот между отсчетами  - результаты оценок параметров сигналов ФМР методом максимального правдоподобия удовлетворительные.

Следует отметить, что разработанная методика, основанная на информации о линии Лоренца, не позволяет получить оценки поглощения сигналов ФМР. Амплитуда и частота в линии Лоренца являются относительными ( ).  При изменении коэффициента поглощения  меняются как параметр a, так и параметр b сигнала ФМР. Это воспринимается как изменение амплитуды и частоты сигнала ФМР. Это недостоверная информация.

Для оценки коэффициента затухания сигнала ФМР мнимой части спектра экспоненциально затухающего сигнала недостаточно. Мнимая часть спектра определяет линию поглощения:

.  (8)

В аналитическом выражении для линии поглощения (8.) параметрами также являются относительная амплитуда a и относительная частота b. Извлечь из (8.) информацию о коэффициенте поглощения достаточно сложно. В ФМР спектроскопии информация о коэффициенте поглощения определяется или по ширине линии поглощения (8.), или по виду логарифма огибающей экспоненциально затухающего сигнала. Такая методика приемлема для одиночного сигнала ФМР. Однако она недостаточна для оценки поглощения при суперпозиции двух или более сигналов ФМР с разными коэффициентами поглощения. В данных модельных исследованиях приведены результаты, полученные при одинаковых коэффициентах затухания двух сигналов ФМР. Исследования возможностей обработки сигналов ФМР с различными коэффициентами затухания будут проведены в следующих работах.

Рассмотрим с помощью модельных расчетов помехоустойчивость метода максимального правдоподобия при обработке суперпозиции сигналов ФМР. При моделировании зададим следующие параметры двух сигналов ФМР:  . Интервал между отсчетами - -Δx = 0,5 . Амплитуда шума меняется, определяя зависимость решения от отношения сигнал/шум. В соответствии с таблицей 1, при различии в частотах , разрешающая способность метода   и классический метод обработки сигналов ФМР не может решить задачу разрешения.

На рис. 4 показан исходный сигнал (линия Лоренца), полученный при отношении сигнал/шум ~14 дБ.

 

Рис. 5.Поверхность функционала

Шумовая дорожка хорошо видна и существенно затрудняет обнаружение второго сигнала ФМР. На рис. 5 показана поверхность функционала правдоподобия. Уровень пьедестала и максимумов, вследствие понижения отношения сигнал/шум, уменьшился до значения ~3 для пьедестала и ~5,5 для максимума. Однако оценки параметров сигналов ФМР достаточно хорошие:  .

Таблица 2. Исследование помехоустойчивости

 (дБ)

Метод максимального правдоподобия

Классический метод

a1

b1

a2

b2

a1

b1

40

10,01

100,0

12,02

102,0

17,4

101,25

26

10,06

100,0

12,11

102,0

17,59

101,25

20

10,12

100,0

12,23

102,0

17,75

101,25

14

10,25

100,0

12,47

102,0

18,08

101,25

10

10,38

100,0

12,7

102,0

18,43

101,5

6

10,63

100,0

13,18

102,0

19,1

101,5

3

9,1

99,75

15,23

101,75

19,78

101,5

0

9,35

99,75

16,1

101,75

20,79

101,5

-3,5

9,77

99,75

17,56

101,75

22,48

101,5

В таблице 2 представлены данные расчетов, позволяющие оценить количественно эффективность обработки. Из таблицы следует, что метод максимального правдоподобия обеспечивает оценку параметров двух сигналов ФМР при отношении сигнал/шум до -3,5 дБ и коэффициентом корреляции между ними - .

Рецензенты:

  • Волхонская Е. В., д.т.н., профессор кафедры теоретических основ радиотехники Балтийской государственной академии рыбопромыслового флота, г. Калининград.
  • Захаров В. Е., д.ф.-м.н., профессор, зав. кафедрой радиофизики и информационной безопасности Балтийского федерального университета им. И. Канта,  г. Калининград