Оптическая информатика. Конспект лекций [Электронный ресурс] (90
.pdf
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 |
|
¶x¶y ¶x¶z |
|
|
|
|
|
¶y¶x |
|
|
|
|
¶y |
|
|
|
|
|
¶y¶z |
|
|
|
|
|
|
¶z¶x |
|
|
¶z¶y |
|
|
|
¶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) |
|
|
|
|
ω |
|
2π |
λ0 |
|
|
|
|
|
|
|
|||
где |
k = ω με = |
υ |
= |
λ |
– волновое число, λ = n – длина волны в среде, λ0– |
длина |
||
волны в воздухе.
Фазовая скорость одной монохроматической гармонической волны:
υ = ω . |
(31) |
k |
|
Для квазимонохроматической волны
ω+ ω, ω >> ω ,
выражение для фазовой скорости совпадает с (31), но вводится также понятие групповой скорости:
|
υg = |
Dω . |
(32) |
||
|
|
Dk |
|
||
В пределе при ω→0 |
|
|
|
|
|
υg |
= lim |
Dω = |
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.
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) |
2πR 2 |
|
||||
|
|
R |
|
||
