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

Научная статья
DOI:
https://doi.org/10.23670/IRJ.2021.9.111.011
Выпуск: № 9 (111), 2021
Опубликована:
2021/09/17
PDF

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

Научная статья

Шашихин В.Н.1, *, Будник С.В.2

1 ORCID: 0000-0002-3718-9623;

1, 2 Санкт Петербургский политехнический университет Петра Великого, Санкт-Петербург, Россия

* Корреспондирующий автор (shashihin[at]bk.ru)

Аннотация

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

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

SUPPRESSION OF CHAOTIC FLUCTUATIONS IN NONLINEAR SYSTEMS WITH AFFINE STRUCTURE

Research article

Shashikhin V.N.1, *, Budnik S.V.2

1 ORCID: 0000-0002-3718-9623;

1, 2 Peter the Great St. Petersburg Polytechnic University, Saint Petersburg, Russia

* Corresponding author (shashihin[at]bk.ru)

Abstract

The current article examines the problem of reducing chaotic modes of nonlinear systems to regular modes. The stabilization of singular points and limit cycles is achieved by introducing feedback on phase coordinates and forming the corresponding spectrum of Lyapunov characteristic indicators of nonlinear systems with an affine control structure. The feedback coefficients are calculated using the modal control method generalized to nonlinear systems, based on the reduction of the system to the canonical Frobenius form. The proposed method of control synthesis is examined using a nonlinear Chua chain model as an example. The results of a computational experiment to determine the spectrum of Lyapunov characteristic indicators and construct a phase portrait of a nonlinear system confirm the operability of the proposed method for suppressing chaotic dynamics.

Keywords: control systems, stabilization, fluctuations.

Введение

Основная часть технических объектов и процессов окружающего нас мира имеет нелинейный характер поведения. Хорошо известно, что в нелинейных системах при определённых условиях возможно возникновение детерминированного хаоса, который может играть как положительную роль, так и отрицательную. В этой ситуации проблема управления хаосом приобретает большое теоретическое и прикладное значение [1].

Среди различных задач управления хаосом можно выделить задачу подавления, имеющую важную практическую роль при нейтрализации шумов в электронных системах и вибраций в различных конструкциях [2], [3]. Задача подавления хаоса заключается в формировании управляющего воздействия, при котором вместо хаотической динамики замкнутая нелинейная система будет иметь устойчивые особые точки или периодические решения [4], [5], [6], [7].

Для стабилизации хаотических систем исторически первыми и наиболее активно развиваемыми в настоящее время являются следующие методы: метод Отта–Гребоджи–Йорке (OGY-метод) [8], метод Пирагаса [9], метод Магницкого [10]. Также для управления хаотической динамикой используются традиционные подходы и методы автоматического управления. Например, в работе [11] стабилизация хаотической системы осуществляется с помощью обратной связи по отклонению, а для выбора коэффициентов регулятора используется критерий Рауса-Гурвица. Развивается подход и на основе аналитического конструирования агрегированных регуляторов [12], основанный на идеи введения совокупности инвариантных многообразий и поэтапной декомпозиции исходной динамической системы [13].

Настоящая работа посвящена синтезу управления нелинейной системой путем введения обратной связи. Искомое управление позволяет решать задачи, как подавления хаоса, так и его усиления за счет формирования требуемого спектра характеристических показателей Ляпунова.

Постановка задачи Пусть возмущенные движения нелинейного динамического объекта описываются векторным дифференциальным уравнением 02-11-2021 11-22-16   (1)

где 02-11-2021 11-22-30 ‑ вектор состояния; 02-11-2021 11-22-38 ‑ вектор управления; 02-11-2021 11-23-08 ‑ векторная функция, 02-11-2021 11-23-17 ‑ вещественные функции, определенные и непрерывные в области 02-11-2021 11-24-02 и имеющие в ней непрерывные частные производные, которые ограничены в замкнутой области 02-11-2021 11-24-25.

Вид траекторий системы (1) определяется характеристическими показателями Ляпунова. Одна из особенностей нерегулярных режимов –неустойчивость траекторий, принадлежащих хаотическому (странному) аттрактору. Количественной мерой этой неустойчивости являются характеристический показатель Ляпунова. Среди всего набора характеристических показателей Ляпунова наиболее важен наибольший (старший) показатель 02-11-2021 11-25-02. Характеристические показатели, упорядоченные по убыванию 02-11-2021 11-25-11, задают спектр Ляпунова нелинейной динамической системы.

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

02-11-2021 11-30-56  – состояние равновесия; (2а)
 02-11-2021 11-31-06– предельный цикл; (2б)
 02-11-2021 11-31-18– странный аттрактор, 02-11-2021 11-31-30. (2в)

Задача подавления хаоса состоит в преобразовании хаотического режима системы (1), который характеризуется спектром (2в), в регулярный режим со спектром (2а) или (2б), то есть обеспечить аттрактор в виде особой точки или предельного цикла.

Для решения этой задачи необходимо найти управление в виде обратной связи по фазовому вектору нелинейной системы (1)

02-11-2021 11-33-12     (3) которое обеспечит в замкнутой системе 02-11-2021 11-33-21    (4) спектр характеристических показателей Ляпунова 02-11-2021 11-33-34    (5) равный желаемому (требуемому) спектру 02-11-2021 11-33-57    (6)

Желаемый спектр (6) определяется требуемым характером регулярного движения системы (4).

Синтез обратной связи для линеаризованной системы

Пусть уравнение (1) описывает отклонения фазовых координат нелинейного объекта от некоторой траектории 29-10-2021 12-33-03, на которой он удерживается с помощью управляющего воздействия 29-10-2021 12-33-13. Используя формулу Тейлора в предположении дифференцируемости компонент функции 02-11-2021 11-35-32 в окрестности заданной траектории 29-10-2021 12-33-03, преобразуем уравнение (1) к квазилинейному виду

29-10-2021 12-34-15

с коэффициентом 29-10-2021 12-33-53

02-11-2021 11-36-29    (7)

Здесь 29-10-2021 12-33-53 - матрица Якоби нелинейной системы (1). Если матрица Якоби (7) не имеет собственных чисел на мнимой оси, тогда для синтеза управления, стабилизирующего нелинейную систему (1), в силу теоремы Гробмана-Хартмана [14], может быть использовано линеаризованное уравнение

29-10-2021 12-33-40   (8)

Желаемые собственные числа матрицы Якоби (7) можно задать с помощью методики синтеза модального регулятора (3) [15]. Модальный регулятор позволяет задать желаемый характер поведения линеаризованной системы за счет выбора требуемых собственных значений матрицы Якоби замкнутой системы. Пусть линеаризованная система (8) имеет вид

02-11-2021 11-36-54   (9)

где 02-11-2021 11-22-30 – вектор фазовых координат; 02-11-2021 11-38-42  – управляющее воздействие; 02-11-2021 11-38-50 - матрица (столбец) входа.

Требуется определить параметры 29-10-2021 12-39-46 линейного закона управления с обратной связью 29-10-2021 12-39-53 обеспечивающего заданные собственные числа 29-10-2021 12-40-04 матрицы замкнутой системы 29-10-2021 12-40-22 Из выражения для матрицы замкнутой системы 29-10-2021 12-40-37 нельзя непосредственно получить значения коэффициента обратной связи l, так как сама матрица 29-10-2021 12-41-09 неизвестна. Поэтому используется замена переменных 02-11-2021 11-40-24 такая, чтобы математическая модель преобразованной системы (9), замкнутой синтезируемым управлением,

02-11-2021 11-40-39    (10)

имела каноническую форму Фробениуса:

02-11-2021 11-42-01

где 29-10-2021 12-41-52 – коэффициенты характеристического полинома матриц  02-11-2021 11-42-08.

Для того чтобы привести матрицу системы к канонической форме, используется матрица T образованная из коэффициентов характеристического полинома матрицы объекта управления A следующим образом:

02-11-2021 11-42-14

Для преобразования 02-11-2021 11-43-11 используется матрица  где 29-10-2021 12-52-52 – матрица управляемости системы. Рассмотрим подробно такое преобразование:

02-11-2021 11-43-21

Если ввести обозначения

02-11-2021 11-43-28

тогда модель системы примет вид, определяемый выражением (10). Учитывая особенности построения матрицы Q, матрица системы имеет каноническую форму Фробениуса, а вектор b будет приведён к простейшему виду:

02-11-2021 11-47-31

где 29-10-2021 12-56-01 – матрица размера , имеющая на p-й «наддиагонали» единицы, а остальные нули; 29-10-2021 12-56-29 – единичный вектор размера n, n-я координата которого равна единице, а остальные – нулю; 29-10-2021 12-56-50 – вектор коэффициентов характеристического полинома матрицы A. После выполненных преобразований матрица замкнутой системы примет вид:

02-11-2021 11-48-12

Для того чтобы эта матрица имела требуемые собственные значения необходимо, чтобы коэффициенты ее характеристического уравнения соответствовали вектору 29-10-2021 12-57-25 – коэффициенты характеристического полинома матриц  02-11-2021 11-48-48. Тогда

02-11-2021 11-51-01

а коэффициенты регулятора определяются соотношением

29-10-2021 12-57-50  (11)

и обеспечивают желаемые собственные значения в замкнутой системы [16].

Синтез обратной связи для нелинейной системы

Рассмотрим методику синтеза регулятора для модели, которая в безразмерных переменных и коэффициентах имеет вид [17]:

02-11-2021 11-51-20    (12)

Будем использовать следующие значения параметров системы 02-11-2021 11-52-48, при которых характеристические показатели Ляпунова системы Чуа (12) равны 02-11-2021 11-52-59 и соответствуют спектру (2в), что свидетельствует о наличии хаоса и, следовательно, аттрактор системы - странный.

Матрица Якоби (7) нелинейной системы зависит от значений фазовых координат точки, в которой она вычисляется

02-11-2021 11-53-21    (13)

Матрица Якоби (13) при заданных параметрах системы зависит только от фазовой координаты 29-10-2021 13-05-43. Выберем точку 29-10-2021 13-05-49 на траектории неуправляемой системы (12), в которой необходимо изменить собственные числа матрицы Якоби. Для этого рассмотрим зависимость собственного числа матрицы Якоби 02-11-2021 11-58-24  на траектории системы (12), начинающейся в точке 29-10-2021 13-06-13 (см. рисунок 1). Из графика видно, что значение собственного числа симметрично относительно точки 29-10-2021 13-06-21 и эта точка является точкой локального экстремума функций 02-11-2021 11-58-38. В качестве точки 29-10-2021 13-05-49 выбираем точку, в которой 29-10-2021 13-06-48 – особая точка системы (12). В точке 29-10-2021 13-05-49 матрица Якоби равна 29-10-2021 13-06-55, а вектор её собственных значений равен.

02-11-2021 12-00-4029-10-2021 13-27-50

Рис. 1 – Зависимость действительной и мнимой части 02-11-2021 11-58-24

 

Исходную систему (12) необходимо изменть так, чтобы в точке 29-10-2021 13-05-49 матрица Якоби 29-10-2021 13-15-04 замкнутой системы имела заданные собственные значения: 02-11-2021 12-01-07 (a – параметр, близкий к единице)

29-10-2021 13-16-11

где  02-11-2021 12-01-31 – коэффициент обратной связи, который находится по формуле (11).

При управлении первой компонентой фазового вектора система (12), замкнутая модальным регулятором, принимает вид

02-11-2021 12-01-42   (14)

Зависимость показателя Ляпунова 02-11-2021 12-03-59 имеет вид (см. рисунок 2). Из графика зависимости 02-11-2021 12-03-59 выбирается 29-10-2021 13-16-55 так, чтобы либо сформировать предельный цикл, либо подавить хаос. Для формирования предельного цикла необходимо выбрать 29-10-2021 13-16-55 такой, что 02-11-2021 12-04-13. Тогда аттрактором замкнутой системы будет предельный цикл. Возьмем 29-10-2021 13-17-17, так как при нём второй показатель Ляпунова минимальный. Величине 29-10-2021 13-17-17 соответствует коэффициент обратной связи 29-10-2021 13-33-47 и спектр характеристических показателей Ляпунова (5) 02-11-2021 12-05-04. Последнее свидетельствует об устойчивости нелинейной системы (14) с синтезированным управлением.

График проекции траектории замкнутой системы с начальным условием 29-10-2021 13-34-05 на фазовую плоскость 29-10-2021 13-34-11 имеет вид (см. рисунок 3).

29-10-2021 13-42-38

Рис. 2 – Зависимость  от 02-11-2021 12-03-59

29-10-2021 13-42-57

Рис. 3 – Предельный цикл системы с обратной связью

 

Возьмём коэффициент 29-10-2021 13-16-55 такой, что  02-11-2021 12-10-24, тогда аттрактором замкнутой системы будет особая точка. Выберем параметр 29-10-2021 13-38-10. Ему соответствуют коэффициент обратной связи 29-10-2021 13-38-23 и спектр отрицательных характеристических показателей Ляпунова 02-11-2021 12-11-00 что свидетельствует об устойчивости нелинейной системы с обратной связью. График проекции траектории замкнутой системы с начальными условиями  на фазовую плоскость 29-10-2021 13-38-53 имеет вид (см. рисунок 4).

29-10-2021 13-43-16

Рис. 4 – Стационарная точка системы с обратной связью

  Заключение

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

Конфликт интересов Не указан. Conflict of Interest None declared.

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

  1. Андриевский Б.Р. Управление хаосом: Методы и приложения. Часть 1. Методы / Б.Р. Андриевский, А.Л. Фрадков // Автоматика и телемеханика. - 2003. - № 5. - С. 3-45.
  2. Анищенко В.С. Лекции по нелинейной динамики / В.С. Анищенко, Т.Е. Вадивасова // Известия вузов. Прикладная нелинейная динамика. - 2010. - Т. 18. - № 3. - С. 186-191.
  3. Малиницкий Г.Г. Структуры. Вычислительный эксперимент: введение в нелинейную динамику / Г.Г. Малиницкий // М.: Эдиториал УРСС, 2002. - 255 с.
  4. Лоскутов А.Ю. Основы теории сложных систем / А.Ю. Лоскутов, А.С. Михайлов // М.: Ижевск: Институт компьютерных исследований, 2007. - 272 с.
  5. Крищенко А.П. Подавление хаотической динамики / А.П. Крищенко, А.В. Кавинов // Дифференциальные уравнения. - 2004. - Т. 40. - № 12. - С. 1629-1635.
  6. Талагаев Ю.В. Стабилизация осциллятора Диффинга-Ван дер Поля оптимальной коррекцией управляющих параметров / Талагаев Ю.В., Тараканов А.Ф. // Письма в ЖТФ. - 2006. - Т. 36. - Вып. 24. - С. 1-9.
  7. Grebogi C. Controlling chaotic dynamical systems / C. Grebogi, Y.C. Lai // Syst. Contr. Lett. -1997. Vol. 31. N 3. P. 307-312.
  8. Ott E. Controlling chaos / E. Ott, C. Grebogi, G. Yorke // Phys. Rev. Lett. - 1990. - Vol. 64. - N 11. - P. - 1196-1199.
  9. Pyragas K. Continuous control of chaos by self-controlling feedback / K. Pyragas // Phys. Lett. A. - 1992. - Vol. 170. - N 6. - P. - 421-428.
  10. Магницкий Н.А. Новые методы хаотической динамики / Н.А. Магницкий, С. В. Сидоров // М.: Едиториал УРСС, 2004. - 320 с.
  11. Yamapi R. Harmonic oscillations, stability and chaos control in a non-linear electromechanical system / R. Yamapi, J.B. Chabi Orou // Journal of Sound and Vibration. - 2003. - Vol. 259. - N 5. - P. - 1253– 1264.
  12. Колесников А.А. Прикладная синергетика: основы системного синтеза / А.А. Колесников // Таганрог: Изд-во ТТИ ЮФУ, 2007. - 384 с.
  13. Кучерова В.Ю. Стабилизация состояний равновесия базовых нелинейных моделей с хаотической динамикой / В.Ю. Кучерова, Л.К. Кильчукова // Вестник АГТУ. Серия: Управление, вычислительная техника и информатика. - 2015. - № 3. - С. - 96-103.
  14. Гробман Д. Гомеоморфизм систем дифференциальных уравнений / Д. Гробман // ДАН СССР. - 1959. - Т. - № 5. - С. 880–881.
  15. Будник С.В. Управление крупномасштабными динамическими системами / С.В. Будник, В.Н. Шашихин //. СПб.: Изд-во Политехпресс, 2020. - 308 с.
  16. Козлов В.Н. Управление энергетическими системами. Ч.1. Теория автоматического управления / В.Н. Козлов, В.Е. Куприянов, В.Н. Шашихин // СПБ.: Изд-во Политехн. ун-та, 2006. - 316 c.
  17. Matsumoto T. Chaotic Attractor from Chua’s Circuit / T.А. Matsumoto // IEEE Transactions on Circuits & Systems. - 1984. - Vol. CAS-31. - N 12. - P. - 1055–1058.

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

  1. Andrievsky B. R.Upravlenie khaosom: Metody i prilozhenija. Chast' 1. Metody [Managing Chaos: Methods and Applications. Part 1. Methods] / B. R. Andrievsky, A. L. Fradkov // Avtomatika i telemekhanika [Automation and telemechanics]. - 2003. - No. 5. - pp. 3-45 [in Russian]
  2. Anishchenko V. S.Lekcii po nelinejjnojj dinamiki [Lectures on nonlinear dynamics] / V. S. Anishchenko , T. E. Vadivasova // Izvestiya vysshikh uchebnykh zavedenii. Fizika’ [Higher Education Institutions Bulletin. Applied nonlinear dynamics]. - 2010. - Vol. 18. - No. 3. - pp. 186-191 [in Russian]
  3. Malinitsky G. G. Struktury. Vychislitel'nyjj ehksperiment: vvedenie v nelinejjnuju dinamiku [Structures. Computational experiment: an introduction to nonlinear dynamics] / G. G. Malinitsky // Moscow: Editorial URSS, 2002. - 255 p. [in Russian]
  4. Loskutov A. Yu. Osnovy teorii slozhnykh sistem [Fundamentals of the theory of complex systems] / Loskutov A. Yu., Mikhailov A. S. // Moscow: Izhevsk: Institute of Computer Research, 2007. - 272 p. [in Russian]
  5. Krishchenko A. P. Podavlenie khaoticheskojj dinamiki [Suppression of chaotic dynamics] / Krishchenko A. P., Kavinov A.V. // Differencial'nye uravnenija [Differential equations]. - 2004. - Vol. 40. - No. 12. - pp. 1629-1635 [in Russian]
  6. Talagaev Yu. V. Stabilizacija oscilljatora Diffinga-Van der Polja optimal'nojj korrekciejj upravljajushhikh parametrov [Stabilization of the Diffing-Van der Pol oscillator by optimal correction of control parameters] / Yu. V. Talagaev, F. Tarakanov // Pisma v ZhTF [Applied Physics Letters]. - 2006. - Vol. 36. - Issue. 24. - p. 1-9 [in Russian]
  7. Grebogi C. Controlling chaotic dynamical systems / C. Grebogi, Y.C. Lai // Syst. Contr. Lett. -1997. Vol. 31. N 3. P. 307-312.
  8. Ott E. Controlling chaos / E. Ott, C. Grebogi, G. Yorke // Phys. Rev. Lett. - 1990. - Vol. 64. - N 11. - P. - 1196-1199.
  9. Pyragas K. Continuous control of chaos by self-controlling feedback / K. Pyragas // Phys. Lett. A. - 1992. - Vol. 170. - N 6. - P. - 421-428.
  10. Magnitsky N. A. Novye metody khaoticheskojj dinamiki [New methods of chaotic dynamics] / N. A. Magnitsky, V. Sidorov // Moscow: Editorial URSS, 2004. - 320 p. [in Russian]
  11. Yamapi R. Harmonic oscillations, stability and chaos control in a non-linear electromechanical system / R. Yamapi, B. Chabi Orou // Journal of Sound and Vibration. - 2003. - Vol. 259. - N 5. - P. - 1253–1264.
  12. Kolesnikov A. A. Prikladnaja sinergetika: osnovy sistemnogo sinteza [Applied synergetics: fundamentals of system synthesis] / Kolesnikov A. A. // Taganrog: Publishing house of Taganrog Technological Institute SFU, 2007. - 384 p. [in Russian]
  13. Kucherova V. Yu. Stabilizacija sostojanijj ravnovesija bazovykh nelinejjnykh modelejj s khaoticheskojj dinamikojj [Stabilization of equilibrium states of basic nonlinear models with chaotic dynamics] / V. Yu. Kucherova, L. K. Kilchukova // Vestnik AGTU. Serija: Upravlenie, vychislitel'naja tekhnika i informatika [Bulletin of ASTU. Series: Management, Computer Engineering and Computer Science]. - 2015. - No. 3. - pp. - 96-103 [in Russian]
  14. Grobman D.Gomeomorfizm sistem differencial'nykh uravnenijj [Homeomorphism of systems of differential equations] / D. Grobman // DAN USSR. - 1959. - Vol. 128. - No. 5. - pp. 880-881 [in Russian]
  15. Budnik S. V. Upravlenie krupnomasshtabnymi dinamicheskimi sistemami [Management of large-scale dynamic systems] / S. V. Budnik, V. N. Shashikhin // Saint Petersburg: Izd-vo Politekhpress, 2020. - 308 p. [in Russian]
  16. Kozlov V. N. Upravlenie ehnergeticheskimi sistemami. Ch.1. Teorija avtomaticheskogo upravlenija [Management of energy systems. Vol.1. Automatic control theory] / Kozlov V. N., Kupriyanov V. E., Shashikhin V. N. // St. Petersburg: Publishing House of the Polytechnic University, 2006. - 316 p. [in Russian]
  17. Matsumoto T. Chaotic Attractor from Chua’s Circuit / T.А. Matsumoto // IEEE Transactions on Circuits & Systems. - 1984. - Vol. CAS-31. - N 12. - P. - 1055–1058.