НАУЧНОЕ ИЗДАНИЕ МГТУ ИМ. Н. Э. БАУМАНА
НАУКА и ОБРАЗОВАНИЕ
Эл № ФС77 - 48211. Государственная регистрация №0421200025. ISSN 1994-0408
электронный научно-технический журнал
Управление хаотической динамикой в системе Спротта <<А>>
# 05, май 2013
Б01:10.7463/0513.0555404
Кавинов А. В.
УДК 517.935.4
Россия, МГТУ им. Н.Э. Баумана kavinov@newmail.ru
Введение
С начала девяностых годов прошлого века внимание исследователей начинает привлекать задача управления хаосом. Суть этой задачи состоит в синтезе такого управления, под воздействием которого динамическая система перестает обладать хаотической динамикой. Обычно замкнутая система имеет при этом один или несколько устойчивык предельных циклов. В качестве примеров способов решения этой задачи можно привести ООУ-метод [1, 2, 3], метод обратной связи с запаздыванием [4, 5], управление с помощью искусственых нейронных сетей [6, 7, 8].
В [9, 10] описан алгоритм решения задачи управления хаосом для аффиннык стационарных систем
х = А(х) + В(х)и, х е П С Ега, и е М1, (1)
эквивалентны« на некотором открытом подмножестве пространства состояний П регулярной системе канонического вида [11]
у(п) = f(у) + д(у)и, у = (у, у,..., у[п-1))т. (2)
В предположении, что система (1), замкнутая управлением и = 0, обладает хаотической динамикой, найдено управление, при котором эта система не имеет хаотической динамики.
Предложенный метод заключался в стабилизации некоторого специально найденного движения системы (2). В результате замкнутая полученным управлением система имела асимптотически устойчивый предельный цикл. Если при этом решения исходной системы обадали следующим свойством: для любой заранее заданной положительной величины е на любом решении с хаотическим характером движения существуют такие две точки х(^1) и х(Ь2), где ¿2 > Ь, что ||х(^) — х(£2)|| < е (при этом ¿2 — t1 ^ 8 > 0), то полученное управление могло быть сделано сколь угодно малым по абсолютной величине.
Описанный метод неприменим к системам, эквивалентным системам квазиканонического вида [12, 13]
У(к) = /(У) + #(УК »>1 = МУ^ Чп-к = ^п-Й^^ (3)
У = У, ... ,У(п-1))Т, П = (П1, ..., Пп-к)Т.
Это вызвано тем, что для системы (3) невозможно, вообще говоря, стабилизировать произвольно заданное движение.
В данной статье рассматривается задача подавления хаотической динамики в системе, не эквивалентной регулярной системе канонического вида ни на каком открытом подмножестве пространства состояний.
Постановка задачи
В [15] описано 19 простых систем третьего порядка с квадратичной правой частью, обладающих хаотической динамикой. Первая из этих систем, обозначенная автором буквой «Л», имеет вид
ж = у,
У = —ж + уг,
¿ = 1 — у2.
Если ввести управление во второе уравнение этой системы, получим систему квазиканонического вида
(ж = У,
У = —ж + уг + и, (4)
¿=1 — У2,
не эквивалентную никакой системе канонического вида ни на каком подмножестве К3 [14].
Построим такое управление, при котором вместо хаотической динамики замкнутая им система (4) будет иметь устойчивые периодические решения.
Создание периодических решений
Пусть функция € С2 является периодической с периодом Т. Рассмотрим управление
и(ж, У,г, ¿) = ж — Уг + г/(£). (5)
Тогда решение задачи Коши (4)-(5) с начальными условиями
ж (¿о) = жо = ^(¿о), У (¿о) = Уо = У(£о), ¿(¿о) = ¿0 (6)
имеет вид
г
ж(*) = У(*) = т ¿(¿) = £(¿) = ¿о + I (1 — п2(т)) ¿т. (7)
го
Заметим, что это решение принадлежит цилиндрической поверхности с образующей, параллельной оси г. Параметрические уравнения
x = n(t) У = n(t) (8)
задают направляющую этой цилиндрической поверхности. Дополнительно потребуем, чтобы для периодической функции n(t) для любого t выполнялось условие
t+T
J (1 - П2(т)) dT = 0. (9)
t
Тогда для любого t
t t+T t
£(t + T) = zo + / (1 - П2(t)) dT + J (1 - П2(т)) dT = zo + / (1 - П2(t)) dT = £(t).
to t to
Это означает, что компоненты решения (7) являются функциями периода T, T и kT соответственно, где k — натуральное число. Будем искать функцию n(t), удовлетворяющую условию (9), в виде
2n(t - to) , , . 2n(t - to) (10)
n(t) = a cos-t--+ b sin-t-, (1°)
где a и b — некоторые константы. Тогда
2п . 2n(t - to) ,2n 2n(t - to)
n(t) = -a— sin-—--+ b— cos
T T T T
t+T
í (1 - n2(t)) dT = T -
2n2(a2 + b2)
T
t
Условие (9)и
^(¿о) = хо, п(^о) = уо (11)
будут выполняться в том случае, если
2П
а = хо, Ь — = уо, Т2 = 2п2 (а2 + Ь2). (12)
Система (12) при |уо| < \[2 имеет решение
а = хо, Ь = Т = (13)
2п V 2 — уо2
Таким образом, для любой начальной точки (хо,уо,го), |уо| < л/2, мы построили управление (5),(10),(13), при котором система совершает периодическое движение по замкнутой кривой, лежащей на эллиптическом цилиндре с образующей, параллельной оси г. Для начальный точек (хо, уо, го) и (хо, уо, г1), отличающихся только значениями третьих координат,
соответствующие им фазовые кривые будут принадлежать одному и тому же цилиндру, причем вторая кривая может быть получена из первой параллельным сдвигом вдоль оси z на величину z\ — z0.
В рассматриваемой системе цилиндр состоит из периодических траекторий, три из которых изображены на рис. 1. Изображенные траектории получены интегрированием системы (4), замкнутой полученным управлением (5), для начальных точек (Тл/2/п, 0, ü), (TV2/n, 0, 4), (TV2/n, 0, 8), при T = 100.
Рис. 1. Фазовые кривые
Стабилизация периодических решений
Рассмотрим теперь задачу Коши с произвольными начальными условиями
x(to ) = x*, y(to) = y*, z(to) = z* (14)
для системы (4), замкнутой управлением
u(x,y,z,t) = x - yz + n(t) - ki(y - rj(t)) - ko(x - nit)), (15)
где функция nit) построена для начальной точки (x0, y0, z0) (6); k\ и k0 — некоторые константы. В точках траектории x(t) = n(t), y(t) = n(t), z(t) = £(t) это управление совпадает с управлением (5). Подставив (15) в исходную систему (4), получим замкнутую систему
x = У,
y = nit) - ki(y - nit)) - ko(x - nit)), (16)
z = 1 - y2.
Первые два уравнения (16) запишем в переменных Дх = ж — п(£), Ду = у — П(£):
аДх
Ду,
^ (17)
аДу
—— = Ду + коДх.
При положительных к1 и к0 эта линейная система с постоянными коэффициентами будет иметь асимптотически устойчивое в целом нулевое решение.
Поэтому Дх и Ду стремятся к нулю при £ ^ Это означает, что расхождение между направляющей (8) цилиндрической поверхности и проекцией на плоскость Ожу решения задачи Коши (4),(15),(14) стремится к нулю при £ ^
Пусть далее к1 и к0 таковы, что характеристическое уравнение системы (17) имеет вещественные отрицательные корни А1 и Л2, А1 > Л2. Тогда общее решение (17) будет иметь вид
Дх = С1вЛ1(4-4о) + С2вЛ2(4-4о), Ду = С1А1вЛ1(4-4о) + С^е^Ч
Значения постоянных С1 и С2, соответствующие решению задачи Коши с начальными условиями
Дх(£о) = ж* — жо, Ду(£о) = у* — уо, (18)
могут быть найдены по формулам
с = —А2Дх(£о) + Ду(£о) с = А1Дх(£о) — Ду(£о) 1 А1 — А2 2 А1 — А2
Теперь рассмотрим вопрос о поведении Дг = г(£) — £(£) при £ ^ где г(£) — третья компонента решения задачи Коши (4), (15), (14). Поскольку
г
г(£) = г* + У(1 — у2(т)) От,
го
то
г г
Дг = г* + у (1 — у2(т)) ат — го ^(1 — п2(т)) ^т =
го го
г г
= — го ^ У (у2(т) — ?)2(т)) ат = — го — J (у(т) — п(т))(у(т) + п(т)) ат =
го го
г
= г* — го — У Ду(т) (Ду(т) + 2т)(т)) ¿т;
го
при этом величина
Ду(т) = С1А1вЛ1(т-го) + С2А2вЛ2(т-го)
бесконечно мала при t ^ а выражение
а(т) = Ау(т) + -г](т) =
= Ci\ieXl(T-'О + С2Х2вХ2(т-'О — 4ТП sin + ^ cos ^^^ =
= CiAie^i(T-'о) + C2A2e^2(T-tO — ^ sin + -у, cos
ограничено при t G [to, как сумма ограниченых функций. Обозначим константу,
ограничивающую а(т) по абсолютной величине при t G [t0, через M.
+t»
Покажем, что интеграл f Лу(т)а(т) dт сходится. Действительно,
to
/ Ау(т)а(т) Лт = у С1\1вХ1(т) а(т) Лт + у С2\2еХ2(т-0) а(т) Лт.
Ьо 10 Ьо
Для первого интеграла в правой части равенства имеем
CiAieAl(T-о)а(т) dт
to
|CiAi|
eXl(T-to)а(т) dт
to
Так как
и
eXl(T-'о)а(т)| < M |eXl(T-to)| = MeXl(T-to)
MeXKT-o) dт = - M M
to
Ai |Ai |
то
интеграл f CiAieXl(T to)а(т) dт сходится, причем
to
CiAieXl(T-^а(т) dт
to
< |CiAi| ■ |M| = |Ci|M.
Аналогично
C2A2eX2(T-^)а(т) dт
to
< CIM.
Таким образом, существует конечный предел
lim Az(t) = lim (z(t) — £(t)) = C, t^+t t^+t
где величина
C = J Ау(т)а(т) dт
to
зависит лишь от начальных условий (14). Следовательно, ^-предельное множество решения задачи Коши (4), (15), (14) будет замкнутой кривой
х(£) = п(£), у(£) = П(£), г (£) = £ (£) + С,
лежащей на цилиндрической поверхности с направляющей (8).
Результаты моделирования процесса управления хаосом при помощи построенного управления приведены на рис. 2. Моделирование производилось для трех начальных точек 51 (1, 1, 1), 52(—1, —1, —20) и 53(10, 3, 0) с одинаковым значением Т = 100; были выбраны значения к1 = 0,3 и ко = 0,02, жо = Тл/2/п, уо = 0. На рисунке синим цветом показана тректория при движении из 51, красным — из 52, зеленым — из 53.
.2 -30
Рис. 2. Моделирование процесса управления хаосом Заключение
В статье описан способ, позволяющий решить задачу управления хаосом в одной из систем Спротта с управлением. Рассматриваемая система примечательна тем, что она не эквивалентна регулярной системе канонического вида ни на каком подмножестве пространства состояний. Очевидно, что описанный способ допускает обобщение на более широкий класс систем. Так, в последнем уравнении системы правая часть может иметь вид произвольного многочлена от ж и у. При этом в процедуре управления изменится лишь вид последнего условия в (12).
Работа выполнена при поддержке грантов РФФИ 11-01-00733 и Программы Президента РФ по государственной поддержке ведущих научных школ (грант №НШ-3659.2012.1).
Список литературы
1. Ott E., Grebogi C., Yorke J.A. Controlling chaos // Phys. Rev. Lett. 1990. Vol. 64, no. 11. P. 1196-1199.
2. Grebogi C., Lai Y.C. Controlling chaos in high dimensions // IEEE Trans. Circ. Syst. I. 1997. Vol. 44. P. 971-975.
3. Grebogi C., Lai Y.C. Controlling chaotic dynamical systems // Syst. Contr. Lett. 1997. Vol. 31, no. 3. P. 307-312.
4. Nitsche G., Dressler U. Controlling chaotic dynamical systems using time delay coordinates // Physica D: Nonlinear Phenomena. 1992. Vol. 58, no. 1-4. P. 153-164. DOI: 10.1016/0167-2789(92)90106-W.
5. Pyragas K. Continuous control of chaos by self-controlling feedback // Phys. Lett. A. 1992. Vol. 170, no. 6. P. 421-428. DOI: 10.1016/0375-9601(92)90745-8.
6. Alsing P.M., Gavrielides A., Kovanis V. Using Neural networks for controlling chaos // Phys. Rev. E. 1994. Vol. 49. P. 1225-1231.
7. Weeks E.R., Burgess J.M. Evolving artifical neural networks to control chaotic systems // Phys. Rev. E. 1997. Vol. 56, no. 2. P. 1531-1540.
8. Lin C.T., Jou C.P. Controlling chaos by GA-based reinforcement learning neural network // IEEE Trans. Neural Netw. 1999. Vol. 10. P. 846-859.
9. Крищенко А.П., Кавинов А.В. Подавление хаотической динамики // Дифференциальные уравнения. 2004. Т. 40, № 12. С. 1629-1635.
10. Кавинов А.В. Стабилизация программных движений и подавление хаотической динамики малыми управлениями // Нелинейная динамика и управление. Вып. 5: Сб. статей /Под ред. С.В. Емельянова, С.К. Коровина. М.: ФИЗМАТЛИТ, 2006. С. 175-184.
11. Краснощёченко В.И., Крищенко А.П. Нелинейные системы: геометрические методы анализа и синтеза. М.: Издательство МГТУ им. Н.Э. Баумана, 2005. 520 с.
12. Крищенко А.П., Клинковский М.Г. Преобразование аффинных систем с управлением и задача стабилизации // Дифференциальные уравнения. 1992. Т. 28, № 1. С. 1945-1952.
13. Фетисов Д.А. Исследование управляемости регулярных систем квазиканонического вида // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2006. № 3. С. 12-30.
14. Ткачев С.Б., Шевляков А.А. Преобразование аффинных систем со скалярным управлением к квазиканоническому виду // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2013. №1. С. 3-16.
15. Sprott J.C. Some simple chaotic flows // Phys. Rev. E. 1994. Vol. 50, no. 2. P. 647-650.
SCIENTIFIC PERIODICAL OF THE BAUMAN MSTU
SCIENCE and EDUCATION
EL № FS77 - 48211. №0421200025. ISSN 1994-0408
electronic scientific and technical journal
Control of chaotic dynamics in the Sprott A system
# 05, May 2013
DOI: 10.7463/0513.0555404
Kavinov A. V.
Bauman Moscow State Technical University 105005, Moscow, Russian Federation kavinov@newmail.ru
Since the early nineties, the problem of controlling chaos attracts researches' attention. The essence of the problem lies in synthesis of control which eliminates chaotic dynamics of a dynamic system. A normally closed system has one or several stable limit cycles. This article deals with the problem of controlling chaos in one of the Sprott systems and also describes the method of elimination of chaotic dynamics by stabilizing a cylindrical invariant sub-manifold in the phase space of the closed system which consists of system's stable periodic solutions. The system under consideration belongs to the class of affine systems which are not equivalent to a regular system of the canonical form on any subset of the state space. The described method is applicable not only to the specific system in question but may be extended.
References
1. Ott E., Grebogi C., Yorke J.A. Controlling chaos. Phys. Rev. Lett., 1990, vol. 64, no. 11, pp.1196-1199.
2. Grebogi C., Lai Y.C. Controlling chaos in high dimensions. IEEE Trans. Circ. Syst. I., 1997, vol. 44, no. 10, pp. 971-975. DOI: 10.1109/81.633886.
3. Grebogi C., Lai Y.C. Controlling chaotic dynamical systems. Syst. Contr. Lett., 1997, vol. 31, no. 3, pp. 307-312.
4. Nitsche G., Dressler U. Controlling chaotic dynamical systems using time delay coordinates. PhysicaD: Nonlinear Phenomena, 1992, vol. 58, pp. 153-164. DOI: 10.1016/0167-2789(92)90106-W.
5. Pyragas K. Continuous control of chaos by self-controlling feedback. Phys. Lett. A, 1992, vol. 170, no. 6, pp. 421-428. DOI: 10.1016/0375-9601(92)90745-8.
6. Alsing P.M., Gavrielides A., Kovanis V. Using Neural networks for controlling chaos. Phys. Rev. E, 1994, vol. 49, pp. 1225-1231.
7. Weeks E.R., Burgess J.M. Evolving artifical neural networks to control chaotic systems. Phys. Rev. E, 1997, vol. 56, pp. 1531-1540.
8. Lin C.T., Jou C.P. Controlling chaos by GA-based reinforcement learning neural network. IEEE Trans. NeuralNetw., 1999, vol. 10, pp. 846-859.
9. Krishchenko A.P., Kavinov A.V. Podavlenie khaoticheskoi dinamiki [Suppression of chaotic dynamics]. Differentsial'nye uravneniia, 2004, vol. 40, no. 12, pp. 1629-1635. (Trans. version: Differential Equations, 2004, vol. 40, no. 12, pp. 1709-1715. DOI: 10.1007/s10625-005-0102-9.)
10. Kavinov A.V. Stabilizatsiia programmnykh dvizhenii i podavlenie khaoticheskoi dinamiki malymi upravleniiami [Stabilization programmed motion and suppression of chaotic dynamics using small controls]. Nelineinaia dinamika i upravlenie. Vyp. 5: Sbornik statei [Nonlinear dynamics and control. Iss. 5: collection of articles]. Moscow, PhysMatLit, 2006, pp. 175-184.
11. Krasnoshchechenko V.I., Krishchenko A.P. Nelineinye sistemy: geometricheskie metody analiza i sinteza [Nonlinear systems: geometrical methods of analysis and synthesis]. Moscow, Bauman MSTU Publ., 2005. 520 p.
12. Krishchenko A.P., Klinkovskii M.G. Preobrazovanie affinnykh sistem s upravleniem i zadacha stabilizatsii [The transformation of affine systems with control and stabilization problem]. Differentsial'nye uravneniia, 1992, vol. 28, no. 1, pp. 1945-1952.
13. Fetisov D.A. Issledovanie upravliaemosti reguliarnykh sistem kvazikanonicheskogo vida [Study of controllability of regular systems of quasi-canonical type]. VestnikMGTU im. N.E. Baumana. Ser. Estestvennye nauki [Bulletin of the Bauman MSTU. Ser. Natural science], 2006, no. 3, pp. 12-30.
14. Tkachev S.B., Shevliakov A.A. Preobrazovanie affinnykh sistem so skaliarnym upravleniem k kvazikanonicheskomu vidu [The transformation of affine systems with scalar control to the quasi-canonical form]. VestnikMGTUim. N.E. Baumana. Ser. Estestvennye nauki [Bulletin of the Bauman MSTU. Ser. Natural science], 2013, no. 1, pp. 3-16.
15. Sprott J.C. Some simple chaotic flows. Physical review E, 1994, vol. 50, no. 2, pp. 647-650.