О существовании периодических решений в системе Лоренца

    image

    Это третий мой топик на Хабре (часть 1 и часть 2), посвященный динамической системе Лоренца. Я продолжаю заниматься исследованием вопроса о существовании периодических решений (циклов) в этой системе. Удалось получить интересный результат при определенном соотношении ее параметров.

    Рассмотрим систему дифференциальных уравнений Лоренца
    image
    где image, r и b — некоторые положительные числа, параметры системы.

    Докажем, что если image, то в системе (1) нет периодических решений (исключая, конечно, положения равновесия).

    Сделаем замену
    image
    где image — некоторая функция от времени.

    Продифференцируем (2), получим
    image

    В левую часть выражения (3) подставим правую часть третьего уравнения системы (1), а в правую часть (3) — правую часть первого уравнения системы (1), учитывая, что image. Получим
    image

    Вместо z в (4) подставим выражение (2), откуда имеем уравнение
    image
    решением которого является функция
    image
    где image — произвольная постоянная.

    Теперь во второе уравнение системы (1) подставим вместо z выражение (2). При этом выразим y из первого уравнения системы (1). Получим
    image
    и
    image

    Подставив (5) и (6) в (7), имеем
    image

    Рассмотрим неавтономный случай, когда image в уравнении (8). Предположим, что в этом случае уравнение (8) имеет периодическое решение с периодом T. Так как производная периодической функции с периодом T есть периодическая функция с периодом T, то левая часть уравнения (8) является периодической функцией с периодом T. Однако правая часть уравнения (8) непериодична, так как image не является периодической функцией. Получили противоречие.

    Таким образом, при image уравнение (8) не имеет периодических решений.

    Рассмотрим теперь случай, когда image. Имеем автономное уравнение второго порядка
    image
    у которого по критерию Бендиксона [1, с. 142-143] нет периодических решений, что и доказывает их отсутствие в системе Лоренца при image.

    Заметим, что в этом случае параметр r может принимать любые значения. Тогда при достаточно больших его значениях в системе Лоренца также будут отсутствовать периодические решения, что кажется весьма неочевидным, поскольку параметр r пропорционален разности температур между нижним и верхним слоем жидкости при свободной конвекции. При увеличении градиента температуры в слое должны возникнуть в жидкости конвективные валы, а здесь жидкость со временем приходит в стационарное состояние (ламинарный режим). Это подтверждается и в численном эксперименте (наблюдались устойчивые фокусы при разных значениях r — рисунок (проекция дуги траектории на плоскость xOy) в начале топика). Скорее всего, это объясняется тем, что система Лоренца достаточно грубо описывает данный процесс, хотя при других соотношениях между image и b (r принимает достаточно большое значение) в системе (1) наблюдается устойчивый предельный цикл [2, с. 291-294].

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

    Рассмотрим другой случай, когда image. В известной мне литературе он исследуется в линейном приближении. Применим второй метод Ляпунова. Составим функцию Ляпунова
    image
    обладающую свойствами:
    image
    image для image
    image при image
    и
    image для image в силу правых частей уравнений системы (1). Тогда по теореме Барбашина-Красовского [3, с. 248-250] любое решение системы (1)
    image при image

    Более общий случай описан в работе [4] (лемма 1.2), где доказывается отсутствие циклов в системе (1) (предельные режимы — положения равновесия) для
    image и image

    Также отметим, что у системы (1) всегда есть частные решения вида
    image
    где image — произвольная постоянная.

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

    Литература


    1. Немыцкий В.В., Степанов В.В. Качественная теория дифференциальных уравнений. — М.: Едиториал УРСС, 2004.
    2. Неймарк Ю.И., Ланда П.С. Стохастические и хаотические колебания. — М: ЛИБРОКОМ, 2009.
    3. Демидович Б.П. Лекции по математической теории устойчивости. — М.: Наука, 1967.
    4. Юдович В.И. Асимптотика предельных циклов системы Лоренца при больших числах Рэлея // Рукопись деп. в ВИНИТИ, №2611-78. 1978.
    5. Robbins K.A. Periodic Solutions and Bifurcation Structure at High R in the Lorenz Model // SIAM Journal on Applied Mathematics, 36(3): 457-472, 1979.
    6. Shimizu T. Analytic Form of the Simplest Limit Cycle in the Lorenz Model // Physica A: Statistical Mechanics and its Applications, 97(2): 383-398, 1979. DOI: 10.1016/0378-4371(79)90113-4.
    7. Покровский Л.А. Решение системы уравнений Лоренца в асимптотическом пределе большого числа Релея. I. Система Лоренца в простейшей квантовой модели лазера и приложение к ней метода усреднения // Теоретическая и математическая физика, 62(2): 272-290, 1985.
    8. Jibin Li, Jianming Zhang. New Treatment on Bifurcations of Periodic Solutions and Homoclinic Orbits at High r in the Lorenz Equations // SIAM Journal on Applied Mathematics, 53(4): 1059-1071, 1993.
    • +19
    • 9,6k
    • 7
    Поделиться публикацией

    Комментарии 7

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

      Вы, кстати теорией катастроф не занимались? Не пробовали в этой системе элементарные катастрофы искать?
        0
        А где можно вкратце почитать?
          0
          Про бифуркационные диаграммы или теорию катастроф?
            +1
            О теории катастроф? Есть хорошая книга Владимира Игоревича Арнольда «Теория катастроф».
              0
              Ох, в какой области он только не работал. Спасибо.
            0
            Только читал.
            Я о периодических решениях уже писал здесь.

        Только полноправные пользователи могут оставлять комментарии. Войдите, пожалуйста.

        Самое читаемое