Автоматика и телемеханика, № 5, 2023
Нелинейные системы
© 2023 г. В.Н. ТХАЙ, д-р физ.-мат. наук (tkhai@ipu.ru)
(Институт проблем управления им. В.А. Трапезникова РАН, Москва)
СТАБИЛИЗАЦИЯ КОЛЕБАНИЙ УПРАВЛЯЕМОЙ
АВТОНОМНОЙ СИСТЕМЫ
Рассматривается гладкая автономная система общего вида, допускаю-
щая невырожденное периодическое решение. Строится глобальное семей-
ство (по параметру h) невырожденных периодических решений, выво-
дится закон монотонного изменения периода на семействе, доказывается
существование редуцированной системы второго порядка. Для нее реша-
ется задача стабилизации колебания управляемой системы, выделенного
значением параметра h. Находится гладкое автономное управление, кон-
струируется притягивающий цикл.
Ключевые слова: автономная система, невырожденное периодическое ре-
шение, глобальное семейство, теорема Ляпунова о центре, управление,
притягивающий цикл, естественная стабилизация.
DOI: 10.31857/S0005231023050033, EDN: AFVGEL
1. Введение
В 1927 г. Ван дер Поль предложил уравнение, описывающее линейный
осциллятор, на который действует малая нелинейная сила, линейная по ско-
рости, и создающая диссипацию в каждой текущей точке траектории (дис-
сипация Ван дер Поля). Уравнение допускает притягивающий цикл. Впо-
следствии Л.С. Понтрягин в [1] нашел достаточные условия для выделения
цикла из семейства периодических решений гамильтоновой системы: приме-
нялись негамильтоновые возмущения. В линейном осцилляторе колебания
изохронные, в нелинейной гамильтоновой системе период колебаний зависит
от постоянной энергии.
В [2] диссипация типа Ван дер Поля вводилась для семейства периоди-
ческих решений, в котором период монотонно зависит от параметра семей-
ства h: для линейного осциллятора получается диссипация Ван дер Поля. По-
строена система, в которой для стабилизации колебания с параметром h вы-
бирается соответствующее значение параметра в управлении. Управляемая
система остается автономной и обладает асимптотически орбитально устой-
чивым циклом. Эти результаты развивались в [3-5], мехатронная система
стабилизация предложена в [6]. В многомерных системах строилась редуци-
рованная система возможного низшего порядка.
29
Другие исследования по стабилизации желаемого режима колебания от-
личаются применением управлений, явно зависящих от времени. Приведем
некоторые из них. Обзор на примере перевернутого маятника дается в [7].
В [8] предлагается раскачивающее управление (swinging control). В [9] ре-
шается задача об орбитальной стабилизации периодических решений мало-
приводных нелинейных систем (с числом независимых приводов на едини-
цу меньше числа степеней свободы неуправляемой консервативной системы).
Синтезированный нелинейный закон управления с обратной связью зависит
от времени. Стабилизация желаемых уровней механической энергии посред-
ством импульсного управления проводится в [10], робастное стабилизирую-
щее управление колебаниями ищется неявным методом Ляпунова в [11].
Редуцированную систему можно выделить, когда описано все множество
периодических решений. В нелинейных системах период колебаний, как пра-
вило, меняется от колебания к колебанию. Множество одночастотных колеба-
ний с монотонным изменением периода на нем от параметра, принимающего
все возможные для множества значения, становится единицей такого мно-
жества и называется глобальным семейством периодических решений. В фа-
зовом пространстве глобальное семейство заполняет связное инвариантное
множество. В системе может быть конечное или счетное число глобальных
семейств.
Задача о глобальном семействе первоначально возникла в многочисленных
приложениях теоремы А.М. Ляпунова о центре (1892 г.). В теореме устанав-
ливается существование локального семейства нелинейных периодических
движений, примыкающей к равновесию консервативной системы. В конкрет-
ных задачах всегда встает вопрос о границе применимости ляпуновского се-
мейства. Прогресс в вопросе о границе начался с исследований А.А. Зевина
(1997 г.). В [12] даны условия на гамильтониан, гарантирующие в компак-
те Ω ∈ R2n продолжение ляпуновского семейства на границу ∂Ω. Результаты
развиты в [13]: найдены строго звездообразные (strictly starshaped) гамиль-
тонианы, для которых условия продолжения [12] выполняются.
Сейчас стало понятно, что вопрос о границе в теореме Ляпунова разреша-
ется знанием о глобальном семействе периодических решений, которое вклю-
чает, как составляющее, локальное ляпуновское семейство.
Обратимая механическая система выделяется свойством пространственно-
временной симметрии, однако в общем случае не допускает первый интеграл.
В ней рассматриваются симметричные периодические движения. Глобальные
семейства в этих системах изучались (см. [4, 5, 14, 15]); понятие глобального
семейства вводилось в [14]. Результаты по глобальному семейству (включая
вопросы его существования, построения, свойства и т.д.) использовались при
нахождения управления с меняющим значение параметром в [4-6].
В данной статье для автономной системы общего вида, не стесненной
дополнительными ограничениями (гамильтоновость, обратимость, консерва-
тивность и т.д.), ставится и решается задача построения глобального семей-
30
ства невырожденных периодических решений. При этом применяется подход
Понтрягина в теории обыкновенных дифференциальных уравнений (см. [16])
введения непродолжаемого решения. Во второй части статьи для редукциро-
ванной системы второго порядка строится управляемая система. Находится
гладкое автономное управление с настраиваемым параметром, выделяющим
стабилизируемое колебание.
2. Невырожденное периодическое решение автономной системы
Рассматривается гладкое автономное уравнение
(1)
Ż = Z(z), z ∈ Rn,
общее решение которого обозначается через z(z01, . . . , z0n, t), где z0 = (z01,
...,z0n) — начальная точка (при t = 0). Необходимое и достаточное условие
существования T -периодического решения системы (1) записывается в виде
равенства
(2)
f ≡ z(z01,...,z0n,T) - z0
= 0.
Пусть уравнение (2) имеет решение z0 = z∗, T = T∗, не совпадающее с рав-
новесием1: Z(z∗) = 0. В силу автономности системы (1) уравнение (2) вместе с
указанным решением всегда обладает семейством решений по параметру γ —
сдвига точки z0 вдоль траектории:
(3)
z0 = z∗(γ), T = T∗.
Вычисляется ранг Ra функциональной матрицы Af (матрицы Якоби) для
функции f с параметром T в точке z∗ при T = T∗. Тогда получается: Ra ≤
≤ n - 1.
Определение 1. Случай Ra = n - 1 называется невырожденным для
периодического решения. При Ra = n - 1 периодическое решение называется
невырожденным, в противном случае, оно называется вырожденным.
Далее исследуется случай Ra = n - 1.
Уравнение (2) может допускать единственное решение вида (3). Тогда ав-
тономное уравнение (1) имеет изолированное периодическое решение — цикл
с периодом T∗. Альтернативой будет семейство решений, в котором период T
меняется от решения к решению, т.е. является функцией некоторого пара-
метра h: z0 = z0(γ, h), T = T (h).
Из равенства (2) в окрестности решения (3) выводится система линейных
равенств
∂fs
∂fs
∂fs
ξs ≡
dz01 + . . . +
dz0n +
dT = 0, s = 1, . . . , n,
(4)
∂z01
∂z0n
∂T
f = (f1,...,fn),
1 замечание рецензента
31
в которой частные производные вычисляются при z0 = z∗, T = T∗. Систе-
ма (4), как и равенство (2), удовлетворяется тождественно по сдвигу γ. Так
как производная
∂fs
∂z(z0,T)
=
= Zs(z0),
∂T
∂T
то это необходимо приводит к линейной зависимости функций Zs(z0), s =
= 1, . . . , n.
В ситуации цикла система (4) выполняется только при ΔT = T - T∗ = 0.
Матрица Af имеет простое нулевое собственное значение.
Далее рассматривается случай, когда матрица Af содержит жорданову
2-клетку из нулевых собственных значений. Здесь имеем скалярный пара-
метр h (см. [17]). Из равенств (4), расматриваемых как система линейных
уравнений с матрицей Af , выводится существование в уравнении (2) локаль-
ного семейства решений по dT . Соответственно, уравнение (1) допускает ло-
кальное h-семейство периодических решений. Для T∗-периодического реше-
ния: T∗ = T (h∗), T′(h∗) = 0 (через (′) обозначается дифференцирование по h).
Параметр h принадлежит окрестности числа h∗.
Далее h называется параметром семейства периодических решений. Функ-
ция T (h) — монотонна, поэтому за h можно выбирать сам период T .
В ситуации семейства решения уравнений в вариациях состоят из частных
производных от периодического решения соответственно по t и h. В резуль-
тате дифференцирования по h получается растущее вместе со временем t
решение, которое выписывается, например, в [18, стр. 416, формула (9.9)].
Вводится определение 2.
Определение 2. Семейство периодических решений уравнения (1) на-
зывается невырожденным, если период T (h) на нем монотонно зависит от
параметра h.
Определение 2 справедливо для рассмотренного локального семейства, в
котором h меняется в окрестности числа h∗. Оно остается справедливым для
любого семейства периодических решений: на нем h меняется в интервале.
Значит, существует интервал изменения параметра h, который отвечает всем
периодическим решениям семейства.
Определение 3. Невырожденное семейство периодических решений,
на котором параметр h принимает всевозможные для решений семейства
значения, называется глобальным семейством.
В фазовом пространстве глобальное семейство представляется связным
множеством точек. Интервал изменения h в глобальном семействе может
быть конечным (энергия маятника от нижнего до верхнего равновесий) или
неограниченным (неограниченные по координате колебания). В уравнении (1)
может существовать одно, нескольно или счетное множество глобальных се-
мейств невырожденных периодических решений. Предметом статьи будет от-
дельное глобальное семейство.
32
Для обозначения невырожденного периодического решения, принадлежа-
щего семейству, далее применяется сокрашение НПР. Тогда глобальное се-
мейство невырожденных периодических решений будет называться глобаль-
ным семейством НПР и обозначается через Σ.
НПР существует для системы размерности n ≥ 2. Колебания математи-
ческого маятника дают пример НПР: x + sin x = 0, n = 2, Ra = 1. При этом
локальные колебания близ нижнего равновесия образуют ляпуновское семей-
ство НПР. Его продолжение по периоду приводит ко всему семейству колеба-
ний маятника (глобальное семейство НПР), связывающему нижнее и верхнее
равновесия. Период колебаний на нем, начиная со значения 2π, монотонно,
вместе с амплитудой колебания, стремящейся к π, стремится к бесконечности.
В линейном осцилляторе реализуется семейство вырожденных (изохрон-
ных) колебаний одного периода: n = 2, Ra = 0.
Замечание 1. В обратимых механических системах НПР могут быть
симметричными, для таких решений понятие НПР уточняется в [4].
3. Локальные свойства НПР
НПР обладает следующими свойствами.
1◦. Свойство локальной продолжимости НПР.
В системе (4) приращение dT не равно нулю и меняется независимо от γ.
При непрерывном двустороннем изменении приращения dT из системы (4)
находятся все кривые семейства НПР. Поэтому уравнение (1) вместе с НПР,
где h = h∗, обладает семейством НПР с монотонной зависимостью перио-
да T (h) от h. При этом параметр h принадлежит некоторой окрестности
числа h∗. Приращение dT ограничивается условием Ra = n - 1.
Семейство НПР также можно выделять интервалом изменения периода T ;
семейство обозначается через Σ(T ).
2◦. Семейство Σ(T ) заполняет двумерную областьΣ(T ).
ОбластьΣ(T ) образована замкнутыми кривыми — НПР уравнения (1),
которые параметризованы сдвигом γ. Эти кривые меняются вместе с пара-
метром семейства h.
Замечание 2. На локальном семействе период является единственным
параметром. При этом T′(h∗) = 0 и период зависит от одного параметра h
(закон зависимости периода нелинейных колебаний от одного параметра [17]).
4. Глобальное семейство НПР
Глобальное семейство НПР (Σ) получается двусторонним продолжением
семейства Σ(T ) по периоду T (по параметру h). Соответственно, продолжение
областиΣ(T ) по периоду T приводит в фазовом пространстве к глобальной
областиΣ, заполненной решениями из Σ.
33
Лемма 1. Глобальное семейство Σ существует и заполняет глобальную
областьΣ. На нем период является монотонной функцией параметра се-
мейства. Для точек областиΣ ранг Ra = n - 1, на ее границе ∂Σ условие
Ra = n - 1 не выполняется.
Доказательство. Рассматривается НПР системы (1). Согласно свой-
ствам 1◦ и 2◦ оно принадлежит семейству Σ(T ) и заполняет областьΣ(T ).
К каждому НПР семейства Σ(T ) применяется свойство 1◦. В результате воз-
никают две ситуации. В первой из них областьΣ(T ) не меняется: глобальное
семейство Σ — построено. Во-второй ситуации областьΣ(T ) расширяется в
том смысле, что содержит все точки прежней области, а также новые точки.
Значит, семейство Σ(T ) получило продолжение по периоду T , как в сторону
увеличения, так и уменьшения T . Процесс продолжения описывается систе-
мой (4). На полученном семействе монотонное изменение периода от пара-
метра семейства сохраняется.
На следующем шаге возникает такая же альтернатива, как на предыдущем
шаге итерации.
В конечном результате получается глобальная область Σ, в которой выпол-
няется условие Ra = n - 1. ОбластьΣ заполняется глобальным семейством
НПР. На границе ∂Σ условие Ra = n - 1 нарушается.
Построение Σ происходит за бесконечное число шагов.
Замечание 3. Для достижения глобального семейства применяется под-
ход Понтрягина к построению непродолжаемого решения в теории обыкно-
венных дифференциальных уравнений.
Лемма 2. Глобальное семейство НПР описывается редуцированной си-
стемой второго порядка.
Доказательство. Равенства (4) выполняются на глобальном семей-
стве Σ тождественно по паре (z0, T ). Ранг матрицы Af равен n - 1, произ-
водные по T линейно зависимы, а вектор df/dT = 0. Поэтому n - 1 линейно
независимых дифференциальных форм из (4), связанных условиями (4), при-
водятся к системе, в которой только одна из полученных форм содержит dT .
Независимость оставшихся форм от dT означает, что решения глобального
семейства описываются системой второго порядка.
Замечание 4. В случае консервативной системы получается редуциро-
ванная консервативная система с одной степенью свободы.
На основе лемм 1 и 2 формулируется теорема 1.
Теорема 1. Пусть уравнение (1) допускает НПР. Тогда оно продолжа-
ется по периоду T на глобальное семейство Σ. На Σ период T(h) монотонно
зависит от параметра семейства h. Семейство Σ заполняет глобальную
областьΣ; Σ описывается редуцированной системой второго порядка. Для
точек областиΣ ранг Ra = n - 1, на ее границе ∂Σ условие Ra = n - 1 не
выполняется.
34
Замечание 5. При подходе к границе ∂Σ производная T′(h) может стре-
мится к нулю, бесконечности или перестает существовать. В случае T′(h) → 0
границей семейства может быть вырожденное периодическое решение. Или
на границе находится центр (примеры 1 и 2, κ < 1). Также семейство может
стать неограниченным по координате (пример 2, κ > 1). Случай T′(h) → ∞
реализуется для центра (пример 2, κ > 1), седла (пример 1), а также — для
неограниченного по координате глобального семейства (пример 2, κ < 1) (см.
также [19]).
Следствие 1. Ляпуновское семейство продолжается на глобальное се-
мейство невырожденных периодических решений (глобальная теорема Ля-
пунова о центре).
Доказательство. Период на локальном ляпуновском семействе моно-
тонно зависит от постоянной энергии (см. [18]). Поэтому семейство состоит
из НПР. По теореме 1 глобальное семейство существует и получается про-
должением любого НПР. Следовательно, теорема Ляпунова о центре имеет
глобальный характер.
Замечание 6. Проблема границы в теореме Ляпунова о центре получает
полное решение описанием глобальной областиΣ и ее границы.
Замечание 7. Из теоремы 1 в частном случае выводятся результаты
[4, 14, 15].
Теорема 1 и замечание 5 иллюстрируются примерами. В частности, на них
показывается, какие движения реализуются на ∂Σ.
Пример 1. x + sinx = 0. В математическом маятнике колебания обра-
зуют глобальное семейство невырожденных колебаний, границами которой
служат равновесия (центр и седло) и сепаратриссы. На семействе период мо-
нотонно возрастает от центра к седлу.
Пример 2. x + xκ = 0,κ > 0. На фазовой плоскости равновесие x = 0,
x = 0 окружено семейством периодических решений. Уравнение допускает
первый интеграл
x2/2 + xκ+1/(κ + 1) = h(const).
Вычислим период колебания с амплитудой xmax:
∫
dx
T =4
√
2(h - xκ+1/κ + 1)
0
Положим x = h1/(κ+1)y. Тогда получается
∫
1-κ
dy
T = 4h2(κ+1)
√
,
2(1 - yκ+1/(κ + 1))
0
где xmax = h1/(κ+1)ymax. Отсюда следует, что при κ < 1 период монотонно воз-
растает вместе с h (вместе с амплитудой колебаний), при κ = 1 (случай линей-
35
ного осциллятора) период колебаний не зависит h, а при κ > 1 период моно-
тонно убывает к нулю с возрастанием амплитуды колебаний. Следовательно,
в случае κ = 1 уравнение обладает семейством невырожденных неограничен-
ных колебаний, одной из границ служит нулевое равновесие, а другая граница
уходит на бесконечность.
5. Управляемая редуцированная система
Пусть уравнение (1) допускает НПР. По теореме 1 оно принадлежит гло-
бальному семейству Σ, которое заполняет глобальную областьΣ. Предпола-
гается, чтоΣ является притягивающей.
Семейство Σ описывается редуцированной системой второго порядка. Для
этой системы ставится задача стабилизации НПР из семейства Σ. Рассмат-
ривается гладкая система.
Таким образом, решается задача нахождения в системе
(5)
x = X(x,y) + εF,
y = Y (x,y) + εG,
гладких управлений F и G с малым коэффициентом усиления регулятора ε
таких, чтобы в (5) реализовался притягивающий цикл, близкий к выбранному
НПР семейства Σ. Предполагается, что семейство Σ, т.е. решения системы (5)
при ε = 0, описываются формулами
(6)
x = ϕ(h,t), y = ψ(h,t).
Находятся управления, решающие задачу стабилизации для любого НПР се-
мейства Σ: для выбранного НПР значение параметра h = h∗.
Для малых ε записывается неоднородная линейная система
δ x = a11δx + a12δy + εF,
(7)
δ y = a21δx + a22δy + εG,
δx = x - ϕ(h,t), δy = y - ψ(h,t),
где однородная часть совпадает с уравнениями в вариациях для решения (6).
В (7) учитываются только линейные члены по δx, δy: добавление нелинейных
членов не меняет получаемых качественных выводов, которые количественно
будут различаться от получаемых на величины o(ε).
Лемма 3. Система (7) приводится к виду
u=ε
ψ(h∗, t)F - ϕ˙(h∗, t)G]/Δ,
(8)
T′(h∗)
v=
u + ε[η(h∗,t)F - ξ(h∗,t)G]/Δ,
T∗
где T∗ — периодические функции ξ(h∗, t) и η(h∗, t) — производные от функций
соответственно ϕ(h,(T/T∗)t) и ψ(h,(T/T∗)t) по h при h = h∗.
36
Доказательство. Преобразованием Ляпунова однородная часть в (7)
приводится к системе с постоянными коэффициентами. Тогда вся преобразо-
ванная система записывается как (8).
Замечание 8. В явном виде преобразование дается в Приложении. Фор-
мулы (8) справедливы для НПР, отвечающего значению параметра h = h∗.
6. Притягивающий цикл управляемой системы
Из первого уравнения системы (8) выводится необходимое условие
∫
(9)
ψ(h∗, t)F -ϕ˙(h∗, t)G]/Δ(h∗
,t)dτ = 0
0
существования T∗-периодического решения управляемой системы. В случае
выполнения (9) интегрированием в (8) второго уравнения находится формула
для переменной v. При этом единственность изолированного T∗-периодиче-
ского решения, т.е. существование цикла выводится из анализа (9).
Условие (9) можно рассматривать как уравнение для определения значе-
ния параметра h∗ для НПР. Поэтому заменой h∗ в функциях ϕ и ψ на h из (9)
получается амплитудное уравнение
∫T∗
(10)
I(h) ≡
ψ
(h, t)F - ϕ˙(h, t)G]/Δ(h, t)dτ = 0.
0
Система (8) задает отображение t : 0 → T на периоде T∗. Неравенством
I′(h∗) = 0 гарантируется единственность неподвижной точки отображения,
а неравенством I′(h∗) < 0 — принадлежность собственного числа левой полу-
плоскости. Оно вычисляется как в [2].
Каждому простому корню h = h∗ уравнения (10) соответствует единствен-
ное T∗-периодическое решение уравнения (8), т.е. цикл управляемой систе-
мы (5). Соответствующее условие притяжения решений управляемой системы
к циклу дается неравенством I′(h∗) < 0.
Таким образом, справедлива лемма 4.
Лемма 4. Простому корню амплитудного уравнения
(10) отвечает
цикл управляемой системы (5). При выполнении неравенства I′(h∗) < 0 цикл
будет притягивающим.
7. Поиск управления
Основой для поиска управления служит амплитудное уравнение (10). Со-
гласно лемме 4 задача стабилизации колебания управляемой системы (5) ре-
шается управлениями F и G, для которых амплитудное уравнение (10) до-
пускает простой корень h∗: I′(h∗) < 0.
37
Переменные x и y и действующие управления F и G входят в систему (5)
равноправно. Равноправность сохраняется в амплитудном уравнении (10).
Получается, что выделенное колебание можно стабилизировать любой парой
(F, 0) или (0, G), а также в общем случае выбирать F = 0, G = 0. С другой
стороны, для периодического решения переменные x и y задают положение и
скорость точки на траектории. Это хорошо демонстрируется в приведенной
системе (8), в которой переменной u задается положение точки, а v — ее ско-
рость. Также видно, что при выборе пары (F, 0) диссипация по переменной x
приводит также к диссипации по y.
В соответствии с постановкой задачи в разделе 5 находится управление,
стабилизирующее любое НПР семейства Σ, выделяемое значением парамет-
ра h. Значит, управление содержит некоторую характеристику K такую, что
при подстановке в выражение для управления вместо числа K функции K(h)
амплитудное уравнение (10) будет выполняться тождественно, т.е.
∫
Φ(K(h), ϕ(h, t), ψ(h, t),ϕ˙(h, t),ψ(h, t))dτ ≡ 0,
(11)
0
Φ = ( ψF -ϕ˙G)/Δ.
Без ограничения общности функция Φ предполагается линейной по K. Тогда
тождество (11) выполняется с функцией
(12)
Φ = [1 - Ka(h,t)]b(h,t).
Явные выражения для a(h, t) и b(h, t) получаются при подстановке функций
ϕ(h, t) и ψ(h, t) вместе с их производными в Φ.
При дифференцировании тождества (11) по h с подставленной в него
функцией (12) получается
∫
dT (h)
dΦ(h)
(13)
[1 - Ka(h, T )]b(h, T )
+
dτ ≡ 0.
dh
dh
0
На НПР одна из производных ϕ˙(h, t) илиψ(h, t) в моменты времени t = 0, T
равняется нулю. Поэтому управления ищутся при условии b(h, 0) = b(h, T ) =
= 0, и при h = h∗ из (13) выводится равенство
∫
dI(h∗)
dK(h∗)
(14)
=
ν, ν = a(h∗,t)b(h∗
,t)dτ.
dh
dh
0
38
Само число K для заданного h∗ вычисляется из амплитудного уравнения
(11):
∫
b(h∗, t)dτ
0
(15)
K(h∗) =
∫
a(h∗, t)b(h∗, t)dτ
0
Зависимость K(h) является характеристикой управления (и одновременно
семейства Σ).
Таким образом, на основе леммы 4 формулируется общая теорема 2.
Теорема 2. Любое НПР глобального семейства Σ, соответствующее
значению параметра h = h∗, стабилизируется управлениями с функци-
ей (12), если для характеристики K(h) при h = h∗ выполняется неравенство
K′(h∗)ν < 0.
Замечание 9. Применяется адаптивное управление с настраиваемым па-
раметром h.
8. Некоторые управления
Теоремой 2 выделяется класс управлений. В частном случае механиче-
ских систем в этот класс входит управление (см. [2-4]), которое естествен-
ным образом, без применения иных управлений, стабилизирует колебание.
В нем применяется нелинейная сила, линейная по скорости, действующая в
текущей точке траектории и приводящая к диссипации. Таким же свойством
обладают и другие управления найденного класса.
Далее приводятся некоторые конкретные управления.
1. Рассматриваются управления
(16)
F = α(1 - Kx2)˙x, G = -β(1 - Ky2
)y,
α, β - const,
c характеристикой K = K(h). Тогда из (11) выводится амплитудное уравне-
ние
∫T∗
I(h) ≡
[1 - K(αx2 + βy2)]ϕ˙(h, t)ψ(h, t)/Δ(h, t)dτ = 0.
0
Отсюда находится функция
∫
ϕ(h, t)ψ(h, t)/Δ(h, t)dτ
0
K(h) =
∫
[αϕ2(h, t) + βψ2(h, t)]ϕ˙(h, t)ψ(h, t)/Δ(h, t)dτ
0
39
Наконец, формула (14) для стабилизации НПР со значением параметра h =
= h∗ приобретает вид
∫
dI(h∗)
dK(h∗)
=
ν, ν = (αx2 + βy2)ϕ˙(h∗,t)ψ(h∗,t)/Δ(h∗,t)dτ.
dh
dh
0
Замечание 10. Для одного уравнения второго порядка x = y, поэтому
в (16) полагается: α = 1, β = 0.
Замечание 11. Результаты для механической системы, подверженной
действию позиционных сил, в частности, потенциальных сил, известны
(см. [2]). Для системы Δ(h, t) = 1.
Замечание 12. В случае Y (x,0) ≡ 0 получается редуцированная обрати-
мая механическая система второго порядка 2. Для системы Δ(h, t) = 1. К ней
в [4] применяются управления (16).
2. Выбираются управления
F = (1 - Kx2)y, G ≡ 0.
Тогда характеристика дается формулой
∫
ψ2(h,t)/Δ(h,t)dτ
0
K(h) =
∫
ϕ2(h,t)ψ2(h,t)/Δ(h,t)dτ
0
В рассматриваемом случае из (14) получается
∫
dI(h∗)
dK(h∗)
=
ν, ν = ϕ2(h∗,t)ψ2(h∗,t)/Δ(h∗,t)dτ.
dh
dh
0
3. Пусть применяются управления
F ≡ 0, G = -(1 - Ky2)˙x.
Для них вычисляется такая характеристика
∫
ϕ2(h,t)/Δ(h,t)dτ
0
K(h) =
,
∫
ψ2(h,t)ϕ˙2(h,t)/Δ(h,t)dτ
0
а формула (14) записывается в виде
∫
dI(h∗)
dK(h∗)
=
ν, ν = ψ2(h∗,t)ϕ˙2(h∗,t)/Δ(h∗,t)dτ.
dh
dh
0
40
Замечание 13. Примеры стабилизации НПР конкретных систем приве-
дены в [2-4, 20]. Далее дается пример стабилизации колебания, принадлежа-
щего вырожденному семейству.
9. Вырожденное семейство
На семействе НПР параметр h отражает монотонное изменение периода T
вместе с h: на семействе T′(h) = 0. На вырожденном семействе неравенство
не удовлетворяется. Для такого семейства теорема 1 о глобальном семействе
не применима. Тем не менее развитый в статье подход применим к вырож-
денному семейству на двумерном многообразии.
Для вырожденного семейства вместо периода можно предложить другую
характеристику семейства решений, а именно — амплитуду колебания. Так-
же, как и период, амплитуда может монотонно меняться на семействе. Реше-
ния семейства параметризуются амплитудой колебания.
Далее приводится примечательный пример стабилизации (в большом)
цикла, рожденного из вырожденного семейства периодических решений ав-
тономной системы.
Пример 3. Рассматривается управляемая система
(17)
x=y,
y + x = ε(1 - Kx2
)y,
где K — постоянная. В ней к линейному осциллятору применяется управле-
ние (16), в котором с точностью до обозначений переменных берутся числа:
α = 1, β = 0.
При ε = 0 колебания даются формулой x = A cos(t+γ), где γ — временной
сдвиг на траектории. Записывается амплитудное уравнение
2π
∫
(1 - KA2 cos2(t + γ))A2 sin2(t + γ)dτ = 0.
0
√
Отсюда находятся два корня: A = 0 и A = 2/
K. Первый корень соответству-
ет началу координат, второй корень — циклу, близкому к колебанию линей-
√
ного осциллятора с амплитудой 2/
K. Амплитуда A выбирается в качестве
параметра семейства колебаний линейного осциллятора, при этом амплитуда
цикла системы (17) задается параметром K в управлении.
Зависимость K(A) дается формулой K = 4/A2, производная K′(A) < 0,
поэтому цикл стабилизируется в малом.
Таким образом, подход к стабилизации колебания, принадлежащего семей-
ству НПР, развивается на семейство вырожденных периодических решений.
Замечание 14. В системе (17) по существу записано уравнение Ван дер
Поля. Цикл системы (17) стабилизируется глобально.
41
10. Заключение
Невырожденное периодическое решение автономной системы может быть
циклом или принадлежать семейству. На НПР период монотонно меняется
вместе с параметром семейства. НПР продолжается по периоду на глобальное
семейство НПР. Глобальное семейство НПР заполняет инвариантное много-
образие и описывается редуцированной системой второго порядка. На гло-
бальном семействе монотонное изменение периода сохраняется.
Задача стабилизации НПР решается для глобального семейства НПР,
заполняющего глобальную область, которая предполагается притягивающей.
Развивается подход, в котором управление строится для всего глобального се-
мейства НПР, а для стабилизации выбранного НПР фиксируется параметр се-
мейства. Применяется гладкое автономное управление, в результате которого
реализуется притягивающий цикл, близкий к выделенному НПР глобального
семейства. Подход применим также к вырожденному семейству периодиче-
ских решений.
ПРИЛОЖЕНИЕ
Приведение однородной системы.
Уравнения в вариациях для решения (6) записываются в виде
∂X(x,y)
∂X(x,y)
δx=
δx +
δy,
∂x
∂y
(Π.1)
∂Y (x,y)
∂Y (x,y)
δy=
δx +
δy,
∂x
∂y
где частные производные вычисляются для x = ϕ(h, t), y = ψ(h, t).
Фундаментальная система решений в (Π.1) дается матрицей
∂ϕ(h, t)
∂ϕ(h, t)
∂t
∂h
(Π.2)
∂ψ(h,t)
∂ψ(h,t)
∂t
∂h
c определителем
∫t
)
(∂X(x,y)
∂Y (x,y)
Δ(h, t) = Δ(h, 0) exp
+
dτ .
∂x
∂y
0
Записывается общее решение системы в (Π.1)
δx = c1ϕ˙(h,t) + c2ϕ′(h,t),
(Π.3)
δy = c1ψ˙(h,t) + c2ψ′(h,t)
42
с постоянными c1 и c2. Разрешение системы (Π.3) относительно c1 и c2 ис-
пользуется для перехода к новым переменным u, v:
(Π.4)
u = -(ψ˙
δx -ϕ˙δy)/Δ, v = [η(h,t)δx - ξ(h,t)δy]/Δ.
При этом полагается
T′(h)
∂ϕ(h, t)
ξ(h, t) =
t ϕ(h, t) +
,
T∗
∂h
T′(h)
η(h, t) =
t ψ(h,t) + ∂ψ(h,t),
T∗
∂h
где функции ξ(h∗, t), η(h∗, t) будут T∗-периодическими. Наконец, из выраже-
ния для v в (Π.4) получается
T′(h)
(Π.5)
v=
tu + (ψ′(h,t)δx - ϕ′
(h, t)δy)/Δ.
T∗
Подстановкой в формулы (Π.4) и (Π.5) первого решения из матрицы (Π.2)
вычисляется: u = 0, v = 1. Соответственно для второго решения записыва-
ется:
u = -( ψϕ′ -ϕ˙ψ′)/Δ = 1, v = T′(h∗)tu + u.
T∗
Таким образом, система (Π.1) приводится к виду
T′(h∗)
u = 0,
v=
u.
T∗
СПИСОК ЛИТЕРАТУРЫ
1. Понтрягин Л.С. О динамических системах, близких к гамильтоновым // Журн.
эксперим. и теорет. физики. 1934. Т. 4. Вып. 9. С. 883-885.
2. Тхай В.Н. Стабилизация колебания управляемой механической системы // АиТ.
2019. № 11. С. 83-92.
Tkhai V.N. Stabilizing the Oscillations of a Controlled Mechanical System // Autom.
Remote Control. 2019. V. 80. No. 11. P. 1996-2004.
3. Тхай В.Н. Стабилизация колебания управляемой механической системы с N
степенями свободы // АиТ. 2020. № 9. С. 93-104.
Tkhai V.N. Stabilizing the Oscillations of an N Degree of Freedom Controlled Mecha-
nical System // Autom. Remote Control. 2020. V. 81. No. 9. P. 1637-1646.
4. Тхай В.Н. Стабилизация колебаний управляемой обратимой механической си-
стемы // АиТ. 2022. № 9. С. 94-108.
Tkhai V.N. Stabilization of Oscillations of a Controlled Reversible Mechanical
System // Autom. Remote Control. 2022. V. 83. No. 9.
5. Тхай В.Н. Режим цикла в связанной консервативной системе // АиТ. 2022. № 2.
С. 90-106.
Tkhai V.N. Cycle Mode in a Coupled Conservative System // Autom. Remote
Control. 2022. V. 83. No. 2. P. 237-251.
43
6.
Тхай В.Н. Мехатронная схема стабилизации колебаний // Изв. РАН. Теория и
системы управления. 2022. № 1. С. 9-16.
7.
Boubaker O. The Inverted Pendulum Benchmark in Nonlinear Control Theory: a Sur-
vey // Int. J. Adv. Robot. Syst. 2013. V. 10. No. 5. 233-242.
8.
Fradkov A.L Swinging Control of Nonlinear Oscillations // Int. J. Control.
1996.
V. 64. Iss. 6. P. 1189-1202.
9.
Shiriaev A., Perram J.W., Canudas-de-Wit C. Constructive Tool for Orbital Sta-
bilization of Underactuated Nonlinear Systems: Virtual Constraints Approach //
IEEE T. Automat. Contr. 2005. V. 50. No. 8. P. 1164-1176.
10.
Kant K., Mukherjee R., Khalil H. Stabilization of Energy Level Sets of Underac-
tuated Mechanical Systems Exploiting Impulsive Braking // Nonlinear Dynam.
2021. V. 106. P. 279-293.
11.
Guo Yu., Hou B., Xu Sh., Mei R., Wang Z., Huynh V.Th. Robust Stabilizing Con-
trol for Oscillatory Base Manipulators by Implicit Lyapunov Method // Nonlinear
Dynam. 2022. V. 108. P. 2245-226.
12.
Zevin A.A. Nonlocal generalization of Lyapunov theorem // Nonlinear Analysis,
Theory, Methods and Applications. 1997. V. 28. No. 9. P. 1499-1507.
13.
Zevin A.A. Global continuation of Lyapunov centre orbits in Hamiltonian systems //
Nonlinearity. 1999. V. 12. P. 1339-1349.
14.
Тхай В.Н. Колебания и равновесия в обратимой механической системе // Вест-
ник СПбГУ. Сер. 1. Матем. Механ. Астрон. 2021. Вып. 4. С. 709-715.
Tkhai V.N. Equilibria and oscillations in a reversible mechanical system // Vestnik
SPbSU. Mathematics. Mechanics. Astronomy. 2021. Vol. 54. No. 4. P. 447-451.
15.
Tkhai V.N. Spatial oscillations of a physical pendulum // Proc.
2022 16th Int.
Conf. on Stability and Oscillations of Nonlinear Control Systems (Pyatnitskiy’s
Conference), IEEE Xplore: 29 June 2022.
16.
Понтрягин Л.С. Обыкновенные дифференциальные уравнени. М.: Наука, 1974.
17.
Тхай В.Н. Закон о зависимости периода нелинейных колебаний от одного пара-
метра // Прикл. матем. механ. Т. 75. Вып. 3. C. 430-434.
18.
Малкин И.Г. Некоторые задачи теории нелинейных колебаний. М.: Гостехиздат,
1956.
19.
Devaney R.L. Blue Sky Catastrophes in Reversible and Hamiltonian Systems //
Indiana University Mathematics Journal. 1977. V. 26. No. 2. P. 247-263.
20.
Тхай В.Н. Стабилизация колебаний автономной системы // АиТ. 2016. № 6.
Tkhai V.N. Stabilizing the Oscillations of an Autonomous System // Autom. Remote
Control. 2016. V. 77. No. 6. P. 972-979.
Статья представлена к публикации членом редколлегии А.М. Красносельским.
Поступила в редакцию 25.10.2022
После доработки 20.01.2023
Принята к публикации 26.01.2023
44