CC BY
106
УДК 681.586.69
В. А. Пахотин, С. В. Молостова, М. А. Никитин, А. С. Чугайнов
РАЗРЕШАЮЩАЯ СПОСОБНОСТЬ ОПТИЧЕСКИХ ПРИБОРОВ
Results of processing of diffraction images of dot sources are presented by a method of the maximum credibility. It is shown that the permission of two dot sources can be better, than the Releevsky permission. Results of modeling calculations are presented.
Ключевые слова: оптические приборы, дифракция, теория оптимального приема, разрешающая способность, критерий Рэлея.
Key words: optic instruments, diffraction, theory of optimum reception, resolution capability, Rayleigh criterion.
Разрешающая способность оптических приборов ограничена явлением дифракции света. В ее результате на входном зрачке оптической системы точечный источник света преобразуется в протяженный объект — дифракционное распределение света. Ширина главного максимума определяется известным выражением.
где X -длина волны, Я -радиус входного зрачка оптической системы [1].
При наличии двух точечных источников света, согласно критерию Рэлея, их главные дифракционные максимумы будут наблюдаться раздельно, если угловое различие точечных источников Дф > 0,61А. Этим условием определяется разрешающая способность микроскопов, телескопов, глаза человека, т. е. всех оптических приборов, световой поток которых ограничен входным зрачком.
В настоящей работе рассмотрена теория оптических приборов с точки зрения метода максимального правдоподобия [2]. При этом решаются две задачи. Первая связана с оценкой углового положения одного точечного источника света ф0 и его амплитуды по дифракционному распределению поля при известных значениях радиуса отверстия и длины волны. Определяется также дисперсия параметра ф0. Во второй задаче оценивается угловое различие между двумя точечными источниками света и их амплитуды по дифракционному распределению светового поля. Принято, что дифракционные распределения света от разных источников складываются аддитивно, без учета фазовых зависимостей.
26
Представлены результаты обработки дифракционных изображений точечных источников методом максимального правдоподобия. Показано, что разрешение двух точечных источников может быть лучше, чем рэлеевское. Представлены результаты модельных расчетов.
9 * 0,61—, R
(1)
Вестник Балтийского федерального университета им. И. Канта. 2012. Вып. 4. С. 26-33.
Рассмотрим явление дифракции точечного источника света на круглом отверстии радиуса Я . Длина волны светового потока равна X. В этом случае напряженность светового поля в зависимости от угла дифракции д определяется известным соотношением
ТТ( \ тт ^ (кК8Ш(Я ~ Я)) тт і \
и (д) = и0 ( • (-----------у) + иш (я) ,
( эш ( - Яо))
(2)
где ио — амплитуда напряженности поля, д0 = у0 — угловое положение точечного источника света, Зх (кЯ 8т (д — д0)) — функция Бесселя первого порядка, Пш (д) — шумовая добавка, которая характеризуется
нормальным распределением с дисперсией S2 и средним значением, равным нулю [1]. Запишем на основании выражения (2) логарифм функции правдоподобия [2]:
1п ( ( ,ио))=--4-1
2& Хи
и (в )-ио
2(кЯэт(в - в'0))
(кЯ эт (в - в0))
(3)
Учитывая, что он является выпуклой функцией, продифференцируем выражение (3) по амплитуде ио и приравняем нулю. В результате получим
І и (в)/(в, воув
и'0(в'°) = —т--------------------
І | / (в, во ) dв
(4)
где /(в, в0 ) =
2^ (кЯ (в - в0))
(кЯэт (в - в'0))
Математическое ожидание от решения выражения (4) будет:
м (и о)=и о
І / (0,00)/ (000 )0
-Т_______________________ ,
т
ІІ / (0,00 ) d0
(5)
Оценим дисперсию параметров. Для этого получим двумерную информационную матрицу Фишера по выражению [3].
( ^ (¡8 (Ь (а )))^|
¿у =-М У , , >П , (6)
где а т =(и', в0) — вектор параметров сигнала.
В результате двойного дифференцирования логарифма функции правдоподобия, согласно (6), получим информационную матрицу Фишера:
27
28
О 2 Т.
V «21
(7)
22
где
«п = 112 (в. о'0 ув; «12 = «21 = | / (в. 00 ) (0 У) <
т (
Л (.Г (в. 00)
Вычисляя диагональные элементы матрицы, обратной к информационной матрице Фишера, получим дисперсию амплитуды Ои и дисперсию углового положения точечного источника Од.
Точечный источник расположен под углом 20 °. Полуширина дифракционного пятна определяется первым нулевым значением и согласно выражению (1) равна 2,33 °, что отражено на рисунке 1.
Рис. 1. Зависимость дифракционного распределения света от точечного источника от угла дифракции
На рисунке 2 показана производная функции по углу дифракции. Модуль производной функции часто используется вместо самой функции для более точной оценки углового положения точечного источника по минимуму.
В качестве примера можно привести пеленгационную характеристику при пеленгации разностным методом [2]. Она часто применяется для слежения за пеленгом. Если ввести нормированный коэффициент а,
корреляции г = ■
12
то можно получить дисперсию амплитуды и
(«11 «22 )
дисперсию углового положения точечного источника света в вице
(1 - г 2)
; =-
(8)
1 «
29
Рис. 2. Производная функции по углу дифракции
По аналогии с сигналами в радиотехнике [2; 3] а11 и а22 имеют смысл обобщенной энергии. Следовательно, дисперсии Би и Бд зависят от отношения «сигнал/шум» и от коэффициента корреляции между функцией и ее производной.
Таким образом, метод максимального правдоподобия позволяет дать оценку углового положения точечного источника света. Эта оценка не зависит от размера дифракционного пятна. Она определяется с точностью, определяемой дисперсией углового положения источника света. Вместо дифракционного пятна (протяженного объекта) в результате данной обработки можно получить точечное изображение, угловое положение которого случайно и характеризуется дисперсией. Эффект дифракции, который можно отождествить с аппаратной функцией оптических приборов, может быть исключен. Контуры протяженных объектов могут быть уточнены до значений дисперсии.
Рассмотрим случай дифракции двух точечных источников света на круглом отверстии радиуса Я . В этом случае распределение света в зависимости от угла дифракции будет представлено в виде суммы двух дифракционных распределений:
и I 01 -,, ' ( 5Ш (0 - »¡)) +и 2-'і (йч (0 - 02)) . )
( 1 1 (БІ11 (-01 )) 2 ( БІП (-02)) ш ( ^ ( )
где q1, д2 — угловые положения точечных источников света, их,и2 — их амплитудные значения. Обозначим
( 2'1 ((1(0-0-)) ; ( 2'1 (ЖЯ11 (И-02)) (1п)
Л ( * 1 ) ( БІП (-01)) ; Л ’ 1 ) (БІП (-02)
и запишем логарифм функции правдоподобия на основании (9).
Дифференцируя и приравнивая дифференциалы к нулю, получим систему уравнений.
и ( 0)/1 ( 0,0--) - и-/12 ( 0,0-)+и2л ( 0,0-- )/2 ( 0,02) и( 0)/2( 0,02)- и1-лЩЩЩ+и2л12 ( 0,0-)'
Черта сверху означает интегрирование по параметру q. Решим эту систему:
Если угловое различие между дифракционными максимумами удовлетворяет соотношению
то коэффициент корреляции близок к нулю и решения уравнения (13) совпадают с классическими решениями:
Если угловое различие меньше, чем определено в условии (13), то, согласно классическим представлениям (критерий Рэлея), светящиеся точки сливаются вместе. Коэффициент корреляции в этом случае будет отличен от нуля. Однако решения (12) и в этом случае дают более точные оценки угловых положений двух светящихся точек. Решения находятся по значению поверхности функционала А в двумерном пространстве углов 0[ ид'2. Функционал А определяется выражением и зависит от амплитуд и угловых положений источников света. Перебирая все значения углов 0[ ид'2, получим полную поверхность функционала. Минимум этой поверхности определяет решения, т. е. оценочные значения амплитуд и[ и и2 и углов 0[ ид'2 двух точечных источников. Решения (12) несмещенные. Они учитывают взаимную корреляцию между дифракционными максимумами двух светящихся точек и дают более точные значения их угловых положений.
Оценим дисперсии углового положения двух светящихся точек и их амплитудных значений напряженности поля. Для этого дважды продифференцируем логарифм функции правдоподобия. Получим элементы информационной матрицы Фишера:
Определяя диагональные элементы матрицы, обратной к информационной матрице Фишера, получим дисперсии амплитуд света от первой и второй точек Би 1. Би2.
и< = и ( в) / ( в. в)/2 ( в. 02) - яи ( в) / ( в. 02) 1 X ( в .01)/2 (в .02)- я2
и,= и ( в)/2 ( в .02)/2 (вд) - яц( в) / ( в .в1) 1 /12 ( в .в1 Л ( в .в2)- я2
(12)
30
2
лв = |в2 - в/ | > 0.61 я,
(13)
и = и ( 0)/, ( вв) и, = и ( в)/2 ( в.в2) 1 /12 ( в .0') ’ 2 /2 ( в .02)
(14)
где «1ц = /12 ( в. 01); *112 = ^21 = /1 ( в. 01) /2 ( в. 02 )= Я1; = /2 ( в. 02).
, А
-11° 1к
а111 а122 - Ш12
а111 а122 - Ш12
Аналогично можно получить дисперсии угловых положений светящихся точек Бд1, Бд2. В этом случае информационная матрица Фишера будет
Определяя диагональные элементы обратной матрицы, получим
Таким образом, дисперсия угловых положений светящихся точек зависит от отношений «сигнал/шум» и взаимной корреляции дифракционных максимумов. Дисперсии будут минимальны, когда Ш2 = 0. Однако, если отношение «сигнал/шум» достаточно большое (более 20 дБ), при приемлемых дисперсиях коэффициент корреляции может достигать значения 0,9 и выше. Это сильно перекрывающиеся дифракционные максимумы.
Предварительные модельные расчеты проведены с целью иллюстрации принципиальной возможности получения более точного решения задачи разрешения двух точечных источников света по их дифракционным изображениям, чем это следует из работы [1]. Принято следующее. Амплитуды напряженности поля от точечных источников света равны 2 и 1,5. Угловое положение первого из них 20 °, второго — меняется в пределах от 26 до 21°. Отношение «сигнал/шум» равно
30 дБ, радиуса отверстия к длине волны Ш / Л = 15. Обработка дифракционного распределения поля проводилось согласно классическим представлениям. Это соответствует методу Фурье с рэлеевским ограничением на разрешение двух светящихся точек. Вместе с тем обработка дифракционного распределения поля проводилось и по методу максимального правдоподобия (табл.).
Как видно из таблицы, благодаря методу максимального правдоподобия можно решить данную задачу вплоть до угла ф2 = 21,5 ° и, в принципе, эта возможность зависит от отношения «сигнал/шум». Метод Фурье позволяет решить задачу практически лишь до угла ф2 = 24 ° вне зависимости от отношения «сигнал/шум».
где
31
Таблица 1
Сравнение работы ММП и метода Фурье при решении дифракционной задачи
Модельное, значение, (рг Метод максимального правдоподобия Метод Фурье
и, и 2 9 92 и и 2 91 92
26,0 1,97 1,47 20,0 26,0 2,0 1,51 20,0 25,9
25,0 1,97 1,48 20,0 25,0 1,99 1,51 20,0 25,1
24,0 1,97 1,5 20,0 24,0 1,9 1,41 19,9 24,1
23,0 1,97 1,52 20,0 23,0 1,93 — 20,8 —
22,0 1,96 1,52 20,0 22,0 2,67 — 20,8 —
21,5 2,22 1,28 19,9 21,7 2,98 — 20,6 —
21,0 3,03 0,4 19,7 21,9 — — — —
На рисунке 3 показана поверхность обратного функционала А = 1/Л, зависящая от значений углов П = в[ и ^2 = в'2. Максимум обратного функционала определяет решение. По одному максимуму оцениваются одновременно два значения угловых положений светящихся точек. Необходимости введения термина «разрешение (разрешающая способность)» в методе максимального правдоподобия нет. В методе же Фурье каждой светящейся точке соответствует свой максимум и возникает задача их разрешения.
Рис. 3. Вид обратного функционала
В настоящей работе представлено решение задач оценки параметров и разрешения подобных сигналов в области дифракции света на круглом отверстии (входном зрачке оптических приборов). Решения
получены методом максимального правдоподобия. Угловое положение точечного источника света в новом методе решения позволяет исключить дифракционное пятно. Точечный источник преобразуется в точечный со случайным распределением углового положения, которое оценивается дисперсией, а не угловым размером. В результате границы светящихся областей будут определяться с большей точностью.
Разрешающая способность при цифровой обработке дифракционных максимумов, полученных от двух светящихся точек методом максимального правдоподобия, может быть увеличена по сравнению с рэлеевской разрешающей способностью. Ограничение на разрешение точек в новом решении связано с отношением «сигнал/шум». Новое решение задачи углового разрешения двух светящихся точек имеет большое значение для цифровых оптических систем: телескопов, микроскопов.
Список литературы
1. Сивухин Д. В. Оптика. М., Наука, 1980.
2. Пахотин В. А., Бессонов В. А., Молостова С. В., Власова К. В. Теоретические основы оптимальной обработки сигналов : курс лекций для радиофизических специальностей. Калининград, 2008.
3. Перов А. И. Статистическая теория радиотехнических систем : учебное пособие для вузов. М., 2003.
Об авторах
Валерий Анатольевич Пахотин — д-р физ.-мат. наук, проф., Балтийский федеральный университет им. И. Канта, Калининград.
E-mail: [email protected]
Светлана Валерьевна Молостова — канд. физ.-мат. наук, доц., Балтийский федеральный университет им. И. Канта, Калининград.
E-mail: [email protected]
Михаил Анатольевич Никитин — д-р физ.-мат. наук, проф., Балтийский федеральный университет им. И. Канта, Калининград.
E-mail: [email protected]
Александр Сергеевич Чугайнов — асп., Балтийский федеральный университет им. И. Канта, Калининград.
E-mail: [email protected]
About authors
Valery Pakhotin — Dr, professor, I. Kant Baltic Federal University, Kaliningrad.
E-mail: [email protected]
Svetlana Molostova — PhD, associate professor, I. Kant Baltic Federal University, Kaliningrad.
E-mail: [email protected]
Mikhail Nikitin — Dr, professor, I. Kant Baltic Federal University, Kaliningrad.
E-mail: [email protected]
33
Alexander Chygainov— PhD student, I. Kant Baltic Federal University, Kaliningrad. E-mail: [email protected]
