Из таблицы 1 следует, что видность уменьшается до 38% оптимального значения, когда интенсивности различаются в 25 раз, т.е., например, когда обладающая высоким коэффициентом отражения измеряемая поверхность расположена напротив непокрытой опорной Физо-пластины. В случае, когда оценка распределения фазы выполнена средствами фазовых измерений, а интенсивности определены ПЗС-камерой, это эквивалентно деградации контраста на 1,38 бита.
Если камера в составе интерферометра 8-разрядная, что имеет место в большинстве систем, это может привести к неприемлемой потере точности измерений, в то время как применение 10- и 12-разрядных камер не приводит к потере точности. В любом случае этот эффект может быть компенсирован увеличением числа измерений с усреднением получаемых результатов, т.е.
, что в данном примере означает примерно семь измерений.Если эти меры не помогают, то необходимо выровнять интенсивности двух пучков, например, путем уменьшения коэффициента отражения поверхности с
коэффициентом отражения. Возможные способы достижения этого выравнивания обсуждаются далее для двух вариантов оптической схемы интерферометра Тваймана-Грина и для одного варианта интерферометра Физо. В схемах Тваймана-Грина предполагается, что светоделитель BS не обладает поглощением и что измеряемая T и опорная R поверхности имеют коэффициенты отражения соответственно . Все три конфигурации представлены на рисунке 3.![]() a) Оптическая схема 1
![]() b) Оптическая схема 2
![]() c) Оптическая схема 3
отражения измеряемой T и референтной (опорной)
R поверхностей
Коэффициенты ослабления, отражения и пропускания светоделителя описываются формулами:
a) Оптическая схема 1
,
b) Оптическая схема 2
,c) Оптическая схема 3
.В оптической схеме интерферометра Физо коэффициент отражения опорной поверхности не должен быть больше 4% во избежание влияния на интерферограмму многократных отражений.
В оптических схемах интерферометра Тваймана-Грина имеет место только двухлучевая интерференция, в то время как оптическая схема интерферометра Физо внутри полости может испытывать многолучевую интерференцию. Этим эффектом можно пренебречь при условии, что коэффициент отражения опорной поверхности (так называемого "пропускающего пробного стекла") не превышает 4%. При сверхточных измерениях и использовании формулы для анализа результатов фазовых измерений необходимо принимать во внимание истинную форму интерференционных полос.
Плоскопараллельная пластина ослабителя в оптических схемах (рисунок 3) должна быть изготовлена из плавленого кварца или стекла ВК7 с хорошим поглощающим покрытием. Из рисунка 3 следует, что ослабитель A является встроенной частью полости и влияние ее (двойное прохождение) на прошедший сквозь нее волновой фронт будет сказываться на результате измерений, порождая погрешности в дополнение к тем, что присущи измеряемой поверхности. Если возможно провести измерение с пустой "полостью", то легко определить влияние ослабителя на измеряемый волновой фронт и тем самым математически скорректировать получаемый результат.
3.2.6 Зависимость контраста от источника излучения
3.2.6.1 Общие сведения
В обобщенном виде интерферометр можно считать состоящим из:
- источника излучения;
- коллиматора;
- интерференционной части;
- фокусирующей оптической системы;
- приемника оптического излучения.
Свет, излучаемый источником, предназначен для формирования интерферограмм, а приемник - для сбора информации, т.е. их записи.
Контраст интерферограммы зависит от следующих параметров:
a) области спектра;
b) степени когерентности излучения;
c) поляризации света.
Источник излучения выбирается в соответствии с назначением прибора.
3.2.6.2 Область спектра
3.2.6.2.1 Излучатель и приемник
Пользователь обязан согласовать спектральный диапазон источника излучения со спектральной чувствительностью приемника, причем первый из них должен быть близок к максимуму второго.
Примечание - При визуальном наблюдении интерферограммы наиболее подходящими источниками считаются источник зеленого (ртутная лампа с фильтром) или желтого (натриевая спектральная лампа, которая обеспечивает близкую к нулю оптическую разность хода) света (оба соответствуют максимуму чувствительности глаза); для операторов с признаками дальтонизма не рекомендуется использовать лазерный диод с излучением красного цвета.
3.2.6.2.2 Факторы, влияющие на пропускание и отражение излучения
Свет проходит через коллиматор, интерференционную часть и фокусирующую оптическую систему, которые расположены между излучателем и приемником. При этом он отражается или проходит через оптические элементы схемы, что может привести к изменению исходного спектрального состава излучения.
Примечание - Например, позолоченное зеркало имеет слабый коэффициент отражения в сине-зеленой области спектра, отражательная способность увеличивается в желтой его части и особенно высока в красном и инфракрасном диапазонах длин волн. При покрытии серебром пик отражения приходится на желтый участок спектра, а при алюминиевом покрытии все видимое излучение отражается очень слабо.
Следовательно, выбору излучателя должно быть уделено особое внимание.
3.2.6.2.3 Требования контроля
Если исследуемый объект предназначен для использования на определенной длине волны, пользователь ограничен в выборе излучателя и приемника и должен использовать эту длину волны. В некоторых случаях возможны измерения на другой длине волны, но это должно быть отражено в протоколе измерений.
3.2.6.3 Временная когерентность
3.2.6.3.1 Прагматический подход
Пучок излучения следует рассматривать как последовательность (цуг) световых волн. Длительность цуга эквивалентна длине когерентности. Цуг распространяется от излучателя к приемнику. Достигая светоделителя, цуг делится на два цуга, один из которых проходит через опорное, а другой - через предметное плечи интерферометрической части. Пройдя каждый по своему оптическому пути, оба цуга достигают оптического элемента, используемого для реконструкции волнового цуга (в качестве такового иногда применяют светоделитель). Эти два цуга при равных оптических длинах путей обоих плечей полностью перекрываются; в случае различия длин один из цугов задерживается по отношению к другому.
Интерференционный контраст максимален при равных длинах плеч и убывает по мере увеличения их разности.
Примечание - Контраст при использовании в качестве излучателя гелий-неонового лазера уменьшается при разности оптических длин порядка нескольких сантиметров (длина когерентности обычно считается равной или близкой к длине лазера).
Контраст, наблюдаемый при использовании недорогого лазерного диода, весьма незначителен при разности хода в несколько сантиметров и периодически меняется. Это объясняется многолучевой эмиссией лазерного диода.
Контраст, наблюдаемый в свете натриевой спектральной лампы, осциллирует по мере изменения разности хода; период осцилляций соответствует разности хода 0,6 мм, что вызвано наличием натриевого дублета.
Контраст, наблюдаемый в свете ртутной спектральной лампы низкого давления, стремится к нулю, когда разность хода приближается к 1 см.
Когерентность излучения ртутной спектральной лампы высокого давления уменьшается достаточно быстро по мере возрастания температуры лампы и спадает почти мгновенно до уровня менее 1 мм. Тем не менее можно увеличить длину когерентности с помощью реостата в цепи источника ее питания, охлаждая лампу в процессе работы.
3.2.6.3.2 Влияние поверхности объекта
Пользователь должен ответить на вопрос: "Соответствует ли когерентность излучения источника шероховатости поверхности объекта?" При измерении дефектов поверхности этот вопрос, как правило, не возникает, поскольку характерный размер шероховатости (вершина - впадина) не превышает обычно нескольких микрометров. С другой стороны, важно согласовать форму волнового фронта с поверхностью объекта. Это означает, что пользователь должен адаптировать форму волнового фронта к форме объекта путем изменения разности хода в точках, меньших, чем длина когерентности источника.
Примечание - Например, наилучший способ измерения качества асферической поверхности заключается в том, чтобы придать волновому фронту радиус почти идеальной сферы, минимизирующей при своем падении на объект расстояния до его поверхности.
3.2.6.3.3 Влияние толщины объекта
При прохождении предметного пучка сквозь измеряемый объект имеют место два эффекта.
Во-первых, возникает задержка цуга при его прохождении по предметному плечу. Рассмотрим объект толщиной T с показателем преломления n. В предметное плечо вводится разность хода 2(n - 1)T. Она может быть компенсирована перемещением зеркала в опорном плече в направлении от светоделителя.
Во-вторых, объект следует рассматривать как толстую (в виде окна) плоскопараллельную пластину. Она является источником собственного изображения, формируемого зеркалом в опорном плече и перенесенного из его начального положения на расстояние (n - 1)T/n в направлении к светоделителю.
В итоге пользователь хочет удлинить опорное плечо, а образец сокращает длину предметного плеча. Таким образом, пользователь стремится расположить опорное зеркало в компромиссном положении, принимая во внимание когерентность излучения и видность интерференционных полос.
3.3.1 Общие сведения
Идеальным считается интерферометр, правильные и точные результаты измерений которого не зависят от влияний окружающей среды: вибраций, температуры и ее градиентов, силы тяжести и распределения усилий в креплениях, удерживающих измеряемый образец. Такой идеальный интерферометр можно представить как фильтр, полностью блокирующий воздействия окружающей среды. В то же время он должен давать всю информацию о форме поверхности с высочайшей точностью в пределах субнанометровой шкалы, причем, при необходимости, в более чем миллионе измеряемых точек. Кроме того, такой идеальный прибор не должен искажать результаты измерений присущими ему инструментальными погрешностями или шумом.
Практически это невыполнимо, однако у разработчика имеется ряд возможностей ослабить влияние на интерферометр внешних воздействий. Кроме того, стабилизируя воздействующие факторы в лаборатории, можно получать результаты измерений с хорошей сходимостью, но не обязательно правильные. Поэтому значимость высококачественной (в смысле стабилизации внешних условий) лаборатории часто переоценивается. Однако при этом не следует забывать и о необходимости высокой квалификации оператора.
3.3.2 Влияние вибраций
При испытаниях и контроле оптических поверхностей (как, например, с использованием интерферометра Физо на рисунке 14 в ИСО/ТО 14999-1:2005) две интерферирующих волны отражаются от последней поверхности V и опорного зеркала M. Эти две отражающих поверхности должны быть смонтированы на опоре, обеспечивающей их стабильность. При этом зазор между ними должен сохраняться неизменным во время фазовых измерений с временной несущей, как описано в 4.4.
Основанием для выдвижения подобных требований является тот факт, что в случае дополнительного, вызванного вибрациями опоры изменения "интерферометрической полости" уравнение (8) (см. 4.4) оказывается неточным, поскольку
где a - амплитуда вибраций,
Если даже шум имеет сложный спектр колебаний, во многих случаях достаточно свести рассмотрение к самой низшей механической собственной частоте
Задаваясь подходящим алгоритмом реконструкции фазы и частотой кадров камеры, возможно вычисление фазовых погрешностей при данной амплитуде a в предположении, что она достаточно мала (порядка
К сожалению, во многих случаях не удается эффективно подавить вибрации и устранить их влияние. Тогда полезно изучить степень влияния вибраций на результаты измерений и принять меры к минимизации этих воздействий путем частичной модернизации измерительной установки. Например, возможно применение электронной стабилизации полости интерферометра с использованием высокоскоростной камеры, мощного лазера и "двухступенчатого" фазового алгоритма, а также эффективна обработка данных с использованием алгоритмов фазовых измерений, нечувствительных к нежелательным и неизвестным фазовым сдвигам между несколькими кадрами камеры.
В данном случае предпочтительны алгоритмы двух групп:
- алгоритмы измерения фазы с пространственной несущей частотой (см. 4.5);
- алгоритмы, не требующие априори известных фазовых шагов.
Последний подход особенно важен в случаях, когда уже заведомо гарантированы подходящие начальные значения фазовых шагов, впоследствии подвергающиеся уточнению путем итерационной оптимизации, т.е. когда виброизоляция достаточно совершенна, но остаточные ее влияния слишком сильны и мешают проведению качественных измерений интерферометром.
Выбираемые при этом фазовые шаги являются в данном случае дополнительными неизвестными, которые подлежат дальнейшему определению вместе со значениями фазы. Только значения фазы являются функциями пространственных координат и не зависят от времени, в то время как фазовый шаг между двумя соседними кадрами камеры является единственной зависящей от времени величиной.
При использовании интерферометра, чувствительного к фазовым изменениям, обусловленным вибрациями, важно знать, как проконтролировать степень влияния вибраций на результаты измерений. Для лучшего понимания описываемой далее процедуры полезно констатировать, что вызываемые вибрациями погрешности измерений
описываются формулой + ... малые члены высших порядков, (7)где
- используемого алгоритма;
- амплитуды вибраций;
- отношения основной частоты
- времени экспозиции камеры.
Заключенный в квадратные скобки основной член в формуле погрешности, обусловленной вибрацией, содержит удвоенную измеряемую фазу, т.е.
. Следовательно, для определения значимости a) ввести в интерферограмму полосы путем наклона объекта;
b) компенсировать центральную часть интерферограммы и считать число "полос", остающихся после обработки результатов измерений фазового сдвига.
Если это число равно удвоенному числу ранее установленных полос в интерферограмме, то
3.3.3 Влияние силы тяжести и крепления образца
В большинстве практически важных случаев влияние этого эффекта успешно оценивается математически (например, [4]). Подобный анализ представляет собой довольно сложную процедуру, связанную с решением дифференциальных уравнений. Очевидно, что столь строгий математический анализ необходим лишь при получении сверхточных результатов измерений.
Несколько основных фактов и правил помогают избежать этих трудностей:
a) жесткость оптических материалов изменяется на порядок (см. таблицу 2). Оптическое стекло ВК7 занимает по этому свойству среднее геометрическое место по отношению к экстремумам;
b) из механики известно, что прогиб двухопорной балки (с закрепленными концами) минимизирован, если точки опоры расположены на расстоянии одной пятой длины балки от ее концов (так называемые точки минимального отклонения отделены друг от друга расстоянием 0,55 длины балки). Это правило оказывается полезным, если горизонтально расположенная эталонная поверхность закреплена таким оптимальным способом;
c) таблица 2 демонстрирует несколько численных значений прогиба однородно нагруженных (сила тяжести) круглых пластин, закрепленных по окружности (наружному контуру). В случае A сравнению подлежат разные материалы; в случае B для одного и того же материала варьируется толщина в соответствии с известным правилом D/d = 8; в случае C представлена толщина, необходимая для получения при указываемом диаметре прогиба ~= 33 нм.
Таблица 2
Существуют различные способы выполнения интерференционных контрольных измерений. Здесь приводится обзор методов и источников систематических и случайных погрешностей, равно как и выполнения прецизионных измерений. Методики анализа интерферограмм делятся на две группы:
a) визуальные и ручные методики, дающие качественную информацию о форме измеряемого волнового фронта (поверхности) и предоставляющие простейшие характеристики (например, величину PV);
b) методики фазовых измерений, обеспечивающие получение полной количественной информации о волновом фронте (поверхности).
4.2.1 Общие сведения
Интерферограмма формируется при взаимодействии измеряемого волнового фронта с опорным волновым фронтом. Форма и ориентация волнового фронта должна быть выбрана так, чтобы получать типовые, легко интерпретируемые интерферограммы (например, прямые полосы).
Сравнивая отклонения полученной интерферограммы от идеальной, можно обнаружить и классифицировать различные оптические аберрации волновых фронтов, равно как и типичные искажения как волновых фронтов, так и поверхностей. Кроме того, очень важно выделить информацию, касающуюся знака аберрации волнового фронта, т.е. либо это искажение является пиком или впадиной, либо измеряемая поверхность выпукла или вогнута относительно опорного волнового фронта. Эта информация о фазе не содержится в статической интерференционной картине или интерферограмме. Для получения ее необходимо найти и зарегистрировать местоположение нулевого порядка интерференционной картины, т.е. то место, где разность хода интерферирующих волновых фронтов равна нулю. Зная это место и внося определенным способом изменение оптической разности хода, можно предсказать знак отклонения волнового фронта (поверхности) путем наблюдения за направлением движения интерференционных полос:
- если разность хода уменьшается, полосы движутся в направлении от местоположения нулевого порядка (для наклонных волновых фронтов - к утолщению воздушного клина);
- если разность хода увеличивается, полосы движутся в направлении к местоположению нулевого порядка (для наклонных волновых фронтов - к сужению воздушного клина);
- если угол наклона между волновыми фронтами увеличивается, расстояние между полосами уменьшается, и они сдвигаются в направлении к полосе нулевого порядка.
Возможно, интерферометр должен быть настроен так, чтобы имелся наклон, тогда разность хода содержит воздушный клин, и можно получить практически прямые полосы.
4.2.2 Пример 1 - Интерферометр Физо
Вогнутость будет вынуждать полосы двигаться к сужению воздушного клина (рисунок 4).
![]() a) Вид сверху
![]() b) Вид сбоку
![]() c) Интерферограмма
1 - вогнутость; 2 - выпуклость; 3 - опорная поверхность;
4 - измеряемая поверхность; 5 - расстояние между полосами
при наилучшей аппроксимации интерференционной картины;
интерференционной картины
и выпуклости, соответственно
При вогнутой поверхности полосы у края интерференционной картины искривляются в сторону утолщения клина, а при выпуклой - к сужению воздушного клина (рисунок 5).
![]() 1 - вогнутая сфера; 2 - выпуклая сфера; 3 - утолщенная часть
клина; 4 - опорная пластина; 5 - измеряемая поверхность;
6 - суженная часть клина
поверхностей соответственно
На рисунке 5 полагается, что угол клина между двумя сравниваемыми волновыми фронтами больше, чем искривления волнового фронта, в результате чего не образуются замкнутые интерференционные полосы. На рисунках углы клина сильно преувеличены.
Если направление клина неизвестно, его можно определить в соответствии с таблицей 3.
Таблица 3
В общем случае, когда полосы не прямые и не круговые, нужно рассматривать и учитывать следующее обстоятельство. Во-первых, для случая прямых полос должно быть определено направление нулевого порядка. Во-вторых, для определения формы поверхности волнового фронта (вогнутая или выпуклая) следует уточнить направление кривизны полос относительно нулевого порядка. Это помогает обнаружению локальных искажений.
В случае сферических волновых фронтов или поверхностей знак оптической разности хода
![]() ![]() 1 - надавливание; 2 - нулевой порядок; 3 - неизменное
расстояние, соответствующее размеру полосы
Рисунок 6 - Движение полос
Во всех предыдущих рассуждениях "поверхность" может быть заменена "волновым фронтом".
4.2.3 Пример 2 - Интерферометр Тваймана-Грина
В таблице 4 отображены несколько методик определения знака погрешности.
Таблица 4
Изменение расстояния без изменения относительного наклона волновых фронтов друг относительно друга практически невозможно.
При испытаниях оптических систем допускается визуальный анализ интерферограмм в простых случаях и при слабых аберрациях, что позволяет грубо оценить источники и природу этих аберраций. Несколько типичных случаев продемонстрированы на рисунках 7 и 8.
![]() a) Простая расфокусировка
![]() b) Простой наклон
![]() c) Дефокусировка и наклон
Интерферограммы демонстрируют юстировку интерферометра
без заметных аберраций
![]() a) Сферическая аберрация
![]() b) Астигматизм
![]() c) Кома
Дополнительная иллюстрация зависимости интерференционных
картин от аберраций
при наличии небольших аберраций
Для измерения параметров волнового фронта по интерференционной картине необходимо количественно оценить ее отклонение от максимально приближенного к идеальному изображения, характерного для некоего "идеального" волнового фронта, подлежащего измерению. В зависимости от типа такого теоретического волнового фронта в интерферограмме можно обнаружить некоторые типы искажений, т.е. искривление волнового фронта, нерегулярность расположения полос и нерегулярность осевой симметрии (см. ИСО 10110-5:2007).
В идеальном случае аппроксимация должна быть выполнена методом наименьших квадратов.
Практически, по возможности, следует работать с прямыми полосами (обладающими минимальным изгибом), что облегчит ручную обработку интерферограмм и избавит от необходимости использования линейки или специального параллелограмма. Если это невозможно (например, при использовании защитных очков), то центральную полосу следует сравнивать не с прямой линией, а с ломаной, соединяющей два конца и центр центральной полосы.
Отклонение от идеальной картины, обозначаемое искажением или дефектом, обычно находится между соседними полосами идеальной интерферограммы (см. рисунок 9).
![]() S - расстояние между интерференционными полосами практически
идеальной интерферограммы;
полосы от идеальной картины
Количественно отклонение оценивается отношением
В том месте, где фронты пересекаются, минимальное расстояние между ними становится отрицательным; при вычислении расстояния между пиком и впадиной следует учитывать знак. Если один из фронтов отображает теоретически плоский фронт, то расстояние "пик-впадина" именуется "деформацией типа пик-впадина".
Искажение в единицах длины волны равно произведению
Искажение = (
Например, при оценке качества поверхности зеркала при нормальном падении света в интерферометрах Физо или Тваймана-Грина "масштабный коэффициент интерферограммы" равен 1/2: расстояние между полосами =
В таблице 5 приведены значения "масштабного коэффициента интерферограммы" для некоторых оптических схем интерферометров.
Таблица 5
Отклонение формы поверхности может быть описано "шагом полос интерференционной картины" - этот термин имеет отношение не к поперечному расстоянию между полосами, а к тому, что число полос, различаемых в интерференционной картине, соответствует числу полуволн, укладывающемуся в отклонение формы поверхности.
О качестве распространившегося волнового фронта сквозь оптический элемент судят по деформации плоского волнового фронта, дважды прошедшего сквозь оптический элемент. Эта деформация измеряется вдоль направления распространения, т.е. перпендикулярно плоскости волнового фронта. Единицей измерения служит шаг полос интерференционной картины при двукратном прохождении света, равный длине волны.
Среднеквадратичное отклонение оценки искажения волнового фронта не может быть определено визуально или простейшим ручным способом обработки данных.
Более полное описание методик анализа интерферограмм содержится во многих книгах (в частности, в [5]).
4.4.1 Общие сведения
Основные методики аналитической оценки интерференционных картин описаны в 4.4.2 - 4.4.4, где используются изменения во времени опорной фазы
измеряемого волнового фронта (без учета знака и многозначности фазового сдвига). Рассмотрение ограничено наиболее общим случаем двухлучевой интерферометрии.Интенсивность I(x, y) интерференционной картины, зарегистрированной на однородной прямоугольной сетке, описывается уравнением (ИСО/ТО 14999-1:2005, подраздел 3.1.1)
где I0(x, y) - почти равномерная интенсивность фона;
V(x, y) - видность (контраст) полос;
- измеряемое распределение фазы;Поскольку в этом интерференционном уравнении содержится три неизвестных I0, V,
.Предположим, что взаимная ориентация как измеряемого, так и опорного волновых фронтов такова, что
Временная модуляция опорной фазы осуществляется путем непрерывного или скачкообразного изменения разности хода в одном из плеч интерферометра с помощью различных устройств и способов:
- движением опорного зеркала, смонтированного на пьезоэлектрическом преобразователе (PZT);
- наклоном стеклянной пластины;
- смещением дифракционной решетки;
- вращением поляризаторов в зависимости от поляризационной чувствительности интерферометрической установки;
- акустооптическим модулятором.
Эти способы и устройства используются для введения разности фаз между пучками либо дискретно, либо непрерывно, в зависимости от условий измерительной задачи, соответственно вводимому частотному сдвигу.
В зависимости от способа изменения фазы применяют три метода:
- гетеродинной интерферометрии
;- интерферометрии с фазовой синхронизацией
;- интерферометрии фазовых шагов и фазового сдвига
.Примечание - Способ, описанный в 4.3, основан на пространственном изменении фазового сдвига
соответственно углу наклона между референтным (опорным) и измеряемым волновыми фронтами, благодаря чему формируются клинообразные прямые полосы с определенной пространственной частотой. Более сложные методы, основанные на пространственном изменении фазы, описаны в 4.5.Разность фаз
между интерферирующими пучками варьируется с помощью непрерывного фазового модулятора (например, вращающейся радиальной дифракционной решетки, акустооптического модулятора или вращающегося двухлучепреломляющего элемента). Поскольку непрерывно возрастающая разность фаз достигает удвоенной частоты (второй гармоники), то возможно применение двухчастотного лазера. Модулированная интенсивность интерференционной картины регистрируется фотоэлектрическим устройством. Для каждой координаты (x, y) интерференционной картины (т.е. в каждой ее точке) измеряется сигнал с частотой биений Один из вариантов полагает размещение неподвижного приемника излучения в точке с координатами (x0, y0), а второй приемник сделать сканирующим по всей интерференционной картине. В каждой точке (x, y) электронным фазометром со сканирующим приемником излучения на входе определяется разность фаз
. Путем последующего интегрирования Другая возможность использования второго приемника состоит в наблюдении интерференции пучков, не воздействующих на датчик. Фазовый сигнал этого приемника фиксируется и служит в качестве реперного в электронной фазометрии.
4.4.3 Интерферометрия с фазовой синхронизацией
Разность фаз интерферирующих пучков модулируется синусоидально введением дополнительного управляемого фазового смещения
(9)Простейшим способом считается применение опорного зеркала с пьезоэлектрическим приводом. Результирующая интенсивность
(10)Уравнение для I(t) в данной точке интерференционной картины (x, y) содержит члены типа
и . Эти члены могут быть разложены в ряды по бесселевым функциям в качестве коэффициентов (11)Частотный спектр детектированного сигнала U содержит компоненты: U0 (постоянную составляющую) и
(12) (13)Компонента сигнала
(14)Компонента
Основные особенности этой разновидности интерферометрии:
a) фазовое детектирование весьма чувствительно к фазовым сдвигам внутри всего интерферометра при выполнении сканирования; среднеквадратичное отклонение погрешности практически не превышает
b) область санирования должна быть определена предварительно (с использованием любого типа интерактивного программного обеспечения), чтобы избежать прерывания фазового сопровождения (процесса сканирования) в случае пересечения среза апертуры интерферометра.
4.4.4.1 Общие сведения
Поскольку уравнение (8) содержит три неизвестных (I0, V,
. Для более трех измерений и получения более высокой точности следует применять аппроксимацию методом наименьших квадратов.Разность фаз двух пучков должна быть дискретной, но предварительно фаза опорного волнового фронта сдвигается так, чтобы R измерений были равномерно распределены в пределах одного или более периодов модуляции:
(15)Для решения результирующей системы уравнений уравнение (8) следует переписать в виде
(16)где L = I0;
умножается на умножается на и r = 1, ..., R.Отсюда следует
(17)где Ir - интенсивность, относящаяся к опорной фазе
(18)Умножение на
В зависимости от необходимого числа фазовых шагов существуют несколько общих схем.
4.4.4.2 Четырехшаговая схема
(19)а видность полос описывается выражением
Значение V(x, y) представляет информацию о надежности измерения
в каждой точке или об отсутствии таковой. Если процедура измерений автоматизирована, то эта информация используется для определения границы апертуры интерферометра или наличия частиц пыли, интенсивность излучения в которых не зависит от (см. 4.4.4.3), то видность равна 0,9Vstep (20).Что касается времени измерений, то измерения фазы с применением шагового двигателя нельзя отнести к числу самых быстродействующих процедур, поскольку двигатель требует после каждой остановки периода стабилизации положения прикрепленного к нему зеркала. В модифицированной схеме фаза сдвигается непрерывно в пределах интервала времени интегрирования приемника.
Фазовое изменение (фазовый шаг)
. Для каждого из участков суммарная интенсивность равна (21)После интегрирования этого выражения зарегистрированная интенсивность равна
(22)где
.Таким образом, единственное отличие этого метода от ступенчатого состоит в уменьшении контраста полос. При рассмотрении различных фазовых ступенчатых схем коэффициент уменьшения контраста будет приниматься во внимание.
4.4.4.4 Трехшаговая схема
Существует несколько вариантов реализации трехшаговой схемы; в большинстве случаев применяются шаги 60°, 90° или 120°:
(23) (24) (26) (27)Контраст V уменьшается до уровня 0,9Vstep, если интегрирование фазы осуществляется в пределах
: и (29) (30)Контраст V уменьшается до уровня 0,83Vstep, если интегрирование фазы осуществляется в пределах
.При фазовом сдвиге большем, чем
(31)при фазовых сдвигах
При проведении точных измерений фазовый сдвиг известен (например, путем калибровки фазосдвигающего устройства). Используя четырехшаговую схему, возможно определить неизвестный, но постоянный фазовый сдвиг
(32)Фаза
(33)При измерениях по этой схеме модуляция интенсивности вычисляется по формуле
(34)Здесь предполагается
. Если фазовый сдвиг выходит за пределы +/- 10%, то и оценка V также превосходит +/- 10%.Эта схема эффективна при неизвестном фазовом сдвиге. Фаза
в каждой точке интерферограммы может быть определена по модулю 4.4.4.6 Калибровка фазосдвигающего устройства
Точность результатов, получаемых с применением описанных схем, зависит от линейности фазосдвигающего устройства. В частности, в случае применения пьезоэлектрического устройства (актуатора) (PZT) линейность во многих случаях зависит от линейности его характеристики преобразования, искажаемой гистерезисом и другими нелинейными эффектами. Следовательно, необходима калибровка фазосдвигающего устройства.
Одним из вариантов служит схема Карре (см. 4.4.4.5). Для вычисления
При наличии пяти результатов измерений интенсивности, разделенных равными фазовыми интервалами, может быть использовано следующее уравнение [8]:
(35)При некоторой априорной информации о
Измеренные значения позволяют установить скорректированные значения напряжений, подаваемых на контроллер фазосдвигающего устройства, управляемого компьютером.
4.5.1 Анализ полос с применением преобразования Фурье
Если пространственная частота
, то не наблюдается замкнутых полос и для получения из интерферограммы значения фазы Уравнение (8) может быть переписано в виде
(36)где
.После преобразования Фурье получаем
(37)где A, C и C* - комплексные фурье-амплитуды.
Путем цифровой фильтрации одна боковая полоса распределения интенсивности фильтруется и подвергается обратному быстрому преобразованию Фурье. В результате определяется величина c(x, y), по которой и вычисляется фаза:
(38)или с использованием комплексного логарифма
![]() Если применяется теорема дискретизации (т.е. на период косинусной функции приходится более двух элементов приемника излучения), неопределенности значений фазы легко устраняются (см. 4.4.4).
При частоте полос на интерферограмме, меньшей максимального наклона волнового фронта, т.е.
, применимо двухкоординатное быстрое преобразование Фурье [10]. Например, замкнутые полосы на интерферограмме очевидно свидетельствуют о нарушении этого условия.При использовании двухкоординатного алгоритма уравнение для F(I) выглядит следующим образом:
(39)Фильтрация комплексной величины C выполняется в двух ортогональных направлениях (с различными фильтрами в
Неоднозначности могут быть устранены при выполнении двух независимых вычислений
Второй вариант процедуры основан на оценке двух интерференционных картин, когда одна из них сдвинута на
по отношению к другой. Соответствующие коэффициенты Фурье c2(x, y) и c1(x, y) можно было использовать для вычисления (40)Считается, что
.4.5.2 Пространственно-синхронный анализ интерференционных полос
В тех местах интерферограммы, где полосы расположены чаще, локальный наклон волнового фронта становится более крутым, чем средний относительный наклон между измеряемым и опорным волновыми фронтами. В этих местах коэффициенты
и в уравнении для F(I) частично перекрываются, и знак становится неопределимым (изменение знака вычисленной фазы). Для преодоления этого затруднения следует увеличить относительный наклон опорного волнового фронта, соответствующий наивысшей пространственной частоте Необходимая информация о локальной фазе
измеряемого волнового фронта может быть извлечена путем умножения соответствующей уравнению (8) интерференционной картины непосредственно на косинусную и синусную функции, обладающие той же линейной зависимостью фазы, что и средняя несущая пространственная частота интерферограммы [11, 12].Поскольку интерференционная картина обычно сканируется, возможно использование однокоординатного преобразования. Пусть f0 является грубой оценкой средней пространственной частоты f интерференционной картины. Тогда
(41)и
(42)Если разность между измеряемым волновым фронтом
и опорным волновым фронтом мала, то третий член в этих выражениях отображает компоненту с низкой пространственной частотой, которая может быть отрезана фильтром нижних частот с функцией окна H(x, y) (предпочтительно типа Хэннинга).Результирующая муаровая картина в полосе низких пространственных частот записывается в виде:
(43)
(44)и
(45)Так как окно свертки может иметь ширину нескольких полос, требуется высокая их частота, чтобы при каждом измерении фазы оно было существенно меньше размеров зрачка. При этом теряется значительно меньше фазовой информации, чем в случае, когда окно свертки перекрывает край зрачка.
Простейший способ получения квадратично увеличенной муаровой картины состоит в формировании одноразмерного опорного наклона применительно к интерферограммам, к которым может быть добавлен значительный наклон (за исключением замкнутых полос). Тогда
(46) (47)и
![]() При другом варианте высокочастотная компонента удаляется из уравнения путем фильтрации на частоте f0. Тогда
(48)и
(49)что приводит к выражению для фазы
(50)Этот метод именуется "алгоритмом синусоидального окна".
В третьем методе используется комплексное синусоидальное окно, чтобы отфильтровать фазовую информацию одного бокового лепестка Фурье-преобразования
(51)в результате чего получается функция, вещественная и мнимая части которой идентичны F1 и F2 в методе синусоидального окна. Фаза в этом случае равна арктангенсу отношения реальной и мнимой частей.
Предыдущие уравнения справедливы в пределах модуля
a) определение фазы по модулю
в выражении для фазы, где имеет место арктангенс этих величин. Фаза вычисляется по модулю Таблица 6
В схеме Карре вместо
(52)и
(53)b) в случае выполнения условия дискретизации и когда разность фаз между соседними пикселями превышает
В интерферометрах, непосредственно измеряющих волновой фронт, его фазу или форму поверхности, регистрация информации осуществляется непосредственно ПЗС-камерой. В каждой точке с координатами xp, yp, соответствующей определенному пикселю на камере, где имеется модулированная фазой интенсивность интерферограммы, регистрируется фаза или "высота". Некоторые объекты заполняют всю ПЗС-матрицу, а другие - лишь ее часть. В первом случае центр матрицы, содержащей зарегистрированную информацию, считается центром или началом обрабатываемых данных. Отношение числа пикселей по вертикали к числу пикселей по горизонтали, отображающих круговые полосы, называется форматом изображения (коэффициентом формы).
Во втором случае началом координат обрабатываемых данных является центр тяжести (или центроид) базы данных. Для объектов с осевой симметрией должно быть показано, как координатная система измеряемого объекта соотносится с системой координат матрицы пикселей. У большинства ПЗС-камер пиксели имеют форму не квадрата, а прямоугольника. Однако всегда желательно, чтобы круглый объект выглядел и обрабатывался как круглый. Как центрирование, так и коэффициент формы пикселей представляются следующим образом. Суммируя x- и y-координаты каждого пикселя изображения, определяют средние значения множества всех xp и yp и для модулированных фазовых данных.
В результате получаются средние значения x0 и y0, являющиеся началом координат в системе координат пикселей. Координаты центрированных данных
xpc = xp - x0 и ypc = yp - y0. (54)
Именно эти центрированные данные в системе координат пикселей и подлежат статистической обработке различными методами.
Проблемным является вопрос выбора координатной системы при картировании криволинейных поверхностей на плоской поверхности ПЗС-матрицы. При плоской поверхности поперечные (x, y) координаты объекта картируются линейно на поверхности ПЗС-матрицы. Для поверхности, обладающей оптической силой (вогнутой или выгнутой), в частности обладающей высокой числовой апертурой, возникает вопрос о способе картирования поверхности на плоскость чувствительной поверхности приемника. В большинстве интерферометров измеряемый объект размещается в плоскости входного зрачка, точнее, в плоскости вершины объекта. Картирование осуществляется лучами, параллельными оптической оси, т.е. коллимированным пучком. Таким образом, координаты x, y измеряемой поверхности один к одному картируются в координатах плоскости зрачка, при этом все прогибы в направлении оптической оси подавляются. Плоскость зрачка картируется линейно в плоскость детектора так же, как и в случае плоской поверхности.
Отсюда следуют два вывода:
- плоская сетка из прямых, перпендикулярных оптической оси, расположенная между объектом и интерферометром, будет изображена в плоскости детектора с подушкообразной дисторсией;
- x, y-координаты в плоскости зрачка являются теми же x, y-координатами, что используются при анализе всех данных о фазе, т.е. поставляют информацию обо всех высотных или топографических отступлениях от плоской поверхности.
4.8.1 Представление и отображение фазовой информации
После преобразования или картирования фазовой информации из системы координат ПЗС-матрицы она может быть представлена двумя различными способами: либо в виде массива отдельных точек, либо после интерполяции в виде аналитических функций с соответствующим числом или набором коэффициентов. Каждый из способов в зависимости от дальнейшего использования полученной информации обладает присущими ему достоинствами и недостатками. Ни один из двух способов не является оптимальным во всех случаях. В ряде случаев часть интерферограммы в области низких пространственных частот лучше отображать аналитическими функциями, а высокочастотную - необработанной совокупностью пикселей.
Необработанные данные о фазе в виде пикселей привлекательны своей "необработанностью", поскольку содержат максимум присущей объекту информации. При любой "выборочной" информации ее объем слишком мал при использовании шкалы пространственных частот меньше двух пикселей (а иногда эта информация в пределах такой шкалы и вовсе отсутствует). Однако в случае дискретизации с использованием теоремы Найквиста получается достаточно полная информация.
Однако способ сохранения информации попиксельно имеет много недостатков, поскольку отсутствует возможность анализа ряда полезных аспектов информации, содержащейся в интерферограмме. При этом для получения ряда требуемых статистических параметров необходимо привлекать дополнительную стороннюю информацию.
Если информация в виде пикселов аппроксимирована 2D-аналитическими функциями, то много полезных данных дает изучение коэффициентов этих функций, при этом нет необходимости в привлечении сторонней дополнительной информации. Более того, при использовании определенных типов аналитических функций по их коэффициентам можно судить о качестве работы системы формирования изображений.
Недостатком аппроксимации является фильтрация части полезной информации, содержащейся, как правило, в области высоких пространственных частот. Если необработанные данные хотя бы частично исчезли (были удалены) из памяти, то аналитическим способом эту информацию восстановить нельзя.
Аналитические функции, описывающие поверхности круглой или прямоугольной формы, приведены в 4.8.2 - 4.8.6, а детальное их описание - в приложении A.
Поскольку большинство оптических систем имеют круглые зрачки, в качестве аппроксимирующих функций, отображающих оптическую поверхность и волновой фронт, обычно используются полиномы Цернике. Для единичного круга эта система полиномов является полной и ортогональной. Кроме того, при определенной нормировке коэффициенты полиномов Цернике являются средними квадратическими отклонениями индивидуальных оптических аберраций, присущих полученным фазовым данным, непосредственно относящимся к формированию изображений оптическими системами.
Полезные свойства полиномов Цернике распространяются только на круглые апертуры, а неправильные результаты получаются в случае применения этих полиномов к некруглым объектам. В случае отчетливой эллиптической апертуры полиномы Цернике применимы, поскольку на зеркало под углом 45° падает круглый пучок. В этом случае допустимая норма астигматизма для эллиптической апертуры эквивалентна режиму падения нормального или преломленного пучка.
В приложении A представлены первые шесть порядков полиномов Цернике, нормированных в двух вариантах в зависимости от применения, при этом использованы заимствованные из [13] обозначения с надстрочными и подстрочными индексами. Это отличается от различных подходов "серийной нумерации", используемых для облегчения компьютерного кодирования.
Одним из достоинств полиномов Цернике является возможность описания ими симметричных аберраций. Эти симметрии становятся очевидными, судя по обозначениям полиномов и их коэффициентам, причем подстрочные индексы указывали на радиальный порядок, а подстрочный индекс - на угловой или азимутальный порядок.
Свойства симметрии нескольких первых порядков очевидны. Следующим преимуществом полинома Цернике является возможность заключения по его порядку о дополнительных свойствах объекта, не различимых в традиционной схеме.
В приложении A описаны два способа нормирования: экстремальный, т.е. типа "вершина-впадина", обозначаемый прописной буквой, и способ среднеквадратичного отклонения - "среднеквадратичный", обозначаемый строчной буквой. В первой публикации, описывавшей полиномы Цернике, для сокращения объектов информации об интерферограмме был применен способ "вершина-впадина", причем коэффициент 1 свидетельствовал о двойной погрешности определения этого диапазона. Этот способ также применим и при интерпретации интерферограмм пороговым методом. Однако "среднеквадратичный" метод более полезен при необходимости определения одной или более аберраций оптической системы и необходим для точной расшифровки данных, получаемых способом "вершина-впадина".
В целях лучшего понимания настоящего стандарта было принято следующее решение, относящееся к нормированию и обозначениям.
Если коэффициенты Цернике обозначаются прописной (заглавной) буквой, как, например,
Преобразование одного способа нормирования в другой осуществляется умножением "прописных" коэффициентов на коэффициент нормирования (см. приложение A) для получения "среднеквадратичных" коэффициентов.
Полезной первоначальной операцией, предшествующей выполнению преобразования Цернике, является удаление из интерференционных данных вершин, наклонов, перекрытий. Среднеквадратическая погрешность представляет собой квадратный корень из суммы квадратов коэффициентов.
Подобно этому, коэффициенты высоких порядков могут быть сведены к нулю путем удаления аберраций высоких порядков. Среднеквадратическая погрешность остающейся части погрешности определения волнового фронта или поверхности является по-прежнему квадратным корнем из суммы квадратов коэффициентов.
Если интерферограмма представлена пикселями, то среднеквадратичная погрешность определения поверхности или волнового фронта равна
, где N - число пикселей в пределах площади открытой апертуры, а z - высота (пик) над или под средней (опорной) плоскостью. Суммирование выполняется по всем N точкам.Если аппроксимация полиномами Цернике используется при "среднеквадратичном" нормировании, то среднеквадратичная погрешность определения поверхности или волнового фронта равна
, где суммирование коэффициентов Цернике выполняется после того, как устранены наклон, фокус и т.д.4.8.3 Аппроксимация полиномами Цернике при эллиптической форме границы интерферограммы
С эллиптическими апертурами в оптических системах приходится часто иметь дело, поскольку на пути коллимированного пучка устанавливается плоское поворотное зеркало, изгибающее оптический путь. Во многих оптических системах часть мощности оптического пучка не фокусируется в плоскости изображений. Однако если плоское зеркало обладает хотя бы слабой преломляющей силой, да еще не при нормальном падении излучения, оно вносит астигматизм, искажающий характеристики оптической системы.
С другой стороны, если плоское зеркало обладает астигматизмом, совпадающим с направлением падения излучения, то этот астигматизм вносит в оптическую систему эффект преломления и дополнительной фокусировки. Таким образом, основной задачей применения полиномов Цернике при анализе эллиптических апертур является определение степени и направленности рефракции и астигматизма, присущих поверхности зеркала.
Процедура при этом такова. Сначала положим, что большая ось a апертуры направлена вдоль оси x. Пусть a и b - длины полуосей эллипса. Теперь умножим все x-координаты точек локализации полос на b/a. При этом зрачок становится круглым, а аппроксимация полиномами Цернике - полезной.
Если оптическая система осуществляет фокусировку преломленных лучей, то коэффициент
После устранения членов
В случае если аппроксимация полиномами Цернике была выполнена применительно к данным круглой апертуры, контурная карта поверхности поверх эллиптической апертуры должна быть получена умножением координат x, полученных при круговой апертуре, на отношение (a/b) с последующим построением полученной картины.
Возможны и другие случаи применения аппроксимации полиномами Цернике к фазовым данным, полученным с эллиптической апертурой. Описанный выше подход также может быть использован при условии отчетливого понимания того, как масштабирование координат в том или ином направлении повлияет на результаты аппроксимации.
4.8.4 Полиномы Цернике-Татляна для круглой апертуры с центральным круглым отверстием
В случае, когда имеется интерферограмма круглой апертуры (например, при контроле телескопа Кассегрена), аппроксимация может быть осуществлена с помощью модифицированных полиномов Цернике (см. приложение A). Эти полиномы являются функцией
4.8.5 Полиномы Лежандра для апертуры прямоугольной формы
В некоторых случаях для аппроксимации апертуры прямоугольной формы применяется 2D-аналитическая функция, которая обычно используется при испытаниях цилиндрической оптики. В этих ситуациях предпочтительно использование полиномов Лежандра.
Часто возможно вычисление ортогональных полиномов для определения требуемых характеристик поверхностей и участков произвольной формы, причем при необходимости с применением весовой функции. Математически это не новшество, но в оптике это не получило распространения. Термин ортогональная функция или "неполином" используется в настоящем стандарте как способ описания абсолютно любой базовой функции. Ею может быть, например, функция
для точного моделирования дефокусировки.В приложении B приведены несколько примеров ортогональных полиномов без указания деталей вычисления. С использованием матриц эти вычисления могут быть очень простыми, обеспечивая возможность получения совокупности полиномов один за другим детерминированным путем (без итераций и аппроксимаций).
Результаты каждых испытаний, калибровки или их серий, выполненных в лаборатории, должны быть оформлены аккуратно, ясно и объективно в соответствии с инструкциями по эксплуатации и испытаниям или методиками калибровки. Результаты обычно заносятся в протокол испытаний или калибровочный сертификат и должны содержать всю требуемую клиентом и необходимую для интерпретации информацию об испытаниях и калибровке.
Каждый протокол испытаний и калибровочный сертификат должен содержать, по крайней мере, следующую информацию до тех пор, пока лаборатория не сочтет ее излишней:
a) наименование (например "Протокол испытаний" или "Калибровочный сертификат");
b) наименование и адрес лаборатории, а также, при необходимости, местоположение пункта проведения испытаний и калибровки, если они проводились вне лаборатории;
c) отличительный идентификационный символ протокола или сертификата (например, серийный номер на каждой странице);
d) фамилию, имя, отчество и адрес клиента;
e) краткое описание метода измерений;
f) краткое описание, включая однозначную идентификацию испытанных или калиброванных образцов;
g) дату проведения испытания или калибровки;
h) результаты калибровки или испытаний с указанием единиц измерений;
i) фамилию, имя, отчество, должность и подпись оператора, подтверждающего правильность полученных результатов.
Твердые копии (CD или дискеты) протоколов испытаний и калибровочных сертификатов также должны включать в себя номер страницы и их общее число (например: стр. 1 из трех).
Рекомендуется, чтобы лаборатории включали требования о том, что протокол или сертификат могут быть воспроизведены только полностью.
В дополнение к изложенному в 5.2 протокол испытаний должен при необходимости дополнительно содержать следующие данные:
a) отклонения, дополнения или исключения (корректировки) метода испытаний, включая сведения об условиях окружающей среды;
b) при необходимости заключение о согласованности (или несогласованности) испытуемого объекта с требованиями или спецификациями;
c) при необходимости заключение, содержащее оценку неопределенности измерений;
d) при необходимости мнения и рекомендации, основанные на результатах измерений;
e) дополнительная информация, зависящая от специфики методики испытаний или требований клиента.
5.4.1 Основы
В дополнение к требованиям 5.2 при необходимости калибровочные сертификаты должны содержать следующие данные:
a) условия (например, окружающая среда), которые могли при калибровке повлиять на результаты измерений;
b) неопределенность измерений и/или утверждение о согласованности с идентичным метрологическим сертификатом, выданным ранее;
c) свидетельство о прослеживаемости результатов измерений к национальным эталонам.
5.4.2 Спецификация
Калибровочный сертификат должен относиться только к величинам и результатам, определяющим функционирование объекта.
Если утверждение о согласовании спецификации сделано без учета результатов измерений и присущих им неопределенностей, лаборатория должна зарегистрировать эти результаты и сохранить их для возможного обсуждения в будущем.
Когда установлено согласование с предыдущими результатами, то неопределенность измерений должна быть принята во внимание.
5.4.3 Юстировка или ремонт
Когда прибор для калибровки был отъюстирован или отремонтирован, то должны быть зарегистрированы результаты калибровки до и после выполнения этих операций.
Если в документ включаются мнения и рекомендации, то он должен содержать их обоснование. Кроме того, в протоколе испытаний их следует выделить.
Мнения и рекомендации, включенные в протокол, могут содержать следующее (но не обязательно ограничиваться этим):
- мнение о заключении, подтверждающем соответствие или несоответствие полученных результатов предъявляемым требованиям;
- рекомендации по использованию полученных результатов;
- рекомендации и указания по усовершенствованию объектов и методик.
В случае передачи результатов испытаний или калибровки электронными средствами (по телефону, факсом, факсимильно и др.) должны быть соблюдены требования настоящего стандарта.
Формат должен быть согласован с полученными данными во избежание непонимания и недоразумений.
Внимание должно быть уделено размещению материала в протоколе или сертификате с целью максимальной доступности восприятия его содержания читателем.
Заголовки должны быть предельно унифицированы для обоих документов.
Модернизация (изменения) протоколов или сертификатов должна быть выполнена в той же форме, что и исходный документ с указанием, что это "Дополнение к протоколу испытаний (калибровочному сертификату) с серийным номером (XXXX)x".
При необходимости выпуска полностью нового протокола или сертификата должно быть однозначно указано, какой документ он заменяет.
Информация, содержащая сведения о формате данных интерферометрического вычислительного процесса, может быть найдена, например, в /template/go.php?url=https://support.zygo.com/resource/manuals.cgi?type=mp или в /template/go.php?url=https://www.fisba.ch/(Metrology-Products-Software).
(справочное)
A.1 Общие сведения
В данной части описаны различные совокупности полиномов, используемые при испытаниях оптических элементов двух наиболее распространенных форм (таблица A.1):
a) круглой (преобразуемой в эллиптическую);
b) квадратной (преобразуемой в прямоугольную).
Каждая совокупность обозначается комбинацией букв:
- форма площади: C (окружность/эллипс); S - (квадрат/прямоугольник);
- основополагающий базис: P - (полярный); R - (прямоугольный);
- форма выражения: Z - (Цернике); P - (полярная); R - (прямоугольная).
Примечания
1 В этом приложении термины "полиномиальный" и "функциональный" иногда взаимозаменяемы. Например, функция Цернике может быть получена путем умножения полинома Цернике на гармоническую функцию (
2 Ортогональные совокупности, описываемые в приложении, являются ортогональными на непрерывных поверхностях определенной формы (круглой или квадратной). При использовании этих совокупностей в дискретизованном виде, что, как правило, и встречается на практике, приводит к неортогональности этих функций. Задача пользователя состоит в том, чтобы определить, в какой степени эта неортогональность искажает результаты расчетов коэффициентов. Однако для низших порядков функции эти искажения столь незначительны, что практически не влияют на показания приборов при современных уровнях точности.
Таблица A.1
Внимание! Любая ортогональная совокупность должна быть применена с осторожностью. Применительно к вещественным (реальным) данным (например, при математической обработке аберраций) отдельные виды ортогональных функций не пригодны для анализа.
Причина, по которой ортогональные функции, описываемые в этом приложении, полезны для практического анализа, заключаются в форме нескольких первых членов (терминов), а именно "вздутие", "острие", "наклон", "рефракция", описывающих изменения в сопряжениях, а также в использовании ряда полезных аберраций нижнего порядка: "кома", "астигматизм", "сферическая аберрация".
Универсально используемая совокупность Цернике (A.3.1), так же как альтернативная совокупность (A.3.3), содержит эти полезные члены. Совокупности (ряды) Лежандра (A.3.2) не включают в себя все эти термины и члены (в частности, речь идет о рефракции). Однако в более высоких порядках все термины и члены становятся полезными и значимыми на практике. К числу редко упоминаемых компонентов относятся Z20, L17 или B29.
A.2 Вычисление коэффициентов
При выполнении аппроксимации множества ортогональных функций {P0, P1, P2 ...} функцией f коэффициенты {a0, a1, a2 ...} ортогональных функций вычисляются по формуле
![]() где
A - площадь интегрирования;
dA - инкремент;
Значения
A.3 Определения, таблицы и иллюстрации, относящиеся к ортогональным функциям
Примечание -
.A.3.1.1 Общие сведения
Речь идет о классических полиномах Цернике, широко используемых в оптике. Это множество полиномов (рядов) ортогонально применительно к круглым поверхностям (преобразуемым в эллиптические).
Определения приведены в A.3.1.2. Значения
Таблица A.2
![]()
Цернике (Z0 - Z35)
Определяемая область: круг с радиусом, равным 1.
Полиномы Цернике
![]() Главный член
![]() Квадрат среднего квадратического значения полинома
![]() Функция Цернике (ортогональные ряды, используемые в оптике)
![]() Примечание - Среднее квадратическое значение функций Цернике дополняется функциями косинуса и синуса. У этих функций (m > 0) средний квадрат величины является половиной соответствующего полинома Цернике. У вращательно-симметричных функций (m = 0) средние квадратичные величины равны соответствующим полиномам.
A.3.2.1 Общие сведения
Существуют 10 полиномов Лежандра. Перемножая каждый из этих полиномов на x и y, получаем 20 полиномов (функций), ортогональных на квадратных поверхностях (преобразуемых в эллиптические).
Определения приведены в A.3.2.2. Значения
Таблица A.3
![]()
Лежандра (2D), ортогонального множеству (L0 - L20)
Определяемая область: центрированный квадрат со стороной длиной 2.
Полином Лежандра 10
![]() Главный член
.1 D квадрат среднего квадратического значения полинома
![]() 2 D полиномы Лежандра (ортогональные ряды, применяемые в оптике
L(p, q, x, y) = L(p, x)·L(q, y).
2 D средние квадратические значения полиномов
||L(p, q, x, y)||2 = 1/[(2p + 1)(2q + 1)].
A.3.3.1 Общие сведения
2D Лежандра функции имеют один недостаток: полиномы нижнего порядка не содержат полезных вращательно-симметричных полиномов, входящих в состав полиномов Цернике, в частности полиномов до Z8 (т.е. до 4-го порядка), включая наиболее желательный рефракционный член, т.е. (x2 + y2).
Этот недостаток может быть компенсирован формированием полиномов с заметной вращательной симметрией.
Невращательная природа симметрии квадратной поверхности делает полиномы более сложными, от 6-го порядка и выше. Однако до 4-го порядка (в том числе) полиномы имеют ту же форму, что и полиномы Цернике, но с отличными коэффициентами.
Однако, как свидетельствуют графики, незаметная вращательная симметрия становится видимой даже при изображении функций высшего порядка.
В A.3.3.2 даны определения; значения
Таблица A.4
Таблица A.5
![]()
множества (B0 - B35) для квадратных поверхностей на базе
множества Цернике (Z0 - Z35)
A.3.2.2 Определения
Обозначения: (n + m) = V: порядок функции;
(n - m) >= 0: m не более n;
(n - m) - четное.
Определяемая область: центрированный квадрат со стороной длиной
Главный член: функция Цернике с идентичным индексом.
(справочное)
B.1 Общие сведения
В данном приложении даны рекомендации по вычислению ортогональных полиномов, описывающих требуемые характеристики поверхностей произвольной формы с учетом необходимой весовой функции. Математически это не новый метод, но ранее он не применялся в оптической практике. В этом приложении используется термин функция, а не полином, поскольку этот метод применим к любым базисным функциям. Например, для точного моделирования дефокусировки можно пользоваться функцией
.Приводятся несколько примеров ортогональных полиномов без деталей вычислений, которые упрощаются с использованием матриц, с помощью выстраиваемых друг за другом полиномов, причем детерминированным путем (без итераций и аппроксимаций).
B.2 Предварительные рассмотрения
Прежде всего нужно решить, будут ли используемы функции/полиномы для сохранения данных или для их анализа (обработки).
a) Сохранение данных
В этом случае важно иметь в виду, чтобы функции были ортогональны или, по крайней мере, принимать во внимание частичную потерю ортогональности, в противном случае реконструкция волнового фронта, хоть и незначительно, но может оказаться ошибочной.
Часто в цифровых интерферометрах приходится иметь дело с дискретными зрачками. В этом случае приходится создавать специальное множество ортогональных функций для каждого типа зрачка.
b) Анализ данных
Потеря ортогональности может оказывать влияние на амплитуду вычисляемых функций.
Если функции "почти" ортогональны (например, функции Цернике на дискретном зрачке), то погрешность будет мала. Для кривизны "Z3", например, погрешность равна 1/n·2; для более высоких порядков она может быть больше.
B.3 Общий метод и определения
Даже не в системе прямоугольных координат определим x и y как основные координаты.
Сформируем ряды {f0(x, y), f1(x, y), ..., fn(x, y)}.
B.3.1 Основная поверхность
"зрачок"
Примечание - Весовая функция (см. B.3.2) может участвовать в этом определении.
w(x, y).
Амплитуда волнового фронта на поверхности зрачка.
B.3.3 Корреляция двух функций
f1(x, y), f2(x, y).
Корреляция:
.Краткая запись: C(1, 2) = C(f1f2).
B.4 Генерирование функций
B.4.1 Начальные условия
В оптике по очевидным причинам предпочтительно начинать ряды функций с {1, x, y} или даже {1, x, y, r2}.
Примечание - ИСО 10110-5 и ИСО 10110-14 определяют "центр" зрачка таким образом, что три первых функции всегда ортогональны, {r2} никогда не может быть ортогональна с {1}, поскольку, как и {1}, эта функция всегда положительна.
B.4.2 Свойства - примеры
B.4.2.1 Пример 1 - Функции (полиномы) Цернике
Zn должен иметь форму
, где n и p - целые числа. На практике косинусный член подвергается растяжению и записывается в виде . Это значительно упрощает вычисление. Полиномы Цернике наиболее просты для генерации.B.4.2.2 Пример 2 - "Квадратно-радиальные" функции
Рассмотрим пример:
, который также может быть переписан в виде![]() Член 2x введен в это выражение лишь для коррекции усредненного наклона основного члена
(степень 3 здесь введена с целью общей нормировки).На поверхности квадратного зрачка усредненный наклон базового члена вероятно будет другим. "Радиально-квадратные" функции будут идентичны функциям Цернике только в пределах базового члена, но в "корректирующих членах" будут отличаться от них.
"Квадратное" множество функций имеет серьезный недостаток, поскольку 4-симметричный квадратный зрачок будет всегда являться источником корреляции между
и членом с круговой симметрией. Но зрачковая (т.е. амплитудная) функция играет роль в распространении волнового фронта и его дифракции. Следовательно, некорректно говорить о циркулярно-симметричных функциях на поверхности циркулярно-асимметричного зрачка.Таким образом, в противоположность круглому зрачку (Цернике) нельзя сначала генерировать циркулярно симметричные функции, а затем добавлять члены
B.4.2.3 Пример 3 - "Квадратные" функции
Запишем, что основной член будет иметь форму {xp, yq} с корректирующими членами той же формы, но с парой индексов, "меньших", чем (p, q), при той же "значимости" (этого свойства вполне достаточно, так как симметричная квадратная поверхность гарантирует, что различные "значимости" будут устраняться при интегрировании по поверхности зрачка).
Основной ряд выглядит следующим образом:
{1, x, y, x2, xy, y2, x3, x2y, xy2, y3, x4, x3y и т.д.}.
Каждый порядок (p + q) содержит (p + q + 1) функций, т.е. так же, как и в функциях Цернике (Z8; Z11; Z12; Z16; Z17 и т.д.).
B.4.3 Генерация множества
Для каждого полинома следует отметить, что он не коррелирован с полиномами других генераций. Это обеспечивает возможность определять как главный, так и корректирующий члены. Если базовые члены выбраны правильно (а это значит, что между ними отсутствует взаимная корреляция), то и вся генерация выбрана правильно.
B.5 Дискретные (т.е. "реальные") поверхности
К дискретным поверхностям относится все сказанное ранее.
B.6 Матричное вычисление
Использование матриц облегчает процесс вычисления.
Вернуться в "Каталог нормативных документов"
Источник информации: https://internet-law.ru/documents/prod/gost-r_gosudarstvennyj-standart/41/gost_39214.html
На правах рекламы:
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||