RU2753627C1 - Способ определения волновых аберраций оптической системы - Google Patents

Abstract

Изобретение относится к области оптического приборостроения и касается способа определения волновых аберраций оптической системы. При осуществлении способа направляют световой пучок с длиной волны λ на оптическую систему и измеряют распределение интенсивности светового пучка в различных плоскостях в пространстве изображений. Определяют положение плоскости П_макс, в которой пиковая интенсивность имеет максимальное значение I_макс, и определяют контур изофоты в распределении интенсивности светового пучка в плоскости П_макс по уровню 0,2⋅I_макс и расстояние ρ максимального удаления точки с интенсивностью I_макс от контура изофоты. Измеряют расстояние L_макс от плоскости П_макс до плоскости П, соприкасающейся с последней поверхностью оптической системы. Далее измеряют распределение интенсивности светового пучка в координатах двух плоскостей П_1рег и П_2рег, расположенных от плоскости П на расстояниях, определяемых в зависимости от величин λ, ρ и L_макс. По полученным данным определяют волновые аберрации оптической системы w₁(x,y) и w₂(x,y) в координатах плоскости П. Значения волновых аберраций w(x,y) оптической системы определяют как w(x,y)=[w₁(x,y)+w₂(x,y)]/2. Технический результат заключается в повышении точности и достоверности определения волновых аберраций. 4 ил., 4 табл.

Description

Изобретение относится к области оптики и может быть использовано для характеристики качества оптических систем или оптических элементов при их изготовлении и эксплуатации.

В настоящее время существуют различные методы измерений волновых аберраций оптических систем, в том числе измерений распределения волновых аберраций w(x,y) по поперечному сечению пучка на выходе оптической системы. Особенности схемного и конструктивного решения оптической системы, связанные с ее назначением и условиями применения, а также свойства материалов оптических элементов системы могут значительно усложнять задачу определения w(x,y), особенно в случае применения методов, использующих распределение интенсивности в поперечном сечении прошедшего систему пучка. Связано это с тем, что оптическая система способна вносить в проходящий через нее световой пучок наряду с аберрационными, еще и амплитудные искажения.

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

В тех случаях, когда на пути распространения светового пучка за оптической системой формируется неоднородное распределение интенсивности по сечению пучка, являющееся следствием суперпозиции аберрационных и амплитудных искажений, определение волновых аберраций w(x,y) из данных по распределению интенсивности в поперечном сечении прошедшего систему пучка неизбежно приводит к значениям, которые отличаются от действительных. Под действительными значениями здесь и далее мы подразумеваем те значения волновых аберраций, которые вносит оптическая система в проходящий через нее световой пучок.

Известен способ определения волновых аберраций оптической системы [Патент РФ №2680657, МПК G01M 11/02, дата приор. 10.04.2018 г., опубл.25.02.2019 г., Бюлл.№6], принятый за прототип. При осуществлении способа направляют некогерентный гомоцентрический световой пучок с длиной волны X на оптическую систему. Определяют положение плоскости П_макс, в которой пиковая интенсивность светового пучка в распределении интенсивности светового пучка имеет максимальное значение I_макс. Измеряют расстояние L_макс от плоскости П_макс до плоскости П, перпендикулярной оси светового пучка и соприкасающейся с последней поверхностью оптической системы. Определяют контур изофоты в распределении интенсивности светового пучка в плоскости П_макс по уровню 0,2 I_макс и расстояние ρ максимального удаления точки О_макс с пиковой интенсивностью I_макс от контура изофоты. Регистрируют распределение интенсивности

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

плоскости П_рег, положение которой выбирают на удалении

от плоскости П_макс. Волновые аберрации w(x,y) определяют из уравнения, связывающего w(x,y) в координатах (х,у) плоскости П, с распределением интенсивности светового пучка в плоскости П_рег:

при этом используют распределение интенсивности светового пучка I(х,у), выраженное в координатах плоскости П путем замены переменных в распределении

Значения w(x,y) находят путем численного решения уравнения (1).

К причинам, препятствующим достижению указанного ниже результата при использовании известного способа, принятого за прототип, относятся следующие.

В способе прототипе уравнение (1) устанавливает в приближении геометрической оптики связь между волновыми аберрации w(x,y) на выходе оптической системы, т.е. в плоскости П, с распределением интенсивности по сечению пучка в плоскости П_рег. Однако, как указывалось выше, существуют оптические системы, которые вносят в проходящий световой пучок, помимо аберрационных, еще и амплитудные искажения, т.е. сразу на выходе формируют неоднородное распределение интенсивности по сечению пучка. Так, амплитудные искажения с параболическим профилем по сечению пучка и перепадом от центра к краю вплоть до нескольких десятков процентов на выходе объектива могут проявляться в светосильных системах, а поперечные градиенты интенсивности - во внеосевых. Эти искажения распределения интенсивности накладываются в плоскости регистрации на те, которые имеют аберрационную природу, в результате определяемые значения волновых аберраций существенно отличаются от действительных значений, и использование способа не обеспечивает получение достоверных значений волновых аберраций.

Сущность изобретения заключается в следующем.

Задачей изобретения является определение волновых аберраций оптической системы, формирующей на выходе неоднородное распределение интенсивности пучка, с высокой точностью и достоверностью, без ужесточения требований к реализации способа и усложнения обработки результатов измерений.

Технический результат, который может быть достигнут при осуществлении заявленного способа, заключается в определении волновых аберраций оптической системы, формирующей на выходе неоднородное распределение интенсивности пучка, с высокой точностью и достоверностью, без ужесточения требований к реализации способа и усложнения обработки результатов измерений.

Указанный технический результат при осуществлении изобретения достигается тем, что в способе определения волновых аберраций оптической системы, включающем направление некогерентного гомоцентрического светового пучка с длиной волны X на оптическую систему, измерение распределений интенсивности светового пучка в плоскостях, перпендикулярных оси светового пучка, прошедшего оптическую систему и сфокусированного в пространстве изображений, определение положения плоскости П_макс, перпендикулярной оси светового пучка в пространстве изображений, в которой пиковая интенсивность светового пучка в распределении интенсивности светового пучка имеет максимальное значение I_макс среди экспериментально зарегистрированных пиковых интенсивностей в распределениях интенсивностей светового пучка в плоскостях, перпендикулярных оси светового пучка, определение контура изофоты в распределении интенсивности светового пучка в плоскости П_макс по уровню 0,2⋅I_макс и расстояние ρ максимального удаления точки О_макс с пиковой интенсивностью I_макс от контура изофоты, измерение расстояния L_макс по оси пучка от плоскости П_макс до плоскости П, перпендикулярной оси светового пучка и соприкасающейся с последней поверхностью оптической системы, с учетом правила знаков, измерение распределения интенсивности I_1рег(ξ₁,η₁) светового пучка в координатах (ξ₁,η₁) плоскости П_1рег, перпендикулярной оси светового пучка в пространстве изображений, которая удалена от плоскости П_макс на расстояние L_1рег по оси пучка, при этом начало координат в плоскостях П, П_1рег, П_макс лежит на оси светового пучка, определение волновых аберраций оптической системы w₁(x,y) в координатах (х,у) плоскости П из уравнения

- I₁(x,y) - распределение интенсивности светового пучка, выраженное в координатах плоскости П, полученное путем замены переменных в распределении интенсивности I_1рег(ξ₁,η₁) по формулам ξ₁=x⋅(L_1рег/L_макс) и η₁=y⋅(L_1рег/L_макс);

- I₁(0,0) - значение функции I₁(x,y) на оси пучка,

путем численного решения его вариационным методом при свободных граничных условиях, аппроксимируя w₁(x,y) степенным полиномом, удовлетворяющим условию w₁(0,0)=0 на оси пучка, в соответствии с заявляемым техническим решением положение плоскости П] _рег выбирают между плоскостью П и плоскостью П_рег, перпендикулярной оси светового пучка в пространстве изображений и удаленной от плоскости П_макс на расстояние L_per=-16⋅ρ²/λ с учетом правила знаков, дополнительно измеряют распределение интенсивности L_2рег(ξ₂,η₂) светового пучка в координатах (ξ₂,η₂) плоскости П_2рег, перпендикулярной оси светового пучка, положение которой выбирают в пространстве изображений за плоскостью П_макс по ходу пучка из условия L_2рег=-L_1рег⋅[L_макс/(L_макс-2L_1рег)], где L_2рег - расстояние по оси пучка от плоскости П_макс до плоскости П_2рег, при этом начало координат в плоскости П_2рег лежит на оси светового пучка, определяют волновые аберрации оптической системы w₂(x,y) в координатах (х,у) плоскости П из уравнения

- I₂(х,у) - распределение интенсивности светового пучка, выраженное в координатах плоскости П, полученное путем замены переменных в распределении интенсивности I_2рег(ξ₂,η₂) по формулам ξ₂=x⋅(L_2рег/L_макс) и η₂=y⋅(L_2рег/L_макс);

- I₂(0,0) - значение функции I₂(х,у) на оси пучка,

путем численного решения его вариационным методом при свободных граничных условиях, аппроксимируя w₂(х,у) степенным полиномом, удовлетворяющим условию w₂(0,0)=0 на оси пучка, значения волновых аберраций w(x,y) оптической системы находят из соотношения w(x,y)=[w₁(x,y)+w₂(x,y)]/2, считая значения волновых аберраций w₁(x,y) и w₂(х,у) оценочными.

Совокупность вышеизложенных признаков изобретения связана причинно-следственной связью с техническим результатом изобретения.

Уравнение (1) связывает волновые аберрации w(x,y) в плоскости П, перпендикулярной оси светового пучка и соприкасающейся с последней поверхностью оптической системы, с распределением интенсивности светового пучка в плоскости регистрации (П_рег в способе прототипе). При этом предполагается, что условия геометрооптического приближения в плоскости регистрации выдерживаются и пучок имеет равномерное распределение интенсивности на выходе оптической системы в плоскости П.

Если же световой пучок в этой плоскости имеет неравномерное распределение интенсивности, порожденное свойствами оптической системы, то наличие амплитудных искажений неизбежно приводит к отклонениям от точного результата в сторону больших или, наоборот, меньших значений в зависимости от положения плоскости регистрации. Выбор положения плоскости регистрации крайне важен при проведении экспериментов по определению волновых аберраций. В самом деле, профиль распределения интенсивности на выходе оптической системы в отсутствие аберраций, обусловленный лишь амплитудными искажениями, сохраняется в поперечных сечениях светового пучка по пути его распространения (если пренебречь влиянием дифракции). Классический тому пример - распространение сфокусированного гауссовского пучка: в отсутствие аберраций во всех его поперечных сечениях как до, так и за плоскостью фокусировки сохраняется гауссовский профиль. Однако если оптическая система вносит аберрационные искажения в распределение интенсивности светового пучка, то степень неоднородности регистрируемого распределения интенсивности будет изменяться на дистанции распространения светового пучка.

Как было обнаружено в наших экспериментах, выраженные в координатах (х,у) плоскости П аберрационные искажения распределения интенсивности в сечении пучка, зарегистрированные в плоскости регистрации, расположенной в пространстве между плоскостью П и П_макс, имеют профиль, обратный по знаку к выраженному в тех же координатах профилю распределения интенсивности, которое зарегистрировано в плоскости регистрации, расположенной за плоскостью П_макс. Например, если под влиянием сферической аберрации в плоскости регистрации перед П_макс распределение интенсивности приобретает параболический профиль с максимумом в приосевой области, то в плоскости регистрации, расположенной за П_макс, также формируется параболический профиль, но с минимумом на оси пучка. Поэтому в тех участках пути распространения светового пучка, где амплитудные и аберрационные искажения распределения интенсивности имеют сходный профиль, они накладываются друг на друга и усиливают степень неоднородности суммарного распределения интенсивности. И, наоборот, в тех участках пути, где профили амплитудных и аберрационных искажений распределения интенсивности противоположны по отношению друг к другу, степень неоднородности распределения понижается. В обоих случаях это влияет на достоверность определения значений волновых аберраций и приводит к получению значений, отличающихся от действительных: в первом случае - в сторону увеличения значений волновых аберраций, во втором - в сторону уменьшения значений.

В заявляемом способе регистрация распределений интенсивности светового пучка осуществляется в двух плоскостях П_1рег и П_2рег, расположенных перед и за плоскостью П_макс по ходу пучка, соответственно. Расстояния удаления плоскостей регистрации от плоскости П_макс обозначим как L_1рег и L_2рег, соответственно.

В способе прототипе было получено предельное расстояние удаления плоскости регистрации изображения от плоскости П_макс для аберрированного пучка, равное 16ρ²/λ по обе стороны от П_макс, определяющее зону пространства, вне которой еще можно рассматривать распространение пучка в геометрооптическом приближении.

Положение плоскости П_рег мы выбираем на расстоянии L_рег=-16⋅ρ²/λ от плоскости П_макс с учетом правила знаков. Тогда с учетом упомянутых выше ограничений положение плоскости П_1рег можно выбирать перед плоскостью П_макс в промежутке между плоскостями П и П_рег. Определяя положение плоскости П_1рег в экспериментах с реальными оптическими системами, следует руководствоваться следующим. С одной стороны, условия геометрооптического приближения выполняются тем лучше и, следовательно, тем меньше проявляются погрешности расчета, связанные с дифракцией, чем более удалены положения плоскостей регистрации от плоскости П_макс. С другой стороны, уменьшается амплитуда неоднородности распределения интенсивности аберрационного происхождения и, соответственно, ее относительный вклад в сравнении с амплитудными искажениями распределения и случайными флуктуациями интенсивности.

Как показали наши эксперименты и последующий анализ, выбранное положение плоскости П_1рег задает положение второй плоскости регистрации П_2рег. Было обнаружено, что оптимальное для точного определения волновых аберраций расстояние L_2рег нужно выбирать из условия:

Уравнения (2) и (3) связывают волновые аберрации искаженного пучка в координатах (х,у) плоскости П с распределениями интенсивности светового пучка в плоскостях П_1рег и П_2рег. Используя эти уравнения, можно определить значения волновых аберраций w₁(x,y) и w₂(x,y) в плоскости П по распределениям интенсивности светового пучка в плоскостях П_1рег и П_2рег, соответственно. С учетом того, что значения функций w₁(x,y) и w₂(x,y) получены с отличиями от действительных в сторону их увеличения или уменьшения, мы будем считать их оценочными значениями волновых аберраций оптической системы.

В заявляемом способе значения волновых аберраций w(x,y) оптической системы находят из соотношения w(x,y)=[w₁(x,y)+w₂(x,y)]/2. При арифметическом усреднении функций w₁(x,y) и w₂(x,y), предполагающем их сложение, отличия от действительных значений будут взаимно скомпенсированы. Функция волновых аберраций w(x,y) имеет такую же, как w₁(x,y) и w₂(х,у), структуру в виде линейной комбинации базовых функций, описывающих отдельные виды аберраций.

Таким образом, регистрация распределений интенсивности светового пучка в двух плоскостях, расположенных по разные стороны от П_макс на определенном удалении, определение по ним значений волновых аберраций, рассматриваемых нами как оценочные, последующее арифметическое усреднение полученных значений позволяет экспериментально отделить вклад присутствующих в них амплитудных искажений и получить действительные значения волновых аберраций с высокой точностью, не уступающей точности определения w(x,y) с равномерным распределением интенсивности на выходе оптической системы. Как будет показано ниже, проведенные нами эксперименты подтвердили высокую степень точности и достоверности определяемых значений волновых аберраций. Способ осуществляется без ужесточения требований к его реализации и усложнения обработки результатов измерений.

На фиг. 1 приведена оптическая схема установки для измерения волновых аберраций линзового объектива, где 1 - некогерентный источник света (лампа накаливания); 2 - конденсор; 3 - узкополосный светофильтр; 4 - микродиафрагма; 5 - многокомпонентный линзовый объектив; 6 - ПЗС-матрица, регистрирующая распределение интенсивности в поперечном сечении светового пучка. Плоскость П - плоскость, в системе координат которой определяют волновые аберрации, соприкасается с последней поверхностью линзового объектива 5, плоскость П_1рег - плоскость, в которой регистрируют распределение интенсивности I_1рег(ξ₁,η₁), плоскость П_макс - плоскость, в которой пиковая интенсивность светового пучка достигает максимального значения I_макс среди экспериментально зарегистрированных пиковых интенсивностей в распределениях интенсивностей светового пучка в плоскостях, перпендикулярных оси светового пучка, О_макс - точка пиковой интенсивности I_макс в плоскости П_макс, П_2рег - плоскость, в которой регистрируют распределение интенсивности L_2рег(ξ₂,η₂). Плоскость П_рег удалена от плоскости П_макс на расстояние L_рег=-16⋅ρ²/λ с учетом правила знаков, где ρ - расстояние максимального удаления точки О_макс от изофоты распределения интенсивности в плоскости П_макс на уровне 0,2⋅I_макс. Плоскости П_1рег, П_2рег, П_рег и П_макс перпендикулярны оси светового пучка, начало координат в плоскостях П, П_1рег, П_2рег, П_рег и П_макс лежит на оси светового пучка. L_1рег, L_2рег и L_макс - расстояния от плоскости П_макс до плоскостей П_1рег, П_2рег и П, соответственно, по оси пучка с учетом правила знаков.

На фиг. 2, А и Б, представлены экспериментально зарегистрированные изображения поперечного сечения пучка в плоскостях П_1рег и П_2рег, отстоящих от плоскости П_макс на расстояния L_1рег=-20 мм и L_2рег=25 мм, соответственно.

На фиг. 2, В и Г, представлены приведенные к координатам в плоскости П экспериментальные распределения интенсивности светового пучка вдоль горизонтального диаметра в плоскостях П1_рет и П_2рет, отстоящих от плоскости П_макс на расстояния L_1рег=-20 мм и L_2рег=25 мм, соответственно.

На фиг. 3 представлены в виде гистограммы нормированные на длину волны среднеквадратические значения волновых аберраций σ_w1/λ, σ_w2/λ и σ_w/λ на выходе оптической системы, соответствующие эксперименту по определению волновых аберраций w₁(x,y), w₂(x,y) и w(x,y). На горизонтальной оси гистограммы приведены соответствующие эксперименту значения L_1рег=-20 мм и L_2рег=25 мм. По вертикальной оси гистограммы высота левого столбца представляет величину σ_w1/λ, - нормированное среднеквадратическое значение волновых аберраций w₁(x,y), вычисленных по распределению пучка в плоскости П_1рег. Высота правого столбца представляет величину σ_w2/λ, - нормированное среднеквадратическое значение волновых аберраций w₂(x,y), вычисленных по распределению пучка в плоскости П_2рег. Высота среднего столбца σ_w/λ гистограммы соответствует нормированному среднеквадратическому значению волновых аберраций w(x,y).

На фиг. 4 представлены в виде трех гистограмм нормированные на длину волны среднеквадратические значения волновых аберраций σ_w1/λ, σ_w2/λ и σ_w/λ на выходе оптической системы, соответствующие трем проведенным экспериментам. В каждом из экспериментов определялись волновые аберрации w₁(x,y), w₂(x,y) и w (х,у). Эксперименты отличались выбором положения плоскостей П_1рег и П_2рег. На горизонтальной оси каждой гистограммы приведены соответствующие эксперименту значения L_1рег и L_2рег. По вертикальной оси каждой гистограммы высота левого столбца представляет величину σ_w1/λ - нормированное среднеквадратическое значение волновых аберраций w₁(x,y), вычисленных по распределению пучка в плоскости П_1рег. Высота правого столбца представляет величину σ_w2/λ - нормированное среднеквадратическое значение волновых аберраций w₂(x,y), вычисленных по распределению пучка в плоскости П_2рег. Высота среднего столбца σ_w/λ каждой гистограммы соответствует нормированному среднеквадратическому значению волновых аберраций w(x,y), рассчитанных для соответствующего эксперимента.

Заявленный способ был экспериментально осуществлен для измерения волновых аберраций многокомпонентного линзового объектива. Для этого была создана экспериментальная установка.

Пример конкретного выполнения способа.

Излучение от лампы накаливания 1 проецируется конденсором 2 на отверстие микродифрагмы 4 диаметром 50 мкм через узкополосный светофильтр 3 с максимумом пропускания на длине волны 0,55 мкм (фиг. 1). За отверстием микродиафрагмы, имитирующей точечный источник излучения, формируется сферический пучок с однородным распределением интенсивности в пределах входной апертуры исследуемого объектива. Такой пучок можно считать гомоцентрическим. Микродиафрагма 4 была установлена на расстоянии 870 мм от первой оптической поверхности объектива 5 со смещением в горизонтальном направлении относительно центра поля зрения объектива. В такой геометрии эксперимента прошедший оптическую систему световой пучок приобретает волновые аберрации в основном типа комы и сферической. Диаметр пучка на выходе оптической системы составлял 20 мм.

ПЗС-матрицу 6 автоколлимационным методом выставляют ортогонально оси светового пучка, прошедшего объектив 5. Затем определяют местоположение плоскости П_макс на оси пучка путем перемещения ПЗС-матрицы вдоль каустики светового пучка, сравнения пиковых значений в распределении интенсивности в разных позициях по ходу пучка и фиксации положения ПЗС-матрицы, в котором пиковая интенсивность излучения имеет максимальное значение I_макс. ПЗС-матрицу при этом ориентируют ортогонально оси пучка. Используя систему отсчета, задаваемую системой координат в плоскости П_макс, в которой начало координат лежит на оси светового пучка, а положительное направление оси аппликат совпадает с направлением распространения светового пучка, измеряют расстояние L_макс от плоскости П_макс до плоскости П, перпендикулярной оси светового пучка и соприкасающейся с последней поверхностью объектива 5. В нашем эксперименте L_макс=-200 мм.

В плоскости П_макс регистрируют распределение интенсивности в поперечном сечении пучка. Программными средствами из комплекта технического обеспечения ПЗС-матрицы определяют контур изофоты на уровне сигнала 0,2 от пикового значения I_макс и находят расстояние ρ максимального удаления точки О_макс пиковой интенсивности I_макс от этой границы. Определяют положение плоскости П_рег, перпендикулярной оси светового пучка и удаленной от плоскости П_макс на расстояние L_рег=-16⋅ρ²/λ с учетом правила знаков. Выбирают расположение плоскости П_1рег в пространстве между плоскостями П и П_рег, этому положению будет соответствовать расстояние L_1рег, определяющее удаление плоскости П_1per от плоскости П_макс. В нашем эксперименте ρ=23 мкм, L_рег=-15 мм, L_1рег=-20 мм. В плоскости П_1per регистрируют ПЗС-матрицей распределение интенсивности I_1рег(ξ₁,η₁) светового пучка. Зарегистрированное изображение поперечного сечения пучка в этой плоскости представлено на фиг. 2,А. Полученное распределение затем преобразуют путем замены координат по формулам ξ₁=x⋅(L_1рег/L_макс) и η₁=y⋅(L_1рег/L_макс) в распределение интенсивности I₁(x,y) в координатах плоскости П. На фиг. 2,В представлено приведенное к координатам плоскости П распределение интенсивности светового пучка вдоль горизонтального диаметра I₁(x,0).

Определяют волновые аберрации оптической системы w₁(x,y) в плоскости П из уравнения (2), в котором I₁(0,0) - значение функции I₁(x,y) на оси пучка, путем численного решения уравнения вариационным методом при свободных граничных условиях. При этом для аппроксимации волновых аберраций в качестве базовых функций

составляющих сумму w₁(x,y), выбирают степенные полиномы 3-й и 4-й степени. Каждая из базовых функций

отвечает определенному типу волновых аберраций и удовлетворяет условию

Так, волновым аберрациям третьего порядка типа комы будут отвечать функции вида

а сферической аберрации - функция вида

Значения аберрационных коэффициентов A_q являются решением системы линейных уравнений [Л.В.Канторович, В.И. Крылов. Приближенные методы высшего анализа. М.-Л.: Гос. Изд-во физ.-мат.лит.1962. - 708 с.]:

где коэффициенты

вычисляют по формуле:

коэффициенты β_m вычисляют по формуле:

Здесь

- это правая часть уравнения (2), т.е.

область интегрирования ∑ - есть поперечное сечение пучка за исключением узкой краевой зоны поперечного сечения, где распределение интенсивности несколько возмущено дифракцией пучка на апертурной диафрагме оптической системы.

С учетом дискретного характера экспериментальных данных и замены координат (х,у) на дискретные переменные по формулам

базовые функции

преобразуют к виду:

- нечетная функция от координаты х_к описывает кому 3-го порядка,

- нечетная функция от координаты у_к описывает кому 3-го порядка,

- четная функция от координат x_k1 и y_n1 описывают сферическую аберрацию 3-го порядка.

Здесь

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

δ - размер пиксела ПЗС-матрицы;

k₁=0,…k_j и n_i=0, …N_i - индексы нумерации пикселей ПЗС-матрицы;

k₁₀ и n₁₀ - индексы оси пучка (или центра поперечного сечения пучка) в кадре.

Аппроксимирующая волновые аберрации в целом функция w₁(x,y), выраженная в дискретных координатах, имеет вид:

Используя указанные базовые функции и выраженное в дискретных координатах распределение интенсивности по сечению пучка I₁(x_k,y_n), находят коэффициенты α_q,m и β_m системы уравнений (5), заменяя операцию интегрирования суммированием по индексам пикселей ПЗС-матрицы. Решая алгебраическую систему уравнений (5), вычисляют искомые аберрационные коэффициенты А_3х, А_3у и А₄. В проведенном эксперименте были получены следующие значения аберрационных коэффициентов: А_3х=-8,42⋅10^-7 мм^-2, А_3у=1,48⋅10^-7 мм^-2 и A₄=8,06⋅10^-8 мм^-3.

Поскольку на длине шага координатной сетки изменения величины волновой аберрации очень малы, то можно пренебречь дискретностью координат и рассматривать волновые аберрации для каждого эксперимента как непрерывные функции координат х и y.

На следующем этапе измеряют распределение интенсивности I_2рег(ξ₂,η₂) светового пучка в координатах (ξ₂,η₂) плоскости П_2рег, положение которой задают в пространстве изображений за плоскостью П_макс по ходу пучка из условия L_2рег=-L_1рег⋅[L_макс/(L_макс-2L_1рег)], где L_2рег - расстояние по оси пучка от плоскости П_макс до плоскости П_2рег. Выбранному значению L_1рег=-20 мм отвечает значение L_2рег=25 мм. На фиг. 2, Б представлено зарегистрированное изображение поперечного сечения пучка в плоскости П_2рег. Распределение I_2рег(ξ₂,η₂) преобразуют путем замены координат по формулам ξ₂=x⋅(L_2рег/L_макс) и η₂=y⋅(L_2рег/L_макс) в распределение интенсивности I₂(х,у) в координатах плоскости П. На фиг. 2,Г представлено приведенное к координатам плоскости П распределение интенсивности светового пучка вдоль горизонтального диаметра I₂(х,0).

Из сравнения распределений интенсивности I₁(x,0) и I₂(х,0), представленных на фиг. 2,В и Г, отчетливо видно, что под влиянием сферической аберрации в плоскости регистрации П_1рег распределение интенсивности приобретает параболический профиль с максимумом в приосевой области, а в плоскости П_2рег - параболический профиль с минимумом на оси пучка.

Измеренные распределения интенсивности используют для определения волновых аберраций оптической системы w₂(x,y) в координатах (х,у) плоскости П из уравнения (3), в котором I₂(0,0) - значение функции I₂(х,у) на оси пучка, путем численного решения уравнения (3) вариационным методом при свободных граничных условиях. При этом для аппроксимации волновых аберраций w₂(x,y) выбирают те же степенные полиномы 3-й и 4-й степени, что были использованы при аппроксимации волновых аберраций w₁(x,y), но с заменой координат (х,у) на дискретные переменные по формулам

- шаг координатной сетки в плоскости П,

k₂=0,…k₂ и n₂=0,…N₂ - индексы нумерации пикселей ПЗС-матрицы;

k₂₀ и n₂₀ - индексы оси пучка (или центра поперечного сечения пучка) в кадре.

В этом случае аппроксимирующая волновые аберрации функция w₂(x,y), выраженная в дискретных координатах, имеет вид:

Искомые аберрационные коэффициенты В_3х, В_3у, В₄ также находят описанным выше методом, аналогично вычислению аберрационных коэффициентов А_3х, А_3у и А₄. В проведенном эксперименте были получены следующие значения аберрационных коэффициентов: В_3х=-2,11⋅10^-7 мм^-2, В_3у=4,23⋅10^-7 мм^-2 и В₄=7,94⋅10^-8 мм^-3.

Пренебрегая, как и выше, дискретностью координат, полученные волновые аберрации будем рассматривать как непрерывные функции координат х и у:

Как уже указывалось выше, значения волновых аберраций w₁(x,y) и w₂(x,y) считают оценочными. На последнем этапе определяют волновые аберрации w(x,y) путем арифметического усреднения найденных функций волновых аберраций: w(x,y)=[w₁(x,y)+w₂(x,y)]/2. В силу линейности операции арифметического усреднения, w(x,y) также имеет вид суммы базовых функций, характеризующих отдельные типы аберраций. Аберрационные коэффициенты для определения волновых аберраций w(x,y), рассчитанные по формулам

имеют значения

а искомые волновые аберрации представлены выражением:

Волновые аберрации, включающие в себя несколько аберрационных компонент, обычно количественно характеризуют среднеквадратическим значением, вычисляемым по формуле:

где S - площадь поперечного сечения светового пучка на выходе оптической системы.

На фиг. 3 представлены в виде гистограммы нормированные на длину волны среднеквадратические значения волновых аберраций

соответствующие эксперименту, в котором значения L_1рег=-20 мм и L_2рег=25 мм. По вертикальной оси гистограммы высота левого столбца представляет величину

- нормированное среднеквадратическое значение волновых аберраций w₁(x,y), вычисленных по распределению пучка в плоскости П_1рег. Высота правого столбца представляет величину

- нормированное среднеквадратическое значение волновых аберраций w₂(x,y), вычисленных по распределению пучка в плоскости П_2рег. Высота среднего столбца

гистограммы соответствует нормированному среднеквадратическому значению волновых аберраций w(x,y). Видно, что значение

заметно отличается от значения

(на десятые доли λ), что говорит о присутствии амплитудных искажений пучка оптической системой в эксперименте. При этом, как и следовало ожидать, отклонение значений

и

от значения

имеет разный знак.

Таким образом, осуществление операций способа позволило устранить влияние амплитудных искажений и получить действительные значения волновых аберраций.

Для подтверждения высокой точности и достоверности полученного результата были проведены дополнительные эксперименты. В них положение плоскости П_1рег выбиралось на пути распространения светового пучка как до, так и после положения, соответствующего описанному выше основному эксперименту, значения L_1рег составляли, соответственно, L_1рег=17 мм для первого дополнительного эксперимента и L_1рег=-25 мм для второго дополнительного эксперимента.

Аналогичным образом, как и в описанном выше основном эксперименте, были получены значения аберрационных коэффициентов А_3х, А_3у и А₄ для двух дополнительных экспериментов. Эти значения приведены в Табл.1 (для наглядности вместе со значениями, соответствующими значению L_1рег=-20 мм).

Выбранному значению L_1рег=-17 мм в первом дополнительном эксперименте отвечает значение L_2рег=20 мм, значению L_1рег=-25 мм во втором дополнительном эксперименте - значение L_2рег=33 мм.

Далее аналогичным образом, как и в описанном выше основном эксперименте, были получены значения аберрационных коэффициентов В_3х, В_3у, В₄ для двух дополнительных экспериментов. Эти значения приведены в Табл.2 (для наглядности вместе со значениями, соответствующими значению L_2рег=25 мм).

Аналогичным образом, как и в описанном выше основном эксперименте, были получены значения аберрационных коэффициентов С_3х, С_3у и С₄ для двух дополнительных экспериментов. Эти значения приведены в Табл. 3 (для наглядности вместе со значениями, полученными в основном эксперименте).

На фиг. 4 представлены в виде трех гистограмм нормированные на длину волны среднеквадратические значения волновых аберраций σ_w1, σ_w2 и σ_w, полученные в основном и двух дополнительных экспериментах для трех положений плоскостей П_1рег и П_2рег. Их числовые значения приведены в таблице 4.

Из фиг. 4 и таблицы 4 видно, что значения σ_w1 и σ_w2 в каждом из выбранных положений П_1рег и П_2рег сильно отличаются - их разница составляет ~0.2λ÷0.3λ, причем во всех случаях σ_w1/λ > σ_w2/λ. Такое различие свидетельствует о том, что световой пучок на выходе оптической системы имеет неоднородное по сечению распределение интенсивности. Сравнение значений σ_w1 для трех разных положений П_1рег и П_2рег показывает, что они также отличаются. То же можно сказать и о значениях σ_w2. В противоположность этому разброс значений σ_w/λ вычисленных для каждого из трех положений П_1рег и П_2рег, незначителен и не превышает величины λ/50.

Поскольку в каждом из трех проведенных экспериментов реализуется своя независимая суперпозиция вкладов амплитудных и аберрационных искажений в неоднородности распределения интенсивности пучка и, в то же время, результаты всех трех измерений волновых аберраций, оцениваемые среднеквадратическим значением, оказываются чрезвычайно близки, то это говорит о высокой степени их точности и достоверности.

Отметим также, что при определении волновых аберраций заявленным способом остаются неизменными приемы настройки схемы измерений и используемая аппаратура регистрации изображений светового пучка, применяются те же самые математические методы и алгоритмы обработки изображений, что и при реализации способа-прототипа. Поэтому можно утверждать, что по сравнению с прототипом ужесточения требований к реализации заявленного способа и усложнения обработки результатов измерений не происходит.

Claims (9)

1. Способ определения волновых аберраций оптической системы, включающий направление некогерентного гомоцентрического светового пучка с длиной волны λ на оптическую систему, измерение распределений интенсивности светового пучка в плоскостях, перпендикулярных оси светового пучка, прошедшего оптическую систему и сфокусированного в пространстве изображений, определение положения плоскости П_макс, перпендикулярной оси светового пучка в пространстве изображений, в которой пиковая интенсивность светового пучка в распределении интенсивности светового пучка имеет максимальное значение I_макс среди экспериментально зарегистрированных пиковых интенсивностей в распределениях интенсивностей светового пучка в плоскостях, перпендикулярных оси светового пучка, определение контура изофоты в распределении интенсивности светового пучка в плоскости П_макс по уровню 0,2⋅I_макс и расстояния ρ максимального удаления точки О_макс с пиковой интенсивностью I_макс от контура изофоты, измерение расстояния L_макс по оси пучка от плоскости П_макс до плоскости П, перпендикулярной оси светового пучка и соприкасающейся с последней поверхностью оптической системы, с учетом правила знаков, измерение распределения интенсивности I_1рег(ξ₁,η₁) светового пучка в координатах (ξ₁,η₁) плоскости П_1рег, перпендикулярной оси светового пучка в пространстве изображений, которая удалена от плоскости П_макс на расстояние L_1рег по оси пучка, при этом начало координат в плоскостях П, П_1рег, П_макс лежит на оси светового пучка, определение волновых аберраций оптической системы w₁(x,y) в координатах (х,у) плоскости П из уравнения
2. 
3. где I₁(x,y) - распределение интенсивности светового пучка, выраженное в координатах плоскости П, полученное путем замены переменных в распределении интенсивности I_1рег(ξ₁,η₁) по формулам ξ₁=x⋅(L_1рег/L_макс) и η₁=y⋅(L_1рег/L_макс);
4. I₁(0,0) - значение функции I₁(x,y) на оси пучка,
5. путем численного решения его вариационным методом при свободных граничных условиях, аппроксимируя w₁(x,y) степенным полиномом, удовлетворяющим условию w₁(0,0)=0 на оси пучка, отличающийся тем, что положение плоскости П_1рег выбирают между плоскостью П и плоскостью П_рег, перпендикулярной оси светового пучка в пространстве изображений и удаленной от плоскости П_макс на расстояние L_per=-16⋅ρ²/λ с учетом правила знаков, дополнительно измеряют распределение интенсивности L_2рег(ξ₂,η₂) светового пучка в координатах (ξ₂,η₂) плоскости П_2рег, перпендикулярной оси светового пучка, положение которой выбирают в пространстве изображений за плоскостью П_макс по ходу пучка из условия L_2рег=-L_1рег⋅[L_макс/(L_макс-2L_1рег)], где L_2рег - расстояние по оси пучка от плоскости П_макс до плоскости П_2рег, при этом начало координат в плоскости П_2рег лежит на оси светового пучка, определяют волновые аберрации оптической системы w₂(x,y) в координатах (х,у) плоскости П из уравнения
6. 
7. где I₂(х,у) - распределение интенсивности светового пучка, выраженное в координатах плоскости П, полученное путем замены переменных в распределении интенсивности I_2рег(ξ₂,η₂) по формулам ξ₂=x⋅(L_2рег/L_макс) и η₂=y⋅(L_2рег/L_макс);
8. I₂(0,0) - значение функции I₂(х,у) на оси пучка,
9. путем численного решения его вариационным методом при свободных граничных условиях, аппроксимируя w₂(x,y) степенным полиномом, удовлетворяющим условию w₂(0,0)=0 на оси пучка, значения волновых аберраций w(x,y) оптической системы находят из соотношения w(x,y)=[w₁(x,y)+w₂(x,y)]/2, считая значения волновых аберраций w₁(x,y) и w₂(x,y) оценочными.

Images