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

Якубовский  Евгений  Георгиевич

инженер  вычислительного  центра,  Национальный  Минерально-Сырьевой  Университет  «Горный»,  РФ,  г.  Санкт- Петербург

E-mail 

CHAOTIC  SOLUTIONS  ORDINARY  DIFFERENTIAL  EQUATIONS

Evgeniy  Georgievich  Yakubovski

Engineering  Computer  Center,  National  University  of  Mineral-Raw  "Mountain",  Russia  St.  Petersburg

АННОТАЦИЯ

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

ABSTRACT

In  mathematics,  we  introduce  the  notion  of  chaotic  solutions  of  ordinary  differential  equations  .  For  small  deviations  of  initial  data  obtained  large  deviation  solutions.  It  turns  out  that  when  we  have  multiple  coordinates  of  position  of  equilibrium  observed  chaotic  solution.

Ключевые  слова :  хаотическое  решение  обыкновенных  дифференциальных  уравнений;  кратные  координаты  положения  равновесия

Keywords :  chaotic  solution  of  ordinary  differential  equations;  multiple  equilibrium  position  coordinates 

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

  Ставится  задача  нахождение  условий,  когда  система  нелинейных  уравнений  имеет  хаотическое  решение.  Т.  е.  перескок  от  одного  положения  равновесия  к  другому.  Эта  система  уравнений  описывает  скачкообразное,  турбулентное  решение  уравнений  гидродинамики,  когда  локальная  скорость  потока  не  предсказуема  во  времени.  Оказывается,  что  это  условие  наступает,  при  кратных  положениях  равновесия  системы  обыкновенных  нелинейных  уравнений,  к  которым  сводятся  решения  нелинейных  уравнений  в  частных  производных  с  помощью  метода  Галеркина. 

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

    (1) 

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

    (2)

Будем  рассматривать  для    уравнения  координату    независимой,  а  остальные  координаты  как  функции  времени,  являющиеся  решением  дифференциального  уравнения.  Система  уравнений  разбивается  на    независимых  неавтономных  уравнений,  причем    уравнение  зависит  от  одной  переменной    и  от  времени.  В  случае,  когда  решение  равно  координатам  положения  равновесия,  правая  часть  дифференциального  уравнения  не  зависит  от  времени.  При  этом  в  точке  положения  равновесия  с  двукратным  корнем  правая  часть  системы  (2)  представляет  неопределенность  ,  и,  следовательно,  определяет  координаты  не  кратного  положения  равновесия  для  системы  (2).  Числитель  соотношения  неопределенности  равен  квадрату  нуля  из-за  двукратного  положения  равновесия  системы  (1).  Знаменатель  равен  нулю  из-за  вырожденности  преобразования.  При  этом  связь  между  решениями  относительно  преобразования    одинакова  со  связью  между  положениями  равновесия  (связь  вырожденна  в  кратном  положении  равновесия,  но  при  этом  для  остальных  произвольных  значениях    получается  единственное  значение  ).  При  этом  справедливо  следующее  равенство

.

Допустим,  матрица  уравнения  (1)  имеет  двукратное  положение  равновесия.  Тогда  матрица  уравнения  (2),  полученного  из  уравнения  (1),  с  помощью  преобразования    равна

.

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

  Но  преобразование  переменных    не  произвольно.  Преобразование  таково,  что  действительным  или  комплексным  корням    должны  соответствовать  действительные  или  комплексные  корни  .  Строить  систему  (2)  необходимо,  так  как  добиваемся,  чтобы  определитель  линеаризованной  системы  уравнений  (2)  равнялся  нулю.

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

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

  В  предлагаемой  статье  определен  простой  критерий  существования  хаотического  решения.  Наличие  хаотических  решений  связано  с  наличием  кратного  положения  равновесия  системы  уравнений. 

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

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

Т.  е.  определитель  матрицы    равен  нулю  в  этом  представлении.  Решаем  линеаризованное  уравнение  (3)

    (3)

Это  уравнение  имеет  решение  в  окрестности  положения  равновесия  ,  так  как  ,  где  собственные  векторы    и  собственные  числа    определятся  из  системы  уравнений

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

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

,

в  моменты  времени,  когда  по  устойчивым  направлениям  достигнуто  положение  равновесия.  При  этом  решение  увеличится  в  разы

  (4)

Но  в  силу  приближенности  решения,  получится  не  бесконечность,  а  переход  к  другому  положению  равновесия.  Точное  решение  при  численном  счете  получается  мгновенно  при  условии    и  имеет  значение,  зависящее  от  малости  величины  .  Т.  е.  происходит  произвольный  скачок.  Остальные  направления    собственного  вектора,  соответствуют  устойчивому  собственному  числу  положения  равновесия.

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

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

Отметим,  что  подобная  картина  наблюдается  и  у  системы  (1),  так  ее  положения  равновесия  достигаются  одновременно  с  положениями  равновесия  системы  (2)  в  силу  преобразования  координат.  Приближение  будет  по  всем  направлениям  к  положению  равновесия,  и  удаление  по  одному  направлению,  даже  если  нет  собственного  числа  равного  нулю.  Все  равно  будет  существовать  направление  удаления  решения  в  силу  одинаковости  приближения  к  положению  равновесия  этих  систем  уравнений. 

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

.

Причем  координаты  минимума  не  зависят  от  начальных  условий  и  решение  достигает  координат  минимума  в  силу  доказанной  теоремы. 

Аналогичная  теорема  справедлива  для  задачи  Коши  в  комплексной  плоскости.  Т.  е.  при  кратных  корнях  при  комплексных  уравнениях  наблюдается  хаотическое  решение  и  в  случае  решения  в  комплексной  плоскости. 

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

1.Магницкий  Н.А.,  С.В.  Сидоров  Новые  методы  хаотической  динамики.  М.:  Едиториал  УРСС,  2004  г,  —  с.  320.