CC BY
27
УДК 519.7: 658.512.2
А. Б. Бушуев, О. К. Мансурова
СИНТЕЗ УПРАВЛЕНИЙ В ДВУМЕРНОЙ СИСТЕМЕ С ГЛАДКИМИ НЕЛИНЕЙНОСТЯМИ*
Рассматривается объект управления — нелинейный компенсационный гомео-стат как математическая модель конфликтных ситуаций между двумя каналами. Для синтеза управлений, разрешающих конфликт, используются потенциальные функции канонических катастроф.
Ключевые слова: канонические катастрофы, гиперболическая омбилика, компенсационный гомеостат.
Рассмотрим класс двухканальных систем управления, каждый из каналов которых задается дифференциальным уравнением первого порядка с правой частью вида
^ (х, у, а) = , (1)
п т г
где аг — коэффициенты (/=0,1,2,...), х=х(0 и у=у(0 — координаты системы, хеЯ1, уеЯ1 (п=0,1,2,...; т=0,1,2,...), I — время.
Целью работы является синтез алгоритма управления, обеспечивающего разрешение конфликта между координатами х(() и у(^) для типовых ситуаций: конкуренции, нейтралитета, союзничества, конфронтации и т.п. [1]. Управление формируется в виде нелинейной функции от координат, в качестве которой выбирается потенциальная функция одной из канонических катастроф Р. Тома [2].
Известны примеры использования канонических катастроф для синтеза управления в нелинейном компенсационном гомеостате [3], математическая модель которого задается уравнением (1). Управление синтезируется на основе несимметричной катастрофы типа „эллиптическая омбилика". Сигнал управления складывается из двух составляющих: „мягкую" часть управления вырабатывает командный генератор, а „жесткая" часть формируется в виде обратной связи от координат управляемого гомеостата.
Возможно использовать симметричную каноническую катастрофу типа „гиперболическая омбилика" при синтезе управления для линейных стационарных систем, матрица объекта управления которых представляется в нормальной жордановой форме [4]. Для каждого собственного движения, задаваемого клеткой жордановой матрицы, формируется замкнутое управление в виде потенциальной функции канонической катастрофы типа „гиперболическая омбилика". В качестве координат потенциальной функции выбираются переменные состояния соответствующей клетки жордановой матрицы. Коэффициенты регулятора выбираются из условия обеспечения асимптотической устойчивости замкнутой системы. Как отмечается в работе [4], такой закон управления обладает робастной устойчивостью в случае неопределенности задания как структуры, так и параметров объекта управления.
Для моделирования типовых ситуаций взаимной работы каналов используются градиентные системы, у которых антиградиент потенциальной функции системы направлен по вектору скорости координат, т. е.
Работа выполнена при финансовой поддержке ФЦП „Научные и научно-педагогические кадры инновационной России" на 2009—2013 гг. (государственный контракт № 16.740.11.0553).
Синтез управлений в двумерной системе с гладкими нелинейностями 39
дН (х, y, а) = _ • дН (х, у, а) =_ у дх ду
где Н (х, у, а) — потенциальная функция системы (1). Коэффициент пропорциональности между правой и левой частями в уравнениях (2) для упрощения принят единичным. Для градиентных систем переходные процессы по координатам х=х(7) и у=у(0 при выборе численных значений коэффициентов а\ системы, обеспечивающих условие устойчивости при задании антисимметричных начальных условий х(0)= - у(0), получаются плавными, одна из координат располагается в первом, а другая — в четвертом квадранте плоскости координат. Та или иная типовая ситуация обеспечивается выбором знаков коэффициентов аг обратных и перекрестных связей в каналах системы (1). При п=т=0 в (1) в качестве слагаемых используются коэффициенты а0, задающие вынужденное движение и рассматриваемые как управление, подаваемое на систему (1), например, для канала с координатой х(^):
'х = ЪЪЪа>хпут + ао. (3)
п+1 т+1 г+1
Выберем в качестве управления потенциальную функцию одной из канонических катастроф и подадим его на систему (3). В качестве координат х и у канонической катастрофы выбираем координаты системы (1) х=х(7) и у=у(0. Тогда уравнение правой части замкнутой системы в каждом канале будет иметь вид
х = ЪЪЪагхпут + аоУ (х, у, X), (4)
п+1 т+1 г+1
где V(х, у, X) — потенциальная функция канонической катастрофы, X — вектор управляющих параметров катастрофы X/ (/=1, 2, ..., к), к — число управляющих параметров катастрофы. Потенциальные функции некоторых канонических катастроф приведены в таблице.
к Потенциальная функция V (х, у, X) Катастрофа
1 3 х --+ X, х 3 1 Складка
2 4 2 х х --+ X1--+ X 2 х 4 2 Сборка
3 3 3 х + у + X1 ху + X 2 х + X 3 у Гиперболическая омбилика
3 х13 _ 3 ху2 + X1 (х2 + у2) + X 2 х + X 3 у Эллиптическая омбилика
Приравняв выражение (4) к нулю, получим нелинейные алгебраические уравнения для определения стационарных состояний системы
ЪЪЪа>хпут + V (х, у, X) = 0,
п+1 т+1 г+1
связывающих параметры аг (4) с управляющими параметрами канонической катастрофы.
Для синтеза регулятора выбираем значения управляющих параметров X, обеспечивающие устойчивость заданного конечного состояния равновесия.
Рассмотрим пример синтеза управлений для компенсационного гомеостата, моделирующего стереотип конкуренции антагонистов [5]. Уравнения антагонистов задаются системой (1):
2 2 х = «2 х + «1 у + а0, у = а.5 у + а4 х + а3.
(5)
Здесь первое слагаемое в правой части первого уравнения получается при /=2, п=2, т=0; второе — при /=1, п=0, т=1; третье — при /=0, п=0, т=0. Для второго уравнения (5) первое слагаемое получается при /=5, п=0, т=2; второе — при /=4, п=1, т=0; третье — при /=3, п=0, т=0.
Сначала рассмотрим собственное движение антагонистов при а0=а3 =0 из начальных условий х(0) = -у(0) = -0,1. Для обеспечения устойчивости заданных состояний равновесия х(ю) = -у(да) = -1 выбираем значения а1=а4 = -3, а2= -а5 =3. График переходного процесса собственного движения приведен на рис. 1, а.
Как видно из рисунка, обе координаты строго антисимметричны, и явление конкуренции не проявляется. При подаче на вход канала у(0 ресурса питания конкурентов в виде ступенчатого воздействия (а3 =1(0) проявляется явление конкуренции: ресурс распределяется неравномерно, поскольку х(ю) = -1,09, у(ю) = 1,19 (рис. 1, б).
а)
х, у 1
0,5 0
-0,5 -1
б)
у(^)
х(0
х, у 1,5
1
0,5 0 -0,5 -1 -1,5
у®
.....1,19
х(0 -1,09
0
/, с Рис. 1
0
А с
Синтезируем регулятор, обеспечивающий разрешение конфликта между конкурентами, т.е. получим управление, которое переводит систему (5) из конкурентного х = -1,09, у = 1,19 в устойчивое стационарное состояние х=у=0. В качестве управления выберем потенциальную функцию катастрофы типа „сборка". Подставив V(х,у, К) в (5), получим
х = а2 х +а1 у + а0
( 4 2 Л
х л х -
--+ К--+ К х
V 4 1 2 2 у
у = а5 у +а4 х + аз
( 4 2
— + К3 — + К4 у
V 4 3 2 4
Найдем стационарные состояния из системы уравнений
„2 >
0 = а2 х +а1 у + а(
( 4 х л х л
--ьК--+ К х
V 4 1 2 2 у
( 4 2
0 = а5у2 +а4х + а3 + К3 + К4у V 4 2
(6)
Выберем следующие соотношения К = —
2а2
К2 =-
а
0
а1 у а0 х
К3=--
2а5
К 4 = —
ал х
а
3
а3 у
(7) Тогда
4 4
ху
получим систему 0 = а0 —, 0 = а^ , в которой имеются четыре пары нулевых корней х=0 и
у=0. Знак управляющего сигнала можно задать знаком коэффициентов а0 и а3 исходя из условий обеспечения устойчивости стационарных состояний х=у=0 при заданных начальных условиях х=-1,09 и у =1,19.
1
2
1
2
Синтез управлений в двумерной системе с гладкими нелинейностями 41
График затухания конфликта приведен на рис. 2 для а0=32 и а3 =-32. Как видно, графики конфликтующих сторон асимптотически сходятся к нулю.
X, V 1
0,5
0
-0,5 -1 -1,5
0 1 2 3 с
Рис. 2
Выбор канонической катастрофы для синтеза управления зависит от вида нелинейной функции правой части уравнения (1) гомеостата. Необходимо выбирать потенциальную функцию катастрофы, старшая степень координаты которой хотя бы на порядок выше старшей степени нелинейной функции (1). Для разрешения конфликта назначением управляющих параметров X необходимо компенсировать первые степени координат состояния системы х и у и обеспечить условие устойчивости стационарных состояний.
1,19 y(t)
-1,09 x(t) - -
список литературы
1. Горский Ю. М. Основы гомеостатики (Гармония и дисгармония живых, природных, социальных и искусственных систем). Иркутск: Изд-во ИГЭА, 1998. 337 с.
2. Постон Т., Стюарт И. Теория катастроф и ее приложения. М.: Мир, 1980. 607 с.
3. Бушуев А. Б. Математическое моделирование конфликтов в техническом творчестве // Науч.-техн. вестн. СПбГУ ИТМО. Программирование, управление и информационные технологии. 2005. Вып. 19. С. 26—32.
4. Жуматаева Ж. Е., Бейсенби М. А. Построение систем управления в классе трехпараметрических структурно-устойчивых отображений // Вестн. Казахского национального технического университета. 2010.
5. Бушуев А . Б. Математическое моделирование процессов технического творчества. СПб: СПбГУ ИТМО, 2010. 306 с.
Сведения об авторах
Александр Борисович Бушуев — канд. техн. наук, доцент; Санкт-Петербургский национальный иссле-
довательский университет информационных технологий, механики и оптики, кафедра систем управления и информатики Ольга Карибековна Мансурова — канд. техн. наук; Национальный минерально-сырьевой университет
„Горный", Санкт-Петербург; доцент; E-mail: [email protected]
Рекомендована кафедрой Поступила в редакцию
систем управления и информатики 13.12.12 г.
