Алгоритм решения краевой задачи для нелинейной управляемой системы и его численное моделирование Текст научной статьи по специальности «Математика»

УДК 531.36 Вестник СПбГУ. Математика. Механика. Астрономия. 2017. Т. 4 (62). Вып. 4

МЯО 93С10

ЛЛГОРИТМ РЕШЕНИЯ КРАЕВОЙ ЗАДАЧИ ДЛЯ НЕЛИНЕЙНОЙ УПРАВЛЯЕМОЙ СИСТЕМЫ И ЕГО ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ

А. Н. Квитко, О. С. Фирюлина, А. С. Еремин

Санкт-Петербургский государственный университет,

Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7—9

Предложен достаточно удобный для численной реализации алгоритм построения дифференцируемой управляющей функции, гарантирующей перевод широкого класса нелинейных стационарных систем обыкновенных дифференциальных уравнений из начального состояния в заданное конечное состояние фазового пространства с учетом ограничений на управление. Получено конструктивное достаточное условие, при котором указанный перевод возможен. Рассмотрена конкретная практическая задача и проведено ее численное моделирование. Библиогр. 16 назв. Ил. 1.

Ключевые слова: граничные условия, стабилизация, управление, нелинейная система.

1. Введение и постановка задачи. Одной из проблем математической теории управления являются вопросы, связанные с разработкой точных или приближенных методов построения управляющих функций и соответствующих им траекторий, соединяющих заданные точки в фазовом пространстве. Этим исследованиям посвящено большое количество публикаций. Некоторые из них приведены в [1-14]. В настоящее время проблема решения граничных задач достаточно подробно изучена для линейных и нелинейных управляемых систем специального вида. Однако теория решения граничных задач для нелинейных управляемых систем общего вида еще недостаточно разработана. Основные усилия авторов данной статьи направлены на разработку удобного для численной реализации и устойчивого к погрешностям вычислений алгоритма решения граничных задач для широкого класса нелинейных управляемых систем обыкновенных дифференциальных уравнений в классе синтезирующих управлений, а также на нахождение конструктивного достаточного условия, гарантирующего существование решения этих задач. Поставленная цель достигнута сведением решения исходной задачи к решению задачи стабилизации линейной нестационарной системы специального вида и решению задачи Коши для вспомогательной системы обыкновенных дифференциальных уравнений.

Простота реализации полученного метода обеспечивается тем, что наиболее трудоемкая его часть, связанная с построением вспомогательной системы и решением задачи стабилизации вспомогательной системы выполняется аналитическими методами и может быть реализована средствами компьютерной алгебры. Традиционные итеративные методы построения программных управлений [1] требуют большого количества вычислений и не обладают устойчивостью к погрешностям вычислений.

Эффективность полученного алгоритма иллюстрируется на численном моделировании конкретной практической задачи.

Объектом исследования является управляемая система обыкновенных дифференциальных уравнений

* = 1'(х, и), (1)

© Санкт-Петербургский государственный университет, 2017

где х = (х1,..., хп) , и = (и\,..., иг) , х € К", и € Кг, г ^ п, время £ € [0,1],

Г € С4"(К" х Кг; К"), f =(/1,...,/п)г, (2)

г(0, 0) = 0, (3)

Я = (В, АВ, А2В,..., АП_1В), гапкЯ = п, А =—(0,0), В =—(0,0), (4)

1|и|| < N. (5)

Задача. Найти пару функций х(£) € С[0,1], и(£) € С[0,1], удовлетворяющих системе (1) и условиям

х(0) = 0, х(1) = X, Х=(Х1,...,Х" )т. (6)

Указанную пару х(£), и(£) будем называть решением задачи (1), (6).

Теорема. Пусть для правой части системы (1) выполнены условия (2)-(4). Тогда Зе > 0 такое, что Ух € К": ||х|| < е существует 'решение задачи (1), (6), которое может быть получено после решения задачи стабилизации линейной нестационарной системы с экспоненциальными коэффициентами и последующим решением задачи Коши для вспомогательной системы обыкновенных дифференциальных уравнений.

Главная идея доказательства теоремы состоит в том, что посредством преобразований зависимых и независимых переменных, решение исходной задачи сводится к решению задачи стабилизации нелинейной вспомогательной системы обыкновенных дифференциальных уравнений специального вида при постоянно действующих возмущениях. Для ее решения находится синтезирующее управление, обеспечивающее экспоненциальное убывание фундаментальной матрицы линейной части вспомогательной системы. На заключительном этапе осуществляется переход к исходным переменным.

2. Построение вспомогательной системы. Функцию х(£), входящую в решение задачи (1), (6) ищем в виде

Хг(£) = Ог(£) + Хi, % = 1, . . . , П. (7)

В новых переменных система (1) и граничные условия примут вид

а = Г(х + а, и), (8)

а(0) = -х, а(1) = 0. (9)

Пару функций а(£), и(£), удовлетворяющих системе (8) и условиям (9), будем называть решением задачи (8), (9).

Задача. Найти пару функций а(£) € С[0,1], и(£) € С[0,1], удовлетворяющую системе (8) и условиям

а(0) = -х, а(£) ^ 0 при £ ^ 1. (10)

Замечание 1. Переходя к пределу в решении задачи (8), (10) при £ ^ 1 получим решение задачи (8), (9).

Сделаем в системе (8) преобразование независимой переменной г:

г = 1 - е-ат, т е [0, +то), (11)

где а > 0 — некоторое фиксированное число, подлежащее определению. Тогда в новой независимой переменной т система (8) и условия (10) примут вид

— = ае-аЧ(х + с,(1), (12) ат

г(0) = -X, г(т) ^ 0 при г ^ то, (13)

г(т) = а(г(т)), а(т) = и(г(т)), т е [о, +то), (14)

где г = (01,..., сп)т, d = (¿1,..., ¿г)т. Пару функций г(т), ^т), удовлетворяющую системе (12) и условиям (13), будем называть решением задачи (12), (13). Имея решение задачи (12), (13), с помощью формул (11), (14) можно получить решение задачи (8), (10).

Введем обозначения:

с = X + в,г, d = в¿d, в, е [0,1], i = 1,...,п,

п г

\к\ ^^ /¿, |т| ^^ т,, к! = /^1! • ... • кп!, т! = т1! • ... • тг!.

¿=1 ¿=1

Используя свойство (2) и разложение правой части системы (1) в ряд Тейлора в окрестности точки (X, 0), систему (12) можно записать в виде

^ = «е-ат/»(х, 0) + ае-ат ]Г 0)с,- + ае^ ]Г

3=1 3 з=1 3

1 / п п „ 2 х п г

+ 5«- ее з^ь*»)-*» + *ее а^ь* »>-**+

\3=1к=1 3 к 3=1 к=1 3 к

3=1 к=1

Е1 д

—.-г-1-1---(х, .. . с^г]"11 . .. (¿™г +

„.,. . к!т! дхк1 ...дхП" дм^1 ... дм^ п 1 Г

|к| + |т| = 1 1

=4п-2

—-—г-т-^-(с, а)с?1 ...<£"<гГ (15)

|к| + |т|= 1 1

=4п-1

для всех i = 1,..., п.

Ограничим область изменения г(т) неравенством

||г(т)|| <С1, т е [0, то). (16)

Сделаем множество преобразований сдвигов функций с,(т), i = 1,..., п. Главная их цель состоит в том, чтобы в правой части системы, полученной в результате этих

преобразований, все слагаемые, не содержащие в явном виде степеней компонент с и ^ в области (5), (16) удовлетворяли оценке 0(е-4"ат||х||) при т ^ то и ||х|| ^ 0. На первом этапе выполним замену еДт) на с(1)(т) по формуле

с>(т )= е(1)(т) - е-ат /i(X, 0), % = 1,...,п. (17)

Пусть

.,. . . Я|к| + М/ + = ^-—ъ---, г= 1,...,п.

дх\г ...дх^ди™1 ...Эй™"

После подстановки (17) в левую и правую части системы (15) с учетом введенного обозначения получим систему

¿=1 и

1 " " „2 /

¿=1к=1 и " я -Р " " о2 х

^ ¿=1к=1 ^

¿=1 ¿=1 к=1 ^

¿=1 к=1 3 к

У ттЦ £>|к|+|т|Л(х,0)(с(11)-е-отЛ(х,0))к1...х ^^ к!т! 1

|к| + |т|=4"-2

х (е"1) - е-ат/¡(х, 0))кп^ . .. + 1

к

|к| + |т|=4"-1

У) Й)(с(11) - О))*1 • • • х

к!т!

х (е"1) - е-ат /¡(х, 0))к" . , % = 1, .. ., п. (18)

Из (13), (17) следует

е(1)(0) = + /¡(х, 0), % = 1,...,п. (19)

Используя (19) и индуктивный переход на к-м шаге получим искомое преобразование вида

е(к-1)(т ) = е(к)(т)+ е-кат ф(к)(х), ф(к) (0) = 0, % = 1,...,п. (20)

Вестник СПбГУ. Математика. Механика. Астрономия. 2017. Т. 4(62). Вып. 4 611

Если преобразования (20) применить (4n — 1) раз, объединить в полученной системе слагаемые, линейные по компонентам вектора c(4n-1) и содержащие коэффициенты e-íaT, i = 1,..., n, а также слагаемые, линейные по компонентам вектора d и содержащие коэффициенты e-íaT, i = 1,..., 2n, то согласно (18)—(20) будем иметь систему и начальные данные, которые в векторной форме можно записать следующим образом:

dc(4n-1)

= Pc(4n"1} + Qd + Ri (с^4™-1), d, х, г) + R2 (с^4™-1), d, х, т) +

+ R3(c(4n-1), d, т) + R4(c(4n-1), d, x, т), (21)

Ri = (Ri,..., R" )T, R2 = (R2,..., R" )T, R3 = (R3,..., R3 )T, R4 = (R4,..., R4 )T.

Функции R1 содержат все слагаемые, которые линейно зависят от компонент вектора c(4n-1) c коэффициентами e-íaT, i ^ n + 1, а также слагаемые, входящие в последнюю сумму правой части полученной системы, которые не содержат степеней компонент вектора d и имеют сумму степеней компонент вектора c(4n-1), равную единице.

Функции R2 содержат все слагаемые, которые линейно зависят от компонент вектора d с множителями e-íaT, i ^ 2n+1, а также слагаемые, входящие в последнюю сумму правой части полученной системы, которые не содержат степеней компонент вектора c(4n-1) и имеют сумму степеней компонент вектора d, равную единице.

Функции в R3 содержат все слагаемые, нелинейные по компонентам векторов c(4n-1) и d.

Функции R4 состоят из суммы слагаемых, не содержащих степеней компонент векторов c(4n-1) и d.

Кроме того, можем записать

P(x) = ae-aT (P1(x) + e-aTP2(x) + • • • + e-(n-1)aTPn_1(x)) ,

Q(x) = (Q1(x) + Q2(x) + • • • + e-(2"-1)aTQ2n-1(x)) , (22)

Pi(x) = ^(x,0), Pi(0) = A, Ql(x) = |^(x,0), Qi(0) = B,

c

(4n-1)(0) = —x + f(x, 0) — ф(2)(x) — ф(3)(x)-----^(4n-1)(x),

Ф(<) = (#,..., #)T, ф(<) (0)=0, i = 1,..., 4n — 1.

(23)

3. Оценка слагаемых правой части вспомогательной системы. Из построения системы (21) следует, что в области (5), (16) имеют место оценки

(25)

||Р,(х)|| ^ 0, (х)|| ^ 0 при i = 2, ...,п - 1, 3 =2,..., 2п - 1; (24)

|^1(г(4п-1), d, х,т)|| < е-(п+1)ат¿1|г(4п-1)|, ¿1 > 0, |^2(г(4п-1), d, X, т)|| < е-(2п+1)ат¿2|d|, ¿2 > 0,

|^з(г(4п-1), ^т)|| < е-ат(||г(4п-1)||2 + |^||2) , ¿з > 0. (26)

Кроме того, из условий (2), (3) и построения R4 следует

|^4(г(4п-1), d, X, т) || < ¿4|х|е-4пат, ¿4 > 0. (27)

612 Вестник СПбГУ. Математика. Механика. Астрономия. 2017. Т. 4(62). Вып. 4

Оценка (27) вытекает из представления д Г

№,0) = — (0,О)х, 6 = (в1Х1,...,впХп)Т, вг (Е [0, 1], 1=1,...,п. дх

Замечание 2. Обозначим через qг1 %-й столбец матрицы р1. Построим матрицу

81 = Ы, Plq1,..., Рк1 -1q1, q2, Р^?,..., Рк2qí,..., Р2г}.

Здесь к3-, 3 = 1,...,г, — максимальное количество столбцов вида qj, Р^Ц, ..., Ркз таких, что все столбцы матрицы 81 линейно независимы. Из условий (4) и (22) следует, что существует такое е1 > 0, что Ух € К" таких, что ||х|| < е1, выполняется гапк81 = п.

Рассмотрим систему

¿с(4"-1)

-= Рс(4п-1} + СИ (28)

ат

4. Вспомогательная лемма.

Лемма. Пусть для системы (1) выполнены условия (2), (4). Тогда Зе > 0: е < е1 такое, что Ух € К", ||х|| < е, существует управление ^т) вида

а(т ) = М(т )с(4"-1), (29)

при котором фундаментальная матрица системы (28) экспоненциально убывает.

Доказательство. При доказательстве используем метод Н. Н. Красовского стабилизации линейных нестационарных систем. Пусть ЬЦ, 3 = 1,...,г,— 3-й столбец матрицы р. Построим матрицу

82 = {Ь1, ь2 ,..., ьк1, Ь2, Ь2,..., Ь^ ,..., ц,..., цг}, (зо)

Ц=рЦ_1--1-±, о = 1,...,г, г = 2,...,Щ.

Здесь к3-, 3 = 1,..., г, — максимальное количество столбцов вида Ь*, Ь1,..., Ьк. таких, что столбцы 82 линейно независимы. Покажем, что выполняется

гапк82 = п. (31)

Пусть ЬЦ, 3 = 1,..., г, — 3-й столбец матрицы ае-атР1. Введем в рассмотрение матрицу

83 = {Ь1, ь 1,..., ь к1, Ь1, ь 2,..., Ь2к2,..., ь 1,..., ь кг}, (32)

Ц=ае-агР1Ц_1--] = \ ,...,г, ¡ = 2,...,^,

¿Я

¿т

где величины к3-, 3 = 1,..., г, определяются так же, как при построении матрицы 82. Условия (2), (22), (24) гарантируют 82 ^ 83 при ||х|| ^ 0. Отсюда следует существование е2 > 0: е2 < е1 такого, что Ух € К": ||х|| < е2 выполняется

гапк82 = гапк83. (33)

Рассуждая методом от противного с учетом замечания 2 и равенства (33), убеждаемся, что Ух е Кп: ||х|| < £2 справедливо равенство

гапкЯ3 = п. (34)

Из равенств (33), (34) следует справедливость условия (31) в области ||х|| < £2. Кроме того, из структуры матрицы Яз следует, что в указанной области имеет место оценка

П»-1! = 0(епат), т ^то. (35)

Используя (34), выполним замену переменной г(4п-1) по формуле

г(4п-1) = Я2(т )у. (36)

В результате получим систему

Согласно [15] матрицы в правой части системы (37) имеют вид

Я^1 ^Р82 - = {ё2,... ,ёк1,(рк1(т),... ,ёк1+...+кг_1+2,... ,ёп,(ркг(т)},

»-1 Р = {§1,..., ..., ё7+1}, 7 = I] 1=1 /¿, п = ^ Г=1 /¿.

Здесь ё; = (-^,..., -^1,..., - <,..., , 0,..., 0)ПХ 1; ё = (0,..., 1,..., 0)ПХ1, где единица стоит на ьм месте; являются коэффициентами разложения вектора

3+1 по векторам М,..., ,..., ц,..., ьк., з =l,..., r, т.е.

Ч + 1 = - Е (т )Ь1-----£ Л; (т )Ь3. (38)

¿=1 ¿=1

Рассмотрим задачу стабилизации системы

= 1 = 1,...,г, (39)

Уk¿ = (у^ , ...,У^) "Т х 1, = (0,..., 1,..., 0)т:х1, где единица стоит на i-м месте,

/1 к*\ Т

<ё к* = ^ к* ) к ¿х1.

Пусть уk¿ = Фазовые переменные системы (39) связаны с функцией ^(т) и ее производными равенствами

..., у 1. = к*-1) + Г кг-2(т-2) + • • • + Г1 (т)^(1) + го(т

Если продифференцировать последнее равенство из (40), то из первого уравнения системы (39) получим

к<) + £ к4-1(т к-1) + ••• + £о(т = ¿¿, i = 1,...,г. (41)

В (40), (41) функции гд^_2(т), ..., го(т), -i(r), ..., £о(т) являются линейными комбинациями функций (т), i = 1,..., kj, j = 1,..., r, и их производных.

Замечание 3. Из структуры матриц P и Q (см. (22) и представление (38)) следует, что функции ^ (т), ..., (т), их производные, а также функции rki-2(т), ..., г0(т), eki-1(r), ..., £0(т) ограничены. Остальные элементы столбцов , j = 1,..., r, удовлетворяют оценке O(e(n-1)aT), т ^ то. Пусть

k

di = Х>*-j(т) - Yki-j)0(ki-j), i = 1,...,r, (42)

j=i

где ya^ -j, j = 1,..., ki, выбраны так, чтобы корни Л^., ..., А^ уравнений

Aki + Yfci-iAfci-1 + • • • + Yo =0, i = 1, • • •, r, удовлетворяли условиям

Afc, = A^ , если i = j, A^ < — (2n + 1)a - 1, j = 1,..., ki, i = 1,...,r. (43)

В силу (40), (43) и замечания 3 закон управления di = S^T-1yki, i = 1,...,r, обеспечивает экспоненциальное убывание решений системы (37). Возвращаясь в (42) к исходным переменным с учетом (36), получим

di = Т-^1 S-fc1 c(4n-1), i = 1,...,r, (44)

Ski = (£ki-1(т) — Yfci-1,..., £о(т) — Yo), Tfci — матрица равенства (40), т.е. yfci = Tfci0, 0 = (0(ki-1),..., 0)T, S-1 — матрица, состоящая из соответствующих строк матрицы S-1. Найденное управление можно записать в виде (29), где М(т) = SaT-1S2fc1 =

= (Ski ^^kl , . . . , Skr T-r1 ^kt )T .

Обозначим через Ф(т), Ф(0) = E (E — единичная матрица) фундаментальную матрицу системы (41), замкнутую управлением (44). Очевидно, что элементами матрицы Ф(т) являются экспоненты с отрицательными показателями и их производные. Рассмотрим систему (28), замкнутую управлением (44):

dc(4n-1)

---= D(r)c(4™-1), D(r) = P(r) + Q(r)M(r). (45)

ат

Введем в рассмотрение блочно-диагональную матрицу T(t), на диагонали которой стоят матрицы Tki, i = 1,..., r. Тогда на основании (36), (40) фундаментальная матрица Ф(т), Ф(т)-1 = E системы (45) имеет вид

Ф(т) = S2(t )T(t )Ф(т )T-1 (0)S-1(0). (46)

Из (46), структуры матриц S2(t), Ф(т) и замечания 3 следует оценка

||Ф(т)|| < Ke-aTе-Лт, Л > 0, K> 0. (47)

Положим £ = £2. Тогда из (47) следует справедливость утверждения леммы. □

Кроме того, на основании (35), (44), (46) и замечания 3 имеем

||Ф(т)Ф-1(*)|| < К1в-Л(т-i)e(n-1)ai, t < т, т е [0, то), K1 > 0,

||М(т)|| = O(enaT), т ^то. ( )

5. Доказательство теоремы. Система (21), замкнутая управлением (44), имеет

вид

аг(4п-1)

—-= К1(с(4п-1),М(г)С(4п-1),Х,Г)+К2(С(4п-1),М(Г)С(4п-1),Х,Г) +

ат

+ Б(т)г(4п-1) + Rз(г(4n-1), М(т)г(4п-1),т) + R4(г(4n-1), М(т)г(4п-1), х,т), (49) Выполним в системе (49) замену переменной

г(4п-1) = ze-3naт, г(4п-1)(0) = z(0). (50)

В результате получим систему

— = С(т)х + е3пат111 (е-3пат2, М(т)е-Зпатг, х, г) + ат

+ е3патR2 (е-3патz, М(т)е-3патz, X, т) + е3патRз (е-3патz, М(т)е-3патz, т) +

+ е3патR4(e-3nатz, М(т)е-3патz, х,т), С(т) = Б(т) + 3паЕ. (51)

Покажем, что все решения системы (51) с начальными данными (50), начинающиеся в достаточно малой окрестности нуля, экспоненциально убывают.

Пусть Ф1(т), Ф1 (0) = Е — фундаментальная матрица системы. Тогда на основании (47), (48), (50) получим

||Ф1(т)|| < Ке-вт, ||Ф1(т)Ф-1 (г)|| < К1е-в(т-4)е(п-1)а4, в = А - 3па. (52)

Выберем а так, чтобы было выполнено

в > 0. (53)

Решение системы (51) с начальными данными (23) при т е [т1, то) принимает вид

1 „3па( .

z(т) = Ф1 (т)Ф-1(тl)z(тl) + Гф^т)Ф-1 (г)е3пах

т1

Rl (е-3па^, М(г)е-3па^, X, г) + (е-3паЧ М(г)е-3паЧ X, г) + + Rз(e-3nаíz, М(г)е-3па^, г) + R4(e-3nаíz, М(г)е-3паz, X, г) ¿г, (54) а при т е [0, т1] —

z(т) = Ф1 (т )г(4п-1)(0) + / 1 Ф1(т)Ф-1 (г)е3па4х

Jо

х |н.1(е-3паЧМ(г)е-3паЧх,г) + R2(е-3паг^,М(г)е-3паЧх,г)+

+ Rз(e-3naíz, М(г)е-3паЧ г) + R4(e-3naíz, М(г)е-3паz, X, г) ¿г. (55)

В свою очередь из (54), (55) с учетом (25)-(27), (29), (50) и (52) в области (5) (при замене г на z), (16) имеем оценки:

|^(т)|| < Ке-втЦФ^ЫЦ^ЫЦ + Г К1е-в(т-4)(1е-^|^(г)|| + Ь4е-^||х||)аг (56)

т1

при т € [т1, то) и

|^(т)|| < Ке-вт||с(4"-1)(0)|| + [Т1 К1е-в(т-4)(Ье-а4|^)|| + Ь4в-а4||х||)^ (57)

Jo

при т € [0,т1]. Константа Ь зависит от области (5), (16).

Применяя к неравенствам (56), (57) известный результат [16], получим

в(т)|| < Ке-мт||Ф-1(т1)|||Ит1)|| + [ К1в-^(т-4)Ь4в-а4ухуО^, т € [ть то),

в(т)|| < Ке-^1т||с(4"-1)(0)|| + / 1 К1в-^1(т-4)Ь4в-а4ухуа4, т € [0,т1],

Jo

(58)

где ^ = в - К1Ье-ат1, = в - К1Ь.

Используя условие (53), зафиксируем т1 > 0 так, чтобы было выполнено неравенство ^ > 0. Далее, ограничим выбор а > 0 условием а < Тогда после интегрирования вторых слагаемых правых частей оценок (58) имеем неравенства

|^(т)|| < Ке-мт|Ф-1(т1 )||Кт1)|| + К2е-атЬ4||х||, т € [т1, то), К > 0, |^(т)|| < Кз||с(4"-1)(0)|| + К4Ь4|х|, т € [0, т 1], Кз,К > 0.

На основании (23), (27) в области ||х|| < е последние две оценки можно записать в виде одного неравенства

|^(т)|| < К5е-ат||х||, т € [0, то), К5 > 0. (59)

В (59) константа К5 зависит от области (5), (16) и константы т1. Воспользовавшись (29), (48), (50), (59), получим оценку ||а(т)||:

||а(т)|| < ||М(т)||е-3"атК5е-ат||х|| < Кев-(2"-1)ат||х||, т € [0, то), Кб > 0. (60)

Положим величину е > 0 в условии теоремы равной е = шт{е, С1/К5, Ж/Кб}. В результате показано, что решение системы (51) с начальными данными (50), (23) при ||х|| < е не покинет область ||z| < С1, экспоненциально убывает и при этом согласно (60) соответствующая функция ^т) удовлетворяет ограничению (5). Кроме того, если это решение подставить в формулы (50), (29) и вернуться к переменным с(т) по формулам (20), (17), то будем иметь решение задачи (12), (13). Далее, возвращаясь к исходным зависимым и независимым переменным по формулам (14), (11), (7) и переходя к пределу при £ ^ 1 согласно замечанию 1, получим решение исходной задачи (1), (6). Теорема доказана. □

Алгоритм решения поставленной задачи состоит из следующих этапов.

1. Построение вспомогательной системы (21). Выполняется аналитическими методами.

2. Решение задачи стабилизации системы (28) аналитическими методами. В результате находится управление вида (29).

3. Решение задачи Коши для вспомогательной системы (21), замкнутой обратной связью (29). Выполняется численными методами.

4. Возврат к исходным зависимым и независимым переменным. Выполняется аналитическими методами.

6. Пример. Решение задачи перевода материальной точки в центральном поле тяготения. В качестве иллюстрации эффективности предложенного метода рассмотрим задачу перевода материальной точки, движущейся по круговой орбите с постоянной скоростью в центральном поле тяготения с помощью реактивной силы. Согласно [12] система уравнений в отклонениях относительно указанного движения по круговой орбите и условия (6) имеют вид

у , (хо +^з)2 , , , ,

Х1=Х2, х2 = -7-^7 + 7-+ аги, х3 = (г0 + х\)а,фи, (Ы )

(го + Х1 )2 (го + Х1 )3

х,(0)=0, х,(1) = X,, i = 1, 2, 3. где Х1 = г - го, Х2 = Г, хз = х - хо, Хо = (^го)1/2, м = т/т, го

(62)

радиус круговой орбиты, г — радиальная скорость; х — обобщенный импульс, аг, а^, —проекции вектора относительной скорости отделяющейся частицы на направление радиуса и поперечное направление соответственно (постоянные величины), т и т —соответственно масса и скорость изменения массы, V = ^оМ, ^о — постоянная всемирного тяготения; М — масса Земли, х = (х , х2,хз)т, м е К1. Система (12) и условия (13) примут вид

а°1 -ат а°2 — ат / / N. Л а°3 -ат/ , - , N ,

ае (С1 + Ж1 + г0)а,ф(1, (63)

ае~атс2, ^ = ае ат(5г(сь с3) + аг£), ат ат ат

д(с1, сз) = --

+

(Хо + сз + хз)2

(го + 01 + Х1 )2 (го + 01 + Х1 )3 ' с,(0) = -X,, с,(т) ^ 0 при т ^ то, i = 1, 2, 3. (64)

Всюду в дальнейшем будем предполагать, что X = (х, 0, 0)т. Для решения задачи (63), (64) необходимо выполнить семь преобразований типа (20):

с2 = г2-е атд(х1), = + ^е 2атд(х1),

1

г2 =-ш2 - -е 6

3ат

дс1

(х1 )^(ж1),

= + w1 = и1 + е-6ат ё,

"4ат w2=v2-e-^cJ

дс1

— 7ат =

и = М2 - е ё,

(65)

Положим аг = 1. Матрицы Р, Р и система, которая является аналогом системы (28), примут вид

аёг

— = Рс + дс1, с = (г;ьм2,с3)т, (66)

р

ат

010

а21 0 а23 000

0

ае-атЪ+ %е-3ата,ф)

где а21 = а23 = Ъ = (г0 + х^Лф.

_ 9д

V

ат

Для решения задачи стабилизации системы (66) строим матрицу Я = {Ь1, Ь2, Ьз},

= Ц, и = - г = 2,3.

После преобразования с = Яу, у = (у1,у2,уз)т получим систему

= <Р1{т)уз + <1, = У1 + <Р2{т)уз, = У2 + <Рз(т)уз- (67)

Замена уз = ф(т) приводит систему (67) к линейному уравнению третьего порядка

ф^(т) - Мт)Ф{2)(т) - ^ + Мт)) Ф(1)(т) - + ^ + Ыт)) Ф(т) = <1.

Т (68) Переменные у1, у2, Уз связаны с переменными ф(2), ф(1), ф равенством

у = ТФ, Т

1 ~<Рз +

-<£>з

0 0 1

, Ф = (ф(2),ф(1),ф)т

Пусть ¿=-(6 + ^з )ф(2) - (И + 2^ + " (б + После

подстановки ¿(т) в (68) получим уравнение, корни характеристического уравнения которого суть -1, -2, -3. Возвращаясь к исходным переменным, имеем

й = Г(т)Т-\т)8-\т)с, Г = - + П + + 6 + + ^ + ^

V ат ат2 ат

(69)

Очевидно, управление (69) обеспечивает экспоненциальное убывание решений системы (66). На заключительном этапе решаем задачу Коши для системы, полученной из системы (63) после замены фазовых координат по формулам (65) и замкнутой управлением (69) в новых переменных, и осуществляем переход к исходным переменным. При этом начальные данные для задачи Коши имеют вид

<п\ - З^1) 1 д9 I- \ I- \ = уЦО) = -хх-----24:'дс^ ~ 9'

1 дд

м2(0) = д(х 1) + - — (х1)д(х1) + д + д, с3(0) = 0. 6 дс1

В процессе численного моделирования интегрировалась вспомогательная система, полученная из системы (63), (69) после замены фазовых координат по формулам (65) с начальными данными и>1(0), «2(0), сз(0) при с = 10 м, г о = 7 • 10б м, а = 0,1, а^, = 0,1 на промежутке [0; 12,5]. Далее был выполнен переход от вспомогательных переменных к исходным. На рисунке представлены графики управляющей функции и соответствующих ей функций изменения фазовых координат Х1(£), Х2(£), хз(£) в исходной независимой переменной

На основании результатов численного моделирования можно сделать следующие выводы.

1. Наибольшие энергетические затраты для управляющего воздействия приходятся на начальный период перевода и напрямую зависят от значения конечного состояния Х1.

а

б

x3(t)

0.6

0.8

t

-600J v

Управляющая функция (а); функции изменения фазовых координат (б).

2. Максимальное отклонение конечного состояния по координате, допустимое для реализации алгоритма, при заданном времени перехода составляет 100 м.

3. Численное моделирование рассмотренного примера легко реализуется на персональной ЭВМ средних возможностей.

Литература

1. Зубов В. И. Лекции по теории управления. 2-е изд. М.: Лань, 2009.

2. Комаров В. А. Синтез ограниченных управлений для линейных неавтономных систем // АиТ. 1984. № 10. С. 44-50.

3. Крищенко А. П. Исследование управляемости и множеств достижимости нелинейных систем управления // АиТ. 1984. №6. С. 30-36.

4. Aeyels D. Controllability for polynomial systems // Lect. Notes Contr. and Inf. Sci. 1984. Vol.63. P. 542-545.

5. Qin H. On the controllability of nonlinear control system // Comput. & Maths. with Appls. 1985. Vol. 10, N6. P. 441-451.

6. Пантелеев В. П. Об управляемости нестационарных линейных систем // Дифференц. уравнения. 1985. Т. 21, №4. С. 623-628.

7. Коробов В. И. Почти полная управляемость линейных стационарных систем // Укр. мат. журн. 1986. Т. 38, №2. С. 163-169.

8. Айсагалиев С. А. К теории управляемости линейных систем // АиТ. 1991. №2. С. 35-41.

9. Бердышев Ю. И. О построении области достижимости в одной нелинейной задаче // Изв. РАН. Теория и системы упр. 2006. №4. С. 22-26.

10. Квитко А. Н. Об одном методе решения граничной задачи для нелинейной нестационарной управляемой системы с учетом результатов измерений // АиТ. 2012. №12. С. 89-110.

11. Kvitko A. N., Taran T. S., Firyulina O. S. Control problem with incomplete information // Proc. IEEE Int. Conf. "Stability and Control Processes" in Memory of V. I. Zubov (SCP). Saint-Petersburg. 2015. P. 106-109.

12. Красовский Н. Н. Теория управления движением. М.: Наука, 1968.

13. Кабанов С. А. Оптимизация динамики систем с коррекцией параметров структуры управления // Вестник СПбГУ. Сер. 1. 2014. Т. 1(59), вып. 2. С. 254-260.

14. He W., GeS.S. Vibration control of a flexible beam with output constraint // IEEE T. Ind. Electron. 2015. Vol.62, N8. P. 5023-5030.

15. Смирнов Е. Я. Стабилизация программных движений. СПб.: Изд-во С.-Петерб. ун-та, 1997.

16. Ба,рбашин Е. А. Введение в теорию устойчивости. М.: Наука, 1967.

Статья поступила в редакцию 6 апреля 2017 г.; рекомендована в печать 22 июня 2017 г.

Сведения об авторах

Квитко Александр Николаевич — доктор физико-математических наук, профессор; a.kvitko@spbu.ru Фирюлина Оксана Сергеевна — кандидат физико-математических наук; o.firulina@spbu.ru Еремин Алексей Сергеевич — кандидат физико-математических наук, доцент; a.eremin@spbu.ru

AN ALGORITHM OF SOLUTION OF A BOUNDARY VALUE PROBLEM FOR A NONLINEAR STATIONARY CONTROL SYSTEM IT MODELLING

Alexander N. Kvitko, Oksana S. Firyulina, Alexey S. Eremin

St. Petersburg State University, Universitetskaya nab., 7—9, St. Petersburg, 199034, Russian Federation; a.kvitko@spbu.ru, o.firulina@spbu.ru, a.eremin@spbu.ru

An algorithm for a control function construction is proposed. It is applicable for a wide class of stationary nonlinear systems of ordinary differential equations. The control transfers system to a given final state and obeys certain restrictions. The algorithm is easy to implement numerically. A sufficient condition of the desired transfer function existance is constructed. A certain practical problem is considered and simulated numerically with the presented method. Refs 16. Fig. 1.

Keywords: boundary condition, stabilization, control, nonlinear system.

References

1. Zubov V. I., Lectures on control theory (2th ed., Lan' Publ., Moscow, 2009) [in Russian].

2. Komarov V. A., "Design of constrained control signals for nonlinear non-autonomous systems", Automat. Rem. Contr. 45(10), 1280-1286 (1984).

3. Krishchenko A. P., "Controllability and attainability sets of nonlinear control systems", Automat. Rem. Contr. 45(6), 707-713 (1984).

4. Aeyels D., "Controllability for polynomial systems", Lect. Notes Contr. and Inf. Sci. 63, 542-545

(1984).

5. Qin H., "On the controllability of nonlinear control system", Comput. & Maths. with Appls. 10(6), 441-451 (1985).

6. Panteleev V. P., "Controllability of time-dependent linear systems", Diff. Equat. 21(4), 623-628

(1985) [in Russian].

7. Korobov V. I., "Almost full controllability of linear stationary systems", Ukranian Math. J. 38(2), 163-169 (1986) [in Russian].

8. Aisagaliev S. A., "On the theory of the controllability of linear systems", Automat. Rem. Contr. 52(2/1), 163-171 (1991).

9. Berdyshev Yu.I., "On the construction of the reachability domain in one nonlinear problem", J. Comput. Sys. Sci. Int. 45(4), 526-531 (2006).

10. Kvitko A. N., "On one method of solving a boundary problem for a nonlinear nonstationary controllable system taking measurement results into account", Automat. Rem. Contr. 73(12), 2021-2037 (2012).

11. Kvitko A. N., Taran T.S., Firyulina O.S., "Control problem with incomplete information", Proc. IEEE Int. Conf. "SStability and Control Processes" in Memory of V. I. Zubov (SCP), Saint-Petersburg, 106-109 (2015).

12. Krasovsky N.N., Theory of motion control (Nauka Publ., Moscow, 1968) [in Russian].

13. Kabanov S. A. "Dynamic system optimization using correction of control structure parameters," Vestnik SPbSU. Series 1, 1(59), issue2, 254-260 (2014) [in Russian].

14. He W., Ge S. S. "Vibration control of a flexible beam with output constraint", IEEE T. Ind. Electron. 62(8), 5023-5030 (2015).

15. Smirnov E. Ya., Stabilization of programmed motion (CRC Press, 2000).

16. Barbashin E.A., Introduction to the theory of stability (Wolters-Noordhoff, 1970).

Для цитирования: Квитко А. Н., Фирюлина О. С., Еремин А. С. Aлгоритм решения краевой задачи для нелинейной управляемой системы и его численное моделирование // Вестник СПбГУ. Математика. Механика. Астрономия. 2017. Т. 4(62). Вып. 4. С. 608-621. https://doi.org/10.21638/11701/spbu01.2017.409

For citation: Kvitko A.N., Firyulina O.S., Eremin A. S. An algorithm of solution of a boundary value problem for a nonlinear stationary control system it modelling. Vestnik SPbSU. Mathematics. Mechanics. Astronomy, 2017, vol. 4(62), issue 4, pp. 608-621. https://doi.org/10.21638/11701/spbu01.2017.409