CC BY
12
УПРАВЛЕНИЕ, МОДЕЛИРОВАНИЕ, АВТОМАТИЗАЦИЯ
DOI: 10.24143/2072-9502-2019-4-7-17 УДК 519.622.2
АНАЛИЗ МИНИМАЛЬНОГО ВРЕМЕНИ СТАБИЛИЗАЦИИ КОЛЕБАТЕЛЬНОГО ПРОЦЕССА ПРИ РАЗЛИЧНЫХ РЕЖИМАХ ЕГО ПРОТЕКАНИЯ
Д. А. Асадова
Институт систем управления Национальной Академии Наук Азербайджана,
Баку, Азербайджан
Представлен численный подход к решению задач оптимального управления процессами установления колебательных систем посредством краевых управлений. В отличие от обзорно рассмотренных работ других авторов по схожей тематике, задача установления колебательного процесса численно исследуется в некоторой временной окрестности требуемого режима, а в математической модели колебательной системы учитывается наличие силы трения. Приводятся схемы вывода основных формул, а также алгоритмы численного их решения методами оптимизации первого порядка. Получены соответствующие аналитические формулы для компонент градиента целевого функционала по управляющим воздействиям. Проводится сравнение минимального времени установления колебательного процесса как при оптимальном управлении им, так и при самоустановлении процесса. Приведены результаты численных экспериментов. Полученные результаты могут быть отнесены ко многим процессам, описываемым гиперболическими уравнениями.
Ключевые слова: колебательные системы, краевые управляющие воздействия, градиент целевого функционала, минимальное время установления колебательного процесса.
Для цитирования: Асадова Д. А. Анализ минимального времени стабилизации колебательного процесса при различных режимах его протекания // Вестник Астраханского государственного технического университета. Серия: Управление, вычислительная техника и информатика. 2019. № 4. С. 7-17. DOI: 10.24143/2072-9502-2019-4-7-17.
Введение
Задачи управления колебательными процессами всегда привлекали внимание исследователей [1-16]. Большой вклад в исследование задач управления граничными условиями для классических гиперболических уравнений второго порядка был внесен В. А. Ильиным и его учениками [2-7]. В их работах установлены необходимые и достаточные условия существования управлений в зависимости от того, в каком соотношении находятся длина колебательной системы I и финальный момент времени T .
В работах В. А. Ильина и Е. И. Моисеева [3], А. А. Никитина [4] получены аналитические решения задач граничного управления для простейшего одномерного волнового уравнения в классах сильных обобщенных решений для произвольного момента времени Г > l. В работе [7] Г. Д. Чабакаури аналитически решена задача о наилучшем приближении к заданной цели.
В классах слабых обобщенных решений Л. Н. Знаменской были получены аналитические выражения граничных управлений для критического момента Г = I в работах [5, 6], а для моментов времени Г > I - в работе [13] авторами М. Gugat, G. Leugering и др.
В работе [17] Ф. П. Васильевым и Р. П. Ивановым был предложен общий метод приближенного решения задач быстродействия для весьма широкого класса линейных управляемых процессов в банаховых пространствах.
В работе [12] предложен численный метод приближенного построения оптимального управления процессом колебаний однородной струны за счет управления на концах. Вычислительные эксперименты проведены для различных режимов колебаний правого конца струны в случае, когда управляющие функции аппроксимируются полиномами Лежандра.
Задачи точечного управления для гиперболического уравнения были рассмотрены в [1, 8]. В [1] изучены задачи финитного управления. В работе Г. Ф. Кулиева [8] получены необходимые и достаточные условия оптимальности в виде вариационного неравенства для задачи точечного управления гиперболическим уравнением общего вида.
В работе [11] для управляемых процессов, описываемых гиперболическим уравнением, были получены конструктивные оценки в классах сильных и слабых обобщенных решений, позволяющие с помощью вариационного метода М. М. Потапова находить устойчивые приближенные решения задач граничного управления процессом перевода колебательной системы в требуемое состояние точно к заданному моменту времени Т. Также для общего случая теоретическое исследование задач быстродействия с граничными управлениями начато в работе [15], в которой отмечается актуальность разработки соответствующих эффективных численных методов. В работах [14, 16] предложены численные методы решения некоторых задач оптимального управления для волнового уравнения.
В отличие от этих работ в данной работе задача установления колебательного процесса численно исследуется в некоторой временной окрестности требуемого режима. Кроме этого, в математической модели колебательной системы учитывается наличие силы трения. Поскольку рассматриваемым в качестве математической модели телеграфным уравнением описываются многие колебательные процессы, такие как неустановившееся движение жидкости в трубопроводах, процессы переноса электроэнергии в системах электропередач, колебательные движения струн, стержней, пластин, полученные результаты могут быть отнесены ко многим процессам, описываемым гиперболическими уравнениями.
Постановка задачи
Пусть движение колебательной системы с учетом сопротивления среды описывается следующим уравнением гиперболического типа [18] в безразмерных переменных:
d2y d2y dy ^ „ Л ^ „ /1Ч
--v + а— = 0, 0 <X < 1, t > 0, (1)
dt дх dt
где а - коэффициент сопротивления среды; y = у(х, t) - функция состояния распределенной колебательной системы, т. е. ее отклонение от состояния равновесия.
Допустим, что до момента времени t = 0 колебательная система находилась в исходном стационарном состоянии:
y(х, 0) = х), yt(х, 0) = х) = const 0 < х < 1, t < 0. (2)
Требуется перевести колебательную систему в конечное заранее заданное стационарное состояние:
y(х, T) = ^ (х), yt(х, T) = ^ (х) = const, 0 < х < 1, t > T, (3)
где T - время, после которого установится новый стационарный режим (3).
Здесь ц(.), £,(.), Ц1О, £i(.) - заданные функции.
Отметим, что в реальности добиться полного установления (стабилизации) колебательного процесса к какому-либо определенному моменту времени T практически невозможно (из-за погрешностей в работе измерительной техники), поэтому будем считать, что режим (3) достигнут, если выполняются условия
|у(Т) -цДх)| <81, (4)
|у,(х, Т) х)| < 82, х е (0,1), Т < t < Т + DT,
где ВТ - определенный малый заранее заданный интервал времени после момента Т установления процесса; 81, 32 - заданные положительные величины, определяемые в процентах от значений некоторого установившегося режима (ц1(х), ^(х)).
В связи с этим будем рассматривать уравнение (1) на множестве О = [0 < х < 1] [0< t < Т+ВТ]. Управление процессом перевода колебательной системы осуществляется за счет краевых условий, т. е. изменением функции состояния на обоих концах колебательной системы:
у(0, г) = щц), у(1, г) = игц), t > 0, (5)
где и1, и2 е L2 [0, Т + ВТ].
Замечание 1. Отметим, что т. к. управляющие функции и^г), и2(г) могут быть разрывными, нельзя говорить о решении начально-краевой задачи (1), (2), (5) в классическом смысле [18]. Поэтому везде в дальнейшем решение краевых задач будем понимать в обобщенном смысле [19].
Определение. Под решением краевой задачи (1), (2), (5), соответствующим управлениям щ(г), и2(г) е Ь2 [0, Т + ВТ], будем подразумевать функцию у = у(х, г) е Ь2(П), имеющую следы у(х, -)е Ь2 [0, Т + ВТ] при всех х е [0, 1]; у(% г)е Ь2 [0, 1] при всех г е [0, Т + ВТ], удовлетворяющую уравнению (1) почти всюду в П, а условиям (2) - в смысле равенства соответствующих следов функции у, и интегральному тождеству
1
И у( х, 0(фй (х, г) - Ф хх (х, г) - аФг) dхdt +1 х)Фг (х, 0)<А -
П 0
1 Т+ВТ
- Г(?(г) + ац(г))Ф(г 0)dr - Г( 11 ■ / ни ■! / ■ —
I (5(х) + ар,(х))Ф(х, 0^ - | (и2 (г)Фх (1, г) - и1(г)Фх (0, г))dt = 0,
0
справедливому для всех функций Ф(х, г) е Ж22(П) [2-7] таких, что Ф(0, г) = 0, Ф(1, г) = 0 при всех 0 < г < Т и Ф(х, Т + ВТ) = 0 и Фг(х, Т + ВТ) = 0 при всех 0 < х < 1, и для которых краевые условия и первое начальное условие выполняются в смысле равенства элементов из Ь2[0, Т + ВТ] и Ь2[0, 1] соответственно, а второе начальное условие - в смысле равенства элементов из Ж (0, 1).
Целью управления в данной задаче является наиболее быстрый перевод колебательной системы из состояния (2) в состояние (3) с учетом того факта, что достижение режима (3) возможно в смысле (4). В связи с этим рассматривается целевой функционал, состоящий из двух слагаемых: первого слагаемого - минимизируемого момента времени Т, и второго слагаемого, отражающего факт установления колебательного процесса в смысле (4):
Т+ВТ 1 . .
J(и, Т) = Т + I |{г1 [у(х, г) - ц1(х)]2 + г2 [у,(х, г) - 41(х)]2 }drdt ^ Шп . (6)
Т 0
Величина ВТ - определенный заранее заданный интервал времени после момента Т установления процесса; Г\, г2 - заданные весовые коэффициенты.
Таким образом, требуется найти такие управляющие функции и1 (г), и2 (г) е Ь2 [0, Т + ВТ] и оптимизируемый параметр Т, при которых функционал (6) принимает минимальное значение.
Численный подход к решению задачи
Для решения задачи оптимального управления процессом установления колебательных систем, рассматриваемой в данной работе, предлагается использовать итерационные методы оптимизации первого порядка, основанные на применении аналитических формул градиента целевого функционала по управляющим параметрам.
0
Поставленная задача является задачей оптимального по быстродействию управления распределенной системой при заданных значениях функций состояния в момент окончания процесса Т, который рассматривается как оптимизируемый параметр, и с управлением в краевых условиях. Для решения данной задачи могут быть применены два подхода. Согласно первому подходу можно рассматривать Т как параметр и использовать двухуровневую оптимизацию: на верхнем уровне для определения оптимального времени проведения переходного процесса Т применить какой-либо из методов одномерной оптимизации, а на нижнем уровне при заданных текущих значениях Т для определения
JT = J (иТ, Т) = тт J (и, Т)
и
решать задачу оптимального управления распределенной системой с фиксированным временем. Согласно второму подходу Т рассматривается как компонента управления, и для нахождения ее оптимального значения применяется процедура одновременной совместной оптимизации Т и и(г). В данной работе применен второй подход.
Формулы градиента целевого функционала, выводимые ниже, могут быть использованы при получении необходимых условий оптимальности для рассмотренной задачи.
Формулы для градиента функционала по управляющим воздействиям
Формулы для компонент градиента функционала по краевым управлениям. Для получения формул градиента функционала по управлениям и{(г),I = 1,2,... дадим им приращения
Диг- (0, т. е. рассмотрим допустимые управления иг(г), (г) + (г) еL2[0, Т + DT].
Пусть у(х, г; ц), у(х, г; ui + Дui) - решения краевой задачи (1), (2), (5), соответствующие
этим управлениям. Обозначим управляющий вектор через и(г) = (и1(^), и2(г)), тогда Ди(г) = = (Ди1(/), Ди2 (/)) е ¿22 [0, Т + БТ ].
Из (1), (2), (5) следует, что Ду(х, г) = у(х, г; и + Ди) - у(х, г; и) является решением следующей краевой задачи:
ДУи -ДУхх + аДу = 0, х е (0, 1), г е (0, Т + ДТ) ; (7)
Ду(х, 0) = 0, Дуг(х, 0) = 0, Ду(0, г) = Ди1(г), Ду(1, г) = Ди2(г). (8)
Формула для приращения функционала (6), с учетом (7), (8), примет следующий вид:
Т+БТ 1
Д„У(и, Т) = J(и + Ди, Т) - J(и, Т) = | |г [у(х, г) + Ду(х, г) -цДх)]2 dxdt -
Т 0
Т+БТ 1 Т +БТ 1
- 1 1Г [у(х, г)-ц1(х)]2 dxdt + | | г2 [ у (х, г) + Ду, (х, г)-^1(х)]2 dxdt -
Т 0 Т 0
Т+БТ 1 Т +БТ 1
- 1 1 г2 [у(х, г)-^1(х)]2 dxdt = | 12г1 [у(х, г) -ц1(х)]Ду(х, t)dxdt +
Т 0 Т 0
Т + БТ 1
+ 1 1 2г2 [уг (х, г)-^(х)]Ду (х, г^г +R, (9)
Т 0
Т + БТ 1 Т+БТ 1
где R = г 1 1 (Ду(х, г))2ёхЛ + г2 11 (Дуг(х, г))2ёхЛ .
Т 0 Т 0
Пользуясь результатами работы [19], можно доказать следующую оценку:
N * М\\Ди\Ь[0,Т+ДТг (10)
где М - положительная постоянная, не зависящая от у(х, г) и и(г).
Пусть функция у(х, г) является обобщенным решением следующей вспомогательной сопряженной краевой задачи:
д2 у д у
1 ~ + 2Г1(у(х, г)-ф1Т (А)) - 2Г2уй (х, г), Т * г * Т + БТ ,
(11)
д 2у дг2
+ а-^ + 2Г1( у (х, г) - Ф1Т (х)) - 2г2 у г (х, г), Т * г * Т + БТ, дх дг
д2у + а^, 0 * г<Т;
дх2 дг
у(0, г) = у(1, г) = 0, 0 * г < Т + БТ; (12)
у(х, Т + БТ) = 0, 0 < х < 1; (13)
уг(х, Т + БТ) = ау(х, Т + БТ) - 2г2 [у (х, Т + БТ) - у (х, Т)], 0 * х * 1. (14)
Учитывая (7), (8), (11)—(14) в (9) и применяя формулу интегрирования по частям, получим
Т+БТ 1 1
" Г "ЧБТ)-£,(х)]Ду(х, Т+1
ДuJ (и, Т) =1 12г [у( х, г) - ц1 (х)]Ду( х, г )dxdt +12г2 [уг (х, Т + БТ) - ^ (х)]Ду( х, Т+БТ ^х
Т 0 0
1 Т+ВТ 1
-1 2г2 [уг (х, Т) - ^1(х)]Ду(х, Т ^х - 1 1 2г2 уЙ (х, г)Ду(х, t)dxdt +Я =
0 Т 0
Т+БТ 1 Т 1
= 1 1 (У « - У хх - аУ г) Ду (^ + 11 (У и - У хх - аУ г) Ду (х, -
Т 0 0 0
I Т+ПТ 1
-1 Уг(х, Т + БТ)Ду(х, Т + ПT)dxdt + N = 1 1 (угг - Ухх - ауг)Ду(х, t)dxdt -
0 0 0
II , _ . — . , уг(х, Т + БТ)Ду(х,Т + ПT)dxdt +
00
Т+БТ 1 Т+БТ Т+БТ
+ Г ГУ х, у )Ду„ ( х, г Лйхйг +
1 у г (х, Т + БТ) Ду (х, Т + БТ )dxdt + N =1 у г (х, Т + БТ) Ду (х, Т + БТ )dxdt -
0
+БТ 1 Т+БТ Т+БТ
1 1 у( х, у )Дуг (х, г )йхйг + 1Ди1(г )у х (0, г )dt - 1Ди2(г )у х (1, г -
0 0 0 0 Т+БТ 1 Т+БТ 1
- 1 1 у( х, у)Духх (х, + 1 1 ау( х, г)Дуг (х, -
0 0 0 0 1 Т+БТ 1
-1 у г (х, Т + БТ) Ду( х, Т + БТ )dxdt + N = 11 у( х, у)( ДуЙ (х, г) - Духх (х, г) + аДуг )dxdt +
0 0 0 Т+БТ Т +БТ
+ 1 Ди1(г)ух (0, г)Л - 1 Ди2(г)ух (1, г)Л + N =
0 0 Т +БТ Т+БТ
= 1 Ди1(г)у х (0, г)Л - 1 Ди2(г)у х (1, г)Л + N. (15)
00
Из (15) и оценки остаточного члена (10) следует, что функционал дифференцируем по и(г) на L2[0, Т + БТ] и компоненты его градиента определяются по формулам
gradUlJ(и, Т) = ух (0, г), г е [0, Т + БТ];
grad и, 3 (и,Т) = - у х (1, г), г е [0, Т + БТ ].
Формулы для компонент градиента функционала по моменту времени Т. Для получения формулы градиента функционала по моменту времени Т будем считать, что параметр Т является компонентой управления. Рассмотрим приращение функционала, соответствующее приращению конечного момента времени Т, т. е. ЛТ:
Т+ВТ+ЛТ 1
Л^(и, Т) = J(и, Т + ЛТ)- J(и, Т) = ЛТ + | |г [у(х, г) -^(х)]2 dхdt-
T+ AT 0
T+DT 1 T+DT+AT 1 T+DT 1
- i IГ [y(x, t) - Ri(х)] dxdt + J Jr2 [yt (x, t) - £,1(x)] dxdt - J Jr2 [yt (x, t) - £,1(x)] dxdt =
T 0 T+AT 0 T 0
T+DT 1 T+DT+AT 1
= AT + J Jr1 [y(x, t) - p,1(x)] dxdt + J J r1 [y(x, t) - p,1(x)] dxdt -
T+AT 0 T+DT 0
T+DT 1 T+DT+AT 1 T+DT 1
- J Jr1 [y(x, t) - p,1(x)] dxdt + J J r2 [yt (x, t) - £,1(x)] dxdt - J Jr2 [yt (x, t) - x)] dxdt =
ч S1WJ J J'2L^fv
T 0 T+AT 0 T 0
T +DT 1 T +DT +AT 1
= AT + J J r1 [y(x, t) - x)]2dxdt + J J r1 [y(x, t) - x)]2dxdt-
T +AT 0 T +DT 0
T+AT 1 T + DT 1 Л T + DT 1
(x,t)-u,(x)]2dxdt + f \r[y(x, t)-u,(x4 2
J J r\y(x, t) - x)]2dxdt + J J r\y(x, t) - p,j(x)]2dxdt + J J r2 [y,(x, t) - £,j(x)]2dxdt
T 0 T+AT 0 ! T +AT 0
DT + AT 1 % T + AT 1 T+DT 1
+ J J r2 [y,(x, t) - x)]2dxdt - I J J r2 [y,(x, t) - x)]2dxdt + J J r2 [y,(x, t) - x)]2dxdt
T + DT + AT 1 / T + AT 1 T+DT 1
'2L/ ^ул^^лш -i j j J J '2 LZ i( ^' - ^^
T+DT 0 $ T 0 T + AT 0 !
= AT + AT|Jr{[y(x, T + DT)-x)]2 -[y(x, T)-x)]2}dx +
ио теор. о среднем знач. $ 0 !
+ ATIJ r2 {[y (x, T + DT) - x)f-[y, (x,T) - x)f }dx"
= AT
$0
1
1
0
1
+ J r1 (y(x, T + DT) + y (x, T) - 2^(x))(y(x, T + DT) - y (x, T))dx +
0
+ Jr2(у,(x, T + DT) + у,(x, T) - x))(у,(x, T + DT) - у,(x, T))dx
0
Разделив обе части на ЛТ и перейдя к пределу при ЛТ ^ 0, получим формулу компоненты градиента функционала по Т
1
gradT J(и,Т) = 1 + |г (у(х, Т + ВТ) + у(х, Т) - (х))(у(х, Т + ВТ) - у(х, Т))Сх +
+1Г2(у,(х, Т + ВТ) + У,(х, Т) - 2^1 (х))(у (х, Т + ВТ) - У,(х, Т))СХ
0
Полученные формулы можно использовать, в частности, при применении градиентной процедуры одновременной совместной оптимизации вектора (и, Т), включающей метод Флет-чера - Ривза для нахождения текущих значений управляющей вектор-функции и(г) = (и^г), и2(г)) и метод скорейшего спуска - для значений оптимизируемого параметра Т:
ик+1 (г) = ик (г) - акдк;
сЛ
~дТк
Tk+1 = Tk -ak^-r, к = 1, 2, ...
/
0
где ч = graduJС"1); як = gradмJ("к)- Ра-1; = -||graduJ("к("к-1)||2; а
мерной минимизации, выбираемый из условия
- шаг одно-
ак = а^тт Л | (ик
а>0
J%(ик(г) - адк) Тк - а) .
ик (г )еЦ[0,Г+DT ]
Результаты численных экспериментов
Отметим (о чем уже говорилось во введении), что одним из основных теоретических результатов, полученных А. Г. Бутковским, В. А. Ильиным и его учениками, которые занимались вопросами граничного управления струной, являлся тот факт, что при управлении струной на двух концах минимальное время успокоения равно длине струны, а при управлении на одном конце - удвоенной длине струны.
В данной работе при численном решении задач граничного управления переходными процессами в колебательных системах без учета каких-либо ограничений на управления и фазовые переменные данные факты были подтверждены. При решении задачи в безразмерных переменных (т. е. при 0 < х < 1) было установлено, что минимальное время установления переходного режима равнялось 1 (т. е. длине системы) при управлении на обоих концах колебательной системы, а при управлении на одном конце (на другом конце сразу осуществлялся переход на требуемый режим) оптимальное время составляло 2.
В безразмерных единицах значение скорости функции состояния при исходном режиме было взято равным 1, т. е. ^ = 1, при конечном режиме - равным 2 (т. е. ^ = 2), длина системы считалась равной 1. Управление осуществлялось изменением функции состояния на концах (допустим, струны, или какой-либо другой распределенной колебательной системы). Как видно из графиков, изображенных на рис., функции управления носят импульсный характер.
0,2 0,-41 0,6 0,8 1 1,2 1,4
200 100 0 -100 а -200 -300 -400 -500 -600
0,2 0,4 0,6 0,8 1 1,2 1,4
г
а б
Графики оптимальных управлений при краевом управлении колебательной системой: оптимальное управление и1 на левом конце (а); оптимальное управление и2 на правом конце (б)
Так, функция управления на левом конце совершает резкий скачок вверх до 130 ед., затем вниз до -70, а потом с амплитудой, уменьшающейся по абсолютной величине, совершаются колебания вверх и вниз где-то до Т ~ 0,4, далее небольшие осцилляции продолжаются до полного установления (стабилизации) к моменту времени Т ~ 1. Значения функции управления на правом конце сначала уменьшаются до -500, затем увеличиваются до 70, далее с небольшой амплитудой функция осциллирует до Т ~ 0,2, после чего совершаются незаметные глазу небольшие колебания вплоть до момента полного установления при Т ~ 1.
Естественно, что эти исследования носят чисто теоретический характер, поскольку при управлении реальными колебательными процессами необходимо учитывать ограничения на функции управления. Однако они дают возможность оценить вид оптимальной управляющей функции и выигрыш во времени по сравнению с тем, когда колебательная система устанавливается без оптимального управления ею. Известно, что при наличии силы трения колебательная система, выведенная из состояния равновесия, со временем установится сама по себе, такое
установление системы обычно называют самоустановлением (самостабилизацией). В работе для сравнения были также проведены вычислительные эксперименты по наблюдению за временем самоустановления системы. Колебательная система была выведена из состояния равновесия (при £0 = 1) тем, что на границе (в данном случае, на левом и правом конце) была сразу переведена на конечный установившийся режим, т. е. на концах поддерживалось значение скорости функции состояния £,1 = 2 и проводилось наблюдение за временем установления, т. е. перехода к режиму (3) на всем протяжении системы при = 2.
В табл. приведены значения времени самоустановления колебательной системы при численном исследовании ее с помощью явных и неявных разностных схем (здесь h - шаг по пространственной переменной, т - шаг по временной переменной).
Численные результаты решения задачи самоустановления переходных режимов колебательных систем с использованием различных схем аппроксимации
h т T
Явная схема Неявная схема
0,1 0,05 4,5 3,7
0,08 4,5 3,7
0,09 4,6 3,7
0,1 4,6 3,7
0,05 0,01 3,8 3,7
0,02 3,8 3,7
0,03 4,6 3,7
0,05 4,6 3,7
0,1 - 3,7
0,01 0,005 3,8 3,7
0,008 3,8 3,7
0,009 4,6 3,7
0,1 - 3,7
Из таблицы видно, что при решении краевой задачи с помощью неявной схемы (которая была использована и при расчете краевой задачи в решении задачи оптимального управления) время самоустановления системы составляло Т ~ 3,7. Это в 3,7 раза превышает время установления системы при оптимальном управлении ею.
Все это показывает явные преимущества оптимального управления установлением (стабилизацией) колебательной системы.
Заключение
В работе было исследовано время установления (стабилизации) колебательной системы, совершающей переход от исходного установившегося состояния к конечному установившемуся состоянию, как при оптимальном управлении процессом установления, так и при самоустановлении системы; проведено сравнение полученных результатов. При решении задачи оптимального управления процессом установления колебательной системы предложен численный метод, в котором используются итерационные методы оптимизации первого порядка, основанные на применении аналитических формул градиента целевого функционала по управляющим воздействиям. В качестве управляющих воздействий рассматриваются краевые управления и минимизируемый параметр Т.
Предложенный численный метод может быть применен при исследовании задач оптимального управления механическими колебаниями струн и стержней, колебаниями электромагнитного поля, распространением акустических волн и многими другими процессами.
Выполненные в работе исследования могут быть далее развиты в направлении разработки методов численного решения задач быстродействия для волнового уравнения с управляющими воздействиями других типов, с управлениями, удовлетворяющими определенным ограничениям, а также для многомерных уравнений колебаний.
СПИСОК ЛИТЕРА ТУРЫ
1. Бутковский А. Т. Структурная теория распределенных систем. М.: Наука, 1977. 320 с.
2. Ильин В. А. Граничное управление процессом колебаний на двух концах в терминах обобщенного решения волнового уравнения с конечной энергией // Дифференциальные уравнения. 2000. Т. 36. № 11. С. 1513-1528.
3. Ильин В. А., Моисеев Е. И. Оптимизация граничных управлений колебаниями струны // Успехи математических наук. 2005. Т. 60. Вып. 6. С. 89-114.
4. Никитин А. А. Оптимальное граничное управление колебаниями струны, производимое силой при упругом закреплении // Дифференциальные уравнения. 2011. Т. 47. № 12. С. 1773-1782.
5. Знаменская Л. Н. Управление упругими колебаниями. М.: Физматлит, 2004. 176 с.
6. Знаменская Л. Н. Управление колебаниями струны в классе обобщенных решений из L2 // Дифференциальные уравнения. 2002. Т. 38. № 5. С. 666-672.
7. Чабакаури Г. Д. Оптимизация граничного управления процессом колебаний на одном конце при закрепленном втором конце // Дифференциальные уравнения. 2001. Т. 37. № 12. C. 1655-1663.
8. Кулиев Г. Ф. Задача точечного управления для гиперболического уравнения // Автоматика и телемеханика. 1993. № 3. С. 80-84.
9. Лионе Ж.-Л. Некоторые вопросы оптимального управления распределенными системами // Успехи математических наук. 1985. Т. 40. Вып. 4. С. 55-68.
10. Айда-заде К. Р., Асадова Д. А. Исследование переходных процессов в нефтепроводах // Автоматика и телемеханика. 2011. № 12. С. 156-172.
11. Иванов Д. А., Потапов М. М. Приближения к оптимальным по времени граничным управлениям для слабых обобщенных решений волнового уравнения // Журнал вычислительной математики и математической физики. 2017. Т. 57. № 4. С. 605-624.
12. Гибкина Н. В., Сидоров М. В., Стадникова А. В. Оптимальное граничное управление колебаниями однородной струны // Радиотехника. 2016. № 2. C. 3-11.
13. Gugat M., Leugering G., Sklyar G. Lp - optimal boundary control for the wave equation // SIAM Journal Control and Optimization. 2005. V. 44. N. 1. P. 49-74.
14. Ervedoza S., Zuazua E. Numerical approximation of exact controls for waves. New York: Springer, 2013. 122 p.
15. Loheac J., Zuazua E. Norm saturating property of time optimal controls for wave-type equations // IFAC - PapersOnLine. 2016. V. 49. N. 8. P. 37-42.
16. Kunisch K., Wachsmuth D. On time optimal control of the wave equation and its numerical realization as parametric optimization problem // SIAM Journal Control and Optimization. 2013. V. 51. N. 2. P. 1232-1262.
17. Васильев Ф. П., Иванов Р. П. О приближенном решении задачи быстродействия в банаховых пространствах при наличии ограничений на фазовые координаты // Журнал вычислительной математики и математической физики. 1971. Т. 11. № 2. С. 328-347.
18. Тихонов А. Н., Самарский А. А. Уравнения математической физики. М.: Наука, 1977. 736 с.
19. Васильев Ф. П. Методы оптимизации. М.: Факториал Пресс, 2002. 824 с.
Статья поступила в редакцию 10.04.2019
ИНФОРМАЦИЯ ОБ АВТОРЕ
Асадова Джамиля Абульфаз кызы - Азербайджан, Л7 1141, Баку; Институт систем управления Национальной Академии наук Азербайджана; доцент; ведущий научный сотрудник отдела разработки численных методов принятия решений в детерминированных системах; j [email protected].
ANALYSIS OF MINIMUM TIME BEFORE STABILIZATION OF OSCILLATORY PROCESS UNDER DIFFERENT MODES
J. A. Asadova
Institute of Control Systems of the Azerbaijan Academy of Sciences, Baku, Azerbaijan
Abstract. The paper presents a numerical approach to solving the problems of optimal control over establishing oscillatory systems by means of boundary controls. In contrast to the reviewed works of other authors on similar topics, the problem of establishing the oscillatory process is numerically studied in a temporal neighborhood of the required mode, the friction force being taken into account in the mathematical model of the oscillatory system. The schemes of deriving
the basic formulas are presented, as well as algorithms for their numerical solution by first-order optimization methods. The corresponding analytical formulas for the components of the gradient of the target functional for control actions are obtained. A comparison is made of the minimum time for the oscillation process to be established both with its optimal control and with the process self-determination. The results of numerical experiments are presented. The results can be attributed to the different processes described by hyperbolic equations.
Key words: oscillatory systems, boundary control actions, gradient of target functional, minimal settling time of oscillatory process.
For citation: Asadova J. A. Analysis of minimum time before stabilization of oscillatory processes under different modes. Vestnik of Astrakhan State Technical University. Series: Management, Computer Science and Informatics. 2019;4:7-17. (In Russ.) DOI: 10.24143/2072-95022019-4-7-17.
REFERENCES
1. Butkovskii A. T. Strukturnaia teoriia raspredelennykh sistem [Structural theory of distributed systems]. Moscow, Nauka Publ., 1977. 320 p.
2. Il'in V. A. Granichnoe upravlenie processom kolebanii na dvukh koncakh v terminakh obobshchennogo resheniia volnovogo uravneniia s konechnoi energiei [Boundary control of oscillation process on two ends in terms of generalized solution of wave equation with finite energy]. Differencial'nye uravneniia, 2000, vol. 36, no. 11, pp. 1513-1528.
3. Il'in V. A., Moiseev E. I. Optimizaciia granichnykh upravlenii kolebaniiami struny [Optimization of boundary controls for string vibrations]. Uspekhi matematicheskikh nauk, 2005, vol. 60, iss. 6, pp. 89-114.
4. Nikitin A. A. Optimal'noe granichnoe upravlenie kolebaniiami struny, proizvodimoe siloi pri uprugom zakreplenii [Optimal boundary control of string vibrations produced by force with elastic fastening]. Differencial'nye uravneniia, 2011, vol. 47, no. 12, pp. 1773-1782.
5. Znamenskaia L. N. Upravlenie uprugimi kolebaniiami [Elastic vibration control]. Moscow, Fizmatlit Publ., 2004. 176 p.
6. Znamenskaia L. N. Upravlenie kolebaniiami struny v klasse obobshchennykh reshenii iz L2 [Control of string vibrations in class of generalized solutions from L2]. Differencial'nye uravneniia, 2002, vol. 38, no. 5, pp. 666-672.
7. Chabakauri G. D. Optimizaciia granichnogo upravleniia processom kolebanii na odnom konce pri zakreplennom vtorom konce [Optimization of boundary control of oscillation process at one end with the second end fixed]. Differencial'nye uravneniia, 2001, vol. 37, no. 12, pp. 1655-1663.
8. Kuliev G. F. Zadacha tochechnogo upravleniia dlia giperbolicheskogo uravneniia [Point control problem for hyperbolic equation]. Avtomatika i telemekhanika, 1993, no. 3, pp. 80-84.
9. Lionc Zh.-L. Nekotorye voprosy optimal'nogo upravleniia raspredelennymi sistemami [Some issues of optimal management of distributed systems]. Uspekhi matematicheskikh nauk, 1985, vol. 40, iss. 4, pp. 55-68.
10. Aida-zade K. R., Asadova D. A. Issledovanie perekhodnykh processov v nefteprovodakh [Studying transient processes in oil pipelines]. Avtomatika i telemekhanika, 2011, no. 12, pp. 156-172.
11. Ivanov D. A., Potapov M. M. Priblizheniia k optimal'nym po vremeni granichnym upravleniiam dlia slabykh obobshchennykh reshenii volnovogo uravneniia [Approximations to time-optimal boundary controls for weak generalized solutions of wave equation]. Zhurnal vychislitel'noi matematiki i matematicheskoi fiziki, 2017, vol. 57, no. 4, pp. 605-624.
12. Gibkina N. V., Sidorov M. V., Stadnikova A. V. Optimal'noe granichnoe upravlenie kolebaniiami od-norodnoi struny [Optimal boundary control of vibrations of uniform string]. Radiotekhnika, 2016, no. 2, pp. 3-11.
13. Gugat M., Leugering G., Sklyar G. Lp - optimal boundary control for the wave equation. SIAM Journal Control and Optimization, 2005, vol. 44, no. 1, pp. 49-74.
14. Ervedoza S., Zuazua E. Numerical approximation of exact controls for waves. New York, Springer, 2013. 122 p.
15. Loheac J., Zuazua E. Norm saturating property of time optimal controls for wave-type equations. IFAC - PapersOnLine, 2016, vol. 49, no. 8, pp. 37-42.
16. Kunisch K., Wachsmuth D. On time optimal control of the wave equation and its numerical realization as parametric optimization problem. SIAM Journal Control and Optimization, 2013, vol. 51, no. 2, pp. 1232-1262.
17. Vasil'ev F. P., Ivanov R. P. O priblizhennom reshenii zadachi bystrodeistviia v banakhovykh pros-transtvakh pri nalichii ogranichenii na fazovye koordinaty [Approximate solution to performance problem in Banach spaces under limits on phase coordinates]. Zhurnal vychislitel'noi matematiki i matematicheskoi fiziki, 1971, vol. 11, no. 2, pp. 328-347.
18. Tikhonov A. N., Samarskii A. A. Uravneniia matematicheskoi fiziki [Equations of mathematical physics]. Moscow, Nauka Publ., 1977. 736 p.
19. Vasil'ev F. P. Metody optimizacii [Optimization methods]. Moscow, Faktorial Press, 2002. 824 p.
The article submitted to the editors 10.04.2019
INFORMA TION ABOUT THE A UTHOR
Asadova Jamila Abulfaz - Azerbaijan, AZ 1141, Baku; Institute of Control Systems of the Azerbaijan Academy of Sciences; Assistant Professor; Leading Researcher of the Department of Development of Numerical Methods for Decision-making in Deterministic Systems; j [email protected].
