УДК 517.977 Вестник СПбГУ. Прикладная математика. Информатика... 2024. Т. 20. Вып. 2
МБС 49Л20
Неклассическое условие оптимальности в гибридной задаче управления гиперболическими и обыкновенными дифференциальными уравнениями с запаздыванием*
А. В. Аргучинцев, В. П. Поплевко
Иркутский государственный университет,
Российская Федерация, 664003, Иркутск, ул. К. Маркса, 1
Для цитирования: Аргучинцев А. В., Поплевко В. П. Неклассическое условие оптимальности в гибридной задаче управления гиперболическими и обыкновенными дифференциальными уравнениями с запаздыванием // Вестник Санкт-Петербургского университета. Прикладная математика. Информатика. Процессы управления. 2024. Т. 20. Вып. 2. С. 255-264. https://doi.org/10.21638/spbu10.2024.210
Рассматривается задача оптимального управления линейной гиперболической системой первого порядка, в которой неоднородность в правой части определяется из управляемой линейной системы обыкновенных дифференциальных уравнений с постоянным запаздыванием. Матрица коэффициентов при фазовых переменных в системе обыкновенных дифференциальных уравнений зависит от функции управления. Целевой функционал линеен. На основе точной (без остаточных членов) формулы приращения целевого функционала задача сведена к задаче оптимального управления системой обыкновенных дифференциальных уравнений. Результат сформулирован в виде неклассического условия оптимальности вариационного типа. Предлагаемая редукция задачи существенно сокращает объем вычислений при использовании численных методов оптимизации. Приведен иллюстративный пример.
Ключевые слова: гибридная задача, гиперболическая система, система с запаздыванием, точная формула приращения, вариационное условие оптимальности, редукция задачи.
1. Введение. История исследования задач управления дифференциальными уравнениями с частными производными и эффектом запаздывания главным образом повторяет и развивает полученные ранее результаты для задач управления системами с сосредоточенными параметрами. Можно выделить достижения в вопросах существования и единственности оптимальных управлений [1], получении условий оптимальности типа классического принципа максимума [2], построении вычислительных методов [3], синтеза [4] и стабилизации [5-7].
В настоящей статье описывается задача оптимального управления линейной гиперболической системой первого порядка, в которой неоднородность в правой части определяется из управляемой линейной системы обыкновенных дифференциальных уравнений с постоянным запаздыванием. Задачи такого вида возникают, в частности, в моделях, описывающих динамику двух взаимодействующих популяций «растительноядный консумент — растения» [8, с. 218-219]. Структура растительноядной популяции животных определяется интегродифференциальным уравнением с гиперболическим дифференциальным оператором в левой части. Для этой популяции важно распределение ее членов по возрасту. В качестве независимых переменных выступают время, в течение которого изучается процесс, и возраст особей. Для модели-
* Исследование А. В. Аргучинцева выполнено за счет гранта Российского научного фонда № 23-21-00296, https://rscf.ru/project/23-21-00296/
( Санкт-Петербургский государственный университет, 2024
рования же динамики популяции растения возраст часто не имеет существенного значения, а поэтому можно ограничиться классическим обыкновенным дифференциальным уравнением.
Особенностью рассматриваемой задачи является предположение о том, что матрица коэффициентов при фазовых переменных в системе обыкновенных дифференциальных уравнений зависит от функции управления. Поэтому, несмотря на линейность целевого функционала, классический принцип максимума Л. С. Понтрягина в данной задаче является только необходимым, но не достаточным условием оптимальности. Для решения подобных задач обычно применяются методы, разработанные для общих нелинейных проблем.
В статье впервые для задач с запаздыванием применена модифицированная методика [9], основанная на использовании точной (без остаточных членов) формулы приращения целевого функционала. Этот подход позволил свести исходную задачу оптимального управления к задаче оптимального управления системой обыкновенных дифференциальных уравнений. Результат сформулирован в виде неклассического условия оптимальности вариационного типа. Доказанная теорема по существу определяет целое семейство необходимых и достаточных условий оптимальности. Редукция задачи дает возможность в значительной степени сократить объем вычислений. При реализации оптимизационных методов в подобного рода задачах основное время тратится на многократные интегрирования исходной и сопряженной систем дифференциальных уравнений в частных производных. Полученный в настоящей работе результат позволяет ограничиться максимум двумя процедурами достаточно трудоемкого решения начально-краевых задач для систем гиперболических уравнений: однократное интегрирование сопряженной системы и, если это необходимо, финальное решение для определения состояния процесса, соответствующего найденному оптимальному управлению. Приведенный в п. 5 пример иллюстрирует такой подход.
2. Постановка задачи. Рассмотрим систему линейных гиперболических уравнений первого порядка
/77* Г/Т _
— + = Ф(М)х + ЯМ) + Ф)у, (1)
(в,г) е п, п = 5 х т, 5 =[в0,в1], т = [г0,^],
в которой х(в, г), у(г) — п- и т-мерные вектор-функции; А(в,г) — диагональная матрица порядка (п х п); С (г) и Ф(в, г) — матричные функции размера (п х т), (п х п) соответственно; — п-мерная вектор-функция. Дополнительно считаем, что диа-
гональные элементы а = г = 1, 2,..., п, матрицы коэффициентов знакопо-
стоянны в п:
> 0, г = 1, 2,...,т1,
а$(в,г)=0, г = т1 + 1,т1 + 2, ...,т2,
аДв, г) < 0, г = т2 + 1, т-2 + 2,. .., п.
Пусть матрица А+(в,г) порядка (т1 х т1) и матрица А-(в,г) порядка (п — тг) х (п - т2 ) — диагональные подматрицы, составленные из положительных и отрицательных элементов матрицы А. Рассмотрим подвекторы х+ = (х1, х2,..., хт1) и х- = (хт2+1,хт2+2,... ,хп), которые отвечают положительным и отрицательным элементам матрицы А.
Начально-краевые условия для системы (1) зададим в точках начала характеристик гиперболических уравнений:
x(s,t0) = x0(s), s E S; x+ (s0,t) = n(t), x-(sut) = ^(t), t E T, (2)
с условиями согласования n(t0) = (x0(s0))+, p(t0) = (x0(si))-.
Функция y(t) определяется из управляемой системы обыкновенных дифференциальных уравнений с запаздыванием
— = B(u(t),t)y(t -а)+ d(u(t),t), t Е Т, (3)
dt
y(t)= y0(t), t E [t0 - a,t0],
где B(u(t),t) — матричная функция размера (m x m); d(u(t),t) — m-мерная вектор-функция; a — заданная положительная константа.
В качестве множества допустимых управлений выберем совокупность ограниченных и измеримых на T r-мерных вектор-функций u(t), удовлетворяющих почти всюду на этом отрезке ограничениям типа включения
u(t) E U С Er, t E T, (4)
где U — компактное множество.
Требуется найти допустимое управление, доставляющее минимум целевому функционалу:
J(u) = j(e(s),x(s,t1)) ds + JJ(b0(s,t),x(s,t)) dsdt. (5)
s n
Здесь и далее угловыми скобками обозначается классическое скалярное произведение в конечномерном евклидовом пространстве соответствующей размерности.
Задача оптимального управления (1)-(5) рассматривается при следующих предположениях:
1) диагональные элементы ai (s,t), i = 1, 2,...,n, матрицы A непрерывно дифференцируемы в прямоугольнике П;
2) вектор-функции x0(s), ^(t) и n(t) непрерывны на S и T соответственно;
3) матричные функции Ф(в, t), B(u,t), C(t), а также вектор-функции /(s,t), d(u(t),t) непрерывны по совокупности своих аргументов в областях их определения.
При указанных предположениях на параметры задачи для любого допустимого управления существует единственное обобщенное решение начально-краевой задачи (1), (2) из класса непрерывных в П функций, причем каждая компонента решения xi, i = 1, 2,..., n, непрерывно дифференцируема вдоль характеристик системы гиперболических уравнений [10].
3. Формула приращения функционала. Задача (1)-(5) является линейной, однако классический принцип максимума Л. С. Понтрягина не будет в ней достаточным условием оптимальности. Причина этого — зависимость от управления матрицы коэффициентов в (3). Обычно для решения подобных задач применяются те же методы, что и для общего нелинейного случая. В настоящей работе используется методика, ранее примененная в [9] для задач при отсутствии запаздывания.
Рассмотрим два произвольных различных допустимых процесса: {u, x} и {м = u + Au, x = x + Ax}. Введем обозначение для дифференциального оператора в левой
части (1): Dax = xt + A(s,t)xs. Здесь Dax = (Dixi,..., Dnxn) — обобщенная производная, каждая компонента которой DjXj непрерывна вдоль соответствующего г-го семейства характеристик. Уравнение в приращениях для исходной начально-краевой задачи примет следующий вид:
DAAx = Ф(й,г)Дх + C (t)Ay,
Ax(s,t0) = Ax+(s0, t) = Ax-(sbt) = 0,
Ayt = B(w(t),t)y(t - a) - B(u(t),t)y(t - a) + Ad(u(t), t), Ay(t0) = 0. (6)
Здесь Ad(u(t),t) = d(w(t),t) - d(u(t),t).
Перепишем правую часть (6) таким образом:
B(w(t), t)y(t - a) - B(u(t), t)y(t - a) + Ad(u(t), t) =
= AB(u(t), t)y(t - a) + B(u(t), t)Ay(t - a) + Ad(u(t), t), где AB(u(t),t) = B(w(t),t) - B(u(t),t). Приращение функционала запишем так:
AJ(u) = J(и) - J(u) = J(ß(s), Ax(s,t1)) ds + JJ(b0(s, t), Ax(s,t)) dsdt. (7) s n
В правую часть формулы (7) добавим нулевые слагаемые
JJ(^(s, t), DaAx - Ф(в, t)Ax) dsdt n
и
y(p(t), Ayt - Aü_B(u,t)y(t - a) - B(u,t)Ay(t - a) - Ad(u,t)) dt, T
где ^(s,t) и p(t) — пока произвольные вектор-функции.
Применяя формулы интегрирования по частям, получим, что
AJ(u) = J (ß(s), Ax(s,t1)) ds + JJ(b0(s, t), Ax(s, t)) dsdt + s n
^У(^(s,t1), Ax(s,t1)) - (^(s,t0), Ax(s,t0)) ds -S
^y"(^(s0,t),A(s0,t)Ax(s0,t)) - (^(s1,t), A(s1,t)Ax(s1,t)) dt -
- JJ (Da^ + + Фт Ax) dsdt - JJ (^(s,t), CAy)dsdt + n n
+ (p(t1), Ay(t1)) - J(pt, Ay(t)) dt -
T
- J(p(t), AB(u(t), t)y(t - a) + B(u(t), t)Ay(t - a) + Ad(u, t)) dt.
T
Далее рассмотрим выражение
т = t — a
(p(t),B(u,t)Ay(t - a)) dt-, т E [t0 - a, ti - a]
T
to
J (р(т + a), B(u(t + a), т + a)Ay(T)) d.T + (8)
to—a
ti —a
+ J (р(т + a), B(u(t + a),T + a)Ay(T)) d,T.
to
Перейдем обратно к переменной t. С учетом того, что первое слагаемое (8) равно нулю, приходим к следующему выражению:
J(p(t), B(u, t)Ay(t - a)) dt =
T 11 —a
J (p(t + a), B(u(t + a),t + a)Ay(t)) dt.
Заметим, что
Я 1
11 (ф(8, САу) ¿зА = - Л (§; IШ, сау)
П П я
Применяя формулу интегрирования по частям, получаем, что
J!(ф(в, г), сАу) ¿заг = ^^ ф(в, г) ¿в, САу) ¿г.
П Т Б
Потребуем, чтобы функции ф(в,г), р(г) были решениями сопряженной задачи
Блф + Аяф = -ФТ(в,г)ф + Ьо(в,г), (в,г) е п, (9)
ф(в,г1) = -в(в), в е ф(въг) = 0, ф-(в0,г)=0, ф+(въг) = 0, г е т,
{-БТ(и (г + а),г + а)р(г + а) - / сТф(в,г) ¿в, г е [0; г1 - а), -/ сТф(в,г) ¿в, г е [г1 - а;г1],
р(г1) = 0. (10)
Тогда формула приращения функционала примет вид
А1 (и) = - !(р(г), АБ(и,г)у(г - а) + А¿(u,г)) ¿г. (11)
Т
Формула (11) является точной (без остаточных членов) для любой пары допустимых процессов, при этом исходная система обыкновенных дифференциальных уравнений (3) интегрируется на возмущенном управлении.
4. Условие оптимальности. На основе полученной формулы приращения (11) можно осуществить редукцию исходной задачи оптимального управления гиперболической системой к задаче оптимального управления системами обыкновенных дифференциальных уравнений:
I(v) = -y"(p(t,u), (B(v(t),t) - B(u(t),t))y(t - a,v) + (12)
T
+ d(v(t), t) - d(u(t), t)) dt ^ min, yt = B(v(t), t)y(t - a) + d(v(t), t), t e T,
y(t) = y0(t), t e [to - a,to], (13)
v(t) e U.
Здесь y(t) — m-мерная функция состояния, u(t) и p(t,u) — фиксированные функции, v(t) — управляющее воздействие, удовлетворяющее тем же ограничениям на управления, что и в исходной задаче. Тем самым доказана справедливость следующего утверждения.
Теорема. Для оптимальности управления u*(t) в задаче (1)-(5) необходимо и достаточно, чтобы управление v* = u*(t) было оптимальным в задаче (12), (13). При этом оптимальное значение функционала определяется формулой
J (v*) = J (u) + I (v*).
Полученный результат можно рассматривать и как вариационное условие оптимальности (в отличие от традиционных конечномерных условий оптимальности). Отметим, что теорема справедлива для любого допустимого управления u(t).
Опишем схему решения исходной задачи.
1. Задается произвольное допустимое управление u = u(t). Вычисляется соответствующее ему решение p = p(t,u) сопряженной задачи (10). Для этого необходимо также найти решение задачи (9), которое не зависит от выбора допустимого управления.
2. Решается вспомогательная задача оптимального управления системой обыкновенных дифференциальных уравнений с линейным целевым функционалом (12).
Таким образом, для решения задачи (1)-(5) необходимо всего лишь два раза проинтегрировать систему дифференциальных уравнений с частными производными (поиск -0(s,t) и при необходимости x* (s,t,v*)).
5. Иллюстративный пример. Рассмотрим пример применения описанного алгоритма. В квадрате S х T = [0; 2] х [0; 2] поставим следующую задачу оптимального управления:
Xlt + Xls = Х2 - y, X2t - X2s = (s + 2) sint, yt = u • y(t - 0.15), u(t) e [0, 3], y(t) = t, t e [-0.15; 0],
xl (0,t)=0, xi(s, 0) = s, x2 (2,t)=0, x2(s, 0) = 0. Поставим задачу на минимум целевого функционала
J(u) = j[3xl(s, 2) + 2x2(s, 2)] ds, u e U.
В принятых обозначениях
$(s,t) ^0 0 ) ' C(t) = ( -l) , B(u, t) = u, d(u,t) = 0. Сопряженная задача принимает вид
фи + 'Фи =0, 'ф-2Ь - ф2з = -Ф1, ф1(s, 2) = -3, ф2(s, 2) = -2, фi(2,t) = 0, ф2(0,t) = 0, ' -p(t + 0.15) • u(t + 0.15) + f ф1 ds, t e [0; 1.85),
Pt
/ф1 ds, t e [1.85; 2],
.s
p(2) = 0.
Выберем допустимое управление u(t) = 0 для всех t e [0, 2]. Тогда решение сопряженной задачи можно записать так:
ф^>={-3; t С;
3
p(t) = --t2+ 6, t £ [0; 2]. Исходная задача может быть сведена к задаче оптимального управления
2
I(v) = J(1.5t2 - 6) • v(t) • z(t - 0.15, v) dt ^ min,
0
zt = v • z(t - 0.15), z(t)= t, t e [-0.15; 0], v(t) e [0; 3].
6. Заключение. Отметим, что редуцированная задача оптимального управления (12), (13) имеет специальную структуру. Это также задача с запаздыванием, в которой матрица коэффициентов при фазовых переменных зависит от состояния. Для решения данной задачи можно использовать широкий арсенал методов, разработанных для оптимизации систем с сосредоточенными параметрами. В частности, билинейность задачи позволяет эффективно применять алгоритмы, разработанные в [11]. В качестве дальнейшего направления работы можно предложить исследование случая квадратичного целевого функционала с помощью методики [12]. Однако
случай квадратичного критерия качества требует применения формул приращения второго порядка и матричных импульсов Р. Ф. Габасова.
Литература
1. Teo K. Optimal control of systems governed by time-delayed, second-order, linear, parabolic partial differential equations with a first boundary condition // Journal of Optimization Theory and Applications. 1979. Vol. 29. N 3. P. 437-481.
2. Sadek I. Optimal control of time-delay systems with distributed parameter // Journal of Optimization Theory and Applications. 1990. Vol. 67. N 3. P. 567-585.
3. Mai T., Nguyen H., Nguyen V., Vu V. Applying Pade approximation model in optimal control problem for a distributed parameter system with time delay // International Journal of Computing and Optimization. 2017. Vol. 4. N 1. P. 19-30. https://doi.org/10.12988/ijco.2017.61227
4. Провоторов В. В., Провоторова Е. Н. Синтез оптимального граничного управления параболической системы с запаздыванием и распределенными параметрами на графе // Вестник Санкт-Петербургского университета. Прикладная математика. Информатика. Процессы управления. 2017. Т. 13. Вып. 2. C. 209-224. https://doi.org/10.21638/11701/spbu10.2017.207
5. Liu D., Wang Q., Xu G. Stabilization of distributed parameter systems with delays in the boundary based on predictors // IEEE Transactions on Automatic Control. 2021. Vol. 66. N 7. P. 33173324. https://doi.org/10.1109/TAC.2020.3035586
6. Гребенщиков Б. Г. Асимптотические свойства и стабилизация одной системы нейтрального типа с постоянным запаздыванием // Вестник Санкт-Петербургского университета. Прикладная математика. Информатика. Процессы управления. 2021. Т. 17. Вып. 1. C. 81-96. https://doi.org/10.21638/11701/spbu10.2021.108
7. Furtat I., Orlov Yy. Synchronization and state estimation of nonlinear systems with unknown time-delays: Adaptive identification method // Cybernetics and Physics. 2020. Vol. 9. N 3. P. 136-143. https://doi.org/10.35470/2226-4116-2020-9-3-136-143
8. Алексеев В. В., Крышев И. И., Сазыкина Т. Г. Физическое и математическое моделирование экосистем. СПб.: Гидрометеоиздат, 1992. 368 c.
9. Аргучинцев А. В. Оптимальное управление гиперболическими системами. М.: Физматлит, 2007. 186 с.
10. Рождественский Б. Л., Яненко Н. Н. Системы квазилинейных уравнений и их приложения к газовой динамике. М.: Наука, 1978. 592 с.
11. Срочко В. А., Аксенюшкина Е. В. Параметризация некоторых задач управления линейными системами // Известия Иркутского государственного университета. Сер. Математика. 2019. Т. 30. С. 83-98. https://doi.org/10.26516/1997-7670.2019.30.83
12. Аргучинцев А. В. Вариационное условие оптимальности в задаче минимизации нормы конечного состояния составной системой гиперболических и обыкновенных дифференциальных уравнений // Вестник Санкт-Петербургского университета. Прикладная математика. Информатика. Процессы управления. 2023. Т. 19. Вып. 4. C. 540-548. https://doi.org/10.21638/11701/spbu10.2023.410.
Статья поступила в редакцию 21 января 2024 г. ^атья принята к печати 12 марта 2024 г.
Контактная информация:
Аргучинцев Александр Валерьевич — д-р физ.-мат. наук, проф.; [email protected] Поплевко Василиса Павловна — канд. физ.-мат. наук, доц.; [email protected]
The non-classical optimality condition in the hybrid control problem of hyperbolic and ordinary differential equations with delay*
A. V. Arguchintsev, V. P. Poplevko
Irkutsk State University, 1, ul. K. Marksa, Irkutsk, 664003, Russian Federation
* The research of A. V. Arguchintsev was supported by the Russian Science Foundation, project N 23-21-00296, https://rscf.ru/project/23-21-00296/
For citation: Arguchintsev A. V., Poplevko V. P. The non-classical optimality condition in the hybrid control problem of hyperbolic and ordinary differential equations with delay. Vestnik of Saint Petersburg University. Applied Mathematics. Computer Science. Control Processes, 2024, vol. 20, iss. 2, pp. 255-264. https://doi.org/10.21638/spbu10.2024.210 (In Russian)
In this paper, we consider an optimal control problem for a system of linear first-order hyperbolic equations in which the inhomogeneity in the right-hand side is determined from the controlled linear system of ordinary differential equations with constant delay. The coefficient matrix at phase variables in the system of ordinary differential equations depends on the control function. The cost functional is linear. On the basis of the exact increment formula (without remainder terms) of the cost functional, the problem is reduced to the optimal control problem of a system of ordinary differential equations. The result is formulated in the form of a non-classical variational optimality condition. The proposed problem reduction significantly reduces the amount of calculations when using numerical optimization methods. An illustrative example is given.
Keywords: hybrid problem, hyperbolic system, delayed system, exact increment formula, variational optimality condition, problem reduction.
References
1. Teo K. Optimal control of systems governed by time-delayed, second-order, linear, parabolic partial differential equations with a first boundary condition. Journal of Optimization Theory and Applications, 1979, vol. 29, no. 3, pp. 437-481.
2. Sadek I. Optimal control of time-delay systems with distributed parameter. Journal of Optimization Theory and Applications, 1990, vol. 67, no. 3, pp. 567-585.
3. Mai T., Nguyen H., Nguyen V., Vu V. Applying Pade approximation model in optimal control problem for a distributed parameter system with time delay. International Journal of Computing and Optimization, 2017, vol. 4, no. 1, pp. 19-30. https://doi.org/10.12988/ijco.2017.61227
4. Provotorov V. V., Provotorova E. N. Sintez optimal'nogo granichnogo upravleniia parabolicheskoii sistemi s zapazdyvaniem i raspredelennymi parametrami na grafe [Synthesis of optimal boundary control of parabolic systems with delay and distributed parameters on the graph]. Vestnik of Saint Petersburg University. Applied Mathematics. Computer Science. Control Processes, 2017, vol. 13, iss. 2, pp. 209-224. https://doi.org/10.21638/11701/spbu10.2017.207 (In Russian)
5. Liu D., Wang Q., Xu G. Stabilization of distributed parameter systems with delays in the boundary based on predictors. IEEE Transactions on Automatic Control, 2021, vol. 66, no 7, pp. 3317-3324. https: //doi.org/10.1109/TAC.2020.3035586
6. Grebenshchikov B. G. Asimptoticheskie svoistva i stabilizaciya odnoi sistemy neitral'nogo tipa s postoyannym zapazdyvaniem [Asymptotic properties and stabilization of a neutral type system with constant delay]. Vestnik of Saint Petersburg University. Applied Mathematics. Computer Science. Control Processes, 2021, vol. 17, iss. 1, pp. 81-96. https://doi.org/10.21638/11701/spbu10.2021.108 (In Russian)
7. Furtat I., Orlov Y. Synchronization and state estimation of nonlinear systems with unknown time-delays: adaptive identification method. Cybernetics and Physics, 2020, vol. 9, no. 3, pp. 136-143. https://doi.org/10.35470/2226-4116-2020-9-3-136-143
8. Alekseev V. V., Kryshev I. I., Sazykina T. G. Fizicheskoe i matematicheskoe modelirovanie ehkosistem [Physical and mathematical modeling of ecosystems]. St. Petersburg, Hydrometeoizdat Publ., 1992, 368 p. (In Russian)
9. Arguchintsev A. V. Optimalnoe upravlenie giperbolicheskimi sistemami [Optimal control of hyperbolic systems]. Moscow, Fizmatlit Publ., 2007, 186 p. (In Russian)
10. Rozhdestvenskiyi B. L., Yanenko N. N. Systemi kvazilineinykh uravnenii i ikh prilozheniya k gazovoi dinamike [Systems of quasilinear equations and their applications to gas dynamics]. Moscow, Nauka Publ., 1978, 592 p. (In Russian)
11. Srochko V. A., Aksenyushkina E. V. Parametrizatsiya nekotorykh zadach upravleniya lineynymi sistemami [Parameterization of some linear systems control problems]. The Bulletin of Irkutsk State University. Series Mathematics, 2019, vol. 30, pp. 83-98. https://doi.org/10.26516/1997-7670.2019.30.83 (In Russian)
12. Arguchintsev A. V. Variacionnoe uslovie optimalnosti v zadache minimizacii normy konech-nogo sostoyaniia sostavnoi sistemoi giperbolicheskikh i obyknovennykh differencialnykh uravnenii [The variational optimality condition in the problem of minimizing the finite state norm by a composite system
of hyperbolic and ordinary differential equations]. Vestnik of Saint Petersburg University. Applied Mathematics. Computer Science. Control Processes, 2023, vol. 19, iss. 4, pp. 540—548. https://doi.org/10.21638/11701/spbu10.2023.410 (In Russian)
Received: January 21, 2024. Accepted: March 12, 2024.
Authors' information:
Alexander V. Arguchintsev — Dr. Sci. in Physics and Mathematics, Professor; [email protected]
Vasilisa P. Poplevko — PhD in Physics and Mathematics, Associate Professor; [email protected]