Модель мультиагента прогнозирования автоматизированных систем управления Текст научной статьи по специальности «Математика»

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

Критерий «Целевая функция» можно пояснить следующим образом: характеристики оптимального объекта (ТП) обеспечивают при заданных ограничениях их значений достижение экстремума заданной целевой функции [5].

Технологической себестоимостью детали называется та часть ее полной себестоимости, элементы которой существенно изменяются для различных вариантов ТП. К таким изменяющимся элементам относятся: Мо - стоимость исходной заготовки; Зо - заработная плата станочника; Ао -амортизационные отчисления от оборудования; Ато - амортизационные отчисления от технологического оснащения; Ро - затраты на ремонт и обслуживание оборудования; И - затраты на инструмент; Пэ - затраты на силовую электроэнергию; Пп - затраты на амортизацию и содержание производственных площадей. Сумма остальных элементов - технологическая себестоимость обработки. Варианты оцениваются с помощью сравнения технологической себестоимости выполнения операции: С=Зо+Ао+Ат.о+Ро+И+Пэ+Пп+Ппр.

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

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

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

Литература

1. Попов П.М. Оптимизация технических решений проектирования и управления на основе экономико-математических методов анализа. Ульяновск: УлГТУ, 2000. 154 с.

2. Основы автоматизации машиностроительного производства: учеб. для машиностроит. спец. вузов / Е.Р. Ковальчук [и др.]; под ред. Ю.М. Соломенцева. М.: Высш. шк., 1999. 2-е изд. 312 с.

3. Ларин С.Н. Основные задачи обеспечения технологичности конструкции изделия в автоматизированных системах // Автоматизация управления. 2004. .№ 4. С. 62-67.

4. Информационная поддержка жизненного цикла изделий машиностроения: принципы, системы и технологии CALS /ИПИ / А.Н. Ковшов [и др.]. М.: Издат. центр «Академия», 2007. 304 с.

5. СALS (Continuous Acquisition and Life Cycle Support -непрерывная информационная поддержка жизненного цикла продукции) в авиастроении / Б.М. Абрамов [и др.]; науч. ред. А.Г. Братухин. М.: Изд-во МАИ, 2002. 260 с.

МОДЕЛЬ МУЛЬТИАГЕНТА ПРОГНОЗИРОВАНИЯ АВТОМАТИЗИРОВАННЫХ СИСТЕМ УПРАВЛЕНИЯ

(Работа выполнена в рамках проекта по аналитической ведомственной программе «Развитие научного потенциала высшей школы (2009-2010 годы)»)

В.Е. Лихачев (Рязанский государственный университет им. С.А. Есенина, vlikhachev@mail.ru)

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

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

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

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

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

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

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

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

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

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

Эвристическими методами решается большой класс задач приближения функций многих переменных. Для таких задач характерно плохо обусловленное множество координат - признаков многомерной функции. Если значения функции определяют некоторые классы Шъ k=1,2,..., m, то говорят о задаче классификации, или распознавания образов. Для задач распознавания образов характерно большое признаковое пространство вектора x = (xn,xl2,...,xln), /=1, 2,..,р, и малое количество строк матрицы A = | |xln||, для которой

известна принадлежность классам шк номеров l. Здесь А - обучающая матрица.

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

Рассмотрим в n-мерном пространстве (n>1) следующую задачу: заданы значения функции F(x) в точках x е Rn; требуется построить полином наилучшего приближения [2]:

Tpn(xii, xi2 ,...,xta), i=1,2,..., p, (1)

для которого

S[F(Si)-Tpn(x,)]2 ^ min. (2)

i=1

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

Для тригонометрического полинома (1) имеет место выражение

TV1,V2,...,v,(x) = £ Ck^knexp^x! +... + knxn)], (3)

-vi iki iVi

где у1,...,уп - натуральные числа; хх,х2,...,хп -комплексные переменные; Ск1,...,кп - постоянные коэффициенты. Однако реализация вычислительного процесса на основании (3) затруднительна. Рассмотрим построение решения (1) в виде Т(Ю = 1, 1е [а, Ь],

pn V

v L

SSCij(x)9ij(t) = U(x,t),

i=ij=i

(4)

где U(x,t) = max в точке t; {фу(0} - ортогональная на [a, b] система функций, для которых имеют место v непересекающихся классов. Предполагается, что значения функции F(x) распределены на [a, b].

Если F(x)=s, s=1,2,..., m, то разделяющей функцией для (4) будет

v L

Ds(x) = S SCij(x)9ij(s), se [a,b]. (5)

i=ij=i

Следовательно, задача приближения функций многих переменных и распознавания образов в случае (4) сводится к определению полиномов: j Cij(S). (6)

Точное решение в этом случае усложняется, так как необходимо построить vL n-мерных функций. В то же время известно, что для коэффициентов Фурье существует такое 8>0, что из p(Cj,c ij) <8 следует p(U(Cij),U(c ц)) <е, где Cj -точные коэффициенты полинома (6); Cij - приближенные коэффициенты; р - обобщенная метрика, причем е(8)^0 при 8^0, что говорит об устойчивости рядов Фурье. Значит, существует некоторая 8-окрестность, в которой изменения коэффициентов Фурье не приводят к существенному изменению самого ряда.

Устойчивость коэффициентов Фурье допускает некоторый произвол в их измерении, что особенно важно для практических задач, так как вектор x может быть измерен только приближенно. Для определения интерполируемых значений коэффициентов Фурье (6) введем следующее определение.

Определение. Функция y(t) е LL

[a,b]

называ-

ется стробирующей [2]; если Ж(хк) = 1ек е [а,Ь], то у(1;к) > е [а,Ь],1ек#1. Нетрудно усмотреть, что стробирующая функция удовлетворяет (4), следовательно,

v L

SS Cij(x)9ij(t) = ¥(t),

i=ij=i

отсюда в силу ортогональности b

J V(t)9ij(t)dt

C = a_

Ja" (t)dt

(7)

(8)

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

Отыщем приближенное решение (2).

1. Предположим, что Cjj=Cjj(xk), i=1,2,...,v, j,k=1,2,...,n, xkeQx. В этом случае можно выделить два алгоритма.

а) Признак с номером k относится к одному из коэффициентов Cjj каждого класса. Тогда, очевидно, L=n и алгоритм упорядочения признаков x/k сводится к отысканию того признака, для которого

p I I

S Cjj(x,k)-Pjj(xJ ^ min , (9)

i=1,2,.,v, j,k=1,2,.,n, x/kE Qx.

Pjj(xlk) - полином наилучшего приближения коэффициентов Фурье стробирующей функции ■y(t) с аргументами, расположенными в k-м столбце матрицы А.

В каждом классе (9) отыскание признака проводится для всех n-1 коэффициентов Фурье с отбрасыванием из рассмотрения тех столбцов, которые ранее были отнесены к какому-либо коэффициенту. Наконец, к коэффициенту с номером n относится свободный столбец.

б) Произвольный из n признаков относится к любому из n коэффициентов для каждого из v классов при выполнении (9). Может случиться так, что часть столбцов матрицы исключается из рассмотрения.

2. Для Cjj строятся полиномы наилучшего приближения, зависящие от двух и более признаков. Упорядочение признаков при этом аналогично (9).

3. Положим

Cjj = Fc(x). (10)

В этом случае Cjj определяется, как в (4). Коэффициенты Фурье в аналогичном преобразова-

SS C*j(x)9ij(t) = U*(x,t)

i=ij=i

(11)

определяются по алгоритмам 1а, 1б.

Для С* 5-окрестность будет несколько больше, так как определяемые через них коэффициенты Су допускают некоторую неточность в измерении.

Аналогичное (10) преобразование можно, в свою очередь, применить для С* (х). Обозначим

оператор q-квадратного вложения (в смысле (10)) ортогонального полинома в аппроксимируемую функцию через Т41. Тогда для 1а, 1б можно записать

F(x) = Ti(x,A)

и для алгоритма 3 F(x) = T2(x,A),

(12)

(13)

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

Рассмотрим тригонометрическую систему функций {sin kx, cos kx}, ортогональную на [-п,п], в качестве нелинейной разделяющей функции. При этом коэффициенты Фурье рассматриваются как факторы решаемой задачи. В этом случае T(x) = t, te [-п,п],

2 ak (х) cos kt + bk sinkt = U(x,t).

k=1

Пусть

(14)

in, te [t -e, t +e], ¥(t)= \ - (15)

[0, te[ti -e, ti + e],t e [-n,n], e> 0.

Интерполируемые значения коэффициентов Фурье определяются в виде

. 1

ak = — [sink(ti +e) -sink(ti -e)], k

. 1

bk = — [cos k(t, - e) - cos k(t, + e)]. k

(16)

(17)

Возьмем в качестве Су(х) линейную комбинацию признаков х1, х2,..., хп, то есть

п

С„(х) = Е qkij + (18)

' к=1 ' '

Коэффициенты линейной формы представления коэффициентов тригонометрического ряда вычисляются с помощью метода наименьших квадратов. Коэффициенты qkij и qij можно опреде-

лить из среднеквадратичного приближения значений Су^^) для соответствующих векторов хi обучающей матрицы А. Если х соответствует ш и

C„(s,q)-

z q... + q..

k=1 kij U

< 8„

(19)

классификация будет устойчивой.

На основании изложенного сделаем следующие выводы.

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

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

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

Литература

1. Штейнбух К. Некоторые аспекты распознавания при моделировании и решении задач // Распознавание образов. Адаптивные системы: тр. Междунар. симпоз. по технич. и био-логич. проблемам управления; гл. ред. В.А. Трапезников. (24— 28 сентября 1968 г., Ереван). М.: Наука, 1971. 270 с.

2. Журавлев Ю.И., Рязанов В.В., Сенько О.В. Распознавание. Математические методы. Программная система. Практические применения. М.: ФАЗИС, 2006. 176 с.

n

АВТОМАТИЗИРОВАННОЕ ПРОЕКТИРОВАНИЕ ПРИРОДНО-ТЕХНИЧЕСКИХ СИСТЕМ УТИЛИЗАЦИИ ОТХОДОВ

С.Н. Костарев, к.т.н. (Пермский государственный технический университет, iums@dom.raid.ru)

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

Ключевые слова: АРМ твердых бытовых отходов, САПР, полигон твердых бытовых отходов, фильтрат.

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

как природно-техническая система утилизации отходов (система «Природа-Техника-Отходы» (ПТО)). Под системой ПТО понимается совокупность действующих на определенной территории инженерных сооружений (технических средств) в пределах участка захоронения отходов (УЗО) и хозяйственной зоны полигона ТБО, транспортных средств, системы коммуникаций, предназначенных для складирования, изоляции и обезвреживания ТБО, а также сфер их проектирования, строительства, реконструкции и содержания.

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