Бифуркационная диаграмма: различия между версиями

Материал из sawiki
Перейти к навигации Перейти к поиску
(Вторая версия страницы)
Строка 9: Строка 9:
 
\begin{gather*}
 
\begin{gather*}
 
\begin{cases}
 
\begin{cases}
   \dot{u}_1=f(u), ~ u \in \mathbb{R} \\
+
   \dot{u}=f(u), ~ u \in \mathbb{R} \\
   \dot{v}_2=g(v), ~ v \in \mathbb{R}.
+
   \dot{v}=g(v), ~ v \in \mathbb{R}.
 
  \end{cases}
 
  \end{cases}
 
\end{gather*}
 
\end{gather*}
  
Зафиксируем некоторое значение вектора параметров $ a = a_0 $ и рассмотрим в
+
Зафиксируем некоторое значение вектора параметров $$ a = a_0 $$ и рассмотрим в
пространстве параметров максимальное связное множество, содержащее $ a_0, $ такое,
+
пространстве параметров максимальное связное множество, содержащее $$ a_0, $$ такое,
 
что во всех его точках вышеприведённая система топологически эквивалентна этой же системе
 
что во всех его точках вышеприведённая система топологически эквивалентна этой же системе
при $ a = a_0. $ Рассматривая такие множества в пространстве параметров, получим так
+
при $$ a = a_0. $$ Рассматривая такие множества в пространстве параметров, получим так
 
называемый параметрический портрет вышеприведённой системы. Параметрический портрет
 
называемый параметрический портрет вышеприведённой системы. Параметрический портрет
 
вместе с характерными для каждого множества параметров фазовыми портретами
 
вместе с характерными для каждого множества параметров фазовыми портретами
 
составляют бифуркационную диаграмму.
 
составляют бифуркационную диаграмму.
 +
 
При качественном анализе динамической системы желательно получить ее би-
 
При качественном анализе динамической системы желательно получить ее би-
 
фуркационную диаграмму, так как в ней в сжатом виде содержатся все возможные
 
фуркационную диаграмму, так как в ней в сжатом виде содержатся все возможные
 
модели поведения данной системы.
 
модели поведения данной системы.
 +
 +
[[Файл:BifDiagNew.png|300px|мини|справа|Бифуркационная диаграмма для примера 1 при значении $$ v_0 = 0.4 $$.]]
 +
 +
== Алгоритм построения бифуркационной диаграммы ==
 +
Определяем отрезок $$ [ \alpha_0, \alpha_{max} ] $$ на котором будем производить исследование.
 +
 +
Запускаем цикл по j от 1 до N, перебирая по равномерной сетке все значения $$ \alpha. $$
 +
 +
Для каждого значения $$ \alpha_j $$ запускаем цикл по i от 1 до N, в котором находим $$ v_i = g(v_{i-1}). $$ Таким образом, получаем $$ v_N. $$
 +
 +
Для каждого такого $$ \alpha_j, $$ при уже найденном $$ v_N, $$ запускаем новый цикл по i от 1 до M, в котором производим ещё M итераций нахождения $$ v_i = g(v_{i-1}). $$ Каждое полученное значение $$ v_i $$ заносим на график как точку с координатами $$ ( \alpha_j, v_i) . $$
 +
 +
Бифуркационная диаграмма получена.
 +
[[Файл:BifDiag3New.png|300px|мини|слева|Бифуркационная диаграмма в фазово-параметрическом пространстве динамической системы примера 2.
 +
Жирная линия обозначает многообразие положений равновесия, причем сплошная линия отвечает устойчивому положению равновесия, а
 +
пунктирная — неустойчивому. Жирные точки обозначают положения равновесия. Приведены
 +
три фазовых портрета, как сечения бифуркационной диаграммы, для значений параметра $$ a_1, a_c, a_2. $$]]
  
 
== Примеры бифуркационных диаграмм==
 
== Примеры бифуркационных диаграмм==
'''Пример.'''
+
'''Пример 1.'''
 +
Пусть задана динамическая система:
 +
 
 +
\begin{gather*}
 +
\dot{v}_{t+1}= \alpha {v_t}^{\frac{3}{2}} (1 - v_t).
 +
\end{gather*}
 +
Тогда бифуркационные диаграммы приведены в правой части страницы.
 +
N - число итераций необходимое для того, чтобы траектория сошлась к некоторому постоянному состоянию,
 +
 
 +
M - число итераций необходимое для того, чтобы найти возможные положения равновесия в данном состоянии,
 +
 
 +
$$ v_0 $$ - точка, в которой происходит поиск.
 +
 
 +
[[Файл:BifDiag2New.png|300px|мини|справа|Бифуркационная диаграмма для примера 1 при значении $$ v_0 = 0.3 $$.]]
 +
 
 +
'''Пример 2.'''
 +
Пусть задана динамическая система:
 +
 
 +
\begin{gather*}
 +
\dot{u}= a + u^2 = f(u; a).
 +
\end{gather*}
  
 +
'''Пример 3.'''
 
Пусть задана динамическая система:
 
Пусть задана динамическая система:
  
 
\begin{gather*}
 
\begin{gather*}
\dot{v_{t+1}}_2= \alpha {v_t}^{\frac{3}{2}} (1 - v_t).
+
\dot{u}= au - u^3, ~ a \in \mathbb{R}.
 
\end{gather*}
 
\end{gather*}
Тогда бифуркационные диаграммы выглядят так:
 
  
[[Файл:BifDiag.png|300px|мини|справа]]
+
[[Файл:BifDiag4New.png|300px|мини|слева|Бифуркационная диаграмма в фазово-параметрическом пространстве
[[Файл:BifDiag2.png|300px|мини|справа]]
+
динамической системы примера 3. Обозначения как для примера 2.]]
  
 
== Список литературы ==
 
== Список литературы ==
 
# Абрамова В.В. "Лекции по динамическим системам и биоматематике", 2023.
 
# Абрамова В.В. "Лекции по динамическим системам и биоматематике", 2023.

Версия 17:34, 23 сентября 2023

Определение

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

Смысл

\begin{gather*} \begin{cases} \dot{u}=f(u), ~ u \in \mathbb{R} \\ \dot{v}=g(v), ~ v \in \mathbb{R}. \end{cases} \end{gather*}

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

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

Бифуркационная диаграмма для примера 1 при значении $$ v_0 = 0.4 $$.

Алгоритм построения бифуркационной диаграммы

Определяем отрезок $$ [ \alpha_0, \alpha_{max} ] $$ на котором будем производить исследование.

Запускаем цикл по j от 1 до N, перебирая по равномерной сетке все значения $$ \alpha. $$

Для каждого значения $$ \alpha_j $$ запускаем цикл по i от 1 до N, в котором находим $$ v_i = g(v_{i-1}). $$ Таким образом, получаем $$ v_N. $$

Для каждого такого $$ \alpha_j, $$ при уже найденном $$ v_N, $$ запускаем новый цикл по i от 1 до M, в котором производим ещё M итераций нахождения $$ v_i = g(v_{i-1}). $$ Каждое полученное значение $$ v_i $$ заносим на график как точку с координатами $$ ( \alpha_j, v_i) . $$

Бифуркационная диаграмма получена.

Бифуркационная диаграмма в фазово-параметрическом пространстве динамической системы примера 2. Жирная линия обозначает многообразие положений равновесия, причем сплошная линия отвечает устойчивому положению равновесия, а пунктирная — неустойчивому. Жирные точки обозначают положения равновесия. Приведены три фазовых портрета, как сечения бифуркационной диаграммы, для значений параметра $$ a_1, a_c, a_2. $$

Примеры бифуркационных диаграмм

Пример 1. Пусть задана динамическая система:

\begin{gather*} \dot{v}_{t+1}= \alpha {v_t}^{\frac{3}{2}} (1 - v_t). \end{gather*} Тогда бифуркационные диаграммы приведены в правой части страницы. N - число итераций необходимое для того, чтобы траектория сошлась к некоторому постоянному состоянию,

M - число итераций необходимое для того, чтобы найти возможные положения равновесия в данном состоянии,

$$ v_0 $$ - точка, в которой происходит поиск.

Бифуркационная диаграмма для примера 1 при значении $$ v_0 = 0.3 $$.

Пример 2. Пусть задана динамическая система:

\begin{gather*} \dot{u}= a + u^2 = f(u; a). \end{gather*}

Пример 3. Пусть задана динамическая система:

\begin{gather*} \dot{u}= au - u^3, ~ a \in \mathbb{R}. \end{gather*}

Бифуркационная диаграмма в фазово-параметрическом пространстве динамической системы примера 3. Обозначения как для примера 2.

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

  1. Абрамова В.В. "Лекции по динамическим системам и биоматематике", 2023.