Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

itmo373

.pdf
Скачиваний:
29
Добавлен:
21.03.2016
Размер:
1.43 Mб
Скачать

h(η)

1.0

τ(ρ)

1.01

2

1

 

2

3

 

3

 

 

 

η

 

 

 

ρ

-1.12 -0.61

0

0.61 1.12

 

а) функция пропускания по зрачку

б) функция рассеяния точки

Рис.9.3.2. Влияние неравномерности пропускания по зрачку на ФРТ.

9.3.3. Безаберационная ОПФ. Предельная пространственная частота

Оптическая передаточная функция вычисляется при помощи выражения для автокорреляции зрачковой функции (9.2.15). Для безаберрационной оптической системы волновая аберрация W = 0 , тогда интеграл автокорреляции будет выглядеть следующим образом:

D(ωx ,ωy )=

1

∫∫dρxdρy =

Ω(ωx ,ωy )

(9.3.6)

 

Ω0

 

 

Ω0 Ω(ωx ,ω y )

 

где Ω ωx ,ωy – область интегрирования, показанная на рис.9.3.3.

Ω(ωx ,ωy )

ρy

 

Ω0 (ρx ωx ,ρy ωy )

 

 

 

 

 

ωy

ρx

 

 

 

0ωx

Ω0 (ρx ,ρy )

Рис.9.3.3. Области зрачков, смещенные относительно друг друга на ωx ,ωy .

Таким образом, безаберрационная ОПФ пропорциональна площади

 

 

, которая является функцией

перекрытия двух зрачков Ω ωx ,ωy

пространственных частот. Из рис.9.3.3 следует, что максимальная каноническая пространственная частота ωmax = 2 . Для более высоких частот

площадь Ω ωx ,ωy становится нулевой (рис.9.3.4).

140

Максимальной канонической пространственной частоте соответствуют предельные реальные пространственные частоты:

 

2A

2Ay

 

νlim x =

x

, νlim y =

 

(9.3.7)

 

λ

 

λ

 

ОПФ

1

идеальная о.с.

безаберрационная о.с.

ω

0

2

 

Рис.9.3.4. Безаберрационная ОПФ.

Таким образом, ОПФ реальной оптической системы при отсутствии аберраций не соответствует ОПФ идеальной оптической системы, и всегда ограничена предельными частотами, обусловленными дифракцией света.

9.4.Критерии качества оптического изображения

9.4.1.Предельная разрешающая способность по Релею

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

Предельная разрешающая способность – это минимальное расстояние

σR между двумя точками, при котором их изображение отличимо от изображения одной точки.

Критерий Релея гласит, что при провале в распределении интенсивности в изображении двух близких точек в 20% точки будут восприниматься как раздельные. Для этого необходимо, чтобы центральный максимум в изображении одной точки приходился бы на первый минимум в изображении другой (рис.9.4.1). Для оптических систем при отсутствии аберраций σR = 0.61

канонических единиц.

Разрешение по Релею удовлетворительно характеризует качество изображения астрономических телескопов, спектральных приборов, для которых предметами являются близко расположенные точки или линии, а также визуальных приборов (предназначенных для работы с глазом).

141

20%

η

σR

Рис.9.4.1. Разрешение по Релею.

9.4.2. Разрешающая способность по Фуко

Критерий Фуко применяется для оценки качества изображения оптических систем, передающих объекты сложной структуры. Разрешающая способность R определяется как максимальная пространственная частота периодического тест-объекта, состоящего из черно-белых штрихов (миры Фуко), в изображении которого еще различимы штрихи. Разрешающую способность обычно определяют для миры единичного (абсолютного) контраста по графику ЧКХ оптической системы (рис.9.4.2). Разрешающая способность R определяется для заданного контраста (обычно для контраста k′= 0.2 ).

k

1

0.2

 

 

ω

 

 

 

0

 

R

R0 = 2

 

 

 

 

Рис.9.4.1. Разрешающая способность по Фуко.

Предельная разрешающая способность R0 для оптических систем

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

9.5. Влияние аберраций на ФРТ и ОПФ

Влияние малых аберраций (волновая аберрация составляет доли длин волн) на ФРТ проявляется в том, что часть энергии из центрального максимума переходит в кольца. В результате в центральном максимуме остается около 6070% вместо 84%, при этом размеры центрального максимума сохраняются, а интенсивность в центре уменьшается (рис.9.5.1).

142

h(η)

1.0

η

1.22

Рис.9.5.1. Влияние аберраций на ФРТ.

Аберрации разных типов по-разному влияют на вид пятна рассеяния (картину Эри). В случае симметричных аберраций (расфокусировка, сферическая) сохраняется радиальная симметрия пятна (рис.9.5.2.а). В случае несимметричных аберраций (кома, астигматизм) симметрия пятна нарушается

(рис.9.5.2.б, рис.9.5.2.в).

ηy

ηy

ηy

ηx

ηx

ηx

а) симметричные аберрации

б) кома

в) астигматизм

Рис.9.5.2. Картины Эри для аберраций различных типов.

При дальнейшем увеличении аберраций сходство ФРТ с безаберрационной полностью теряется, и ее форма определяется картиной поперечных аберраций (точечной диаграммой). Практически вся энергия из центрального максимума перекачивается в кольца (в центральном максимуме остается меньше 40% энергии). Однако при этом сохраняется дифракционный узор с шагом 0.5 в канонических координатах.

9.5.1. Число Штреля

Поскольку при малых аберрациях часть энергии из центрального максимума перекачивается в кольца, уменьшается интенсивность в центральном максимуме. Обозначим значение ФРТ в ее максимуме при отсутствии аберраций h0 (0), а при наличии аберраций h(0) (рис.9.5.3).

143

h(η)

1.0

h(0) h0 (0)

η

Рис.9.5.3. Число Штреля.

Число Штреля (критерий Штреля, Strehl ratio) показывает влияние аберраций на ФРТ:

St =

h(0)

(9.5.1)

h (0)

 

0

 

 

Значение числа Штреля находится в пределах 0 St 1, энергия в кольца перекачивается в таком же соотношении. Если St =1 – оптическая система безаберрационная, если St 0.8 – система практически безаберрационная, поскольку уменьшение центрального максимума на 20% почти незаметно.

9.5.2. Критерий Релея для малых аберраций

Критерий или допуск Релея заключается в том, что если величина волновой аберрации (при условии, что в системе присутствует только

сферическая аберрация) не превосходит λ4 (рис.9.5.4), то число Штреля

St 0.8 . Отсюда Релей распространил свой критерий и на другие типы аберраций.

ρ

W

λ4

Рис.9.5.4. Величина волновой аберрации.

Релеевский допуск на остаточные аберрации:

Wmax < λ

4

(9.5.2)

 

 

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

144

сколько изменится число Штреля при Wmax = λ6 или Wmax = λ3 . Релеевский допуск точного ответа на этот вопрос не дает.

9.5.3. Формула Марешаля. Допуск Марешаля для малых аберраций

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

Рассмотрим вывод формулы Марешаля. Функция рассеяния точки, в соответствии с выражением (9.2.10):

h(ηx ,ηy )=

 

F 1[f (ρx , ρy )]2 =

 

+∫∫

f (ρx , ρy )e2πi(ηx ρx +ηy ρ y )dρxdρy

 

2

 

 

 

 

 

 

(9.5.3)

 

 

 

 

 

 

−∞

 

 

 

Значение ФРТ в ее центральном максимуме:

h(0)= h(ηx = 0,ηy = 0)=

 

+∫∫

f (ρx , ρy ) dρxdρy

 

2

 

 

 

 

(9.5.4)

 

 

−∞

 

 

 

Воспользовавшись выражением для зрачковой функции (9.2.3), получим:

h(0)=

 

∫∫e2πiW (ρx ,ρy )dρx dρy

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(9.5.6)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ω0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В случае малых аберраций W <

1

4

, следовательно

 

W

 

3

,

 

W

 

4

, ... <<1. Тогда

 

 

 

 

при разложении функции e2πiW

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в ряд, можно оставить только три члена, а

 

 

 

x

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

остальные отбросить e

 

=1+ x +

 

x

 

+K , отсюда:

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e2πiW 1+ 2πiW +

1

(2πiW )2 =1+ 2πiW 2π 2W 2

 

 

 

 

 

(9.5.7)

2

 

 

 

 

 

Тогда можно записать приближенное выражение для ФРТ:

h(0)

 

∫∫[1+ 2πiW 2π 2W 2 ]dρxdρy

 

2

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

Ω0

 

 

 

 

(9.5.8)

=

 

∫∫dρxdρy + 2πi ∫∫W (ρx , ρy )dρxdρy 2π 2 ∫∫W 2 (ρx , ρy )dρxdρy

 

2

 

 

 

 

 

 

 

 

 

Ω0

 

Ω0

 

Ω0

 

 

 

Введем обозначение для среднего значения волновой аберрации по зрачку:

 

 

 

Ω

∫∫

 

 

 

W

=

1

 

W (ρx , ρy )dρxdρy

(9.5.9)

 

 

0 Ω0

и среднего квадрата волновой аберрации:

145

 

 

2 =

1

∫∫

W 2

(ρx , ρy )dρxdρy

(9.5.10)

W

 

Ω

 

 

 

 

 

0 Ω

 

 

 

 

 

 

 

 

0

 

 

 

Тогда выражение (9.5.8) запишется в виде:

h(0)≈ Ω02 1+ 2πiW 2π 2W 2 2 (9.5.11)

Модуль комплексного числа z = a +ib вычисляется как сумма квадратов вещественной и мнимой частей z 2 = a2 +b2 , следовательно:

h(0)≈ Ω02 1

 

+ 4π 4 (W 2

)2 + 4π 2

 

 

= Ω02 1 4π 2 [W 2

(

 

)2 ] (9.5.12)

+ 2π

W 2

W 2

W

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Значение ФРТ в максимуме при отсутствии аберраций определяется выражением:

h (0)= Ω2

 

 

 

(9.5.13)

0

0

 

 

 

 

 

 

 

 

Тогда формула Марешаля:

 

St =

h(0)

 

14π 2 [W 2

 

(

 

)2 ]

(9.5.14)

 

W

h (0)

 

 

 

 

 

 

 

 

 

 

 

0

 

[W 2

 

 

)2 ]

 

 

 

 

Величина

(

 

называется дисперсией волновой аберрации по

W

зрачку (дисперсия – это разность среднего квадрата и квадрата среднего значения):

DW = (W W )2 = (W 2 2WW + (W )2 )=W 2 2WW + (W )2 =W 2 (W )2 (9.5.15)

Формула Марешаля показывает, что важна не сама волновая аберрация, а ее изменение (деформация волнового фронта) по зрачку.

Средний квадрат деформации волнового фронта – это квадратный корень из дисперсии:

DW =Wскв

(9.5.16)

Формула Марешаля дает возможность приблизительно оценить число Штреля, если известен средний квадрат деформации волнового фронта:

St 1 4π 2 D

=1 4π 2W

2

=1 (2πW

скв

)2

 

 

(9.5.17)

W

скв

 

 

 

 

 

Если St 0.8 14π 2 D

,

то, следовательно,

D

1

, а допуск

 

 

W

 

 

 

 

W

200

 

 

 

 

 

 

 

 

 

Марешаля на средний квадрат деформации волнового фронта:

Wскв

 

1

(9.5.18)

14

 

 

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

146

9.5.4. Влияние аберраций на ОПФ. Геометрически-ограниченные и дифракционно-ограниченные оптические системы

При наличии аберраций ОПФ оптической системы становится меньше, чем ОПФ безаберрационной системы. На графике ЧКХ можно показать, как аберрации влияют на форму кривой контраста (рис.9.5.4). Кривые ЧКХ в присутствии аберраций могут иметь сложную форму, но они никогда не превышают кривую безаберрационной ЧКХ.

k

1

дифр. ограниченная о.с.

геом.

безаберрационная о.с.

 

ограниченная о.с.

 

0.2

 

ω

0

0.5

1

2

 

Рис.9.5.4. Влияние аберраций на ЧКХ.

Дифракционно-ограниченные оптические системы имеют рабочий интервал частот, превышающий половину от предельной, то есть ω >1 (рис.9.5.4). Качество изображения в таких системах определяется в основном явлениями дифракции и непосредственно зависит от отношения апертуры к

длине волны Aλ . Остаточные аберрации должны оцениваться по критерию

Марешаля (9.5.18). К дифракционно-ограниченным системам относятся, в частности, измерительные системы, проекционные оптические системы для микроэлектроники и системы, работающие с глазом.

К геометрически-ограниченным относятся оптические системы, рабочий интервал частот для которых не превосходит ω = 0.5 в канонических частотах (рис.9.5.4). Качество изображения таких систем определяется картиной поперечных аберраций и непосредственно не зависит от длины волны и апертуры. Степень коррекции геометрически-ограниченных систем оценивается поперечными аберрациями. К таким системам относятся, в частности, кино-, фото-, и телевизионные объективы.

147

Приложения.

Приложение А. Дифференциальные операторы математической теории поля

Математический аппарат, применяемый в теории электромагнитного поля – это теория скалярного и векторного поля. Предмет этого раздела математики – скалярные (А.1) и векторные (А.2) функции от трех пространственных переменных (радиус-вектора r точки в пространстве):

U (x, y, z)

F(x, y, z)

=U (r)

Fx (x, y, z)

=Fy (x, y, z) = Fxi + Fy j + Fzk

Fz (x, y, z)

(А.1)

(А.2)

В соотношениях теории поля используются дифференциальные операторы

B

) и по пространственным

дифференцирования по времени (например, B =

t

 

 

координатам. Операторы дифференцирования по пространственным координатам могут быть векторами или скалярами, с которыми можно производить все известные из векторной алгебры действия, в частности, векторное произведение (А.3), скалярное произведение (А.4) и смешанное произведение (А.5):

×E

(А.3)

D

(А.4)

×( U)

(А.5)

Дифференциальные операторы 1-го порядка

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

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

=

y

=

x i +

y j +

z k

(А.6)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

Применяя оператор к скалярному или векторному полю, можно получить следующие скалярные и векторные величины:

U

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

U

 

=

U i +

U j +

U k

(А.7)

U = gradU =

 

 

U

y

 

x

y

z

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

148

 

F

 

Fy

 

 

F

 

 

 

F = divF =

 

 

 

x

 

+

 

 

 

+

 

z

 

 

(А.8)

 

x

 

y

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

j

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

×F = rotF =

 

x

 

y

 

 

(А.9)

 

 

 

 

 

z

 

 

 

 

F

 

 

 

F

y

 

 

F

 

 

 

 

 

 

 

x

 

 

 

 

 

 

z

 

 

Результаты выражений (А.7) и (А.9) – векторы, а результат выражения (А.8) – скаляр.

Дифференциальные операторы 2-го порядка

Оператор 2 называется оператором Лапласа (является скаляром):

2 = =

2

+

 

2

+

 

2

(А.10)

x2

y2

z2

 

 

 

 

Применение этого оператора к скалярному полю дает скалярную величину (А.11), а к векторному – векторную (А.12):

2

2U

 

2U

 

2U

U =

x2

+

y2

+

z2

2 Fx2F = 2 Fy2Fz

(А.11)

(А.12)

Основные математические тождества теории поля

( F)= grad (divF)

(А.13)

×( ×F)= rot(rotF)

(А.14)

( F)= 2F + ×( ×F)

(А.15)

×( U )= rot(gradU )= 0

(А.16)

( ×F)= div(rotF)= 0

(А.17)

Все рассмотренные соотношения широко используются в оптике для описания светового поля, вывода уравнений и законов геометрической оптики.

149

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]