12
ТЕХНИЧЕСКИЕ НАУКИ
УДК 004.056(075.8)
А. А. Шелупанов, А. П. Зайцев
Оценка ПЭМИ электронных устройств
Рассмотрена физическая природа побочных электромагнитных излучений электронных устройств, приведены соотношения для расчета электрической и магнитной составляющих излучаемого электромагнитного поля. Указано на целесообразность применения для исследований автоматизированных программно-аппаратных комплексов.
Основные уравнения электромагнитного поля. Электромагнитное поле определяется во всех точках двумя векторными величинами — электрическим и магнитным полем.
Для расчета электромагнитного поля наиболее пригодны уравнения электродинамики в интегральной и дифференциальной формах [1].
Электромагнитное поле характеризуют четыре векторные величины:
E — напряженность электрического поля (В/м);
D — электрическая индукция (вектор электрического смещения) ((Кл/м2);
H — напряженность магнитного поля (А/м);
B — магнитная индукция (Тл).
Основными уравнениями электромагнитного поля являются уравнения Максвелла. Первое уравнение Максвелла соответствует вихрям магнитного поля:
rot H = 5 + ddD. (1)
dt
Физический смысл этого уравнения можно толковать следующим образом: магнитное поле возбуждается совместным действием тока проводимости с плотностью 5 и изменением во времени электрического поля (вектора электрического смещения). Вектор 5 указывает направление движения зарядов и по абсолютному значению равен пределу
5 = lim —, (2)
AS ® 0 AS
где AI — ток через площадку AS , перпендикулярную 5. Плотность тока проводимости 5 = gE, где g — удельная проводимость.
Второе уравнение Максвелла выражает скорость изменения магнитной индукции B через пространственную производную (rot) напряженности электрического поля E :
rot E = - Щ. (3)
dt
Физический смысл второго уравнения Максвелла состоит в том, что электрическое поле может возбуждаться не только электрическими зарядами, но и изменениями во времени магнитного поля (вектора магнитной индукции B ).
Второе уравнение Максвелла можно рассматривать как обобщенный закон электромагнитной индукции.
К основным уравнениям Максвелла относят также следующие два уравнения в дифференциальной форме:
div D = р; (4)
div B = 0.
Согласно первому уравнению расходимость электрической индукции равна объемной плотности заряда р — величине, определяемой предельным соотношением
Р= lim Aq, (5)
AV®0 AV
где Aq — заряд, содержащийся в элементарном объеме A V .
Исходные уравнения путем преобразований приводятся к дифференциальной форме записи:
4. й dD 5 rot H = —— + S, dt
*5 oB
rot E =--,
dt
div D = p,
div B = 0, (6)
в = mH, D = eE, S = gE.
Дальнейшим преобразованием исходную систему уравнений (6) можно привести к форме, в которой переменными будут только напряженности электрического и магнитного полей:
rot H = ereo dE + S, dt
* 5 dH
rot e = -mrmo -of'
div E = p (er = const), (7)
er e0
div H = 0 (mr = const),
S = g E (при Естор = 0).
Система уравнений (7) может являться исходной при изучении и расчетах электромагнитного поля.
Элементарный электрический излучатель. Диполь, момент которого изменяется во времени, называют элементарным излучателем. Различают электрический и магнитный излучатели: электрический и магнитный диполи. Диполь, момент которого изменяется по синусоидальному закону, называют гармоническим.
Рассмотрим элементарный электрический излучатель. Для этого представим отрезок проводника l, ориентированный вдоль координатной оси z, по которому течет ток I = Im cos ю t (рис. 1).
В [1] показано, что при условии постоянства амплитуды тока вдоль всего участка можно условно полагать сосредоточение равных по абсолютной величине и противоположных по знаку колеблющихся зарядов (рис. 2) с комплексными амплитудами
qm = ± ^. (8)
ю
Это значит, что рассматриваемый отрезок с током можно представить как диполь, момент которого Pm = Zolq совершает гармонические колебания с частотой ю и имеет комплексную амплитуду
Pm =-JImlZo. (9)
ю
+ q
q
" i ' i
Рис. 1 Рис. 2
q
I
z
z
q
Изображенный на рис. 1 элемент тока (колеблющийся диполь) рассматривается в качестве элементарного излучателя.
Компоненты поля, создаваемого диполем Герца в произвольной точке пространства M(r,J, a) (рис. 3), определяются по приведенным выше формулам и при переходе от комплексов к векторам поля принимают вид:
Ha = kkIm [-Lcos(m£ - kr) - sin(rnt - kr)]sin J;
4pr kr
Er = k1Im [^sin(a)t - kr) + cos(rnt - kr)]cos J;
2pwer2 kr (10)
k2lI 1 1
Ej =-m[( „ 0 - 1)sin(rnt - kr) + — cos(rnt - kr)]sin J;
4гов£ r k2r2 kr
Hr = HJ = Ea=
Ближняя зона (зона квазистационарности). Границы этой зоны определяются условиями r>>l (l — длина элемента тока или плечо вибратора) и kr<<1, или r<<1/k. В силу равенства k=2n/Á второе условие принимает вид r<<Á/2n (условие квазистационарности). Для ближней зоны (на расстояниях от вибратора существенно меньших длины волны) формулы (10) можно упростить, отбрасывая малые члены в квадратных скобках и пренебрегая фазовым сдвигом kr:
Ha = l1"2 sin J cos wt; Er = Pm 3 cos J sin wt;
4nr 2per
II (11)
Ej = Pm ч sin J sin wt; pm = —m.
4яег3 w
Поле согласно (11) не имеет волнового характера, так как выражения (11) получены в пренебрежении излучением в ближней зоне вследствие его незначительности. Пространственное распределение в этом случае свойственно статическому диполю. Выражения (11) содержат одну составляющую вектора напряженности магнитного поля элемента тока и две составляющие вектора напряженности электрического поля вибратора, характеризующиеся в каждый момент времени как «стационарные» величины. Из (11) следует, что величины E и H сдвинуты по фазе на угол 90°.
Дальняя зона. Рассмотрим поле на расстояниях, значительно превышающих длину волны, когда
r >> l и kr >> 1 . В этом случае можно пренебречь членами порядка 1 k2 r2 и 1/kr . Тогда уравнения (10) принимают вид [1]:
klI
Ha =--m sin J sin(wt - kr);
4 яг
klwOl (12)
Er = 0; Ej =--m sin Jsin(wt - kr).
4Я r
В (12) введено отношение амплитуд Em и Hm , которое равно W0 = = J^ и называется
Hm ' e
волновым сопротивлением неограниченной среды. Для вакуума W0 = í— = 120я [Ом].
\£0
Уравнения (12) соответствуют полю излучения. Оно представляет собой сферическую волну. Векторы E и H расположены перпендикулярно к направлению распространения волны, взаимно перпендикулярны и синфазны. Излучение максимально в экваториальной плоскости (J= 900) и отсутствует в осевом направлении (J = 0).
Элементарный магнитный излучатель. В теории электромагнитного поля доказывается [35], что замкнутый виток (рис. 4, а) с постоянным током на превышающих его размеры расстояниях создает такое же магнитное поле как если бы на его месте находился магнитный диполь (рис. 4, б) с моментом m = ¿^ImS.
При гармоническом токе витка I = Im cos wt переменный магнитный диполь характеризуется комплексной амплитудой момента m = Z0Imm S . Такой виток называют элементарным магнитным излучателем или магнитным диполем Герца.
Решение уравнений Максвелла для магнитного диполя Герца в комплексной форме имеет вид [ 1 ]
E = -a0 ImS (! + jk)ej(wt-kr) sin J;
4pr r (13)
H = mR,2(i + jk) cos J + J0 ^(-i- + jk - k2)sin J]ej(wt-kr).
4я r2 r r r2 r
{
f + q m
- q'
Рис. 4
Из (13) определяется запись компонент электромагнитного поля:
I k2SW0 1
m [— sin(wt - kr) + cos(wt - kr)] sin J;
E = —
4pr L kr
Hr = ImkS [-L cos(wt - kr) - sin( wt - kr)] cos J;
r 2nr2 lkr Imk2Sr, 1
HJ = / lv,2 2 4ppr k2r2
1
m-[( 0 0 - 1)cos(wt - kr) -—— sin(wt - kr)] sin J;
kr
Er = EJ = Ha= Из (14) получаем компоненты ближнего поля:
Ea = ImmS<W sin J sin wt; Hr = mm ^ cos J cos wt;
4pr2 m„
r 2Щг 3
Hj = ' m 3 cos J cos wt; mm = ImmS 4яцг 3
(14)
(15)
z
б
а
и поля излучения:
I k2SW
Ea = —-cos(w t - kr) sin J; Hr = 0;
4яг
I k2S
HJ = -cos(w t - kr) sin J. (16)
4яг
В дальней зоне элементарный магнитный излучатель создает волновое поле, которое отличается от поля элементарного электрического излучателя только ориентацией (рис. 5). Диаграмма направленности магнитного излучателя не отличается от диаграммы направленности элементарного электрического излучателя.
Рассмотренные соотношения могут служить основой для разработки программ расчета компонентов электромагнитного поля, создаваемого побочными электромагнитными излучениями электронных устройств в любой точке пространства.
В компьютерных системах исследованию на ПЭМИ подлежат следующие устройства [2]:
• видеоподсистема;
• накопители на жестких и гибких дисках;
• устройства CD, CD-R, CD-RW, DVD, DVD-RW;
16
ТЕХНИЧЕСКИЕ НАУКИ
• клавиатура;
• последовательные порты;
• принтеры.
Применяемые средства исследования могут быть различными — ручными или автоматизированными. Укажем на целесообразность применения автоматизированных программно-аппаратных комплексов, таких как «Легенда» и «Сигурд», выполненных на базе анализаторов спектра «К&З» и «№К».
Эти комплексы имеют управляющую и расчетную программы, способны автоматически опознавать опасные сигналы по форме их огибающих, заданных тест-программами. Оба комплекса в результате специальных исследований определяют опасные зоны ^ , Г и т{ и формируют отчетный протокол. Эти исследования могут быть дополнены исследованиями по методу реальных зон.
Упрощенная схема измерения ПЭМИ программно-аппаратным комплексом «Легенда» показана на рис. 6.
Исследуемое ТС Анализатор
Рис. 6
Внешний вид основного рабочего экрана управляющей программы «Легенда» представлен на рис. 7.
Рис. 7
На рис. 7 цифрами обозначены следующие элементы основного рабочего экрана (главного окна):
1 — окно спектра (отображается спектр сигнала для выбранного диапазона частот);
2 — таблица исследований (отображаются все программы исследований, которые будут выполнены в автоматическом режиме);
3 — окно эталона (осциллограмма найденного в полуавтоматическом режиме эталона тестового сигнала);
4 — окно сигнала (осциллограмма сигнала, который будет сравниваться с эталоном);
5 — рабочая таблица промежуточных результатов.
Литература
1. Никольский В.В. Теория электромагнитного поля. Учеб. пособие. — М.: Высшая школа, 1964. — 384 с.
2. Бузов Г.А., Калинин С.В., Кондратьев А.В. Защита от утечки информации по техническим каналам: Учеб. пособие. — М.: Горячая линия — Телеком, 2005. — 416 с.
Шелупанов Александр Александрович
ГОУ ВПО «Томский государственный университет систем управления и радиоэлектроники»,
зав. кафедрой комплексной информационной безопасности электронно-вычислительных систем, д.т.н., профессор. Эл. адрес: [email protected].
Зайцев Александр Петрович
ГОУ ВПО «Томский государственный университет систем управления и радиоэлектроники», к.т.н., профессор кафедры комплексной информационной безопасности электронно-вычислительных систем. Эл. адрес: [email protected].