Авторефераты по всем темам  >>  Авторефераты по техническим специальностям Куприянов Александр Сергеевич

МОДЕЛИРОВАНИЕ РЕАКЦИИ СЛОЖНОЙ АДАПТИВНОЙ СИСТЕМЫ НА ИМПУЛЬСНОЕ ВОЗДЕЙСТВИЕ

Специальность 05.13.18 - Математическое моделирование, численные методы и комплексы программ

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата технических наук

Санкт-Петербург 2012

Работа выполнена в ФГБОУ ВПО Петербургский государственный университет путей сообщения на кафедре Высшая математика.

доктор физико-математических наук, профессор,

Научный консультант:

Фоменко Виктор Николаевич

Официальные оппоненты: доктор технических наук, профессор, профессор кафедры Математика и моделирование ФГБОУ ВПО Петербургский государственный университет путей сообщения Дегтярев Валентин Григорьевич кандидат физико-математических наук, доцент кафедры Прикладная математика ФГБОУ ВПО Санкт-Петербургский государственный политехнический университет Мясникова Екатерина Марковна Ведущая организация ФГБОУ ВПО Балтийский государственный технический университет Военмех имени Д. Ф. Устинова

Защита состоится 13 декабря 2012 г. в 13 ч. 30 мин. на заседании диссертационного совета Д218.008.06 на базе Петербургского государственного университета путей сообщения по адресу: 190031, Санкт-Петербург, Московский пр., д. 9, ауд. 1-217.

С диссертацией можно ознакомиться в библиотеке Петербургского государственного университета путей сообщения.

Автореферат разослан 13 ноября 2012 г.

Ученый секретарь диссертационного совета, кандидат технических наук, профессор Кудряшов Владимир Александрович

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность работы. В настоящее время математическое моделирование применяется практически во всех областях науки и техники.

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

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

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

В диссертации затронуты следующие области знания: обработка изображений, цифровая обработка сигналов (ЦОС), классификация объектов.

Обработке изображений посвящены работы Д. Форсайта, Ж. Понса. ЦОС находит очень широкое применение в промышленных и бытовых приборах, выпускаются специализированные процессоры (DSP-процессоры), оптимизированные для ЦОС. Среди фундаментальных работ можно отметить книги А.Б. Сергиенко, Т.К. Муна, В.С. Стирлинга, А.В. Опенгейма, Р.В. Шафера, Б. Боасхаша. Одним из разделов ЦОС является анализ временных рядов; этой теме посвящены работы Льюис-Бека, Ф. Хёппнера, Х.С. Ганди.

Классификации объектов посвящены работы таких ученых, как К.А. Куликовски, Ю. Шурман, В.Н. Вапник, К. В. Воронцов.

Система обработки результатов обследования, используемая до настоящего времени, имела ряд недостатков. Выделение объекта на полутоновом изображении до недавнего времени осуществлялось в полуавтоматическом режиме, с участием оператора. Размеры объекта оценивались по высоте и ширине описанного прямоугольника, точность такой оценки была низкой. Параметры временного ряда (зависимости площади зрачков от времени) оценивались на основе позиций и значений минимумов и максимумов самого временного ряда и его производной. Такие оценки являются очень неустойчивыми и зависят от незначительных колебаний, шума сигнала. Наконец, решение о принадлежности временного ряда к нормальной выборке принималось на основании нормы (групповой или индивидуальной):

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

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

Предметом исследования являются: модель адаптивной системы, методы классификации объектов, а также алгоритмы и методы цифровой обработки сигналов.

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

На защиту выносятся:

1. Практически реализованная модель реакции зрачководвигательной системы человека на кратковременную световую вспышку.

2. Классификация пациентов по полученным оценкам параметров построенной модели.

3. Методика построения временного ряда по последовательности цифровых изображений.

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

5. Программная реализация модели реакции зрачководвигательной системы человека.

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

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

Научная новизна:

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

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

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

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

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

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

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

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

- повысить точность измерения объекта (зрачка) на изображении в 10-(в зависимости от диаметра зрачка) раз по сравнению с ранее использовавшимся методом;

- автоматически выделять объект на изображении;

- оценивать его форму;

- моделировать пупилломоторную реакцию;

- оценивать степень рассогласованности сокращения зрачков;

- классифицировать обследуемых с помощью факторного анализа, в частности, повысить вероятность выявления обследуемых в состоянии наркотической интоксикации (на имеющейся выборке этот показатель составил 100%).

На основании предложенных в диссертации методик и алгоритмов разработан программный комплекс, который применяется в НИИ Гигиены, профпатологии и экологии человека при Федеральном медико-биологическом агентстве (НИИ ГПЭЧ).

Внедрение результатов. Программное обеспечение, разработанное в ходе диссертационного исследования, применияется в НИИ ГПЭЧ. Кроме того, разработанное программное обеспечение зарегистрировано в Реестре программ для ЭВМ под названием Программа для моделирования временных рядов дифференциальными уравнениями, свидетельство № 2012613235.

Апробация результатов. Основные научные и практические результаты настоящей диссертационной работы доложены, обсуждены и одобрены:

- на совместном заседании секции № 4 Проблемной комиссии ФМБА России, Санкт-Петербургского отделения Всероссийского общества токсикологов и Ученого Совета ФГУП НИИГПЭЧ ФМБА России Вопросы обеспечения химической безопасности в Российской Федерации (СанктПетербург, 2007);

- на конференции Вопросы обеспечения химической безопасности в Российской Федерации (Санкт-Петербург, 2007);

- на LXXII научно-практической конференции студентов, аспирантов и молодых ученых Транспорт: проблемы, идеи, перспективы. Неделя науки - 2012 в Петербургском государственном университете путей сообщения, (Санкт-Петербург, 2012).

Публикации. По результатам выполненных диссертационных исследований опубликованы девять печатных работ, из них три статьи в изданиях, рекомендованных ВАК, и одна статья - в сборнике материалов международной конференции SWorld.

Структура и объем работы. Диссертационная работа представлена введением, четырьмя главами и списком литературы, при этом составляет 1страницы основного текста, содержит 53 рисунка и 5 таблиц.

Библиографический список составляют 68 наименований.

СОДЕРЖАНИЕ РАБОТЫ

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

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

Там же дан обзор современных достижений в области обработки изображений, анализа сигналов и классификации объектов. Рассмотрены:

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

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

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

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

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

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

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

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

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

В диссертации решающее значение придается двум аспектам геометрических характеристик объекта: точности и информативности. Точность характеристики при прочих равных условиях зависит от количества пикселей (то есть, от шага квантования), участвующих в подсчете; так например, оценка площади объекта будет гораздо точнее (шаг квантования - меньше), чем оценка ее протяженности, так как в первом случае в подсчете участвуют все пиксели объекта, а во втором случае - только крайние. Информативность характеристики показывает, насколько хорошо она описывает свойство объекта применительно к данной предметной области; как правило, специфические характеристики (например, характеристики, описывающие форму объекта), более информативны, чем общие.

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

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

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

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

Моделирование временных рядов можно провести одним из следующих способов:

1) аппроксимацией временного ряда аналитически заданной функцией (в частности, решением дифференциального уравнения);

2) линейным разложением временного ряда на составляющие (одним из таких разложений является факторный анализ);

3) непосредственным вычислением различных характеристик временного ряда (например, положений и значений его экстремумов).

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

Рисунок 2 - Моделирование временного ряда путем минимизации невязки Так как невязка как функция параметров обычно многомодальна, может быть использован алгоритм глобальной оптимизации на основе алгоритма Нелдера-Мида. В рассматриваемом случае алгоритм был модифицирован:

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

Моделирующая функция должна удовлетворять четырем критериям:

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

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

x''+20x'+0 (x - Dn) = v0 (t -0).

В качестве переменной х выступает площадной диаметр зрачка. Параметры данной модели имеют физическую интерпретацию, в отличие от использовавшейся ранее кусочно-линейной модели. Предполагается, что в момент времени 0 система получает мгновенную скорость v0. Зрачок имеет некий ненулевой диаметр Dn (он вычисляется по первым отсчетам до вспышки, а не рассчитывается с помощью алгоритма минимизации). Сама пупилломоторная система представляет собой систему с собственной частотой 0 и коэффициентом демпфирования . Искусственная задержка 0 связана с тем, что скорость прохождения импульса через нервные волокна конечна.

Наконец, параметр v0 характеризует глубину реакции системы. Ниже показаны три типичных случая - зрачковые реакции здорового обследуемого, пупиллограмма наркозависимого в состоянии абстинентного синдрома и наркозависимого в состоянии интоксикации. Видно, что отличаются не только начальные диаметры, но и время задержки, начальная скорость и глубина реакции. В первом случае параметры модели таковы: Dn=5.87, 0 =0.256, 0=1.76, =1.25, v0=5.74, во втором - Dn=2.82, 0 =0.287, 0=2.09, =1.30, v0=3.46, в третьем Dn= 7.60, 0 =0.306, 0=2.88, =1.0, v0=2.34.

Рисунок 3 - Типичная пупиллограмма здорового обследуемого и ее аппроксимация Рисунок 4 - Типичная пупиллограмма наркозависимого в состоянии абстинентного синдрома и ее аппроксимация Рисунок 5 - Типичная пупиллограмма наркозависимого в состоянии интоксикации и ее аппроксимация Данная модель имеет ряд существенных преимуществ перед использовавшейся ранее кусочно-линейной аппроксимацией. Прежняя модель была весьма нестабильна и зависела от случайных шумов, присутствующих в пупиллограмме, особенно параметр 0. Поэтому операторам нередко приходилось повторять обследование, а неопытные операторы могли делать неправильные выводы.

Была также оценена невязка (среднеквадратичное отклонение) между пупиллограммами и моделями, средняя невязка составляет 0.09 мм. Таким образом, предложенная модель хорошо аппроксимирует пупиллограммы.

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

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

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

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

Метод факторного анализа позволяет представить величины xi в виде p xi = mi +i ik fk + ui, i =1,K,n, a k=где mi и - математическое ожидание и стандартное отклонение данного i временного сечения (всего рассматриваются n дискретных временных точек), fk - некоррелированные нормированные (Mfk = 0; Dfk =1) случайные величины (лобщие факторы), аik - факторные нагрузки, ui - индивидуальные факторы (лостатки). Если индивидуальные факторы невелики, то приведенное представление означает, что временные сечения заполняет некую область в nмерном пространстве, являющуюся частью линейного многообразия, размерность которого равна количеству общих факторов p.

Факторные нагрузки определяются соотношениями aik = k vik.

Здесь{k}, (k = 1,K, p) - наибольшие собственные значения матрицы:

Du1 /1 0 L 0 0 Du2 / L ~ R = R 0 0 O 0 0 L Dun / n , где vik - компоненты нормированного собственного вектора этой матрицы с собственным значением k, R - матрица корреляций, Dui - дисперсии индивидуальных факторов (линдивидуальности).

Как критерий точности представления величин через общие факторы p n обычно рассматривается сумма дисперсий временных сечений D a fk, i ik i =1 k =вычисленных без учета индивидуальных факторов. Эта величина сравнивается n с полной дисперсией. Близость отношения этих величин (относительная i i=суммарная общность) к единице свидетельствует о хорошей точности модели.

В качестве иллюстрации в диссертации рассмотрены пупиллограммы как временные ряды. Каждая пупиллограмма представлена 155 временными сечениями; в качестве величины xi выступает площадной диаметр зрачков. В результате факторного анализа были выделены три общих фактора. В нашем иссследовании относительная суммарная общность оказалась равной 0,996, что свидетельствует о высокой точности предложенной модели пупиллограмм.

Непосредственно измеряемые параметры временного ряда. Помимо моделирования временного ряда, для описания временного ряда могут быть использованы непосредственно вычисляемые характеристики. К ним относятся:

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

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

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

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

В случае обработки данных пупиллометрического обследования в качестве параметров объекта используются параметры модели (способ нахождения таких параметров с помощью аналитически заданной функции описан в главе 3) и некоторые дополнительные параметры (неравенство зрачков, вытянутость зрачков). При этом в результате применения модели Важные экстремумы - термин предложенный в Harith Suman Gandhi. Important extrema of time series: theory and applications. 2004.

повысилась вероятность выявления наркозависимых в состоянии интоксикации.

Так, с помощью ранее использовавшейся кусочно-линейной модели только 98% обследуемых в состоянии интоксикации были классифицированы как обследуемые в состоянии интоксикации, то с помощью параметров модели этот показатель возрос до 100% (на той же выборке). Кроме того, параметры кусочно-линейной аппроксимации были чувствительны к небольшим изменениям пупиллограммы, это приводило к необходимости повторного обследования. Предложенная модель не имеет этого недостатка.

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

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

n u = uiu C ij j, i, j p = fk2.

f k =Здесь Cij - элементы матрицы, обратной по отношению к матрице ковариации индивидуальных факторов ui. В силу принятого выше предположения о нормальном распределении факторов статистики имеет - распределение с p и n степенями свободы соответственно.

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

n n fk = iaik ; ui =,u )wj, wi Cov(ui j i=1 j=где величины wi находятся как решения системы уравнений:

n (Cov(x, xj ))wj = xi - mi.

i j=Пусть Pu и Pf доверительные вероятности для обоих критериев, а F(r ) (x) функция - распределения с числом степеней свободы r. Тогда допустимая область критерия значимости определяется системой неравенств < Ku u, f < K f, где правые части удовлетворяют уравнениям F(n) (Ku ) = Pu ; F( p) (K ) = Pf.

2 f Данный метод не требует какой-либо априорной информации о характере временного ряда, поэтому его можно применять к широкому кругу задач.

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

первая - обследуемые, не принимающие наркотики (703 пупиллограммы), вторая - наркозависимые в состоянии интоксикации (30 пупиллограмм), третья - наркозависимые в состоянии абстинентного синдрома (20 пупиллограмм).

Первая выборка была принята за норму, по ней была рассчитана матрица ковариации. Далее были выделены три фактора при индивидуальных общностях не менее 0.99. К выборке был применен описанный выше критерий значимости с доверительной вероятностью 0.95. В результате по критерию u были отбракованы 137 кривых, а по критерию были отбракованы f кривые, по объединению критериев - 159. То есть, можно сказать, что с данной доверительной вероятностью в первой выборке 544 нормальных обследуемых и 159 обследуемых с отклонениями.

Далее данная методика была применена ко второй и третьей выборкам. В обоих случаях только по критерию u были отбракованы практически все пупиллограммы (30 из 30 в состоянии интоксикации и 19 из 20 с абстинентным синдромом). Это связано с тем, что остатки ui для данных выборок имеют не такое распределение, как остатки для первой выборки, и, как следствие, 2 величина u не имеет -распределения.

2 Рисунок 6 - Критерии u и для пупиллограмм здоровых (о) и f наркозависимых обследуемых в состоянии интоксикации (х) и абстиненции (+) Критерий позволил дифференцировать между собой выборки с f патологиями. В частности, оказалось, что значение статистики лежит в f интервале (0,84 17,11) для лиц с интоксикацией и в интервале (4,28 32,57) для пациентов с абстинентным синдромом.

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

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

Классификационный подход на основе особых точек. Автором предложен метод кластеризации временных рядов, основанный на комбинации:

- метода выделения важных точек и оценки их важности;

- расстояния Левенштейна;

- кластеризации слиянием.

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

В заключении изложены основные результаты исследования.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ 1. Разработана модель зрачководвигательной реакции.

2. Реализован алгоритм подбора оптимальных параметров модели зрачководвигательной реакции.

3. Разработана схема классификации обследуемых на основе параметров модели их зрачководвигательной реакции. Точность такой схемы классификации составила 100% для обследуемых в состоянии абстинентного синдрома, интоксикации и 91% для обследуемых, не принимающих психотропные вещества, на имеющейся экспериментальной выборке.

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

5. Разработан алгоритм выделения темных объектов, близких по форме к кругу, на светлом фоне. Данный алгоритм нашел применение для обработки данных медицинского обследования. Применение такого алгоритма позволило отказаться от ручного выделения зрачка и ускорить обследование в 2-2,5 раза.

6. Разработан алгоритм оценки размера и формы объекта на изображении.

Данный алгоритм также нашел применение для обработки данных медицинского обследования.

7. На основе результатов пп. 1-6 реализован прототип соответствующего программного обеспечения.

Список публикаций в изданиях, рекомендуемых ВАК:

1. А.С. Куприянов / Методы обработки и анализа пупиллограмм. // Известия вузов. Приборостроение. - 2009. - №8 (52). С. 58-63.

2. А.С. Куприянов, В.Н. Фоменко / Математические модели зрачковых реакций глаз человека (пупиллограмм). // Известия ПГУПС, выпуск 4, 2010.

С. 220 -230.

3. В.В. Гарбарук, В.Н. Фоменко, А.С.Куприянов / Применение факторного анализа для исследования временных рядов. // Известия ПГУПС, выпуск 3, 2011. С. 143 - 152.

Список статей, опубликованные в материалах международных конференций:

4. В.В. Гарбарук, В.Н. Фоменко, А.С. Куприянов / Классификация случайных процессов на основе факторного анализа // Сборник научных трудов SWorld. Материалы международной практической конференции Современные проблемы и пути их решения в науке, транспорте, производстве и образовании - 2011. - Выпуск 4. Том 11. - Одесса: Черноморье, 2011. - ЦИТ: 411-0164.

С. 25-26.

Список публикаций по теме диссертации:

5. Ф.А. Цимбал, А.С. Куприянов. / О повышении точности измерения размеров зрачка. // Актуальные проблемы химической безопасности в Российской Федерации. Сборник трудов научно-практической конференции под общ. ред. В.Р. Рембовского и А.С. Радилова. Санкт-Петербург, 2007.

С. 218Ц220.

6. Ф.А. Цимбал, А.С. Куприянов. / О построении математической модели зрачковой реакции. // Актуальные проблемы химической безопасности в Российской Федерации. Сборник трудов научно-практической конференции под общ. ред. В.Р. Рембовского и А.С. Радилова. Санкт-Петербург, 2007.

С. 220-222.

7. Ф.А. Цимбал, А.С. Куприянов, А.Л. Куцало. / О программной реализации экспресс-диагностики наркотической интоксикации. // Вопросы обеспечения химической безопасности в Российской Федерации. Сборник трудов совместного заседания под общей редакцией А.Н. Гребенюка и др.

Санкт-Петербург, 2007, изд-во Фолиант. С. 99-101.

8. Ф.А. Цимбал, А.С. Куприянов, Д.С. Слухова, И.И. Караченцев. / Об анализе главных компонент пупиллограмм. // Вопросы обеспечения химической безопасности в Российской Федерации. Сборник трудов совместного заседания под общей редакцией А.Н. Гребенюка и др. СанктПетербург, 2007, изд-во Фолиант. С. 102-104.

9. А.С. Куприянов / Классификация адаптивных систем по заданному эталону. // Транспорт: проблемы, идеи, перспективы. Материалы LXXII всероссийской научно-технической конференции студентов, аспирантов и молодых ученых. Санкт-Петербург, ПГУПС, 2012. С. 147-148.

Подписано к печати.11.2012 г. Печ. л. - Печать - ризография Бумага для множит. апп. Формат 60x84 1/ Тираж 100 экз. Заказ № СР ПГУПС 190031, С.-Петербург, Московский пр., Авторефераты по всем темам  >>  Авторефераты по техническим специальностям