Разработка метода выделения влияющих на организм человека возмущений геомагнитного поля для создания информационно-измерительных средств медицинского назначения
На правах рукописи
Алексанян Грайр Каренович РАЗРАБОТКА МЕТОДА ВЫДЕЛЕНИЯ ВЛИЯЮЩИХ НА ОРГАНИЗМ ЧЕЛОВЕКА ВОЗМУЩЕНИЙ ГЕОМАГНИТНОГО ПОЛЯ ДЛЯ СОЗДАНИЯ ИНФОРМАЦИОННО-ИЗМЕРИТЕЛЬНЫХ СРЕДСТВ МЕДИЦИНСКОГО НАЗНАЧЕНИЯ Специальность 05.11.17 – Приборы, системы и изделия медицинского назначения
АВТОРЕФЕРАТ
диссертации на соискание ученой степени кандидата технических наук
Новочеркасск 2012 2
Работа выполнена на кафедре «Информационные и измерительные системы и технологии» ФБГОУ ВПО «Южно-Российский государственный технический университет (Новочеркасский политехнический институт)»
Научный консультант: - доктор технических наук, профессор Н.И. Горбатенко (ЮРГТУ (НПИ), г.Новочеркасск)
Официальные оппоненты: - доктор технических наук, профессор Н.Н. Чернов (ТТИ ЮФУ, г.Таганрог) - кандидат технических наук, доцент А.В. Литвин (ДГТУ, г.Ростов-на-Дону ) Ведущее предприятие: Федеральное государственное унитарное предприятие «Научно-исследовательский институт специальных информационно-измерительных систем», г. Ростов-на-Дону
Защита состоится « 27 » сентября 2012 г. в 1420 ч. на заседании диссертационного Совета Д212.208.23 при Южном федеральном университете по адресу: 347928, г. Таганрог, ул.Шевченко, 2., ауд.Е-306.
С диссертацией можно ознакомиться зональной библиотеке Южного федерального университета
Автореферат разослан «_» 2012 г.
Ученый секретарь И.Б. Старченко диссертационного совета Актуальность проблемы. Человек в процессе своей жизни подвергается влиянию гелио-метеорологических и геофизических факторов. Одним из таких факторов является геомагнитное поле (ГМП), возмущения которого оказывают большое влияние как на больных, так и на практически здоровых людей.
Во время геомагнитных возмущений (ГВ) возрастает вероятность возникновения проблем с сердечно-сосудистой системой, изменяется чувствительность к медикаментозным воздействиям. В дни усиления солнечной активности возрастает как нижнее, так и верхнее артериальное давление, увеличивается свертываемость крови. Число нарушений сердечного ритма во время магнитной бури (МБ) увеличивается в 2,5 раза в состоянии покоя и в 1, раза при физической нагрузке. Число инфарктов миокарда возрастает на (13 14) %, а инсультов головного мозга на 7 %, наблюдается уменьшение вариабельности как у больных, так и у клинически здоровых людей до (22-26) % по сравнению со спокойными в геомагнитном отношении днями. Количество несчастных случаев и травматизма во время МБ увеличивается в 1,5 раза. Также установлена сильная корреляция геомагнитной активности с числом ошибок в работе операторов диспетчерских служб, водителей автотранспорта, машинистов железных дорог. Достоверно выявлено влияние ГВ на выделения таких гормонов, как кортизол, дофамин, норадреналин, мелатонин.
Для прогнозирования МБ используются результаты измерения параметров ГВ наземными службами наблюдения за Солнцем и космическими аппаратами, работающими в земной магнитосфере. Достоверность известных способов прогнозирования МБ составляет около 20 %. Для составления прогноза анализируется выброс коронарного вещества, ритмы вариации солнечной активности, движение солнечного ветра, направления межпланетного магнитного поля. Но породят ли перечисленные факторы масштабные вариации магнитосферы Земли остается неизвестным до наступления возмущений. В этой связи создание портативных технических средств оценки и прогнозирования ГМП в бытовых условиях, позволяющих людям, чувствительным к геомагнитным воздействиям, принимать необходимые меры является весьма актуальным.
Работа выполнена в соответствии с приоритетным направлением развития науки, технологий и техники РФ «Информационно-телекоммуникационные системы» (утверждено указом Президента РФ от 30.03.02 г. и переутверждено 21.05.06 г.);
с научным направлением Южно-Российского государственного технического университета (НПИ) «Теория и принципы построения информационно-измерительных систем и систем управления» (утверждено решением ученого совета университета от 25.01.03 г. и переутверждено 01.03.06 г.).
Цель работы. Разработка метода и устройства оценки и прогнозирования влияющих на организм человека геомагнитных возмущений. Для достижения поставленной цели в рамках диссертационной работы решались следующие основ ные задачи по разработке:
– математической модели регулярных вариаций ГМП;
– метода выделения влияющих на организм человека возмущений ГМП на основе разработанной математической модели регулярных вариаций ГМП;
– структур, алгоритмов функционирования устройства оценки и прогнозирования состояния ГМП;
– системы оценки и прогнозирования состояния ГМП, реализующей предложенный метод и позволяющей получать информацию о возмущениях ГМП, оказывающих негативное влияние на здоровье метеочувствительных людей.
Методы исследований: методы теории измерений, статистического анализа, теории вероятности, математического моделирования с использованием лицензированных пакетов прикладных программ STATISTICA, MATLAB, MathCAD, LabVIEW, Micro-Cap.
Научная новизна работы, состоит в следующем:
1 Предложен новый метод оценки состояния ГМП, который основан на выделении возмущений ГМП с использованием математической модели солнечно суточных вариаций ГМП.
2 Разработана математическая модель солнечно-суточных вариаций ГМП на основе быстрого преобразования Фурье, описывающая изменение напряженности ГМП в течение суток.
3 Разработана математическая модель солнечно-суточных вариаций ГМП с использованием ортогональных полиномов Чебышева.
4 Предложен и разработан метод прогнозирования состояния ГМП, заключающийся в определении функции плотности вероятности нормального распределения и расчета вероятности отнесения прогнозного значения к тому или иному уровню возмущенности по формуле Байеса.
Обоснованность и достоверность научных положений и выводов подтверждаются корректностью допущений, принимаемых при математическом моделировании и при разработке устройств;
согласованием теоретических положений с результатами экспериментальных исследований;
непротиворечивостью предложенных моделей, критическим обсуждением основных результатов работы с ведущими специалистами в области моделирования и магнитоизмерительной и медицинской техники.
Практическая значимость работы Разработана информационно-измерительная система оценки и прогнозирования состояния ГМП, реализующая предложенный выше метод.
Система регистрирует северную и восточную составляющие вектора напряженности магнитного поля Земли, рассчитывает горизонтальную составляющую и на основе заданного алгоритма осуществляет прогнозирование состояния ГМП. Данная система позволяет медицинским организациям, службам ЧС, а также непосредственно населению получать ценную информацию о возмущенности ГМП с целью уменьшения количества несчастных случаев, минимизации негативных последствий и планирования работ.
Кроме того, разработан ряд программных продуктов, обеспечивающих:
выделение возмущений ГМП, построение математической модели на основе ортогональных полиномов Чебышева, определение количества значимых слагаемых в уравнении модели регулярных вариаций ГМП, краткосрочное прогнозирования состояния ГМП.
Внедрение результатов работы Результаты работы используются (получен ряд актов внедрений) в научно исследовательской деятельности и в учебном процессе ЮРГТУ (НПИ) (г.Новочеркасск), в СП ООО «А-мобайл» (г.Сухум), ООО «Стальмаркет» (г.Краснодар).
Основные положения, выносимые на защиту:
– метод выделения влияющих на здоровье человека возмущений ГМП;
– математические модели регулярных составляющих ГМП;
– результаты теоретических и экспериментальных исследований оценки и прогнозирования состояния ГМП;
– структуры, алгоритмы, конструктивные и схемотехнические решения информационно-измерительной системы оценки и прогнозирования влияющих на организм человека возмущений ГМП.
Апробация работы. Основные положения и научные результаты исследований докладывались на научно-технических конференциях и семинарах. В том числе: «Мехатроника. Современное состояние и тенденции развития» (г. Новочеркасск, 2009 г.);
59-я науч.-техн. конф. профессорско-преподават.состава, науч. работников, аспирантов и студентов ЮРГТУ(НПИ) «Результаты исследований 2010» (г. Новочеркасск, 2010 г.);
XI Междунар. науч.-практ. конф.
«Моделирование.Теория, методы и средства» (г.Новочеркасск, 2011-2012гг.);
24 Междунар. науч.-практ. конф. «Математические методы в технике и технологиях» (г.Пенза, г.Киев, 2011 г.);
«Студенческая научная весна» (г.
Новочеркасск, 2009-2011 г.г.);
Всероссийская научная школа «Микроэлектронные информационно-управляющие системы и комплексы» (г. Новочеркасск, 2011 г.);
Микропроцессорные, аналоговые и цифровые системы: проектирование и схемотехника, теория и вопросы применения (г. Новочеркасск, 2012 г.);
Методы и алгоритмы прикладной математики в технике, медицине и экономике (г. Новочеркасск, 2012 г.);
на научных семинарах кафедры «Информационные и измерительные системы и технологии» ЮРГТУ (НПИ).
Публикации. По теме диссертации опубликовано 24 печатные работы, в том числе 3 статьи в рецензируемых научных журналах, входящих в перечень ВАК и свидетельства о регистрации программного продукта.
Структура и объем работы. Диссертация состоит из введения, четырех глав, заключения, списка литературы из 128 наименований и приложений. Общий объем работы 204 страницы, включая 45 страниц приложений, 87 иллюстраций и таблиц.
СПИСОК ИСПОЛЬЗУЕМЫХ СОКРАЩЕНИЙ ГМП – геомагнитное поле;
ГВ – геомагнитные возмущения;
МПЗ – магнитное поле Земли;
ГБ – геомагнитная буря;
ПГМП – переменное геомагнитное поле;
ССВ – солнечно-суточные вариации;
ГГМП – горизонтальная компонента геомагнитного поля;
ММ – математическая модель;
МП – магнитное поле;
АЦП – аналого-цифровой преобразователь;
УИ – устройство измерения;
УОП – устройство оценки и прогнозирования.
СОДЕРЖАНИЕ РАБОТЫ
Во введении обоснована актуальность темы, изложены цели и задачи диссертационной работы, научная новизна и практическая ценность.
В первой главе «Анализ влияния геомагнитного поля на организм человека и постановка задачи исследования» осуществлен обзор и анализ методов и устройств измерения напряженности магнитного поля Земли.
Выполнен обзор литературных и библиографических источников по ключевым моментам диссертационной работы. Рассмотрены механизмы и особенности влияния ГМП на организм человека, определены основные «группы риска».
Проанализированы существующие на данный момент технические средства оценки и прогнозирования ГБ.
Геомагнитное поле непрерывно регистрируют на магнитных обсерваториях постоянно действующей глобальной сети и на спутниках (рисунок 1). Получение и публикация геомагнитной информации осуществляется посредством мировых центров по сбору и обработке геомагнитной информации, которые получают данные от международных служб космической погоды и центров по сбору первичной информации. Наиболее распространенными средствами оценки и прогнозирования состояния ГМП являются устройства оптического наблюдения (УОН) и летательные Центры по сбору первичной космические аппараты информации (ЛКА),которые УОН осуществляют прогнозы Международная на 3, 7, 27, 45 и 55 дней.
служба космической погоды Мировые УОН и ЛКА дают ЛКА центры данных прогнозы, касающиеся Международная лишь солнечной служба геомагнитных УРГМП активности. Но для индексов возникновения ГБ Сеть магнитных необходимо выполнение Евроазиатский центр обсерваторий и космической погоды нескольких условий:
спутников образования магнитного Рисунок 1 – Глобальная сеть регистрации ГМП облака или области сжатия в солнечном ветре;
большая и продолжительная южная компонента межпланетарного МП;
траектория движения магнитного облака должна оказаться вблизи магнитосферы Земли и вступить с ней во взаимодействие. Каждое из этих условий возникает с некоторой вероятностью, а достоверность прогноза, основанного на сопоставлении солнечных явлений магнитным бурям, составляет около 20 %. Неотъемлемой частью мониторинга состояния ГМП являются устройства регистрация геомагнитного поля (УРГМП).
Они осуществляют измерение напряженности МПЗ и определяют его состояние.
Геомагнитные бури приводят к резкому изменению всех элементов земного магнетизма и оказывает негативное влияние на здоровье людей.
Таким образом, существует необходимость в разработке новых методов и устройств оценки и прогнозирования ГБ, способных оповещать население и специализированные организации о наступлении бури для принятия необходимых мер. Учитывая это, сформулированы задачи исследования и определены пути ее решения.
Во второй главе «Теоретический анализ составляющих ГМП и разработка метода прогнозирования его возмущений» предложен метод оценки и прогнозирования состояния ГМП, основанный на выделении возмущений ГМП посредством ММ регулярных составляющих МПЗ.
Напряжённость МПЗ Н можно представить, как сумму магнитных полей, имеющих различные причины возникновения:
Н = Н0+Нm+На+НВ+Н11+НГ+НЛ+НС+ НР+НЕД (1) где Н0 – напряжённость МП, создаваемого однородным намагничиванием земного шара;
Нm – напряжённость материкового МП, связанного с неоднородностями строения внутренних слоёв;
На – напряжённость аномального МП, обусловленного намагниченностью верхних слоёв земной коры;
НВ – вековые вариации МП;
Н11 – 11-летние вариации МП, связанные с 11-летним периодом солнечной активности;
НГ – годовые вариации МП, обусловленные изменением оси вращении Земли;
НЛ – лунно-суточные вариации МП, связанные с положением Луны по отношению к горизонту;
НС – солнечно-суточные вариации МП, зависящие от положением горизонта к Солнцу;
НР – рекуррентные возмущения МП, связанные с 27-дневным периодом обращения Солнца;
НЕД – единичные возмущения МП, связанные с кратковременными вспышками на Солнца. Из анализа выражения (1) видно, что все составляющие напряженности Н, кроме последнего, носят регулярный характер. И только последнее слагаемое НЕД описывает возмущения МПЗ. Причем (Н0+Нm+На) образуют постоянное ГМП, неизменное во времени. Слагаемые НВ, Н11, НГ, НЛ, НС и НР описывают ПГМП.
Для использования предложенного выше метода необходимо разработать ММ регулярных составляющих ГМП. На рисунке 2 приведена обобщенная схема применения ММ.
Входными параметрами Координаты точки наблюдения модели являются время и координаты точки ВУ Математическая НМ НЕД наблюдения. Таким образом, Время модель ПГМП каждому моменту времени в модели соответствует Низм УИ определенное значение напряженности НМ для Рисунок 2 – Схема определения НЕД с помощью заданного места наблюдения.
математической модели ПГМП Устройство измерения (УИ) осуществляет регистрацию напряженности НИЗМ МПЗ. Вычитающее устройство (ВУ) на основании полученных с модели НМ и с УИ НИЗМ выделяет единичные возмущения НЕД. Для реализации схемы, приведенной на рисунке 2 была разработана ММ ПГМП.
При этом рассматривался временной ряд напряженности Н t ПГМП.
Согласно выражению (1) он описывается аддитивной моделью и имеет структуру:
Н t Н TR Н C Н SS Н, t t t t где Н TR Н C – тренд-циклическая составляющая;
Н SS – сезонная составляющая;
t t t Н – случайная составляющая.
t Вековые, годовые и одиннадцатилетние вариации представлены в модели как тренд-циклические компоненты, солнечно-суточные вариации, лунно-суточные вариации и рекуррентные возмущения представлены как сезонные составляющие.
Возмущения ГМП, в которых и заключены ГБ – случайная компонента. На первом этапе к исходным данных ПГМП был применен спектральный анализ. Он показал, что в ПГМП из всех перечисленных компонент наибольший вклад вносят лишь ССВ, с периодом 24 часа.
Пусть p – период последовательности Н SS, так что Н SS Н SS для всякого t. t p t t Необходимо оценить Н SS по наблюдениям напряженности ГМП Нt. Обозначим t через Н t скользящее среднее с периодом p, построенное по ряду Н t. Для упрощения обозначений начнем нумерацию величин Н t с единицы, так что ряд из скользящих средних есть: Н 1, Н 2,..., Н k. Соответственно изменим нумерацию исходного ряда так, чтобы величине Н t соответствовал член Н t. При этом приходится отбросить целую часть от деления [p/2] первых членов ряда, для которых значения Н t не определены. Предположим, что k m 1p, где m – положительное целое число. При этом общая длина N исходного ряда равна N m 2p.
Для каждого сезона i, 1 i p, рассмотрим все относящиеся к нему разности: H (ti ) Н (ti ), H (ti p ) Н (ti p ),,…, H (timp ) Н (timp ),.
Каждое из этих отклонений H i от H i можно рассматривать как результат влияния сезонных изменений. Усреднение этих разностей дает нам оценку сезонной компоненты Н SS : t 1 m1 H il p Н il p, для i = 1…p.
Н SS m 1 l i * Скорректированная оценка сезонной составляющей Н SS определяется по i формуле:
* 1p Н SS Н SS Н SS, для i = 1…p i p i 1 i i На основании исходных данных напряженности Н t ПГМП (статистика лет) были выделены ССВ для каждого месяца каждого года, а затем усреднены по соответствующему месяцу. Для определения меры взаимосвязи выделенных ССВ ГГМП между собой по формуле (2) были рассчитаны коэффициенты парной корреляции r.
Н Н срс Н id Н срс n ic r i, (2) Н Н срс Н Н срd n n 2 ic id i i 1 i где Н ic, Н id i – i-е значение сравниваемых временных рядов, i=1,2…n;
Н срс, Н срd – среднее значение временных рядов.
Рассчитанные значения r лежат в пределах 0,7–0,9, что говорит о сильной взаимосвязи между месяцами всего года.
При помощи кластерного (рисунок 3) и факторного (метод главных компонент) (рисунок 4) анализов выделенные ССВ ГГМП были разделены на три обособленные группы (рисунок 3). Главные компоненты упорядочиваются в соответствии со значимостью, затем используется критерий восстановления характеристик и малозначимые компоненты отбрасываются. Использование метода главных компонент позволило сократить количество кривых (рисунок 5), с помощью которых можно описать изменения ССВ ГГМП в течение года, что обеспечивает оптимальную классификацию ССВ ГГМП.
Факторные нагрузки, Фактор 1 и Фактор RR Метод клавных компонент 1, 12 1 0, 0, Декабрь, Январь, Февраль 0, Фактор 3 0, Март, Октябрь, Ноябрь 0, Апрель, Май, Июнь, Июль, Август, Сентябрь -0, -0, -0, -1,1 -1 -0,9 -0,8 -0,7 -0,6 -0,5 -0, Фактор Рисунок 3 – Кластеризация ССВ ГГМП, в Рисунок 4 – Классификация ССВ ГГМП результате получены три обособленные группы К первой группе относятся Н,нТл ССВ в декабре, январе, феврале, ко второй – в апреле, в мае, в июне, в июле, в августе, к третьей – в марте, 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 в октябре и в ноябре. В качестве t, часы - меры расстояния использовалось Декабрь, Январь, Февраль - евклидово расстояние RR, -15 Март, Октябрь, Ноябрь определяемое по формуле:
- RR ( Н ic Н id ) 2, - Апрель, Май, Июнь, Июль, Август, сентябрь i - где c и d – порядковый номер Рисунок 5 – Изменение ССВ ГГМП месяца, c, d =1…12.
в течение года Минимальный коэффициент парной корреляции при уровне значимости =5 % в пределах первой группы составил 0,97, второй – 0,98, а третьей – 0,92.
Проводились исследования по оценке степени взаимосвязи полученных групп фактическим значениям ГМП в магнитоспокойные дни (К-индекс не более 2) на статистике 14 лет. Коэффициент парной корреляции при уровне значимости =5 % лежит в диапазоне 0,90,98.
Разработан алгоритм краткосрочного прогнозирования состояния ГМП. Он реализуется с помощью формул (3) –(6):
1 Hi, Н (3) n i n, Н 2 H i, (4) i ( n1), Н ПР Н 2 (1 ) Н 1, (5) H Н ПРi Н ПР ( 1), (6) ПРi i 1, где – параметр сглаживания;
– горизонт прогнозирования;
n – период формирования прогноза;
Н ПР – прогноз на шагов;
t – время. На рисунке приведено пояснение краткосрочного прогноза состояния ГМП при =3. Было установлено, что среднеквадратическая ошибка краткосрочного прогноза ГМП достигается при =0,2, и при n 5. Алгоритм прогнозирования основан на взвешенной оценке предыдущего прогноза и текущего фактического значения. Новый прогноз получается в результате корректировки предыдущего прогноза с учетом его ошибки.
Алгоритм работает следующим образом.
1) Вычисляем по формуле (3) Н 1 – среднее значение напряженности на интервале ·n.
Рисунок 6 – Схема алгоритма краткосрочного прогнозирования состояния ГМП 2) При =1 прогноз осуществляется по формуле (5), где Н 2 = Низм. Начиная с =2, Н 2 рассчитывается по формуле (4).
3) Затем Н 2 принимается за фактическое значение и по формуле (5) строится предварительный прогноз Н ПР на основе Н 1 и Н 2.
4) Зная предыдущее прогнозное значение H ПРi 1, рассчитывается прогноз Н ПРi на -й шаг.
Точность прогноза оцениваем величиной доверительного интервала Н(t) I(1, 2) (рисунок 7) для заданной f(x) вероятности его осуществления P(Bi).
Н(tn+) Достоверность – оценка вероятности Н осуществления прогноза в заданном ПРi Н(tn) доверительном интервале, т.е P(1 Н(tn+) 2). Вероятность P(Bi) Н(t1) определяем как частоту события Аi в n независимых опытах, где i – количество градаций возмущенности ГМП, i = 1,2,…,7. Были вычислены t1 t2 t3 t4 tn tn+ t n P(Bi), которые составили: P(B1)=0,11, Рисунок 7 – Схема алгоритма краткосрочного P(B2)=0,619, P(B3) = 0,13, P(B4) = 0,08, прогнозирования ГМП P(B5) = 0,04, P(B6) = 0,02, P(B7) = 0, Величина доверительного интервала I(1, 2) для рассчитанных P(Bi) определялась по формулам (7) и (8):
P( Bi ) (1 P( Bi )) 1 P( Bi ) t, (7) n P( Bi ) (1 P( Bi )) 2 P( Bi ) t. (8) n Для отнесения прогнозного значения к i-й градации возмущенности ГМП используется вероятностный метод, заключающийся в том, что вероятность попадания прогнозного значения в i-ую группу рассчитывается по формуле Байеса (9):
P( Bi ) PB ( A) PA ( Bi ), i (9) P( B1 ) PB ( A) P( B2 ) PB ( A)... P( Bn ) PB ( A) 1 2 n где P ( Bi ) – априорная вероятность гипотезы;
PB (A) – условная вероятность i события А при выполнении гипотезы P( Bi ).Условная вероятность рассчитывается по формуле (10):
( Н x) k гр2 ПРi e 2 dx, PB ( A) (10) 2 k i гр где – среднеквадратическое отклонение, kгр1 и kгр2 – нижняя и верхняя границы i-й градации возмущенности ГМП.
После расчета по формулам (7) – (10) выбирается градация возмущенности, вероятность которого максимальна, и делается вывод о прогнозируемом состоянии ГМП на -й шаг. Полученная оценка сохраняется даже после реального уменьшения амплитуды возмущений в течение времени mT, где T-период главной фазы бури, m – число, m=1...p. Это необходимо для эффективного прогноза наложенных друг на друга ГБ, следующих одна за другой.
Среднесрочный прогноз 0, P состояния ГМП основан на 0, 27-дневной повторяемости ГБ.
Рассчитаны вероятности 0, кратной повторяемости Кп на 0, статистике 14 лет. Полученная 0, зависимость приведена на 0, рисунке 8.
Проведен анализ КП распределения плотности 1 2 3 4 5 Ри вероятности ошибки сунок 8 – Зависимость вероятности повторения ГБ от прогнозирования. Рассчитаны кратности их появления параметры полученного распределения. Результаты позволяют определять ошибки прогнозирования ГБ первого и второго рода.
Разработана информационно-измерительная система оценки и прогнозирования состояния ГМП, реализующая описанный выше метод (рисунок 9).
БД Г ПИ ДП УС WС Uпит Нх ПО ПО БКП U ОиП РиО БОП Iк WК Nх АЦП ДТ МК ПК Ny БД 2 АЦП Нy Рисунок 9 – Структурная схема системы оценки и прогнозирования состояния ГМП На рисунке 9 БД1, БД2 – блоки датчиков, имеющих одинаковые структуры;
ПИ –плата измерения;
HX, HY – северная и восточная составляющие вектора напряженности МПЗ;
IK – ток компенсации;
ДП – датчик поля, содержащий магниторезистивный мост, обмотку компенсации WК и катушку сброса WС;
Г– генератор импульсов сброса;
УС – усилитель сброса;
ДТ – датчик тока;
БКП – блок компенсации поля;
МК – микроконтроллер;
БОП – блок оцифровки и передачи;
ПК – персональный компьютер;
ПО ОиП – программное обеспечение оценки и прогнозирования;
ПО РиО – программное обеспечение регистрации и отображения. Система оценки и прогнозирования работает следующим образом.
При воздействии ГМП на магниторезистивный мост в ДП на его выходе появляется напряжение U0, которое поступает на БКП. Усилитель сброса создает размагничивающие импульсы тока в катушке сброса WСБ магнитоизмерительного преобразователя. В БКП происходит преобразование напряжения в компенсирующий ток, который, проходя через ДТ поступает на компенсирующую обмотку WК. В результате чего суммарное магнитное поле, создаваемое ГМП и током в компенсирующей обмотки, стремится к нулю, а напряжение на ДТ пропорционально измеряемому МП. Аналого-цифровой преобразователь формирует из полученного напряжения цифровой код (Nх и Ny), который передается на МК. МК передает цифровой код с обоих каналов в ПК, где происходит вычисление ГГМП Н по следующему выражению:
Н ( Н х )2 ( Н y )2.
В ПО ОиП, а также ПО РиО осуществляется оценка, прогнозирование и отображение состояния ГМП с применением предложенного метода. Вся полученная и обработанная информация накапливается в памяти ПК.
В третьей главе «Разработка математических моделей регулярных составляющих геомагнитного поля» с помощью спектрального анализа были получены математические модели изменения напряженности Нi ГГМП в зависимости от времени суток t. Целью анализа является разложение ряда на функции синусов и косинусов различных частот, для определения тех, появление которых особенно существенно и значимо.
В общем виде нелинейное уравнение регрессии может быть представлено:
m H i (a j cos(2 к ) t b j sin( 2 к ) t ), j где аj, bj – коэффициенты регрессии модели Фурье;
к – частота циклов за единицу времени.
В результате спектрального анализа были получены ММ для выделенных в главе 2 ССВ ГГМП. Таким образом, модель изменения напряженности ССВ ГГМП НЗИМ в декабре, январе и феврале имеет вид:
Н (t ) ЗИМ. 0,7 sin( 2 0,042t ) 1,79 sin( 2 0,083t ) 1,42 cos(2 0,125t ) В марте, октябре, ноябре ССВ ГГМП НОСН описывается:
Н (t ) ОСН. 5,43 cos(2 0,042t ) 4,32 sin( 2 0,083t ) В апреле, мае, июне, июле, августе и сентябре ССВ ГГМП НЛЕТ описываются уравнением:
Н (t ) ЛЕТ 6,26 sin( 2 0,042t ) 7,32 sin( 2 0,042t ) 8,65 cos(2 0,083t ) Проведена проверка значимости полученных коэффициентов уравнения регрессии при помощи t-критерия Стьюдента, незначащие коэффициенты были отсеяны. Проведена проверка адекватности всех трех моделей экспериментальным данным с применением критерия Фишера.
Построена ММ изменения напряженности ССВ ГГМП Нi в течение суток c использованием ортогональных полиномов Чебышева, позволяющих значительно упростить вычисления, так как при каждом повышении степени полинома, добавление новых слагаемых не изменяет вычисленных ранее коэффициентов:
H 0 0 ( x) 11 ( x)... n n ( x), где 0 ( xi ) 1 – полином нулевого порядка;
1 ( xi ) ( x 1 ) 0 ( x) 1 – полином первого порядка;
2 ( xi ) ( x 2 ) 1 ( x) 2 – полином второго порядка;
n1 ( x) ( x n1 ) n ( x) n1 n1 ( x) – полином порядка;
n-го x ( x ), n1, i i n i 2 – коэффициенты xi n1 ( xi ) n ( xi ) H (x ) n1 i n i ( x ) n ( xi ) n1 ( xi ) 2 n i связи.
Были проведены исследования необходимой и достаточной степени применяемого полинома Чебышева для аппроксимации графика изменения напряженности магнитного поля Н в течение суток.
Было установлено, что выделенные посредством разработанных ММ возмущения ГМП НВОЗМ имеют нормальный закон распределения.
Четвертая глава «Разработка и исследование информационно измерительной системы оценки и прогнозирования влияющих на организм человека возмущений геомагнитного поля» посвящена созданию системы, реализующей обобщенную структурную схему, предложенную в главе 2, а также предложенный метод выделения влияющих на организм человека возмущений ГМП. В состав системы входят: двухкомпонентный датчик, который регистрирует северную Нх и восточную компоненты Нy вектора напряженности магнитного поля Земли, со встроенными катушкой компенсации и катушкой сброса;
блок компенсации поля;
генератор и усилитель импульсов сброса;
датчик тока;
плата сбора информации Е–24;
состоящая из измерительного блока с четырьмя 24 разрядными АЦП и 8-ми битным процессором AVR фирмы ATMEL;
персональный компьютер.
Внешний вид разработанной системы приведен на рисунке 10. В качестве датчиков магнитного поля использованы магниторезистивные преобразователи фирмы KMZ Phillips. Оцифровка сигнала с платы измерения осуществляется посредством модуля E-24 фирмы L Card. Связь аппаратной части с компьютером осуществляется через интерфейс RS-232. На персональном компьютере в среде графического Рисунок 10 – ИИС оценки и прогнозирования программирования Labview состояния ГМП реализовано программное обеспечение для обработки данных и отображения результатов. В программном обеспечении реализованы: расчет ГГМП, выделение возмущений посредством предложенного метода, оценка текущего уровня геомагнитной активности по амплитуде возмущений, определение уровня ГБ, краткосрочный и среднесрочный прогнозы, расчет ошибок прогнозирования, определение точности и достоверности прогноза, сохранение результатов измерений и обработки с привязкой к местному времени, отображение всех рассчитанных параметров на экране персонального компьютера.
Основные технические характеристики: количество одновременно измеряемых проекций магнитного поля: 2;
диапазон измеряемого магнитного поля: ±200 мкТл;
погрешность измерения: не более ±4 %;
количество градаций возмущенности магнитного поля Земли: 7;
потребляемая мощность платы измерения: 1,5 Вт.
Все проведенные исследования экспериментально подтвердили высокую эффективность применения предложенной модели и разработанной системы для оценки и прогнозирования состояния ГМП.
Проведены клинические испытания разработанной системы. Установлено, что для более 60 % пациентов с заболеваниями сердечно-сосудистыми наблюдается зависимость их состояния от параметров ГМП.
Разработанное устройство принято к внедрению в учебный процесс кафедры «Информационные и измерительные системы и технологии», в СП ООО «А мобайл» (г.Сухум), ООО «Стальмаркет» (г.Краснодар).
Разработан ряд программных продуктов, предназначенных для обработки результатов измерений и последующего прогнозирования.
Разработанная информационно-измерительная система может быть использована в качестве дополнительного модуля в системах оценки функционального состояния организма человека.
ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ 1 На основании проведенного обзора и анализа установлено влияние ГМП на организм человека и малая эффективность прогнозов ГБ, существующих в настоящее время. Показано, что перспективным направлением повышения достоверности прогноза ГБ является его осуществление на основе измерения и обработки параметров вектора напряженности МПЗ. Установлена необходимость разработки методов прогнозирования ГМП и портативных устройств, позволяющих людям, чувствительным к геомагнитным воздействиям, и специализированным медицинским организациям, заблаговременно получать информацию о наступлении ГБ и принимать необходимые меры.
2 Предложен и разработан новый метод оценки состояния ГМП, основанный на выделении возмущений ГМП посредством математической модели ССВ ГГМП. Возмущения ГМП определяются как разность между измеренным значением ГГМП и значением напряженности, полученной с помощью математической модели. Это позволяет осуществлять оценку и последующий прогноз состояния ГМП, основываясь только на измерении напряженности МПЗ.
3 Построены и исследованы математические модели регулярных вариаций ГМП. Установлено, что ССВ ГГМП имеют сильную зависимость от времени года.
Поэтому математическая модель изменения ССВ ГГМП в течение календарного года может быть построена на основе системы из трех уравнений. Перове уравнение описывает январь, февраль, декабрь, второе – март, октябрь, ноябрь, а третье – апрель, май, июнь, июль, август и сентябрь. Разработанные математические модели описывают с высокой степенью корреляции ССВ и позволяют реализовать предложенный метод выделения возмущений ГМП.
4 Разработаны алгоритмы и созданы программы, позволяющие по выделенным возмущениям ГМП осуществлять краткосрочный прогноз (до часов) состояния ГМП.
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 Междунар. науч. конф. : в 10 т. / Сарат. гос. техн. ун-т - Пенза : Пенз. гос. технол. акад., 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 : материалы регион. науч.-техн. конф. студентов, аспирантов и молодых ученых вузов Ростовской области / Юж.-Рос. гос. техн. ун-т (НПИ) - Новочеркасск : ЮРГТУ, 2011. - С. 188-189.
13 Алексанян Г.К. Спектральный анализ Фурье временных рядов годовых вариаций геомагнитного поля / Г.К. Алексанян // Моделирование. Теория, методы и средства : материалы XI Междунар. науч.-практ. конф., Новочеркасск, 31 марта г. / Юж.-Рос. гос. техн. ун-т (НПИ) - Новочеркасск : ЮРГТУ, 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 Междунар. науч.-практ. конф., г. Новочеркасск, апр. 2012 г. / Юж.-Рос. гос. техн. ун-т (НПИ). – Новочеркасск: ЮРГТУ(НПИ), 2012. – С.4-7 с.
21 Свид. об. оф. рег. прогр. для ЭВМ Российская Федерация Выделение регулярных вариаций геомагнитного поля (GMP 1.0) / Г.К. Алексанян, Н.И.
Горбатенко [и др.];
Роспатент. – № 2012611370;
заявл. 27.02.2012.
22 Свид. об. оф. рег. прогр. для ЭВМ Российская Федерация. Построение математической модели регулярных вариаций геомагнитного поля (MAT_GMP 1.0) / Г.К. Алексанян, Н.И. Горбатенко [и др.];
Роспатент. –№ 2012611337;
заявл. 27.02.2012.
23 Свид. об. оф. рег. прогр. для ЭВМ Российская Федерация Определение количества значимых слагаемы уравнения модели (SLAG_MOD 1.0) / Г.К. Алексанян, Н.И. Горбатенко [и др.];
Роспатент. – № 2012611384;
заявл. 27.02.2012.
24 Свид. об. оф. рег. прогр. для ЭВМ Российская Федерация Программа краткосрочного прогнозирования активности геомагнитного поля (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.
Формат 6084 1/16 Бумага офсетная. Ризография.
Усл. печ. л. 1,0. Уч.-изд. л. 1,6. Тираж 110 экз. Заказ 48-5146.
Южно-Российский государственный технический университет (Новочеркасский политехнический институт) Адрес университета:
346428, г.Новочеркасск, ул.Просвещения, Отпечатано в ИД «Политехник» 346400, г. Новочеркасск, ул. Первомайская, 166. тел 22-70-