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

Воробьев Теория электромагн поля и СВЧ (Кривець)

.pdf
Скачиваний:
1062
Добавлен:
23.03.2016
Размер:
2.65 Mб
Скачать

101

Выражение (3.2) описывает так называемые незатухающие волны. Если волна затухает или нарастает в

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

В электронной технике обычно имеют дело с гармоническими волновыми процессами, которые

описывают

гармоническими

функциями

аргумента

ω(t ± z / vф )=ωt ± kz . Здесь ω = 2π f

круговая частота

временных

колебаний;

k =ω / vф

фазовая

постоянная

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

волновое

число; vф

– фазовая

скорость (скорость распространения фронта волны). Эти параметры волны связаны с временным периодом T и

пространственным периодом

λ

(длина

волны)

однотипными соотношениями:

ω = 2π /T ,

k = 2π / λ ,

vф = 2π / k .

 

 

 

Перейдём к рассмотрению основных видов волновых уравнений. Определить структуру и другие характеристики поля, непосредственно используя уравнения Максвелла, затруднительно, поскольку в каждом из них присутствуют по два неизвестных

→ → →

параметра: H , D и E , B соответственно. Поэтому для описания поля в однородной линейной среде при

отсутствии токов (δ = 0 ) и свободных зарядов ( ρ =0)

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

(δ 0 ) и свободных зарядов ( ρ 0) используются неоднородные волновые уравнения Даламбера.

102

Уравнения Гельмгольца можно получить из первого и второго уравнений Максвелла в дифференциальной форме записи (1.9), (1.10):

rot H = ε

 

(δ = 0 ); rot E = −µ

 

 

E

 

H .

 

 

 

 

 

a

t

 

 

a

t

Подействуем оператором ротора на левую и правую части уравнений:

 

 

 

 

 

 

 

 

 

 

 

 

 

rotrot H = εa

 

 

 

 

rot E

;

rotrot E = −µa

 

 

 

rot H

.

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

Подставим

 

 

 

значения

 

 

 

 

 

 

 

 

и

 

 

 

 

 

 

 

 

из

 

 

 

 

 

rot E

 

 

 

 

 

 

 

 

 

 

rot H

 

вышеприведенных уравнений Максвелла:

 

 

 

 

 

 

 

 

 

 

∂ ∂

 

 

 

 

 

 

 

 

 

 

 

∂ ∂

 

 

rotrot H = −εa

 

 

 

 

 

 

 

µa H

 

;

rotrot E

= −µa

 

 

 

 

 

 

εa E .

 

 

 

 

 

 

 

 

 

 

 

t t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t t

 

 

 

Применив

 

тождество

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rotrot K

= graddiv K 2

K ,

можно записать

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

µ

 

2

,

 

 

 

 

 

 

 

graddiv H 2

H = −ε

a

a

 

H

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

µ

 

 

2

 

 

 

 

 

 

 

 

 

 

graddiv E2

E = −ε

a

a

 

 

E .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С учётом того, что div H

0 (δ = 0 ) и div E = 0 ( ρ =0),

а εa µa имеет размерность, обратную квадрату скорости

распространения волны, то уравнения Гельмгольца примут вид:

 

1

2

 

1

2

 

 

2

H

 

H

= 0 , 2

E

 

E

= 0 .

(3.3)

v2

t2

v2

t2

 

 

 

 

 

 

В отличие от уравнений Максвелла каждое из полученных дифференциальных уравнений второго

→ →

порядка зависит только от одной переменной H или E и в

103

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

2 Hx, y,z

+

2 Hx, y,z

+

2 Hx, y,z

1

 

2 Hx, y,z

= 0 .

x2

y2

z2

v2

 

t2

 

 

 

 

 

Аналогичным образом запишутся и уравнения для

вектора E :

2 Ex, y,z

+

2 Ex, y,z

+

2 Ex, y,z

1

 

2 Ex, y,z

= 0 .

x2

y2

z2

v2

 

t2

 

 

 

 

 

Воспользовавшись первым и вторым уравнениями Максвелла в комплексной форме записи (1.14), можно получить уравнения Гельмгольца для гармонических колебаний:

 

 

2

(3.4)

Ep2

E = 0 ; 2

H p2

H = 0 ,

где p = jω εka µa – коэффициент распространения волны.

Неоднородное волновое уравнение Даламбера для

векторного потенциала A .

Воспользуемся первым уравнением Максвелла в дифференциальной форме (1.9) и домножим его на µa :

 

µ rot H = µ δ+ µ ε

 

 

E .

 

 

 

 

 

 

 

 

 

 

 

a

 

 

a

 

a a

t

 

 

Внесём µa

под оператор rot , тогда с учётом того,

что

→ →

 

 

 

 

 

 

 

 

 

µa H = B = rot A , получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 E

 

 

 

rotrot A

= µa δ+

 

 

 

.

 

 

v2

t

 

 

 

 

 

 

 

 

 

Применив

тождество

 

 

 

 

 

rotrot K

= graddiv K 2

K ,

можно записать

104

 

 

 

 

 

 

 

1

E .

 

graddiv A2

A = µa

δ+

(3.5)

v2

 

 

 

 

t

 

Чтобы перейти от

вектора

к вектору

E

A ,

воспользуемся вторым уравнением Максвелла (1.10) и выражением магнитной индукции через векторный потенциал (1.5).

rot E = − B

= −rot A .

 

t

t

Если равны роторы от двух функций, то равны и их функции с точностью до градиента от некоторой скалярной функции (так как rotgradϕ 0 ), т.е.

 

 

 

A

gradϕ .

(3.6)

E = −

 

t

 

 

(Это несложно проверить: если на (3.6) подействовать оператором ротора, то получится второе уравнение

Максвелла, а для стационарного поля связь E = −gradϕ ). Подставив в (3.6) в (3.5), имеем:

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

graddiv A2

A = µa

δ

 

 

 

 

 

A + gradϕ

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

v

 

t

 

 

 

t

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 2

 

 

 

1 ϕ

 

 

2

 

A

 

grad

div A+

 

 

 

 

 

 

 

A+

 

 

 

 

 

2

= µa

δ . (3.7)

v

2

t

 

 

v

2

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для перехода

от (3.7)

к

стационарному

уравнению

Пуассона 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A = −µa

δ

(

 

 

 

0 )

 

необходимо выполнение

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

равенства

105

1

ϕ .

 

div A = −

(3.8)

 

 

v2 t

 

Последнее выражение

называется

калибровкой

Лоренца.

 

 

 

Окончательно получаем уравнение Даламбера для векторного потенциала:

 

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

2

A

 

A

= −µa δ .

(3.9)

 

 

 

v2

t2

 

 

 

 

 

 

 

 

 

 

 

Уравнение Даламбера для векторного потенциала

можно

 

записать

через

 

четырёхмерный

Лапласиан

( 2 =

2

1 2

 

 

 

 

 

 

 

 

 

 

 

 

):

 

 

 

 

 

 

 

 

 

v2

t2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 A = −µa δ .

(3.10)

Неоднородное волновое уравнение Даламбера для скалярного потенциала ϕ .

Воспользуемся третьим уравнением Максвелла в дифференциальной форме (1.11) и подставим в него

значение E из (3.6):

 

 

 

 

 

ρ

 

 

 

ρ

div E =

; div

A

+ gradϕ

= −

εa

t

εa

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

ρ

 

 

 

div A

+ graddivϕ = −

 

.

 

 

εa

 

t

 

 

 

С учётом калибровки Лоренца (3.8) получим уравнение Даламбера для скалярного потенциала, которое также можно записать через четырёхмерный Лапласиан:

2ϕ = −

ρ

.

(3.11)

 

 

εa

 

106

Для большей наглядности сгруппируем полученные волновые уравнения в виде табл. 3.1.

Таблица 3.1 – Волновые уравнения

Вид

Электрические

 

 

Магнитные

 

уравнения

составляющие поля

составляющие поля

 

Гельмгольца

 

 

1

2

 

 

 

 

 

 

 

 

1

2

 

 

2 E

E = 0 ;

 

2 H

 

 

H = 0 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t2

 

 

 

 

v2

t2

 

 

 

 

 

 

 

v2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

E = 0

 

 

 

 

 

 

2 H = 0

 

Даламбера

2

1 2ϕ

 

 

ρ

 

 

 

 

 

 

2

 

ϕ

 

 

 

 

2

= −

 

 

;

2

A

1

 

 

 

 

A

= −µa δ

;

 

v

2

t

εa

 

 

 

 

 

v2 t2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2ϕ = − ρ

 

 

 

 

 

2 A = −µa δ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

εa

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

3.2 Параметры плоской волны в однородной среде

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

107

Для описания распространения электромагнитных волн используется понятие фазового фронта

поверхности, проходящей через точки с одинаковыми фазами.

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

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

По форме фазового фронта, кроме плоской волны,

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

(рис. 3.2).

H

y

E

 

 

 

x

z

 

Рисунок 3.2 – Схематическое изображение фронта волны, радиально расходящейся от источника

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

величины которого E и H в каждый момент времени во всех точках плоскости (x0y), перпендикулярной

108

направлению распространения волны (ось z ), принимают

 

 

 

 

постоянные значения ( E

= E

= H

= H

= 0 ), т.е.

E и

 

x

y

 

x

y

 

 

перпендикулярны

и

зависят

только

от

H взаимно

координаты z

и не зависят от x

и

y . Обычно векторы

E и H изменяются вдоль оси z по гармоническому закону (закону синуса или косинуса).

Определим закон распространения плоской электромагнитной волны, т.е. найдём волновую функцию плоской волны. Для этого решим волновые уравнения Гельмгольца в комплексной форме (3.4), которые для плоской волны запишутся следующим образом:

 

d

2

2

 

 

 

E

p2

E = 0 ;

d

H

 

p2

H = 0 ,

(3.12)

 

dz2

dz2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где p – коэффициент распространения;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p = jω εka µa

=α + jβ ;

(3.13)

α – коэффициент затухания; β – коэффициент фазы.

Решение уравнения Гельмгольца будем искать в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E = Em ek z .

 

 

(3.14)

В результате подстановки (3.14) в (3.12) получаем

характеристическое уравнение

 

 

 

 

 

 

 

 

 

 

k2 p2

= 0,

 

 

откуда k = ±p . Тогда

 

 

 

 

 

 

 

 

 

 

ep z = E+ eα z

ejβ z

+ Eeα z

e jβ z .(3.15)

E = E+ ep z + E

m

 

 

m

m

 

 

 

m

 

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

109

умножить на e jωt и взять действительную часть, используя тригонометрическую формулу Эйлера

Re(e j(ωtβ z) )= cos(ωt β z).

В результате получим

E = Em+ eα z cos(ωt β z +ϕ+ )+

(3.16)

+Emeα z cos(ωt + β z +ϕ).

Здесь ϕ+ и ϕ– фазы волны, зависящие от выбора

начала отсчёта.

Аналогичные выражения можно получить из второго уравнения (3.12) для H -компоненты поля.

Если принять, что на пути распространения волны нет преград, то амплитуды отражённой волны Em, Hmдолжны

равняться нулю и решения для плоской электромагнитной волны примут вид:

E = Em+ eα z cos(ωt β z +ϕ+ );

(3.17)

H = Hm+ eα z cos(ωt β z ϕ+ ).

 

Следует отметить, что для незатухающей волны (отсутствуют потери в среде) графики решений (3.17) будут по характеру идентичны рис. 3.1. Если среда с потерями (α 0), то амплитуды гармонических функций

будут затухать по экспоненциальному закону eα z .

Основными параметрами, характеризующими распространение электромагнитной волны, являются:

коэффициент распространения p ;

коэффициент фазы β ;

коэффициент затухания α ;

фазовая скорость vф ;

110

длина волны λ ;

волновое сопротивление Zв ;

глубина проникновения волны .

Коэффициент распространения p является

комплексной величиной, которая характеризует изменение амплитуды и фазы бегущей электромагнитной волны и для плоских однородных волн при заданной частоте ω определяется только параметрами среды (εa , µa и γ ).

Коэффициент распространения p (3.13) в общем виде

может быть записан как через коэффициенты затухания и фазы (α , β ), так и через тангенс угла потерь ( tgδ ) с

учётом того, что εka = εa (1jtgδ ):

p =α + jβ = jω µaεa (1jtgδ ) =

 

j

ω

(3.18)

=

µε (1jtgδ ).

 

 

c

 

Коэффициент фазы β показывает изменение фазы

волны при прохождении 1 м пути и равен мнимой части коэффициента распространения p .

Коэффициент затухания α определяет уменьшение амплитуды волны при прохождении 1 м пути и равен действительной части коэффициента распространения p .

Найдём α и β из общего выражения для

коэффициента распространения (3.18), проделав ряд несложных тригонометрических операций:

 

2

 

 

 

 

 

 

2

( 1+tg2δ ); (3.19)

 

 

 

 

 

 

p

 

=

 

α2 + β

2 2

=α

2 + β2 =

ω2 εµ

 

 

 

 

 

 

 

 

 

c

 

 

 

 

 

p2

= −

ω2

µε (1jtgδ );

 

 

 

 

 

 

 

 

c2