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

Кураев А.А., Попкова Т.Л., Синицын А.К. Электродинамика и распространение радиоволн

.pdf
Скачиваний:
342
Добавлен:
02.05.2014
Размер:
5 Мб
Скачать
☆

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

К настоящему времени предложены более совершенные приближенные методы учета дифракционных полей, из которых следует упомянуть метод дифракционных лучей, в котором учитываются «ползущие» вокруг поверхности препятствия волны, связанные с затеканием токов на теневую сторону (П.Я. Уфимцев «Метод краевых волн в физической теории дифракции». М., Сов. Радио, 1962 г.), а также более сильный метод параболического уравнения, в котором учитывается поперечная диффузия дифракционных лучей (В.А. Фок, Дифракция радиоволн вокруг земной поверхности. АН СССР, 1946 г.).

251

ГЛАВА XIX

ОТРАЖЕНИЕ И ПРЕЛОМЛЕНИЕ ВОЛН НА ГРАНИЦЕ РАЗДЕЛА ДВУХ СРЕД

19.1. Отражение и прохождение волн при нормальном падении на плоскую границу раздела

Пусть плоская электромагнитная волна, проходящая в среде 1 падает нормально на плоскую границу, разделяющую среду 1 и среду 2 (рис. 19.1). Среда 1 характеризуется параметрами ε1 , µ1 , среда 2 - ε2 , µ2 (значок «а» у εа иµа здесь и далее опустим для сокращения записи). Непосредственно на

границе раздела существуют три вида волн: падающая (Е0 , Н0 ), отраженная (Е− , Н− ), проходящая (Е+ , Н+ ). Мы будем предполагать, что как среда 1, так и среда 2 однородны и изотропны, поэтому волны в обеих средах распространяются прямолинейно и нормально к границе, как это изображено на рис. 19.1.

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 19.1

 

 

r

 

 

Направления распространения

волн

определяются векторами

=

1

r r

в отраженной волне изменена по от-

S0

[E, H], поэтому ориентация E− , H−

 

 

2

 

 

 

 

 

ношению к падающей, так что направления S00 иS0−

противоположны.

 

 

 

Запишем компоненты трех волн, придерживаясь условной ориентации

Е, Н, заданной на рис. 19.1 (действительная относительная ориентация, т.е.

знак компонент (Е− , Н− ) будет определен соотношением ε1 , µ1 ; ε2 , µ2 ).

 

1.

Падающая волна:

 

 

 

 

 

 

r

. 0

r

 

 

= ω ε1µ1

,

 

 

 

 

H

0 = A y0e−jk1z , k1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

(19.1)

r

. 0

r

 

e−jk1z , W0

= µ

 

ε

 

 

 

E0 = A W0 x

0

1

1

 

 

 

 

1

 

1

 

 

 

 

2.

Отраженная волна:

 

 

 

 

 

 

252

r

. − r

 

 

 

 

 

 

 

 

H

− = −A y0 e jk1z ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(19.2)

 

 

 

 

 

 

 

 

 

r

. −

r

 

 

 

 

 

 

 

E− = A W0 x

e jk1z

 

 

 

 

 

 

 

1

0

 

 

 

 

 

 

 

3. Проходящая волна:

 

 

 

 

 

 

r

. + r

 

 

 

 

 

 

 

 

H

+ = A y0e−jk2z , k 2 = ω ε2µ2 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

(19.3)

 

 

 

 

 

 

 

 

 

r

. +

r

e−jk2z , W0

= µ

 

ε

 

 

 

E+ = A W0 x

2

2

 

 

 

2

0

2

 

 

 

 

Запишем уравнение связи между компонентными амплитудами трех волн на границе раздела двух сред (z = 0). Предположим, что на границе раздела поверхностных токов нет; тогда тангенциальные составляющие как Е, так и Н полного (суммарного) поля при переходе границы должны быть непрерывными (см. главу II первой части). В рассматриваемом случае все компоненты Е и Н являются тангенциальными, и поэтому

r

0 +

r

−

 

=

r

 

,

r

r

r

 

.

H

H

 

H

+

E0

+ E

= E+

 

 

 

z =0

 

 

z =0

 

 

z =0

z

=0

 

(19.4)

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя в (19.4) (19.1) ÷(19.3), имеем

 

 

. 0

. −

 

. +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A − A

= A ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(19.5)

 

 

 

 

 

 

 

 

 

 

 

 

 

. 0 .

−

 

.

+

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

+ A

 

0

 

 

 

 

 

 

 

 

A

W1

= A W2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Введем для удобства дальнейших расчетов коэффициент отражения и коэффициент прохождения

. −

.

−

 

 

ρ. = Eτ (0)

=

 

A

 

− коэффициент отражения,

 

0

. 0

.

 

 

Eτ (0)

 

 

A

 

 

 

 

 

 

 

 

(19.6)

. +

.

+

 

 

τ. = Eτ (0)

=

A W20

− коэффициент прохождения.

 

. 0

.

0

 

 

Eτ (0)

 

 

A W0

 

 

 

 

 

 

1

 

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

253

. − .

0 .

 

. +

.

0 .

W0

 

A

 

= A

 

ρ,

 

A

= A

τ

1

.

(19.7)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

Подставляя (19.7) в (19.5), получим уравнения для ρ.

и τ. :

 

. .

0

 

 

 

 

 

 

 

 

1−ρ = τ

W1

,

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

W

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(19.8)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

. .

 

 

 

 

 

 

 

 

 

 

1+ρ = τ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решая (19.8), находим:

 

.

 

 

 

0

 

 

 

 

 

 

 

 

 

 

τ =

 

2W2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

W

0

+ W

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(19.9)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

0

 

 

0

 

 

 

 

 

 

 

 

ρ =

 

W2

− W1

 

 

 

 

 

 

 

 

 

 

W0

+ W0

 

 

 

 

 

 

 

 

 

2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим

 

предельные

случаи для

 

 

(19.9),

т.е.

случаи, когда

 

ρ.

 

 

= 0 и

 

ρ.

 

=1.

 

 

 

 

 

 

 

 

 

 

 

ρ.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1. Случай согласования двух сред, т.е.

 

 

= 0 . Из (19.9) следует, что это

 

 

 

 

 

 

возможно при условии

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W20 = W10 или

µ2

=

µ1 .

 

 

 

 

 

 

 

 

(19.10)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ε

2

 

ε

1

 

 

 

 

 

 

 

 

 

При этом τ =1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2. Случай полных отражений от границы раздела при нормальном па-

дении, т.е.

 

ρ.

 

=1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из анализа формулы для ρ.

следует, что | ρ.

|→1 при:

 

 

 

 

 

 

 

W0

 

 

 

W0

 

 

 

 

 

 

 

 

а)

 

2

→ 0 (ρ = −1)или б)

 

1

→ 0 (ρ =1).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W0

 

 

 

W0

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Условие а) выполняется, когда среда 2 по своим свойствам приближа-

ется к идеальному

проводнику.

Действительно, W20 =

µ

σ → 0 при σ → ∞.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ε − j

ω

При этом τ = 0 , т.е. поле во второй среде отсутствует.

 

 

 

 

 

 

 

 

 

Условие б) соответствует случаю, когда волна распространяется в сре-

де с большой оптической плотностью (например, в диэлектрике с очень большим ε ), а среда 2 имеет относительно большое волновое сопротивление. В этом случае получаем, что ρ=1, а τ = 2 . Последнее, однако, не говорит о

254

том, что поле в среде 2 есть. Действительно, используя формулу (19.7), име-

.

+ . 0 0

. 0

0

 

0

 

ем: A

= A τ

W1

= A 2

W1

→ 0 при

W1

→ 0 .

W0

W0

W0

 

 

 

 

 

 

2

 

2

 

2

 

Тот факт, что τ = 2 лишь означает, что электрическое поле на границе раздела по амплитуде вдвое превышает амплитуду электрического поля падающей волны ( Е(0)=2 Е0(0)), магнитное же поле на границе равно нулю

( Н(0)=0).

При ρ. =1 в области 1 образуется картина чисто стоячего поля: ампли-

туда отраженной волны равна амплитуде падающей волны. В общем же случае имеем:

r

&

0 r

− jk z

 

.

jk z

&

0 r

 

 

 

.

2 jk z

 

− jk z

r

&

H

− jk z

Н1

= А у0

е

 

− ρ е

 

 

= А у0

1 −

ρе

 

 

е

 

= у0 А (z)e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

&

0 r

0

.

 

2 jk z

 

− jk z r

&

E

 

 

− jk z

 

 

 

 

 

 

E1

= А х0W1

1 +

ρе

 

 

е

 

= х0

А (z)e

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из приведенных формул видно, что при любом ρ. распределение амплитуды электрического поля А&E (z) и распределение амплитуды магнитного поля А&H (z) вдоль z сдвинуты по фазе на λ1 4 .

В минимуме поле обращается в нуль только при ρ. =1; нули

E1 (z) и H1 (z) смещены на λ1 4 .

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

Введем следующие предположения.

1.Падающая на границу раздела волна – линейно поляризованная и ее электрический вектор параллелен границе раздела («горизонтальная» поляризация).

2.1 и 2 среды изотропные и однородные.

Последнее предположение позволяет считать, что плоскость поляризации отраженной и проходящей волн совпадают с плоскостью поляризации падающей, а направления распространения падающей, отраженной и проходящей (преломленной) волн лежат в одной плоскости. На этом основании чертеж, изображающий направления распространения трех волн и ориентацию векторов Е, Н в них, будет иметь вид, показанный на рис. 19.2. Система x, y, z – расчетная, в ней наиболее просто задается граница раздела сред. В системе x1, y1, z1 наиболее просто описывается распространение падающей волны ( Е0, Н0). Такой же смысл имеет система x2, y2, z2 в отношении отра-

255

женной волны, а система x3, y3, z3 – в отношении преломленной волны. Указанные системы z1, z2, z3 развернуты по отношению к расчетной в плоскости чертежа (в «плоскости падения») соответственно на угол ϕ1 , ϕ2 , ϕ3 .

Правила пересчета направляющих векторов развернутых систем координат в расчетной системе состоят в следующем (рис. 19.3):

r

 

 

r

 

 

 

 

 

 

 

 

x

0i

= x0 ,

 

 

 

r

 

 

r

 

 

r

cos ϕi

 

 

 

 

y

0i

= y0

− z0

sin ϕi ,

(19.11)

r

 

 

r

sin ϕi

 

r

 

 

z0i

= y0

+ z0

cos ϕi ,

 

zi = ysin ϕi

+ z cos ϕi .

 

 

 

 

Запишем поле падающей волны сначала в системе координат z1

 

r

0

. 0

r

 

−jk z

 

 

 

H

= A

у01е

 

1 1 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(19.12)

 

 

 

 

 

 

 

 

 

 

 

r

0

. 0

0 r

 

 

 

 

 

 

 

 

−jk1z1

 

 

E

 

= A W1

х01е

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 19.2 Рис. 19.3

Пользуясь правилами (19.11), перепишем (19.12) в расчетной системе z

r

 

. 0

r

cos

ϕ1

r

sin ϕ1 )e−jk1

 

 

H0

= A

(у0

− z0

(y sin ϕ1 +z cos ϕ1 ),

 

 

 

 

 

 

 

 

 

 

 

 

(19.13)

 

 

 

 

 

 

 

 

 

 

 

r

0

. 0

r

0

е

−jk

(y sin ϕ +z cos ϕ )

 

 

E

 

= A x0 W1

1

 

1

1

 

 

 

 

 

 

 

 

 

 

 

Аналогичным образом запишем и поля отраженной и проходящей волн в расчетной системе координат:

256

r

.

−

r

cos ϕ2

r

sin ϕ2 )е−jk1

 

H− = A

(y0

− z0

(y sin

r

.

−

 

r

 

е−jk1

(y sin ϕ2 +z cos ϕ2 ),

 

E− = A

 

W0 x

0

 

 

 

 

1

 

 

 

 

 

r

.

+

r

cos ϕ3

r

sin ϕ3 )е−jk2

 

H+ = A

(y0

− z0

(y sin

r

.

+

 

r

 

е−jk2

(y sin ϕ3 +z cos ϕ3 )

 

E+ = A W0 x

0

 

 

 

 

l

 

 

 

 

 

ϕ2 +z cos ϕ2 ),

(19.14)

ϕ3 +z cos ϕ3 ),

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

текают следующие уравнения, формулирующие законы Снеллиуса:

 

1)

k1 sin ϕ1

= k1 sin ϕ2 ,

(19.15)

2)

k1 sin ϕ1

= k2 sin ϕ3 .

(19.16)

Уравнение (19.15) имеет бесконечное множество решений из-за периодичности функции sin . Выберем решение, соответствующее физическому смыслу задачи. Решение ϕ1 = ϕ2 очевидно, не подходит, поскольку отраженная волна распространяется в среде 1, а не в среде 2. Приемлемым, повидимому, будет решение

ϕ1 = π−ϕ2 = γ

 

 

(19.17)

Обозначим далее ϕ1 = ϕ - угол падения. Тогда (19.17) формулирует пер-

вый закон Снеллиуса: угол падения (ϕ) равен углу отражения (γ).

 

Обозначая ϕ3 = ϑ - угол преломления, из (19.16) получим формулиров-

ку второго закона Снеллиуса:

 

sin ϑ

= k1 =

ε1µ1

= n1

(19.18)

sin ϕ

k2

ε2µ2

n2

 

 

 

 

В (19.18) по аналогии с принятыми в оптике терминами обозначено

 

k1

=

n1

, где n1 – показатель преломления среды 1, n2 – показатель преломле-

 

k2

 

 

 

n2

 

ния среды 2.

 

 

 

 

В соответствии с первым законом Снеллиуса имеем:

 

 

 

 

cosϕ2 = −cos γ = −cosϕ

(19.19)

Это условие учтем в дальнейшем при использовании (19.14).

257

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

H

yI

z=0

= H

yII

z=0

,

 

 

(19.20)

 

 

 

 

 

ExI

z=0

= ExII

z=0 .

 

 

 

Перепишем условия (19.20), используя (19.14) и (19.19)

 

. 0

.

−

 

 

. +

 

 

 

 

 

− A

 

 

ϕ = A

cosϑ,

 

A

cos

 

 

 

 

 

 

 

 

 

 

(19.21)

. 0

 

−

 

. +

 

 

.

0

0

 

 

 

 

+ A

 

 

= A W2

.

 

A

W1

 

 

 

 

 

 

 

 

 

 

 

Введем, как и ранее, коэффициенты отражения и прохождения.

 

−

(0)

 

. −

 

 

+

(0)

.

+

0

.

 

A

 

.

 

A

 

ρг =

Ex

=

,

τг =

Ex

=

 

 

W2

E0

(0)

. 0

E0

(0)

.

0

 

W0

 

x

 

 

A

 

 

x

 

 

A

 

1

 

 

 

 

 

 

 

 

 

 

 

 

Теперь имеем

.

− .

0 .

.

+ .

0 .

W0

A

= A

ρ

; A

= A

τг

1

W0

 

 

г

 

 

 

 

 

 

 

 

 

2

(19.22)

(19.23)

Подставляя (19.23) в уравнения (19.21), получаем следующую систему

для определения ρ. г , τ. г :

.

.

W0

cos ϑ

,

.

.

1−ρг = τг

W0

cos ϕ

1+ ρг = τг .

 

 

1

 

 

 

 

 

 

2

 

 

 

 

Решая (19.24) получаем

.

 

 

0

ϕ

 

 

τг =

 

 

2W2 cos

 

 

 

 

 

 

,

W

0

cos ϕ + W

0

cos ϑ

 

 

 

 

 

2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

0

 

0

cos ϑ

 

ρг =

W2

cos ϕ − W1

.

 

0

 

0

 

 

W2

cos ϕ + W1

 

 

 

cos ϑ

(19.24)

(19.25)

Проанализируем полученные формулы для ρ. Γ , τ. Γ в одном из наиболее важных в теории распространения волн в случае двух идеальных диэлектриков. В этом случае

258

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ε1'' = ε'2' = 0, µ1 = µ2

= µ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W

0

ε

 

 

sin ϑ

 

ε

 

 

 

 

 

 

 

 

 

 

 

 

 

,

=

 

,

 

 

 

 

 

 

 

(19.26)

1 =

ε

2

sin ϕ

ε

1

 

 

 

 

 

 

 

W

0

1

 

 

2

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ε1

 

 

 

 

 

 

 

 

cosϑ = 1−sin

2

ϑ =

1−

sin

2

 

 

 

 

 

 

 

ε

 

 

ϕ.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Учитывая (19.26), получаем

 

 

 

 

 

 

 

 

 

 

ρ. г

=

ε1 cos ϕ −

 

ε2

− ε1 sin 2 ϕ .

 

 

 

 

(19.27)

 

 

ε1 cos ϕ +

 

ε2

− ε1 sinε ϕ

 

 

 

 

 

 

 

Рассмотрим два случая: ε2 > ε1 и ε2 < ε1 .

 

 

1)

ε2

> ε1 , т.е. волна падает из менее плотной среды (например, воздуш-

ной) на более плотную. Представим ρ.

г = ρге−jΦΓ

и проанализируем зависимо-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 19.4

сти ρг , Φг от ϕ:

 

 

 

 

а) ϕ = 0 → ρ. г =

ε1

−

ε2

< 0, т.е. ρг <1, Φг =1800.

 

ε1

+

ε2

 

б) ϕ → 900 , ρ. г → −1

Φг

=1800 , ρг =1

Таким образом, зависимости ρΓ , ΦΓ от ϕ имеют вид, представленный на рис. 19.4.

2) ε2 < ε1 , т.е. волна распространяется внутри более плотной среды.

а) ϕ = 0 → ρ. г

=

ε1

−

ε2

> 0, т.е. ρг <1, ΦΓ = 0.

 

 

ε1

+

ε2

 

б) ϕ = 900 , ρ. г

= −1,

Φг =1800 , ρг =1.

259

Очевидно, что существует ϕ = ϕкр , при котором ρ. впервые становится

равным 1

 

.

 

. Условием этого, очевидно, является уравнение

 

ρ =1

 

 

 

 

 

ε2 −ε1 sin2 ϕкр = 0 , т.е.

 

sin ϕкр =

ε2

(19.28)

 

ε

 

 

1

 

При ϕ > ϕкр выражение для ρ. г ввиду того, что вторые члены в числите-

Рис. 19.5

ле и знаменателе оказываются мнимыми, принимает вид:

ρ. г = a − jb = cе− jα

=1 е−j2α ,

 

 

a + jb

cеjα

 

 

 

 

 

a = ε1 cos ϕ,

b =

ε1 sin 2 ϕ − ε2 .

 

 

Таким образом, при ϕ > ϕкр

 

ρ.

 

 

г

=1, Φг = 2α.

 

 

 

 

 

 

 

Графики ρг , Φг имеют вид, изображенный на рис. 19.5.

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

260