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

m1012

.pdf
Скачиваний:
1
Добавлен:
15.11.2022
Размер:
10.66 Mб
Скачать

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

Тогда из уравнений Максвелла для волн, распространяющихся в свободном пространстве, получают волновое уравнение [16]:

 

 

2

E(x, y, z,t)

 

2

E(x, y, z,t)

1

.

 

c

t

 

 

 

2

 

Решением этого уравнения для монохроматического будет синусоидальное скалярное поле:

E(x, y, z) a(x, y, z)cos[ t (x, y, z)],

где = 2 f, f – частота света; или

E(x, y, z,t) Re[a(x, y, z)e

i t

],

 

(2.1)

света

(2.2)

(2.3)

где a(x, y, z) – комплексная амплитуду, так и фазу волны:

амплитуда (фазор) определяет как

a(x, y, z) a(x, y, z)e

i ( x, y, z )

.

(2.4)

 

Комплексная амплитуда a(x, y, z)

содержит всю информа-

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

Процесс регистрации световых волн, их распространение через однородное пространство, оптические элементы, голограмму и т.д. будем рассматривать в рамках теории линейных, простран- ственно-инвариантных оптических систем [28, 47, 67, 79] не в координатной, а в пространственно-частотной области.

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

Однако операция, связывающая координатную и частотные области, – преобразование Фурье – отражает физическую сущность действия оптических систем и позволяет (в отличие от рассмотрения в координатной области) с единой точки зрения рас-

31

смотреть методы голографической интерферометрии, муара, спекл-интерферометрии.

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

где функция

 

 

 

 

 

 

 

 

 

 

a(x, y)

 

A( fx

, f y )e

i 2 ( f

x f

y)

dfxdf y

,

(2.5)

 

x

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A( fx , f y ) называется фурье-образом функции a(x, y) ,

или ее спектром. Последнее выражение

форме

F a(x, y) A( fx , f y )

чаще представляют в операторной

. Таким образом, сложное распреде-

ление амплитуд в некоторой плоскости может быть представлено в виде суммы плоских волн, распространяющихся в разных направлениях и имеющих разные амплитуды [47, 67]. Величины fx , f y называются пространственными частотами и позволяют

определить направления распространения этих плоских волн по направляющим косинусам:

K

 

f

,

K

 

f

 

,

K

 

f

,

( f

 

 

1 ( f

2

f

2

) ).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

x

 

 

y

 

y

 

 

z

z

 

 

z

 

 

x

 

y

 

(2.6)

Регистрируемая в плоскости голограммы выходная функция b(x, y) , соответствующая входной a(x, y), имеет спектр B( fx , f y ), определяемый произведением спектра входной функции на ча- стотно-передаточную функцию оптической системы H ( fx , f y ) :

B( f

, f

y

) A( f

, f

y

) H ( f

, f

y

),

x

 

x

 

x

 

 

(2.7)

а сама функция

b(x, y

ции B( fx , f y ) :

 

b (x, y)

A( fx

)

является обратным фурье-образом функ-

F 1 B( fx , f y )

(2.8)

, f y )H ( fx , f y )e i 2 ( fx x f y y )dfxdf y ,

т.е. суммой всех ее фурье-компонент.

32

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

ивосстановления одной из фурье-компонент (плоской волны). Результат для произвольной волны на основании принципа суперпозиции получают простым суммированием. С этих позиций

ирассмотрим основные этапы проведения экспериментов, перечисленные выше.

2.3. Регистрация информации

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

энергии и пропорционально квадрату амплитуды

a

2

(x, y, z).

 

Квадрат этой амплитуды обозначают через

I (x, y, z) a

2

(x, y, z)

 

и называют интенсивностью.

Степень воздействия света на фотодетектор пропорциональ-

на экспозиции Е:

2

( – время экспозиции) и не зависит

E I a

от фазы. Поэтому все известные фотодетекторы, включая и глаз человека, называют квадратичными.

Для того чтобы сведения о фазе сигнала не были потеряны, они кодируются в голографии в виде амплитудных изменений на фотодетекторе при интерференции объектной волны (волновой

вектор которой K0 ) с когерентной ей опорной – волновой вектор

KR (рис. 2.1, а).

aR

В качестве опорной волны выберем плоскую с амплитудойconst (не зависящей от координат x, y, z). Плоская волна,

произвольно ориентированная относительно осей координат, имеет фазу, определяемую скалярным произведением

k r k

x k

y

y k

z,

x

 

z

 

где направление волнового вектора

k

сов-

падает с

направлением распространения плоской волны,

а

k 2 Пусть волна распространяется в плоскости yoz, а

k

составляет

угол

R

с осью z (см.

рис. 2.1, а). В этом случае

k r ky y kz z,

где

ky

k sin R ,

kz

k cos R , а пространствен-

ные частоты равны:

fy

ky / 2

fz

kz / 2 Тогда в плоскости

хоy амплитуда опорной волны имеет вид:

 

 

 

 

 

 

 

 

33

 

 

 

 

 

 

 

a(x, y) a e

i(k

y

)

,

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

где R – фаза волны на прямой y = 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

K

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а)

 

 

 

 

 

 

б)

 

 

 

0 + 2 R

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K

0

 

 

 

 

 

 

 

 

 

K

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

K

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

K

R

H

 

 

 

 

H

0

 

K

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 2.1. Общая схема:

а – запись голограммы Н; б – восстановление изображения

(2.9)

Если объектная волна также плоская с волновым вектором,

имеющим проекции: k0

k sin

и k

0

k cos (см. рис. 2.1, а),

 

 

 

y

 

 

 

 

0

 

 

 

 

z

 

 

 

 

0

 

 

то амплитуда этой волны в плоскости xoy имеет вид:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

y

)

 

 

 

 

 

 

 

 

a0 (x, y) a0e

i(k

.

 

 

 

(2.10)

 

 

 

 

y

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Интенсивность суммарной волны I(x, y) равна:

 

 

 

I (x, y) (a a

R

)(a a

R

)*

 

 

 

 

 

 

 

 

 

 

 

0

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

a2 a2

a a

 

ei[(ky ky0 ) y ( R 0 )] a a

e i[(ky ky0 ) y ( R 0 )]

 

0

R

0

 

R

 

 

 

 

 

 

 

 

 

 

 

0 R

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I (x, y) a2

a2

2a a

 

cos[

2

(sin

 

sin

) y (

 

)] , (2.11)

R

 

k

R

0

R

 

0

 

 

 

 

 

 

 

 

 

 

 

0

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т.е. представляет собой периодическую структуру с частотой f:

 

 

 

f

(sin R sin 0 ) / .

 

 

 

(2.12)

При a0 = aR картина имеет наибольший контраст:

 

 

 

 

I (x, y) 2a2[1 cos(2 f y )],

 

 

(2.13)

где = R 0 – начальная фаза.

 

 

 

 

 

 

 

 

 

 

 

 

34

Типичная для галоидосеребряных сред кривая зависимости амплитудного пропускания обработанной фотоэмульсии от экспозиции на среднем участке может быть аппроксимирована прямой [42]. Если регистрации интерференционной картины соответствует линейный участок (линейная регистрация), то амплитудное пропускание фотопластинки будет линейно связано с интенсивностью I(x, y). Тогда результат интерференции двух плоских волн a0

и aR после фотообработки можно представить в виде транспаранта, функция амплитудного пропускания которого равна t:

t t

0

t

cos(2 f

y

).

 

1

 

 

(2.14)

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

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

2.4. Восстановление объектной волны

Если процесс записи голограммы основан на интерференции объектной и опорной волны, то ее восстановление – это дифракция опорной волны на записанном транспаранте. Направим плоскую опорную волну на записанный транспарант под тем же углом ( ky k sin R ), что и при записи (рис. 2.1, б):

aR

a e

i(k

y k

z )

y

z

 

 

 

R

 

 

 

.

(2.15)

Считая, что голограмма расположена в плоскости z = 0, найдем распределение комплексных амплитуд непосредственно за голограммой с амплитудным пропусканием (2.14):

a

R

eiky yt a

t

eiky y a

t cos(2 f

y

y )eiky y

 

 

 

 

 

R 0

 

 

 

 

R 1

 

 

a

t

eiky y

1

a

t eiky y [ei (k p y ) e i (k p y ) ]

(2.16)

 

 

 

 

R 0

 

2

 

R 1

 

 

 

 

 

a

t

eiky y a

 

t ei ei(ky k p ) y a

 

t e i ei(ky k p ) y ,

 

 

 

 

R 0

 

 

 

 

R 1

 

 

R 1

 

где kp k(sin R sin 0 ) .

35

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

k

y

k sin

0

,

 

 

 

 

 

k

y

k

p

k sin

 

 

 

 

 

k

y

k

p

k sin

 

 

 

 

 

, .

(2.17)

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

sin

R

sin

0

 

 

 

 

 

 

sin

R

(sin

R

sin

) sin

 

 

0

 

sin

R

(sin

R

sin

) sin

 

 

0

 

0 R ,

0 ,

0 2 R .

(2.18)

Одна из волн – продолжение падающей aR, вторая – с волно-

вым вектором

k

0

– волна, воспроизводящая предметную, кроме

 

того, при ky

kp

k

образуется и третья волна с волновым векто-

 

 

 

 

 

 

ром

k .

При наличии более сложного распределения на транспа-

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

2.5. Дифракция плоских волн на периодических и квазипериодических структурах

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

Если пропускание транспаранта по оси y носит синусоидальный характер (2.14), то, положив в (2.16) ky 0 ( R 0), получим

a

(x, y,0) at

 

 

1

at e

i 2 f

y

 

1

at e

i 2 f

y

,

 

 

y

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

0

 

2

1

 

 

 

2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.19)

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

36

щая, и не несет никакой информации о пространственном распределении пропускания транспаранта, а две другие распростра-

няются симметрично относительно нее под углами

0

(

0

 

 

 

 

arcsin( f

y

)).

 

 

Вслучае произвольного пропускания его можно представить

ввиде суммы фурье-компонент. Тогда прошедшая такой транспарант волна имеет вид:

0

(x, y,0) at

0

 

 

1

1

 

i 2

f

 

x

 

 

1

1

i 2

f

 

x

 

a

 

 

 

 

 

 

at e

kx

 

 

at e

kx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

2

 

 

 

 

 

 

 

k

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

i 2

f

 

y

 

 

1

 

1

i 2

f

 

y

 

 

 

 

 

 

 

 

 

 

 

 

at e

 

 

ky

 

 

 

 

 

 

at e

 

 

 

ky

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

2

 

 

 

 

 

 

 

 

k

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.20)

т.е. является суммой плоских волн, имеющих пространственные частоты, кратные пространственным частотам транспаранта fx , f y .

Спектр функции a(x, y) можно наблюдать в задней фокальной плоскости линзы, формирующей действительное изображение транспаранта (рис. 2.2). Этот спектр является в общем случае непрерывной функцией ( fx , f y ); если пропускание транспаран-

та – периодическая функция x, y, то спектр дискретный.

Рис. 2.2. Схема наблюдения изображения при изменении спектра сигнала: G – объект; Л – линза; Н – изображение

Различные фурье-компоненты, прошедшие через линзу, складываясь, дают в плоскости изображения копию решетки. Помещая в фокальную плоскость различные диафрагмы, можно, удаляя отдельные спектральные компоненты, сознательно воздействовать на регистрируемый спектр изображения [28].

37

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

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

а)

б)

в)

г)

д)

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

38

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

2.6. Получение интерферограмм

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

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

стота исследуемого объекта

f (x, y)

изменилась на достаточно

малую величину f (x, y), a будет иметь вид:

то амплитуда восстановленной волны

a ~ a e

i 2 f y y

a e

i 2 ( f y f ) y

a (1 e

i 2 fy

)e

i 2 f y y

.

 

 

 

 

0

 

0

 

0

 

 

 

 

(2.21)

Поскольку f f y , то последнее выражение можно интер-

претировать как

восстановленную волну, амплитуда которой

промодулирована с частотой f . Это и есть интерферограмма.

Определение

деформированного состояния исследуемого

объекта на основе связи модуляции амплитуды

a(x, y)

с переме-

щениями U(x, y), V(x, y), W(x, y) или с деформациями (x, y) является основной задачей при расшифровке интерферограмм в методах голографической интерферометрии.

Сама связь между fx , f y и деформациями x , y может быть получена на основании теорем о преобразовании Фурье.

39

mx (x, y),

Пусть мации по функции a

исследуемый объект испытывает однородные дефоросям x и y (причем достаточно малые) так, что вид (x, y) не меняется, а изменяется только ее масштаб:

a(x, y) a(m x, m

y).

x

y

 

(2.22)

Если

F ( f

, f

y

)

x

 

 

– фурье-образ функции

a(x, y),

то в соответ-

ствии с теоремой функции a(mx x, m

y

масштабов [28] для преобразования фурье- y) получим

F

( f

, f

 

)

1

 

F( f

 

/ m

, f

 

y

m m

 

x

y

m

x

 

 

y

 

x

 

 

 

 

 

 

x

 

 

 

 

 

/

m

y

 

)

.

(2.23)

Если интерферограмма зарегистрировала оба состояния объекта, имеющего амплитуды a(x, y) и a(mx x, my y), фурье-образы

которых равны соответственно

F

( f

, f

y

)

m

x

 

 

и

1

F( f

 

/ m

, f

 

/ m

),

m m

x

y

y

x

 

y

 

x

 

 

 

 

 

 

то при восстановлении в соответствии с результатами, найденными ранее при наличии фильтрации (п. 2.5), получим восстановленные волны с амплитудами, промодулированными соответственно с частотами fx , f y , равными разности пространствен-

ных частот этих состояний:

f

 

f

 

 

f

x

 

f

(1

1

 

) f

 

 

m 1

f

 

 

 

,

 

 

 

 

 

 

 

 

x

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

x

 

m

 

x

 

 

m

 

 

x

m

 

x

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

x

 

 

 

 

 

 

x

 

 

 

x

 

 

 

 

 

 

f

 

f

 

(1

1

) f

 

 

y

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

m

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.24)

(2.25)

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

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

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

my (x, y) уже нельзя считать постоянными, а следовательно, непостоянны и частоты fx (x, y), f y (x, y).

40

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