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

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

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

Em =

60PΓGm F(θ,ϕ)

(20.10)

 

r

 

Этой формулой мы и будем пользоваться в дальнейшем, обозначая

PГ = P0 .

20.2.Напряженность поля излучателя, поднятого над плоской

иоднородной поверхностью Земли на высоту h1

Для определенности будем иметь в виду вертикальный элементарный электрический излучатель и рассмотрим схему расположения излучателя и приемника, изображенную на рис. 20.1. Волна от источника (0) к приемнику (А) распространяется двумя путями – r1 и r2. В последнем случае волна отражается в области В от плоской земной поверхности. Направление на приемник по пути r1 определяется углом θ1 , по пути r2 −θ2 . Высоты подъема антенн передатчика и приемника соответственно h1 и h2. Распространение волны по пути r2 можно условно представить, как распространение по прямолинейному пути из точки О1 (разумеется, учитывается коэффициент отражения в области В, а свойства среды предполагается теми же, что над поверхностью земли).

Для расчета интерференционного поля в точке А введем следующие предположения: h1<<r1, h1<<r2. Тогда r1 ≈ r2 , r1 || r2 (рис. 20.2). Соответственно

 

 

 

Рис. 20.1

 

 

Рис. 20.2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

θ2 1800 −θ1 ,

r2

= r1 + 2h1 cosθ

(θ = θ1 )

(20.11)

Учитывая эти упрощения, запишем поле в точке А

 

. . . .

.

.

. .

 

 

(20.12)

E = E1

+ E2 , E1

= Em1 е−jkr1 , E2 = Em2 ρе−jkr2

Здесь

.

 

 

 

 

 

. .

 

ρ = ρе−jΦ - коэффициент отражения в точке В. Запишем

E1 , E2 че-

рез мощность генератора Р0

271

.

=

60P0 Gm

F(θ)е−jkr1 еjωt

,

 

 

E1

 

 

 

 

r1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E2

=

60P0G m

 

 

 

 

 

F(1800 − θ)ρе−jΦе−jkr2 еjωt .

.

 

 

 

 

 

 

 

 

 

 

r2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Используя

(20.11) и

 

предполагая, что F(θ) = F(1800 −θ),

(r2 ≈ r1 = r)

 

 

 

 

 

 

 

.

 

60P0G m

e

j(ωt−kr)

(1+ ρe

−j(2h1k cos θ+Φ)

)F(θ)

E =

r

 

 

 

 

 

 

 

 

 

 

 

Запишем теперь модуль комплексной амплитуды

Em =

60P0G m

F(θ)[1+ ρ2 + 2ρcos(2h1k cos θ + Φ)]1/ 2 .

 

r

 

(20.13)

получим

(20.14)

(20.15)

Формула (20.15) является основой для дальнейшего анализа влияния отражений от земной поверхности на диаграмму направленности излучателей. В случае горизонтального вибратора в формуле (20.15) следует угол θ заменить на угол ϕ .

Таким образом, как следует из сравнения формул (20.10) и (20.15), за счет отражений от земной поверхности в выражении Em появляется интерфе-

ренционный множитель I (θ) = 1+ ρ2

+ 2ρ cos(2h k cosθ +Φ ) 1/ 2 .

 

1

 

20.3. Поле вертикального вибратора над плоской отражающей поверхностью

Пусть свойства среды над отражающей поверхностью определяются как ε. 1 ,µ. 1 , свойства среды под отражающей поверхностью - ε. 2 ,µ. 2 . Соответст-

 

 

 

 

 

 

.

 

.

 

венно k1 = ω ε.

1 µ.

1 , k2

= ω ε.

2 µ.

2

, W10 = µ.

1 ,

W20 = µ.

2 .

 

 

 

 

 

 

ε1

ε2

Коэффициент отражения ρ. В для вертикально поляризованной волны выражается следующим образом:

ρ. B =

W10

cosθ− W20

cosϑ

(20.16)

W0

cosθ+ W0

cosϑ

 

 

 

1

2

 

 

Здесь ϑ - угол преломления. Используя второй закон Снеллиуса, исключим cos ϑ

272

| I(θ)| определяются условием
I(θ). Нули этой функции определяются ус-
(20.18)
θ = θmax
ности излучателя изменяется:

cosϑ =

1− k12 sin2 θ

 

 

(20.17)

 

 

 

 

k22

 

 

 

Подставляя (20.17) в (20.16), получим

ρ.

= k2 W10

cosθ− W20

k22

− k12 sin2

θ

B

k

W0

cosθ+ W0

k2

− k2 sin 2

θ

 

 

 

2

1

2

2

1

 

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

1. Отражение происходит от границы среды, которая может быть охарактеризована как проводник (ДВ, СВ над сушей и морем, КВ над морем). В этом случае W20 → 0, ρB =1, ΦB = 0 . Соответственно

Em =

60P0G m

F(θ)

2[1+ cos(2h1k cos θ)]=

60P0 G m

F(θ)2cos(h1k cos θ).

 

r

 

 

r

 

В отличие от поля в свободном пространстве здесь появляется интерференционный множитель I(θ)= 2cos(h1k cosθ). Теперь диаграмма направлен-

F1 (θ)= F(θ) I(θ).

Проанализируем функцию ловием

h1k1 cosθ0 = π2 , 32π,...

Таким образом,

cosθ0

=

2n +1

 

λ

, n = 0, 1, 2,...

4

 

 

 

 

 

 

h1

Соответственно максимумы

cosθmax = 2nh1 λ .

В направлениях θ = θ0 излучение отсутствует, в направлениях

напряженность поля оказывается вдвое большей, чем в свободном пространстве (| I(θ)| =2). Диаграмма направленности элементарного излучателя при h1 λ =3 показана на рис. 20.3. Число лепестков диаграммы направленности N легко определить, исходя из формулы (20.18)

273

Рис. 20.3

N = n max +1, n max ≤ 2λh1 .

Таким образом, чем больше относительная высота подъема излучателя h1 λ , тем больше интерференционных лепестков имеет диаграмма направленности.

2. По своим свойствам отражающая среда (земля) близка к диэлектри-

ку, т.е. ε'2 >> σ2 ω (КВ, УКВ над сушей). В этом случае ρ. В = f (θ). Напомним характер этой зависимости (рис.20.4 – идеальный диэлектрик, рис. 20.5 – ди-

электрик с потерями ρ. В = ρBe−jΦB ).

 

Рис. 20.4

 

Рис. 20.5

Здесь θ0 - угол Брюстера. Как следует из приведенных зависимостей

ρВ, ρВ =1

достигается только при θ → 900 ,

т.е. только для скользящих лучей.

При θ > θ0 ΦB = π, при θ < θ0

ΦB = 0 для идеальной отражающей поверхно-

сти (в

отсутствии потерь).

Поэтому

рассмотрим раздельно случаи

θ> θ0 и θ < θ0

1)θ >θ0 , ΦB ≈π, ρB ≤1,

274

Em =

60P0 Gm

F(θ)

1+ ρB

2 − 2ρB cos(2h1k1 cos θ),

 

r

 

 

 

I(θ)= 1+ ρB 2 − 2ρB cos(2h1k1 cos θ)

Проанализируем положения минимумов и максимумов интерференционного множителя I(θ).

Imax соответствует cos(2h1k1 cosθmax )= −1, т.е. cosθmax =

2n +1

 

λ

.

4

 

 

 

 

 

h1

Соответственно направления минимального излучения находятся из условия

cosθmin = n λ . 2 h1

При этом

Emax =

60P0Gm

F(θ

 

)[1+ ρ

 

(θ

 

 

 

 

 

max

B

max

)],

 

m

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(20.19)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

60P0Gm

F(θmin )[1−ρB (θmin )]

 

 

Emmin =

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Заметим, что для n = 0, θmin

= π 2 иρB (θmin )=1 и, следовательно,

Emmin = 0 .

2) θ < θ0 , ΦB ≈ 0, ρB <1,

 

 

 

 

 

 

 

 

 

Em =

60P0G m F(θ)

1+ ρB

2 + 2ρB cos(2h1k1 cos θ).

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Соответственно

 

 

 

 

 

 

 

 

 

 

 

cosθmax = n

λ

,

cos

θmin =

2n +1

λ

.

 

 

h1

 

 

2

h1

 

 

 

 

4

 

 

 

 

т.е. эти условия по сравнению с предыдущими противоположны.

Поля Emmin и Emmax описываются теми же формулами (20.19). Типичная диаграмма направленности вертикального диполя для рассматриваемого случая

 

ε21

>>

σ2

, ε1'

< ε2'

 

приведена на рис. 20.6.

 

 

 

ω

 

 

 

 

275

Z

θ

 

 

 

 

F(θ)

 

0

2

X

 

Рис. 20.6

 

20.4. Поле горизонтального вибратора, поднятого над плоской поверхностью земли на высоту h1

В рассматриваемом случае поле излучателя рассматривается в экваториальной плоскости и поэтому в расчетных формулах следует заменить θ на ϕ. Напомним одновременно, что излучение элементарного электрическо-

го излучателя в экваториальной плоскости изотропно, т.е. Em ≠ f (ϕ)и F(ϕ)=1. Таким образом, в рассматриваемом случае

Em =

60P0G m

1+ ρг

2 + 2ρг cos(2h1k1 cos ϕ + Φг ).

 

r

 

 

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

.

−

Φ

 

=

W0

cos ϕ − W0

cos ϑ

=

k

2 W20

cos ϕ − W10

k 22

− k12

sin 2

ϕ

.

ρг = ρге

j

 

Γ

2

1

cos ϑ

k

 

W0

cos ϕ + W0

k 2

− k 2

sin 2

ϕ

 

 

 

 

 

W20

cos ϕ + W10

 

2

 

 

 

 

 

 

 

 

 

 

 

2

1

2

1

 

 

 

Рассмотрим те же предельные случаи, что и в предыдущем разделе.

1. Отражающая среда близка по своим свойствам к проводнику

 

 

 

σ

 

 

 

 

 

 

 

 

 

W0

.

 

 

ε'2

<<

 

2

. В этом случае

 

2

→ 0 и ρг = −1, т.е. ρг =1,

Φг = π. Соответственно

 

 

W0

 

 

 

ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

Em =

60P0Gm

2sin(h

1k1 cos ϕ).

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

Направления максимального излучения, ϕmax

находятся из условия

 

 

 

cosϕmax =

2n +1

 

 

λ

,

n = 0, 1, 2,...

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

h1

 

 

 

 

Направления, в которых излучение отсутствует, определяются как

276

cosϕmin = n

 

λ

.

 

2

 

h1

Диаграмма направленности излучателя в рассматриваемом случае показана на рис. 20.7. Число лепестков диаграммы направленности определяется как

2h

1

 

1

 

N = nmax +1 ≤

 

+

 

.

λ

 

2

 

 

 

 

На рис. 20.7 изображен случай, когда h1 λ =2,3.

Рис. 20.7 2. Отражения происходят от поверхности среды, близкой по своим

свойствам к идеальному диэлектрику. В этом случае графики зависимостей ρг (ϕ)и Φг (ϕ) имеют вид, изображенный на рис. 20.8. В предельном случае, ко-

Рис. 20.8

гда σ2 ε2′ω → 0, ΦΓ =π . Поэтому

Em =

60P0G m

1+ ρг

2 − 2ρг cos(2h1k1 cos ϕ).

 

r

 

 

Направления минимального излучения определяются как

277

cos ϕmin

= n

 

λ , Emmin =

60P0 Gm [1−ρг (ϕmin )].

 

2

 

h1

r

Направления максимального излучения определяются следующим образом:

cos ϕmax = 2n4+1 hλ1

При n = 0 ϕmin = π2

, n = 0,1,2,...

Emmax =

60P0Gm [1+ ρг (ϕmax )].

 

 

r

 

и ρг =1,

 

π

= 0 .

поэтому Em 2

Типичная диаграмма направленности в рассматриваемом случае приве-

Рис. 20.9

дена на рис. 20.9. В отличие от случая вертикального вибратора возможно излучение в вертикальном направлении.

278

20.5.Распространение радиоволн над неоднородной

инегладкой отражающей поверхностью

Свойства отражающей поверхности на протяжении радиотрассы (в пространстве между передатчиком и приемником) обычно не остаются одинаковыми, т.е. отражающая поверхность обычно неоднородна, и поэтому возникает проблема определения ρ&β, ρ&r в интерференционных формулах

типа (19.25), (19.36). Заранее очевидно, что при отражении играют роль свойства ограниченного участка отражающей поверхности, границы которого определяются пересечением эллипсоида существенной области и отражающей

.

Рис. 20.11

.

Рис. 20.10

поверхности (рис. 20.10). Этот участок имеет вид вытянутого в направлении радиотрассы эллипса с поперечными размерами (рис. 20.10).

 

 

R = 1

1

 

a

s

(nλr)2

,

 

n

2

 

 

 

 

 

 

 

Rn - диаметр существенной области,

bs = cosasϕ .

В КВ диапазоне существенная область при отражении обычно в поперечнике занимает десятки метров, в продольном направлении – десятки километров ( h1 ≈ h2 ). На таком протяжении свойства отражающей поверхности могут меняться. Обычно используется среднее значение коэффициента отра-

−

жения ρ по существенной области. Значение ρ может быть определенно

279

экспериментально при изменении h1 или h2 с использованием формул типа

(20.19):

Emmin

=

1

−

ρ

;

ρ

=

1

− Emmin / Emmax

.

 

 

 

 

 

 

Emmax

 

1

+

ρ

1

+ Emmin / Emmax

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

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

Пусть неровность имеет высоту h над отражающей «средней» плоскостью. И

Рис. 20.12

пусть плоская волна падает на эту систему под углом θ . Луч 1 отражается от вершины неровности, луч 2 – от «средней» плоскости. Разность их хода в соответствии с чертежом (рис 20.12). ∆r = r2 − r1 = 2h cos θ, Фазовое запаздывание

луча 2 относительно луча 1, тогда составляет ∆ϕ = 2kh cosθ, k = 2π/ λ . Это запаздывание несущественно, если ∆ϕ < π/ 4 .

Это условие приводит к следующему критерию (критерию Релея):

h < 16cosλ θ

Поскольку в знаменателе справа стоит cosθ , то очевидно, что для скользящих лучей (θ →π 2) критерий Релея выполняется для весьма боль-

ших hλ . А этот случай, как указывалось выше, как раз и характерен для наземных радиотрасс, когда h1 ~ h2 .

280