13. Anchurin S.A., Maksimov V.N., Morozov E.S., Golovan' A.S., Shilov V.F. Blok inertsi-al'nykh datchikov [Unit inertial sensors], Nano- i mikrosistemnaya tekhnika [Nano- and Microsystem technology], 2011, No. 1, pp. 50-53.
14. Elwenspoek M., Wiegerink R. Silicon micro accelerometers, Mechanical microsensors, 2005, pp. 230-236.
15. Lysenko I.E. Modeling of the micromachined angular rate and linear acceleration sensors LL-type with redirect of drive and sense axis, World Applied Sciences Journal, 2013, No. 27 (6), pp. 759-762.
16. Lysenko I.E. Integral'nyy sensor uglovykh skorostey i lineynykh uskoreniy [Integrated sensor of angular velocities and linear accelerations], Inzhenernyy vestnik Dona [Engineering journal of Don], 2010, No. 3. Available at: http://ivdon.ru/magazine/ (free access).
17. Lysenko I.E., Lysenko A. V. Integral'nye sensory uglovykh skorostey i lineynykh uskoreniy LR-tipa na osnove uglerodnykh nanotrubok // Inzhenernyy vestnik Dona [Engineering journal of Don], 2012, No. 4. Available at: http://ivdon.ru/magazine/ (free access).
18. Lysenko I.E., Ezhova O.A. Kriterii ravenstva sobstvennykh chastot kolebaniy chuvstvitel'nykh elementov mikromekhanicheskikh giroskopov-akselerometrov [The criteria of equality of the natural frequencies of oscillation of the sensitive elements of micromechanical gyroscopes-accelerometers], Inzhenernyy vestnik Dona [Engineering journal of Don], 2014, No. 2. Available at: http://ivdon.ru/ru/magazine/archive/n2y2014/2475 (free access).
19. Lysenko I.E. Funktsional'no integrirovannye mikro- i nanomekhanicheskie sensory uglovykh skorostey i lineynykh uskoreniy [Functionally integrated micro - and nanomechanical sensors of angular velocities and linear accelerations]. Taganrog: Izd-vo YuFU, 2013, 167 p.
20. Lysenko I.E., Kulikova I.V., Polishchuk E.V., Khayrulina V.A. Modelirovanie elementov mikrosistemnoy tekhniki v programme ANSYS [Modeling elements of Microsystems technology in the ANSYS program]. Taganrog: Izd-vo TTI YuFU, 2007, 42 p.
21. Abramov I.I. Problemy i printsipy fiziki i modelirovaniya pribornykh struktur mikro- i nanoelektroniki. Ch. II. Modeli poluklassicheskogo podkhoda [Problems and principles of physics and modeling of device structures of micro- and nanoelectronics. Part II. Model semiclassical approach], Nano- i mikrosistemnaya tekhnika [Nano- and Microsystem technology], 2006, No. 9, pp. 26-36.
Статью рекомендовал к опубликованию д.т.н., профессор Б.Г. Коноплев.
Лысенко Игорь Евгеньевич - Южный федеральный университет; e-mail: [email protected]; 347928, г. Таганрог, пер. Некрасовский, 44; тел.: 88634371603; кафедра конструирования электронных средств; д.т.н.; зав. кафедрой.
Ежова Ольга Александровна - e-mail: [email protected]; кафедра конструирования электронных средств; аспирант.
Lysenko Igor Evgenievich - Southern Federal University, e-mail: [email protected]; 44, Nekrasovsky, Taganrog, 347928, Russia; phone: +78634371603; the department of electronic apparatuses design; professor.
Ezhova Olga Aleksandrovna - email: [email protected]; the department of electronic apparatuses design; postgraduate student.
УДК 519.711 Б01 10.23683/2311-3103-2017-4-232-245
С.В. Фролов, Т.А. Фролова
МЕТОДЫ РЕШЕНИЯ ГЛОБАЛЬНОЙ ЗАДАЧИ УПРАВЛЕНИЯ ДЛЯ РАСПРЕДЕЛЁННЫХ ИЕРАРХИЧЕСКИХ СИСТЕМ С ИНТЕРВАЛЬНОЙ НЕОПРЕДЕЛЁННОСТЬЮ
В большинстве практических случаев для неопределённых параметров технологических процессов нахождение характеристик законов распределения и функций принадлежности является затруднительным. Обычно для неопределённых параметров известны только интервалы их значений. Целью исследования является разработка теоретических
основ построения системы управления сложным химико-технологическим процессом в условиях интервальной неопределённости параметров. Сложный химико-технологический процесс совместно с системой управления представляет собой сложную химико-технологическую систему. Первый или нижний уровень управления иерархической структуры составляют аппараты или машины совместно с локальными системами управления, построенными на основе программируемых микроконтроллеров. Элементы нижнего уровня иерархии объединяются по технологическому принципу в блоки, которые вместе с системой управления блоком образуют второй уровень иерархии. Третий или верхний уровень иерархии составляет совокупность элементов второго уровня и система управления этой совокупностью. Управление подпроцессами второго уровня иерархии осуществляют координирующие станции, которые получают координирующие сигналы от системы управления третьего уровня, информационные сигналы с нижнего уровня и подают управляющие сигналы на нижний уровень иерархии. Координирующая станция верхнего управляет системами второго уровня. Поставлены глобальные задачи управления распределённой трёхуровневой иерархической системой с интервальными параметрами: задача получения удовлетворительных решений и задача оптимизации. Задача получения удовлетворительных решений сводится к поиску управлений, при которых для любых значений неопределённых параметров из заданных интервалов выполняются установленные технологические требования. Глобальная задача оптимизации дополняется условием достижения критерием оптимального значения в условиях неопределённости. Сформулированы принципы, предложены методы и алгоритмы декомпозиции глобальной задачи управления, позволяющие для технологических процессов с интервальной неопределенностью параметров решение глобальной задачи свести к параллельно решаемым задачам управления. Изложенный подход дает возможность на стадии разработки для каждой задачи применить принцип модульности. Этот принцип заключается в том, что при формализации, постановке и алгоритмизации каждой задачи может участвовать свой независимый разработчик, которому не обязательно знакомиться с содержанием других задач. Принцип модульности позволяет изменять или модифицировать задачу, не вызывая каких-либо последствий в других задачах. Результаты работы являются теоретической основой для проектирования новых и модернизации существующих систем управления производствами, характеризующихся множеством стадий, энерго- и материалоёмкостью, наличием большого количества параметров и возмущений.
Химико-технологический процесс; интервальная неопределённость; иерархическая система управления; декомпозиция задачи управления; задача оптимизации.
S.V. Frolov, T.A. Frolova
METHODS FOR SOLVING THE GLOBAL CONTROL PROBLEMS FOR DISTRIBUTED HIERARCHICAL SYSTEMS WITH THE INTERVAL
UNCERTAINTY
Finding the characteristics of distribution and membership functions is difficult in most practical cases for uncertain parameters of technological processes. Usually only the intervals of uncertain parameters values are known. A complex chemical-technological process together with a control system is considered as a complex chemical-technological system. The first or the lower level of control of the hierarchical structure is composed of machines or machines in conjunction with local control systems based on programmable microcontrollers. Elements of the lower level of the hierarchy are coupled according to the technological principle into blocks, which together with the block management system form the second level of the hierarchy. The third or upper level of the hierarchy is the set of the elements of the second level and the control system. Control of the second-level sub-processes of the hierarchy is carried out by coordinating stations that receive coordinating signals from the third-level control system, information signals from the lower level, and provide control signals to the lower level of the hierarchy. The coordinating station of the top controls the systems of the second level. The aim of the study is the formulation of theoretical foundations for the development of the control system for complicated chemical processes with interval parameter uncertainties. Global control problems are formulated for a three-level hierarchical system with interval uncertainty of parameters: the problem of obtaining satisfactory solu-
tions and the optimization problem. The problem of obtaining satisfactory solutions reduces to the search control in which for all values of uncertain parameters from the specified intervals the technological requirements are fulfilled. The global optimization problem is supplemented for achieving optimal value criterion under uncertainty. The principles are formulated; methods and algorithms are proposed for decomposition of the global control problems, which allow reducing a global control problem to a set of control problems, which could be solved in parallel for technological processes with the interval uncertainty of parameters. The presented approach makes it possible to apply the modularity principle for each task at the development stage. This principle is that each of the tasks can be formalized, set up and algorithmized by an independent developer who does not need to get acquainted with the content of other tasks. The modularity principle allows modifying a task without causing any consequences in other tasks. The results are the theoretical basis for the design of new and modernization of existing control systems for production, which is characterized by a plurality of stages, energy and materials consumption, the presence of a large number of parameters and disturbances.
Chemical-technological process; interval uncertainty; hierarchical control system; decomposition of control problem; optimization problem.
Введение. При управлении современными технологическими процессами существует проблема высокой размерности и неопределенности их параметров [1]. Проблему размерности решает иерархическое построение систем управления, декомпозиция задачи управления и параллельная организация вычислительного процесса [2, 3].
Другая характеристика - неопределенность неизбежно присутствует практически в любой задаче проектирования и управления (принятия решения) [4, 5]. Под неопределенностью будем понимать неоднозначность исходных данных, относящихся к технологическому процессу, на основе которых вырабатываются управляющие воздействия.
Неопределенность может вызываться стохастическим поведением объекта (действие неконтролируемых внутренних и внешних возмущений), погрешностями датчиков, исполнительных механизмов, информационных каналов, неточностью математической модели объекта и т.п.
Для работы в условиях неопределенности известны три основные методологии: теория нечетких множеств [6], стохастический подход [5] и интервальный анализ [4, 5, 7].
В большинстве практических случаев для параметров технологических процессов определение характеристик законов распределения и функций принадлежности является затруднительным. Обычно для неопределённых параметров известны только интервалы их значений [4, 5, 7].
Актуальной задачей является разработка новых методов управления сложными технологическими процессами с интервальной неопределённостью параметров на основе распределённых систем, построенных по иерархическому принципу.
1. Принципы построения и функционирования иерархических систем управления технологическими процессами. В терминологии общей теории систем [8, 9] технологические процессы, которые обеспечивают превращение исходного сырья в целевой продукт путем прохождения материальных потоков через определенное число функционально различных ступеней переработки, относятся к классу сложных систем. На основе общепринятой терминологии, эти системы называются химико-технологическими системами (ХТС) [10, 11].
Элементом ХТС считается аппарат или машина, в которых протекает химико-технологический процесс, качественно (или) и количественно преобразующий физические переменные входных материальных и энергетических потоков в физические переменные выходных материальных и энергетических потоков. Таким образом, каждый элемент ХТС осуществляет преобразование, которое представляется типовым технологическим оператором (ТТО).
Процесс, составляющий отдельную единицу ХТС, в общем случае формализуется как физико-химическая система (ФХС) [11]. Оператор ТТО является суперпозицией целого ряда элементарных технологических операторов, и имеет смешанную детерминированно-стохастическую природу.
В соответствии с [10] ХТС совместно с системой управления (АСУ ТП) образует сложную ХТС (СХТС). Эта СХТС представляет собой трехуровневую иерархическую систему (рис. 1). Под иерархической системой будем понимать систему, состоящую из взаимосвязанных подсистем, каждая из которых в свою очередь является иерархической по своей структуре.
СХТС ОпСо КСс
—I— Общая вычислительная сеть — ус > С "'
I МХТС! ОПС1 1 КС1 МХТС, Опй 1 КС; 1
Локальная вычислительная сеть Локальная вычислительная сеть
| |лтс| |лтс [ггтс]... [Йт 3 Гр| |лтс| |ПТС| |ПТС лтс| _ |лтс |птс| _ _ |лтс
Тп ти ... Тш, Т.-1 ...
СФХС,] СФХСц СФХСШ1 СФХС,! СФХС9.
Рис. 1. Иерархическая трехуровневая структура СХТС
Первый или нижний уровень иерархической структуры СХТС составляют ТТО и локальные системы управления. Средством нижнего уровня управления ТТО является локальная технологическая станция (ЛТС) [12-17], которая непосредственно соединяется с датчиками, исполнительными устройствами и решает автономные задачи управления ТТО. В современной терминологии ЛТС представляет собой программируемый контроллер или микропроцессорный регулятор [18],
снабженный средствами для работы в составе сети. Совокупность ТТО (Т) (рис. 1)
и локальной системы управления, включающей ряд ЛТС, будем называть сложной ФХС (СФХС).
Сложная ХТС (СХТС) разделяется на отдельные блоки, выделенные по технологическому принципу. Блок I (I = 1,N) состоит из совокупности
(СФХС,СФХС,2,..,СФХС,,,..СФХС,М), образующей подпроцесс Р1, и системы управления подпроцессом Р^. Каждый блок называется малой ХТС (МХТС) [11], а множество (МХТС, МХТС2, МХТС„ МХТС^) образует второй уровень иерархии СХТС.
В состав системы управления второго уровня иерархии входит операторская станция (ОпС,) и координирующая станция (КС,) (1 = 1,N [12-17]. Операторская станция выполнена в виде промышленного компьютера и выполняет функции диспетчерского управления, визуализации процесса, сбора и архивации данных, обеспечивает связь между оператором и ТТО, а также служит задачам по программированию, настройке и контролю за работой ЛТС. Координирующая стан-
ция предназначается для решения сложных прикладных вычислительных задач и является промышленным компьютером большой вычислительной мощности. Связь между первым и вторым уровнем иерархии в рамках МХТС осуществляется через локальную вычислительную сеть (ЛВС), к которой подключены ЛТС, ОпС,- и КС, (7 = 1,Ы).
Третий или верхний уровень иерархии СХТС составляет совокупность (ЫХТС, МХТС, - - -, МХТС,,..., МХТС^) и система управления этой совокупностью. В состав системы управления верхнего уровня СХТС входят связанные по общей вычислительной сети ОпС0 и КС0. Станции верхнего уровня выполняют функции по управлению МХТС, которые аналогичны задачам, решаемым ОпС,- и КС,- на втором уровне иерархии. Связь между вторым и третьим уровнями осуществляется через устройства связи (УС).
Таким образом, СХТС соответствует определенный технологический процесс (ТП), обеспечивающий превращение исходного сырья в готовый для реализации конечный продукт, т.е. исходные и конечные продукты СХТС являются товарными.
Пусть ,-ый подпроцесс, которому соответствует совокупность (СФХС,!, СФХС,.2,..., СФХС ,..СФХС,м) задается отображением:
Р : и х X х В х Е ^ ^ 7 = 1, Ы, где и - множество управляющих сигналов и,-,
приходящих с КС,; Хг - множество контролируемых возмущений х,-; Вг - множество входов Ь,-, посредством которых подпроцесс Р, связывается с другими подпроцессами; - множество неопределенностей, связанных со всеми факторами f ^ е Е), влияющими на выход у,- подпроцесса Р,-; ^ - множество выходов у,- подпроцесса Р,.
Будем считать, что множество неопределенных параметров определяется интервальным вектором, т.е. [У ] = Е.
Будем рассматривать ,-ую координирующую станцию КС, (рис. 2) как локальную систему управления подпроцессом Pi.
Рис. 2. Иерархия процесса решения глобальной задачи СХТС
К системе КС, поступают сигналы двух видов: координирующий сигнал g,- (& еО,) от вышестоящей системы управления КС0 и информационный сигнал обратной связи z,- ^ е2) от подпроцесса Р,-. Множества О, и образуют соответственно множества координирующих и сигналов обратной связи. Выходом КС, является локальное управление и,-. выбираемое из множества управлений и, (и,-е и, ). Считается, что КС, реализует отображение:
КС,: О, х 2, ^ и,
или
ц = КС, ^
7 = .
Множества управляющих сигналов и, информационных сигналов 2 и выходов У в масштабах всей СХТС представляются в виде декартова произведения:
и = и хиг х.хи х...хи„,
2 = 2 х 2 2 х . х 2,. х . х 2Ы, 7 = у хУ2 х...ху х...хУ„. Составляющие g,- выходного сигнала g станции КС0 являются координирующими для станций КСЬ..., КС,-,..., КС№ где
Я = ^ g2,—, & Я ы) , О = О1 х О2 х...х х.х °, § е О Системы управления можно подразделить на жесткие и адаптивные. Жесткой системой будем называть систему, имеющую неизменяемый оператор, который устанавливает связь между векторами входных параметров и выходным вектором управлений. Адаптивной системой будем считать систему, в которой оператор, связывающий вход и выход, меняется в зависимости от входных параметров. В рассматриваемом случае имеется отображение вида:
К : О х 2 ^ А,
или
КС, = К, (§,, г,),
I = 1, N,
где А - множество операторов управления. Для жесткой системы множество
операторов управления состоит из одного элемента: А, = {КС,}.
Станция КС0 осуществляет отображение: КС: Щ ^ О, в котором
Щ = № х № х . х Щ х ... х , где Щ - множество поступающих со станции
КС, информационных сигналов, с помощью которых осуществляется обратная связь. Сигнал обратной связи z,- , поступающий на вход локальной системы управления КС,-, содержит информацию относительно поведения подпроцесса Р, и связан функциональной зависимостью с управляющим сигналом и,-, интервальным вектором неопределенных параметров и выходом у,:
р: иi х Е х У ^ 2.
, = 1Ы.
Поступающий по каналам обратной связи информационный сигнал М направляется в вышестоящую систему управления КС0, содержит в себе информацию относительно поведения нижестоящих систем управления и задается отображением: р: О х 2 х и ^ Щ, где Щ - множество информационных сигналов
М (М е Щ). Информационный сигнал М ^^,М,.••,Ч,.••,М), (М еЩ) является функцией координирующего сигнала Я, информационного сигнала обратной связи г = (г,., г,. • •, , (г е 2) от нижестоящих систем и управляющего воздействия и = (ц,., и,,., и^, (и е и), поступающего на нижестоящие системы.
Предполагается, что существует связанная с системой управления КС0 решаемая задача D0( М) и множество решений О, таких, что для любого w е W выход g = КС0( М ) ^ е О) есть решение задачи D0( М). Таким образом, сигнал обратной связи М конкретизирует задачу D0( М). Аналогично координирующий сигнал g и сигнал обратной связи z конкретизирует совокупность задач
{О^zl),.--А^,г,),^,гN)} Задача О,(е,,г,) связана с нижестоящей системой управления КС,- так, что для каждой пары zi) из множества ^ х ^ вектор управлений и,-, определяемый как и = ^ (&■, г1), является решением задачи О (^, г).
Для упрощения дальнейших выкладок определим Р(и,й), где D - произвольно решаемая задача (совокупность задач), как предикат:
Р(и, й) = и есть решение задачи (совокупности задач) D.
Пусть с СХТС связана глобальная задача й. Задачи {й1 (^, г1 ), .--А (^, г ), .--,0^ (^, )}, которые решаются нижестоящими системами управления КС, = 1, N), координируемы относительно глобальной задачи
й , если справедливо следующее предложение:
(3е),(3и):[р(и,{й1(е1, г1 ),...Д ^, г, (gN, ZN)}) и р(и,й)]. (1)
Координируемость относительно глобальной задачи й означает, что система управления (координатор) КС0 может влиять с помощью сигналов
е = (е, е,gi,на нижестоящие системы управления КС, (/ = N) так, чтобы их результирующее воздействие и = (и, и,Ц,и^) на процесс в целом давало решение глобальной задачи й .
Координирующие сигналы g из О могут содержать прогноз связующих входов подпроцессов Р,, которые будут иметь место в связи с подачей управляющих воздействий и =(ц, и2,..., Ц ,..., и^). Прогноз связующих входов представляется
в виде оценочной области Б8 = В8 х В8 х ... х В8 х ... х В8, которой принадлежит вектор связующих входов Ь8 е В8, (ь8 = (ь8,Ь8,...,Ь8,...,Ь^). Оценочную область в8 будем представлять в виде интервального вектора, т.е. [ь]8 = в8
и [Ь]8 = В8,1 =
Связующие сигналы Ь = (ь,Ь2,.••,Ь,,.,Ь^), поступающие на входы подпроцессов Р,, могут быть функционально связаны с управляющими воздействиями и, контролируемыми возмущениями х и неопределенными параметрами Г : :их Xх Г ^В, где множество связующих входов В является подмножеством множества выходов У (В с У).
Для решения проблемы координации сформулируем принцип оценки взаимодействий. Форма принципа оценки взаимодействий выражается предложением: №(Уи): {Ии^.г),...Д(е,^^)}) ири,хДОе В8(2)
^р(и,О )}
Координатор КС0 определяет область Bg , в пределах которой составляющие вектора связующих входов bg могут варьироваться. Этот принцип утверждает,
что управляющее воздействие u решает поставленную глобальную задачу D всякий раз, когда u является решением задач {Dj (g, z), •••, D (g, Z),•••, D^ (gw, zN)} и фактические значения составляющих вектора связующих входов [b] = R(u, x, [f ]), появляющиеся при применении управляющего воздействия u, попадают внутрь
прогнозируемой оценочной области
2. Постановка глобальной задачи управления СХТС. Глобальная задача
СХТС D может быть задачей получения удовлетворительных решений (задачей "удовлетворения") [19, 20] или задачей оптимизации.
Глобальная задача "удовлетворения" формулируется следующим образом. Пусть технологические требования ф (s = 1, S) задают отображения:
фх: U х X х Y ^ R, 5 = 1S
Технологическому процессу СХТС соответствует отображение: P: U х F ^ Y, где F - множество всех возможных значений неопределенных параметров; P - оператор технологического процесса.
Дано множество допустимых управлений U. Требуется найти такое u е U для заданного x е X, при котором справедливо:
^ ф;(U x as, (3)
[у]= P(u,x,[f ]), (4)
5 = 1, S,
где [y] = ([y ],[[ ],„., [у. ],„., [y^ ] - интервальный вектор выходов подпроцессов P;
[f ] - интервальный вектор неопределенных параметров; параметры as, (5 = 1, S)
определяют уровни технологических требований.
Решение u считается удовлетворительным, если оно приводит к значениям
(Vs): min ф(и, x, у), не опускающихся ниже определенного уровня as.
уе[у]
Для глобальной задачи оптимизации СХТС вводится понятие целевой функции:
J: U х X х Y ^ Q ,
или (5)
J = J(u, x, у),
где Q - множество значений критерия целевой функции J.
Глобальная задача оптимизации СХТС формулируется следующим образом. Необходимо найти такие и* е U, у* е Y , при которых для заданного x е X,
любого и е U, у е Y справедливо: j(u*,x,y*)> max J(u,x,y), и выполняются
уф]
условия (3)-(4).
В частном случае для задач "удовлетворения" и оптимизации множество неопределенностей F может быть пустым множеством: F = 0.
3. Решение проблемы декомпозиции глобальной задачи СХТС. Для решения глобальной задачи D необходимо решить проблему ее декомпозиции. Под декомпозицией задачи D понимается определение таких задач D0( W) и
{01 ,г1), ..^Р,. ,г, \...,РЫ,гм)}, которые могут быть поставлены соответственно перед вышестоящей КС0 и нижестоящими системами управления КС,-| ■ = 1, N | так, чтобы имело место выполнение следующих условий.
Должен быть выполнен принцип координации.
Должен быть найден эффективный алгоритм поиска координирующего сигнала g такой, чтобы задачи {Р1 ^,z1),..., Р,,,zг■),..., ОNN,zN)} были координируемы относительно глобальной задачи Р, т.е. должно быть справедливо предложение (1).
Иерархия решения глобальной задачи СХТС Р показана на рис. 2. Координатор КС0 получает информацию Щ от локальных блоков КС,- , путем
решения задачи й0( Щ) выбирает координирующий сигнал g(k) и посылает нижестоящим системам управления КС,- ( ■ = 1, N) сигнал
g(k) = g(k),...,^. Получив координирующий сигнал g(k) и сигнал
обратной связи z,-, локальная система управления КС,- определяет управление и ^ ), являющееся решением задачи Р.г.), и посылает информацию координатору KC0 в виде сигнала . На основе полученной информации щ= ^,ЩЩ,от всех КС/ -■ = 1,N1 и выбранной стратегии координации: Б0: О х W ^ О координатор КС0 посылает координирующий сигнал g(k+l) = (g(k+1),g!,k+1),...,g(k+1),...,gjk+1)) системам управления КС/ - = 1,ж). Итерационный процесс g(k+1) = Б0 (ц,(к), ))) продолжается до тех пор пока не будет
достигнуто решение глобальной задачи Р .
Таким образом, решение глобальной задачи СХТС Р представляет собой итерационный процесс, в котором вычислительные операции распределяются между всеми системами управления КС0, КСЬ..., КС/,..., КСдт .
Для того, чтобы определить задачи D0(w),
{РД&,г1),...,Р,(g¡,гi),...,Рж^,гж)} проводятся следующие рассуждения. Будем считать, что [Ь; ) есть формируемый КС0 на к-ой итерации интервальный вектор связующих входов --ой МХТС, т.е.
[Ь,Г) - ). (6)
Вектором неопределенных параметров --го подпроцесса можно считать интервальный вектор [^ ], который определяется как декартово произведение:
[у, ]- ^) х ^
или
[у,] = [Ь,Г х[г] (7)
Вектор выходных параметров [у. ] выражается посредством оператора математической модели подпроцесса Р.
[у, ] = М,(и,, х,, [у,]), (8)
Критерий оптимизации (5) можно представить в виде аддитивной функции:
N
J = IJ • (9)
Суть предлагаемого метода решения глобальной задачи оптимизации СХТС
D заключается в фиксации оценочной области B, в пределах которых может варьироваться вектор связующих входов b. Оптимизация на каждой к-ой итерации (рис. 3) сводится к оптимизации отдельных МХТС:
max J<kk = max min J{k> +... + max min J(k> +... + max min J<k>.
u^eU ufeUi yie[yi ] ufeU, y,e[y, ] uJe^n У n e[y n ]
Оптимальное значение критерия J* зависит согласно (6)-(8) от интервального вектора связующих входов [b]: J* = J*([bj), где [b]=([b1], [b2 ],..., [b, ],..., [bN ]).
В этом случае оптимальное значение критерия (9) находится в результате поиска в допустимой области B такого интервального вектора связующих входов
[b ]*, при котором:
4b]*)=max Jb]).
С точки зрения иерархии рассматриваемого процесса управления (рис. 2) роль сигналов обратной связи W, играют величины J *, роль сигналов g,- выполняет интервальный вектор связующих входов [b]. Таким образом, глобальная задача
оптимизации СХТС D разделяется на задачи D0( W), {DiCgi, zi), .,D, (g,, z,),..., D
NVfeN'
Z N )}.
Для формулирования задач оптимизации и "удовлетворения" вводятся следующие обозначения. Пусть l-ый компонент вектора входных переменных b,-
МХТС,- обозначен как Ъг.(/ = 1, N; l = 1, Li), а r-ый компонент вектора выходных
переменных у/ МХТС/ как у . (j = 1, Н; г = 1, ^ ), где Ь, ^ - размерности вектора выходных переменных соответственно МХТС,-. и МХТС/. Тогда элемент матрицы связи технологических потоков СХТС записывается как 2^. Этот эле-
¡Г 1
мент принимает значения: ^ = 1, если г-ый компонент вектора у/ МХТС/ является 1-ым компонентом вектора Ь,- МХТС,, т.е. верно равенство = у ■. Элемент
¡Г I
2.. имеет нулевое значение: 2 = 0, когда х,, Ф у,,-.
.} . ¡. и
Для того, чтобы значения составляющих вектора выходных параметров попадали в оценочную область к) должны выполняться условия:
£(£Ы- о, сю)
.У 1Ч Srj^, IJ
i=i V I=1 I=1
(r = 1, rj, j = 1, N),
(r = 1, Rj, J = 1, N),
где Ьц, Ьц - верхние и нижние границы интервального вектора [Ь, ] связующих
входов (Ь, ]е Б(к), г =
В этом случае неравенства (10), (11) включаются в ограничения для задачи Р; & 1, г; ), (/ = 1, N) в виде:
(и 1,х 1,У1 ^ ая,,
Здесь
q j = \,2RJ.
N L, L, N L, , _ N
^.x, )=-e e< ^=-e e j, ц
1=1 1=1 1=1 1=1 1=1
N L, L, N L,
v(qJ+1)(u j, x, y j)=EE у* E=EE№ J
1=1 ;=1 ;=1 1=1 ;=1
(' = 1, Rj),
где bii, bn (l = 1, Ц) - верхние и нижние границы интервального вектора
[ь, i i=in.
Задача оптимизации D (g, z. ), которая ставится перед системой управления KC,, формулируется следующим образом. Необходимо для заданного x. е Xi, вектора связующих входов [ь] найти управление u * е U., при котором:
J*(u*Д,)= max min J(u,д,), при выполнении условий:
y"ein, Mu' 'x, ,y)-(12)
q = 1,2,...,2^,., min w (u,,x,,y)> а ,
г l т q, \ i? i / q, '
y,e[y,] (13)
q, = 2Д. +1,..., Qi,
[y,.] = M,(u, , x, ,[v, ]), (14)
[vi]=[bi]x[f], (15)
[b]=([bil [b2 ],., [ь, ],., [ь„ ]).
Задача D0( w), которая решается системой управления KC0, формулируется в следующем виде. Необходимо найти такой вектор [ь]* и соответственно такие
значения •*([ь]*), •2*(ь]*),..,•*([ь]*),...,•*(ь]), при которых достигает максимума
критерий (9). Поиск оптимального интервального вектора [ь]* сводится к определению его верхних и нижних границ, которые представляются в виде ь =(ь*,ь2,...,ь* ,...,ьn), ь* =(ь*,ь*,.,ь*,.,ь*¥),те. задача оптимизации, решаемая KC0, записывается в виде:
j* ь ь )=f^B ej ([ь]).
При решении глобальной задачи "удовлетворения" D задача D (g, Z) сводится к отысканию такого вектора u„ при котором для заданного x„ выполняются условия (12)-(15). Задача D0( W), в этом случае заключается в генерации вектора
связующих входов [ь], (верхних и нижних границ: Ь, Ь). В качестве сигналов обратной связи ^ служит информация о нарушении ограничений (12),(13), сигналов обратной связи g,- является вектор [ь].
Предложенный метод решает проблему декомпозиции глобальной задачи
СХТС Р. Во-первых, выполняется принцип координации: принцип оценки взаимодействий (2). Во-вторых, задача 00( М ), относящаяся к системе управления КСо, обеспечивает поиск координирующего сигнала - интервального вектора связующих входов [ь]. Например, при решении глобальной задачи оптимизации при поиске оптимального вектора связующих входов [ь] могут быть использованы известные методы нелинейного программирования [21], которые позволяют определить оптимальные верхние и нижние границы вектора [ь].
Таким образом, предложенный метод позволяет решение глобальной задачи СХТС Р свести к параллельно решаемым на станциях КС0, КС,- (/ = 1, Н) задачам Оо( М X {р1(&, (gí, z¡), ...,й„ , z N )}.
Изложенный подход дает возможность на стадии разработки для каждой задачи р (^,z ) (/ = 1,N) применить принцип модульности [22]. Этот принцип заключается в том, что при формализации, постановке и алгоритмизации каждой задачи р (^, z ) может участвовать свой независимый разработчик, которому не
обязательно знакомиться с содержанием других задач р (^, z ) (у ф г). Принцип модульности позволяет изменять или модифицировать задачу р.,), не вызывая каких-либо последствий в других задачах р (^, z ) (у ф г).
Заключение. Сформулированные в работе принципы и постановки задач, предложенные методы и алгоритмы декомпозиции глобальной задачи управления для технологических процессов с интервальной неопределенностью являются теоретической основой для проектирования новых и модернизации существующих производств, характеризующихся множеством стадий, энерго- и материалоёмкостью, наличием большого количества параметров и возмущений.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Островский Г.М., Зиятдинов Н.Н., Лаптева Т.В. Оптимизация технических систем.
- М.: КНОРУС, 2012. - 432 с.
2. Фролов С.В. Развитие отечественных систем промышленной автоматизации // Автоматизация и производство. - 2001. - № 2. - С. 22-25.
3. Шильяк Д. Децентрализованное управление сложными системами. - М.: Мир, 1994.
- 576 с.
4. Холоднов В.А., Лебедева М.Ю. Приближенный метод оптимизации химико-технологических систем при интервальной параметрической неопределенности // Известия Санкт-Петербургского государственного технологического института (технического университета). - 2016. - № 33 (59). - С. 97-100.
5. Островский Г.М., Лаптева Т.В., Зиятдинов Н.Н. Проектирование оптимальных химико-технологических систем в условиях неопределенности // Теоретические основы химической технологии. - 2014. - Т. 48, № 5. - С. 527-537.
6. Пегат А. Нечеткое моделирование и управление. - М.: БИНОМ. Лаборатория знаний, 2013. - 798 с.
7. Фролова Т.А., Фролов С.В., Туляков Д.С. Решение интервальных математических моделей технологических процессов // Электронное научно-техническое издание «Наука и образование». - 2012. - № 09. - С. 343-360. - UPR: http://technomag.bmstu.ru/doc/454499.html (дата обращения: 17.02.2017).
8. Растригин Л.А. Современные принципы управления сложными объектами. - М.: Сов. радио, 1980. - 232 с.
9. Саймон Г. Науки об искусственном: пер с англ. - 2-е изд. - М.: Едиториал УРСС, 2004.
- 144 с.
10. Кафаров В.В., Мешалкин В.П. Анализ и синтез химико-технологических систем. - М.: Химия, 1991. - 432 с.
11. Кафаров В.В., Дорохов И.Н. Системный анализ процессов химической технологии. Основы стратегии. - М.: Наука, 1976. - 499 с.
12. Харазов В.Г. Интегрированные системы управления технологическими процессами: учебное пособие для вузов. - СПб.: Профессия, 2013. - 655 с.
13. Харазов В.Г. Интеллектуальные приборы и системы управления // Известия Санкт-Петербургского государственного технологического института (технического университета). - 2014. - № 26 (52). - С. 92-94.
14. Харазов В.Г. Проблемы и пути развития интегрированных АСУ ТП // Известия Санкт-Петербургского государственного технологического института (технического университета). - 2013. - № 18 (44). - С. 089-092.
15. Шенброт И.М., Антропов М.В., Давиденко К.Я. Распределенные АСУ технологическими процессами. - М.: Энергоатомиздат, 1985. - 240 с.
16. Шенброт И.М. Распределенные АСУ ТП - АСУ нового класса // Приборы и системы управления. - 1983. - № 12. - С. 5-6.
17. Шенброт И.М., Алиев В.М. Проектирование вычислительных систем распределенных АСУ ТП. - М.: Энергоатомиздат, 1989. - 88 с.
18. Шишов О.В. Программируемые контроллеры в системах промышленной автоматизации: учебник. - М.: ИНФРА-М, 2017. - 365 с.
19. Месарович М., Мако Д., Такахара Я. Теория иерархических многоуровневых систем.
- М.: Мир, 1973. - 344 с.
20. Месарович М., Такахара Я. Общая теория систем: математические основы. - М.: Мир, 1978. - 311 с.
21. Химмельблау Д. Прикладное нелинейное программирование. - М.: Мир, 1975. - 534 с.
22. Йодан Э. Структурное проектирование и конструирование программ. - М.: Мир, 1979.
- 415 с.
REFERENCES
1. Ostrovskiy G.M., Ziyatdinov N.N., Lapteva T. V. Optimizatsiya tekhnicheskikh sistem [Optimization of technical systems]. Moscow: KNORUS, 2012, 432 p.
2. Frolov S.V. Razvitie otechestvennykh sistem promyshlennoy avtomatizatsii [The development of domestic industrial automation systems], Avtomatizatsiya i proizvodstvo [Automation and manufacturing], 2001, No. 2, pp. 22-25.
3. Shil'yak D. Detsentralizovannoe upravlenie slozhnymi sistemami [Decentralized control of complex systems]. Moscow: Mir, 1994, 576 p.
4. Kholodnov V.A., Lebedeva M.Yu. Priblizhennyy metod optimizatsii khimiko-tekhnologicheskikh sistem pri interval'noy parametricheskoy neopredelennosti [Approximate method for optimization of chemical-technological systems with interval parametric uncertainties], Izvestiya Sankt-Peterburgskogo gosudarstvennogo tekhnologicheskogo instituta (tekhnicheskogo universiteta) [Proceedings of Saint-Petersburg state technological Institute (technical University)], 2016, No. 33 (59), pp. 97-100.
5. Ostrovskiy G.M., Lapteva T.V., Ziyatdinov N.N. Proektirovanie optimal'nykh khimiko-tekhnologicheskikh sistem v usloviyakh neopredelennosti [Design of optimal chemical-technological systems in conditions of uncertainty], Teoreticheskie osnovy khimicheskoy tekhnologii [Theoretical foundations of chemical engineering], 2014, Vol. 48, No. 5, pp. 527-537.
6. Pegat A. Nechetkoe modelirovanie i upravlenie [Fuzzy modeling and control]. Moscow: BINOM. Laboratoriya znaniy, 2013, 798 p.
7. Frolova T.A., Frolov S.V., Tulyakov D.S. Reshenie interval'nykh matematicheskikh modeley tekhnologicheskikh protsessov [The solution of interval mathematical models of technological processes], Elektronnoe nauchno-tekhnicheskoe izdanie «Nauka i obrazovanie» [Electronic scientific and technical periodical "Science and education"], 2012, No. 09, pp. 343-360. Available at: http://technomag.bmstu.ru/doc/454499.html (accessed 17 February 2017).
8. Rastrigin L.A. Sovremennye printsipy upravleniya slozhnymi ob"ektami [Modern principles of management of complex objects]. Moscow: Sov. radio, 1980, 232 p.
9. Saymon G. Nauki ob iskusstvennom [Science of artificial]: translation from English. 2 nd ed. Moscow: Editorial URSS, 2004, 144 p.
10. Kafarov V.V., Meshalkin V.P. Analiz i sintez khimiko-tekhnologicheskikh system [Analysis and synthesis of chemical-technological systems]. Moscow: Khimiya, 1991, 432 p.
11. Kafarov V.V., Dorokhov I.N. Sistemnyy analiz protsessov khimicheskoy tekhnologii. Osnovy ctrategii [System analysis of chemical engineering processes. The basics of the strategy]. Moscow: Nauka, 1976, 499 p.
12. Kharazov V.G. Integrirovannye sistemy upravleniya tekhnologicheskimi protsessami: uchebnoe posobie dlya vuzov [Integrated systems of control of technological processes: textbook for universities]. St. Petersburg: Professiya, 2013, 655 p.
13. Kharazov V.G. Intellektual'nye pribory i sistemy upravleniya [Intelligent devices and control systems], Izvestiya Sankt-Peterburgskogo gosudarstvennogo tekhnologicheskogo instituía (tekhnicheskogo universiteta) [Bulletin of Saint-Petersburg state technological Institute (technical University)], 2014, No. 26 (52), pp. 92-94.
14. Kharazov V.G. Problemy i puti razvitiya integrirovannykh ASU TP [Problems and ways of development of an integrated process control system], Izvestiya Sankt-Peterburgskogo gosudarstvennogo tekhnologicheskogo instituta (tekhnicheskogo universiteta) [Bulletin of Saint-Petersburg state technological Institute (technical University)], 2013, No. 18 (44), pp. 089-092.
15. Shenbrot I.M., Antropov M.V., Davidenko K.Ya. Raspredelennye ASU tekhnologicheskimi protsessami [Distributed automatic control system of technological processes]. Moscow: Energoatomizdat, 1985, 240 p.
16. Shenbrot I.M. Raspredelennye ASU TP - ASU novogo klassa [Distributed - ASU new class], Pribory i sistemy upravleniya [Devices and control systems], 1983, No. 12, pp. 5-6.
17. Shenbrot I.M., Aliev V.M. Proektirovanie vychislitel'nykh sistem raspredelennykh ASU TP [Design of computer systems distributed control system]. Moscow: Energoatomizdat, 1989, 88 p.
18. Shishov O.V. Programmiruemye kontrollery v sistemakh promyshlennoy avtomatizatsii: uchebnik [Programmable controllers in industrial automation systems: tutorial]. Moscow: INFRA-M, 2017, 365 p.
19. Mesarovich M., Mako D., Takakhara Ya. Teoriya ierarkhicheskikh mnogourovnevykh system [Theory of hierarchical multilevel systems]. Moscow: Mir, 1973, 344 p.
20. Mesarovich M., Takakhara Ya. Obshchaya teoriya sistem: matematicheskie osnovy [General systems theory: mathematical foundations]. Moscow: Mir, 1978, 311 p.
21. Khimmel'blau D. Prikladnoe nelineynoe programmirovanie [Applied nonlinear programming]. Moscow: Mir, 1975, 534 p.
22. Yodan E. Strukturnoe proektirovanie i konstruirovanie program [Structural design and construction programs]. Moscow: Mir, 1979, 415 p.
Статтью рекомендовал к опубликованию д.т.н., профессор С.П. Тарасов.
Фролов Сергей Владимирович - Федеральное государственное бюджетное образовательное учреждение высшего образования «Тамбовский государственный технический университет»; e-mail: [email protected]; 392000, г. Тамбов, ул. Советская, 106; тел.: +79204817586; кафедра биомедицинской техники; д.т.н.; профессор; зав. кафедрой.
Фролова Татьяна Анатольевна - e-mail: [email protected]; тел.: +79202307269; кафедра компьютерно-интегрированных систем в машиностроении; к.т.н.; доцент.
Frolov Sergey Vladimirovich - Tambov State Technical University; e-mail: [email protected]; 106, Sovetskaya street, Tambov, 392000, Russia; phone: +79204817586; the department of biomedical engineering; dr. of eng. sc.; professor; head of department.
Frolova Tatyana Anatolievna - e-mail: [email protected]; phone: +79202307269; the department of computer-integrated systems in mechanical engineering; cand. of eng. sc.; associate professor.