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

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

Лушавин Андрей Петрович

ПРОГНОЗИРОВАНИЕ ХАРАКТЕРИСТИК ТЕХНОЛОГИЧЕСКИХ

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

Специальность: 05.13.01 - Системный анализ, управление и обработка информации (промышленность)

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

1 7 МАЙ 2012

Самара-2012

005044127

005044127

Работа выполнена на кафедре «Вычислительная техника» Федерального государственного бюджетного образовательного учреждения высшего профессионального образования «Самарский государственный технический университет»

Научный Заслуженный деятель науки РФ,

руководитель: доктор технических наук, профессор

Семенов Владимир Семенович

Официальные Батишев Виталий Иванович,

оппоненты: доктор технических наук, профессор

ФГБОУ ВПО «Самарский государственный технический университет», г. Самара, зав. кафедрой «Информационные технологии»

Небоженко Виктор Александрович,

кандидат технических наук ОАО «Институт «Гипровостокнефть», г. Самара, главный специалист отдела «Информационные технологии»

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

организация: бюджетное образовательное учреждение

высшего профессионального образования «Самарский государственный архитектурно-строительный университет» г. Самара

Защита состоится « 28 » мая 2012 г. в 11 часов на заседании диссертационного совета Д 212.217.03 ГОУВПО «Самарский государственный технический университет» по адресу: г. Самара, ул. Галактионовская, 141, корпус 6, аудитория 33.

С диссертацией можно ознакомиться в библиотеке ФГБОУ ВПО «Самарский государственный технический университет» по адресу: г. Самара, ул. Первомайская 18.

Отзывы на автореферат (в двух экземплярах, заверенные печатью) просим направлять по адресу: 443100, г. Самара, ул. Молодогвардейская, 244, СамГТУ, главный корпус, на имя ученого секретаря диссертационного совета Д 212.217.03.; факс: (846)278-44-00

Автореферат разослан «27» апреля 2012 г.

Ученый секретарь

диссертационного совета Д 212.217.03 к.т.н., доцент

убанов Н.Г.

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

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

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

Систематизированное научно обоснованное развитие прогнозирования осуществляется с первой половины 50-х годов прошлого века. Ключевые вопросы разработки и исследования методов прогнозирования разрабатывались российскими учёными Матросовым Н.Д., Вишнёвым С.М., Панасюк Б., Мосиным В.О., Яновым Я.И., Адирим И.Г., Сутягиным В. В развитие методологии прогнозирования большой вклад внесли отечественные и зарубежные учёные И.В. Бестужев-Лада, Л. Клейн, В. Гольдберг.

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

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

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

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

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

Дня достижения поставленной цели необходимо решение следующих задач.

1. Исследование временных рядов, характеризующих технологические процессы нефтяной, газовой, химической и нефтехимической промышленности;

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

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

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

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

Основными методами исследования являются методы теории множеств, теории вероятностей, математической статистики, линейной алгебры, а также теории нейронных сетей и нечеткой логики.

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

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

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

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

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

На защиту выносятся следующие основные положения.

1. Алгоритм описания и прогнозирования технологических процессов, имеющих нестационарные и непериодические характеристики, основанный на сочетании методов вейлет-преобразования и нечеткой логики.

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

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

Реализация результатов работы Разработанная методика анализа временных рядов внедрена в системе технологического контроля газокомпрессорной станции (Сызранское ЛПУМГ ООО «Газпром трансгаз Самара»). Полученные научные результаты использованы в учебном процессе на кафедре «Вычислительная техника» ГОУ ВПО Самарский государственный технический университет.

Апробация работы. Основные положения работы и результаты исследования докладывались и обсуждались на одной Международной и четырёх Всероссийских научно-технических и научно-практических конференциях, а именно: Пятая Всероссийская межвузовская научно-практическая конференция «Компьютерные технологии в науке, практике и образовании». (Самара, 2006), 65-я Всероссийская научно-техническая конференция по итогам НИР университета за 2007 год. Самарский гос. арх.-строит. ун-т. (Самара, 2007), 66-я Всероссийская научно-техническая конференция по итогам НИР университета за 2008 год. Самарский гос. арх.-строит. ун-т.(Самара, 2008), 67-я Всероссийская научно-технической конференции «Традиции и инновации в строительстве и архитектуре». / СГАСУ,- Самара, 2010, Международная научно-техническая конференция «Информационные, измерительные и управляющие системы» (ИИУС-2010)./ СамГТУ - Самара, 2010.

Публикации. По теме диссертации опубликовано 12 печатных работ, в том числе 4 работы из списка журналов, рекомендованных ВАК РФ.

Структура и объём работы: Работа состоит из введения, четырёх глав, заключения, списка литературы и трёх приложений. Основная часть изложена на 158 страницах, содержащих 48 рисунков и 5 (пять) таблиц.

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

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

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

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

2. Анализ временных рядов необходимо проводить в два этапа. На первом требуется выполнить отделение регулярной составляющей от случайных, а на втором - исследовать нерегулярные составляющие.

Во второй главе выполнен анализ методов описания и прогнозирования временных рядов и предложен алгоритм анализа и прогнозирования технологических процессов, имеющих нестационарные и непериодические характеристики, отличающийся сочетанием методов вейвлет-преобразования и нечеткой логики, что позволяет более точно оценить характеристики процессов. Рассмотрены классические методы Фурье, ARIMA и экспоненциального сглаживания, а также методы математической статистики. Самый известный их них - ARIMA использует алгоритмы, разработанные Г.Е.П. Боксом (G.E.P. Box) и Г.М. Дженкинсом (G.M. Jenkins) [4]. Существует целая иерархия моделей Бокса-Дженкинса, которую можно представить в виде следующей схемы:

AR(p) + MA(q) -> ARMA{p,q) -»■ ARMA(p,q)(P,Q) -» ARIMA(p,q,r)(P,Q,R) ->.......> (])

где AR(p) - авторегрессионная модель порядка p.

Компьютерные эксперименты, выполненные с помощью системы STATISTIC А, результаты которых приведены на рисунке 1, показали ограниченность применения метода ARIMA для прогнозирования характеристик технологических процессов общего вида.

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

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

Кхказз, «ойеі: (0.1,1X0.1.1) виколої 14

ев» Мх*(0.1.1X0.1.1) ЗвзопаІ к? 12

Іпрі/1- 5Ш*ЗТХ ЙагШогідл 1 ЬйгіогідЛ ІЮ

(а) (Ь)

Рисунок 1- Результаты прогнозирования, полученные в системе прогнозирования 8ТАТ18'ПСА при применении метода АЮМА для рядов с разными характеристиками.

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

а)

б)

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

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

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

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

-1, если 0 < х < — 2

<р(х) = 1, если — < х < 1 (2)

2

О, если х < 0 или х > 1

Функция е>(*) = £г(Я) называется масштабирующей, если она может быть представлена в виде: = ^- л) п этом ХКГ <*>. Наб

1 пег

{К)называется масштабирующим фильтром. Произведём преобразование Фурье масштабирующего соотношения = т/л'"1 £сли

яеЛ

принять, что ИЛ<0) = 5 тогда = при этом функция

Иа (®) будет называться частотной функцией масштабирующей функции <?(х).

Подробное изложение теоретической методики кратномасштабного вейвлет -разложения приведено в работе [3]. Здесь приводятся лишь

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

Частотные функции удовлетворяют свойству унитарности

матрицы:

Я0(г) Я0(р2) ... Н0(р"-'2) Я,(г) Н,(рг) ... ЯДр^'г)

(,Я„.,(7) Я^Дрг) ...

(3)

_____ _

где 2 = е , р-е

Для функции Хаара масштабирующее соотношение порядка N имеет вид: <Р(х) = <р(Ых) + <р{Ых -1) +... + <р(Их -Ы + \) = £ А„<з(Л6с - и); (4)

я=0

где К = — = К-\ = 1/л/лГ.

Найдем вейвлет-функции ч/\х),...,щ"-\х). Пусть при этом

Ъ' с где ортонормированный базис в -ортонормированные

базисы подпространств вейвлетов, причём имеет место разложение в прямую сумму ортогональных подпространств V, =У0®1У„ Ф...ФИ-'/4. Пространство К входит в К, следовательно функции также

раскладываются по базису (^т, М) пространства = Ддя

любой функции «'олМ = <Р(* - ") имеем:

Л, = +■ (5)

Тогда условие ортогональности к К0 принимает вид:

Так как («^М!-ортонормированный базис, то получаем: где коэффициенты (в»)-фильтры вейвлетов .....V (дг).

Условие ортогональности вейвлетов ,¥') = 0 следующее : Х^» £» ~ °> в результате получаем систему уравнений:

1 Е«^» = к= и,...,лг -1,ЛГ- 1Д * у

(8)

Она имеет множество решений. В простейшем случае, соответствующем минимальному набору ненулевых коэффициентов, решение находится неоднозначно, т. е. решением будет любой набор из А'-У ортогональных векторов С'= в плоскости Д*, заданной

уравнением:

= 0:

(9)

Решением будет матрица порядка N с ортогональными строками, первая строка которой есть вектор А = (Ло.А.....А-,) = у^1-1.....Например, можно

использовать матрицу Адамара (в этом случае должно выполняться условие существования таких матриц, а именно, её порядок N должен быть кратен 4, кроме того значения N1, N/12 или N/20 должны быть степенью 2). Другим примером служит матрица, которая может быть получена из (9)

А =

і -1

і і

і

о о

-3

1 о о о

-N + 1

(10)

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

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

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

&<* ЬГ'о« 1гивл ХооЬ 2й$пск>*. Це1р □ ^ а а ч а л у \ & & о

26 24 22 20

и 50 100 160 200 2£

10 9

0 50 100 150 200 260 300 350

(ряд Б) - коэффициенты его вязкости, измеряемые один раз в час.

Рисунок 3. Графики исследованных в работе рядов - ряд С (вверху) и ряд Б (внизу).

Очевидно, что нижний ряд (ряд Б) имеет очень сложную динамику с наихудшими оценками по показателю Хёрста и энтропии Шеннона, что существенно влияет на возможность прогнозирования его значений. Это подтверждается и результатами, представленными в таблице 1. Таблица 1 Сравнительные данные для рядов Б и С

Ряд Показатель Хсрста Энтропия Шеннона

с 1,1870 -7.565444693555317е+005

1,2383 -7.564407934301824е+005

0,9539

И 0,4396 -1.147871 165128392е+005

0,4065 -1.145013563577244е+005

-0,3181

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

1. Выбирать уровень разложения п (для примера -4).

2. Выполнить разложение до уровня п (при 4 - на квартили);

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

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

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

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

Так, если в ряду О, приведенном на рисунке 3 и состоящем из 310 отсчетов, изменить только одно значение (например, двадцатое с величины 8,6 на 9,0), как показано на рисунке 4, то это приведет к заметному изменению графика вейвлет-коэффициентов (см. рис. 4). При этом вид исходной кривой О практически не изменился. Статистические показатели коэффициентов приведены в таблице 2.

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

Э50

ВД

Рисунок 4. Изменение графиков ряда О (вверху) и вейвлет - коэффициентов (внизу) при изменении значения одного отсчета (место изменения на значениях ряда показано стрелкой).

Таблица 2. Статистические характеристики (ряд Б) но уровням разложения

Таким образом, предлагаемая методика, основанная на частотных

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

8,6 9 8.6 9 8,6 9

1 0,3342 0,3339 0,1817 0,1804 2,275 2,275

2 0,2815 0,2819 0,1087 0,1099 2,275 2,275

3 0.2794 0,2803 0,1067 0,1080 2,275 2,275

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

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

Е*в ЇИ, Ігмоп Хось Ыв*>

* А / & & о

> ?

10 ■

о 60 100 150 200 250 300

10 20 ЭО 40 60

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

статистики

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

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

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

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

Рисунок 6. Исходные и предсказанные значения по двум моделям прогнозирования (MLP и ANFIS). Сплошная линия - исходный ряд, круги - MLP- прогноз, звёздочки - ANFIS -

система.

Среднеквадратическая ошибка для MLP- модели составила 0,5474, а для ANFIS сети - 0,2242, т.е. уменьшилась в 2,44 раза. При этом структура «чистого» персептрона (MLP) была 62-13-1 (двухслойный персептрон с 62 входами). По результатам экспериментов было выявлено, что гибридная сеть ANFIS показывает худшие результаты на длительных прогнозах для временного ряда с множеством точек перегиба и является наиболее эффективной на монотонных участках. Характеристики различных методов прогнозирования приведены в таблице 3.

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

ff = lg(K/S)lg(n/2), (И)

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

значениями расчёта показателя Херста Н считаются следующие: Я= 0,5-наблюдения независимы и эффективное прогнозирование невозможно, /£>0,5 - ряд персистентный, //<0,5 — ряд имеет свойство антиперсистентности. Таблица 3. Сравнительные характеристики различных способов прогнозирования

Метод Время вычисления Ошибка SSE % к наивному прогнозу

«Наивный» - 1,6953 100%

Экспонен циальное Сглаживание (без тренда) 1.9835 117%

Экспоненциальное сглаживание (линейный тренд) - 2,8859 170%

Экспоненциальное сглаживание (экспоненциальный тренд) - 2,7840 164%

Экспоне нциал ьное сглаживание (демпфированный тренд) - 2,0207 119,2%

ARIMA - 2,3369 137,85%

MLP (62-19-1) 225 с. 0,2504 14,77%

MLP-ANFIS(150 функций, 50 эпох) 36 с. 0,2175 12,83%

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

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

Рисунок 7. Схема общего алгоритма анализа и прогнозирования временного ряда

Заключение

Основные результаты и выводы, полученные в работе

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

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

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

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

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

Основное содержание диссертации опубликовано в следующих работах: Статьи в журналах, рекомендованных ВАК:

1. Семёнов, B.C. Выявление незначительных отклонений значений и составление частотной карты временного ряда с помощью кратномасштабного вейвлет-преобразования./ B.C. Семёнов, А.П. Лушавин// Вестник Самарского государственного технического университета. Серия «Технические науки». - Самара, 2010. - №4(27). -С. 235-238.

2. Лушавин, А.П. Численные данные и анализ применения различных методов прогнозирования временного ряда./ А.П. Лушавин// Естественные и технические науки. - Москва, 2011 - № 2, - С. 3 80-3 81.

3. Лушавин, А.П. Анализ улучшения аппроксимирующих свойств гибридной нейронной сетью временных рядов. Анализ эффективности методов прогнозирования временных рядов методами, не основанными на аппарате матемптической статистики./А.П. Лушавин // Вестник Ростовского государственного университета путей сообщения. - Ростов-на-Дону, 2011.— №2-С. 89-91.

4. Лушавин, А.П. Программа приведения временного ряда к стационарному виду посредством предобработки гибридной сетью./ А.П. Лушавин // Вестник Кузбасского государственного технического университета-Кемерово, 2011.-№4-С. 75-77.

Публикации в других изданиях:

5. Лушавин, А.П. Сравнительные эксперименты и анализ применения различных методов прогнозирования временного ряда./ А.П. Лушавин // Научный обозреватель. Научно-аналитический журнал. - Уфа, 2011.-№1 -С. 58.

6. Лушавин, А.П. Улучшение точности прогнозирования монотонных участков временного ряда посредством гибридной нейронной сети. Сравнительные компьютерные эксперименты./ А.П. Лушавин // Актуальные проблемы гуманитарных и естественных наук.- Москва, - 2011, 2010,-№2.-С. 31-33.

7. Лушавин, А.П. Причины использования нейросетевых алгоритмов в оптимизационных задачах./ А.П. Лушавин // Труды пятой Всероссийской межвузовской научно-практической конференции 16 ноября 2006г. Самарский государственный технический университет. - С. 110-113.

8. Семёнов, В.С Построение системы принятия решений./ B.C. Семёнов, А.П. Лушавин // Материалы 65-й Всероссийской научно-

техническая конференция по итогам НИР университета за 2007 год. Самарский государственный архитектурно-строительный университет. - С. 524-525.

9. Семёнов, B.C. Использование вейвлетов в нейронных сетях./ B.C. Семёнов, А.П. Лушавин// Материалы 66-й Всероссийской научно-техническая конференция по итогам НИР университета за 2008 год. Самарский государственный архитектурно-строительный университет. - С. 247.

10. Лушавин, А.П. Особенности временных рядов в описании неэкономических показателей./ А.П. Лушавин // Традиции и инновации в строительстве и архитектуре. Материалы 67-й Всероссийской научно-технической конференции./ Самарский государственный архитектурно-строительный университет. - Самара, - 2010. - С. 128-129.

11. Лушавин, А.П. Аппроксимация и сглаживание временного ряда посредством гибридной нейронной сетиУ А.П. Лушавин // Традиции и инновации в строительстве и архитектуре. Материалы 67-й Всероссийской научно-технической конференции./ Самарский государственный архитектурно-строительный университет.- Самара, - 2010. - С. 163-164.

12. Лушавин, А.П. Преимущества выбора показателя Херста в качестве критерия фрактальности временного ряда./ А.П. Лушавин // Международная научно-техническая конференция. Информационные, измерительные и управляющие системы (ИИУС-2010)/ Самарский государственный технический университет - Самара, - 2010.-С. 139-140.

Личный вклад автора в работах, написанных в соавторстве заключается в следующем: в [1] приведены примеры компьютерных экспериментов, исследована степень погрешности, которая получается при применении кратномасштабного вейвлет- анализа при исследовании временных рядов; в [8], [9] исследованы возможности улучшения аппроксимационных свойств искусственных нейронных сетей.

Автореферат отпечатан с разрешения диссертационного совета Д 212.217.03 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Самарский государственный технический университет» (протокол №6 от 23 апреля 2012)

заказ № 370. Тираж 100 экз.

Отпечатано на ризографе. ФГБОУ ВПО «СамГТУ» Отдел типографии и оперативной печати 443100 г. Самара ул. Молодогвардейская,244