УДК 004.85
doi:10.18720/SPBPU/2/id24-196
Мрыхин Захар Владимирович.,
студент бакалавриата
СТАБИЛИЗАЦИЯ ПРОГРАММНЫХ ДВИЖЕНИЙ САУ С ИНТЕРВАЛЬНЫМИ ОГРАНИЧЕНИЯМИ
Россия, Санкт-Петербург, Санкт-Петербургский политехнический университет Петра Великого,
mrykhin.zv@edu.spbstu.ru
Аннотация. В работает решается задача стабилизации программных движение САУ с интервальными ограничениями на управление. Производятся преобразования и погружение модели в оптимизационную задачу, которая решается точным алгоритмом.
Ключевые слова: стабилизация программных движений, система автоматического управления, оптимизация, квадратичный функционал, исследование операций.
Zakhar V. Mrykhin,
Graduate student
STABILIZATION OF PROGRAM MOVEMENTS OF ACS WITH INTERVAL RESTRICTIONS
Peter the Great St. Petersburg Polytechnic University, St. Petersburg, Russia,
mrykhin.zv@edu.spbstu.ru
Abstract. The article solves the problem of stabilizing the software movement of the ACS with interval control restrictions. Transformations are performed and the model is immersed in an optimization problem, which is solved by an exact algorithm.
Keywords: stabilization of program movements, automatic control system, optimization, quadratic functional, operations research.
Введение
Стабилизация программных движений для объектов управления является важной задачей ТАУ и исследования операций. В данной работе вычисляется минимум квадратичного функционала с интервальными ограничениями, что обеспечивает стабилизацию программных движений объекта управления.
1. Постановка задачи
Пусть объект САУ задаётся следующим разностным уравнением:
**+i = Нхк + Fuk, (1)
где к — дискретный момент времени, хи — вектор состояния в к-й момент времени, ии — управление, действующее в к-й момент времени, а Н и Е — матрицы, характеризующие САУ. Таким образом, зная текущее положение, а также задав управление, действующее на неё, можно получить, где объект будет находиться в следующий момент.
Можно записать, как будет выглядеть задача в общем случае:
и
: а^шт \ ф =
тпд - X
Ли х,
к+1
Хк+1 = Нхк + ЕИк ( ,
без ограничений на управление и соответственно
2
и* = аг^тт |ф =
хпд - х хк хк+1
Хк+1 = Нхк + ¥ик'и- ^ < иг (,
(2)
(3)
с ограничением на управление.
В этой задаче мы ищем оптимальное управление, которое будет минимизировать функционал нормы разности между х™ и хк+1. х™ — вектор программных движений, то состояние системы, к которому она должна стремиться. и' и и+ — соответственно нижняя и верхняя границы управления.
2. Решение задачи
Задача без ограничений будет решаться не сложно. Преобразуем функционал к следующему виду:
ф =
гпд _ г
хк хк+1
= (ХкЯ хк +1У (^к^ хк +1) =
= (хпд - Нхк - Ещ )Т (хпд - Нхк - Ещ).
Можно сделать замену: Z = (х£д -Нхк - Еик), получим ф = Минимум такого функционала находится в начале координат:
2 = 0 < = > х^д - Нхк - Еик = 0.
Выразим ик, использовав псевдообратную матрицу Е = (ЕГЕ)-1 ЕТ:
и*
Е(хТ - Нхк).
(4)
Решим теперь задачу, имеющую ограничения.
гпд _ г
хк хк+1
хп ^Нх^т Еи ^
ф=
Сделаем замену:
Z = хпд - Нхк - Еик» щ = Е (хпд - Нхк - Z).
(5)
Выразим ограничения на 2. Подставим в интервальные ограничения на ик:
*
и < Р(х^ - Их, - Z) < и+.
Выразим отсюда 2:
Ри+ - х^ + Ихк < Z < Риг - х^ + Ихк,
где
р = (РТР)-1 РТ.
Обозначим Z~ = Ри - х^ + Ихк, Z + = Ри" - х^ + Ихк. Функционал будет иметь вид:
р =
Z
= ZTZ.
Таким образом, мы свели исходную задачу вычисления управления к задаче:
z * = а^тт (=zTz
z -<z<z"
(6)
■к , , ^ = ри х 1 называемой задачей оп-
где
Z = хпд - Ихк - Рик» ик = Р(хпд - Ихк - ^ z - = Ри + - хрд + Ихк, z + = Ри - - хрд + Щ
тимизации квадратичного функционала с интервальными ограничения ми. Решив её, мы сможем вычислить оптимальное управление.
У функционала с интервальными ограничениями может быть 2 слу чая условного минимума:
- условный минимум находится внутри интервальных ограничений;
- условный минимум находится на одном из ограничений.
2
Рис. 1. Квадратичный функционал и плоскости, являющиеся границами соответствующих ограничений-неравенств
Рис. 2. Линии уровня и прямоугольник, который задают ограничения
Исходя из этого рисунка видно, что в первом достаточно просто вычислить положение вершины параболоида. В нашем случае это соответствует уравнению:
Z = х? - Нхк - Еик = 0, решение которого является уравнение (4), рассмотренное выше.
Более интересным является случай, когда вершина параболоида находится не внутри области, удовлетворяющей ограничениям задачи (см. рис. 3)
Решить данную задачу можно следующим способом:
- Zг = 0, если удовлетворяет интервальному ограничению;
- Z* = , если 0 < Z Z + и Z = Z+i, если Z < Z+i< 0,
то есть производится проецирование начала координат на ближайшее ограничение.
Рис. 3. Линии уровня, проецирование и интервальные ограничения
Этот подход верен из-за того, что полученная таким образом точка будет находится ближе всего к началу координат, то есть обладает минимальной нормой из всего допустимого множества.
Рассмотрим другой подход к решению задачи стабилизации, а именно погружение модели в задачу математического программирования следующим образом.
Пусть имеется объект управления:
хк+1 = Ихк + Рик.
Перенесём слагаемое из правой части в левую и проделаем замены:
хк+1- Рик = Ихк = А2 = Ъ
где А = (Е -Р), Ь = И2к,
2к — вектор расширенного пространства переменных. Тогда задача будет иметь вид:
Z* = агтп{( = ZTZ\Z-г <2г <Z+г, AZ = Ъ);
Л2 = Ь — линейное многообразие, Z г < Д < Z + — интервальные
ограничения, многомерный параллелепипед, ф = 2[2 — минимизируемая функция, по сути, её минимизация сводится к поиску ближайшей к началу координат точке, которая будет лежать на линейном многообразии и удовлетворять интервальным ограничениям.
Предлагаемый приближённый алгоритм схематически может быть изображён на рисунке 4.
Рис. 4. Схема метода
Алгоритм заключается в последовательных проекциях текущего приближения (можно начать с начала координат) на линейное многообразие и интервальные ограничения. Однако можно заметить, что если остановиться на линейном многообразии, то решение не будет удовлетворять интервальным ограничениям, а если остановиться на интервальных ограничениях, то решение не будет удовлетворять линейному мно-
гообразию. Поэтому предлагается ввести точность £, создав новые ограничения следующим образом:
Z-. + *. <^ <Z+..
Теперь стоит делать проекции на новые ограничения и после проекций на линейное многообразие стоит проверять удовлетворение изначальных интервальных ограничений. Таким образом, получаем следующий алгоритм:
1) Сужение области Z . + < ^ < Z+■ -
2) Проекция начала координат на линейное многообразие
^ = ( АтА) 1 АтЪ.
3) Если точка удовлетворяет ограничениям, то победа.
4) Проекция точки на суженные интервальные ограничения (см. ниже).
5) Проекция точки на линейное многообразие
^ = ( АтА)-1Ат (Ъ + AZг_l) + Zг_l = ( АтА)- 1АтЪ +
6) Если точка удовлетворяет исходным ограничениям, то победа, иначе возвращение к шагу 4.
Выполнить 4-й пункт вышеупомянутого алгоритма можно по аналогии описанного ранее метода: проходим по всем интервальным ограничениям и если они не удовлетворяются, то меняем координату точки на ближайшую, удовлетворяющую ограничениям.
Пример:
Интервальные ограничения 3 < х < 5, -1 < у < -4.
Точка 2г-1 = (0; -2). Смотрим на первую координату — «0». Не удовлетворяет. Ближайшее ограничение к нулю — «3». Значит меняем координату на «3». Вторая координата «-2». Удовлетворяет ограничениям, значит не меняем координату.
Получили точку = (3; -2).
Визуально представить данный алгоритм можно следующим образом:
Рис. 5. Визуализация приближенного метода
Таким образом, возможно получить решение, удовлетворяющее всем ограничениям с точностью е.
Заключение
В данной работе были представлены два метода стабилизации программных движений — точный и численный, основанный на последовательных проекциях на линейное многообразие и на интервальные ограничения.
Список литературы
1. Козлов В.Н. Проекционный метод синтеза ограниченных оптимальных управлений динамических систем энергетики: монография. - СПб.: ПОЛИТЕХ-ПРЕСС, 2019. - DOI: 10.18720/SPBPU/2/i19-277.
2. Козлов В. Н. Системный анализ, оптимизация и принятие решений: учеб. пособие / В.Н. Козлов. - СПб.: Изд-во Политехи. ун-та, 2011. - 244 с.
3. Потапов А.П. Линейная алгебра и аналитическая геометрия. Теория, задачи и упражнения: учебное пособие. - СПб.: Изд-во Политехн. ун-та, 2013. -DOI: 10.18720^РВРи/2^20-146.
4. Козлов В.Н., Куприянов В.Е., Шашихин В.Н. Теория автоматического управления: учебное пособие. - СПб.: Изд-во Политехн. ун-та, 2008.
УДК 519.217.8
ёо1:10.18720/БРВРШМ24-197
Овечкин Николай Николаевич,
студент
ВЫВОД МОДЕЛИ СМО С ПОСТОЯННЫМИ КОЭФФИЦИЕНТАМИ
Россия, Санкт-Петербург, Санкт-Петербургский политехнический университет Петра Великого,
оуесИкт.01 @тЬох.ги
Аннотация. В статье рассматривается вывод системы массового обслуживания с постоянными коэффициентами. Вывод строится на предположении, что процесс работы системы массового обслуживания можно рассматривать как марковский процесс, а поток событий ординарный и стационарный. Поток событий является потоком события без последействия. Также поток является простейшим и стационарным. Результатом проделанной работы является модель, описывающая порт как систему массового обслуживания с постоянными коэффициентами. Модель порта представлена системой массового обслуживания с ограниченной очередью.
Ключевые слова: система массового обслуживания, порт, закон Пуассона, марковский процесс, очередь, поток событий.