Неподвижные точки системы: различия между версиями

Материал из sawiki
Перейти к навигации Перейти к поиску
Строка 1: Строка 1:
== Динамические системы с дискретным временем ==
+
Пусть задана динамическая система:
Пусть задана динамическая система с дискретным временем
 
  
\begin{equation} \label{sist1}
+
\label(sist1)
N_{t+1}=f(N_{t}), N_{t}\in\mathbb{R}, f:\mathbb{R} \rightarrow \mathbb{R}
+
\[
\end{equation}
+
u \mapsto f(u) = f(u;r),[[Участник:Kirich23|Kirich23]] ([[Обсуждение участника:Kirich23|обсуждение]])u \in \mathcal{U} \subset X,[[Участник:Kirich23|Kirich23]] ([[Обсуждение участника:Kirich23|обсуждение]])r \in \mathcal(R), [[Участник:Kirich23|Kirich23]] ([[Обсуждение участника:Kirich23|обсуждение]])f:\mathcal{U} \rightarrow \mathcal{U},
 +
\]
 +
где множество $$X \subset \mathcal{R}^n$$
  
'''Определение 1.'''
+
Определение 1.
 +
'''Неподвижными точками'''
 +
[https://sawiki.cs.msu.ru/index.php/%D0%94%D0%B8%D0%BD%D0%B0%D0%BC%D0%B8%D1%87%D0%B5%D1%81%D0%BA%D0%B0%D1%8F_%D1%81%D0%B8%D1%81%D1%82%D0%B5%D0%BC%D0%B0#.D0.9E.D0.BF.D1.80.D0.B5.D0.B4.D0.B5.D0.BB.D0.B5.D0.BD.D0.B8.D0.B5_.D0.B4.D0.B8.D0.BD.D0.B0.D0.BC.D0.B8.D1.87.D0.B5.D1.81.D0.BA.D0.BE.D0.B9_.D1.81.D0.B8.D1.81.D1.82.D0.B5.D0.BC.D1.8B динамической системы]
 +
(\ref(sist1)) называются такие точки пространства состояний $$u^{∗}$$, что $$f(u^{∗}) = u^{∗}$$.
  
Решения задачи (\ref{sist1}), не изменяющиеся с течением времени $$t$$ называются неподвижными точками отображения (\ref{sist1}).
+
Определение 2.
Неподвижные точки определяются как решение уравнения $$N=f(N)$$.
+
'''пространством состояний''' (или
 +
[https://ru.wikipedia.org/wiki/%D0%A4%D0%B0%D0%B7%D0%BE%D0%B2%D0%BE%D0%B5_%D0%BF%D1%80%D0%BE%D1%81%D1%82%D1%80%D0%B0%D0%BD%D1%81%D1%82%D0%B2%D0%BE фазовым пространством)
 +
называется множество всевозможных состояний $$u_t$$.
  
Заметим, что нередко функцию $$f(N)$$ представляют в виде $$f(N)=Nf(N)$$, чтобы подчеркнуть существование тривиальной неподвижной точки $$N^{*}=0$$.
+
Кроме термина '''"неподвижная точка"''' используют иногда термины '''"стационарная точка"''' или '''"положение равновесия"'''.
В этом случае остальные неподвижные точки — решения уравнения $$f(N)$$.
 
  
 
== Устойчивость неподвижных точек ==
 
== Устойчивость неподвижных точек ==
Строка 18: Строка 23:
 
[[Файл:Fig2staticpoints.jpg|мини|300px|$$r=0.9, b=7$$, система имеет единственный корень $$v=0$$.]]
 
[[Файл:Fig2staticpoints.jpg|мини|300px|$$r=0.9, b=7$$, система имеет единственный корень $$v=0$$.]]
  
'''Определение 2.'''
+
Известно, что никакая реальная система не может находиться все время
 +
в одном и том же состоянии, так как помимо всего прочего подвержена внешним
 +
воздействиям. Что произойдет, если немного возмутить состояние системы так, что
 +
ее состояние окажется в некоторой окрестности положения равновесия? Траектории
 +
могут покинуть эту окрестность, остаться в этой окрестности, или приблизиться к
 +
положению равновесия. Естественно назвать положение равновесия неустойчивым
 +
в первом случае и устойчивым в двух других. Чтобы формализовать обсуждаемое
 +
понятие устойчивости, введем следующее определение.
  
Неподвижная точка $$N^{*}$$ отображения (\ref{sist1}) называется '''устойчивой по Ляпунову''', если для любого $$\varepsilon > 0$$ существует такое $$\delta > 0$$, что
+
'''Определение 3.'''
для любых начальных данных $$N_{0}$$ из $$\delta$$-окрестности точки $$N^{*}$$ вся траектория системы $$N_{t}$$, $$t = 0, 1, 2, ...$$ содержится в $$\varepsilon$$-окрестности точки $$N^{*}$$.
 
  
Если, кроме того, $$\displaystyle{\lim_{t\to\infty}} N_{t} = N^{*}$$, то точка $$N^{*}$$ называется '''асимптотически устойчивой'''.
+
Неподвижная точка $$u^{*}$$ отображения (\ref{sist1}) называется '''устойчивой по Ляпунову''', если для любого $$\varepsilon > 0$$ существует такое $$\delta > 0$$, что
 +
для любых начальных данных $$u_{0}$$ из $$\delta$$-окрестности точки $$u^{*}$$ вся траектория системы $$u_{t}$$, $$t = 0, 1, 2, ...$$ содержится в $$\varepsilon$$-окрестности точки $$u^{*}$$.
  
'''(Определение 2 в символьном виде)'''
+
Если, кроме того, $$\displaystyle{\lim_{t\to\infty}} u_{t} = u^{*}$$, то точка $$u^{*}$$ называется '''асимптотически устойчивой'''.
  
$$N^{*}=f(N^{*})$$ '''устойчивая по Ляпунову''', если $$\forall \varepsilon > 0 ~ \exists \delta(\varepsilon) ~ \forall N_0 \in \mathcal{U}_{\delta}(N^{*})$$ верно $$|N^{*}-N_{t}|< \varepsilon, ~\forall t \geqslant 0$$.
+
То есть:
 +
$$u^{*}=f(u^{*})$$ '''устойчивая по Ляпунову''', если $$\forall \varepsilon > 0 ~ \exists \delta(\varepsilon) ~ \forall N_0 \in \mathcal{U}_{\delta}(u^{*})$$ верно $$|u^{*}-u_{t}|< \varepsilon, ~\forall t \geqslant 0$$.
  
Если, кроме того, $$|N(t; N_0) − N^{*}| \rightarrow 0$$, при $$t \rightarrow \infty$$, то положение равновесия $$N^{*}$$ называется '''асимптотически устойчивым'''.
+
Если, кроме того, $$|u(t; u_0) − u^{*}| \rightarrow 0$$, при $$t \rightarrow \infty$$, то положение равновесия $$u^{*}$$ называется '''асимптотически устойчивым'''.
  
Асимптотически устойчивые неподвижные точки иногда называют аттракторами, а неустойчивые неподвижные точки иногда называют репеллерами.
+
Асимптотически устойчивые неподвижные точки иногда называют
 +
[https://ru.wikipedia.org/wiki/%D0%90%D1%82%D1%82%D1%80%D0%B0%D0%BA%D1%82%D0%BE%D1%80 аттракторами],
 +
а неустойчивые неподвижные точки иногда называют
 +
[репеллерами].
  
 
'''Теорема 1.'''
 
'''Теорема 1.'''
  
Пусть $$N^{*}$$ — неподвижная точка отображения (\ref{sist1}), т.е. $$N^{*}$$ = $$f(N^{*})$$, и пусть $$f$$ обратима в малой окрестности $$N^{*}$$. Тогда $$N^{*}$$ асимптотически устойчива, если
+
Пусть $$u^{*}$$ — неподвижная точка отображения (\ref{sist1}), т.е. $$u^{*}$$ = $$f(u^{*})$$, и пусть $$f$$ обратима в малой окрестности $$u^{*}$$. Тогда $$u^{*}$$ асимптотически устойчива, если
$$|f^{'}(N^{*})| < 1$$,
+
$$|f^{'}(u^{*})| < 1$$,
 
и неустойчива, если
 
и неустойчива, если
$$|f^{'}(N^{*})| > 1$$.
+
$$|f^{'}(u^{*})| > 1$$.
Если $$ |f^{'}(N^{*})| = 1$$, то требуются дополнительные исследования.
+
Если $$ |f^{'}(u^{*})| = 1$$, то требуются дополнительные исследования.
  
 
''Доказательство.''
 
''Доказательство.''
  
Пусть $$|f^{'}(N^{*})| < 1$$ и пусть $$N$$ принадлежит малой окрестности $$N^{*}$$.
+
Пусть $$|f^{'}(u^{*})| < 1$$ и пусть $$u$$ принадлежит малой окрестности $$u^{*}$$.
 
Так как
 
Так как
 
\[
 
\[
\displaystyle{\lim_{t\to\infty}} \frac{|f(N)-f(N^{*})|}{|N-N^{*}|}=|f^{'}(N^{*})|,
+
\displaystyle{\lim_{t\to\infty}} \frac{|f(u)-f(u^{*})|}{|u-u^{*}|}=|f^{'}(u^{*})|,
 
\]
 
\]
поэтому существует такая окрестность $$N^{*}$$, что
+
поэтому существует такая окрестность $$u^{*}$$, что
 
\[
 
\[
\frac{|f(N)-f(N^{*})|}{|N-N^{*}|} \leqslant a,
+
\frac{|f(u)-f(u^{*})|}{|u-u^{*}|} \leqslant a,
 
\]
 
\]
для всех $$N$$ из этой окрестности; здесь $$a$$ — некоторое число, такое что $$|f^{'}(N^{*})| \leqslant a < 1$$.
+
для всех $$u$$ из этой окрестности; здесь $$a$$ — некоторое число, такое что $$|f^{'}(u^{*})| \leqslant a < 1$$.
Таким образом, $$f(N)$$ остается в той же окрестности, что и $$N$$, и, кроме того, ближе к
+
Таким образом, $$f(u)$$ остается в той же окрестности, что и $$u$$, и, кроме того, ближе к
неподвижной точке $$N^{*}$$, по крайней мере, на множитель $$a$$. Отсюда следует, что
+
неподвижной точке $$u^{*}$$, по крайней мере, на множитель $$a$$. Отсюда следует, что
\[|f(f(N)) − f(f(N^{*}))| \leqslant a |f(N) − f(N^{*})| \leqslant a^2|N N^{*}|,
+
\[|f(f(u)) − f(f(u^{*}))| \leqslant a |f(u) − f(u^{*})| \leqslant a^2|u u^{*}|,
 
\]
 
\]
 
или, по индукции,
 
или, по индукции,
\[|f^{k}(N) − N^{*}|  \leqslant a^{k}|N N^{*}|,
+
\[|f^{k}(u) − u^{*}|  \leqslant a^{k}|u u^{*}|,
 
\]
 
\]
где $$f^{k}$$ обозначает $$k$$-ую суперпозицию отображения $$f$$. Таким образом мы доказали,что последовательность $$f^{k}(N)$$ будет сходиться к $$N^{*}$$, то есть является асимптотически устойчивой.
+
где $$f^{k}$$ обозначает $$k$$-ую суперпозицию отображения $$f$$. Таким образом мы доказали, что последовательность $$f^{k}(u)$$ будет сходиться к $$u^{*}$$, то есть является асимптотически устойчивой.
  
 
Вторая часть утверждения доказывается сходным образом.
 
Вторая часть утверждения доказывается сходным образом.
Строка 114: Строка 130:
 
Таким образом, по теореме 1, точка $$v^{*}_1=0$$ будет асимптотически устойчивой, при $$r<1$$ и неустойчивой при $$r>1$$.
 
Таким образом, по теореме 1, точка $$v^{*}_1=0$$ будет асимптотически устойчивой, при $$r<1$$ и неустойчивой при $$r>1$$.
  
== Динамические системы с непрерывным временем ==
 
 
'''Определение 3.'''
 
 
Точка $$a \in \mathbb{R}^n$$ называется неподвижной точкой динамической системы $$\dot x_i = f_i(x_1,...,x_n)$$, где $$(x_1,...,x_2) \in D \subset \mathbb{R}^n, i=\overline{1,n}, f=(f_1,...,f_n)$$, если $$f(a)=0$$.
 
 
 
Пусть задана динамическая система с непрерывным временем
 
 
\begin{equation} \label{sist3}
 
\dot u=f(u), u \in U \subseteq \mathbb{R}^{n}, f: U \rightarrow \mathbb{R}^{n}.
 
\end{equation}
 
 
'''Определение 4.'''
 
 
Положение равновесия динамической системы (\ref{sist3}) называется
 
гиперболическим, если число собственных значений $$n_0$$ якобиана вектор-функции $$f(u)$$ равняется нулю,
 
т.е. не существует собственных чисел, расположенных на мнимой оси матрицы Якоби.
 
Гиперболическое положение равновесия называется гиперболическим седлом,
 
если $$n_+$$ $$n_- \neq 0$$, где $$n_+$$, $$n_-$$, $$n_0$$ — количество собственных
 
значений якобиана (с учётом их кратности) с положительной, отрицательной и равной нулю
 
вещественной частью соответственно.
 
 
'''Теорема 2. (А. М. Ляпунов, А. Пуанкаре)'''
 
 
Пусть $$u^*$$ — гиперболическое положение равновесия (\ref{sist3}). Тогда, если $$n_+=0$$, то положение равновесия $$u^*$$ асимптотически устойчиво, если $$n_+>0$$, то неустойчиво.
 
 
Устойчивость гиперболического положения равновесия определяется собственными числами матрицы Якоби, вычисленными в точке $$u^*$$.
 
В случае, когда вещественная часть одного или нескольких собственных значений
 
равна нулю, необходимы дополнительные исследования.
 
 
== Пример поиска и исследования на устойчивость неподвижных точек для системы с непрерывным временем==
 
[[Файл:Fig7staticpoints.jpg|мини|300px|Параметрический портрет системы (\ref{sist44})]]
 
[[Файл:Fig8staticpoints.jpg|мини|300px|Поведение траекторий вблизи положения равновесия $$(1, 16.5)$$ при значениях $$\mu = 1, \nu=33$$. Построен неустойчивый узел.]]
 
[[Файл:Fig9staticpoints.jpg|мини|300px|Поведение траекторий вблизи положения равновесия $$(5, 2.5)$$ при значениях $$\mu = 5, \nu=0.6$$. Построен неустойчивый фокус.]]
 
Рассмотрим динамическую систему
 
 
\begin{equation} \label{sist44}
 
\begin{cases}
 
\dfrac{du}{d\tau}= \dfrac{\mu \nu u^2}{1+u}-uv,\\
 
\dfrac{dv}{d\tau}= v(-\mu + u),
 
\end{cases}
 
\mu, \nu > 0.
 
\end{equation}
 
 
Найдем неподвижные точки получившейся системы (\ref{sist4}). Для этого необходимо решить
 
 
\begin{equation} \label{sist4}
 
\begin{cases}
 
0=\dfrac{\mu \nu u^2}{1+u}-uv,\\
 
0=v(-\mu + u).
 
\end{cases}
 
\Longrightarrow
 
\begin{cases}
 
uv=\dfrac{\mu \nu u^2}{1+u},\\
 
uv=v \mu.
 
\end{cases}
 
\end{equation}
 
 
Из второго уравнения системы (\ref{sist4}):
 
 
Пусть $$v=0$$.
 
Тогда $$\dfrac{\mu \nu u^2}{1+u}=0, \mu,\nu>0$$, значит $$u=0$$.
 
   
 
Пусть $$v \neq 0$$.
 
Тогда $$u = \mu$$ и $$v = \dfrac{\mu^2 \nu}{1+\mu}$$.
 
 
Других решений система (\ref{sist4}) не имеет.
 
 
Получили 2 неподвижные точки:
 
 
1) $$u=0,v=0$$. Обозначим как точку $$A$$.
 
 
2) $$u=\mu, v=\frac{\mu^2 \nu}{1+\mu}$$. Обозначим как точку $$B$$.
 
 
Чтобы исследовать найденные неподвижные точки на устойчивость, выпишем [https://ru.wikipedia.org/wiki/%D0%AF%D0%BA%D0%BE%D0%B1%D0%B8%D0%B0%D0%BD якобиан]
 
системы (\ref{sist4}) и воспользуемся теоремой 2.
 
 
\[
 
    J(u,v) = \begin{bmatrix}
 
        \frac{2 \mu \nu u}{1+u} - \frac{\mu \nu u^2}{(1+u)^2} - v & -u \\
 
  v & u-\mu
 
\end{bmatrix}
 
\]
 
 
Рассмотрим значения якобиана в точке $$A=(0,0)$$.
 
 
\[
 
J(A) = \begin{bmatrix}
 
        0 & 0 \\
 
  0 & -\mu
 
\end{bmatrix}\ 
 
\]
 
 
Собственными значениями являются числа $$\lambda_1=0, \lambda_2=-\mu$$. В этом случае теорема 2 не применима. Устойчивость точки $$A=(0,0)$$ без дополнительных исследований установить нельзя.
 
 
Рассмотрим значения якобиана в точке $$B=\left(\mu,\frac{\mu^2 \nu}{1+\mu}\right)$$.
 
 
\[
 
J(B) = \begin{bmatrix}
 
        \frac{\mu^2 \nu}{\mu+1}-\frac{\mu^3 \nu}{(\mu+1)^2} & - \mu \\
 
  \frac{\mu^2 \nu}{\mu+1} & 0
 
\end{bmatrix}\
 
\]
 
 
Найдём собственные значения.
 
 
\[
 
\det(J(B)- I \cdot \lambda)
 
=
 
\begin{bmatrix}
 
        \frac{\mu^2 \nu}{\mu+1}-\frac{\mu^3 \nu}{(\mu + 1 )^2} - \lambda & - \mu \\
 
  \frac{\mu^2 \nu}{\mu+1} & - \lambda
 
\end{bmatrix}\
 
= -
 
\lambda
 
\left(
 
\frac{\mu^2 \nu}{\mu+1} - \frac{\mu^3 \nu}{(\mu+1)^2 }
 
- \lambda
 
\right) +
 
\frac{\mu^3 \nu}{\mu+1}
 
=0.
 
\]
 
 
\[
 
(\mu+1)^2 \lambda^2 - \mu^2 \nu \lambda + (\mu + 1) \mu^3 \nu =0.
 
\]
 
 
Далее решаем квадратное уравнение относительно $$\lambda$$.
 
Дискриминант:
 
\[
 
D = (\mu^2 \nu)^2 - 4(\mu + 1)^2(\mu+1) \mu^3 \nu =  \mu^4 \nu^2 - 4\mu^3 \nu (\mu^3 + 3\mu^2 + 3\mu +1).
 
\]
 
 
Тогда корни уравнения:
 
\[
 
\lambda_{1,2} =
 
\frac{\mu^2 \nu \pm \sqrt{\mu^4 \nu^2 - 4 \mu^3 \nu (\mu^3 + 3\mu^2 + 3\mu +1)}}{2(\mu+1)^2}
 
\]
 
 
Так как параметры $$\mu, \nu>0$$, то положения равновесия являются гиперболическими. Возможны 2 случая:
 
 
'''1) $$D \geqslant 0$$.'''
 
 
Этот случай возникает при $$\nu \geqslant 4(\mu^2 + 3\mu + 3 + \frac{1}{\mu})$$. Собственные значения $$\lambda_{1,2}$$ будут вещественными и положительными, таким образом положение равновесия — узел.
 
 
По теореме 2 узел будет неустойчивым.
 
 
'''2) $$D < 0$$.'''
 
 
Этот случай возникает при $$\nu < 4(\mu^2 + 3\mu + 3 + \frac{1}{\mu})$$. Собственные значения $$\lambda_{1,2}$$ будут комплексными с положительной вещественной частью, значит положение равновесия — фокус.
 
 
По теореме 2 фокус будет неустойчивым.
 
 
Проанализируем параметрический портрет системы.
 
 
В области $$I$$ точка $$B$$ будет неустойчивым узлом.
 
В области $$II$$ точка $$B$$ будет неустойчивым фокусом.
 
  
== формулировки основных теорем из функционального анализа о неподвижных точках ==
+
== Связанные теоремы ==
  
'''Теорема 3. (Банах)'''
+
'''Теорема 2. (Банах)'''
  
 
Пусть $$(X, d)$$ — полное метрическое пространство с метрикой $$d$$.
 
Пусть $$(X, d)$$ — полное метрическое пространство с метрикой $$d$$.
Строка 286: Строка 144:
 
точке $$\xi$$.
 
точке $$\xi$$.
  
'''Теорема 4. (Брауэр)'''
+
'''Теорема 3. (Брауэр)'''
  
 
Любое непрерывное отображение замкнутого шара в себя в конечномерном евклидовом пространстве имеет неподвижную точку.
 
Любое непрерывное отображение замкнутого шара в себя в конечномерном евклидовом пространстве имеет неподвижную точку.
Строка 292: Строка 150:
 
Более подробно, рассмотрим замкнутый шар в ''n''-мерном пространстве <math>B^n\subset \mathbb R^n</math>. Пусть <math>f \colon B^n\to B^n</math> — некоторое непрерывное отображение этого шара в себя (не обязательно строго внутрь себя, не обязательно биективное, т.е. даже не обязательно сюръективное). Тогда найдется такая точка <math>x\in B^n</math>, что <math>f(x)=x</math>.
 
Более подробно, рассмотрим замкнутый шар в ''n''-мерном пространстве <math>B^n\subset \mathbb R^n</math>. Пусть <math>f \colon B^n\to B^n</math> — некоторое непрерывное отображение этого шара в себя (не обязательно строго внутрь себя, не обязательно биективное, т.е. даже не обязательно сюръективное). Тогда найдется такая точка <math>x\in B^n</math>, что <math>f(x)=x</math>.
  
'''Теорема 5. (Шаудер — Тихонов)'''
+
'''Теорема 4. (Шаудер — Тихонов)'''
  
 
В локально выпуклом топологическом векторном пространстве любое непрерывное отображение $$f : K\to K$$ выпуклого компактного множества $$K$$ в себя имеет неподвижную точку.
 
В локально выпуклом топологическом векторном пространстве любое непрерывное отображение $$f : K\to K$$ выпуклого компактного множества $$K$$ в себя имеет неподвижную точку.

Версия 22:47, 16 декабря 2023

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

\label(sist1) \[ u \mapsto f(u) = f(u;r),[[Участник:Kirich23|Kirich23]] ([[Обсуждение участника:Kirich23|обсуждение]])u \in \mathcal{U} \subset X,[[Участник:Kirich23|Kirich23]] ([[Обсуждение участника:Kirich23|обсуждение]])r \in \mathcal(R), [[Участник:Kirich23|Kirich23]] ([[Обсуждение участника:Kirich23|обсуждение]])f:\mathcal{U} \rightarrow \mathcal{U}, \] где множество $$X \subset \mathcal{R}^n$$

Определение 1. Неподвижными точками динамической системы (\ref(sist1)) называются такие точки пространства состояний $$u^{∗}$$, что $$f(u^{∗}) = u^{∗}$$.

Определение 2. пространством состояний (или [https://ru.wikipedia.org/wiki/%D0%A4%D0%B0%D0%B7%D0%BE%D0%B2%D0%BE%D0%B5_%D0%BF%D1%80%D0%BE%D1%81%D1%82%D1%80%D0%B0%D0%BD%D1%81%D1%82%D0%B2%D0%BE фазовым пространством) называется множество всевозможных состояний $$u_t$$.

Кроме термина "неподвижная точка" используют иногда термины "стационарная точка" или "положение равновесия".

Устойчивость неподвижных точек

$$r=30, b=3$$, система имеет два различных неотрицательных корня.
$$r=0.9, b=7$$, система имеет единственный корень $$v=0$$.

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

Определение 3.

Неподвижная точка $$u^{*}$$ отображения (\ref{sist1}) называется устойчивой по Ляпунову, если для любого $$\varepsilon > 0$$ существует такое $$\delta > 0$$, что для любых начальных данных $$u_{0}$$ из $$\delta$$-окрестности точки $$u^{*}$$ вся траектория системы $$u_{t}$$, $$t = 0, 1, 2, ...$$ содержится в $$\varepsilon$$-окрестности точки $$u^{*}$$.

Если, кроме того, $$\displaystyle{\lim_{t\to\infty}} u_{t} = u^{*}$$, то точка $$u^{*}$$ называется асимптотически устойчивой.

То есть: $$u^{*}=f(u^{*})$$ устойчивая по Ляпунову, если $$\forall \varepsilon > 0 ~ \exists \delta(\varepsilon) ~ \forall N_0 \in \mathcal{U}_{\delta}(u^{*})$$ верно $$|u^{*}-u_{t}|< \varepsilon, ~\forall t \geqslant 0$$.

Если, кроме того, $$|u(t; u_0) − u^{*}| \rightarrow 0$$, при $$t \rightarrow \infty$$, то положение равновесия $$u^{*}$$ называется асимптотически устойчивым.

Асимптотически устойчивые неподвижные точки иногда называют аттракторами, а неустойчивые неподвижные точки иногда называют [репеллерами].

Теорема 1.

Пусть $$u^{*}$$ — неподвижная точка отображения (\ref{sist1}), т.е. $$u^{*}$$ = $$f(u^{*})$$, и пусть $$f$$ обратима в малой окрестности $$u^{*}$$. Тогда $$u^{*}$$ асимптотически устойчива, если $$|f^{'}(u^{*})| < 1$$, и неустойчива, если $$|f^{'}(u^{*})| > 1$$. Если $$ |f^{'}(u^{*})| = 1$$, то требуются дополнительные исследования.

Доказательство.

Пусть $$|f^{'}(u^{*})| < 1$$ и пусть $$u$$ принадлежит малой окрестности $$u^{*}$$. Так как \[ \displaystyle{\lim_{t\to\infty}} \frac{|f(u)-f(u^{*})|}{|u-u^{*}|}=|f^{'}(u^{*})|, \] поэтому существует такая окрестность $$u^{*}$$, что \[ \frac{|f(u)-f(u^{*})|}{|u-u^{*}|} \leqslant a, \] для всех $$u$$ из этой окрестности; здесь $$a$$ — некоторое число, такое что $$|f^{'}(u^{*})| \leqslant a < 1$$. Таким образом, $$f(u)$$ остается в той же окрестности, что и $$u$$, и, кроме того, ближе к неподвижной точке $$u^{*}$$, по крайней мере, на множитель $$a$$. Отсюда следует, что \[|f(f(u)) − f(f(u^{*}))| \leqslant a |f(u) − f(u^{*})| \leqslant a^2|u − u^{*}|, \] или, по индукции, \[|f^{k}(u) − u^{*}| \leqslant a^{k}|u − u^{*}|, \] где $$f^{k}$$ обозначает $$k$$-ую суперпозицию отображения $$f$$. Таким образом мы доказали, что последовательность $$f^{k}(u)$$ будет сходиться к $$u^{*}$$, то есть является асимптотически устойчивой.

Вторая часть утверждения доказывается сходным образом. $$\blacksquare$$

Пример поиска неподвижных точек с помощью графического метода

Аттрактор в точке $$v_2^*$$ при параметрах $$r=10, b=2, v_0=1.$$
Аттрактор в точке $$v_2^*$$ при параметрах $$r=3, b=0.9, v_0=4.5$$.

Пусть задана система $$v_{t+1} = \frac{r v_t}{(1+ v_t)^b}$$.

Графически неподвижные точки — это точки пересечения графика функции $$f(N)$$ и биссектрисы первого координатного угла $$N_{t+1} = N_{t}$$ (нас интересуют только неотрицательные решения). Для нахождения неподвижных точек заданной системы рассмотрим возможные пересечения графика функции $$f(v)=\frac{rv}{(1+v)^b}$$ с прямой $$g(v)=v$$.

Заметим, что система \begin{equation} \label{sist2} \begin{cases} f(v)=\frac{rv}{(1+v)^b},\\g(v)=v. \end{cases} \end{equation} при любых $$v \geqslant 0, ~r,b > 0$$ имеет хотя бы один корень, а именно $$v=0$$. Второй корень $$v=r^{1/b}-1$$ система может иметь только при значении параметра $$r> 1$$. Больше двух неотрицательных корней система (\ref{sist2}) не имеет.


Аттрактор в точке $$v_1^*$$ при параметрах $$r=0.9, b=0.7, v_0=4.$$
Репеллер в точке $$v_1^*$$ при параметрах $$r=2, b=0.9, v_0=2.$$

Пример исследования неподвижных точек на устойчивость для системы с дискретным временем

Исследуем на устойчивость неподвижные точки $$v^{*}_1=0$$ и $$v^{*}_2=r^{1/b}-1$$ для системы (\ref{sist2}).

\begin{equation} f_v(v^{*})=\frac{r((1+v)^b-bv(1+v)^{b-1})}{(1+v)^{2b}}. \end{equation}

Исследуем сначала точку $$v^{*}_2=r^{1/b}-1$$:

Подставим $$v_2^{*}$$ в выражение (\ref{sist2}) и с учетом наложенных ограничений $$r>1, b>0$$ для существования точки получим

\[ f_v(v_2^{*})=b(r^{-1/b}-1)+1. \]

Согласно теореме 1 точка $$v^{*}_2=r^{1/b}-1$$ будет асимптотически устойчивой, при $$r^{-1/b}<1$$ и неустойчивой при $$r^{-1/b}>1$$. Отметим, что точка $$v^{*}_2=r^{1/b}-1$$ (отличная от нуля) существует только при значении параметра $$r > 1$$, таким образом $$r^{-1/b}$$ принимает значения меньше единицы при любых $$r> 1$$, значит, точка $$v^{*}_2$$ асимптотически устойчива всегда, если она существует.

Теперь исследуем на устойчивость точку $$v^{*}_1=0$$:

Подставим $$v_1^{*}$$ в выражение (\ref{sist2}) и получим \[ f_v(v_1^{*})=r. \] Таким образом, по теореме 1, точка $$v^{*}_1=0$$ будет асимптотически устойчивой, при $$r<1$$ и неустойчивой при $$r>1$$.


Связанные теоремы

Теорема 2. (Банах)

Пусть $$(X, d)$$ — полное метрическое пространство с метрикой $$d$$.

Пусть задано отображение $$f : X \rightarrow X$$ и существует число $$a, 0 \leqslant a < 1$$, такое, что для любых $$x, y \in X: d(f(x), f(y)) \leqslant a \cdot d(x, y)$$.

Тогда существует единственная точка $$\xi \in X$$ такая, что $$f(\xi) = \xi$$, и начиная с любой точки $$x_0 \in X$$, последовательность итераций $${f_n(x_0)}_{n=1,2,...}$$ сходится к точке $$\xi$$.

Теорема 3. (Брауэр)

Любое непрерывное отображение замкнутого шара в себя в конечномерном евклидовом пространстве имеет неподвижную точку.

Более подробно, рассмотрим замкнутый шар в n-мерном пространстве \(B^n\subset \mathbb R^n\). Пусть \(f \colon B^n\to B^n\) — некоторое непрерывное отображение этого шара в себя (не обязательно строго внутрь себя, не обязательно биективное, т.е. даже не обязательно сюръективное). Тогда найдется такая точка \(x\in B^n\), что \(f(x)=x\).

Теорема 4. (Шаудер — Тихонов)

В локально выпуклом топологическом векторном пространстве любое непрерывное отображение $$f : K\to K$$ выпуклого компактного множества $$K$$ в себя имеет неподвижную точку.

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

1. Братусь А.С., Новожилов А.С., Платонов А.П. Динамические системы и модели биологии, 2011.

2. Абрамова В.В. Лекции по курсу "Динамические системы и биоматематика", 2023.

3. Banach S. Sur les opérations dans les ensembles abstraits et leur applications aux équations integrales, Fund. Math, 1922

4. Шашкин Ю.А. Неподвижные точки, М.: Наука, 1989

5. J. Schauder, Der Fixpunktsatz in Funktionalräumen, Studia Math. 2, 1930