УПРАВЛ1ННЯ У ТЕХН1ЧНИХ СИСТЕМАХ
УПРАВЛЕНИЕ В ТЕХНИЧЕСКИХ СИСТЕМАХ
CONTROL IN TECHNICAL SYSTEMS
УДК 681.5.013
Дорофеев Ю. И.
Канд. техн. наук, доцент, доцент кафедры системного анализа и управления Национального технического
университета «Харьковский политехнический институт», Харьков, Украина
ДЕСКРИПТОРНЫЙ ПОДХОД К СИНТЕЗУ ОГРАНИЧЕННОГО РОБАСТНОГО ГАРАНТИРУЮЩЕГО УПРАВЛЕНИЯ ЗАПАСАМИ С ИСПОЛЬЗОВАНИЕМ ПАРАМЕТРИЗОВАННОЙ ФУНКЦИИ ЛЯПУНОВА
Решена задача синтеза робастного гарантирующего управления запасами для сетей поставок с параметрической структурной неопределенностью в условиях действия неизвестного, но ограниченного внешнего спроса и наличия несимметричных ограничений на значения состояний и управлений. Закон управления строится в виде линейной динамической обратной связи по сигналу невязки между наличными и страховыми уровнями запаса ресурсов. Для подавления влияния внешних возмущений, моделирующих изменения спроса, одновременно с обеспечением робастной устойчивости замкнутой системы применен метод инвариантных эллипсоидов, который был усовершенствован посредством использования дескрипторного описания системы и построения параметризованной функции Ляпунова, что позволяет уменьшить степень влияния неопределенности значений транспортных запаздываний на результат синтеза управления. С помощью математического аппарата линейных матричных неравенств задача синтеза управления представлена в виде последовательности задач полуопределенного программирования, для решения которых применяются свободно распространяемые специализированные программные пакеты. В рамках предложенной методики возможен выбор оптимальных значений страховых запасов ресурсов, поскольку рассмотренное решение задает, фактически, алгоритмическую зависимость между уровнем страховых запасов и оптимальным значением критерия качества. Рассмотрен численный пример.
Ключевые слова: управление запасами, робастное управление, метод инвариантных эллипсоидов, дескрипторный подход, параметризованная функция Ляпунова, линейное матричное неравенство.
НОМЕНКЛАТУРА
«+» - псевдообращение Мура-Пенроуза;
0пхп - нулевая матрица соответствующей размерности;
А - матрица динамики расширенной модели сети поставок;
В - матрица влияния управлений расширенной модели сети поставок;
В{ - матрицы влияния управлений исходной модели сети поставок;
С - матрица выходов расширенной модели сети поставок;
Б - множество допустимых значений внешних возмущений;
Е - матрица влияния возмущений исходной модели сети поставок;
О - матрица влияния возмущений расширенной модели сети поставок;
1п - единичная матрица соответствующей размерности;
© Дорофеев Ю. И., 2015
БОТ 10.15588/1607-3274-2015-4-13
К (к) - нестационарная матрица коэффициентов обратной связи;
и - множество допустимых значений управляющих воздействий;
- положительно определенная диагональная весовая матрица критерия качества;
Wu - положительно определенная диагональная весовая матрица критерия качества;
X - множество допустимых значений состояний; Стах - вектор верхних граничных значений внешнего спроса;
Ст1П - вектор нижних граничных значений внешнего спроса;
С (к) - вектор внешних возмущающих воздействий; к - номер дискретного интервала времени; т - размерность вектора управляющих воздействий; п - количество узлов сети поставок;
и тах - вектор граничных значений управляющих воздействий;
и (к) - вектор управляющих воздействий; q - размерность вектора внешних возмущений; г - количество узлов сети, интервалы запаздывания которых варьируются;
Хтах - вектор граничных значений состояний; х(к) - вектор состояний исходной модели сети поставок;
е - малая положительная константа; ^ (к) - вектор состояний расширенной модели сети поставок;
Л - максимальное значение интервалов запаздывания материальных потоков в сети;
ЛМН - линейное матричное неравенство; ПФЛ - параметризованная функция Ляпунова.
ВВЕДЕНИЕ
Задача управления запасами возникает в системах производства-хранения-распределения ресурсов, когда с целью удовлетворения потребительского спроса создаются запасы материальных ресурсов. Примерами могут служить производственные, транспортные системы, системы распределения ресурсов (воды, электроэнергии) и т. п., в состав которых входят поставщики сырья, производственные узлы, хранилища ресурсов и продавцы конечной продукции, тогда как потребители, формирующие заказы на поставку ресурсов, рассматриваются в качестве источников внешних возмущений. Существуют различные типы топологии подобных систем, которые определяются взаимным размещением произ-водствен-ных звеньев, промежуточных складов и потребителей. Если некоторые виды сырья или полуфабрикатов используются в нескольких процессах, проходящих одновременно, система приобретает эшелонированную структуру, вследствие чего рассматриваемые системы называют сетями поставок.
Предполагается, что каждый узел сети поставок в реальном времени принимает заказы от узлов, являющихся потребителями его продукции, а также от внешних потребителей, и формирует заказы узлам, которые являются для него поставщиками ресурсов. Управление запасами заключается в определении моментов времени и размеров заказов на их восполнение. Выбор стратегии управления запасами определяется характером внешнего спроса. На практике зачастую нет оснований для того, чтобы рассматривать спрос в качестве случайных, либо гармонических, либо убывающих с течением времени внешних возмущений - какая-либо дополнительная информация, кроме той, что внешний спрос является ограниченным, отсутствует.
Спецификой задачи управления запасами является наличие запаздываний по управлению, обусловленных задержками в пополнении запасов относительно моментов формирования заказов. В процессе функционирования сети поставок величины запаздываний могут отличаться от своих номинальных значений. В результате возникает необходимость обеспечения робастности системы управления запасами относительно возможных вариаций значений внешнего спроса и величин запаздывания.
Другой особенностью рассматриваемой задачи является необходимость учета структурных ограничений на объемы хранилищ и размеры заказов. В теории управления традиционно рассматривают ограничения, заданные в какой-либо норме. Однако, для задач управления запасами характерно требование неотрицательности значений переменных, что приводит к необходимости учета несимметричных ограничений на значения управляющих воздействий и состояний.
Одним из наиболее распространенных подходов к синтезу ограниченного стабилизи-рующего управления является построение квадратичной функции Ляпунова и использование математического аппарата ЛМН. Основным недостатком такого подхода является консерватизм полученных результатов [1], который проявляется в том, что с практической точки зрения получаемые границы робастности оказываются неоправданно заниженными. Причиной является использование единой функции Ляпунова для всех возможных вариантов реализации неопределенности модели системы и внешних воздействий.
Для уменьшения степени консерватизма используют параметризованную функцию Ляпунова. Однако, при этом полученные ЛМН содержат произведение нестационарной матрицы динамики системы на матрицу, которая участвует в построении функции Ляпунова. Для преодоления указанного недостатка применяется деск-рипторный подход, который позволяет добиться разделения матрицы динамики системы и матрицы Ляпунова, что ведет к существенному уменьшению степени консерватизма результатов синтеза управления.
Целью работы является синтез робастной стратегии управления запасами гарантированной стоимости, которая строится на основе дескрипторного подхода и может использоваться для определения в каждый момент времени размеров заказа ресурсов с учетом несимметричных ограничений на их значения в виде функции от уровней запаса ресурсов в узлах системы, которые позволяют удерживать состояния в ограниченном компактном множестве при любых допустимых возмущениях.
1 ПОСТАНОВКА ЗАДАЧИ
Для математического описания управляемой сети поставок используется дискретная модель в пространстве состояний, уравнения которой описывают изменение уровня запасов каждого вида ресурсов с течением времени. В качестве переменных состояний рассматрива-ют-ся наличные уровни запаса ресурсов. Управляющими воздействиями являются размеры заказов на поставку ресурсов, формируемые узлами в текущем периоде. Размеры спроса на ресурсы, поступающие из внешней среды, выступают в качестве внешних возмущений.
Для описания запаздываний используется модель дискретной задержки, поскольку предполагается, что номинальные значения длительности транспортировки и переработки ресурсов известны и кратны некоторому периоду дискретизации. Также предполагается, что структура сети известна, а состояния доступны непосредственному измерению. Тогда математическая модель сети поставок задается разностным уравнением с запаздыванием:
х(к +1) = х(к) + Б{и(к - 0 + Еа (к). (1)
р-К8К 1607-3274. Радюелектронжа, iнформатика, управлiння. 2015. № 4 е-ЕЗБЫ 2313-688Х. Яаёю Шео^оп^, Сотр^ег Баепое, СоПго1. 2015. № 4
Предполагается, что для каждого узла сети заданы максимально допустимые уровни запаса ресурсов и максимальные размеры заказов. Тогда в процессе функционирования системы должны выполняться ограничения:
х(к) е X = {х е Я п :0 < х < хтах ), и (к) е и = {и е Я т :0 < и < и тах).
(2)
Будем предполагать, что внешние возмущения удов -летворяют ограничениям:
С (к) е Б = {с е Я 9 :0 < сСт1п < С < сСг
}.
Выполним преобразование модели (1) к стандартному виду без запаздывания на основе расширения вектора состояний [2] путем включения в него векторов, определяющих размеры ранее заказанных ресурсов, находящихся в процессе транспортировки и переработки:
хТ (к), иТ (к -1), иТ(к - 2),..., иТ(к - Л)]. Тогда
«к) =
уравнения расширенной модели примут вид:
Е,(к +1) = А^(к) + Ви (к) + ОС (к), х(к) = СЩ), (3)
где матрицы А е Я, В е Я*хт, О е Я *х9, с е Я пх*, N = п + тЛ имеют соответствующую блочную структуру [3].
Рассмотрим построение матрицы динамики А расширенной модели в том случае, когда величины запаздывания
управляемых потоков Л,, , = 1, п отличаются от своих номинальных значений. В этом случае матрица становится нестационарной и в каждый момент времени к > 0 может принимать какое-либо значение из множества
А(8) = {аеЯ: А = Ао +^=10/(к)А(г), Э(к)е©), (4)
где Ь = 2; 8,- (к), , = 1, Ь - набор параметров, которые описывают структурную неопределенность модели и удовлетворяют следующим требованиям:
© = {8 е ЯЬ : 0г (к) > 0, £^8, (к) = 1)
(5)
Таким образом, модель сети в условиях неопределенности значений запаздывания управляемых потоков может рассматриваться как выпуклый многогранник, который задается списком вершин
{А(1), В, О, С) А(2), В, О, С)..., (а(Ь), В, О, С)), и может быть представлена в виде модели с параметрической структурной неопределенностью следующего вида:
Е,(к +1) = А(8)^(к) + Ви (к) + ОС (к),
х(к) = С^(к), А(8) е О = Со{(1), ..., А(Ь)), (6)
где Со{- ) - выпуклая оболочка; А(,), , = 1,Ь - , -я вершина выпуклого множества О.
Запишем критерий качества в случае бесконечного временного горизонта в виде:
J«>(к) = ЕГ=0(())Т((к)-5*)+иТ(к)WuU(к) 1, (7)
*
где 0 X W5е Я *х* ,0 X Wu е Я тхт; 5 = £
*Т *Т ..., х ]
Л+1
*
составной вектор, у которого компоненты вектора х определяют размеры страховых запасов ресурсов в узлах сети и вычисляются на основании верхних граничных значений внешнего спроса с учетом запаздываний и продуктивной модели Леонтьева:
= (( - П) 1 С, С, =
, = 1,9,
|л,ст _
I 0, , = 9 +1, п,
(8)
где П е Япхп - продуктивная матрица, значение элемента Пу которой равно количеству единиц ресурса ,, необходимому для производства единицы ресурса ].
Первое слагаемое в выражении (7) определяет размеры штрафов за отклонение текущих уровней запаса ресурсов от страховых, второе - стоимость производства и транспортировки ресурсов.
Для системы (6) с параметрической неопределенностью (4), (5) рассматривается задача синтеза робастной стратегии управления запасами, которая для любого допустимого спроса С (к) е Б у к > 0 обеспечивает:
1) полное и своевременное удовлетворение спроса, то есть выполнение первого из ограничений (2) на значения состояний;
2) робастную устойчивость замкнутой системы при ограничениях (2) на значения управлений;
3) гарантированную стоимость управления, которая означает, что значение критерия качества (7) не превысит некоторого граничного значения.
В качестве дополнительного условия выдвигается требование снижения степени влияния изменений неопределенного вектора параметров 6(к) на результат синтеза управления.
2 ОБЗОР ЛИТЕРАТУРЫ
Анализ различных подходов к управлению запасами можно найти в работе [4] и обширной библиографии к ней. Среди многообразия моделей управления запасами выделяют два основных типа: модель критических уровней и модель периодической проверки. В первом случае предполагается непрерывный контроль за состоянием запасов и размещение заказов фиксированного размера при снижении текущих запасов до некоторых критических уровней. Второй тип модели предполагает проверку уровня запасов через равные промежутки времени и размещение заказа, размер которого определяется в соответствии с выбранной стратегией.
В данной работе рассматривается модель периодической проверки, а задача управления запасами сформулирована как задача подавления влияния неслучайных ограниченных внешних возмущений, методы решения которой рассмотрены в работе [5].
Одним из подходов к данной проблематике в теории робастного управления является концепция инвариантных множеств [6]. Среди различных форм инвариантных множеств особо выделяются эллипсоиды вследствие их
Т
х
простой структуры и прямой связи с квадратичными функциями Ляпунова. При этом для синтеза оптимального регулятора требуется решить эквивалентную задачу поиска наименьшего по некоторому критерию инвариантного эллипсоида замкнутой динамической системы. В работе [5] на основе техники ЛМН устанавливается достаточное условие устойчивости замкнутой системы -это существование квадратичной функции Ляпунова, построенной на решениях системы.
После того, как были развиты вычислительные методы, основанные на идеях выпуклой оптимизации, и для их реализации были разработаны соответствующие алгоритмы и программное обеспечение [7], техника ЛМН используется в качестве общего метода анализа и синтеза динамических систем как в непрерывном, так и в дискретном времени. Однако, применение указанного подхода приводит к консервативным результатам. Это означает, что максимально допустимая величина неопределенности, при которой сохраняется робастность системы, определяется наихудшим элементом семейства. Иными словами, подход рассчитан на наихудшую возможную неопределенность, реализация которой на практике может быть крайне маловероятной.
Для преодоления указанного недостатка предлагается использовать так называемую ПФЛ (см., например, [8]). Исследования показали, что условия устойчивости, полученные на основе ПФЛ являются менее консервативными, чем те, при получении которых использовалась функции Ляпунова, не зависящая от параметра. Необходимо отметить, что большинство работ, авторы которых используют ПФЛ, посвящены анализу и синтезу динамических систем в непрерывном времени. Среди работ, в которых предлагается применять ПФЛ для систем в дискретном времени, следует отметить [9].
В настоящее время подход на основе использования ПФЛ стал мощным инструментом для анализа и синтеза линейных систем с неопределенностью различного рода. В работе [10] предложено расширение указанного подхода на основе дескрипторного описания системы, которое первоначально было предложено для исследования устойчивости и синтеза управления в системах с запаздыванием. Дескрипторный подход обладает двумя основными преимуществами: во-первых, введение дополнительных ослабляющих переменных позволяет значительно уменьшить консерватизм результатов; во-вторых, подход может быть применен не только для анализа устойчивости систем с неопределенностью, заданной в виде многогранника, но и для решения задачи синтеза оптимального регулятора.
Однако, полученные результаты напрямую неприменимы к задаче управления запасами в сетях поставок, поскольку рассматриваемая модель системы не содержит внешних возмущений, авторы не учитывают ограничения на значения состояний и управлений, а также не вводят критерий качества, позволяющий оценить стоимость полученного управления.
В результате возникает необходимость расширения метода инвариантных эллипсоидов для задач управления запасами на основе дескрипторного подхода с использованием параметризованной функции Ляпунова.
3 МАТЕРИАЛЫ И МЕТОДЫ
Выполним аппроксимацию множества D значений внешнего спроса эллипсоидом наименьшего объема, который задается уравнением
Е(/, Qd) = | й е Яч : (й(к) - /)ТQ-1(((к) - /)< ^ (9)
Матрица Qd е Я 4x4 и вектор й * е Я ч координат центра эллипсоида определяются в результате решения задачи полуопределенного программирования аналогично тому, как это сделано в работе [11].
Будем строить закон управления в виде линейной динамической обратной связи по сигналу рассогласования между наличными и страховыми уровнями запаса ресурсов
и(к) = К (к )(;(к) - 4*). (10)
Тогда расширенную модель замкнутой системы для управления (10) можно представить в виде
4(к +1) = Лг (к, 0)(^(к) - 4*)+ Л(9)4* + Ой (к),
х(к) = С4(к), Лг (к, 0) = Л(0) + БК(к), Л(0) е О. (11)
Задача синтеза управления эквивалентна решению минимаксной задачи
(
и (к) = ащ тт
и (к )еи
Л
тах J ж (к)
й (к )еЕ (й *,Qd), Л(0)еО
. (12)
Выполним преобразование системы (6) с помощью дескрипторной системы следующего вида
0 N х N
0 N xN 0 N xN
4(к +1)
У (к +1).
0 N xN ^
Л(9) - ^
4(к)
У(к)
0 N х N
Б
и (к) +
0 N xN
О
й (к).
(13)
Пусть
>!(к, 9) 0 N xN' Р2(к ) Рз(к )
где Р) е Я вектора
Р(к, 9) =
NxN, I = 1,2,3, Р1 (к, 9) - аффинная функция от
параметров
ь
9 такая, что Р1г (к) ^ 0,1 = 1ь,
Рх(к,9) = РТ(к,9) = ^^=191Ри(к)
Рз(к) = Р,Т(к).
Определим модифицированную параметризованную функцию Ляпунова, которая построена на решениях системы (11), в следующем виде
V («к) - 4*, 9)= («к) - 4* )т р (к, 9)(к) - 4*). (14)
Динамическая система (11) с параметрической неопределенностью (4), (5) является робастно устойчивой, если
V(4(к)-4*,9)>0 и ау(4^)-4*,е)= V(4^+1)-4*,е)-
-V(4(к)-4*,9)<0 У9 е© при (4(к)-4*)*0.
+
+
Введем
обозначения
a (9) =
NxN
I
N
G _
a(9) - in l(k) =
NxN
G
A(9) =
4(k) y(k)
E _
0 N xN
A(9) - In_
IN 0 NxN 0 NxN 0 NxN
B _
0 NxN
B
, Ч * _ .ч* 1 с *
.Ч _
Вычислим первую по к разность ПФЛ (14) в силу системы (11) с учетом преобразования (13):
лг() - Ч*,е)=((к+1) - 1 ^(к^к+1) - 1 *)--(- 1*)т Ещ,о(к) - !*)=
(- 1*)т ^к,е)К(к) -1*)=
4(k+1) - Ч
0Nx1
P(k,9)
4(k+1)-ч
0Nx1
y(k) - Ч
A(9)4(k) + Bu(k) + Gd(k) - y(k) y(k)-4
A(9)4(k) + Bu(k)+Gd(k) - y(k) -((k) - ff EP(k,0)((k) - 4*)=
>I<
y(k) - 4
P(k,9)x
P(k,9)x
A(9)(4(k) - 4*)+Bu(k) + Gd(k) - y(k) + 4* +(((9) - In ) y(k)-4
A(9)(4(k) - 4*)+Bu(k)+Gd(k) - y(k) + 4* +(((9) - In)
-(4(k) - ljT EP(k,0)((k) - 4*)= ()() - 4*)+Bu(k) + Gd(k) + ^4(9)4*)T(((k,9)+ PT(k,9)):
x (l(0)((k) - Ч*)+Bu(k) + Gd(k) + A(0)i;*)-(4(k) - Г) EP(k,9)((k) - Ч*)) - 4*)T PT(k,9)ET(((k) - ЧТ)
Потребуем, чтобы значение ПФЛ (14) с течением времени убывало с некоторой гарантированной скоростью, которая определяется значением показателя качества (7):
Мч<
,(k)-Ч )+Bu(k)+Gd(k)+A(9)4 ) \P(k,9)+PT(k,9 x ((9)((k) - Ч*)+Bu(k)+Gd(k)+A(9)4 *)--((k) - 4*)T Epk,9)((k) - 4*)-(((k) - 4*)T PT(k,9)ET (4(k) - Ч*)
<-i((k) - Ч*)® - Ч*)+uT(k)WUu(k)
Если последнее неравенство выполняется, то можно показать, что ПФЛ (14) Vk > 0 определяет верхнее граничное значение критерия качества (7):
max J„ (k) < Vk(k) б).
d(k)eE(d ,Qd),A(9)eQ
Тогда задача (12) эквивалентна задаче минимизации
ПФЛ (14) u(к) = arg min V\t(k)б), для решения
и (k)eU V '
которой применим метод инвариантных эллипсоидов [4]. Кратко изложим идеи метода.
Эллипсоид, описываемый уравнением
*, Q(k, 6)) = { % е r" : (¡(к) - )T Q-1 (k, 6)\((k) - ) < ), (15)
называется инвариантным по состоянию для системы (11), если из условия ¡(0) е £(%*, q(0,6)) следует, что
¡(к) е e(%*,Q(k,6)) Vk > 0. Другими словами, любая траектория системы, начавшись в инвариантном эллипсоиде, остается в нем для любого момента времени k > 0.
Эллипсоид (15) может рассматриваться в качестве аппроксимации множества достижимости замкнутой системы (11), то есть позволяет характеризовать влияние внешних возмущений и неопределенности параметров модели на траекторию замкнутой системы. Тогда минимизация в некотором смысле инвариантного эллипсоида (15) соответствует робастному управлению системой (11).
Сравнение выражений (14) и (15) позволяет сделать вывод, что если выполняется тождество
P[(k, 6) = Q _1(k, 6), то эллипсоид (15) представляет собой множество, находящееся внутри поверхности уровня ПФЛ (14).
Тогда задача робастной стабилизации заключается в вычислении в каждый момент времени k > 0 матрицы K (k) такой, чтобы регулятор (10) обеспечивал минимизацию по некоторому критерию эллипсоида (15) при ограничениях (2). Выберем в качестве критерия сумму квадратов полуосей эллипсоида, то есть след матрицы Q(k, 6). Тогда целевая функция оптимизационной задачи будет иметь вид
Xf=1trace (p- 4k))
— mm,
(16)
то есть является нелинейной. Однако, задача (16) эквивалентна следующей задаче:
trace (И, (k)) — min
при ограничениях
И, (k) IN ' In Рц (k)
(17)
h 0, И, (k) ^ 0,
: RNxN, i _ 1,l.
Н, (к )<
Указанная эквивалентность доказывается с помощью неравенств Н^ (к) У Гц (к), / = 1, Ь и леммы Шура [8].
Тогда результат решения задачи синтеза ограниченного робастного гарантирующего управления запасами в сетях поставок на основе параметризованной функции Ляпунова с помощью дескрипторного подхода может быть представлен в виде следующей теоремы.
Теорема. Рассмотрим систему (11) с параметрической неопределенностью (4), (5) и ограничениями (2), и
пусть матрицы S(к) и У (к) получены в результате решения оптимизационной задачи (17) при ограничениях
T
на матричные переменные р (к) = РцТ (к) е Я Иг (к) е Я *х*, , = 1,1, Р2(к) е Я* х*, Б (к) е Я тхт:
Г (к) е Я
тх*.
2 ('') е Я(+3?+т)х(7*+3д+т), , = II;
2 (/) = 1 (■/7)> , = 1, Ь -1, / = г +1, Ь и скалярный параметр а( к )
2 (11) 2 (12) 2 (21) 2 (22)
2 (Ь1) 2 (12)
2(1Ь) 2(2Ь)
2
(ЬЬ)
< 0,
а(к) > 0,
Р2(к) ВВТ ВБ (к) Б (к)
Б т(к) В Т
^ 0,
(18)
Р1г (к) ^ 0, И, (к) ^ 0,
с
С Т Рц (к )
И, (к) I*
. I* Р, (к). ^ 0, , = 1"!,
^ 0,
(19)
ГТ(к)еет (((к) -4*) Гт(к)
Г (к)
Б (к)
Г т(к )и тах (((к) -(*) Г Т(к)
Г (к)
Б(к)
^ 0,
^ 0,
(20)
Где Qx е Япхп - матрица эллипсоида наименьшего объема, который аппроксимирует множество X; ет = [1,1,...,1]Т е Ятх1;
Если задача минимизации линейной функции (17) при ограничениях (18)—(20), которая является задачей полуопределенного программирования, имеет решение, то: 1) для любого начального состояния х(0) > Хтах, и (к) = 0 тх1 У к < 0 и У4(9) е а также внешнего возмущения d(к) е Е(С *, Qd ) система (11) с параметрической неопределенностью (4), (5) является робастно устойчивой при ограничениях (2); 2) среди всех линейных управлений вида (10) регулятор с матрицей
К (к) = Б _1(к )7(к)
(21)
доставляет минимум по критерию следа матрицы инвариантному эллипсоиду (15) для замкнутой системы (11) в момент времени к.
Из-за ограниченности объема статьи доказательство теоремы опускаем.
Первое из ЛМН (18) гарантирует стабилизацию замкнутой системы (11) при УС (к) е Е (С *, Qd) и У4(9) е Последнее из ЛМН (19), а также ЛМН (20) обеспечивают выполнение заданных ограничений (2). Отметим, что именно наличие ЛМН (20) приводит к необходимости использования нестационарной обратной связи, поскольку матрицы неравенств зависят от текущего значения вектора состояний 4 (к).
4 ЭКСПЕРИМЕНТЫ
В качестве примера рассмотрим сеть поставок, которая изучалась в работе [2]. Модель сети описывается
графом О = ({1,2,3}, {(2,1), (2,3), (3,1)}) . Заданы значения
времени выполнения заказа в узлах сети: Т = = 2, Т3 = 1; и времени транспортировки ресурсов между уз -лами сети: Т21 = Т31 = Т23 = 1.
Представим управляемые потоки и1 и и3, описывающие сборочные процессы, в виде гипердуг, а также добавим поток и 2, который описывает поставки сырья извне (см. рис. 1). Дуги С! и С 2, изображенные пунктирными линиями, представляют внешний спрос. Значение времени транспортировки и количество единиц продук-
-р(к) /°)Тр2(к)+ГТ(к)ЕТ' 0*ху #Р(к)+ГТ(к В 0№у Щ Гтт
рТ(к)А0 +ВГк) 4р (к}-2(2(к)+Р^(к)) Р(к( -1* Р^^кО -2Р2(к) 0Ыху 0ЫхЫ
0ЫхЫ ( -I*) Р2(к) / \Т ( -I*) ) 0ЫхЫ
ОТР2(к) 0уху 03»? ОТР2(к) Оуху
ОТР2(к) 0зх* 0 0уу Чуху К ОТР2(к) Оуху
РТкА +ВГк) -2р(к) Р^кОРТкО Рк Р^
0УхУ 0уху фТР2(к) ОкЦд Оуху
Щ 0Щ 0^ 0ЫхЫ 0Ыху 0ЫхЫ
ТО 0лш 0лх? 0¥х? 0ЫхЫ 0Ыху 0ЫхЫ ^'и
р-ЕЗБЫ 1607-3274. Радюелектронжа, шформатика, управлшия. 2015. № 4 е-ЕЗБЫ 2313-688Х. Ид^о Ше^гоп^, Сошриег Баепое, Сопйо1. 2015. № 4
(0)[1] и2 Д а2
(1)[2]
к I,
= 2,0,) = 1,12
\(1)[1]
и, /
и3 1 "V
(1)[2]
Рисунок 1 - Графическое представление модели сети поставок ции Пу, которое требуется в соответствии с технологическим процессом, указаны для каждого управляемого потока в круглых и квадратных скобках, соответственно. Возле каждого узла в круглых скобках указаны значения времени выполнения заказа Т,. Специфика рассматриваемой системы в том, что на узел 1 действует только внешний спрос; на узел 2 действует как внешний, так и внутренний спрос со стороны узлов 1 и 3; на узел 3 -только внутренний спрос со стороны узла 1.
По формуле Л, = шах7у,г + 7, ), у = 1,3, у Ф )} определим величины запаздывания управляемых потоков для всех узлов сети, в результате получим л = 3. Тогда размерность расширенной модели сети равна N = 12.
Заданы максимальные вместительности хранилищ уз -лов сети и объемы транспортировок хшах = [120,672,240]т, ишах = [25,130, 55]т, а также граничные значения внешнего спроса а™п = [7, 6]Т, ашах = [20,18]т
Пусть время транспортировки ресурсов между узлами 2 и 3 в процессе функционирования сети может увеличиваться на один период, т.е. 7^3 е { 2}. Тогда величина запаздывания управляемых потоков узла 3 принимает значение из множества л3 е {2, 3}. В результате
получим А(9) е о = Со{4(1), А(2)}, где
А(1) =
0 0 0 "0 0 0 1 0 0 "1 0 0~
В)= 0 1 0 , В2»= 0 1 0 , В1= 0 0 0 В2)= 0 0 0
0 0 1 0 0 0 0 0 0 0 0 1
(варьируемые элементы матриц подчеркнуты).
В соответствии с (8) уровни страховых запасов равны
х* = [60, 336,120]Т Следуя [12], вычислим параметры эллипсоида, аппроксимирующего множество значений внешнего спроса: = diag(84,5, 72,0) ,
а * =[13 ,5,12,0] ; а также матрицу эллипсоида Qx, аппроксимирующего множество X допустимых значений
состояний. В качестве начального состояния выбраны
*
значения страховых запасов х(0) = х , а значения диаго-
нальных весовых матриц равны ки ]УУ. = 0,1, у = 13.
5 РЕЗУЛЬТАТЫ
Моделирование осуществлялось в течение 15 периодов. Результаты моделирования получены путем численного решения последовательности задач (17) при ограничениях (18)-(20) в условиях скачкообразно изменяющегося внешнего спроса с помощью свободно распространяемого пакета СУХ [12]. Результаты представлены на рис. 2 - рис. 4, где а - значения наличного и страхового уровней запаса; б - значения внешнего спроса и объемов заказа.
6 ОБСУЖДЕНИЕ
Вначале наблюдается переходный процесс, обусловленный тем, что в момент начала моделирования каналы транспортировки ресурсов не были загружены. В 4, 9 и 14 периодах значение матрицы динамики выбиралось
равным в остальных периодах - следствием чего являются скачкообразные изменения объемов заказов в указанных периодах, особенно заметные для 1 и 3 узлов. В процессе моделирования фазовая траектория замкнутой системы не выходит за пределы инвариантных эллипсоидов, размеры которых зависят от выбранных значений весовых матриц к и к, определяющих критерий качества системы.
13 В В) В1 1 13 В В2) В2)1 во * 20
0x3 0м 0м 03x3 , А2)= 03x3 03x3 03x3 03x3 10
0x3 13 03x3 03x3 03x3 13 03x3 03x3
0x3 0м 13 03x3 03x3 03x3 13 03x3 0 0
5 10
Номер дискретного интервала к
а
&
п_
5 10
Номер дискретного интервала к
б
Рисунок 2 - Графики переходных процессов для узла 1
а
70
60
50
40
30
15
30
25
20
15
10
5
0
0
5 10
Номер дискретного интервала k
а
8 CL
5
________1________ 1 1 ___ --- от
1 1 1
1
1
1 1 1
1 1 1
1 1 1
f----------- 1
о
15
5 10
Номер дискретного интервала к б
Рисунок 3 — Графики переходных процессов для узла 2 130
5 10
Номер дискретного интервала k а
| 30
£
> 20
5 10
Номер дискретного интервала k б
Рисунок 4 - Графики переходных процессов для узла 3
Результаты моделирования показали, что полученная стратегия управления запасами обеспечивает отсутствие дефицита ресурсов в сети, гарантированную стоимость управления, а также робастную устойчивость замкнутой системы при заданных ограничениях на состояния и управления. ВЫВОДЫ
В статье предложен подход к решению задачи синтеза ограниченного робастного гарантирующего управления запасами для сетей поставок с параметрической структурной неопределенностью в условиях действия неизвестного, но ограниченного внешнего спроса и наличия несимметричных ограничений на значения состояний и управлений. Закон управления строится в виде линейной динамической обратной связи по сигналу невязки между наличными и страховыми уровнями запаса ресурсов. Подход основан на использовании метода инвариантных эллипсоидов и параметризованной функции Ляпунова. Использование математического аппарата ЛМН позволило сформулировать задачу синтеза управления в виде последовательности задач полуопределенного программирования, для решения которых применяются свободно распространяемые специализированные программные пакеты. С помощью дескрип-торного описания системы удалось добиться уменьше -ния консерватизма полученных результатов, поскольку матрицы неравенств не содержат произведения нестационарной матрицы динамики системы и матрицы Ляпунова, а при вычислении матрицы коэффициентов обратной связи используются матрицы, не зависящие явно от неопределенных параметров модели. БЛАГОДАРНОСТИ
Работа выполнена в рамках госбюджетной научно-исследовательской темы Национального технического университета «Харьковский политехнический институт» «Развитие теории и методов синтеза децентрализованного робастного управления распределенными сетями поставок в условиях неопределенности» (номер государственной регистрации 0111U002285). СПИСОК ЛИТЕРАТУРЫ
1. Щербаков П. С. Приближенные методы в параметрической роба-стности линей-ных систем управления : дис. ... д-ра физ.-мат. наук : 05.13.01 / Щербаков Павел Сергеевич. - М. : , 2004. - 215 с.
2. Blanchini F. Feedback control of production-distribution systems with unknown demand and delays / F. Blanchini, R. Pesenti, F. Rinaldi, W. Ukovich // IEEE Transactions on Robotics and Automation. - 2000. - Vol. 16, No. 3. - P. 313-317.
3. Дорофеев Ю. И. Построение математических моделей управляемых сетей поставок с учетом запаздываний потоков / Ю. И. Дорофеев, А. А. Никульченко // Системш дослщження та шформацшш технологи. - 2013. - № 1. - С. 16-27.
4. Лотоцкий В. А. Модели и методы управления запасами / В. А. Лотоцкий, А. С. Мандель. - М. : Наука, 1991. - 188 c.
5. Поляк Б. Т. Управление линейными системами при внешних возмущениях: техника линейных матричных неравенств / Б. Т. Поляк, М. В. Хлебников, П. С. Щербаков. - М. : ЛЕ-НАНД, 2014. - 560 с.
6. Blanchini R. Set-theoretic methods in control / R. Blanchini, S. Miani. - Boston: Birkhauser, 2008. - 494 p.
7. Чурилов А. Н. Исследование линейных матричных неравенств. Путеводитель по программным пакетам / А. Н. Чурилов, A. B. Гессен. - СПб: С.-Петерб. гос. ун-т, 2004. - 148 с.
340
320
300
280
260
240
0
15
80
60
40
20
0
15
60
50
40
10
0
0
15
8. Feron E. Analysis and synthesis of robust control systems via parameter-dependent Lyapunov functions / E. Feron, P. Apkarian, P. Gahinet // IEEE Transactions on Automatic Control. - 1996. -Vol. 41. - P. 1041-1046.
9. Oliveira M. C. A new discrete-time robust stability condition / М. C. Oliveira, J. Bernussou, J. C. Geromel // System & Control Letters. - 1999. - Vol. 37. - P. 261-265.
10. Zhang W. Robust stability test for uncertain discrete-time systems: a descriptor system approach / W. Zhang, H. Su, Y. Liang, Z. Han // Latin American Applied Research. - 2011. - Vol. 41, No. 4. -P. 359-364.
Дорофеев Ю. I.
Канд. техн. наук, доцент, доцент кафедри системного аналiзу i управлшня Нащонального техшчного ушверситету «Харгавський полггехшчний шститут», Харгав, Укра!на
ДЕСКРИПТОРНИЙ П1ДХ1Д ДО СИНТЕЗУ ОБМЕЖЕНОГО РОБАСТНОГО ГАРАНТУЮЧОГО КЕРУВАННЯ ЗАПАСАМИ З ВИКОРИСТАННЯМ ПАРАМЕТРИЗОВАНО! ФУНКЦП ЛЯПУНОВА
Виршено задачу синтезу робастного гарантуючого керування запасами для мереж поставок з параметричною структурною невизначешстю в умовах ди невщомого, але обмеженого зовншнього попиту та наявносп несиметричних обмежень на значення сташв i керуючих дш. Закон керування будуеться у виглядi лшшного динамiчного зворотного зв'язку за сигналом нев'язки мiж гсшвковими i страховими рiвнями запасу ресурав. Для подавлення впливу зовншшх збурень, що моделюють змши попиту, одночасно iз забезпеченням робастно! стшкосп замкнуто! системи застосований метод iнварiантних елiпсоiдiв, який був удосконалений за допо-могою використання дескрипторного опису системи та побудови параметризовано! функцп Ляпунова, що дозволяе зменшити стушнь впливу невизначеностi значень транспортних зашзнювань на результат синтезу керування. За допомогою математичного апарату лiнiйних матричних нерiвностей задача синтезу управлiння представлена у виглядi послiдовностi задач нашввизначеного програмуван-ня, для вирiшення яких застосовуються вiльно поширюванi спецiалiзованi програмш пакети. В рамках запропоновано! методики можливий вибiр оптимальних значень страхових запаав ресурсiв, оскiльки розглянуте ршення задае, фактично, алгоритмiчну за-лежшсть мiж рiвнем страхових запасiв i оптимальним значенням критерiю якостi. Розглянуто чисельний приклад.
Ключовi слова: керування запасами, робастне керування, метод iнварiантних елшсодав, дескрипторний пiдхiд, параметрезована функщя Ляпунова, лiнiйна матрична нерiвнiсть.
Dorofieiev Yu. I.
PhD, Associate professor, Associate professor of department of system analysis and control, National Technical University «Kharkiv Polytechnic Institute», Kharkiv, Ukraine
DESCRIPTOR APPROACH TO THE SYNTHESIS OF CONSTRAINED ROBUST GUARANTEED INVENTORY CONTROL USING PARAMETER-DEPENDENT LYAPUNOV FUNCTION
The problem of robust guaranteed inventory control synthesis for supply networks with parametric uncertainty under action of unknown, but bounded external demand and availability of non-symmetrical constraints on the values of states and controls is solved. The control law is constructed in the form of a linear dynamic feedback with respect to deviation between cash and safety stock levels of resources. In order to suppress the influence of external perturbations, modeling changes in demand, while ensuring robust stability of the closed system the method of invariant ellipsoids is used, which has been improved through the use descriptor system approach and building parameter-dependent Lyapunov function, which reduces the impact degree of uncertainty of the transport time-delays on the result of the synthesis control. With the help of the linear matrix inequalities technique the control synthesis problem is presented as a series of semidefinite programming, for solving which is used open source software. In the framework of the proposed method is possible to choose the optimal values of safety stock levels of resources, as the resulting solution determines, in fact, an algorithmic relationship between the levels of safety stocks and the optimal value of the quality criterion. The numerical example is considered.
Keywords: inventory control, robust control, invariant ellipsoids method, descriptor approach, parameter dependent Lyapunov function, linear matrix inequality.
REFERENCES
1. Shherbakov P. S. Priblizhennye metody v parametricheskoj robastnosti linejnyh sistem upravlenija : dis. ... d-ra fiz.-mat. nauk : 05.13.01 / Shherbakov Pavel Sergeevich. Moscow, 2004, 215 р.
2. Blanchini F., Pesenti R., Rinaldi F., Ukovich W. Feedback control of production-distribution systems with unknown demand and delays, IEEE Transactions on Robotics and Automation, 2000, Vol. 16, No. 3, pp. 313-317.
3. Dorofieiev Yu. I., Nikulchenko A. A. Postroenie matematicheskih modelej upravljaemyh setej postavok s uchetom zapazdyvanij potokov, Sistemni doslidzhennja ta informacijni tehnologii, 2013, No. 1, pp. 16-27.
4. Lotockij V. A., Mandel' A. S. Modeli i metody upravlenija zapasami. Moscow, Nauka, 1991, 188 p.
5. Poljak B. T., Hlebnikov M. V., Shherbakov P. S. Upravlenie linejnymi sistemami pri vneshnih vozmushhenijah: tehnika linejnyh matrichnyh neravenstv. Moscow, LENAND, 2014, 560 p.
6. Blanchini R., Miani S. Set-theoretic methods in control. Boston, Birkhauser, 2008, 494 p.
11. Дорофеев Ю. И. Синтез системы оптимального управления запасами с дискрет-ным ПИД-регулятором с использованием ЛМН / Ю. И. Дорофеев // Збiрник наукових праць Харювського ушверситету Повпряних Сил. - Х. : ХУПС, 2014. - Вип. 4(41). -С. 34-41.
12. Grant M. CVX: MATLAB software for disciplined convex programming, version 2.0 beta / M. Grant, Boyd S // [Electronic resource]. - Access mode: http://cvxr.com/cvx.
Статья поступила в редакцию 15.10.2015.
После доработки 19.10.2015.
7. Churilov A. N., Gessen A. B. Issledovanie linejnyh matrichnyh neravenstv. Putevoditel' po programmnym paketam. Sankt-Peterburg, SPbGU, 2004, 148 p.
8. Feron E., Apkarian P., Gahinet P. Analysis and synthesis of robust control systems via parameter-dependent Lyapunov functions, IEEE Transactions on Automatic Control, 1996, Vol. 41, pp. 1041-1046.
9. Oliveira M. C., Bernussou J., Geromel J. C. A new discrete-time robust stability condition, System & Control Letters, 1999, Vol. 37, pp. 261-265.
10. Zhang W., Su H., Liang Y., Han Z. Robust stability test for uncertain discrete-time systems: a descriptor system approach, Latin American Applied Research, 2011, Vol. 41, No. 4, pp. 359-364.
11. Dorofieiev Yu. I. Sintez sistemy optimal'nogo upravlenija zapasami s diskretnym PID-reguljatorom s ispol'zovaniem LMN, Zbirnik naukovih prac' Kharkivs'kogo universitetu Povitrjanih Sil. Kharkiv, KhUPS, 2014, Vip. 4(41), pp. 34-41.
12. Grant M., Boyd S. CVX: MATLAB software for disciplined convex programming, version 2.0 beta [Electronic resource]. Access mode: http://cvxr.com/cvx.