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

Научная статья
DOI:
https://doi.org/10.23670/IRJ.2017.62.066
Выпуск: № 8 (62), 2017
Опубликована:
2017/08/18
PDF

Антоновская О.Г.

ORCID: 0000-0002-5688-7996, Кандидат физ.-мат. наук, Нижегородский государственный архитектурно-строительный университет

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

Аннотация

В настоящей статье предлагается методика определения момента окончания переходного  процесса при изменении параметров непрерывной динамической системы на основе применения прямого метода Ляпунова, основанная на оценивании областей притяжения состояний равновесия. За оценку области, которую траектория системы не покинет с течением времени, принимается окрестность состояния равновесия, в которой первая производная функции Ляпунова отрицательна. Функция Ляпунова выбирается в виде положительно определенной квадратичной формы, для которой знакоотрицательность ее первой производной в силу линеаризованной системы обеспечивается с заданным запасом.

Ключевые слова: непрерывная динамическая система, макроструктура пространства состояний, метод функций Ляпунова.

Antonovskaya O.G.

ORCID: 0000-0002-5688-7996, PhD in Physics and Mathematics, Nizhny Novgorod State University of Architecture and Civil Engineering

APPLICATION OF LYAPUNOV’S QUADRATIC FUNCTIONS IN SOLVING APPLIED DYNAMIC PROBLEMS

Abstract

In the following paper, we propose a technique for determining the moment of the end of the transient process when the parameters of a continuous dynamic system are changed on the basis of the application of the direct Lyapunov method, based on the estimation of the regions of attraction of equilibrium states. For an estimate of the region the trajectory of the system does not leave with time, we take a neighborhood of the equilibrium state where the first derivative of the Lyapunov function is negative. The Lyapunov function is chosen in the form of a positive definite quadratic form, where the negative sign of its first derivative is provided with a given margin due to a linear system.

Keywords: continuous dynamic system, macrostructure of the state space, the Lyapunov function method.

При исследовании динамики систем одной из наиболее важных является проблема определения длительности переходных процессов в системе. В ходе решения задачи определения длительности переходных процессов в существенно нелинейных системах при переключении параметров в заданном диапазоне возникает задача точного определения момента окончания переходного процесса в системе [1], [2], [3]. Дело в том, что критерием его окончания может оказаться попадание траектории на заданное, быть может, неограниченное множество, содержащее состояние равновесия, соответствующее рабочему режиму системы [1], [3]. И тот факт, что траектория попала на заданное множество, не гарантирует того, что она в дальнейшем этого множества не покинет. Таким образом, возникает задача получения оценок областей притяжения, т.е. таких областей, которые траектория с течением времени не покинет. Для решения этой задачи может быть использован второй (прямой) метод Ляпунова или метод функций Ляпунова, когда изучается поведение некоторой вспомогательной функции вдоль траекторий системы [4], [5], [6].

Центральное место во втором (прямом) методе Ляпунова занимает проблема построения функции Ляпунова. При этом функцию Ляпунова нелинейной динамической системы часто ищут в классе положительно определенных квадратичных форм, исходя из того условия, что построенная квадратичная форма является функцией Ляпунова для соответствующей линеаризованной системы [5, С. 120-132]. Кроме того, может ставиться вопрос о построении квадратичных функций Ляпунова с некоторыми заданными свойствами, которые определяются особенностями задачи. В частности, в динамических задачах, когда интерес представляют собой ухе не только качественные, но и количественные характеристики системы, возникает необходимость ограничений на первую производную функции Ляпунова вдоль траекторий системы  [7,8].

В настоящей работе предлагается способ получения аналитических оценок окрестности асимптотически устойчивого состояния равновесия непрерывной динамической системы, которую траектория системы с течением времени не покинет, основанный на  применении прямого метода Ляпунова и использующий методику работ [9,10].

Пусть нелинейная непрерывная динамическая система задана системой дифференциальных уравнений

30-08-2017 10-12-07  (1)

имеющей асимптотически устойчивое состояние равновесие 30-08-2017 10-13-12. И пусть функции fi имеют непрерывные частные производные до второго порядка включительно.  Запишем (1) в виде [4, C. 74-75]

30-08-2017 10-14-37   (2)

где 30-08-2017 10-15-16 - постоянные, равные значениям частных производных функций fi при 30-08-2017 10-15-57 – остаточные члены второго порядка в разложении по формуле Тейлора

30-08-2017 10-16-45  (3)

В силу непрерывности вторых производных для остаточных членов (3) на любом ограниченном множестве, содержащем состояние равновесия, можно указать оценку

30-08-2017 10-17-43  (4)

где 30-08-2017 10-18-23

Выберем функцию Ляпунова в классе положительно определенных квадратичных форм 30-08-2017 10-19-14   (5)

так, чтобы она была функцией Ляпунова соответствующей (2) линеаризованной системы

30-08-2017 10-22-28   (6)

и удовлетворяла при этом условию [6,7]  30-08-2017 10-22-43, где

30-08-2017 10-24-56  (7)

первая производная (5) в силу линеаризованной системы, а 30-08-2017 10-25-52 – собственные числа матрицы системы 30-08-2017 10-27-03. В [8] показано, что выбор коэффициентов (5), удовлетворяющей этому условию, можно осуществить с помощью простых явных соотношений.

Оценим первую производную такой функции Ляпунова для нелинейной системы (2) подобно [4].

Следует отметить, что, согласно [10], в точках поверхности уровня 30-08-2017 10-27-52имеют место неравенства

30-08-2017 10-28-28  (8)

где Aii – алгебраическое дополнение элемента Kii в матрице , а det K  – определитель матрицы K. Но тогда, согласно (4),

30-08-2017 10-31-23   (9)

а первая производная (5) в силу нелинейной системы (2) будет иметь вид 30-08-2017 10-32-13   (10)

где 30-08-2017 10-33-03 есть первая производная (5) в силу уравнений линейной системы. Так как на сечении 30-08-2017 10-33-25, то, учитывая ограничения (8), (9), получаем, что

30-08-2017 10-35-05  (11)

Значит, неравенство 30-08-2017 10-35-52 будет обязательно выполняться, если 30-08-2017 10-44-26   (12) поскольку 30-08-2017 10-45-11. То есть при 30-08-2017 11-02-38  (13)

первая производная в силу нелинейной системы будет отрицательной. Таким образом, доказана следующая теорема.

Теорема. Область

30-08-2017 11-03-20  (14)

где

30-08-2017 11-04-08  (15)

является оценкой снизу области притяжения асимптотически устойчивого состояния равновесия 30-08-2017 11-04-57  системы дифференциальных уравнений (2).

Замечание 1. Оценка (14) может быть использования при определении точного момента окончания длительности переходного процесса в динамической системе, математическая модель которой представляет собой нелинейную систему дифференциальных уравнений вида (2),  при численном исследовании поведения ее траекторий. Пусть критерием окончания переходного процесса является попадание траектории системы (2) на заданное множество 30-08-2017 11-05-50, содержащее устойчивое состояние равновесия системы.  В случае, если множество (14)-(15) 30-08-2017 11-06-28, его можно принять за окрестность состояния равновесия, которую траектория системы не покинет с течение времени. Если это не так,  рассмотрим  множество

30-08-2017 11-07-11   (16)

 для  поверхности уровня 30-08-2017 11-07-58, вписанной в 30-08-2017 11-05-50. Например, в случае, когда  есть «σ -полоса»

30-08-2017 11-09-41  (17)

значение 30-08-2017 11-20-50  (18) Скажем, при n=2 30-08-2017 11-21-52   (19)

Так как в этом случае 30-08-2017 11-22-38, то  на множестве (16), и его можно принять за искомую окрестность состояния равновесия. То есть в качестве окрестности состояния равновесии (2), которую траектория системы не покинет с течением времени, можно взять множество

30-08-2017 11-22-55  (20)

где 30-08-2017 11-25-27, а момент окончания переходного процесса в системе определять следующим образом.

За начальную точку траектории выберем точку, координаты которой соответствуют координатам состояния равновесия при значении параметра, с которого происходит переключение, по отношению к координатам состояния равновесия 30-08-2017 11-31-55, соответствующего значению параметра, на которое переключение происходит. Вычисления будем проводить до тех пор, пока траектория системы не попадет на множество 30-08-2017 11-05-50, содержащее 30-08-2017 11-31-55. Запомнив момент попадания траектории на 30-08-2017 11-05-50, продолжим вычисления до тех пор, пока либо траектория не попадет в область (19), либо не покинет множество 30-08-2017 11-05-50. В первом случае будем считать моментом окончания переходного процесса запомненный ранее момент попадания в 30-08-2017 11-05-50, во втором – продолжим вычисления до следующего момента попадания траектории на 30-08-2017 11-05-50, с последующим повторением всей процедуры.

Замечание 2. Доказанная теорема может быть использована для получения оценок области притяжения состояния равновесия нелинейной системы (2) [9]. Пусть мы ищем оценку области притяжения состояния равновесия 30-08-2017 11-31-55 системы (2), принадлежащую ограниченному множеству 30-08-2017 11-36-25. Поскольку поверхности уровня 30-08-2017 11-36-39, вписанной в 30-08-2017 11-36-25, соответствует значение 30-08-2017 11-36-51 [9], то искомой аналитической оценкой области притяжения является множество

30-08-2017 11-39-17  (21)

Замечание 3. 1. Пусть все корни характеристического уравнения, соответствующего состоянию равновесия 30-08-2017 11-31-55, действительны и различны, причем 30-08-2017 11-46-34.  В этом случае всегда существует линейное невырожденное преобразование координат, приводящее систему (10) к каноническому виду (столбцы матрицы преобразования координат являются собственными векторами матрицы A, соответствующими собственным значениям 30-08-2017 11-47-25 соответственно).

Рассмотрим  систему дифференциальных уравнений вида

30-08-2017 11-47-59   (22)

Следует отметить, что согласно [8], квадратичную функцию Ляпунова для системы (22) можно выбрать в виде

30-08-2017 11-50-16 Тогда, в силу доказанной теоремы, получим, что в этом случае 30-08-2017 11-54-04  (25)
  1. Пусть все корни характеристического уравнения, соответствующего состоянию равновесия 30-08-2017 11-31-55, таковы, что 30-08-2017 11-55-49, т.е. действительны и различны. И пусть 30-08-2017 11-55-59. В этом случае всегда существует линейное невырожденное преобразование координат, приводящее линеаризованную систему к каноническому виду [4].

Рассмотрим систему дифференциальных уравнений

30-08-2017 11-57-09

Заметим [7], что в этом случае квадратичную функцию Ляпунова, удовлетворяющую условию  30-08-2017 11-58-13,  можно выбрать в виде (23), где

30-08-2017 11-59-15

Тогда, в силу доказанной теоремы, в этом случае 30-08-2017 11-59-57.

Список литературы / References

  1. Горюнов В. И. Техническая полоса захвата одноконтурного синтезатора частоты / В. И. Горюнов, В. Н. Ерусланов, Н. И. Лобашов // Техника средств связи. Сер. Техника радиосвязи. – 1990. –  Вып. 2. – С. 88-94.
  2. Горюнов В. И. К анализу быстродействия синтезатора при переключениях по диапазону / В. И. Горюнов, В. Н. Ерусланов, М. Н. Зайцева // Техника средств связи. Сер. Техника радиосвязи. – 1986. – Вып. 3. – С. 44-49.
  3. Антоновская О. Г. К анализу формы и длительности переходных процессов при переключениях синтезатора с делителем частоты и пропорционально-интегрирующим фильтром по диапазону / О. Г. Антоновская, В. И. Горюнов, Н. И. Лобашов // Динамика систем: Межвуз. сб. Горький: Изд-во ГГУ, 1989. –  С.59-72.
  4. Ляпунов А. М. Общая задача об устойчивости движения. / А. М. Ляпунов.. – М.-Л.:  Изд-во техн.-теор. лит., 1950. –  472 с.
  5. Барбашин Е. А. Функции Ляпунова. / Е. А. Барбашин. –  М. Наука, 1970. – 240 с.
  6. Четаев Н. Г. Устойчивость движения. / Н. Г. Четаев. –   М.: Наука, 1965. –  207 с.
  7. Антоновская О. Г. О построении квадратичной функции Ляпунова с заданными свойствами / О. Г. Антоновская // Дифференциальные уравнения. – 2013. –    Т.49. –   № 9. –   С. 1220-1224.
  8. Антоновская О. Г. О построении квадратичной функции Ляпунова с заданными свойствами / О. Г. Антоновская // Дифференциальные уравнения. – 2016. –   Т.52. –    № 3. –   С. 276-281.
  9. Антоновская О. Г. Об одном способе оценки области притяжения устойчивого решения системы дифференциальных уравнений / О. Г. Антоновская // Международный научный журнал «Инновационная наука». – 2015. –  №9. –  С. 11-15.
  10. Антоновская О. Г. Об одном способе оценки размеров области притяжения неподвижной точки нелинейного точечного отображения произвольной размерности / О. Г. Антоновская, В. И. Горюнов // Известия вузов. Математика. –  2016. –  № 12. –  С. 12-18.

Список литературы на английском языке / References in English

  1. Goryunov V. I. Tehnicheskaya polosa zahvata odnokonturnogo sintezatora chastoty [Technical capture strip of one-contour frequency synthesizer] / V .I. Goryunov, V. N. Eruslanov, N. I. Lobashov // Tehnika sredstv svyasi. Ser. Tehnika radiosvyasi [Techniques of communication. Series Techniques of radio-communication]. – 1990. Vyp. 2. P. 88–94. [in Russian]
  2. Goryunov V.I. K analizu bystrodeystviya sintezatora pri pereklucheniyah po diapazonu [On synthesizer speed analysis when switching in a range] / V .I. Goryunov, V. N. Eruslanov, M. N. Zaytseva // Tehnika sredstv svyasi. Ser. Tehnika radiosvyasi [Techniques of communication. Series Techniques of radio-communication]. – 1986. Vyp. 3. P. 44–49. [in Russian]
  3. Antonovskaya O. G. K analizu formy I dlitelnosti perehodnyh processov pri pereklyucheniyah sintezatora s delitelem chastity I proportsionalno-integrtiruyushim filtrom po diapazonu [On the analysis of form and duration of transient processes when switching over a range in the synthesizer with counter and proportional-integrating filter] / O. G. Antonovskaya, V. I. Goryunov, N. I. Lobashov // Dinamika system: Mezhvuz. sbornik [System dynamics: Inter-university Collection]. – Gorky: Izd-vo GGU, 1989. P. 59–72. [in Russian]
  4. Lyapunov A. M. Obshchaya zadacha ob ustoychivosti dvizheniya [General problem of the stability of motion]. / A. M. Lyapunov. – Moscow–Leningrad : Izdat. Nthn.-teor. Lit., 1950. 472 p.
  5. Barbashin E.A. Funktsii Lyapunova [Lyapunov functions]. /E. A. Barbashin. – Moscow: Nauka, 1970. 240 p.
  6. Chataev N. G. Ustoychivost’ dvizheniya [Stability of motion]. / N. G. Chetaev. – Moscow: Nauka, 1966. 207 p.
  7. Antonovskaya O. G. On the construction of quadratic Lyapunov function with given properties / O. G. Antonovskaya // Differential equations. – 2013. – V. 49. –  № 9. –   P. 1187-1191.
  8. Antonovskaya O. G. Determination of coefficients of a quadratic Lyapunov function with given properties / O. G. Antonovskaya // Differential equations. – 2016. – V. 52. – № 3. –   P. 275-261.
  9. Antonovskaya O. G. Ob odnom sposobe otsenki oblasti prityazheniya ustoychivogo resheniya sistemy differentsialnyh uravneniy [On a method of evaluation of attraction domain for fixed point of the stable solution of differential system] / O. G. Antonovskaya // Mezhdunarodnyy zhurnal ”Innovatsionnaya nauka” [International magazine ”Innovation Science”].– 2015. – № 9. – P. 11-15. [in Russian]
  10. Antonovskaya O. G. On a method of evaluation of attraction domain for fixed point of nonlinear point mapping of arbitrary dimension / O. G. Antonovskaya, V. I. Goryunov // Russian Mathematics.– 2016. – V. 60. – № 12. – P. 10-14.