автореферат диссертации по приборостроению, метрологии и информационно-измерительным приборам и системам, 05.11.17, диссертация на тему:Разработка метода выделения влияющих на организм человека возмущений геомагнитного поля для создания информационно-измерительных средств медицинского назначения

На правах рукописи

Алексанян Грайр Кареыович

РАЗРАБОТКА МЕТОДА ВЫДЕЛЕНИЯ ВЛИЯЮЩИХ НА ОРГАНИЗМ ЧЕЛОВЕКА ВОЗМУЩЕНИЙ ГЕОМАГНИТНОГО ПОЛЯ ДЛЯ СОЗДАНИЯ ИНФОРМАЦИОННО-ИЗМЕРИТЕЛЬНЫХ СРЕДСТВ МЕДИЦИНСКОГО НАЗНАЧЕНИЯ

Специальность

05.11.17 - Приборы, системы и изделия медицинского назначения

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

0.4 ОКТ 2012

Новочеркасск 2012

005052478

Работа выполнена на кафедре «Информационные и измерительные системы и технологии» ФБГОУ ВПО «Южно-Российский государственный технический университет (Новочеркасский политехнический институт)»

Научный руководитель: - доктор технических наук,

профессор Н.И. Горбатенко (ЮРГТУ (НПИ), г.Новочеркасск)

Официальные оппоненты:

- доктор технических наук, профессор H.H. Чернов (ТТИ ЮФУ, г.Таганрог)

- кандидат технических наук, доцент A.B. Литвин (ДГТУ, г.Ростов-на-Дону)

Ведущее предприятие: Федеральное государственное

унитарное предприятие «Научно-исследовательский институт специальных

информационно-измерительных систем», г. Ростов-на-Дону

Защита состоится « 27 » сентября 2012 г, в 1420 ч. на заседании диссертационного Совета Д212.208.23 при Южном федеральном университете по адресу: 347928, г. Таганрог, ул.Шевченко, 2., ауд.Е-306.

С диссертацией можно ознакомиться зональной библиотеке Южного федерального университета

Автореферат разослан О$ 2012 г.

Ученый секретарь ' И.Б. Старченко

диссертационного совета

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

Во время геомагнитных возмущений (ГВ) возрастает вероятность возникновения проблем с сердечно-сосудистой системой, изменяется чувствительность к медикаментозным воздействиям. В дни усиления солнечной активности возрастает как нижнее, так и верхнее артериальное давление, увеличивается свертываемость крови. Число нарушений сердечного ритма во время магнитной бури (МБ) увеличивается в 2,5 раза в состоянии покоя и в 1,5 раза при физической нагрузке. Число инфарктов миокарда возрастает на (13-14)%, а инсультов головного мозга на 7%, наблюдается уменьшение вариабельности как у больных, так и у клинически здоровых людей до (22-26) % по сравнению со спокойными в геомагнитном отношении днями. Количество несчастных случаев и травматизма во время МБ увеличивается в 1,5 раза. Также установлена сильная корреляция геомагнитной активности с числом ошибок в работе операторов диспетчерских служб, водителей автотранспорта, машинистов железных дорог. Достоверно выявлено влияние ГВ на выделения таких гормонов, как кортизол, дофамин, норадреналин, мелатонин.

Для прогнозирования МБ используются результаты измерения параметров ГВ наземными службами наблюдения за Солнцем и космическими аппаратами работающими в земной магнитосфере. Достоверность известных способов прогнозирования МБ составляет около 20 %. Для составления прогноза анализируется выброс коронарного вещества, ритмы вариации солнечной активности, движение солнечного ветра, направления межпланетного магнитного поля. Но породят ли перечисленные факторы масштабные вариации магнитосферы Земли остается неизвестным до наступления возмущений. В этой связи создание портативных технических средств оценки и прогнозирования ГМП в бытовых условиях, позволяющих людям, чувствительным к геомагнитным воздействиям, принимать необходимые меры является весьма актуальным.

Работа выполнена в соответствии с приоритетным направлением развития науки, технологий и техники РФ «Информационно-телекоммуникационные системы» (утверждено указом Президента РФ от 30.03.02 г. и переутверждено 21.05.06 г.); с научным направлением Южно-Российского государственного технического университета (НПИ) «Теория и принципы построения информационно-измерительных систем и систем управления» (утверждено решением ученого совета университета от 25.01.03 г. и переутверждено 01.03.06 г.).

Цель работы. Разработка метода и устройства оценки и прогнозирования влияющих на организм человека геомагнитных возмущений. Для достижения поставленной цели в рамках диссертационной работы решались следующие основные задачи по разработке:

- математической модели регулярных вариаций ГМП;

- метода выделения влияющих на организм человека возмущений ГМП на

основе разработанной математической модели регулярных вариаций ГМП;

- структур, алгоритмов функционирования устройства оценки и прогнозирования состояния ГМП;

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

Методы исследований: методы теории измерений, статистического анализа, теории вероятности, математического моделирования с использованием лицензированных пакетов прикладных программ STAT1STICA, MATLAB, MathCAD, Lab VIEW, Micro-Cap.

Научная новизна работы, состоит в следующем:

1 Предложен новый метод оценки состояния ГМП, который основан на выделении возмущений ГМП с использованием математической модели солнечно-суточных вариаций ГМП.

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

3 Разработана математическая модель солнечно-суточных вариаций ГМП с использованием ортогональных полиномов Чебышева.

4 Предложен и разработан метод прогнозирования состояния ГМП, заключающийся в определении функции плотности вероятности нормального распределения и расчета вероятности отнесения прогнозного значения к тому или иному уровню возмущённое™ по формуле Байеса.

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

Практическая значимость работы

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

Кроме того, разработан ряд программных продуктов, обеспечивающих: выделение возмущений ГМП, построение математической модели на основе ортогональных полиномов Чебышева, определение количества значимых

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

Внедрение результатов работы

Результаты работы используются (получен ряд актов внедрений) в научно-исследовательской деятельности и в учебном процессе ЮРГТУ (НПИ) (г.Новочеркасск), в СП ООО «А-мобайл» (г.Сухум), ООО «Стальмаркет» (г.Краснодар).

Основные положения, выносимые на защиту:

- метод выделения влияющих на здоровье человека возмущений ГМП;

- математические модели регулярных составляющих ГМП;

-результаты теоретических и экспериментальных исследований оценки и прогнозирования состояния ГМП;

- структуры, алгоритмы, конструктивные и схемотехнические решения информационно-измерительной системы оценки и прогнозирования влияющих на организм человека возмущений ГМП.

Апробация работы. Основные положения и научные результаты исследований докладывались на научно-технических конференциях и семинарах. В том числе: «Мехатроника. Современное состояние и тенденции развития» (г. Новочеркасск, 2009 г.); 59-я науч.-техн. конф. профессорско-преподават.состава, науч. работников, аспирантов и студентов ЮРГТУ(НПИ) «Результаты исследований 2010» (г.Новочеркасск, 2010г.); XI Междунар. науч.-практ. конф. «Моделирование.Теория, методы и средства» (г.Новочеркасск, 2011-2012гг.); 24 Междунар. науч.-практ. конф. «Математические методы в технике и технологиях» (г.Пенза, г.Киев, 2011г.); «Студенческая научная весна» (г. Новочеркасск, 2009-2011 г.г.);Всероссийская научная школа «Микроэлектронные информационно-управляющие системы и комплексы» (г.Новочеркасск, 2011г.); Микропроцессорные, аналоговые и цифровые системы: проектирование и схемотехника, теория и вопросы применения (г. Новочеркасск, 2012 г.); Методы и алгоритмы прикладной математики в технике, медицине и экономике (г.Новочеркасск, 2012г.); на научных семинарах кафедры «Информационные и измерительные системы и технологии» ЮРГТУ (НПИ).

Публикации. По теме диссертации опубликовано 24 печатные работы, в том числе 3 статьи в рецензируемых научных журналах, входящих в перечень ВАК и 4 свидетельства о регистрации программного продукта.

Структура н объем работы. Диссертация состоит из введения, четырех глав заключения, списка литературы из 128 наименований и приложений. Общий объем работы 204 страницы, включая 45 страниц приложений, 87 иллюстраций и 33 таблиц.

СПИСОК ИСПОЛЬЗУЕМЫХ СОКРАЩЕНИЙ

ГМП - геомагнитное поле; ГВ - геомагнитные возмущения; МПЗ - магнитное поле Земли; ГБ - геомагнитная буря; ПГМП - переменное геомагнитное поле; ССВ - солнечно-суточные вариации; ГГМП - горизонтальная компонента геомагнитного поля; ММ - математическая модель; МП - магнитное поле; АЦП -

аналого-цифровой преобразователь; УИ - устройство измерения; УОП -устройство оценки и прогнозирования.

СОДЕРЖАНИЕ РАБОТЫ Во введении обоснована актуальность темы, изложены цели и задачи диссертационной работы, научная новизна и практическая ценность.

В первой главе «Анализ влияния геомагнитного поля на организм человека и постановка задачи исследования» осуществлен обзор и анализ методов и устройств измерения напряженности магнитного поля Земли. Выполнен обзор литературных и библиографических источников по ключевым моментам диссертационной работы. Рассмотрены механизмы и особенности влияния ГМП на организм человека, определены основные «группы риска». Проанализированы существующие на данный момент технические средства оценки и прогнозирования ГБ.

Геомагнитное поле непрерывно регистрируют на магнитных обсерваториях постоянно действующей глобальной сети и на спутниках (рисунок 1). Получение и публикация геомагнитной информации осуществляется посредством мировых центров по сбору и обработке геомагнитной информации, которые получают данные от международных служб космической погоды и центров по сбору первичной информации. Наиболее распространенными средствами оценки и

УОН

Центры по сбору первичной информации

ЛКА

Международная служба космической погоды

УРГМП

Международная служба геомагнитных индексов

Сеть магнитных обсерваторий и спутников

(УОН) и летательные космические аппараты (ЛКА),которые осуществляют прогнозы на 3, 7, 27, 45 и 55 дней. УОН и ЛКА дают прогнозы, касающиеся лишь солнечной

активности. Но для возникновения ГБ

необходимо выполнение нескольких условий: Рисунок 1 - Глобальная сеть регистрации ГМП образования магнитного

облака или области сжатия в солнечном ветре; большая и продолжительная южная компонента межпланетарного МП; траектория движения магнитного облака должна оказаться вблизи магнитосферы Земли и вступить с ней во взаимодействие. Каждое из этих условий возникает с некоторой вероятностью, а достоверность прогноза, основанного на сопоставлении солнечных явлений магнитным бурям, составляет около 20%. Неотъемлемой частью мониторинга состояния ГМП являются устройства регистрация геомагнитного поля (УРГМП). Они осуществляют измерение напряженности МПЗ и определяют его состояние. Геомагнитные бури приводят к резкому изменению всех элементов земного магнетизма и оказывает негативное влияние на здоровье людей.

Евроазиатский центр космической погоды

—

Мировые

центры

данных

Таким образом, существует необходимость в разработке новых методов и устройств оценки и прогнозирования ГБ, способных оповещать население и специализированные организации о наступлении бури для принятия необходимых мер. Учитывая это, сформулированы задачи исследования и определены пути ее решения.

Во второй главе «Теоретический анализ составляющих ГМП и разработка метода прогнозирования его возмущений» предложен метод оценки и прогнозирования состояния ГМП, основанный на выделении возмущений ГМП посредством ММ регулярных составляющих МПЗ. Напряжённость МПЗ Н можно представить, как сумму магнитных полей, имеющих различные причины возникновения:

Я= Н0+Нт+Нг+Ни+Нц+Н1+Н»+Нс+ НР+Нт (1)

где Я0 - напряжённость МП, создаваемого однородным намагничиванием земного шара; Нт - напряжённость материкового МП, связанного с неоднородное гями строения внутренних слоев; Яа - напряжённость аномального МП, обусловленного намагниченностью верхних слоёв земной коры; На - вековые вариации МП; Ни - 11-летние вариации МП, связанные с 11-летним периодом солнечной активности; Яг - годовые вариации МП, обусловленные изменением оси вращении Земли; Нл - лунно-суточные вариации МП, связанные с положением Луны по отношению к горизонту; Яс - солнечно-суточные вариации МП, зависящие от положением горизонта к Солнцу; - рекуррентные возмущения МП, связанные с 27-дневным периодом обращения Солнца; Нт -единичные возмущения МП, связанные с кратковременными вспышками на Солнца. Из анализа выражения (1) видно, что все составляющие напряженности Я, кроме последнего, носят регулярный характер. И только последнее слагаемое Нт описывает возмущения МПЗ. Причем (Я0+Ят+Яа) образуют постоянное ГМП, неизменное во времени. Слагаемые Я0, Я,,, Я,, Нл, Яс и ЯР описывают ПГМП.

Для использования предложенного выше метода необходимо разработать ММ регулярных составляющих ГМП. На рисунке 2 приведена обобщенная схема применения ММ.

Координаты точки наблюдения

Г

Время ■

Математическая модель ПГМП

УИ

Н.,

Ф

;у

Ни

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

наблюдения. Таким образом, каждому моменту времени в модели соответствует

определенное значение

напряженности Ям для заданного места наблюдения.

Рисунок 2 - Схема определения //щ с помощью математической модели ПГМП Устройство измерения (УИ) осуществляет регистрацию напряженности ЯЮм МПЗ. Вычитающее устройство (ВУ) на основании полученных с модели Ям и с УИ Яизм выделяет единичные возмущения НЕЛ. Для реализации схемы, приведенной на рисунке 2 была разработана ММ ПГМП.

При этом рассматривался временной ряд напряженности Я, ПГМП. Согласно выражению (1) он описывается аддитивной моделью и имеет структуру:

я, +я(,+я«,+я5,

где Нщ +НС - тренд-циклическая составляющая; Я ЛЛ - сезонная составляющая; Н, - случайная составляющая.

Вековые, годовые и одиннадцатилетние вариации представлены в модели как тренд-циклические компоненты, солнечно-суточные вариации, лунно-суточные вариации и рекуррентные возмущения представлены как сезонные составляющие. Возмущения ГМП, в которых и заключены ГБ - случайная компонента. На первом этапе к исходным данных ПГМП был применен спектральный анализ. Он показал, что в ПГМП из всех перечисленных компонент наибольший вклад вносят лишь ССВ, с периодом 24 часа.

Пусть р - период последовательности Я, так что = для всякого /.

Необходимо оценить по наблюдениям напряженности ГМП Я,. Обозначим

через Н, скользящее среднее с периодом р, построенное по ряду Я,. Для

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

скользящих средних есть: Н\,#:>,...,#< .Соответственно изменим нумерацию

исходного ряда так, чтобы величине Н, соответствовал член Я,. При этом приходится отбросить целую часть от деления \р/2] первых членов ряда, для

которых значения Я, не определены. Предположим, что А: =(и +1)/), где т -положительное целое число. При этом общая длина N исходного ряда равна N = (т + 2)р.

Для каждого сезона /, 1 < / < р, рассмотрим все относящиеся к нему разности: Л(/,)- Я(г,).

л

Каждое из этих отклонений Я, от Я, можно рассматривать как результат влияния сезонных изменений. Усреднение этих разностей дает нам оценку сезонной компоненты Яи :

2 <»+1 л

Якк =-- Н "'г, для I = 1 ...р.

т +1 ы

Скорректированная оценка сезонной составляющей Я», определяется по формуле:

Я», =Я ¿Ячч , для /'= 1 ...р р

На основании исходных данных напряженности Я, ПГМП (статистика 13 лет) были выделены ССВ для каждого месяца каждого года, а затем усреднены по соответствующему месяцу. Для определения меры взаимосвязи выделенных ССВ

ГГМП между собой по формуле (2) были рассчитаны коэффициенты парной корреляции г.

(2)

гДе #„,//„,, - ¿-е значение сравниваемых временных рядов, ¿=1,2...«; Нср., ЯрГ среднее значение временных рядов.

Рассчитанные значения г лежат в пределах 0,7-0,9, что говорит о сильной взаимосвязи между месяцами всего года.

При помощи кластерного (рисунок 3) и факторного (метод главных компонент) (рисунок 4) анализов выделенные ССВ ГГМП были разделены на три обособленные группы (рисунок 3). Главные компоненты упорядочиваются в соответствии со значимостью, затем используется критерий восстановления характеристик и малозначимые компоненты отбрасываются. Использование метода главных компонент позволило сократить количество кривых (рисунок 5), с помощью которых можно описать изменения ССВ ГГМП в течение года, что обеспечивает оптимальную классификацию ССВ ГГМП.

«ж

'.'•liVT

last..-* Я.

Ё 0.2

"■ГУ) 0 / Декабрь. Январь. ' Февраль

/31о° . / \

Март. Октябрь.

Ноябрь

Апрель, Май, Июнь.

\ < . / Июль. АвгусУ, Сентябрь У

Рисунок 3 - Кластеризация ССВ ГГМП, в результате получены три обособленные группы

н,нтл-

>■' Г ) 4 j\ С 17 13 19 20 21 22 23 24

J Декабрь. Январь. Февраль

V ' Март. Октябрь. Ноябрь

Апрель. Май. Июнь. Июль Август, сентябрь

Рисунок 5 - Изменение ССВ ГГМП в течение года

Рисунок 4 - Классификация ССВ ГГМП

К первой группе относятся ССВ в декабре, январе, феврале, ко второй - в апреле, в мае, в июне, в июле, в августе, к третьей - в марте, в октябре и в ноябре. В качестве меры расстояния использовалось евклидово расстояние RR, определяемое по формуле:

где с и d - порядковый номер месяца, с, d=\..A2.

Минимальный коэффициент парной корреляции при уровне значимости а=5 % в пределах первой группы составил 0,97, второй - 0,98, а третьей - 0,92.

Проводились исследования по оценке степени взаимосвязи полученных групп фактическим значениям ГМП в магнитоспокойные дни (К-индекс не более 2) на статистике 14 лет. Коэффициент парной корреляции при уровне значимости а=5 % лежит в диапазоне 0,9^0,98.

Разработан алгоритм краткосрочного прогнозирования состояния ГМП. Он реализуется с помощью формул (3) -(6):

н„ (3)

2Ж, (4)

X /-тч/нп.-!

= а-Я'2 + (1-а)-Я;„ (5)

1)- (6)

1-1-1,-1

где а - параметр сглаживания; т - горизонт прогнозирования; п - период формирования прогноза; прогноз на т шагов; ? - время. На рисунке 6

приведено пояснение краткосрочного прогноза состояния ГМП при т=3. Было установлено, что среднеквадратическая ошибка краткосрочного прогноза ГМП

, Дляг--1 достигается при а=0,2, и

при "-5- Алгори™

| прогнозирования основан

на взвешенной оценке

"7Г................ Т предыдущего прогноза и

Для "'т- текущего фактического

I-~-1 гп ; ; /г значения. Новый прогноз

я,- я,., я,^ я,_ я,я,., я, и, ; получается в результате | I { ; корректировки

I--1

им ИЛ яи, н„ я„,яя„, им я,., я,.; ян н

-_! I_ 1 ) предыдущего прогноза с

Ллярз " я;, учетом его ошибки.

и'-.<_ я;. ; 1 Алгоритм работает

: следующим образом. : 1) Вычисляем по -1- формуле (3) Н'и- среднее " " " г значение напряженности

на интервале т-л.

Рисунок 6 — Схема алгоритма краткосрочного прогнозирования состояния ГМП

2) При т=1 прогноз осуществляется по формуле (5), где Н'_2 = Нтм. Начиная с т=2, Щ2 рассчитывается по формуле (4).

3) Затем Н'_2 принимается за фактическое значение и по формуле (5) строится предварительный прогноз#,'„,,.на основе и Н'_2.

4) Зная предыдущее прогнозное значение Я„ на т-й шаг.

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

Рисунок 7 - Схема алгоритма краткосрочного прогнозирования ГМП

Точность прогноза т) оцениваем величиной доверительного интервала 7(^1, Са) (рисунок 7) для заданной вероятности его осуществления Р(В,). Достоверность - оценка вероятности осуществления прогноза в заданном доверительном интервале, т.е Р& < Я(/т„(1) < С2). Вероятность Р(В.) определяем как частоту события А/ в п независимых опытах, где / -количество градаций возмущенности ГМП, / = 1,2,...,7. Были вычислены Р(В,), которые составили: /,(5|)=0,11, Р(В,)=0,619, Р(Вз) = 0,13, Р(В4) = 0,08, Р(В5) = 0,04, Р(В6) = 0,02, Р{В1) = 0,001

Величина доверительного интервала С,2) для рассчитанных Р(В,) определялась по формулам (7) и (8):

Р(В,У(\-Р{В,))

Р(В,)(1-Р(В,))

(7)

(8)

Для отнесения прогнозного значения к /-Й градации возмущенности ГМП используется вероятностный метод, заключающийся в том, что вероятность попадания прогнозного значения в /-ую группу , рассчитывается по формуле Байеса (9):

Р(В ) • Р. (А)

Р(В) =_ '___(9)

Л Р(В])-РВ1(А) + Р(Вг).РвХА) + ...+ Р(В1,)-Ре(АУ где Р(В1) - априорная вероятность гипотезы; Рв (А) - условная вероятность события А при выполнении гипотезыР(В:).Условная вероятность рассчитывается по формуле (10):

1 .......

Г

сЬс,

(10)

ал/2к л „

где о - среднеквадратическое отклонение, к1р1 и к,г2 - нижняя и верхняя границы /-Й градации возмущенности ГМП.

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

фазы бури, т - число, т=\...р. Это необходимо для эффективного прогноза наложенных друг на друга ГБ, следующих одна за другой.

Среднесрочный прогноз состояния ГМП основан на 27-дневной повторяемости ГБ. Рассчитаны вероятности кратной повторяемости Кп на статистике 14 лет. Полученная зависимость приведена на рисунке 8.

Проведен анализ распределения плотности Рисунок 8 - Зависимость вероятности повторения ГБ вероятности ошибки

от кратности их появления прогнозирования. Рассчитаны

параметры полученного распределения. Результаты позволяют определять ошибки прогнозирования ГБ первого и второго рода.

Разработана информационно-измерительная система оценки и прогнозирования состояния ГМП, реализующая описанный выше метод (рисунок 9).

БОП

ПО ОиП

N.

АЦП 1

МК

АЦП 2 N1,.

ПО РиО

ПК

Рисунок 9 - Структурная схема системы оценки и прогнозирования состояния ГМП

На рисунке 9 БД1, БД2 - блоки датчиков, имеющих одинаковые структуры; ПИ -плата измерения; Нх, Ну - северная и восточная составляющие вектора напряженности МПЗ; /к - ток компенсации; ДП - датчик поля, содержащий магниторезистивный мост, обмотку компенсации 1УК и катушку сброса Г-генератор импульсов сброса; УС - усилитель сброса; ДТ - датчик тока; БКП -блок компенсации поля; МК — микроконтроллер; БОП - блок оцифровки и передачи; ПК - персональный компьютер; ПО ОиП - программное обеспечение оценки и прогнозирования; ПО РиО - программное обеспечение регистрации и отображения. Система оценки и прогнозирования работает следующим образом.

При воздействии ГМП на магниторезистивный мост в ДП на его выходе появляется напряжение t/0, которое поступает на БКП. Усилитель сброса создает размагничивающие импульсы тока в катушке сброса IVqí магнитоизмерительного преобразователя. В БКП происходит преобразование напряжения в компенсирующий ток, который, проходя через ДТ поступает на компенсирующую обмотку IVK. В результате чего суммарное магнитное поле, создаваемое ГМП и током в компенсирующей обмотки, стремится к нулю, а напряжение на ДТ пропорционально измеряемому МП. Аналого-цифровой преобразователь формирует из полученного напряжения цифровой код (NK и Ny), который передается на МК. МК передает цифровой код с обоих каналов в ПК, где происходит вычисление ГГМП Н по следующему выражению:

Я = 7(ЯЛ)2+(Я,.)2.

В ПО ОиП, а также ПО РиО осуществляется оценка, прогнозирование и отображение состояния ГМП с применением предложенного метода. Вся полученная и обработанная информация накапливается в памяти ПК.

В третьем главе «Разработка математических моделей регулярных составляющих геомагнитного поля» с помощью спектрального анализа были получены математические модели изменения напряженности H¡ ГГМП в зависимости от времени суток t. Целью анализа является разложение ряда на функции синусов и косинусов различных частот, для определения тех, появление которых особенно существенно и значимо.

В общем виде нелинейное уравнение регрессии может быть представлено:

Н, = X(а,'cos(2tt•yK)-t+bl ■ sin(2rc■ yt)• /),

7=1

где üj, bj- коэффициенты регрессии модели Фурье; ук - частота циклов за единицу времени.

В результате спектрального анализа были получены ММ для выделенных в главе 2 ССВ ГГМП. Таким образом, модель изменения напряженности ССВ ГГМП //зим в декабре, январе и феврале имеет вид:

Я (0ЗИМ = -0,7 sin( 2л ■ 0,0420 + 1,79 sin( 2% ■ 0,0831)--l,42cos(2rc- 0,1250

В марте, октябре, ноябре ССВ ГГМП //осн описывается:

Н(')осч = 5>43 cos(' °>042 0 + 4>32 sin(' °>083 0

В апреле, мае, июне, июле, августе и сентябре ССВ ГГМП Н]ш описываются уравнением:

Н (Олит = 6'26 sin( 271 • °'042 0 ~ 7>32 sín(2п ' °'042 0 + + 8,65cos(2n-0,0830

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

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

Я = а0\|/0 О) + а,\|/,(*) + ... + а„\|/„ (х), где у0(х,) = 1 - полином нулевого порядка; у,(лг,) = (х + р,)■ уДх) + у- полином

полином второго порядка;

Р,

ПОЛИНОМ /7-ГО

a -I^.U) _

' Ikoof

порядка; коэффициенты

первого порядка; \|/,(л:,) = (л: + Р2)\|/|(х) + 72

У„+1(*) = (* + Р„+1)-У„« + + У^У^О) -

1*.к(*,)Г 2>,¥„-,(*>„(*,) 1к(*,)Г ' - Хк-.ООГ

связи.

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

Было установлено, что выделенные посредством разработанных ММ возмущения ГМП ЯВОзм имеют нормальный закон распределения.

Четвертая глава «Разработка и исследование информационно-измерительной системы оценки и прогнозирования влияющих на организм человека возмущений геомагнитного поля» посвящена созданию системы, реализующей обобщенную структурную схему, предложенную в главе 2, а также предложенный метод выделения влияющих на организм человека возмущений ГМП. В состав системы входят: двухкомпонентный датчик, который регистрирует северную Я, и восточную компоненты Я вектора напряженности магнитного поля Земли, со встроенными катушкой компенсации и катушкой сброса; блок компенсации поля; генератор и усилитель импульсов сброса; датчик тока; плата сбора информации Е-24; состоящая из измерительного блока с четырьмя 24-разрядными АЦП и 8-ми битным процессором А УК фирмы АТМЕЬ\ персональный компьютер.

Внешний вид разработанной системы приведен на рисунке 10. В качестве

датчиков магнитного поля

щ

Рисунок 10- ИИС оценки и прогнозирования состояния ГМП

использованы магниторезистивные преобразователи KMZ52 фирмы Phillips. Оцифровка сигнала с платы измерения осуществляется

посредством модуля Е-24 фирмы L-Card. Связь аппаратной части с компьютером осуществляется через интерфейс RS-232. На персональном компьютере в среде графического программирования Labview

реализовано программное обеспечение

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

Основные технические характеристики: количество одновременно измеряемых проекций магнитного поля: 2; диапазон измеряемого магнитного поля: ±200 мкТл; погрешность измерения: не более ±4 %; количество градаций возмущенности магнитного поля Земли: 7; потребляемая мощность платы измерения: 1,5 Вт.

Все проведенные исследования экспериментально подтвердили высокую эффективность применения предложенной модели и разработанной системы для оценки и прогнозирования состояния ГМП.

Проведены клинические испытания разработанной системы. Установлено, что для более 60 % пациентов с заболеваниями сердечно-сосудистыми наблюдается зависимость их состояния от параметров ГМП.

Разработанное устройство принято к внедрению в учебный процесс кафедры «Информационные и измерительные системы и технологии», в СП ООО «А-мобайл» (г.Сухум), ООО «Стальмаркет» (г.Краснодар).

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

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

ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ

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

2 Предложен и разработан новый метод оценки состояния ГМП, основанный на выделении возмущений ГМП посредством математической модели ССВ ГГМП. Возмущения ГМП определяются как разность между измеренным значением ГГМП и значением напряженности, полученной с помощью математической модели. Это позволяет осуществлять оценку и

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

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

4 Разработаны алгоритмы и созданы программы, позволяющие по выделенным возмущениям ГМП осуществлять краткосрочный прогноз (до 3 часов) состояния ГМП.

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

6 Разработана информационно-измерительная система медицинского назначения, предназначенная для оценки и прогнозирования состояния ГМП в бытовых условиях. В нем реализованы метод выделения возмущений ГМП посредством математической модели регулярных составляющих и алгоритмы прогнозирования активности ГМП. Устройство обеспечивает измерение напряженности ГМП с погрешностью не более ±4% и позволяет людям, чувствительным к геомагнитным возмущениям, а также медицинским и др. учреждениям получать информацию о зарождении ГБ и принимать необходимые меры.

ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ

Статьи в журналах, входящих в Перечень ВАК

1 Алексанянн Г.К. Влияние геомагнитного поля на элементы мехатронных систем / Г.К.Алексанян, М.В. Ланкин, В.И. Дубров, А.Г. Толстокожин // Известия вузов. Сев.-Кавк. регион. Техн. науки. - 2009. - Спец. выпуск: Мехатроника. Современное состояние и тенденции развития.- С. 177-179.

2 Алексанян Г.К. Разработка математической модели активности магнитного поля Земли / Г.К.Алексанян, М.В. Ланкин, В.И. Дубров, А.Г. Толстокожин // Известия вузов. Сев.-Кавк. регион. Техн. науки. - 2009. - Спец. выпуск.-С.41-46.

3 Алексанян Г.К. Разработка метода выделения влияющих на организм челоека возмущений геомагнитного поля / Г.К.Алексанян // Современные проблемы науки и образования. - 2012. -№ 3; URL: www.science-education.ru/103-6125 (дата обращения: 03.05.2012).

Статьи в других изданиях и тезисы докладов на конференциях

4 Алексанян Г.К. Разработка математической модели прогнозирования геомагнитного поля / Г.К.Алексанян, М.В. Панкин, В.И Дубров, А.Г. Толстокожин // Студенческая научная весна - 2009: Сборник научных трудов аспирантов и студентов ЮРГТУ(НПИ) / Юж.-Рос. гос. техн. ун-т (НПИ). - Новочеркасск: ЮРГТУ, 2009.-С. 165-166.

5 Алексанян Г.К. Анализ изменения суточных вариаций геомагнитного поля в течение года / Г.К.Алексанян // Результаты исследований - 2010: материалы 59-ой научно-технической конференции профессорско-преподавательского состава, научных работников, аспирантов и студентов ЮРГТУ (НПИ) // Юж.-Рос. гос. техн. ун-т (НПИ). - Новочеркасск: ЮРГТУ, 2010. - С. 231-233.

6 Алексанян Г.К. Математическая модель суточных вариаций геомагнитного поля на основе быстрого преобразования Фурье / Г.К.Алексанян // Студенческая научная весна - 2010: Сборник научных трудов аспирантов и студентов ЮРГТУ(НПИ) // Юж.-Рос. гос. техн. ун-т (НПИ). - Новочеркасск: ЮРГТУ, 2010.-С. 184-185.

7 Алексанян Г.К. Определение частотных составляющих солнечно-суточных вариаций геомагнитного поля / Г.К.Алексанян // Математические методы в технике и технологиях - ММТТ-24 : сб. тр. XXIV Междунар. науч. конф.: в Ют./ Сарат. гос. техн. ун-т - Пенза : Пенз. гос. технол. акад., 2011. - Т. 7. Секц. 11 : Информатизация и моделирование технических систем и процессов. - С. 142-143.

8 Маркарян Н.Р. Калибровка мобильного устройства прогнозирования геомагнитных возмущений / Н.Р.Маркарян, М.В.Ланкин, Г.К.Алексанян // Микроэлектронные информационно-управляющие системы и комплексы : сб. тез. и статей Всерос. науч. шк., Новочеркасск, 5-7 сент.2011г. / Юж.-Рос. гос. техн. ун-т (НПИ) - Новочеркасск : ЛИК, 2011. - С. 72-73.

9 Алексанян Г.К. Разработка метода прогнозирования геомагнитного поля / Г.К.Алексанян, М.В.Ланкин, Н.Р.Маркарян, М.В.Кондыбенко // Микроэлектронные информационно-управляющие системы и комплексы : сб. тез. и статей Всерос. науч. шк., Новочеркасск, 5-7 сент. 2011г. / Юж.-Рос. гос. техн. ун-т (НПИ) - Новочеркасск : ЛИК, 2011. -С. 29-32.

10 Алексанян Г.К. Программа построения математической модели регулярных вариаций геомагнитного поля с испльзованием ортогональных полиномов Чебышева / Г.К.Алексанян, М.В. Ланкин // Микроэлектронные информационно-управляющие системы и комплексы : сб. тез. и статей Всерос. науч. шк., Новочеркасск, 5-7 сент.2011г. / Юж.-Рос. гос. техн. ун-т (НПИ) - Новочеркасск : ЛИК, 2011. -С. 27-29.

11 Алексанян Г.К. Кластерный анализ вариаций горизонтальной составляющей геомагнитного поля/ Г.К.Алексанян // Микроэлектронные информационно-управляющие системы и комплексы : сб. тез. и статей Всерос. науч. шк., Новочеркасск, 5-7 сент.2011г. / Юж.-Рос. гос. техн. ун-т (НПИ) - Новочеркасск : ЛИК, 2011. -С. 25-26.

12 Маркарян Н.Р. Система для калибровки мобильного устройства прогнозирования геомагнитных бурь / Н.Р.Маркарян, М.В.Ланкин, Г.К.Алексанян //

Студенческая научная весна - 2011 : материалы регион, науч.-техн. конф. студент! аспирантов и молодых ученых вузов Ростовской области / Юж.-Рос. гос. техн. yi (НПИ) - Новочеркасск : ЮРГТУ, 2011. - С. 188-189.

13 Алексанян Г.К. Спектральный анализ Фурье временных рядов годов: вариаций геомагнитного поля / Г.К. Алексанян // Моделирование. Теория, методы средства : материалы XI Междунар. науч.-практ. конф., Новочеркасск, 31 марта 20 г. / Юж.-Рос. гос. техн. ун-т (НПИ) - Новочеркасск : ЮРГТУ, 2011. - С. 34-36.

14 Алексанян Г.К. Определение структуры математической моде горизонтальной составляющей геомагнитного поля / Г.К.Алексанян, М.В. Ланкин, Бертрам // Моделирование. Теория, методы и средства : материалы XI Междун; л науч.-практ. конф., Новочеркасск, 31 марта 2011 г. / Юж.-Рос. гос. техн. ун-т (НПИ) Новочеркасск : ЮРГТУ, 2011. - С. 31-36.

15 Алексанян Г.К. Исследование возмущений геомагнитного поля магнитоспокойные дни / Г.К. Алексанян // Моделирование. Теория, методы средства : материалы XI Междунар. науч.-практ. конф., г. Новочеркасск, 27 февра 2012 г. / Юж.-Рос. гос. техн. ун-т (НПИ). - Новочеркасск: ЛИК, 2012. - С. 4-6.

16 Алексанян Г.К. Разработка информационно-измерительной систел непрерывного мониторинга геомагнитного поля / Г.К. Алексанян, Ю.Р. Кревченко Моделирование. Теория, методы и средства : материалы XI Междунар. науч.-пра] конф., г.Новочеркасск, 27 февраля 2012 г. / Юж.-Рос. гос. техн. ун-т (НПИ). Новочеркасск: ЛИК, 2012. - С. 6-8.

17 Алексанян Г.К. Разработка математической модели краткосрочно прогнозирования геомагнитного поля методом экспоненциального сглаживания / Г. Алексанян, Ю.Р. Кревченко // Моделирование. Теория, методы и средства материалы XI Междунар. науч.-практ. конф., г.Новочеркасск, 27 февраля 2012 г. Юж.-Рос. гос. техн. ун-т (НПИ). - Новочеркасск: ЛИК, 2012. - С. 11-12.

18 Алексанян Г.К. Определение начальных параметров математической модели краткосрочного прогноза состояния геомагнитного поля / Г.К. Алексанян, Н.Р. Маркарян // Методы и алгоритмы прикладной математики в технике, медицине и экономике: материалы XII Междунар. науч.-практ. конф., г.Новочеркасск, 12 апр. 2012 г. / Юж.-Рос. гос. техн. ун-т (НПИ).-Новочеркасск: ЮРГТУ (НПИ), 2012. - С. 4-6.

19 Алексанян Г.К. Факторный анализ солнечно-суточных вариаций горизонтальной составляющей геомагнитного поля / Г.К. Алексанян // Методы и алгоритмы прикладной математики в технике, медицине и экономике: материалы XII Междунар. науч.-практ. конф., г. Новочеркасск, 12 апр. 2012 г. / Юж.-Рос. гос. техн. ун-т (НПИ). - Новочеркасск: ЮРГТУ(НПИ), 2012. - С. 7-9.

20 Алексанян Г.К. Информационно-измерительная система оценки и прогнозирования влияющих на организм человека возмущений геомагнитного поля / Г.К. Алексанян, Ю.Р. Кревченко // Микропроцессорные, аналоговые и цифровые системы: проектирование и схемотехника, теория и вопросы применения: материалы XII Междунар. науч.-практ. конф., г. Новочеркасск, 22 апр. 2012 г. / Юж.-Рос. гос. техн. ун-т (НПИ). - Новочеркасск: ЮРГТУ(НПИ), 2012.-С.4-7 с.

21 Свид. об. оф. per. прогр. для ЭВМ Российская Федерация Выделение регулярных вариаций геомагнитного поля (GMP 1.0) / Г.К. Алексанян, Н.И. Горбатенко [и др.]; Роспатент. - № 2012611370; заявл. 27.02.2012.

22 Свид. об. оф. per. прогр. для ЭВМ Российская Федерация. Построение математической модели регулярных вариаций геомагнитного поля (MAT_GMP 1.0) / Г.К. Алексанян, Н.И. Горбатенко [и др.]; Роспатент. -№ 2012611337;заявл. 27.02.2012.

23 Свид. об. оф. per. прогр. для ЭВМ Российская Федерация

Определение количества значимых слагаемы уравнения модели (SLAG_MOD 1.0) / Г.К. Алексанян, Н.И. Горбатенко [и др.]; Роспатент. - № 2012611384; заявл. 27.02.2012.

24 Свид. об. оф. per. прогр. для ЭВМ Российская Федерация Программа краткосрочного прогнозирования активности геомагнитного поля (FORC_GMP 1.0) / Г.К. Алексанян, Н.И. Горбатенко [и др.]; Роспатент. -№ 2012611382;заявл. 27.02.2012.

Личный вклад автора в опубликованных в соавторстве работах: [1] -

анализ влияние геомагнитного поля на технические системы; [2] - идея применения метода сезонной декомпозиции к данным геомагнитного поля, проведение экспериментальных исследований; [4] - идея построения структурной схемы; [8] - разработка структурной схемы системы; [10] -разработка структурной схемы; [9] -идея метода прогнозирования; [10] - идея использования полиномов Чебышева, разработка алгоритма; [12] - проведение исследований; [14] -идея применения кластерного анализа, проведение экспериментальных исследований; [16] - разработка структурной, функциональной схемы системы, выбор магнитоизмерительного преобразователя; [17] - разработка математической модели, проведение исследований; [18] -проведение исследований; [20] -разработка информационно-измерительной системы; [21-24] - разработка алгоритма и написание листинга программ.

Алсксанян Грайр Каренович

РАЗРАБОТКА МЕТОДА ВЫДЕЛЕНИЯ ВЛИЯЮЩИХ НА ОРГАНИЗМ ЧЕЛОВЕКА ВОЗМУЩЕНИЙ ГЕОМАГНИТНОГО ПОЛЯ ДЛЯ СОЗДАНИЯ ИНФОРМАЦИОННО-ИЗМЕРИТЕЛЬНЫХ СРЕДСТВ МЕДИЦИНСКОГО

НАЗНАЧЕНИЯ

Автореферат

Подписано в печать 08.08.2012. Формат 60x84 '/|6 Бумага офсетная. Ризография. Усл. печ. л. 1,0. Уч.-изд. л. 1,6. Тираж 110 экз. Заказ №48-5146.

Южно-Российский государственный технический университет (Новочеркасский политехнический институт)

Адрес университета: 346428, г. Новочеркасск, ул. Просвещения, 132. Отпечатано в ИД «Политехник» 346400, г. Новочеркасск, ул. Первомайская, 166. тел. 22-70-97

СПИСОК ИСПОЛЬЗУЕМЫХСОКРАЩЕНИЙ.

ВВЕДЕНИЕ.

1 АНАЛИЗ ВЛИЯНИЯ ГЕОМАГНИТНОГО ПОЛЯ НА ОРГАНИЗМ ЧЕЛОВЕКА И ПОСТАНОВКА ЗАДАЧИ ИССЛЕДОВАНИЯ.

1.1 Анализ геомагнитного поля как объекта исследования.

1.2 Влияние возмущений геомагнитного поля на организм человека.

1.3 Обзор и анализ методов и устройств измерения напряженности магнитного поля Земли.

1.4 Обзор и анализ методов и устройств прогнозирования состояния геомагнитного поля.

1.5 Выводы и постановка задачи исследования.

2 ТЕОРЕТИЧЕСКИЙ АНАЛИЗ СОСТАВЛЯЮЩИХ ГЕОМАГНИТНОГО ПОЛЯ И РАЗРАБОТКА МЕТОДА ПРОГНОЗИРОВАНИЯ ЕГО ВОЗМУЩЕННОСТИ.

2.1 Метод выделения возмущений геомагнитного поля.

2.2 Определение структуры математической модели регулярных вариаций геомагнитного поля.

2.3 Разработка алгоритмов прогнозирования геомагнитного поля.

2.4 Механизм влияния возмущений геомагнитного поля на организм человека.

2.4 Выводы по главе 2.

3 РАЗРАБОТКА МАТЕМАТЧЕСКИХ МОДЕЛЕЙ РЕГУЛЯРНЫХ СОСТАВЛЯЮЩИХ ГЕОМАГНИТНОГО ПОЛЯ.

3.1 Разработка математической модели ССВ ГГМП на основе быстрого преобразования Фурье.

3.2 Разработка математической модели ССВ ГТМП с использованием ортогональных полиномов Чебышева.

3.3 Исследование метода выделения влияющих на организм человека возмущений ГМП с применением разработанной математической мод ели.

3.4 Разработанные информационные средства для построения математических моделей и проведения исследований.

3.5 Выводы по главе 3.

4 РАЗРАБОТКА И ИССЛЕДОВАНИЕ ИНФОРМАЦИОННО

ИЗМЕРИТЕЛЬНОЙ СИСТЕМЫ ОЦЕНКИ И ПРОГНОЗИРОВАНИЯ ВЛИЯЮЩИХ НА ОРГАНИЗМ ЧЕЛОВЕКА ВОЗМУЩЕНИЙ ГЕОМАГНИТНОГО ПОЛЯ.

4.1 Информационно-измерительная система оценки и прогнозирования влияющих на организм человека возмущений геомагнитного поля.

4.2 Программное обеспечение информационно-измерительной системы.

4.3 Результаты испытаний информационно-измерительной системы оценки и прогнозирования влияющих на организм человека возмущений геомагнитного поля.

4.4 Результаты метрологических испытаний.

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

4.6 Выводы по главе 4.

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

Во время геомагнитных возмущений (ГВ) возрастает вероятность возникновения проблем с сердечно-сосудистой системой, изменяется чувствительность к медикаментозным воздействиям. В дни усиления солнечной активности возрастает как нижнее, так и верхнее артериальное давление, увеличивается свертываемость крови. Число нарушений сердечного ритма во время магнитной бури (МБ) увеличивается в 2,5 раза в состоянии покоя и в 1,5 раза при физической нагрузке. Число инфарктов миокарда возрастает на (13-14)%, а инсультов головного мозга на 7%, наблюдается уменьшение вариабельности сердечного ритма на (22-26) %. Количество несчастных случаев и травматизма во время МБ увеличивается в 1,5 раза. Также установлена сильная корреляция геомагнитной активности с числом ошибок в работе операторов диспетчерских служб, водителей автотранспорта, машинистов железных дорог. Достоверно выявлено влияние ГВ на реакции секреции и выделения таких гормонов, как кортизол, дофамин, норадреналин, мелатонин.

Для прогнозирования МБ используются результаты измерения параметров ГВ наземными службами наблюдения за Солнцем и космическими аппаратами, работающими в земной магнитосфере. Достоверность известных способов прогнозирования МБ составляет около 20 %. Для составления прогноза анализируется выброс коронарного вещества, ритмы вариации солнечной активности, движение солнечного ветра, направления межпланетного магнитного поля. Но породят ли они масштабные вариации магнитосферы Земли остается неизвестным до наступления возмущений. В этой связи особенно необходимым является создание портативных технических средств оценки и прогнозирования ГМП в бытовых условиях, позволяющих людям, чувствительным к геомагнитным воздействиям, принимать необходимые меры.

Цель работы. Разработка метода и устройства оценки и прогнозирования влияющих на организм человека геомагнитных возмущений. Для достижения поставленной цели в рамках диссертационной работы решались следующие основные задачи по разработке:

- математической модели регулярных вариаций ГМП;

- метода выделения влияющих на организм человека возмущений ГМП на основе разработанной математической модели регулярных вариаций ГМП;

- структур, алгоритмов функционирования устройства оценки и прогнозирования состояния ГМП;

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

Методы исследований: методы теории измерений, статистического анализа, теории вероятности, математического моделирования и анализа с использованием пакетов прикладных программ STATISTICАб.О, MATLAB2008, MathCAD14, LabVIEW8.5,Micro-Cap9.0.

Научная новизна работы, состоит в следующем:

1 Предложен новый метод оценки состояния ГМП, который основан на выделении возмущений ГМП с использованием математической модели солнечно-суточных вариаций ГМП.

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

3 Разработана математическая модель солнечно-суточных вариаций ГМП с использованием ортогональных полиномов Чебышева.

4 Предложен и разработан метод прогнозирования состояния ГМП, заключающийся в определении функции плотности вероятности нормального распределения и расчета вероятности отнесения прогнозного значения к тому или иному уровню возмущенности по формуле Байеса.

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

Практическая значимость работы

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

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

Внедрение результатов работы

Основные результаты работы применяются в научно-исследовательской деятельности и в учебном процессе ЮРГТУ (НПИ) (г.Новочеркасск), в СП ООО «А-мобайл» (г.Сухум), ООО «Стальмаркет» (г.Краснодар). Со всех перечисленных организаций получен ряд актов внедрений результатов работы.

Основные положения, выносимые на защиту:

- метод выделения влияющих на здоровье человека возмущений ГМП;

- математические модели регулярных составляющих ГМП;

- результаты теоретических и экспериментальных исследований оценки и прогнозирования состояния ГМП;

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

Апробация работы. Основные положения и научные результаты исследований докладывались на научно-технических конференциях и семинарах. В том числе: «Мехатроника 2009. Современное состояние и тенденции развития» (г. Новочеркасск, 2009 г.); 59-я науч.-техн. конф. профессорско-преподават.состава, науч. работников, аспирантов и студентов ЮРГТУ (НПИ) «Результаты исследований 2010» (г. Новочеркасск, 2010 г.); XI Междунар. науч.-практ. конф. «Моделирование.Теория, методы и средства» (г.Новочеркасск, 2011-2012гг.); 24 Междунар. науч.-практ. конф. «Математические методы в технике и технологиях» (г.Пенза, г.Киев, 2011 г.); «Студенческая научная весна» (г. Новочеркасск, 2009-2011 г.г.);Всероссийская научная школа «Микроэлектронные информационно-управляющие системы и комплексы» (г.Новочеркасск, 2011г.); Микропроцессорные, аналоговые и цифровые системы: проектирование и схемотехника, теория и вопросы применения (г.Новочеркасск, 2012г.); Методы и алгоритмы прикладной математики в технике, медицине и экономике (г. Новочеркасск, 2012 г.); на научных семинарах кафедры «Информационные и измерительные системы и технологии» ЮРГТУ (НПИ).

Публикации. В рамках диссертационного исследования опубликованы 24 печатные работы, в том числе 3 статьи в научных журналах, входящих в перечень ВАК и 4 свидетельства о регистрации программы для ЭВМ.

Структура и объем работы. Диссертация состоит из введения, четырех глав, заключения, списка литературы из 128 наименований и приложений. Общий объем работы 204 страниц, включая 45 страниц приложений, 87 иллюстраций и 33 таблиц.

4.6 Выводы по главе 4

1) На основании проведенных исследований разработаны алгоритмы и созданы программы, позволяющие по выделенным возмущениям ГМП осуществлять краткосрочный прогноз (до 3 часов) состояния ГМП достоверностью не менее 70%.

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

3) Разработана информационно-измерительная система, предназначенная для оценки и прогнозирования состояния ГМП в бытовых условиях. В ней реализованы метод выделения возмущений ГМП посредством математической модели регулярных составляющих и алгоритмы прогнозирования активности ГМП. Оно обеспечивает измерение напряженности ГМП с погрешностью не более 4%.

ЗАКЛЮЧЕНИЕ

В ходе выполнения диссертационной работы получены результаты, приведенные ниже:

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

2 Предложен и разработан новый метод оценки состояния ГМП. Он основан на выделении возмущений ГМП посредством математической модели ССВ ГГМП. Возмущения ГМП определяются как разность между измеренным значением ГГМП и значением напряженности, полученной с помощью математической модели. Это позволяет осуществлять оценку и последующий прогноз состояния ГМП, основываясь только на измерении напряженности МПЗ.

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

4 Разработаны алгоритмы и созданы программы, позволяющие по выделенным возмущениям ГМП осуществлять краткосрочный прогноз (до 3 часов) состояния ГМП.

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

6 Разработана информационно-измерительная система медицинского назначения, предназначенная для оценки и прогнозирования состояния ГМП в бытовых условиях. В нем реализованы метод выделения возмущений ГМП посредством математической модели регулярных составляющих и алгоритмы прогнозирования активности ГМП. Устройство обеспечивает измерение напряженности ГМП с погрешностью не более 4% и позволяет людям, чувствительным к геомагнитным возмущениям, а также медицинским и др. учреждениям получать информацию о зарождении ГБ и принимать необходимые меры.

1. Ботт, М. Внутренне строение Земли / М. Ботт. М.: Мир, 1974. - 376с.

2. Несс, Н.Ф. Геомагнитный хвост / Н.Ф. Несс // Физика атмосферы-М.: Мир, 1972.-С. 122-161.

3. Заболотная Н.А. Индексы геомагнитной активности: Справочное пособие. Изд. 2-е, перераб. / Н.А.Заболотная М.: ЛКИ, 2007. - 88 с.

4. Дубров, А.П. Геомагнитное поле и жизнь. Краткий очерк по геомагнитобиологии / А.П. Дубров. Ленинград.: Гидрометеоиздат, 1974. —175 с.

5. Яновский, Б.М. Земной магнетизм / Б.М. Яновский. Л., Изд-во Ленингр. ун-та., 1978. - 592 с.

6. Акасофу, С.И. Солнечно-земная физика 4.2. / С.И. Акасофу, С. Чепмен. М.: Мир, 1975. - 509 с.

7. Акасофу, С.И. Полярные и магнитосферные суббури / С.И. Акасофу. -М.: Мир.: 1971.-314 с.

8. Беньклва, Н.П. Шевнин А.Д. Геомагнитные поля и их вариации // Электромагнитные поля в биосфере (2 тома) / Под ред. Н.В. Красногорской. М.: Наука, 1984. С. 40-44.

9. Mayaund, P.N. Atlas of indices К. / P.N.Mayaund/ IAGA Bull-1967. -№21. -113 p.

10. Бреус, Т.К. Магнитные бури: медико-биологические и геофизические аспекты / Т.К. Бреус, С.И. Рапопорт. М.: Советский спорт, 2003.-192 с.

11. Ораевский, В.Н. Влияние геомагнитной активности на функциональное состояние организма / В.Н. Ораевский, Т.К. Бреус и др. // Биофизика. 1995. Т.43. Вып.4. - С. 819 - 826.

12. Рапопорт, С.И. Магнитные бури как стресс / С.И. Рапопорт, Т.Д. Большакова и др.. // Биофизика. 1995. -Т.43. Вып.4. - С. 632-639.

13. Комаров, Ф.И. Хронобиологические аспекты природы и характера воздействия магнитных бурь на функциональное состояние организма людей / Ф.И. Комаров, С.И. Рапопорт и др.. // Хронобиология и хрономедицина. -М.: Триада, 2000. С. 299 - 317.

14. Ашмарин, И.П. Нейрохимия в таблицах и схемах / И.П. Ашмарин, Н.Д. Ещенко и др.. М.: Экзамен, 2007. - 143 с.

15. Гурфинкель, Ю.И. Влияние геомагнитных возмущений на капиллярный кровоток больных ишемической болезнью сердца / Ю.И. Гурфинкель, В.В. Любимов и др.. // Биофизика. 1995. Т.4. Вып.4. -С. 793-800.

16. Беспятых, А.Ю. Мелатонин: теория и практика А.Ю. / Беспятых, В.Я. Бродский и др.. М.: ИД «МЕДПРАКТИКА», 2009. - С. 79-90.

17. Плеханов, Г.Ф. Дестабилизация неравновесных процессов как основа общего механизма биологического действия магнитных полей / Г.Ф. Плеханов // Реакция биологических систем на магнитное поле. -М.: Наука, 1976. С.59-81.

18. Парин, В.В. Космическая кардиология / В.В.Парин, P.M. Баевский, и др.. JL: Медицина, 1967. - 206 с.

19. Grigoriev, A.I. Advances in Space Biology and Medicine / A.I.Grigoriev, A.D. Egorov // S.I.Bonting. ed. Greenwich. London. JAI Pres. Inc. 6, 1991. V.l. P. 1-35.

20. Баевский, P.M. Прогнозирование состояний на грани нормы и патологии / P.M. Баевский. М.: Медицина, 1979. - 295 с.

21. Баевский, P.M. Математический анализ сердечного ритма в оценке особенностей адаптации организма к условиям космического полета / P.M. Баевский, Г.А. Никулина, и др.. // Вестник АМН СССР. 1984. -№4.- С. 62-69.

22. Van Ravenswaaij-Arts, С.М. Heart rate variability (review). / C.M. Van Ravenswaaij-Arts, A.A., Kollee, // Annals of Internal Medicine. 1993. - № 1184. P. 436—447.

23. Музалевская, Н.И. Характеристика возмущенного геомагнитного поля как раздражителя / Н.И. Музалевская // Проблемы космической биологии. Влияние некоторых космических и геофизических факторов на биосферу Земли. 1973. - № 4 - С. 123-142.

24. Царфис, П.Г. Биоритмология и хронотерапия / П.Г. Царфис, И.Е. Орнанский. -М.: Высшая школа, 1989 159 с.

25. Селье, Г. Очерки об адаптационном синдроме / Г. Селье.-М.: Биомедгиз, 1961. 234 с.

26. National Geophysical Data Center Solar—Terrestrial Physics Электронный ресурс. /Электрон. дан.- Боулдер: 2010-Режим ROCTyna:http://www. ngdc. noaa.gov.

27. Институт земного магнетизма, ионосферы и распространения радиоволн Электронный ресурс. / Электрон, дан. Троицк: 2009. Режим доступа: http://izmiran.ru.

28. Danmarks meteorologiske Institute Электронный ресурс. / Электрон. дан.-Копенгаген: 2009. Режим доступа: http://dmiweb.dmi.dk.

29. British geophysical survey Электронный ресурс. / Электрон, дан-Эдинбург: 2009. Режим доступа: http://www.geomag.bgs.ac.uk.

30. World Data Center for Geomagnetism Электронный ресурс. / Электрон, дан- Киото: 2009. Режим доступа: http://swdcwww.kugi.kyoto-и. ac.jp/index. html

31. National Geophysical Data Center Solar-Terrestrial Physics Электронный ресурс. / Электрон, дан- Боулдер: 2010-Режим flOCTyna:http://www.ngdc.noaa.gov

32. World Data Center for geophysics Электронный ресурс. / Электрон, дан Beijing: 2009. Режим доступа: http://gp.wdc.cn/

33. Service International of geomagnetic index Электронный ресурс. / Электрон, дан. Сэнт Мор: 2009. Режим доступа: http://isgi.cetp.ipsl.fr.

34. International Ursigram and world days service. Synoptic Codes for Solar and Geophysical Data. 1975.-P. 198.

35. Бараночников, МЛ. Микромагнитоэлектроника / M.JI. Бараночников. M.: ДКМ Пресс, 2001. - 544 с.

36. Applied Physics Systems Электронный ресурс. / Режим доступа: http:// www. appliedph ysics. com.

37. Walker Ldj Scientific: Magnetics and magnetizing systems Электронный ресурс. / Режим доступа: http:// www. walker scientific, com

38. ООО НПО ЭНТ Электронный ресурс. / Электрон, дан. Санкт Петербург: 2009. Режим доступа: http://www.nvlab.spb.ru.

39. TDK Corporation of America Электронный ресурс. / Режим доступа: http:// www, tdk-components. com.

40. Любимов, B.B. Малогабаритные, экономичные и дешевые компонентные вариометры для нужд науки и медицины / В.В. Любимов. -М.: ИЗМИР АН, 1992. 21 с.

41. Любимов, В.В. Феррозондовые диагностические магнитометры, созданные в ИЗМИР АН в период с 1989 по 1994 гг. / В.В. Любимов. М.: ИЗМИРАН, 1994.- 19 с.

42. Любимов, В.В. Диагностический магнитометр-индикатор магнитной бури / В.В. Любимов, А.А. Заруцкий. // Приборы и техника эксперимента. М.: Наука, 1996. - №2. - С. 171.

43. Любимов, В.В. Приборы «ИМПЕДАНС» для геофизики и медицины / В.В. Любимов // Приборы и техника эксперимента. М.: Наука, 2000. - №2. - С.157-158.

44. Универсальный регистратор данных. Модель IDL-07. Проспект Научно-производственной фирмы ИМПЕДАНС. М.: Наука, 1998.-21 с.

45. Магниторезистивные датчики фирмы HONEYWELL Электронный ресурс. / Режим доступа: http:// www.honeywell.com

46. Магниторезистивные датчики фирмы PHILIPS Электронный ресурс. / Режим доступа: http://www.philips.ru.

47. Ермолаев, Ю. И. О некоторых статистических взаимосвязях солнечных, межпланетных и геомагнитосферных возмущений в период 19762000 годов. / Ю.И. Ермолаев, М.Ю. Ермолаев // Космические исследования. 2003. Т41, №6. С. 115.

48. Ермолаев, Ю. И. «Буря мглою небо кроет.» / Ю.И. Ермолаев // Земля и Вселенная. 2004.- №5.-С. 13-22.

49. Веселовский, И. С. Солнечные и гелиосферные явления в октябре -ноябре 2003 г.: причины и следствия / И.С. Веселовский, М.И. Панасюк, и др.. // Космические исследования. 2004. - №5. - С. 453-508

50. Сергиенко, А. Б. Цифровая обработка сигналов / А. Б. Сергиенко. -Спб.: Питер, 2006. 751 с.

51. Павлейно, М.А. Спектральные преобразования в MatLab / М.А. Павлейно, В.М. Ромаданов СПб.: Питер, 2007. - 160 с.

52. Шахтарин, Б.И. Методы спектрального оценивания случайных процессов / Б.И. Шахтарин, В.А. Ковригин. М.: Гелиос АРВ, 2005. - 248 с.

53. Вадутов, О.С. Математические основы обработки сигналов. Практикум: учебное пособие / О.С. Вадутов Томск: Изд-во Томского политехнического университета, 2007. - 100 с.

54. Зелык, Я.И. Спектральные методы анализа вращения солнечных структур / Я.И.Зелык, H.H. Степанян // Изв. Крымской Астрофиз. Обе. -2007. -№1.~ С. 56-69.

55. Теребиж, В.Ю. Анализ временных рядов в астрофизике / В.Ю. теребиж. М.: Наука, 1992. - 392 с.

56. Электронный учебник по статистике // StatSoft. Москва,2001., WEB: http://www.statsoft.ru/home/textbook/default.htm (дата обращения: 20.03.2011).

57. Серебренников, М.Г. Выявление скрытых периодичностей / М.Г. Серебренников, A.A. Первозванский. М.: Наука, 1965. - 244 с.

58. Алексанян, Г.К. Разработка математической модели активности магнитного поля Земли / Г.К.Алексанян, М.В. Ланкин, и др. // Известия вузов. Сев.-Кавк. регион. Техн. науки. 2009. - Спец. выпуск. - С. 41-46.

59. Дуброва, Т.А. Статистические методы прогнозирования в экономике: Учеб. Пособие, практикум, тесты, программа курса / Т.А. Дуброва, М.Ю. Архипова. М.: Моск. гос. ун-т экон., стат. и инфор., 2004. -136 с.

60. Тюрин, Ю.Н. Статистический анализ данных на компьютере / Ю.Н.Тюрин, A.A. Макаров. М.: ИНФРА-М, 1998. - 528 с.

61. Боровиков, В.П. Прогнозирование в системе STATISTICA в среде Windows. Основы теории и интенсивная практика на компьютере / В.П. Боровиков, Г.И. Ивченко. М.: Финансы и статистика, 2000. - 384 с.

62. Боровиков, В.П. STATISTICA. Искусство анализа данных на компьютере: Для профессионалов / В.П. Боровиков. СПб.: Питер, 2003. -688 с.

63. Сошникова, Л.А. Многомерный статистический анализ в экономике: Учеб. пособие для вузов / Л.А. Сошникова, В.Н. Тамашевич, и др.. М.: ЮНИТИ-ДАНА, 1999. - 598 с.

64. Ким, Дж.-О. Факторный анализ: статистические методы и практические вопросы / Дж.-О Ким, Ч.У. Мьюллер. М.: Финансы и статистика, 1989. -215 с.

65. Шуметов, В.Г. Факторный анализ: подход с применением ЭВМ / В.Г. Шуметов, Л.В. Шуметова. Орел.: ОрелГТУ, 1999. - 88 с.

66. Енюкова, И.С. Факторный, дискриминантный и кластерный анализ / И.С. Енюкова. М.: Финансы и статистика, 1989. - 215 с.

67. Лоули Д. Факторный анализ как статистический метод / Д. Лоули, А. Максвелл. М.: Мир, 1964. - 224 с.

68. Львовский, E.H. Статистические методы построения эмпирических формул: Учеб. пособие для втузов / E.H. Львовский. М.: Высш. шк., 1988. -239 с.

69. Нуссбаумер, Г. Быстрое преобразование Фурье и алгоритмы вычисления сверток / Г. Нуссбаумер. М.: Радио и связь, 1985. - 248 с.

70. Рабинер, Л. Теория и применение цифровой обработки сигналов / Л. Рабинер, Б. Гоулд. М.: Мир, 1978. - 848 с.

71. Хуанга, Т.С. Быстрые алгоритмы в цифровой обработке изображений / Т.С. Хуанга. М.: Радио и связь, 1984. - 221 с.

72. Снеддон, И. Преобразование Фурье / И. Снедцон. М.: 1955. - 668с.

73. Марпл-мл., С.Л. Цифровой спектральный анализ и его приложения / С.Л. Марпл-мл. М.: Мир, 1990. - 584 с.

74. Баскаков, С.И. Радиотехнические цепи и сигналы: Учебник для вузов / С.И. Баскаков. М.: Высшая школа, 2000. - 462 с.

75. Хемминг, Р.В. Цифровые фильтры / Р.В.Хемминг. М.: Недра, 1987.-221 с.

76. Сергиенко А.Б. Цифровая обработка сигналов / А.Б. Сергиенко -СПб.: Питер, 2003. 608 с.

77. Васильев, К.К. Теория электрической связи: учебное пособие / К.К. Васильев, В.А. Глушков, и др.. Ульяновск: УлГТУ, 2008. - 452 с.

78. Рангайян, P.M. Анализ биомедицинских сигналов. Практический подход / P.M. Рангайян. М.: ФИЗМАТЛИТ, 2007 - 440 с.

79. Лайонс, Р. Цифровая обработка сигналов / Р. Лайонс.- М.: ООО «Бином-Пресс», 2006. 656 с.

80. Cooley, J.W. An Algorithm for the Machine Computation of Complex Fourier Series / J.W. Cooley, J.W. Tukey //Mathematics Computation, 1965. -Vol.\9.-P. 297-301.

81. Fredrick, J. On the Use of Windows for Harmonic Analysis with the Discrete Fourier Transform / J. Fredrick // Proc. IEEE, 1978 Vol. 66, №1. - P. 51-83.

82. Кривицкий Б.Х. Справочник по теоретическим основам радиоэлектроники / Б.Х. Кривицкий. М.: Энергия, 1977. - 504 с.

83. Центр прогнозов космической погоды ИЗМИР АН Электронный ресурс. / Электрон, дан. Троицк: 2010. Режим доступа: http://forecast.izmiran.ru/

84. Елисеева, И.И. Общая теория статистики / И.И. Елисеева, М.М. Юзбашев. М.: Финансы и Статистика, 2002. - 480 с.

85. Стюард, И.Г. Введение в фурье-оптику / И.Г. Стюард М.: Мир, 1985.- 184 с.

86. Черногоров, Л.Ф. Современные методы спектрального анализа квазипериодических и волновых процессов в ионосфере: особенности и результаты экспериментов / Л.Ф. Черногоров // Геомагнетизм и астрономия. -2008. -Т.48. -№5. С. 681-702.

87. Zalyubovsky I. The Earth Atmosphere Geospace System: Main Properties, Processes and Phenomena /I. Zalyubovsky, L. Chernogor // Space Research in Ukraine. The Report Prepared by the Space Research Institute of NASU-NSAU. Kyiv. 2008. - P. 19-29.

88. Зелык, Я.И. Спектральные методы анализа вращения солнечных структур / Я.И. Зелык, Н.Н. Степанян, и др.. // Изв. Крымской Астрофиз. Обе. 103. 2007. -№1. - С. 56-69.

89. Jones Т.Е. The detection of atmospheric waves produced by the total solar eclipse 11 August 1999 / T.B. Jones, D.M. Wright et al. II J. Atmos. Sol.-Terr. Phys.-2004. Vol. 66, -№5. -P. 363-374.

90. Охоцимский, Д.Е. Основы механики космического полета: Учеб. пособие / Д.Е. Охоцимский, Ю.Г. Сихарулидзе. М.: Наука, 1990. - 448 с.

91. Госсард, Э.Э. Волны в атмосфере / Э.Э. Госсард, У.Х. Хук. М.: Мир, 178,-532 с.

92. Афонский, А.А. Цифровые анализаторы спектра, сигналов и логики / А.А. Афонский, В.П. Дьяконов. М.: СОЛОН-Пресс, 2009. - 248 с.

93. Дьяконов, В.П. Современная осциллография и осциллографы / В.П. Дьяконов. М.: СОЛОН-Пресс, 2004. - 320 с.

94. Глинченко, А.С. Цифровая обработка сигналов / А.С. Глинченко. Красноярск: Изд-во КГТУ, 2001.-383 с.

95. Рабинер, Л. Теория и применение цифровой обработки сигналов / Л. Рабинер, Б. Гоулд. М.: Мир, 1978. - 848 с.

96. Гольденберг, Л.М. Цифровая обработка сигналов. Учебное пособие для вузов / Л.М. Гольденберг. М.: Радио и связь, 1990. - 256 с.

97. Блейхут, Р. Быстрые алгоритмы цифровой обработки сигналов / Р. Блейхут. М.: Мир, 1989. - 448 с.

98. Злобин, С.JI. Двумерное быстрое преобразование Хартли в цифровой обработке изображений / С.Л. Злобин, А.Я. Стальной // Доклады 6-ой Международной конференции «Цифровая обработка сигналов и её применение». 2004. - Том 2. - С. 114-116.

99. Брейсуэлл, Р. Преобразование Хартли. Теория и приложения / Р. Брейсуэлл. М.: Мир, 1990. - 175 с.

100. Магнус, Я.Р. Эконометрика. Начальный курс / Я.Р. Магнус, П.К. Катышев и др.. М.: Дело, 2004. - 576 с.

101. Лютар, 3. Статистическое оценивание / 3. Лютар. М.: Статистика, 1976. - 598 с.

102. Горбатенко, Н.И. Планирование эксперимент: учеб.пособие / Н.И. Горбатенко, М.В. Ланкин и др. / Юж.-Рос. гос. техн. ун-т (НПИ). -Новочеркасск: Оникс+, 2007. 120 с.

103. Львовский, E.H. Статистические методы построения эмпирических формул: Учеб. пособие для втузов / E.H. Львовский. М.: Высш. шк., 1988.-239 с.

104. Львовский, E.H. Пассивный и активный эксперимент при исследовании механических характеристик бетона / E.H. Львовский. Кишинев: Картя Молдовеняскэ, 1988 - 176 с.

105. Дьяконов, В.П. MATLAB. Обработка сигналов и изображений. Специальный справочник / В.П. Дьяконов. СПб.: Питер, 2002. - 608 с.

106. Дьяконов, В.П. MATLAB 6: учебный курс / В.П. Дьяконов СПб.: Питер, 2001.-592 с.

107. Курбатова, Е.А. MATLAB 1. Самоучитель / Е.А. Курбатова. М.: Вильяме, 2006. - 256 с.

108. Макаров, Е.Г. Инженерные расчеты в Mathcad. Учебный курс / Е.Г. Макаров. СПб.: Питер-Юг, 2005. - 448 с.

109. Электронный учебник по статистике/VStatSoft. Москва,2001., WEB: http://www.statsoft.ru/home/textbook/default.htm (дата обращения: 20.03.2011).

110. Немчинов, B.C. Полиномы Чебышева и математическая статистика / B.C. Немчинов. М.: С.-х. акад. им. К. А. Тимирязева, 1946. -139 с.

111. ГОСТ 8.207-76. Прямые измерения с многократными наблюдениями. Методы обработки результатов наблюдений. М.: Изд-во стандартов, 1977. - 10 с.

112. Ивановский, Р.И. Теория вероятностей и математическая статистика. Основы, прикладные аспекты, с примерами и задачами в среде Mathcad! Р.И. Ивановский. СПб.: БХВ-Перербург, 2008. - 528 с.

113. Болышев, JI.H. Таблицы математической статистики / JI.H. Болышев, Н.В. Смирнов. М.: Наука, 1983. - 416 с.

114. Мартынов, Г.В. Критерии омега-квадрат / Г.В. Мартынов. М.: Наука, 1978. - 80 с.

115. Суранов, А.Я. LabVIEW 8.20: Справочник по функциям / А .Я. Суранов. М.: ДМК Пресс, 2007. - 536 с.

116. Рабинович, С.Г. Погрешности измерений / С.Г. Рабинович. JL: Энергия, 1978.-262 с.

117. Сергеев, А.Г. Метрология: Учеб.пособие для вузов / А.Г. Сергеев, В.В. Крохин М.: Логос, 2001. - 408 с.