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

Приведенные в [1—4] примеры иллюстрируют методику синтеза ОУ с архитектурой рис. 1, б.

На основании вышеизложенного можно сделать следующие выводы.

Для минимизации систематической составляющей напряжения смещения нуля операционного усилителя масштабные коэффициенты А:/|3, хп, кп4 и крЪ, хр, кр4 применяемых в них функциональных узлов следует выбирать в соответствии с равенством (7). Это позволит за счет введения слабой токовой ассиметрии уменьшить на один-два порядка статическую погрешность

ОУ, обусловленную конечной величиной р транзисторов и его радиационной (или температурной) зависимостью.

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

Статья подготовлена в рамках госконтракта П507 ФЦП "Научные и научно-педагогические кадры инновационной России на 2009—2013 годы".

СПИСОК ЛИТЕРАТУРЫ

1. Каскодный операционный усилитель с малым напряжением смещения нуля [Текст]: заявка на патент РФ, МПК8 Н 03 F 3/45, Н 03 F 3/34 / Прокопенко H.H., Серебряков А.И., Будяков П.С. № 2009119092/09; заявл. 20.05.2009 (186(B)).

2. Каскодный дифференциальный усилитель с малым напряжением смещения нуля [Текст]: заявка на патент РФ, МПК8 Н 03 F 3/45, 3/34 / Прокопенко H.H., Серебряков А.И., Будяков П.С. №

3. Каскодный дифференциальный усилитель с малым напряжением смещения нуля [Текст] : заявка на патент РФ, МПК8 H 03 F 3/45, 3/34 / Прокопенко H.H., Будяков П.С., Серебряков А.И. № 2009119952/09; заявл. 26.05.2009 (194).

4. Дифференциальный усилитель с малым напряжением смещения нуля [Текст]: заявка на патент РФ, МПК8 H 03 F 3/45, 3/34/ Прокопенко H.H.,

№

заявл. 28.05.2009 (188).

УДК621.39:574.5

В.В. Александров, Б.В. Александров

МЕТОДИКА ИЗМЕРЕНИЯ И АНАЛИЗА СПЕКТРА ЭНЕРГИИ ИЗЛУЧЕНИЯ ТЕЛЕФОННЫХ АППАРАТОВ СОТОВОЙ СВЯЗИ

Посвящается памяти профессора кафедры теоретических основ электротехники СПбГПУ докт. техн. наук И.Ф. Кузнецова

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

медицинские аспекты проблемы крайне обострены и далеки от разрешения [1,2].

Лаборатория электрофизики естественных экосистем Санкт-Петербургского государственного политехнического университета в течение десяти лет проводила экологическую оценку действия ЭМИ различных источников на магнитное поле среды обитания вблизи телефонных аппаратов (ТА) сотовой связи. Потребитель, как правило, используя ТА, близко подносит его к

височной области головы и слуховому нерву и создаются условия непосредственного влияния энергии излучения на нейрональную систему мозга. Попытки исключить или, по крайней мере, снизить влияние ЭМИ теоретически обоснованы в Фонде новых медицинских технологий "Айрэс", на основании которого созданы матричные преобразователи-аппликаторы самых различных модификаций, которые применяются для электронной аппаратуры. Такие преобразователи наклеиваются на изделие с тыльной стороны в районе антенны. Цель исследования:

выявить изменение естественного геомагнитного поля В2 в объёме пространства в ближней зоне ТА различных конструкций и модификаций;

оценить изменчивость В2 при использовании различных режимов функционирования ТА;

испытать матричные преобразователи-аппликаторы конструкции Фонда новых медицинских технологий "Айрэс";

определить пороги индукции ДВ2 при снижении влияния излучения для различных способов получения информации потребителем.

Для решения поставленных задач предварительно обследовались электромагнитные поля (ЭМП) помещения и применявшейся аппаратуры сотовой связи (ТА) и оптимально размещалась установка, включившая в себя: источники ЭМИ (ТА), аквариумы с рыбами (два опытных, один — контрольный), ПЭВМ, оснащенные \¥еЬ-камерами, а также комплект контроля переменных магнитного и электрических полей в лабо-

Рис 1. Преобразователь магнитной индукции ПМИ НЕВА-5 для регистрации многосуточного цикла данных и их картирования

ратории. В комплект входили: преобразователь магнитной индукции ПМИ НЕВА-5 (рис. 1) и измеритель переменного ЭП (НЕВА-Эн); измеритель квазипостоянного МП тесламетр НЕВА-4 (рис. 2); конструкции и аттестации Государственного научного центра ЦНИИ им. акад. АН. Крылова (фирмаЭТН, Санкт-Петербург) [1].

Краткие технические характеристики приборов представлены в [1].

Для проведения экспериментов использовалась операционная система Windows ХР на трех компьютерах следующего наполнения:

Intel Celeron764 МГц, 128 Мб, HDD - 40 ГГб, web-camera Creative;

Intel Celeron 1100 МГц, 128 Мб, HDD - 30 ГГб, web-camera Logitech;

Intel Celeron433 МГц, 128 Мб, HDD - 4 ГГб, CDRW40.

Микроструктура квазипостоянного ЭМИ ТА

Оценка состояния магнитного поля в аудиториях СПбГПУ. Для работы выбиралась горизонтальная сетка 0,5x0,5м и двауровня по вертикали: на высоте 0,7 и 1,5 м от пола. Характерные примеры диаграмм распределения уровня магнитного поля Bz в помещении представлены на рис. 3 и 4 (опыты 2008 года).

Магнитным фоном для исследуемого помещения является величина48 мкТл (± 2,5 %).

Значения вектора магнитной индукции изменяются (мкТл) в пределах:

с включенным электропитанием 20,8—86,4 (0,7 м над полом) и 23,4—56,0 (1,5 мнад полом);

Á

Рис. 2. Тесламетр НЕВА-4 для измерения постоянной и переменной составляющих вектора индукции магнитного поля В2

а)

В, мктл

>>,м

X, м

б)

Д мктл

>>,м

В, мктл:

0-5  5—10  10-15  15-20 В20-25 «25-30

30-35  35—40 П40-45 П45-50 И50-55 И55-60

60-65 "65-70 "70-75 "75-80 "80-85 Ш5-90

х, м

Рис. 3. Распределение магнитного поля в аудитории на уровне 0,7 м от пола с включенным электропитанием

а)

В, мктл

б)

>>,м

В, мктл

5 X, М

В, мктл:

Я 0-5 "5-10 И 10-15 Ш 15-20  20-25  25-30 Я 30-35 В 35-40  40-45 О 45-50 в 50-55  55-60

Рис. 4. Распределение магнитного поля в аудитории на уровне 1,5 м от пола без электропитания

у, м

5 X, М

без электропитания 23,5-57,2 (0,7 м над по- го помещения, определяли источники магнитных

лом) и 29,1—56,8 (1,5 м над полом). аномалий.

Сопоставляя диаграммы распределения Электрическим фономдля лаборатории можно

уровня магнитного поля в помещении и план это- условно принять напряженность 0,005-0,01 кВ/м

при отключенном и включенном электропитании. Значения вектора напряженности ЭП Е2 (кВ/м) изменяются в пределах:

при отключенном электропитании 0,005— 0,050 (0,7 м) и 0,005-0,100 (1,5 м);

при включенном электропитании 0,005— 0,200 (0,7 м) и 0,005-0,160 (1,5 м).

Таким образом, на основании тщательно измеренных значений В2 и Е2 и исходя из реальных возможностей помещения были выбраны площадки (на столах) для размещения предполагаемых биологических объектов (экосистем), действующих источников облучения (телефоны) и \¥еЬ-камеры регистраторов.

Оценка состояния индукции В/ в пределах ра -бочей зоны. Рабочая площадка разделяется сеткой с шагом 5x5 см (5"= 3600 см2). В каждом узле сетки получают значение вертикальной составляющей В2.

Квазистационарное магнитное поле В2 изменялось, как правило, в пределах 40—65 мкТл. Для Северо-Западного региона фоновое значение составляет 48 мкТл. Получив карту распределе-

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

Пространственная модель микроструктуры Е2 и В2 в ближней зоне телефонного аппарата сотовой связи. В пределах рабочей зоны выявляют площадку с фоновыми значениями магнитной индукции. Её разделяют сеткой с шагом 1x1 см. На площадку помещают, например, сотовый телефон Nokia 3110 с нейтрализатором AIRES и изучают структуру электрического и магнитного полей в ближней зоне сотового телефона в стационарном режиме.

Очевидно, что рассчитать пространственную структуру изменения таких полей, учитывая массу побочных источников, весьма сложно, а моделирование также даст лишь идеализированную картину. Поэтому было произведено измерение электрического и магнитного полей вблизи ТА в режиме звонка на реальных расстояниях от оси ориентации приборов (рис. 5 и 6).

а) Д мкТл ]————————————————

05 —......................

........... .........................—........................... ..........

0.3 ...................—.........................-........ .....-

0,2------..... м—-......--+--..........^——лу---^—

°л I Гтт"Т"[ i 1 i i i м м

О 10 20 ЭО 40 50 60 70 S0 90 100 110 120 130 140 £, ММ

б) В, мкТл Чг—I-!-!-!-!-\-\-\-\-\-\-\-\-j-j-

\ : : : : : :::::: : : : :

0 5 JA-J-----:-----L—J-----:.....—----J—-A—X—4—4-----

1 \ M M M M M / I I i

игн---------;-----r — 1 — 1.....r —-----:-----r—1 — —* —r\—

•\ : : : : : : / : - : :. ; j \

u......N-f-i-l.....Ы \ if bfi----m\

^.........\ * * v —4-—j-—i—-i"*--;—--i.........

.....rrrrri N • •

J 1 i i 1 1 i i i i i i i i 1 I ^

0 10 20 30 40 50 60 70 80 SO 100 110 120 130 140 £, MM

i-9 Г*" N ;-з ;-з N V

—' N

: : : : : : : : : : : :

! ! ! '/..... '+"1.....

.....L____J_____I..... j : : i i J \

! ! !    ! ! J : : A \

г_...[~~]........4 r r ^ 4

: : * f f 1 j j 1 1 1

Рис. 5. График изменения напряжённости электрического поля вдоль сотового телефона с нейтрализатором

Отличительная черта объёмных структур — убедительная разница их для ТА в обычном режиме работы и для ТА с нейтрализатором: его наличие обусловливает двухвершинное распределение (рис. 7) электрического и магнитного полей, причем максимальная амплитуда Е2 примерно в два раза меньше, чем для ТА без нейтрализатора, хотя амплитуды по В2 весьма близки (рис. 6).

Структура переменной компоненты В2. При

изучении процесса влияния ЭМИ сотового телефона на среду обитания мы столкнулись с проблемой электромагнитной совместимости технических средств, биологических объектов и изменчивости геофизических факторов (в частности, ГМП). Выяснилась очень интересная подробность при анализе и интерпретации колебаний МП в лаборатории, которые мы регистрировали как фоновую характеристику в реальном амплитудно-временном масштабе. Тем более это стало важным теперь, когда нужны нюансы при оценке слабых воздействий. На рис. 8 представлены реализации колебаний МП в лаборатории, характеризующие помехи, обусловленные работой ТА. Телефонные аппараты помещали возле магнитометрического датчика и точно отмечали время воздействия ЭМИ в режиме ожидания или

в режиме звонка. Применяли попеременно аппараты без "нейтрализатора" и с "нейтрализатором" в течение 400 с. В режиме звонка ЭМИ ТА без "нейтрализатора" индуцирует нарастающий сигнал, который в итоге смещает средний условный фоновый уровень магнитных колебаний (рис. 8, б участок а) почти в шесть раз в направлении "+", при этом амплитуда колебаний также возрастает примерно в два раза. Самым неприятным здесь оказывается последующее состояние датчика МП, у которого отключается автокомпенсация постоянной составляющей компоненты Bz, и условный фоновый уровень становится новым. При автоматической компьютерной регистрации МП это может привести к появлению неконтролируемых ошибок. При использовании ТА с "нейтрализатором" картина вполне удовлетворительная: соответствующий участок звонка (рис. 8, в, участок Р) выделяется также ростом амплитуды колебаний примерно в 1,5—2 раза, но условный фоновый уровень не смещается и автокомпенсация датчика не нарушается. Таким образом достигается гармонизация в электромагнитной совместимости между аппаратурой и естественным фоном среды обитания. Мы полагаем, что это является серьёзным достижением защитного устройства AIRES.

Ez, kB/m

0,28-0,3

0,26-0,28

0,240,26

0,22-0,24

0,2-0,22  0,18-0,2  0,16-0,18  0,140,16

0,12-0,14

0,1-0,12  0,08-0,1  0,06-0,08

0,С40,06

0,02-0,М

0-0,02

0,000-0,020  0,020-0,040  0,040-0,060  0,060-0,080  0,080-0,100

0,100-0,120 то, 120-0,140 И0,140-0,160 И0,160-0,180 И0,180-0,200  0,200-0,220 И0,220-0,240 «0,240-0,260 «0,260-0,280 «0,280-0,300

Рис. 6. Электрическое поле вокруг сотового телефона Nokia 3110 без нейтрализатора во время звонка

Рис. 7. Электрическое (а) и магнитное (б) поля вокруг сотового телефона Nokia 3110

с нейтрализатором во время звонка

Нужны глубокие экспериментальные проработки и статистические данные для получения окончательных рекомендаций и оценок.

Спектральная характеристика ЭМИ сотового ТА. Преобразователь магнитной индукции "Нева-5" регистрирует флуктуационную часть ДВ2 при компенсации постоянного геомагнитного поля ГМП, и частотный спектр ЭМИ заведомо не известен. Следовательно, если регистрация идет в широком спектре частот, причём наша методика допускает разрешение 0,001 с, результирующий эффект будет зависеть именно от этого, и диапазон анализируемого излучения будет в пределах 0—500 Гц. Таким образом, при излучении в диапазоне 900/1800 МГц мы заведомо фильтруем высокочастотную часть спектра и получаем оценочные результаты.

Измерения проводились в четырёх режимах функционирования сотового телефона (ожидание, звонок, разговор и приём текстового сообщения).

Чтобы получить объективную информацию и исключить ошибки, измерения проводились пять дней для двух телефонов по три раза, и каждое измерение длилось по 4 мин (телефонный аппарат находился 1 мин в режиме ожидания, 1 мин шел звонок, 1 мин длился разговор и снова шла 1 мин ожидания). Приём текстовых сообщений проводился с начала второй минуты эксперимента. Измерения проводились на нулевом расстоянии от ПМИ "Нева-5".

Для каждого измерения была построена спектрограмма Gif, /), где /—частота флуктуации ABZ, регистрируемая ПМИ "Нева-5"; t— время,

à)

с 20 40 so so ¡o i: и ií ib 20 22 25 2S 30 :-2 ¿4 зе за 40

OÜOOOOOOOOOOOOOO

б)

-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-r~

0 20 40 60 SO 10 12 14 15 IS 20 22 24 26 2B 30 32 34 Jó 3S 40 CQOOOOOQOOOOOOOO

e)

jß

«ULluiJil, lU i111 vli) J

0 30 40 Я1 30 100 (20 140 1Í0 183 200 !2fl 3+0 260 230 ЭОС 320 340 МО JSO 400

Рис. 8. График изменения вертикальной составляющей переменного магнитного поля около сотового телефона в режиме ожидания

в течение которого длился процесс ЭМИ от источника; Gif, г) — функция спектральной плотности энергии магнитного поля. Затем эти спектры усреднили и построили общую спектрограмму.

Обработка реализаций проводилась с помощью пакетов анализа Microsoft Excel. Каждая запись разбивалась на одинаковые временные отрезки для вычисления спектров на основе метода быстрого преобразования Фурье. Полученные спектры собирали вместе, и для этого массива производилось построение временных спектрограмм (рис. 9, а, б, в):

по горизонтали (осьх) откладывалось время t (0-240 с), ? (0-120 с);

по вертикали (ось .у) откладывалась частота/ (0-500 Гц);

по ортогонали к плоскости чертежа в виде градаций цвета откладывалась функция спектральной плотности Gif, t).

Результаты измерений представлены на спектрограммах флуктуаций плотности энергии магнитного поля (рис. 9) для разных ТА и в режиме приёма сообщений.

На рис. 10 представлены спектральные функции флуктуирующего магнитного поля, излучаемого аппаратами сотовой связи в течение 0—2600 с для различных режимов работы. Как видно, в зависимости от режима (ожидание или приём текстового сообщения) минимум спектра смещается в более высокочастотную область излучения.

Кроме того, представлены спектрограммы плотности энергии В2 в иных условиях измерения и для других ТА (рис. 11) с учётом защитных устройств Фонда "Айрэс".

Известно, что тонкие измерения, которые проводились нами на протяжении многих лет, могут быть связаны с погрешностями. Поэтому

Рис. 9. Микроструктура спектральной плотности энергии магнитного поля Gif, t) в ближней зоне технического устройства сотовой телефонной связи: a — Samsung Duos D-880

весьма необходима достоверность, обусловленная геомагнитной обстановкой в Северо-Западном регионе.

С этой целью выбираются особые дни, когда космическая погода благоприятна. В частности, измерения в 2008 году проводились в течение пяти дней: 01.10.08, 03.10.08, 06.10.08, 09.10.08, 13.10.08 (рис. 12). В эти дни геомагнитная обстановка в Северо-Западном регионе была спокойной (только 03.10.08 наблюдались геомагнитные возмущенная). При этом по данным ИЗМИРАН (Институт земного магнетизма и распространения радиоволн Российской академии наук), пятнообразователь-ная активность Солнца прошедшего периода была на очень низком уровне. Вспышечная активность была на очень низком уровне весь

период. В остальные дни геомагнитная обстановка была спокойной и слабо возмущенной [www. izmiran. ru].

Экспериментальными исследованиями, проведенными в лаборатории электрофизики водных экосистем Учебно-научного центра "Мониторинг и реабилитация природных систем" СПбГПУ установлено, что кроме основного несущего сигнала на частотах 900/1800 МГц ТА в режимах звонка и разговора генерирует переменное электрическое поле в диапазоне 5—2000 Гц и переменное магнитное поле в диапазоне 5—500 Гц. Эти поля фиксируются аппаратными методами в ближней зоне непосредственно над поверхностью корпуса ТА на расстоянии нескольких миллиметров (в диапазоне 0—30 мм) и имеют сложное распределение в пространстве по корпусу

1 UlulAkK. jk . НА

1 ООО 1500 „ ¡ПС 2Я

.

-1 V- ... jJ I_:_|.

РРЦЛ1 '

t : 00 400 60 о\ 8 QÜJ «fflWr >00

Рис. 10. Спектрограммы флуктуации магнитной индукции Bz во времени для телефонных аппаратов сотовой связи а — Samsung duos D-880 в режимах функционирования

телефона, как это представлено на соответствующем графике (см. рис. 5).

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

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

а)

0,001

0,0001

О 50 100 150 200 250 300 350 400 450 J] Гц

б)

0,001

0,0001

О 50 100 150 200 250 300 350 400 450 /, Гц

Рис. 11. Микроструктура спектральной плотности энергии Bz в режиме ожидания

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

Как показывают эксперименты, проведенные в СПбГПУ, применение нейтрализатора электромагнитных аномалий "Айрэс" позволяет изменить как распределение

низкочастотных электрических и магнитных полей по всей протяженности ТА, так и величины напряженностей этих полей, реструктурируя их, что особенно важно для ликвидации привнесенного низкочастотными акустическими колебаниями конфликта в цитострук-туру биоформы.

X ** 2- ¿1$-Э в -п й-и ЯП 14-Н

^ % |у $ 3 2> 3 & ЯА

2 ъ н ч ? % Я. щ

3 % я % } ч Н 1 ш

5 % 3 ъ 3 а > №

л & & 3 ь X г'Н 15

л й, г А <5 Ш

7 & л £ Й- г £

Л 0 0 0 0 о 0 0 л

Р а 0 0 -V % 1 ? 0 л 6

а 0 1 ь ь 1 & Т £ X ь н

<1 1 н ъ 5 с ж

и 13 ж $ з: 3 % ъ Ъ н ь

л * ^ . и £ $

/4 15 1 ъ £ ъ ъ 1 5 X а Я А, 3 А 1 4 №

Рис. 12. Значения /¿-индексов Трехчасовая характеристика цифрового вида геомагнитной обстановки СЗ региона, данные ИЗМИРАНс 1.10.08 по 15.10.08

Основные частоты мозговых колебаний, биоритмы головного мозга человека: 1,5—2, 3,7 и 4,7 Гц и т. д. по своим частотным характеристикам находятся непосредственно в зоне генерации зафиксированных В.В. Александровым с сотрудниками сопутствующих паразитных излучений ТА (от 5 до 500 Гц). Таким образом, отсутствие корреляции между естественными для мозга человека колебаниями и генерируемыми техногенным источником (в данном случае ТА) провоцирует конфликт, проявляющийся на уровне работы различных структур мозга.

Кроме рассмотренных результатов, связанных с искусственным ЭМИ, следует обращать особое внимание на естественное эволюционное поведение ГМП, аппроксимированное поведением гидромагнитного динамо.

Аппроксимация геомагнитного поля плоским геоцентрическим диполем была постулирована К. Гауссом, но современное представление о происхождении ГМП и его свойствах, изменяющихся на протяжении миллиардов лет, обосновано сравнительно недавно, оно базируется на постулировании эффекта самовозбуждения магнитных полей вследствие движения проводящей жидкости или газовой плазмы — динамо-эффекта [3].

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

Теория динамо-эффекта приводит также к возможности самообращения магнитной оси (переполюсовке магнитного поля Земли) и дол-гопериодическим колебаниям ГМП (вековым вариациям), что отражает реальные свойства земного магнитного поля.

Наше мнение заостряет внимание специалистов на том, что биообъекты, адаптированные в течение миллиардов лет ктаким свойствам ГМП, в современных условиях попадают в поле воздействия ЭМИ высоких частот и теряют естественные физиологические свойства как целостных организмов, так и его частей. Этим вполне объяснима угнетающая роль ЭМИ ТА при экспериментах [1,2].

Проблемы экологической безопасности для промышленного Северо-Западного региона России требуют законодательных решений.

Авторы выражают признательность коллегам, студентам, инженерам, внесшим большой вклад в полученные результаты: В.М. Данилову, И.Н. Серо-

ву, ДА. Уставникову, А.Я. Сергееву, Ю.М. Балагуле, Е.А. Кичигину, М.А. Тришину, Э.А. Ярошевичу, ТВ. Акинчиц, В А Тартыжеву, М.Ю. Александровой.

СПИСОК ЛИТЕРАТУРЫ

1. Александров В.В. Экологическая роль элект- 2. BioEM2009. Davos, Switzerland. Technical

ромагнетизма. СПб.: Изд-во Политехи, ун-та. program. 2009. June 14—19. 58 p. 2006. 716 с. 3. Физический энциклопедический словарь.

М.: Сов. энцикл. 1983. С. 161.