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

Зайцев_Технмческие средства защиты информации

.pdf
Скачиваний:
2180
Добавлен:
01.06.2015
Размер:
11.03 Mб
Скачать

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

Рассмотрим некоторые особенности и свойства электромагнитных каналов.

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

Электромагнитное поле представляет собой особый вид материи. Оно, как и вещество, обладает не только энергией, но также массой, количеством движения и моментом количества движения. Поле может превращаться в вещество, как и вещество – в поле. Электромагнитное поле воздействует с определенной силой на заряженные частицы.

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

Для расчета электромагнитного поля наиболее пригодны уравнения электродинамики в интегральной и дифференциальной формах [35].

Электромагнитноеr поле характеризуется четырьмя векторнымиr величинами: E – напряженность электрического поля (В/м); D – электрическая индукция (вектор электрического смещения) ((Кл/м2 ); H – напряженность

магнитного поля (А/м); Br – магнитная индукция (Тл).

Определение поля в некоторой области пространства требует указания этих векторов в любой ее точке. В общем случае взаимосвязь векторов

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

 

D Er;

(1.1)

B H ,

(1.2)

где ε= εr ε0 – диэлектрическая проницаемость среды; ε0 =8,855 1012 – диэлектрическая проницаемость вакуума (Ф/м); εr – относительная диэлектрическая проницаемость среды, в которой находятся заряды; μ=μr μ0 – абсолютная магнитная проницаемость среды; μ0 =4π 107 – магнитная проницаемость вакуума (Гн/м); μr – относительная магнитная проницае-

мость среды.

Безразмерные величины εr и μr для воздушной среды близки к единице. Например, для воздушной среды при температуре 0° εr =1,0006.

28

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

r

r

дD

 

 

rot H =δ+

 

.

(1.3)

дt

 

 

 

 

Физический смысл этого уравнения можно толковать следующим об-

разом: магнитное поле возбуждается совместным действием тока проводи-

мости с плотностью δr и изменением во времени электрического поля (век-

тора электрического смещения D ). Величина дD называется плотностью

 

 

 

 

дt

 

тока смещения. Вектор δr указывает направление движения зарядов и по

абсолютному значению равен пределу

 

 

 

 

 

r

I

 

(1.4)

 

 

δ= lim

,

 

 

 

 

 

S0

S

 

где I – ток через площадку S , перпендикулярную δr

. Плотность тока

 

r

γE, где γ – удельная проводимость.

 

проводимости δ=r

 

r

дD

 

 

 

 

Сумму δ+

 

называют плотностью полного тока.

 

дt

 

 

 

 

 

 

Второеrуравнение Максвелла выражает скорость изменения магнитной

индукции B через пространственную производную (rot) напряженности

электрического поля Er:

дBr

 

 

r

.

(1.5)

rot E =−

дt

 

 

 

Физический смысл второго уравнения Максвелла состоит в том, что электрическое поле может возбуждаться не только электрическими зарядами, но иrизменениями во времени магнитного поля (вектора магнитной индукции B ).

Если изобразить в пространстве произвольную поверхностьr S с контуром L (рис. 1.4), то можно определить поток вектора rot E через эту по-

верхность.

Согласно (1.5) имеем:

r

r

=−

дB

r

 

(1.6)

rot EdS

 

dS.

 

S

r

S дt

 

 

 

Векторный символ

обозначает произведение элемента поверхности

dS

 

dS на единичный вектор нормали к ней n0 . r

где v

 

Применяя теорему Стокса ( rotvdV =vdl ,

– любой вектор) и

S L

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

29

где dlr вектор

r

r

 

д

 

r r

 

(1.7)

Edl

=−

 

BdS

,

 

L

 

 

дt S

 

 

 

– произведение элемента линии dl на касательный к ней единичный rτ0 .

Уравнение (1.7) представляет собой второе уравнение Максвелла в интегральной форме.

Если поверхность S (рис. 1.5) опирается на проводящий контур L (например, проволочный), то выражение (1.7) можно записать как

 

 

 

 

 

e =−дФ,

(1.8)

где

циркуляция

вектора

Er

дt

 

в этом случае

есть не что иное, как ЭДС

e =

r r

 

 

 

 

 

Edl , наводимая в контуре изменяющимся потоком вектора магнитной

 

L

д

r r

 

 

 

 

, где Ф – магнитный поток. В итоге для рас-

индукции, а

 

BdS =−

 

 

 

 

 

дt S

dt

 

 

сматриваемого случая имеем хорошо известный закон электромагнитной

индукции: e =−.

 

rot Er

dt

rotn Er

 

rotn Er

rot E

B

 

 

L

S

 

r

 

τ0

L

 

Рис. 1.4

Рис. 1.5

Второе уравнение Максвелла можно рассматривать как обобщенный закон электромагнитной индукции.

Интегральная форма первого уравнения Максвелла может быть получена интегрированием обеих частей уравнения (1.3) по произвольной поверхности S с контуром L и применением теоремы Стокса:

r

r

 

д

 

r r

r r

 

Hdl

=

 

DdS

+ δdS.

(1.9)

 

L

 

 

дt S

 

S

 

30

Интеграл δrdSr=I – поток вектора δ через поверхность S – является

 

 

 

S

 

током проводимости, пересекающим эту поверхность,

а составляющая

 

д

 

r r

 

 

 

DdS =Iсм – ток смещения. Сумма I +Iсм называется полным током.

 

 

 

дt S

К основным уравнениям Максвелла относят также следующие два

 

 

 

уравнения в дифференциальной форме:

 

 

 

 

divD ;

(1.10)

 

 

 

divB =0.

(1.11)

Согласно первому уравнению расходимость электрической индукции равна объемной плотности заряда ρ – величине, определяемой предельным

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

ρ= lim

q

,

(1.12)

V

V 0

 

 

где q – заряд, содержащийся в элементарном объеме

V .

Интегрированием обеих частей уравнения (1.10) по некоторому объему V и применением к левой части формулы Остроградского-Гаусса получим

DdS =q.

(1.13)

S

 

Здесь S – поверхность, ограничивающая объем V, а q = ρdV – полный

V

заряд в этом объеме.

Равенство (1.13) является интегральной формой уравнения Максвелла (1.10) и является формулировкой теоремы Гаусса: поток электрической индукции через замкнутую поверхность равен заключенному внутри ее заряду. rИнтегральную форму уравнения (1.11) получают интегрированием

divB по объему V и применением формулы Остроградского-Гаусса:

BdS =0.

(1.14)

S

 

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

Интегральная форма:

 

r

r

 

d

 

 

r

r

 

r

r

 

 

H dl

=

 

 

 

DdS +

δdS,

 

 

 

 

 

 

 

L

 

r

 

dt S

 

 

S

 

 

 

 

 

r

 

 

d

 

 

r

r

 

r

r

 

 

Edl

=−

 

 

 

 

BdS

,

DdS

=q,

(1.15)

dt

L

r

r

 

 

S

 

S

 

 

 

 

=0.

 

 

 

 

 

 

 

 

 

 

BdS

 

 

 

 

 

 

 

 

 

 

S

31

Дифференциальная форма:

r

=

дD

 

r

rot H

дt

 

,

 

 

r

 

 

 

 

 

 

r

=−

дB

,

rot E

дt

 

r

 

 

 

 

 

 

 

 

 

 

divD ,

 

 

 

divBr=0,

 

 

(1.16)

BrHr,

 

 

 

 

v

v

 

 

 

 

 

D E,

 

 

 

 

r

r

 

 

 

 

 

δ=γE.

Преобразованием (исключением D иBv) систему уравнений (1.16)

можно привести к форме, в которой переменными будут только напряженности электрического и магнитного полей:

r

 

 

дE

r

r

 

 

дHv

r

ρ

 

rot H

rε0

дt

, rotE

=−μrμ0

дt

, divE =

 

(εr =const),

εrε0

r

 

 

 

r

r

 

r

(1.17)

 

(μr =const),

 

 

 

 

divH =0

δ=γE

(при

Eстор =0).

 

Системы уравнений (1.15)…(1.17) являются исходными при изучении электромагнитного поля.

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

вить как ψ=ψm cos(ωt ψ) . Тогда всякий вектор поля V разлагается на

компонентыv , изменяющиеся по аналогичному закону:

V =av1V1m cos(ωt 1) +av2V2m cos(ωt 2 ) +av3V3m cos(ωt 3), (1.18)

где av1, av2 , av3 – орты некоторой системы координат q1, q2 , q3 .

Величина ω=2πf называется круговой частотой гармонических колебаний; ψm и Vim – амплитуды, ϕψ и ϕi – начальные фазы.

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

e j(ωt) =cos(ωt ) + jsin(ωt )

видно, что скаляр ψ (см. выше) и вектор V можно выразить как вещественные части величин

32

v

 

 

j(ωtψ )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&

 

 

 

 

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ=ψme

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.19)

v&

v

 

 

e

j

(ωt)

v

 

e

j(ωt

)

v

 

 

 

e

j(ωt

)

 

 

 

 

 

 

,

V

=a V

 

 

 

 

 

1

+a V

 

 

2

 

+a V

 

3

 

 

1 1m

 

 

 

 

 

 

 

2

2m

 

 

 

 

 

 

3

3m

 

 

 

 

которые называются их комплексами. В (1.19)

a1,ar2 , ar3 – орты некоторой

 

 

 

 

 

 

 

 

 

 

 

v

 

 

 

 

 

 

 

&

 

 

 

 

 

 

 

системы координат. Таким образом ψ=Re ψ, V =ReV.

 

 

 

 

 

 

Выделим в комплексе Vv&

 

 

 

 

 

 

 

 

&

 

 

 

 

 

 

 

 

 

 

 

множитель

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v&

 

v

 

 

e

jϕ

1

v

 

 

e

jϕ

 

v

 

 

e

jϕ

3

,

 

 

(1.20)

 

V

=a V

 

 

 

+a V

 

 

2 +a V

 

 

 

 

 

 

 

1 1m

 

 

 

 

2

2m

 

 

 

3

3m

 

 

 

 

 

 

 

который называют комплексной амплитудой. Через комплексную амплитуду можно выразить комплекс Vv& как Vv& =Vv&me jωt . Дифференцирование комплекса по времени соответствует его умножение на jω.

Если комплекс Vv& удовлетворяет некоторому линейному дифференциальному уравнению, то данному уравнению удовлетворяют его вещественная и мнимая части.

Сучетом приведенных выше соотношений уравнения Максвелла (1.17)

вкомплексных значениях принимают форму:

rot Hv&m v&m + jωεrε0Ev&m;

 

rot Ev&m =− jωμrμ0Hv&m;

 

div Hv&m =0;

(1.21)

div Ev&

=

ρ&m

;

 

 

 

m

 

εrε0

 

δv&m Ev&m.

Уравнения (1.21) могут быть упрощены, если учесть, что

εr =1,0006.

Рассмотрим некоторую область V (рис. 1.6), в которой распределен заряд (ρ≠0) и присутствует ток (δ≠0 ). В некоторой точке М существует

электрическое поле, потенциал которого φ есть решение уравнения Пуассона

д2ϕ

+

д2ϕ

 

+

д2ϕ

=−

ρ

дx2

дy2

дz2

ε

 

 

 

и выражается формулой

 

1

 

1

 

 

 

 

ϕ=

 

 

ρdV ,

 

 

 

 

 

 

 

 

4πεV r

 

 

 

 

а также магнитное поле, характеризуемое векторным потенциалом ределяемым из решения уравнения 2 A =−μδr как

(1.22)

(1.23)

A , оп-

33

Ar= μ 1δrdV , 4πV r

r д

r д

r д

 

где = x0

 

+ y0

 

+ z0

 

– оператор Гамильтона.

дx

дy

дz

 

 

 

 

 

 

r

 

 

 

 

 

 

M

0 δr

ρ

Рис. 1.6

Для решения системы уравнений (1.21) необходимо определить для электромагнитного поля электрический ϕ и магнитный A запаздывающие

потенциалы:

ϕ(t) =

1

 

 

1ρ(t

r

)dV ,

 

 

 

 

 

 

 

 

4πεV r

 

v

(1.24)

r

 

μ

 

 

1 r

r

 

 

 

 

 

 

 

 

 

 

A(t) =

 

 

 

δ(t

 

 

)dV ,

 

 

 

 

 

 

 

 

 

4πV r

 

v

 

 

 

где r – расстояние до точки наблюдения M ; v – фазовая скорость бегущей

волны, связанная с постоянной распространения волны в неограниченном

пространстве k соотношением

 

k =

 

ω

. Величины ρ и δr связаны между

 

 

 

собой уравнением

 

 

 

 

 

 

 

v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

dρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

divδ=−

 

 

 

 

.

 

 

(1.25)

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В комплексной форме выражения запаздывающих потенциалов при-

нимают вид:

 

 

 

 

 

 

 

 

 

 

 

 

 

&ejkr

 

 

&

 

 

=

 

 

1

 

 

 

dV ,

 

 

 

 

 

 

 

 

 

 

 

 

ϕm

 

 

 

 

 

 

 

ρ

 

r

 

 

 

 

 

 

 

4πεV

 

 

 

 

 

 

(1.26)

r

 

 

 

μ

 

 

 

 

 

r

 

jkr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&

 

 

 

 

 

 

 

 

&e

 

 

 

 

 

 

 

A

=

 

 

 

 

 

δ

 

 

r

 

 

dV.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4πV

 

 

 

 

 

 

 

 

 

Если рассматривать поле, создаваемое одним лишь колеблющимся за-

рядом q m V cosωt =qm cosωt,

 

расположенным в пространстве

V , то

согласно (1.26) комплексная амплитуда потенциала этого поля будет

 

 

 

 

 

 

 

 

 

 

&

 

 

 

 

e

jkr

 

 

ϕm =

 

qm

 

 

 

 

 

,

(1.27)

&

 

 

 

 

 

4πε

 

 

 

 

r

 

 

 

а сам потенциал равен:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

qm

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ=

 

 

 

 

 

cos(ωt kr).

(1.28)

 

4πεr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

34

=−μSrддIt .

В этом случае поле имеет форму сферической волны, расходящейся из точки, в которой расположен заряд, со скоростью v.

С учетом параметров А и ϕ напряженности магнитного и электриче-

ского полей можно выразить как

 

 

 

r

r

 

 

1

H =

 

 

rot A;

μμ0

 

 

(1.29)

 

 

 

r

r

=−

дA

gradϕ,

E

дt

 

 

 

 

где

dϕdx

gradϕ= dϕ

dy

dϕ

dz

1.2.1.2. Элементарный электрический излучатель

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

Электрический излучатель соответствует элементу электрического тока. В этом легко убедиться, если рассмотреть производную по времени от момента электрическогоr диполя. Так как электрический момент (векторная

величина) pr =qlr, то ддpt = ддqt lr=Ilr, при этом положительное направление тока I совпадает с p .

По аналогииr производная по времени от момента замкнутого витка с током rmr =−μI S магнитного диполя соответствует элементу магнитного

тока ддmt

Рассмотрим элементарный электрический излучатель. Для этого представим отрезок проводника l, ориентированный вдоль координатной оси z

ипо которому течет ток I =Im cosωt (рис. 1.7).

В[35] показано, что при условии постоянства амплитуды тока вдоль всего участка можно условно полагать сосредоточение равных по

35

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

 

 

 

 

qm

jIm

.

 

(1.30)

 

 

 

 

&

 

 

ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I

 

z

 

 

 

 

 

q

+ q

z

 

 

 

 

 

 

 

q

 

+ q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

l

 

 

Рис. 1.7

 

 

 

 

 

 

 

Рис. 1.8

 

Это значит, что рассматриваемый отрезок с током можно представить

как диполь, момент которого

pr&

= zr lq

совершает гармонические колеба-

 

 

 

m

 

0

 

 

 

 

 

 

ния с частотой ω и имеет комплексную амплитуду

 

 

 

 

 

pr&

 

=− j

Iml

zr .

 

(1.31)

 

 

 

 

 

 

 

 

 

m

 

 

ω

0

 

 

Изображенный на рис. 1.7 элемент тока (колеблющийся диполь) рассматривается в качестве элементарного излучателя и называется диполем Герца.

Расположив диполь в сферической системе координат (рис. 1.9) получают комплексную амплитуду векторного потенциала элемента тока:

 

r

 

=(rr

r

μlIm ejkr .

(1.32)

 

A&

 

cosϑ−ϑ sinϑ)

 

m

 

0

0

4πr

 

 

 

 

 

 

 

 

Z

 

 

 

 

 

 

 

 

ϑ

 

 

 

 

 

Z

 

 

 

 

 

 

 

 

 

 

r

Hα

 

r

 

M

 

 

ϑ

 

ϑ0

 

 

 

I

r

 

 

 

 

 

0

 

 

Eϑ

 

 

 

 

M

 

 

 

r

α

 

 

 

 

 

 

ϑ0

Рис. 1.9

Компоненты поля, создаваемого диполем Герца в произвольной точке пространства M (r,ϑ,α) , определяются по приведенным выше формулам и

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

36

Hα =

klIm

[

 

1

 

cos(ωt kr) sin(ωt kr)]sinϑ;

4πr

kr

 

 

 

 

 

 

 

 

 

E

=

 

klIm

 

 

 

[

 

1

sin(ωt kr) +cos(ωt kr)]cosϑ;

2πωεr2

 

r

 

 

 

kr

(1.33)

E

=

 

k2lIm

 

 

[(

 

1

1)sin(ωt kr) +

1

cos(ωt kr)]sinϑ;

 

 

 

 

k2r2

 

ϑ

 

 

4πωεr

 

 

 

 

 

kr

Hr =Hϑ =Eα =0.

Ближняя зона (зона квазистационарности). Границы этой зоны определяются условиями r >> l (l – длина элемента тока или плечо вибратора) и kr << 1, или r <<1/k. В силу равенства k = 2π/λ второе условие принимает вид r <<λ/2π (условие квазистационарности). Для ближней зоны (на расстояниях от вибратора существенно меньших длины волны) формулы (1.33) можно упростить, отбрасывая малые члены в квадратных скобках и пренебрегая фазовым сдвигом kr:

H

 

=

lIm

sinϑcosωt;

E

=

 

pm

cosϑsinωt;

 

4πr2

2πεr3

 

α

 

 

 

 

r

 

(1.34)

 

 

 

pm

 

 

 

 

 

lIm

 

E

 

=

 

sinϑsinωt;

p

=

.

 

 

4πεr3

 

 

 

ϑ

 

 

 

m

 

 

ω

 

Поле согласно (1.34) не имеет волнового характера, так как выражения (1.34) получены в пренебрежении излучением в ближней зоне вследствие его незначительности. Пространственное распределение в этом случае свойственно статическому диполю. Выражения (1.34) содержат одну составляющую вектора напряженности магнитного поля элемента тока и две составляющие вектора напряженности электрического поля вибратора, характеризующиеся в каждый момент времени как «стационарные» величины. Из (1.34) следует, что величины E и H сдвинуты по фазе на угол 90°.

Дальняя зона. Рассмотрим поле на расстояниях, значительно превышающих длину волны, когда r λ и kr 1 . В этом случае можно

пренебречь членами порядка 1k2r2 и 1kr . Тогда уравнения (1.33)

принимают вид [35]:

Hα =− klIm sinϑsin(ωt kr); 4πr

Er =0;

(1.35)

E

=−klW 0Im sinϑsin(ωt kr).

 

ϑ

4πr

 

 

 

37