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

УДК 517.925

МЕТОДЫ ИССЛЕДОВАНИЯ ЗАДАЧ СИНХРОНИЗАЦИИ ПЕРИОДИЧЕСКИХ КОЛЕБАНИЙ В НЕЛИНЕЙНЫХ СИСТЕМАХ УПРАВЛЕНИЯ © М.Г. Юмагулов

Ключевые слова: вынужденные колебания; динамические системы; бифуркации; синхронизация.

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

Одно из фундаментальных свойств нелинейных динамических систем - это явление синхронизации, которое заключается в установлении определенных соотношений между параметрами системы и возникновении вследствие этого в системе незатухающих колебаний определенных периодов [1, 2]. Многие задачи приводят к необходимости исследования вопроса о синхронизации автоколебательной системы, имеющей собственную частоту щ, внешним периодическим сигналом частоты V. Здесь при определенных отношениях частот в = у/уо и амплитуды а внешнего сигнала возможно возникновение колебаний разных периодов. В частности, в системе могут возникнуть колебания, период которых является числом кратным периоду Т внешнего сигнала. Такое явление называют синхронизацией на субгармониках, а возникающие периодические решения - субгармоническими колебаниями нелинейной системы. Задачи о синхронизации на субгармониках изучались многими авторами [2, 3]. Здесь особо важны исследования, направленные на поиск признаков такой синхронизации, приближенное построение возникающих решений, анализ их устойчивости.

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

Рассматривается зависящая от скалярного или векторного параметра у динамическая система, описываемая дифференциальным уравнением

X = Е

в котором Е(х,у) и д(Ь,у) - это гладкие по совокупности переменных функции, при этом функция д(Ь, у) является Т -периодической по £: д(Ь + Т, у) = д(Ь, у) . Предполагается, что автономная система

X = Е

имеет точку равновесия х*(у), которую без ограничения общности можно считать нулевой,

т. е. Е(0,у) = 0 . Предполагается также, что при некотором у = уо внешний периодический

сигнал д

2751

при у = уо имеет пару простых собственных значений ±Уо1, где Уо > 0, и не имеет других чисто мнимых собственных значений.

Систему (2) можно представить в виде х! = А(у)х + а(х,у), где нелинейность а(х,у) является однородной по х и начинается с квадратичных слагаемых. Уравнение (1) можно рассматривать как динамическую систему, имеющую собственную частоту Уо и на которую воздействует внешний периодический сигнал д(Ь, у) частоты у = 2п/Т.

Пусть д - натуральное число. Будем говорить, что значение уо параметра у является точкой синхронизации на субгармониках периода дТ системы (1) , если существуют число ео > 0 и определенная при е € [0, ео) непрерывная функция у(е) такая, что у(0) = уо и при этом для каждого е € (0,ео) уравнение (1) при у = у(е) имеет ненулевое периодическое решение х = х(Ь,е) минимального периода дТ такое, что функция х(Ь,е) непрерывно зависит от е, при этом ||х(£,е)||^ 0 при е ^ +0 равномерно по £ .

Отметим, что синхронизация на субгармониках периода дТ является типичным явлением в нелинейной системе (1), если отношение частот Уо и у является рациональным, а именно, если выполнено равенство

Уо = Р

у д

при некоторых натуральных р и д.

Пусть условие (3) выполнено. Для исследования задачи о синхронизации на субгармониках периода дТ предлагается следующая схема. На первом этапе от дифференциального уравнения (1) осуществляется переход к операторному уравнению

У = В

где

яТ

ад = е^>, «у,у) = / е*Т-'>А<">а0ф),у) Л,

о

яТ

и

(у) = У ея(Т 3'>л(^>д(в,у) йв ;

здесь х(£) - это решение задачи Коши для дифференциального уравнения (1) при начальном условии х(0) = у. Задача о дТ -периодических решениях уравнения (1) равносильна задаче о решениях уравнения (4) в следующем смысле: если функция х(£) является дТ -периодическим решением уравнения (1), то вектор у = х(0) будет решением уравнения (4) и, наоборот, если вектор у является решением уравнения (4) то решение х(£) задачи Коши для дифференциального уравнения (1) при начальном условии х(0) = у будет дТ -периодическим решением уравнения (1).

На втором этапе исследование операторного уравнения можно проводить на основе разработанных в [4, 5] операторных методов исследования задач о многопараметрических бифуркациях динамических систем.

Так как матрица А(уо) имеет пару простых собственных значений ±Уо* и не имеет других чисто мнимых собственных значений, то в силу условия (3) матрица В (уо) имеет полупростое собственное значение 1 кратности 2. Тогда найдутся две пары линейно независимых векторов е, д и е*, д* таких, что векторы е и д являются собственными для матрицы В (уо), отвечающими собственному значению 1, а векторы е* и д* являются собственными для транспонированной матрицы В*(уо), отвечающими собственному значению 1. При этом они могут быть выбраны исходя из соотношений:

(е е*) = (д, д*) = 1 (е д*) = (д, е*) = 0.

2752

В рассматриваемой задаче естественным будет предположение, что параметр у является двумерным, т. е. у = (а,в), где а и в - скалярные параметры; положим также уо = = (ао,во)- Тогда уравнение (4) примет вид

у = В (а, в)у + Ь(х,а,в)+ и(а,в), у € Ем.

Определим матрицу

Б = Г «

= [

Здесь и'а и и'в - векторы, полученные дифференцированием вектора и(а, в) по а и в соответственно.

Теорема 1. Пусть ёе! Б = 0 . Тогда значение уо = (ао,во) параметра у = (а,в) является точкой синхронизации на субгармониках периода дТ системы (1).

Рассмотрим теперь вопрос о построении функции у(е), участвующей в определении понятия точки синхронизации на субгармониках. Так как в данном случае имеем у = (а, в) и уо = (ао,во), то речь идет о двух непрерывных скалярных функциях а(е) и в(е), таких, что а(0) = ао и в(0) = во-Теорема 2. Пусть

(Ь(е,ао,во),е*) = 0, (Ь(е, ао,во),д*) = 0 .

Тогда в условиях теоремы 1 уравнение (1) имеет субгармонические решения при у = = (а(е),в(е)), где функции а(е) и в(е) представимы в виде

а(е) = ао + е2а2 + аз(е), в(е) = во + е2в2 + вз(е);

здесь числа а2 и в2 являются ненулевыми, а функции а3(е) и вз(е) - это непрерывные функции, удовлетворяющие соотношениям: а3(е) = о(е2) , вз(е) = о(е2) при е 0 .

При доказательстве теорем 1 и 2 были получены конечные формулы для чисел а2 и в2-Из теоремы 2, в частности, следует, что явление синхронизации на субгармониках периода дТ системы (1) имеет, как правило, “направленный” характер: например, если а2 > 0 и в2 > 0, то субгармонические колебания возникают при а> ао и в> во-

ЛИТЕРАТУРА

1. Гукенхеймер Дж., Холмс Ф. Нелинейные колебания, динамические системы и бифуркации векторных полей. Москва; Ижевск: Институт компьютерных исследований, 2002.

2. Анищенко В.С., Вадивасова Т.Е. Лекции по нелинейной динамике. Москва; Ижевск: НИЦ «Регулярная и хаотическая динамика», 2011.

3. Розо М. Нелинейные колебания и теория устойчивости. М.: Наука, 1971.

4. Юмагулов М.Г. Операторный метод исследования правильной бифуркации в многопараметрических системах. // Доклады АН. 2009. Т. 424. № 2. C. 177-180.

5. Вышинский А.А., Ибрагимова Л.С., Муртазина С.А., Юмагулов М.Г. Операторный метод приближенного исследования правильной бифуркации о многопараметрических динамических системах // Уфимский математический журнал. 2010. Т. 2. № 4. С. 3-26.

БЛАГОДАРНОСТИ: Работа выполнена при поддержке гранта Минобрнауки РФ, соглашение 14.B37.21.0358 «Развитие новых направлений спектральной теории и теории функций, их приложения в задачах математической физики и нелинейной динамики».

Yumagulov M.G. METHODS OF RESEARCH PROBLEMS SYNCHRONIZATION OF PERIODIC OSCILLATIONS IN NONLINEAR CONTROL SYSTEMS

In this paper we study the problem of synchronization of forced oscillations in nonlinear control systems. A new scheme of study of synchronization problems in the sub-harmonics, based on the operator

2753

method for the study of problems of multi-parameter bifurcations is proposed. The scheme allows you to define the new conditions for signs of synchronization and asymptotic formula for the resulting solutions. Key words: forced oscillations; dynamical system; bifurcation; synchronization.

УДК 517.92

ПРИЗНАКИ СУБГАРМОНИЧЕСКОЙ БИФУРКАЦИИ ДЛЯ ФУНКЦИОНАЛЬНО-ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ ЗАПАЗДЫВАЮЩЕГО ТИПА © М.Г. Юмагулов, Д.А. Якшибаева

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

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

Постановка задачи.

Рассмотрим систему функционально-дифференциальных уравнений запаздывающего типа, зависящую от векторного параметра А с Т -периодической по Ь правой частью:

Г

х'

о

где 0 <г ^ Т, Q(А, Ь, т) — квадратная N х N матрица, элементы которой при каждом А являются функциями ограниченой вариации по Ь € [0, г] и при каждом Ь € € [0,г] непрерывно дифференцируемы по А ; х1 = (х(Ь — $1),---,х(Ь — я3)), Я] € [0,Т),] = = 1,...,з; нелинейность а(Ь,А,х^ равномерно по А и Ь удовлетворяет соотношению \\а(Ь, А,х^\\ = О (||х||2) при ||х|| ^ 0 . Интегралы (1) понимаются в смысле Лебега-Стилтьеса, || • \| — норма векторов в евклидовом пространстве Ям .

Уравнение (1) при всех значениях А имеет нулевое решение. В статье исследуется задача о бифуркации периодических решений уравнения (1) в окрестности точки х = =0.

Наряду с (1) будем рассматривать линейную систему

Г

x'

о

Критическими для уравнения (1) будут те значения Ао параметра А, при которых один или несколько мультипликаторов системы (2) по модулю равны 1. Изменение

2754