УДК 62-50
В.Д. ДМИТРИЕНКО, д-р техн. наук, проф. НТУ "ХПИ" (г. Харьков), В.И. НОСКОВ, канд. техн. наук, гл. инженер ООО "Преобразователь"
(г. Запорожье),
Н.В. МЕЗЕНЦЕВ, ст. преподаватель НТУ "ХПИ" (г. Харьков)
СИНТЕЗ РЕГУЛЯТОРОВ МЕТОДОМ АКОР А.А. КРАСОВСКОГО ПРИ НЕЛИНЕЙНО ВХОДЯЩИХ УПРАВЛЕНИЯХ И СЛУЧАЙНЫХ ВОЗМУЩАЮЩИХ ВОЗДЕЙСТВИЯХ
Предлагается новая модификация метода аналитического конструирования регуляторов по критерию обобщенной работы А.А. Красовского для случая нелинейного вхождения управлений в систему уравнений, описывающую объект, а также при наличии случайных возмущающих воздействий.
Ключевые слова: метод аналитического конструирования регуляторов по критерию обобщенной работы, случайные возмущающие воздействия.
Постановка проблемы и анализ литературы. Метод аналитического конструирования регуляторов по критерию обобщенной работы (АКОР), разработанный А.А. Красовским, позволяет выполнять синтез регуляторов для объектов, описываемых системами обыкновенных нелинейных дифференциальных уравнений высокого порядка [1 - 3], чем он выгодно отличается от многих других методов синтеза оптимальных регуляторов [4 -7]. При этом математическая модель объекта управления в самом общем случае имеет вид:
дх,
dt
+ ft (*i,
’t) = Z Фу (xi’---’ xn > t)uj +
j=i
+ Z^k (X1>---> Xn >t )^(t X k=1
(1)
где ( = 1, п) - фазовые координаты объекта управления; /, Ф., пк (/ = 1, п, ] = 1, т , к = 1, г) - непрерывные функции, описывающие объект управления; £,к (/) - белые шумы (последовательности статистически независимых
5-импульсов, случайных по площади и разделенные сколь угодно малыми, но конечными промежутками времени); Uj (] = 1, т) - управления, определяемые из условия минимума функционала
J = M[Уз(xi(t2),Xn(t2), t2)] + M
jQ(X1,... xn,t ')dt
x
n
+
Здесь М - символ математического ожидания; Уз, Q - заданные непрерывные функции; р, д - положительные четные числа, удовлетворяющие условию 1/ р +1/ д = 1; kj (] = 1, т) - положительные коэффициенты; V - функция, являющаяся решением уравнения
Управления и j, определяемые из условия минимума функционала обобщенной работы (2), являются оптимальными как для случая объекта (1), так и для случая детерминированного объекта (£,к (/) = 0, к = 1, г). Это свойство определяемых управлений и возможность синтеза регуляторов для объектов, описываемых системами нелинейных обыкновенных дифференциальных уравнений высокого порядка, является несомненным достоинством метода аналитического конструирования регуляторов по критерию обобщенной работы. Однако линейное вхождение управлений в систему дифференциальных уравнений объекта (1) является определенным ограничением для применения рассматриваемого метода аналитического конструирования регуляторов. Известные способы преобразования математических моделей с нелинейно входящими управлениями к моделям с линейно входящими управлениями не всегда помогают, поскольку такие преобразования часто ухудшают качество синтезируемых систем управления. В связи с этим были разработаны методы синтеза регуляторов по критерию обобщенной работы для различных математических описаний объектов с нелинейно входящими управлениями [8 - 11]. Однако во всех указанных работах синтез регуляторов рассматривался в условиях отсутствия случайных возмущений.
Целью статьи является расширение области применения метода аналитического конструирования регуляторов по критерию обобщенной работы на объекты, описываемые системами нелинейных обыкновенных дифференциальных уравнений с нелинейно входящими управлениями и при наличии случайных возмущающих воздействий.
Исследования и результаты. В работе [11] доказана одна из наиболее общих теорем аналитического конструирования регуляторов по критерию
(3)
при граничном условии
V(/ = Г2) = Vз .
(4)
обобщенной работы с нелинейно входящими управлениями для случая
р = д = 2.
Теорема 1. Пусть объект описывается системой уравнений
Жх
Ж
+ / , (хЪ • • •, хп, О = X ф j■ (x1, • • •, Хп, 1)у\(и1 j■)^12 (u2ij )
и=1
тогда оптимальными в смысле минимума функционала
(
2 т п 2
1 = V [ Х1 (^2 X •, Хп О2Х ^2] + | ^ + ХХ|
^ j=l , =1 11
2 2
и1у и2,
к 2 к 2
1и 2и
dt
(5)
(6)
являются управления
1 2 (ЗV 1 2 //-\
и1 и = -Ку ф и ¥,>■ /2и1 и;
и2и = -к2,у ^ фи V■ V,* /2и2и =
(7)
где V- решение уравнения (3) при граничном условии (4); Уи, (,= 1,п,
и = 1, т) - непрерывные функции от управлений, соответственно ии - и и2у-.
Обобщим эту теорему на случай, когда объект управления будет описываться системой уравнений
ЖХ;
dt
+ /, (хъ •, хп, t) = X фи(^ •, хп, (и1 и )¥12 (и2и) +
и=1
(8)
+£п *(xl, •, хп, t )^к(t).
к=1
Теорема 2. Пусть объект описывается системой уравнений (8), тогда оптимальными в смысле минимума функционала
I = М[Vз(Хl(t2), ..., Хп&), t2)]+ М
+ М
т п t2
XX XI
и=1 ,=1 ^
^2 2 ^
и1ги и2и
к2 к2
1и 2и
dt
(9)
являются управления (7), в которых К есть решение уравнения (3) при граничном условии (4).
+
Доказательство.
Полная производная функции V в силу уравнений объекта (8) и уравнения в частных производных (3) равна
ау дУ -П дУ ах, Л дУ ^ 12^ ^ /■
+ X—— = — +^— X Фу Ъ Ъу + X^\к (') - /,
к=1
а д' дх а д' дх,
,=1 , ,=1 ,
(10)
_ ^ дУ ^ 12 Vі дУхр в ^
=-(ї+X д- X Ф-ъ +X дт X ^ ^(').
і=1 , -=1 і=1 , к=1
Интегрируя выражение (10) в интервале времени ['і, '2], получим
ь
У[х1(і2), х„(/2),/2]-V[х1(/1), х„(ґ1),ґ1 ] = -|<2& +
■XXJ:дУ Ф»ъ- —+ZX.ffУ^ « ('Х*.
і=1 к=1 '1 і
(11)
і=1 -=1 '1 ^
Учитывая, что по условию V [х1 (ї2 ), хп ^ І2 ]= V, [х1(^Х Хп ({2І к ],
из (11) можно получить
М [V, (х^), Хп &)Л )]+ М
12
I Qdt
= М[V(Х1Й), Хп)] +
+ М
п т t2 ;ат/
XX X Ь- ф, ъ!,- , і=1 ,=111 СХ!
+М
п г Ч ят/
X X I дЧ* § * «л
1=1 к=1 ?!дх
(12)
Покажем, используя методику А.А. Красовского [1], что последний член в выражении (12) не зависит от управлений (7). Учитывая принятую модель белых шумов в уравнениях (1) и (8), можно записать
^ (') =X Вка 8('-%а),
(13)
где ^к (t) - обобщенный стационарный случайный процесс; ВкЖ -
независимые случайные центрированные величины; т1, т2, ... - случайные моменты времени, которым соответствуют "импульсы" 5-функции.
Согласно соотношению (13) последний член выражения (12) можно преобразовать к виду
Я = М
XXX
і=1 у=1 а
с дV
Вка | т- Пі*5^ - та
V т а-
(14)
- =1
а
Вследствие воздействия 5-функций 8(t — Tj) (d = const) фазовые координаты объекта управления (8) за интервал времени от Td — 0 до Td + 0 получают приращения
Ax,d = \^,kBkd8(t - Td)dt, i =1, n.
хій - і 'Іік _0
С точностью до бесконечно малых более высокого порядка приращения Axid
не зависят от управлений uU]■ и и2у ( і -1, л, ] -1, от ). Таким образом,
изменения фазовых координат, вызванные воздействием возмущений в виде 5-функций, с точностью до бесконечно малых величин более высокого порядка не зависят от управлений и и2і]-, а это значит, что и возмущения
д V ___ ____
функций У\хъ хп, /], —, цік[х,хп,{] ( і -1, п, к -1,т), вызванные
дхі
воздействием возмущений в виде 5-функций с точностью до бесконечно малых более высокого порядка не зависят от управлений и и2і]- ( і -1, п,
] -1, т). Поэтому и Я с точностью до бесконечно малых величин не зависит от управлений.
Функционал (9) с учетом выражений (12) и (14) преобразуем к виду
I = M[V(Xl(tl), Xn(tl), tl)]+ M
n m t2 p)j/
EE h-ф- ^ j
i=1 j=1 tj OXi
+M
m n t2
E E I
j=1 i=1 tj
22 Uj U-
к 2 к 2
lj
2j
dt
(І5)
+ R .
Поскольку выполняются следующие соотношения:
EE
i=1 j=l
U^ + ки OV Фі/Vi/Vjj I2uij
Kj J Oxi J J J j
=EE
i=1 J=1
2
n m U, .. n m О—
2
- IV .2 lj
U1j O— 1 2 U1J
ё + IT + W 1/
-t? O—
= 2EE-rf+EE^- Ф-/
i=l j=l к1/ ,=l j=l OX,
1.2 .
EE
i=1 J=1
7— + к2,/|— Ф-УЖ I2u
OV
к
2j
Or,
2ij
=EE
i=1 j=1
22 U 2ij OV 1 2 U 2/j
к-+^ф j+кг
2ij i 2ij
+
2
2
п т М~ п т р>Т/
—21 +1 IIV фу у/ у/
,—1 /—1 к2/ ,—1 ./=1 Сл,-
то функционал (15) можно записать в виде
I — МV(хД/Д ..., Хп(/Д /)] +
+М
1 /2 пт
2121
2 /1 , =1 / =1
^М1/ , СК ! 2
*“+к1/^“ ф/У/У// 1/
Л2
а
+ М
1 /2 п т
1111
2 ^ ,-1 ,-1
' /\ ,—1 /=1
Если управления и1гу и и2у
и2/ , СК 1 2 л
+ к2„^ ф/V/ V/ и2,/
V кЪ/ / Сх, /
2
а
определяются соотношениями (7), то
подинтегральные выражения в функционале (16) обращаются в нуль и он принимает минимальное значение. Таким образом, управления (7) действительно являются оптимальными в смысле минимума функционала (9) и, следовательно, теорема 2 доказана.
Заметим, что вид оптимальных управлений (7) при случайных возмущениях описанного типа одинаков для детерминированных и стохастических систем. При этом величины управлений как и в наиболее общих теоремах по аналитическому конструированию регуляторов по критерию обобщенной работы А.А. Красовского [1] не зависят от уровня шумов.
Из общих уравнений (5), (8) объекта управления и его управлений (7) могут быть получены различные частные системы уравнений, важные для различных областей применения. В качестве примера рассмотрим задачу синтеза системы управления асинхронным приводом, математическая модель которого в осях и и V может быть записана следующим образом:
—1 + апУ1 + Я1зУз = А 008 га/; Л
2 + a--у- + а24У4 = А 8т га/;
“У
Оу3
Л
ау4
Л
+ я31У! + а33у3 + а345У40 = 0;
+ а42У2 + 0^44У4 + 0^35У3^ — 0;
ао р г
“Г - Т (Я5о(У2У3 - У1У4 ))“Мс (0 ) - I Пк4к (/) = 0 а J к—1
+
где у, у - потокосцепления по оси и; t - время; а 1,а1з,Я22, •",Я50 -постоянные коэффициенты; А - амплитуда (первое управление) питающего напряжения частоты га (второе управление); у, у - потокосцепления по оси
V; ^ - угловая скорость вращения ротора электродвигателя; Р - число пар полюсов двигателя; 3 - момент инерции двигателя и механизма, приведенный к валу двигателя; Мс - момент сопротивления нагрузки; цк - характеристики объекта (постоянные коэффициенты или непрерывные функции); ^ (/) -белые шумы.
Используя соотношения (7) и систему уравнений объекта (17), несложно получить выражения для управляющих воздействий
т.е. имеем известный закон при частотном управлении: отношение амплитуды сигнала к его частоте равно константе.
Выводы. Таким образом, доказательство теоремы 2 расширило область возможного применения метода АКОР на объекты, описываемые системами обыкновенных дифференциальных уравнений вида (8), т.е. на объекты с нелинейно входящими управлениями и при наличии случайных возмущающих воздействий. В дальнейшем предполагается обобщение полученных результатов и разработка численных методов для синтеза конкретных систем управления.
Список литературы: 1. Красовский А.А. Системы автоматического управления полетом и их аналитическое конструирование. - М.: Наука, 1973. - 560 с. 2. Красовский А.А., Буков В.Н., Шендрик В.С. Универсальные алгоритмы оптимального управления непрерывными процессами. -М.: Наука, 1977. - 272 с. 3. Справочник по теории автоматического управления / Под ред. Красовского А.А. - М.: Наука, 1987. - 712 с. 4. Мирошник И.В. Теория автоматического управления. Нелинейные и оптимальные системы. - СПб.: Питер, 2006. - 272 с. 5. Дорф Р., Бишоп Р. Современные системы управления. - М.: Лаборатория Базовых Знаний, 2004. - 832 с. 6. Методы классической и современной теории автоматического управления. Учебник в 5-ти томах. Т.4: Теория оптимизации систем автоматического управления / Под ред. К.А. Пупкова и ЕД. Егунова. - М.: Изд-во МГ ТУ им. Н.Э. Баумана, 2004. - 744 с. 7. Нестационарные системы автоматического управления: анализ, синтез, оптимизация / Под ред. К.А. Пупкова и ЕД. Егунова. - М.: Изд-во МГТУ им. Н.Э. Баумана, 2007. - 632 с. 8. Эволюционные методы компьютерного моделирования / Верлань А.Ф., Дмитриенко В.Д., Корсунов Н.И.,
Из соотношений (18), (19) находим
A k
— = —- = const, ю k2
ШорохB.A. - К.: Наук. думка, 1992. - 25б с. 9. НосковВ.И., Д^тр/пенко В.Д., Зaпoлoвeкuй Н.И., Леонов С.Ю. Моделирование и оптимизация систем управления и контроля локомотивов. -Харьков: ХФИ "Транспорт Украины", 2003. - 248 с. 10. ДaнuленкoA^., Дттрыенко В.Д., Зaпoлoвeкuй Н.И. Математические модели оптимальных систем управления тяговым асинхронным приводом тепловозов // Электронное моделирование. - 1991. - Т. 13. - № 2. - С. 40 - 44. 11. Дмumpuенкo В.Д., Носков В.И. Мезенцев Н.В. Оптимизация функционала обобщенной работы нри нелинейно входящих управлениях. Праці Луганського відділення Міжнародної Академії Автоматизації. - Луганськ: Луганське відділення Міжнародної Академії Автоматизації. - 2005. №° 1. - С. 17 - 22.
УДК б2-50
Синтез регуляторів методом АКУР О.А. Красовського при управліннях, що нелінійно входять, та випадкових впливах, що збурюють / Дмитрієнко В.Д., Носков В.І., Мезенцев М.В.
// Вісник НТУ "ХПІ". Тематичний випуск: Інформатика і моделювання. - Харків: НТУ "ХПІ". -2009. - № 13. - С. 53 - 60.
Пропонується нова модифікація методу аналітичного конструювання регуляторів за критерієм узагальненої роботи О.А. Красовського для випадку нелінійного входження управлінь у систему диференціальних рівнянь, що описують об'єкт, а також при наявності випадкових впливів, що збурюють. Бібліогр.: 11 назв.
Ключові слова: метод аналітичного конструювання за критерієм узагальненої роботи, випадкові впливи, що збурюють
UDC 62-50
Synthesis of regulators method of AKOR of A.A. Krasovski at nonlinear incoming controls and casual revolting influences / Dmitrienko B.D., Noskov V.I., Mezentsev N.V. // Herald of the National Technical University "KhPI". Subject issue: Information Science and Modelling. - Kharkov: NTU "KhPI". -2009. - № 13. - P. 53 - 60.
New modification of method of the analytical constructing of regulators is offered on the criterion of the generalized work of A.A. Krasovski for the case of the nonlinear including of controls in the system of equalizations, describing an object, and also at presence of casual revolting influences. Refs.: 11 titles.
Key words: method of the analytical constructing of regulators on the criterion of the generalized work, сasual revolting influences.
noemynrn.a в pедaкцuю 25.05.2009