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

ЭМ волны в материальных средах

.pdf
Скачиваний:
27
Добавлен:
23.02.2015
Размер:
579.82 Кб
Скачать

k , E (ω,k ) = ωcμ H (ω,k );

(2.2а)

(k ,εE (ω, k ))= i4πρst (ω,k );

(2.3а)

(k , μH (ω, k ))= 0 .

(2.4а)

2.3. Проектирование уравнений Максвелла на продольное и перпендикулярное на-

правления относительно вектора k

2.3.1. Общий вид уравнений

Спроектируем уравнения Максвелла на плоскость, перпендикулярную вектору k и на направление, параллельное ему. Все векторы представим в виде суммы двух проекций на эти направления, например:

E = E + E ; H = H + H

(2.18)

и т. д., где E k , E k .

Тогда уравнения (2.1а)-(2.4а) преобразуются к виду

k , H

 

= ωε

E

 

+i

4π

j st

ωε

E

+i

4π

j st ;

(2.19)

 

 

c

 

 

 

 

c

 

c

 

 

c

 

 

k , E

= ωμ H

 

+ ωμ H ;

 

 

 

 

(2.20)

 

 

 

c

 

 

 

 

c

 

 

 

 

 

 

(k ,εE )= i4πρst ;

 

 

 

 

 

 

(2.21)

 

(k , μH )= 0 ,

 

 

 

 

 

 

 

 

(2.22)

Каждое из уравнений (2.19)-(2.20) разобьем на поперечную и продольную части. Уравнение (2.19) распадается на два таких уравнения:

k , H

 

= −ωε

E

+i 4π

j st ;

(2.23)

 

 

c

 

c

 

 

 

 

 

 

 

ωε

E +i

4π

j st = 0 ,

 

(2.24)

 

c

 

c

 

 

 

или, что то же самое,

ωεE = i4π j st .

(2.24а)

22

E , посколь-

2.3.2. Анализ проекций уравнений Максвелла на поперечное и продольное направления

Анализ уравнений (2.19)-(2.24) позволяет сделать такие физические выводы. Век-

тор H всегда перпендикулярен вектору k : H = H , поскольку из уравнения (2.22) вид-

но, что H 0 . В соответствии с тем, что H 0 , уравнение (2.20) принимает вид

k , E

 

= ωμ H

 

,

(2.25)

 

 

c

 

 

 

 

 

 

 

откуда видно, что векторы k , E и H образуют правую тройку взаимно ортогональ-

ных векторов.

Уравнение (2.24а) представляет собой особый интерес. Выясним на конкретном примере, что следует понимать под j st и E , и как абсолютное значение этих векторов зависит от направления в пространстве. Направление в пространстве задается направлени-

ем вектора k , которое первоначально удобно отсчитывать от направления тока, возбуж-

дающего поле. Если речь идет о j st , надо понимать, что k j st , значит, электромагнит-

ные процессы рассматриваются в направлении, являющемся продолжением направления прямолинейного тока. В случае, когда говорят о j st , это значит, что интересующее на-

правление перпендикулярно направлению тока, и там вектор k перпендикулярен вектору тока: k j st . Направление векторов поля E , H , E , H в любом случае рассматри-

вается по отношению к вектору k .

Для простоты рассмотрим поле электрического диполя Герца. В направлении, сов-

падающем с направлением оси диполя, j st принимает максимальное значение, а j st = 0 .

Из (2.24а) следует, что и E в этом направлении максимально, а E = 0 . В направлении,

перпендикулярном оси тока, проекция тока j st отсутствует ( j st =0). Из того же равенства

(2.24а) следует, что здесь и E = 0 . Отлично от нуля только E . Пусть отсчет угла ϑ сферической системы координат начинается от направления линейного тока. Тогда с уве-

личением угла ϑ от нуля до π2 величина проекции E изменяется от максимального

значения до нуля, а проекция E , напротив, увеличивается от нуля до максимального значения. Это полностью коррелирует с представлением о диаграмме направленности электрического диполя Герца. В дальней зоне рассматривается только поле

ку с увеличением расстояния от точки наблюдения до источника поле E убывает гораздо быстрее, чем E .

23

Особый интерес представляет собой случай отсутствия сторонних источников. При отсутствии сторонних токов и зарядов имеет место равенство

ωεE = 0 .

(2.26)

В плазме, где под действием поля E происходит направленное движение электро-

нов, возможна ситуация, когда ε = 0 . Тогда E может быть отличным от нуля. Это так

называемые электростатические волны, в процессе распространения которых наблюдают-

ся не поперечные, а продольные колебания вектора электрического поля E .

2.3.3. Выражение векторных спектральных функций E (ω,k ) и H (ω,k ) через спектральную функцию j st (ω,k )

Напомним, что здесь везде идет речь о векторных спектральных функциях. Из совместного рассмотрения уравнений (2.23) и (2.25) можно выразить спектральные функции

E (ω,k ) и H (ω,k ) через спектральную функцию

j st (ω,k ). Исключим вектор H

из уравнения (2.23), подставив сюда его значение из (2.25). Получаем:

 

 

c

k k , E

 

= −ωε E

 

+

4πi

j st .

 

 

 

 

c

 

 

ωμ

 

c

 

 

 

 

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

 

k (k , E )k2 E

= −ω2εμ2

E

+ 4πiωμ2 j st .

(2.27)

 

 

 

 

c

 

 

 

c

 

Скалярное произведение в левой части равенства обращается в ноль. Из (2.27) получаем вектор E (ω,k ) в следующем виде:

E (ω,k )=

4πi

ωμ j st

 

 

 

2

 

 

 

.

(2.28)

c

 

2

 

 

 

 

ω εμ

k2

 

 

 

 

 

 

c2

 

 

 

 

 

 

 

 

 

Пользуясь соотношением (2.25), получим также и выражение для H (ω,k ):

H (ω, k )=

c

 

 

 

k , j st

 

 

 

 

 

 

= 4πic

 

 

 

 

 

 

 

2

2

 

2 .

(2.29)

ωμ

k , E

c

 

 

 

 

ω εμ k

 

 

 

24

Полученные выражения (2.28) и (2.29) представляют собой спектральные функции разло-

жения в интеграл Фурье векторных функций E (r,t ) и H (r ,t ). Этот интеграл для од-

ной из них (например, E (r,t )) с учетом (2.6), (2.14) и (2.28) принимает следующий вид:

∞ ∞

E (r,t )= 4πi ∫ ∫

−∞ −∞

 

 

 

ωμ j st

 

 

e

i(ωtkr )

dωdk .

(2.30)

με

(

ω

ω2

(kc)2

 

 

 

 

 

 

 

)

 

εμ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подчеркнем, что в подынтегральном выражении переменные интегрирования ω и k являются независимыми. Чтобы получить E (r,t ) в замкнутой аналитической форме,

выражение (2.30) необходимо дважды проинтегрировать – сначала по переменной ω, а

потом по k или наоборот. Проведем первое интегрирование по параметру ω.

2.4. Первое интегрирование в процессе определения векторных функций E (r,t ) и

H (r ,t )

Запишем внутренний интеграл J в таком виде

J =

 

ω j st (ω, k )

 

 

e

i(ωtkr )

dω .

(2.31)

 

 

 

 

 

 

 

 

 

 

kc

 

 

kc

 

 

 

−∞

 

 

 

 

 

 

 

 

ε (ω)

ω

 

 

ω +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

εμ

 

εμ

 

 

 

 

Интегрирование по действительной оси в (2.31) можно заменить интегрированием в комплексной плоскости ω. Подынтегральное выражение удовлетворяет условиям лем-

мы Жордана, которая при условии t > 0 позволяет замкнуть контур интегрирования по действительной оси дугой бесконечного радиуса в верхней полуплоскости. При этом контур интегрирования становится замкнутым. Известно, что интеграл по замкнутому контуру в комплексной плоскости равен сумме вычетов подынтегральной функции.

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

точках ω01 и ω02 комплексной плоскости ω

:

 

 

 

 

ω01 =

ck

 

; ω02 = −

ck

.

(2.32)

εμ

 

 

 

εμ

 

Выясним, куда попадут эти полюсы при наличии потерь в диэлектрике. Ранее было показано, что комплексная диэлектрическая проницаемость имеет вид ε =ε′+iε′′, где

25

 

 

 

′′

 

4πσ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ε

=ε , а

ε

= − ω < 0 . Ее можно представить и в другой форме:

ε =

ε

 

e

iδ

, где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δ = arctg

4πσ

. Это значит, что при наличии потерь полюс ω01

=

 

ck

=

 

 

 

ck

 

ei(δ 2)

ωε

εμ

 

 

 

 

 

 

 

 

 

 

 

 

 

ck

 

 

ck

 

 

 

 

 

 

ε

 

 

μ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

находится в верхней полуплоскости, а полюс ω02 = −

= −

 

 

ei(δ 2)

– в нижней.

 

 

 

 

 

 

 

 

 

 

 

 

εμ

ε

μ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

J = 2πi

ω(k ) j st (ω(k ), k )

e

i(ω(k )tkr )

=

ε (ω(k ))2ω(k )

 

 

=πi

j st (ω(k ),k )

e

i(ω(k )tkr )

.

(2.33)

ε (ω(k ))

 

 

 

Далее мы везде будем подразумевать, что ω =ω(k ), и писать просто ω. Подстав-

ляя (2.33) в (2.30), получаем после первого интегрирования следующее выражение для

E (r,t ):

E (r,t )= −4π2

j st (ω,k )

ei(ωtkr )dk .

(2.34)

 

−∞

ε (ω)

 

Аналогичным образом находим H (r,t ):

k , j st (ω, k )

H (r ,t )= −4π2c

 

ei(ωtkr )dk .

−∞ ωε (ω)μ(ω)

Еще раз подчеркнем, что в выражениях (2.34) и (2.35) параметры ω и ляются независимыми, а связаны соотношением

(2.35)

k уже не яв-

k =

ω

ε (ω)μ(ω),

(2.36)

 

c

 

 

где ε (ω)μ(ω) – показатель преломления.

26

Выражение (2.36) представляет собой дисперсионное уравнение среды. Если ε и μ не зависят от частоты, имеет место линейная дисперсия. При этом фазовая скорость

v

= ω

=

c

(2.37)

 

φ

k

 

εμ

 

 

 

 

 

меньше фазовой скорости света в вакууме и не зависит от частоты.

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

2.5. Второе интегрирование при определении векторных функций E (r,t ) и

H (r ,t ). Процесс интегрирования и анализ результата

Для упрощения дальнейших математических выкладок выберем декартову систему координат, ось x которой совпадает с направлением вектора k . При этом проекции век-

тора k на оси y и z равны нулю ( ky = 0, kz = 0 ), а kx = k . Тогда выражение (2.34) су-

щественно упрощается:

E (x,t )= −4π2

j st (ω,k )

ei(

ω

kx)dk .

 

 

t

 

(2.38)

ε (ω)

 

 

 

−∞

 

 

 

 

 

Введем обозначение A(k )= −4π2 j st (ω, k )ε (ω). При этом равенство (2.38)

принимает вид

 

 

 

 

 

 

 

 

E (x,t )=

A(k )ei(ωtkx)dk .

 

 

(2.39)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

Рассмотрим волновой пакет (рис.2.1). Полагаем, что

A(k )

принимает постоянное

значение A0

в узкой полосе 2δk значений k в окрестности k0 , а вне этой полосы A(k )

равно нулю.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

A ,

 

k k

 

 

 

 

<δk

k

 

δk

 

 

 

 

 

A(k )= 0

 

 

0

 

 

(2.40)

 

 

 

 

 

 

 

 

 

 

 

 

k

 

k

 

+δk

 

 

0,

 

k k

 

 

 

 

>δk

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

0

 

0

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис.2.1. Волновой пакет

27

В соответствии с (2.40) выражение (2.39) преобразуется к виду

 

E

(x,t )=

k

+δk

A ei(ωtkx)dk .

 

(2.41)

 

0

 

 

 

 

 

 

k0

δk

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Учитывая узость полосы 2δk ,

которая проявляется в том, что

δk <<1,

разложим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k0

 

ω(k ) в ряд Тейлора в окрестности точки k0 :

 

 

 

 

 

 

 

 

 

ω(k )=ω(k0 )

+ dω

 

k=k

 

(k k0 )+ 1 d 2ω

 

k=k

 

(k k0 )2 +.....

(2.42)

 

 

 

 

 

dk

 

 

 

0

 

 

 

2 dk2

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

k0

+δk

i ω(k0 )+dω

 

k =k0

(kk0 ) tikx

 

 

 

 

 

 

E

(x,t )= A

e

dk

 

 

 

dk .

(2.43)

 

0

 

 

 

 

 

 

 

 

 

k0

δk

 

 

 

 

 

 

 

Обозначим ω(k0 )=ω0 . Введем новую переменную α = k k0 . При этом dk = dα ,

нижний и верхний пределы интегрирования принимают значения δk и δk . После замены переменной интеграл (2.43) принимает вид

 

E (x,t )=

A ei[ω0tk0 x]

δk

i dω

 

k =k

tx α

dα .

 

 

 

 

 

 

 

 

 

e

dk

 

 

 

0

 

 

 

 

 

 

0

δk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dω

 

k=k0 t

 

 

δk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin

 

 

 

 

x

 

E

(x,t )= A ei[ω0tk0 x]2δk

 

 

dk

 

 

 

 

 

 

.

(2.44)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

dω

 

 

 

 

 

 

 

δk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dk

 

 

k=k0 t x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Анализируя выражение (2.44), можем сделать такие физические выводы. На часто-

те ω0 распространяется волна с модулированной амплитудой. Фазовый множитель этой

волны представлен экспонентой ei[ω0tk0 x]. Постоянному значению фазы соответствует равенство

ω0t k0 x = const .

(2.45)

Время течет непрерывно, а потому ω0t непрерывно увеличивается. Для того, чтобы ра-

венство (2.45) оставалось справедливым, должно увеличиваться и x . Надо узнать, с какой скоростью должно бежать вдоль оси x постоянное значение фазы. Иными словами, с какой скоростью должно увеличиваться x , чтобы скомпенсировать увеличение за счет вре-

28

мени первого слагаемого в (2.45). Для этого продифференцируем равенство (2.45) по вре-

мени t . При этом с учетом (2.36) получаем фазовую скорость dxdt в таком виде:

v

= dx

= ω0 =

c

.

(2.46)

 

φ

dt

k0

εμ

 

 

 

Закон модуляции амплитуды (пространственно-временная форма огибающей волны) определяется множителем

 

 

dω

 

 

δk

 

 

 

 

 

 

 

 

sin

 

k=k0 t x

= sinη .

 

A(x,t )

 

dk

 

 

 

(2.47)

 

 

 

 

 

 

 

 

 

 

dω

 

 

δk

η

 

 

 

 

 

 

dk

k=k0 t x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

dω

 

k=k0 t x = const .

(2.48)

 

dk

 

 

 

 

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

vгр =

dω

 

 

k=k0 ,

(2.49)

 

 

dk

 

 

 

 

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

2.6. Соотношение неопределенности

Проанализируем функцию (2.47). Из ее вида (рис.2.2) следует, что основная энергия сигнала заключена между двумя нулями по обе стороны от максимума, когда аргумент η принимает значения ±π , иными словами, когда выполняется соотношение

η 2π .

(2.50)

29

В формуле (2.47) η определяется следующим выражением:

 

η = dω

 

k=k0 t x

δk .

(2.51)

 

dk

 

 

 

 

Величина η может изменяться за счет изменения времени t или координаты x . Зафикси-

руем координату x . В выбранной точке пространства в соответствии с (2.50) должно вы-

полняться условие

dω

 

 

k=k0

tδk 2π . Поскольку

dω

 

 

k=k0 δk =

ω ( ω – прираще-

 

 

 

dk

dk

 

 

 

 

 

 

 

ние ω при малом изменении k ), приходим к соотношению неопределенности

 

 

 

 

 

ω t 2π .

(2.52)

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

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

Зафиксируем теперь момент времени. Изменение η при этом происходит за счет

изменения координаты x . Соотношение неопределенности принимает вид

 

xδk 2π .

(2.53)

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

30

3.ПЛОСКИЕ ВОЛНЫ В АНИЗОТРОПНЫХ СРЕДАХ

3.1.Общие положения

Анизотропной называется среда, диэлектрическая проницаемость которой пред-

 

ε

ε

12

ε

 

 

 

11

 

13

 

. Это симметричный тензор ( εik = εki ), i, k =1, 2,3.

ставлена тензором εˆ =

ε21

ε22

ε23

 

 

ε 31

ε32

ε33

 

 

 

 

 

Всегда можно выбрать систему координат, в которой этот тензор приобретает диагональ-

 

εxx

0

0

 

 

 

 

 

 

εyy

 

 

 

εxx εyy εzz , то это двухосный кристалл. В

ную форму: εˆ =

0

0

. Если

 

0

0

ε

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

zz

 

 

 

случае εxx =εyy =ε ,

εzz ε ,εzz =ε

 

имеем одноосный кристалл. Его диэлектриче-

ская проницаемость приобретает вид:

 

 

 

 

 

 

 

 

ε

0

0

 

 

 

 

 

 

0

ε

 

 

 

 

 

 

εˆ =

0 .

 

 

 

 

 

 

0

0

 

 

 

 

 

 

 

εII

Рассмотрим закономерности распространения однородной плоской волны в такой среде. Вектор электрического поля плоской волны представим в виде:

θ1

E (r ,t ) = E0e

iωtikr

.

1

 

 

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

2

границу раздела двух изотропных сред луч пре-

θ2

ломляется. Угол падения θ0 =θ1 (отражения) и

 

Рис.3.1. Преломление луча

угол преломления θ2 связаны

соотношением:

n1 sinθ1 = n2 sinθ2 , где ni ( i =1,2 ) – показате-

на границе раздела двух сред

 

ли преломления первой и рис.3.1. Преломление

луча на границе раздела двух сред второй среды, определяемые как ni =

εiμi . Направ-

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

связан с волновым числом простым соотношением: n = k

 

ω

. Необходимо выяснить,

 

 

 

 

c

 

что собой представляет показатель преломления n в анизотропной среде. Но сначала про-

31