CC BY
7цесса ТОР. Применение цифровой съемки процесса технического обслуживания позволяет провести дополнительную оценку в режиме оффлайн.
4. Данная методика позволяет в дальнейшем при проведении более масштабного эксперимента может привести к получению корректирующих коэффициентов норм времени при проведении ТОР автомобилей, с целью внедрения их на производстве.
СПИСОК ЛИТЕРАТУРЫ
1. Малышев А.И. Экономика автомобильного транспорта, Учебник для вузов. [Текст] /А.И.Малышев/ - М.: Транспорт, 1983. - 336 с.
2. Пятов С.Г. Автоматизация оценки надежности непараметрическими методами статистики». Материалы региональной НТК «Инновации и эффектив-ность производства». [Текст] /С.Г. Пятов/- ТюмГНГУ. - Сургут, 2006.
УДК 681.3.062
3. Пятов С.Г. Разработка норм времени при проведении технического обслуживания автотранспортных и технологических машин. [Текст] /с.Г.Пятов/ /Материалы Всероссийской НПК «Модернизация экономики и общества: региональный аспект», г. Белокуриха, 20-21 октября. Барнаул: Изд-во АлтГТУ, 2011. - С. 37-38.
Статья подготовлена при выполнении государственного контракта №16.523.11.3009, финансируемого Министерством образования и науки РФ.
д.Т.н., д. Э. н. Осипов Ю.М. - академик Международной академии информатизации, академик Академии навигации и управления движением, заведующий отделением кафедры ЮНЕСКО «Новые материалы и технологии» Томского университета систем управления и радиоэлектроники.
Пятов С. Г., аспирант Отделения кафедры ЮНЕСКО «Новые материалы и технологии» Томского университета систем управления и радиоэлектроники
АЛГОРИТМ ПРОГНОЗИРОВАНИЯ ИЗМЕНЕНИЯ БЫСТРОИЗМЕНЯЮЩЕЙСЯ ВЕЛИЧИНЫ И ЕГО МОДЕЛЬ
С.И. Клевцов
Предложена модель и алгоритм краткосрочного прогнозирования изменения быстроизменяющейся физической величины, характеризующей состояние технического объекта. Модель основана на использовании сглаживающих временных рядов. Определены особенности настройки и использования модели. Приведенные соотношения и алгоритм краткосрочного прогнозирования ориентированы на выполнение обработки данных в микроконтроллере микропроцессорного модуля в реальном времени. Моделирование показало эффективность предлагаемого подхода к прогнозированию изменения параметров объекта.
Ключевые слова: прогнозирование, алгоритм, временной ряд, сглаживание
Введение
Необходимость опережающего мониторинга состояния технического объекта для предотвращения и управления протеканием нештатных ситуаций в целях их ликвидации с минимальными последствиями делает постановку и выполнение задачи определения прогнозируемых значений параметров, а также прогнозируемой оценки состояния отдельных параметров технического объекта или их совокупности в реальном времени актуальной и своевременной. Современное развитие микропроцессорной техники позволяет реализовать эту задачу с использованием микропроцессорных модулей (ММ) в составе информационных систем сложных технических объектов, таких как системы активной безопасности и инерциальной навигации автомобиля. Оценки, реализуемые в ММ,
обеспечивают высокую динамику реакции на ситуации. Для выполнения прогнозной оценки значений параметров технического объекта с помощью ММ целесообразно использовать простые модели, позволяющие получить результат в режиме реального времени без существенной нагрузки на микроконтроллер модуля, например, модели, в основе которых лежат методы анализа и прогнозирования временных рядов [1, 2].
В настоящей статье рассмотрены модель и алгоритм прогнозирования значений быстропеременной физической величины на основе использования многократного экспоненциального сглаживания ее временного ряда. Эффективность модели и алгоритма оценивается на основе обработки данных об изменении значений проекций вектора ускорения автомобиля, измеренных с помощью
трехосевого акселерометра в реальных условиях движения.
Моделирование процесса прогнозирования изменения параметра технического объекта
При построении алгоритма прогнозирования на основе модели многократного экспоненциального сглаживания временного ряда предполагалось, что шаг дискретизации при съеме данных можно установить небольшим, таким, что изменение параметра в течение нескольких последовательных шагов будет незначительным. Это позволит повысить точность прогнозирования.
Пусть значения контролируемого параметра технического объекта у() измеряются микроконтроллерным модулем в дискретные моменты времени с постоянным шагом Ь. В
результате получим массив У = |у| измеренных значений параметра у = у(11) в точках 10Д-|,...Дп; 1 = 1М + И . Необходимо определить значение параметра у в точках Т + кИ, к = 1,К , где Т - текущее значение времени относительно которого делается прогноз на к - шагов или на промежуток времени х = кИ, называемый временем упреждения [1], а К - число, определяющее диапазон прогнозирования.
В качестве прогнозирующего ряда использовалась модель сглаживающего временного ряда первой степени вида [1, 3],
У 1) = 81 + 821,
(1)
характеризующаяся низкой сложностью и обладающая хорошими адаптивными свойствами.
Тогда прогнозное значение у в точках
Т + х :
У Т + х) =(2 + ах£т "(1 + |х)эТ2], (2)
где вт и в^ - экспоненциальные средние,
которые определяются на момент Т по формулам:
вт =аут +Рвт_1, в^авт +рвТ] 1, (3) а - постоянная сглаживания, которую необходимо подбирать, Р =1-а.
Для того, чтобы запустить процесс расчета необходимо задать начальные значения
во и вИ:
Р- 82,0 , з02] = 81,0 - — 82,0 (4)
Здесь 81,0 и 82,0 - начальные оценки коэффициентов 81 и 82 в (1).
Исследование возможности использования модели адаптивного временного ряда для прогнозирования параметра технического объекта проводилось на основе данных, снимаемых с акселерометра, установленного на автомобиле. Зависимость ускорения от времени представляет собой быстроизменяю-щуюся знакопеременную функцию с существенным разбросом значений. Шаг съема данных составлял И = 0,015625сек .
В [3] показано, что процедура прогнозирования непрерывной функции с использованием сглаживающего временного ряда имеет период адаптации, который зависит от точности задания начальных оценок коэффициентов аппроксимации 81,0 и 82,0 в выражении
(1). При этом характерно, что неточное задание начального приближения коэффициентов 8 и 8 влияет только на длительность
участка адаптации и форму прогнозирующей зависимости у() на этом участке и практически не влияет на погрешность прогнозирования за его пределами.
Исследования показали, что в рамках выбранной модели временного ряда погрешность прогнозирования за пределами участка адаптации в значительной степени определяется выбором постоянной сглаживания.
На рисунке 1 представлены результаты прогнозирования на время упреждения И экспериментального временного ряда, характеризующего изменение ускорения во времени, при различных значениях постоянной сглаживания а.
0.22 0.2
2
\ /ь 1 Щ
1 / Г
/ V
3
4
вп = 81,0 -
' 8 8.1 8.2 8.3 8.4 8.5 8.6 8.7 8.8 8.9 Время, сек
Рисунок 1 - Результаты прогнозирования при различных значениях а: график 1 - экспериментальная зависимость ускорения от времени; график 2 -прогноз при а = 0,55; график 3 - при а=0,35; график 4 - при а=0,15
0.24
0.18
0> 0.16
£014
0.12
0.1
Анализ графиков показывает, что при увеличении а прогнозирующая кривая точнее воспроизводит форму экспериментальной зависимости. Однако для каждой зависимости существует определенное значение а, определяющее границу возможностей модели в части точности прогнозирования. Погрешности, соответствующие вышеприведенным результатам и иллюстрирующие указанный вывод, представлены на рисунке 2.
' 8 8.1 8.2 8.3 8.4 8.5 8.6 8.7 8.8 8.9 9 Время, сек
Рисунок 2 - Погрешности прогнозирования при различных значениях а: график 1 -погрешность при а = 0,15; график 2 - при а=0,35; график 3 - при а=0,55
Результаты моделирования показывают, что на участках с резким ростом контролируемого параметра (рисунок 1, временной интервал от 8,4 до 8,5 сек.) меньшую погрешность прогнозирования дает модель с высоким значением постоянной сглаживания и, наоборот, в случае небольших колебаний параметра снижение погрешности прогнозирования связано с уменьшением постоянной сглаживания (например, участок от 8,0 до 8,4 сек., рисунок1).
С помощью модели можно осуществлять прогнозирование на время упреждения г = кИ > И. Результаты прогнозирования для различных г при постоянной сглаживания а=0,35 представлены на рисунке 3.
Анализ результатов моделирования показал, что графики прогноза для различных значений времени упреждения в целом повторяют форму экспериментальной кривой с запаздыванием, которое тем больше, чем больше значение времени упреждения. Для участков с длительным по времени (более 4-5И) ростом или уменьшением ускорения, например, для участка, определяемого временным интервалом от 8,4 до 8,5 сек. (рисунок 1), время запаздывания прогноза пропорционально времени упреждения (рисунок 4).
0.24
0.22
0.2
■ 0.18
д> 0.16 и о
о 0.14
0.12
0.1
0.08
Л А
2 Л_ (г?. ТУ Й л \
1\ 1 1 : #Г V */ \» / ♦ У V ^ V ч • у \!\ •5 У/А.
\ /\ 1 / \! •V ♦ » и« V;
\
Л»*» ., Г
7 Л* V- У А / 1 / г 4
Л* V 3
8.2
8.3
8.4
8.5 8.6 8.7 Время, сек
8.9
Рисунок 3 - Результаты прогнозирования при а=0, 35 и различных значениях г : график 1 - экспериментальная зависимость ускорения от времени; график 2 -прогноз для г = 3И; график 3 -для г = 5И; график 4 - для г = 1011
14
12
о о.
= „Ю _ £
5 16
5 «
а = 0,05
а = 0,15 Я >
/ Ж /
у *' / / л " а-0,65
У ' ^ / /у»
•■''¿у 0,55
0 2 4 6 8 10 12 Время упреждения ( в единицах II)
Рисунок 4 - Зависимость времени запаздывания прогноза от времени упреждения в единицах Л при различных значениях постоянной сглаживания
Погрешности прогнозирования на участках с длительным по времени ростом или уменьшением ускорения максимальны для графика, соответствующего г = 10И, а на участках с незначительно меняющимся трендом погрешности прогнозирования для того же графика в целом меньше по сравнению с другими прогнозирующими кривыми, соответствующими т< 10И (рисунок 5).
Следует отметить, что время запаздывания прогноза уменьшается с увеличением постоянной сглаживания а. Однако, использование а > 0,5 - 0,55 для снижения времени запаздывания не эффективно (рисунок 6).
0.06
0.04
9
0.02
0
-0.02
-0.04
0.06 0.04
л
8 0.02 х
ф 0 а
° -0.02
сц
пз
I -0.04
§ -0.06
о
ю
< -0.08 -0.1
1
А /1 Л М
\ 1 1 л \ \ / Д ' 11 \ ' г' V»' 4\ я ; * / / Щ
V I \У ¥ \ \ V 1 ) ; 1 у1\А Ч Л А //] - /
V у 1 V у \ & Т/. >'
1 \ / / \ \/ » А
V \ V / /V
у
8.2
8.3
8.4
8.5 8.6 8.7 Время, сек
8.8
8.9
Рисунок 5 - Погрешности прогнозирования при а=0,35 и различных значениях г : график 1 - прогноз для г = 3И; график 2 - прогноз для г = 5И; график 3 - прогноз для г = 1011
*__
т-41
-- .........
' - а — ------ -
0,05 0,1 0,15 0,2 0,25 0,3 0,35 0,4 0,45 0,5 0,55 0,6 0,65 0,7 Постоянная сглаживания
Рисунок 6 - Зависимость времени запаздывания прогноза от постоянной сглаживания при различных значениях времени упреждения
Алгоритм реализации процесса прогнозирования изменения параметра технического объекта
Исследование особенностей модели позволили синтезировать обобщенный алгоритм ее настройки и выполнения прогнозирования:
1. В соответствии с соотношением (1) определяются начальные оценки а1,о и §2,0 значений коэффициентов а и а2 . Для этого используются значения параметра у0 ) и
у1 (^ ) , зафиксированные на начальном этапе настройки модели прогнозирования.
2. Устанавливается начальное значение постоянной сглаживания из диапазона [0,15...0,25]. Величина а, выбранная из указанного диапазона, может служить точкой отсчета для подбора наиболее подходящего значения а с точки зрения дальнейшей минимизации погрешности прогнозирования контролируемого параметра у в процессе настройки модели.
3. После определения а, а^ и а2,о вычисляются Бо и вИ по формуле (4).
4. Для каждого Т е [^Дз,...,^] при фиксированном а выполняется прогнозирование на время упреждения г , т.е. прогнозирование на момент времени Т + г. Прогнозное значение ~(Т + т) вычисляется по формулам (3) и (2) соответственно. В процессе прогнозирования выделяется участок адаптации, когда прогнозное значение (Т + т) значительно отличаются от измеряемых значений у(Т + т) . Участок прогнозирования и, соответственно, окончание участка адаптации, можно определить по поведению прогнозирующей кривой, которая на этом участке стремиться приблизиться к зависимости у(Т + т). При этом погрешности прогнозирования в среднем практически не меняются или меняются достаточно медленно по сравнению с динамикой, фиксируемой на участке адаптации.
5. После определения участка прогнозирования необходимо минимизировать погрешности прогнозирования с помощью изменения параметра сглаживания а . Увеличивая или уменьшая а, в зависимости от предполагаемой динамики изменения зависимости у(Т) , определим минимальное (ориентировочно) значение погрешности на участке прогнозирования. Параметр сглаживания а необходимо уменьшать, если прогнозируется медленно меняющаяся зависимость и увеличивать в противном случае. Полученное значение а может быть использовано для дальнейшего прогнозирования значений у(Т) с помощью модели сглаживающего временного ряда у (Т) .
Выводы
Таким образом, разработана модель и алгоритм краткосрочного прогнозирования изменения быстроизменяющейся физической величины, характеризующей состояние технического объекта. Определены особенности настройки и использования модели, в частности, для обеспечения хорошей точности прогнозирования необходимо выбрать постоянную сглаживания, соответствующую динамике прогнозируемого процесса, для устранения или максимального сокращения участка адаптации временного ряда и расширения участка прогнозирования следует точно задавать начальные значения коэффициентов аппроксимации исходной зависимости. Также определены ограничения модели.
9
Приведенные соотношения и схема краткосрочного прогнозирования ориентированы на выполнение обработки данных в микроконтроллере микропроцессорного модуля в реальном времени. Моделирование показало эффективность предлагаемого подхода к прогнозированию изменения параметров объекта. Результаты работы могут быть полезны при разработке программно-алгоритмических средств систем мониторинга и обеспечения безопасности технических, например, систем активной безопасности и инерциальной навигации автомобиля.
Результаты исследований изложенные в данной статье получены при финансовой поддержке Минобрнауки РФ в рамках реализации проекта "Создание высокотехнологичного производства по изготовлению информационно-телекоммуникационных комплексов спутниковой навигации ГЛОНАСС/GPS/
Galileo" по постановлению правительства №218 от 09.04.2010
СПИСОК ЛИТЕРАТУРЫ
1. Лукашин Ю. П. Адаптивные методы краткосрочного прогнозирования временных рядов. [Текст] /Ю.П. Лукашин/ - М.: Финансы и статистика, 2003. -416с.
2. Бокс Дж., Дженкинс Г. Анализ временных рядов, прогноз и управление. [Текст] /Дж.Бокс, Г. Дженкинс/ - М.: Мир, 1974. - 197с.
3. Клевцов С.И. Особенности выбора параметров настройки модели сглаживающего временного ряда для осуществления краткосрочного прогнозирования изменения физической величины. [Текст] /С.И. Клевцов// Известия ЮФУ. Технические науки. 2011. №5(118). - С. 133-138.
Доцент Клевцов С.И. тел. 8-8634-32-80-25, sergkmps@mail.ru - каф. микропроцессорных систем Технологического института Южного федерального университета в г.Таганрог
УДК 536.51: 621.791.923; 621.762.5
ИССЛЕДОВАНИЕ И МОДЕЛИРОВАНИЕ ТЕМПЕРАТУРНЫХ ПОЛЕЙ НА ФАЗОВЫХ ГРАНИЦАХ В ПОРОШКОВЫХ СМЕСЯХ, ПРЕТЕРПЕВАЮЩИХ ПЛАВЛЕНИЕ И ХИМИЧЕСКОЕ
ПРЕВРАЩЕНИЕ
В.В. Иванайский, А.В. Ишков, К.В. Селиверстов, С.Ф. Дмитриев
В статье описаны комплексные исследования по измерению и моделированию температур на границах фаз в сложных плавящихся и термореагирующих порошковых смесях. Для регистрации температуры в процессе индукционной наплавки предложено использовать ХА-микротермопарный метод и метод термоиндикации при помощи СВС-составов. Методом численного расчета нестационарного уравнения теплопроводности смоделирован профиль температур при электростимулированном СВС композиционного контактного материала.
Ключевые слова: порошковый материал, термометрия, индукционная наплавка, микротермопара, термоиндикатор, СВС, моделирование.
Введение
Современные исследования многокомпонентных порошковых смесей, при нагревании которых происходит плавление отдельных компонентов, образование сплавов и покрытий, химические реакции и синтез новых веществ, невозможны без знания величин температуры и ее распределения на границах фаз. Такая информация необходима для управления процессами индукционной наплавки (ИН), сварки, самораспространяющегося высокотемпературного синтеза (СВС) и многими другими технологиями.
Важной составляющей процесса упрочнения деталей методом индукционной
наплавки является нагрев частиц твердого сплава и флюса в наплавочной шихте [1]. В то же время регистрация с приемлемой точностью температуры в процессе ИН имеет определенные технические и методические трудности. Во-первых, толщина покрытий, формируемых методом ИН колеблется в пределах 1 -3 мм, поэтому использование для этого обыкновенных промышленных и исследовательских термопар становится невозможным. Во-вторых, в процессе индукционного нагрева поверхность основного металла излучает большее количество энергии (как в виде теплоты, так и лучистым способом) чем формируемое покрытие, поэтому пирометри-
