CC BY
68
УДК 621.38:681.33:535.4
АНАЛИЗ И СИНТЕЗ ОПТИМАЛЬНОГО ОПТИЧЕСКОГО СЛЕДЯЩЕГО ФАЗОМЕТРА С ВРЕМЕННЫМ ДИСКРИМИНАТОРОМ
С.М.Слободян
Инновационный центр "ТЕСТ". г.Томск E-mail: IC_Test@ inbox.ru
Проведен анализ следящего оптического фазометра с временным дискриминатором видеосигнала, для которого получены оптимальные соотношения параметров при обработке видеосигнала малоконтрастного интерференционного поля, то есть, в условиях гауссова распределения и марковской модели сигнала и фона. Для пуассоновской модели распределения сигнала и фона приведены выражения, определяющие оптимальные размер и скорость следящего сканирования малоконтрастной интерференционной картины по условию минимума погрешности слежения за полосой интерференции.
Средства интерферометрии находят широкое применение в ряде оптических и радиотехнических устройств и систем [1-7]. Интерес к интерферометричес-ким средствам, обеспечивающим прием малоконтрастных низкой интенсивности изображений в условиях воздействия помех, особенно возрос в последние годы в связи с появлением персональных компьютеров высокой мощности и новых методов, например [2, 3, 5, 6], вейвлет-анализа сигналов и изображений, делающих возможным обработку картин интерференции со скоростью, достаточной для получения в реальном времени информации о динамике изменения трехмерной структуры фазового распределения изображений наблюдаемого средством пространства [2-4]. Задача получения оценок динамики изменения пространственно-временных характеристик волнового фронта оптического излучения значительно сложнее подобной [5, 6] радиоволновой задачи. Это объясняется тем, что фотоприемные (ФП) устройства оптического диапазона не осуществляют непосредственно регистрацию комплексной амплитуды волнового поля. Известные методы регистрации комплексной амплитуды, определяя способы получения пространственных и временных характеристик оптического поля, являются косвенными и основаны [1, 2, 4] на обширной гамме различных модификаций основного принципа интерференционного взаимодействия отдельных частей когерентного или некогерентного оптического поля. Название методов, например, голог-рафическая, двухэкспозиционная, спекл-, сдвиговая, муаровая, гетеродинная, интерферометрия интенсив-ностей отражают особенности аналитических подходов, дающих возможность регистрации параметров исследуемого явления, с точностью и разрешением, определенными размерностью длины волны используемого излучения. Так, интерферометрические устройства контроля формы поверхностей по измерению фазы фронта потенциально позволяют определять комплексную амплитуду оптического поля и другие информативные параметры как апериодических, так и колебательных сверхмалых ~(10-2.. ,10-3) А перемещений поверхности [4]. В этом случае изменения фазового фронта наблюдаются как результат сложения волны исследуемого поля с полем другого опорного когерентного оптического излучения. Регистрация этих изменений фазы фронта заключается в последующем изучении структуры, контраста, динамики по-
ведения и положения интерференционных полос (ИП) в плоскости изображения. Малый контраст изображения ИП при слабых световых сигналах является одним из основных факторов ограничения эффективности средств, учитываемый при проектировании систем адаптивной оптики и фазометрических систем контроля в реальном времени динамической деформации волнового фронта.
Эта проблема порождает необходимость удовлетворения противоречивых требований к фазо-метрической системе, а именно, обеспечения максимальной эффективности фильтрации помех при высоких динамической и статической точности в измерении и воспроизведения фазы волнового фронта анализируемого излучения. Наиболее эффективное решение этой проблемы возможно только на основе синтеза системы, реализующей минимум среднего квадрата случайной погрешности, обусловленной воздействием искажающих факторов - помех, при одновременном удовлетворении технических требований к качеству ее других характеристик, например, обеспечения наибольшей чувствительности. Вместе с тем, вопросы синтеза оптических фазометров с определенными характеристиками практически не разработаны и в литературе недостаточно освещены. Основные результаты, полученные в этом направлении и изложенные, например, в [1], не отражают всей полноты требований, предъявляемых на практике к показателям качества таких систем.
Задача синтеза оптимального оптического фазометра для контроля параметров волнового фронта излучения с математической точки зрения может быть рассмотрена как экстремальная задача. Она сводится к исследованию на условный минимум функционала, определяющего средний квадрат случайной погрешности (расстройки или невязки) системы при некоторых ограничивающих условиях, которые соответствуют требованиям чувствительности, статической и динамической точности оптического фазометра. Так, для рассмотренных в [8, 9] следящих фазометров размер и скорость следящего сканирования могут быть оптимизированы на основе общепринятой теории статистических решений по критерию минимума средней квадратической погрешности, максимизи-
рующеи точность системы слежения и зависящей от энергетических, например, отношения сигнал/шум и неэнергетических ее параметров. Рассмотрим их влияние подробнее путем исследования эффективности применения в фазометрах принципа временного дискриминирования совместно с интерференционным методом контроля фазы волнового поля. В фазометр этого типа заложен принцип слежения за подвижным объектом [10] с время-импульсным способом формирования сигнала ошибки [11], основанного на том, что смещение центров следящего сканирования и изображения ИП приводит к временному сдвигу видеоимпульса относительно опорного уровня - виртуального центра сканирования. Измеряя временной сдвиг видеоимпульса от ИП относительно сигнала сканирования для формирования определенного величиной этого сдвига аналогового или цифрового сигнала ошибки положения ИП, можно совместить центры сканирования и изображения ИП.
В представленном ниже варианте слежения за смещением ИП, в течение периода или интервала следящего сканирования Тск, ИП пересекает элемент ФП один раз, и выходной сигнал ФП устройства будет представлять собой результат свертки распределений плотности потока фотонов в изображении ИП размером 2г0 и функции чувствительности апертуры ФП, с учетом формы его элемента размером 2га. Для общности рассматриваем случай замкнутой области связной структуры фрактала изображения ИП некоторого порядка интерференции. При распределении энергетической освещенности в изображении пятна ИП, в том числе нулевого порядка кольцевой структуры, в плоскости ФП E(x',y) и функции совместного распределения прозрачности оптической апертуры и чувствительности элемента ФП ра(х,у) амплитуда сигнала на его выходе определится двумерным интегралом свертки
4 = еЕ0пг02а 2(1 + а 2)-1exp[ -(т? + п02)/(1 + я2)], (1)
где е - чувствительность ФП (в общем случае функция Я); й=га/г0, т0=х0/г0, п0=у0/г0 - относительные размер элемента ФП и координаты; Е0 - энергетическая освещенность сигнала (в точке х0=0; у0=0), с максимальным значением сигнала на выходе ФП
4 = еЕ0пг02а 2/(1 + а2). (2)
При равномерной прозрачности апертуры ФП круглой формы радиусом га
1, для IXI <7 Га2 - у2,1 у I <л/Га2 - X2;
0, для | х | >-у/г2 - у2, | у | >^1г2 - х2,
величина сигнала на основе (1) и (2) определяется выражением
„2 4п ч
Ра(х,у)=
и а =еЕ0ПГ0 х
Ра (X, у ) =
(здесь П01 = (у -у,)/п>; Ф(и) = [ехр(-г2) сИ - ин-
0
теграл вероятности) с максимальным значением сигнала
4м = еЕ0ПГ2[1 - ехр(-а2)]. (3)
Для квадратной формы апертуры ФП со стороной 2гке равномерной прозрачности
1, при IXI < к, I у I < к; 0, при | х | > к, | у | > к, сигнал на выходе ФП при а=гке/г0 равен
4 = 0,25еЕ0ПГ02 [Ф(а + т) + Ф(а-т0)]х х[Ф (а + П0) + Ф (а - П0)],
с максимальным значением 4м = еЕ0пГ02 [Ф(а)]2. (4)
Если учтем, что с точностью не хуже 0,7 % Ф(±и)«±{1-ехр[(-4/п)и2]}1/2, то получим максимальный сигнал, соответствующий наибольшей интенсивности
4м ~бЕ0ПГ02[1 - ехр[(-4/п) а2]}. (5)
Часто представляет интерес среднее значение сигнала, учитывающее случайность положения изображения ИП между тактами сканирования. Для одномерной процедуры сканирования с шагом 25х среднее значение сигнала на выходе ФП составит
1ср = | 4(Х0,у0)^(Х0)Су0,
(6)
а-п0
| [ехр(-П021)]
Ф(д/а2 - (П01 + П0)2 - т0) + +Ф(7а2 - (П01 + П0)2 + т» )
с1пл,
где ^(х0)=1/(25х) - равномерная плотность распределения. Для двумерного точечного с шагом 25хх25у следящего сканирования среднее значение сигнала ФП найдется из
5, 5Х
1ср = | | 4(Х0,у0)^(^(у0)Сх0Су0, (7)
-5У -5х
где ^(х0,у0)=1/(45х,5у) Для линейного изображения ИП: м(х0)=1/(25х), Цу0)=1; 5у=га. Для круглой и квадратной форм апертур ФП с гауссовской прозрачностью выражения (6) и (7) соответственно принимают вид:
4р = 4мФ(Ьх1)(4п /2Ъх\)ехр(-т0/\Д + а2),
4р.1 = 1смФ (Ъх1)п/(2Ъх1),
где Ъх1=5х/(г0лЛ+а2). При фоновой плотности потока ЕФ, величина которого косвенно учитывает темновую составляющую шума, сигнал ФП с апертурой круглой формы равномерной или гауссовской прозрачности равен: ¡Ф=еЕФпг2а, а на выходе ФП с апертурой квадратной формы равномерной прозрачности: 1Ф=еЕФ4г2ке.
На практике заданными являются яркость фона ВФ и сила света излучения или 10 - интенсивность интерференционной картины. При этом освещенность изображения полосы в плоскости ФП равна: Е0=т/0/(А/2г0)2 а фона Еф=птВф/(Д/2/ 2, где т -коэффициент пропускания среды и объектива; Д - диаметр входного зрачка и/ - фокусное расстояние объектива.
При следящем сканировании %(/) ИП (рис. 1 и 2) относительно опорного среднего положения Х}(?)=хио(1к+0,5ТС1) центра траектории сканирования, производимого пьезоприводом типа [12], выходной сигнал устройства представляет собой периодическую последовательность видеоимпульсов. Их огибающая промодулирована случайным образом по ширине и временному положению видеоимпульсов с внутренним заполнением случайным марковским потоком однофотонных импульсов, распределение интенсивности которых отражает изменение плотности потока энергии в сканируемом поле ИП. Для
упрощения формализации анализа точностных характеристик, примем, что для большинства ФП, используемых для регистрации интерференции, справедлива пуассоновская модель сигнала. Тогда, выходной сигнал ФП представляет собой пуассоновс-кую последовательность наносекундных импульсов одноэлектронной природы с плотностью распределения вероятности интервала их следования At р ^) = (1/ топ )exp(-A t/тст),
где тж=1/—(^ - средний период следования выходных импульсов, величина которого обратно пропор-
Рис. 1. Структурная схема формирования сигнала рассогласования ВД Огр.
Кл 1
Кл 2
СВ
СВ
ИНТ. ВЛ
ИНТ. ВЛ
->
1
1
1
......'
ИП
Поле сканирования
Ах
Рис. 2. Динамика переходных процессов при работе следящего фазометра с ВД
циональна их интенсивности пи(г), определяемой плотностью светового потока в изображении полосы на чувствительной поверхности ФП. При слежении за изображением ИП достаточно большой освещенности, сигнал ФП имеет вид последовательности видеоимпульсов. Форма их огибающей определяется видом распределений плотности светового потока в изображении ИП, прозрачности оптической апертуры и элемента ФП и их геометрией. Для соответствия принятой для анализа пуассоновской модели, будем считать, что полоса пропускания оптимального фильтра (рис. 1) блока обработки выходного сигнала ФП, в необходимой для реализации модели мере, согласована с частотой следования выходных импульсов. Среднее значение интенсивности следования импульсов видеосигнала одноэлект-ронной природы, как математическое ожидание пу-ассоновской случайно распределенной величины, равно -сумме средних интенс-ивностей темнового шума и^), внешнего (фона пФ(г) и сигнального пС(г) потоков импульсов: п(г)=п()+пФ(г)+пС(г), определяемых принятыми (аддитивной, мультипликативно й, аппликативно й) моделями изображения. Здесь пФ(г)=аЕФ3Фп/д и пС(г)=аЕ03Фп/д вычисляются с учетом особенностей сканирования, ранее рассмотренных [8, 10]; д = 1,610-19 Кл - заряд электрона;
- площадь элемента ФП.
Сигнал временного дискриминатора (ВД) на интервале тд вычисления величины рассогласования положения видеосигнала и центра следящего сканирования, отражающий расстройку системы
или ее разбаланс, представляется в виде [10]
_
ид= кд\ [ ^ (г) - ад + т)]п(г ул,
0
где К-д - коэффициент пропорциональности, влияющий на крутизну характеристики преобразования ВД; Яи+1(г,т) - нормированная знакопеременная функция, описывающая алгоритм работы ключевых каскадов ВД. Она принимает значение "1" при замкнутом состоянии ключа и "0" -при разомкнутом состоянии. В рассматриваемом случае, учитывая дискретные свойства пуассоновского потока импульсов, с достаточной степенью точности анализа операцию интегрирования можно заменить суммированием на том же временном интервале усреднения. Тогда среднее значение сигнала рассогласования положения видеосигнала относительно центра следящего сканирования, определяемое как математическое ожидание, и его дисперсия найдутся по формулам
а
т = Кд ¿[Я (о - ад+т)]п(() м,
а2и = Кд ¿[Я (г) + Я+1(1+т)]п (г) а г.
Здесь учтены свойства знакопеременной функции Я(0;1): Я(г)-Я+1(г)^0 и Я2(г)=Я(г). Когда скорость следящего сканирования относительно высока по сравнению с динамикой смещения отслеживаемой ИП, то величина Ат интервала рассогласования положения видеосигнала и центра траек-
тории следящего сканирования длиной £ск лежит в области малых значений, и выражение для сигнала дискриминатора, с использованием методики [11], принимает простой вид: ти=К-г0-тд-Дг/ З^-Т^. Откуда крутизна дискриминатора
вд=ди/дДт=ХдПо.тд/Зск- Ta для t-JT=1 составит: ^КмА*.
Дисперсия сигнала рассогласования для тд= ТЛ и длительности видеосигнала слежения, определяемого на уровне e-1, TBC=2r0a(1+a)TJ£CK, найдется из выражения
= кд -Тел[пфт + no(2roa(1 + a)/4К)]. (8)
Нормируя на крутизну ВД среднюю квадрати-ческую погрешность рассогласования при следящем сканировании (время усреднения равно интервалу следящего сканирования) aa=aJQi, найдем оптимальное значение скорости и размера следящего сканирования из условия минимума погрешности слежения дос/д&ск=0:
Зскор! = ro a(1 + a)n0 / ТлПф.
Здесь n0 и -Ф определяются с учетом выражений (2-5), отражающих влияние на интенсивность сигнала формы и вида распределений чувствительности элемента ФП и оптической апертуры. Учитывая это для соответствующих вариаций форм апертур ФП, получим оптимальный размер траектории следящего сканирования £ск
£ скор! = П0шах Г)(1 + a) X
x([a2/(1 + a2)];[1 - exp(-a)];[1 - exp( -4 a/ n)]}/ a пФ шах.
Поведение полученной из (8) и нормированной на значение при а = 1 средней квадратической погрешности для однократного следящего сканирования в зависимости от отношения сигнал/фон для гауссовых апертуры и элемента ФП и изменении соотношения 2r0/£CK opt, принимающего значения, кратные длине волны A-(4r0/xA) или в угловой мере - кратные n(2r0/£c=2r0/%n), представлено на рис. 3. Абсолютное значение средней квадратической временной ошибки слежения за ИП для треугольной формы сиг-ала сканирования составит: Дт„ oM=0,25a(1+a)(-c/n>)-1 TJATJ2 Та, где Ту - интервал времени усреднения сигнала ВД, равный [8] дли-1 »» »» тельности формирования одной ступеньки управляющего напряжения основного канала, а отношение T/2 Та равно числу усредняемых значений отсчетов сигнала ВД на такте управления пьезоприводом.
Численное значение верхней границы временной погрешности (наихудшее значение при оценке точности) метода определим для равных значений и-те-сивности сигнальных и фоновых импульсов (-с=пФ) и элемента пьезопривода [12], выполненного на основе пьезопакета ПП-4 0 30x18x5,5 мм3 с суммарной емкостью 0,5 мкФ ±5 %, резонансной частотой /0=40 кГц, допустимым напряжением ±300 В и чувствительностью 5 10-3 мкм/В. Для Ту = 10-2 с, ТсЛ = 2,5-10-5 с и а = 1 значение верхней границы временной погрешности равно:
supАтслош =к42 -10 5c.
В долях длительности такта интервала управления (~10—2 с) она составит или порядка 0,4 %. С увеличением отношения сигнал/фон повышается и точность слежения. Для отношения сигнал/фон ~10 вносимая шумовыми флуктуация-ми погрешность оценки снижается на порядок, и верхняя граница временной погрешности слежения не превышает 0,04 % длительности одного такта интервала управления пьезоприводом [12].
0 2 4 6 8 10
Максимальное отношение сигнал/фон
Рис. 3. Влияние отношения сигнал/фон и размера сканирования на погрешность оценки координат ПИ
Основной задачей интерферометрии является наиболее полное определение параметров поля интерференции, формируемого в процессе преобразования волнового фронта интерферометром измерительного средства. Интерференционным картинам малого контраста свойственна сложная структура с локализованными фрактальными областями, отражающими свойства объектов и фона в поле наблюдения. Динамика изменения поля наблюдения обусловливает пространственные и временные вариации локализованных фрактальных областей, определяющих информационный формализм видоизме-
СПИСОК ЛИТЕРАТУРЫ
1. Адаптация в информационных оптических системах / И.Н. Матвеев, А.Н. Сафронов, И.Н. Троицкий, Н.Д. Устинов. Под ред. Н.Д. Устинова. — М.: Радио и связь, 1984. — 344 с.
2. Евдокимов Н.В., Комолов В.П., Комолов П.В. Интерференция динамического хаоса гамильтоновых систем: эксперимент и возможности радиофизических приложений // Успехи физических наук. —2001. — Т. 171. — № 7. — С. 775—795.
3. Дремин И.М. Дальние корреляции частиц и вейвлеты // Успехи физических наук. —2000. —Т. 170. —№ 11. —С. 1235—1244.
4. Слободян С.М. Диссекторные гетеродинные системы // Зарубежная радиоэлектроника. — 1986. — № 6. — С. 62—72.
5. Кравченко В.Ф. Синтез антенн на основе нового класса атомарно-фрактальных функций // Доклады РАН. — 2002. — Т. 385. — № 1. —С. 754—761.
6. Werner D.H., Ganguly S. An overview of fractal antenna engineering research // IEEE Antennas and Propagation Magazine. — 2003. —V. 45. — № 1. —P. 38 —57.
7. Журович К.А., Кириллов В.П., Михайлов Ю.А., Склизков Г.В., Судаков О.А. Метод измерения пространственно-частотной характеристики устройства регистрации на основе ПЗС-мат-
нения общей структуры полного изображения. Это определяет необходимость разработки новых методов для большей полноты анализа поля интерференции, чтобы определять совокупно с пространственно-частотными характеристиками поля в целом [2, 4] его локальные особенности [3] с возможностью оценки динамики локальных и глобальных вариаций изображения поля интерференции.
Наиболее широко для анализа поля интерференции используется основанное на габоровском представлении аналитических свойств оптических сигналов Фурье-преобразование [13]. Для этого сигнал задается в области переменных часто неограниченной протяженности в виде ортогонального разложения по частотным базисным функциям. Частотные компоненты разложения поля изображения определяются путем интегрального преобразования полной реализации сигнала на бесконечном интервале. Значительное усреднение затрудняет выделение локальных особенностей общей структуры поля интерференции. При обработке малого контраста картины интерференции зашумлен-ные ее локальные вариации практически теряются.
В настоящее время для анализа, обработки и оценки параметров волновых полей идет интенсивная апробация нового метода [2, 3, 5, 6, 14] многомасштабного анализа сигналов - "вейвлет-анализа". Описанный в [8] принцип следящего фазометриро-вания с оптимальными размером и скоростью следящего сканирования обеспечивает получение как локальных фазовых изменений, так и их глобальное распределение по полю интерференции при существенно сниженном объеме данных, необходимых для получения достаточной информации о поле интерференции. Это позволит осуществлять вейвлет-преобразование изображения фазового поля в реальном времени при меньших вычислительных ресурсах средств обработки.
рицы // Приборы и техника эксперимента. — 2000. — № 5. — С. 79—83.
8. Слободян С.М. Принцип следящего оптического фазометрирова-ния с временным дискриминатором // Известия Томского политехнического университета. —2004. —Т. 307. — № 3. —С. 68—73.
9. Слободян С.М. Следящий оптический фазометр // Известия Томского политехнического университета. — 2003. — Т. 306. — № 6. —С. 102—106.
10. Пустынский И.Н., Слободян С.М. Диссекторные следящие системы. — М.: Радио и связь, 1984. —136 с.
11. Митяшев Б.Н. Определение временного положения импульсов при наличии помех. — М.: Советское радио, 1962. —199 с.
12. Слободян С.М. Многомернокоординатный привод микроуправления // Известия Томского политехнического университета. — 2003. —Т. 306. — № 5. — С. 92—95.
13. Takeda M., Ino H., Kabayashi S. Fourier transform method of fringe-pattern analysis for computer —based tomography and interferome-try // J. Opt. Soc. Amer. —1982. —V. 72. — № 1. — P. 156—160.
14. Mallat S.G. A theory for multiresolution signal decomposition: the wavelet representation // IEEE Trans. on Pattern Analysis and Machine Intel. —1989. —V. 11. —№ 7. —P. 674—693.
