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

Оптическая информатика. Конспект лекций [Электронный ресурс] (90

.pdf
Скачиваний:
9
Добавлен:
15.11.2022
Размер:
450.74 Кб
Скачать

11

Рис. Взаимное расположение векторов электрической (E) и магнитной (H) напряженности и направления распространения света (S).

S =E×H ,

(17)

S – вектор Умова-Пойнтинга, это вектор плотности потока энергии электромагнитного поля, по модулю равен количеству энергии, переносимой через единичную площадь, нормальную к S, в единицу времени, а своим направлением вектор определяет направление переноса энергии волны.

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

Возьмем уравнение для ротора

электрического поля Ñ´Ε = -

B и

 

 

 

 

t

материальное уравение B = μ H , тогда:

 

 

 

 

Ñ´Ε = -

B = -μ

H ,

 

 

t

 

t

 

Векторно домножим это уравнение на Ñ:

 

 

H

Ñ´Ñ´Ε = -μ

Ñ´

 

= -μ

 

t

 

 

t

 

 

2

 

 

2E

= -μ

 

 

(εE)

= -με

 

 

t

2

t

2

 

 

 

 

 

 

 

D

 

(Ñ´H)

= -μ

 

 

 

=

 

 

 

t

t

 

12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Оператор Лапласа D = Ñ2 = Ñ ×Ñ =

 

2

 

+

 

2

 

 

+

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

y2

 

 

z2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2U

2U

 

 

 

2U

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

DU = Ñ U =

x2

+ y2

+ z2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

DF = Ñ2F = (Ñ2 Fx , Ñ2 Fy , Ñ2 Fz ) = Ñ2 Fx i + Ñ2 Fy j + Ñ2 Fz k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

div gradU = Ñ ×(ÑU ) = Ñ2U

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

j

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rot gradU = Ñ´ (ÑU ) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

y

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U

 

U

 

 

U

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U

 

 

 

U

U

 

 

 

 

U

 

 

 

U

 

 

U

 

 

 

 

 

 

 

 

 

=

 

 

 

-

 

 

 

 

i +

 

 

 

 

 

 

 

 

-

 

 

 

 

 

 

 

 

j +

 

 

 

 

 

 

 

-

 

 

 

 

 

 

k = 0

 

 

 

 

 

 

y z

 

 

 

 

x z

 

z

 

 

 

 

 

y

x

x

 

 

 

 

 

 

 

 

 

 

 

z y

 

 

 

 

 

 

x

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

grad divF = Ñ ×(Ñ ×F) =

 

¶ ¶

 

 

 

 

 

 

 

F

 

 

 

Fy

 

 

 

 

F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

,

 

 

 

 

 

×

 

 

 

 

 

x

+

 

 

 

 

 

+

 

z

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

y

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

y

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 F

 

2 Fy

 

2 F

 

 

 

 

 

2 F

 

 

 

2 Fy

 

 

 

 

2 F

 

 

 

 

 

 

2 F

 

 

2 Fy

 

 

2 F

 

 

 

 

=

 

 

x +

 

 

 

 

 

+

z

i +

 

 

 

 

x

+

 

 

 

 

 

 

+

 

 

 

z

j +

 

 

 

 

 

x

+

 

 

 

 

+

 

 

z k =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

2

 

xy xz

 

 

 

 

 

yx

 

 

 

 

y

 

 

 

 

 

yz

 

 

 

 

 

 

zx

 

 

zy

 

 

 

z

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= Ñ2F + Ñ´ (Ñ´ F)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rot rotF = Ñ´(Ñ´ F) = Ñ ×(Ñ ×F) - Ñ2F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

div rotF = Ñ ×(Ñ´ F) =

 

F

 

 

 

 

Fy

 

 

 

 

 

 

 

F

 

 

 

 

 

F

 

 

Fy

 

 

 

F

 

 

 

 

 

 

 

 

 

 

 

z

-

 

 

 

 

 

 

+

 

 

 

 

 

 

x

-

 

 

z

+

 

 

 

 

 

-

 

x

= 0

 

 

 

 

 

y

 

 

 

 

 

 

 

y

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

x

z x

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из Ñ´Ñ´Ε = -με

2E

, используя Ñ´(Ñ´ F) = Ñ ×(Ñ ×F) - Ñ2F , получаем:

 

 

t 2

 

 

 

 

 

 

 

 

Ñ ×(Ñ ×Ε) - Ñ2Ε = -με

2E .

 

 

 

 

 

 

 

t 2

 

В неоднородной среде ε(x,y,z) волновое уравнение выглядит следующим

образом ( D = ε E ):

 

 

 

 

 

 

 

 

 

2

2E

 

Ñε

(18)

 

 

Ñ E - με

t

2

= -Ñ E ×

.

 

 

 

 

 

ε

 

13

Если же среда однородная, то градиент диэлектрической проницаемости равен нулю Ñe=0, и тогда правая часть тоже будет равна нулю ( Ñ ×D = 0 ):

 

Ñ2E - με

2E = 0 ,

(19)

 

 

t 2

 

и аналогично

Ñ2H - με

2H = 0 .

(20)

 

 

t2

 

Уравнения (19) и (20) являются векторными, имеющими по три скалярных компоненты каждый. Уравнение (19) можно представить в декартовых координатах следующим образом:

 

- με

2E

Ñ2Ex

t2x

 

 

2E

 

- με

Ñ2Ey

 

y

t

2

 

 

 

 

- με

2E

Ñ2Ez

t

2z

 

 

 

= 0

= 0

(21)

= 0

Т.е. каждая компонента электромагнитного поля, если обозначить ее через

ψ(x,y,z,t), удовлетворяет скалярному волновому уравнению:

 

 

 

 

 

 

 

 

 

 

 

 

2

 

1 2Y(x, y, z,t)

 

 

 

,

(22)

 

 

 

 

 

 

 

 

 

 

Ñ

 

 

Y(x, y, z,t) -

 

 

t 2

 

 

= 0

 

 

 

 

 

 

 

 

 

 

υ 2

 

 

где υ =

1

 

– скорость волны в среде.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

με

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Учитывая, что скорость волны в воздухе

c =

 

1

 

 

,

показатель преломления

 

 

 

 

 

 

 

 

 

 

μ

ε

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

данной среды n =

c

 

με

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

υ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

μ

ε

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для постоянных гармонических по времени –

монохроматических - полей

временная зависимость имеет определенный вид:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ψ(x, y, z,t) = ψ (x, y, z) exp(iωt + iϕ ) ,

(23)

где ω –

частота излучения, ϕ – начальная фаза.

 

 

 

 

 

 

 

В частности,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ex (x, y, z,t) Ey (x, y, z,t) Ez (x, y, z,t)

= Ex (x, y, z) exp(iωt + iϕ x ),

 

= Ey (x, y, z) exp(iωt + iϕ y ),

(24)

= Ez (x, y, z) exp(iωt + iϕ z ).

 

14

Тогда, беря производную по времени, уравнения Максвелла (13)-(16)

можно переписать:

 

Ñ´ E = -iωμH ,

(25)

Ñ ´ H = iωεE ,

(26)

Ñ × H = 0 ,

(27)

Ñ ×E = 0 ,

(28)

где E(x, y, z) = iEx (x, y, z) + j Ey (x, y, z) + kEz (x, y, z) и т. д.

 

Откуда можно получить (будет подробно рассмотрено в третьей лекции)

волновое уравнение Гельмгольца:

 

 

 

 

 

 

 

(Ñ2 + k 2 )E(x, y, z) = 0 ,

(29)

 

 

 

 

 

 

 

(Ñ2 + k 2 )H(x, y, z) = 0 ,

(30)

 

 

 

 

ω

 

λ0

 

 

 

 

 

 

 

где

k = ω με =

υ

=

λ

– волновое число, λ = n – длина волны в среде, λ0

длина

волны в воздухе.

Фазовая скорость одной монохроматической гармонической волны:

υ = ω .

(31)

k

 

Для квазимонохроматической волны

ω+ ω, ω >> ω ,

выражение для фазовой скорости совпадает с (31), но вводится также понятие групповой скорости:

 

υg =

Dω .

(32)

 

 

Dk

 

В пределе при ω→0

 

 

 

 

 

υg

= lim

Dω =

.

(33)

 

 

ω →0

Dk d k

 

Фактически, групповая скорость представляет собой фазовую скорость огибающей.

15

Волновое уравнение. Частные решения волнового уравнения. Интегральные преобразования в оптике.

В волновой теории света предполагается, что свет распространяется в виде волны. Оптическая волна описывается скалярной комплексной волновой функцией пространственных (x,y,z) и временной t координат:

V(x,y,z,t)

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

Распространение светового поля в пространстве описывается волновым

уравнением

 

 

1

2

 

 

 

 

 

 

 

 

Ñ2 -

 

 

V (x, y, z, t) = 0

(1.1)

c2

 

t 2

 

 

 

 

 

 

 

 

 

где c - скорость света в свободном пространстве,

 

 

Ñ2 =

2

+

2

+

2

 

 

 

x 2

y2

z2

 

 

- оператор Лапласа или лапласиан.

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

В свободном пространстве свет распространяется с одинаковой скоростью c0 = 3×108 м/сек. При этом свет может иметь различную длину волны l.

Оптические длины волн делятся на 3 группы: инфракрасные, видимые и ультрафиолетовые. Длина волны видимого диапазона 390-760 nm, т.е. ~ 0.5

мкм (5×10−7 м).

Для каждой длины волны l определена частота электромагнитных колебаний n, удовлетворяющая следующему соотношению nl=c.

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

exp[ij(x, y, z)] = arg[U (x, y, z)]

16

Источником такого монохроматического света является, например, лазер. Лазеры бывают разных цветов: синие, зеленые, красные. Наиболее распространен He-Ne лазер с длиной волны 633 nm (красный).

Для такой монохроматической волны временная зависимость может быть выражена в явном виде

V (x, y, z, t) = U (x, y, z) exp(−i2πνt) .

(1.2)

Подставим (1.2) в (1.1):

 

 

 

1

2

 

[U (x,

 

Ñ2

-

 

 

U (x, y, z) exp(-i2pnt) = exp(-i2pnt)Ñ2

2

 

t

2

 

 

 

c

 

 

 

 

 

1

 

 

 

 

 

 

2

 

 

= exp(-i2pnt) Ñ2

[U (x, y, z)]-U (x, y, z)

 

(- i2pn)

 

= 0

 

2

 

 

 

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ν

 

1

 

 

 

 

 

2pn

2

 

=

 

 

 

 

 

c

λ

 

 

Ñ2[U (x, y, z)]+U (x, y, z)

 

 

 

= Ñ2

+

 

 

 

 

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y, z)]-U (x, y, z)

2p 2

U (x, y, z)

l

1 2 [exp(-i2pnt)] = c2 t 2

= 0

обозначив k = 2π/λ - волновое число света (для оптических волн ~ 107 м),

получим уравнение Гельмгольца:

[Ñ2 + k 2 ]U (x, y, z) = 0

(1.3)

описывающее распространение пространственной части.

U (x, y, z) = A(x, y, z) exp[ij(x, y, z)]

называется комплексной амплитудой. При этом величина A(x, y, z) = U (x, y, z) называется амплитудой волны, фазой.

С комплексной амплитудой U(x,y,z) можно связать конкретную физическую величину, которая измеряется оптическими приборами. Такой величиной является оптическая интенсивность света:

I (x, y, z) = U (x, y, z) 2 =U (x, y, z)U (x, y, z)

Интенсивность монохроматической волны не зависит от времени. Волновым фронтом называются поверхности равной фазы:

ϕ(x, y, z) = const = 2πq , q - целое.

Нормаль волнового фронта в точке (x,y,z) параллельна вектору градиента

фазы

17

¶j(x, y, z) Ñj(x, y, z) =

x

,

¶j(x, y, z) ¶j(x, y, z)

y

,

z

,

 

 

 

который указывает направление наибольшей скорости изменения фазы.

Частные решения уравнения Гельмгольца: плоская и сферическая волны

Плоская волна

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

Плоская волна имеет следующую комплексную амплитуду:

 

 

F ( x, y, z) = Aexp(ikx) = Aexp[i(kx x + k y y + kz z)].

(1.4)

где A

-

комплексная константа,

называемая

комплексной

оболочкой,

k = (kx , k y , kz ) - волновой вектор.

 

 

 

 

 

Так

как фаза arg{F (x)}= arg{A}+ kx ,

 

то волновой фронт определяется из

уравнения

{

}

,

q - целое,

и будет представлять собой

kx = kx x + k y y + kz z = 2pq - arg A

 

параллельные плоскости, перпендикулярные волновому вектору k (поэтому и называется плоская волна). Расстояние между двумя плоскостями:

{

A

}

и

{

}

kx = kx x + k y y + kz z = 2pq - arg

 

kx'= kx x '+k y y +' kz z =' 2p q(+1 -) argA

 

 

k(x - x' )

 

= 2p

 

 

(x - x '

)

=

2p

= l - длина волны

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

Таким образом, волновая функция (1.4) является периодической в пространстве с периодом λ. Интенсивность плоской волны I (x) = A 2 = const везде в пространстве. Ясно, что эта волна является абстракцией.

Подставим (1.4) в (1.3):

(Ñ2 + k 2 )A exp[ik (ax + by + gz)] = A exp[ik (ax + by + gz)]{(i)2 k 2 (a2 + b2 + g 2 ) + k 2} = 0

равенство выполняется при условии a2 + b2 + g 2 =1 g = ±1 - a2 - b2 .

Следовательно, плоская волна представима в виде:

 

 

 

 

 

 

 

 

 

18

 

F ( x, y, z) = Aexp[ik (αx + βy ± z

 

 

)]= F ( x, y, z = 0)Φ(α, β , z) = F ( x, y)Φ(α, β , z) =

 

1 − α 2 − β 2

 

= Aexp[ik (αx + βy)]exp[ik ± z

 

 

]= F

 

 

 

0

 

 

 

 

 

 

 

1 − α 2 − β 2

( x, y, z)

 

 

 

 

 

 

α ,β

 

 

 

 

где

Φ(α, β , z) = Aexp(± ikz

 

)

(1.5)

1 − α 2 − β 2

- уравнение распространения плоских волн в пространстве вдоль оси z.

Знак "+" перед квадратным корнем выбирается в случае, если волна приходит к плоскости z слева, знак "-" , если волна приходит справа. В дальнейшем будем выбирать знак "+".

Из (1.5) следует условие на параметры α и β:

α2 + β2 < 1

(1.6)

В случае, если α2 + β2 > 1 , то вместо (1.5) будем иметь

Φ(α, β , z) = exp(kz

 

)

(1.7)

α 2 + β 2 − 1

Амплитуда (1.7) описывает затухающие по экспоненте световые волны, которые не распространяются далее в пространстве по оси z, и мы исключим их из рассмотрения.

Базис плоских волн

Наша цель теперь состоит в получении связи между некоторой амплитудой света на плоскости z: U(x,y,z) и амплитудой света на плоскости z=0: U0(x,y) [1].

y

y

U0(x,y)

U(x,y,z)

 

x

x

z=0

z

 

Рис. 4

Представим функцию U(x,y,z) в плоскости z в виде разложения по базису плоских волн:

 

 

 

19

 

 

U ( x, y, z) = ∫ ∫ C(α, β )Fα ,β ( x, y, z) dα d β =

 

 

−∞

.

(1.8)

 

 

= ∫ ∫ C(α, β ) exp[ik(αx + βy)]exp[ikz

 

]dα d β

 

1 − α 2 − β 2

 

−∞

Такое разложение выглядит как результат обратного преобразования Фурье. Действительно, физически, преобразование Фурье - это разложение поля по плоским волнам, а переменные α и β, пропорциональные углам Эйлера, задающим направление распространения плоской волны, являются частотными переменными, С(α,β) - веса соответствующих плоских волн.

При z=0 имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U ( x, y, z = 0) = U0 ( x, y) = ∫ ∫ C(α, β ) exp[ik (αx + βy)]dα d β .

 

(1.9)

 

 

 

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

Следовательно

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C(α, β ) = ∫ ∫U0 ( x, y) exp[ik(αx + βy)]d x d y .

 

(1.10)

 

 

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

Далее, подставив (1.10) в (1.9) , получим

 

 

 

 

 

U(x, y, z) =

∞ ∞

U (x',y ')exp− ik(αx +'βy )' dx

'yd '

expikz

 

 

]

expik(αx + βy)

]

αd βd=

1− α2

− β2

 

 

 

∫ ∫ ∫ ∫

0

[

 

 

 

]

 

 

[

 

 

[

 

 

 

 

−∞ −∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

[

 

 

 

 

]

exp(ik[α (x

x ' )+ β y( − y ])' dα

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

∫ ∫

 

∫ ∫

 

 

 

 

βd

 

dx 'yd ='

 

(1.11)

 

U (x',y ')

 

 

 

2

2

 

 

 

 

 

 

expikz 1− α −β

 

 

 

 

 

 

−∞

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= ∫ ∫U0(x',y ')H (x x ',y y',z )dx 'dy '

 

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

[

 

 

]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∫ ∫

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H (x x ',y y ',z )=

 

 

 

 

 

 

 

exp(ik[α (x x ' )+ β y( − y '])

 

dα dβ

где

−∞

 

expikz

1− α2 − β2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Интеграл (1.11) является конечной целью. Он связывает комплексную

амплитуду

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U0(x,y)

в плоскости

z=0

 

и комплексную амплитуду

 

U(x,y,z) в плоскости z

(рис. 4).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

20

Сферическая волна

Рассмотрим комплексную амплитуду сферической волны

 

exp(ik

 

 

)

 

 

 

 

x 2 + y2 + z2

R =

 

 

 

 

x

2 + y2

+z 2

S (x, y, z) = A

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

x 2

+ y2 + z2

 

 

 

 

 

 

 

 

 

A

exp(ikR )

.

(2.1)

 

 

R

 

которая также является частным решением уравнения Гельмгольца (1.3), (что мы покажем на практике), R = x 2 + y2 + z2 - расстояние от начала координат до текущей точки. Волновой фронт представляет собой поверхности постоянной фазы kR = 2πq − arg{A} , считая для простоты, что arg{A} = 0 , можно записать:

R = λq , q - целое.

Это уравнение концентрических сфер, разделенных по радиальному

расстоянию на величину λ (длина волны). Интенсивность I ( x, y, z) =

 

 

A

 

2

-

 

 

 

 

 

 

R2

 

 

 

 

обратно пропорциональна квадрату расстояния.

Разложение сферической волны по плоским волнам имеет вид (будет показано на практике)

 

 

exp(ikR)

 

 

1

 

 

 

 

 

 

 

 

 

( x, y, z) dα d β =

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

−∞∫ ∫ ik 1 − α 2

β 2

α ,β

 

 

 

 

 

 

 

.

(2.8)

 

 

 

exp(ikz

 

 

 

 

 

 

 

 

)

 

 

 

 

 

1 − α 2 − β 2

exp[ik (αx + βy)]dα d β

 

 

= −2π ∫ ∫

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ik 1 − α

2

β

2

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

Ранее мы

получили

 

связь комплексной амплитуды поля U0(x,y)

в

плоскости z=0

с комплексной амплитудой U(x,y,z) в плоскости z (рис. 4):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U(x, y, z) = ∫ ∫U0(x',y ')H (x x ',y y ',z )dx 'yd '

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H (x, y, z) = ∫ ∫ exp[ikz

 

 

]exp[ik (αx + βy)]d α d β

 

где

1 − α2 − β2

(2.9)

−∞

Учитывая связь правых частей формул (2.8) и (2.9), можно показать (будет на практике), что

 

z exp(ikR )

1

 

H (x, y, z) = −

 

ik

 

 

(2.10)

R 2

 

 

 

R

 

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