Повторить

Теоремы Ляпунова и Четаева. Устойчивость по первому приближению

Прогресс изучения0 из 6

Отмечай теоремы кнопками ✓ (изучено) и ↻ (повторить)

Все теоремы темы6
  1. Теорема Ляпунова об устойчивости
  2. Теорема Ляпунова об асимптотической устойчивости
  3. Теорема Четаева о неустойчивости
  4. Квадратичная функция Ляпунова для асимптотически устойчивой линейной системы
  5. Устойчивость по первому приближению
  6. Критический случай первого приближения
теорема

Теорема Ляпунова об устойчивости

Формулировка

Пусть f∈C1f\in C^1, f(0)=0f(0) = 0, и в окрестности нуля существует положительно определённая функция V∈C1V\in C^1, производная которой в силу системы x′=f(x)x' = f(x) неположительна: V˙≤0\dot V\le 0. Тогда нулевое решение устойчиво по Ляпунову.

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

Развернуть

Пусть ε>0\varepsilon > 0 так мало, что шар ∥x∥≤ε\|x\|\le\varepsilon лежит в окрестности, где определена VV. На сфере ∥x∥=ε\|x\| = \varepsilon непрерывная положительная функция VV достигает минимума m>0m > 0. Так как V(0)=0V(0) = 0 и VV непрерывна, найдётся δ∈(0,ε)\delta\in(0, \varepsilon) с V(x)<mV(x) < m при ∥x∥<δ\|x\| < \delta.

Пусть ∥x(t0)∥<δ\|x(t_0)\| < \delta. Пока решение остаётся в шаре ∥x∥≤ε\|x\|\le\varepsilon, функция V(x(t))V(x(t)) не возрастает, так что V(x(t))≤V(x(t0))<mV(x(t))\le V(x(t_0)) < m. Поэтому решение не может достичь сферы ∥x∥=ε\|x\| = \varepsilon, где V≥mV\ge m.

Значит, решение остаётся в компакте ∥x∥<ε\|x\| < \varepsilon. Решение, не покидающее компакт внутри области, продолжается на всю полуось t≥t0t\ge t_0 (иначе оно должно было бы приблизиться к границе области). Итак, ∥x(t)∥<ε\|x(t)\| < \varepsilon при всех t≥t0t\ge t_0. ■\blacksquare

Источники:Эльсгольц, Дифференциальные уравнения и вариационное исчисление, гл. 4Филиппов, Введение в теорию дифференциальных уравнений, гл. 4
теорема

Теорема Ляпунова об асимптотической устойчивости

Формулировка

Если в условиях теоремы Ляпунова об устойчивости производная V˙\dot V отрицательно определена, то нулевое решение асимптотически устойчиво.

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

Развернуть

По теореме Ляпунова об устойчивости нуль устойчив. Возьмём δ0\delta_0, отвечающее некоторому ε\varepsilon, и решение с ∥x(t0)∥<δ0\|x(t_0)\| < \delta_0; тогда ∥x(t)∥<ε\|x(t)\| < \varepsilon при t≥t0t\ge t_0.

Функция V(x(t))V(x(t)) убывает и ограничена снизу нулём, поэтому имеет предел L≥0L\ge 0. Пусть L>0L > 0. Так как V(0)=0V(0) = 0, найдётся η>0\eta > 0 с V<LV < L при ∥x∥<η\|x\| < \eta, и, значит, η≤∥x(t)∥≤ε\eta\le\|x(t)\|\le\varepsilon при всех t≥t0t\ge t_0. На этом кольце непрерывная отрицательная функция V˙\dot V не превосходит −γ<0-\gamma < 0. Тогда

V(x(t))≤V(x(t0))−γ(t−t0)→−∞,V(x(t))\le V(x(t_0)) - \gamma(t - t_0)\to -\infty,

что противоречит V≥0V\ge 0. Итак, V(x(t))→0V(x(t))\to 0.

Отсюда x(t)→0x(t)\to 0: если бы ∥x(tk)∥≥η′>0\|x(t_k)\|\ge\eta' > 0 для последовательности tk→∞t_k\to\infty, то V(x(tk))≥min⁡η′≤∥x∥≤εV>0V(x(t_k))\ge\min_{\eta'\le\|x\|\le\varepsilon}V > 0. ■\blacksquare

Источники:Эльсгольц, Дифференциальные уравнения и вариационное исчисление, гл. 4Филиппов, Введение в теорию дифференциальных уравнений, гл. 4
теорема

Теорема Четаева о неустойчивости

Формулировка

Пусть f∈C1f\in C^1, f(0)=0f(0) = 0, и в шаре ∥x∥≤r\|x\|\le r определена функция V∈C1V\in C^1 такая, что:

  1. в любой окрестности нуля есть точки, где V>0V > 0;
  2. V˙>0\dot V > 0 во всех точках шара, где V>0V > 0.

Тогда нулевое решение системы x′=f(x)x' = f(x) неустойчиво.

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

Развернуть

Возьмём ε=r\varepsilon = r и любое δ>0\delta > 0. По условию 1 есть точка x0x^0 с ∥x0∥<δ\|x^0\| < \delta и V(x0)=α>0V(x^0) = \alpha > 0. Пусть x(t)x(t) — решение с x(t0)=x0x(t_0) = x^0.

Предположим, что ∥x(t)∥<r\|x(t)\| < r при всех t≥t0t\ge t_0. Пока V(x(t))>0V(x(t)) > 0, функция V(x(t))V(x(t)) возрастает, поэтому V(x(t))≥αV(x(t))\ge\alpha при всех t≥t0t\ge t_0. Значит, решение лежит в множестве

K={∥x∥≤r, V(x)≥α}.K = \{\|x\|\le r,\ V(x)\ge\alpha\}.

Оно компактно, и на нём V˙>0\dot V > 0 по условию 2, поэтому V˙≥β>0\dot V\ge\beta > 0 на KK. Тогда V(x(t))≥α+β(t−t0)→+∞V(x(t))\ge\alpha + \beta(t - t_0)\to+\infty, что противоречит ограниченности непрерывной функции VV на шаре.

Итак, при сколь угодно малом δ\delta найдётся решение, выходящее на сферу ∥x∥=r\|x\| = r. Нулевое решение неустойчиво. ■\blacksquare

Источники:Эльсгольц, Дифференциальные уравнения и вариационное исчисление, гл. 4Филиппов, Введение в теорию дифференциальных уравнений, гл. 4
лемма

Квадратичная функция Ляпунова для асимптотически устойчивой линейной системы

Формулировка

Если все собственные значения матрицы AA имеют отрицательные вещественные части, то матричное уравнение Ляпунова

ATB+BA=−EA^TB + BA = -E

имеет единственное решение BB, и это симметричная положительно определённая матрица:

B=∫0∞etATetA dt.B = \int_0^\infty e^{tA^T}e^{tA}\,dt.

Для квадратичной формы V(x)=xTBxV(x) = x^TBx производная в силу системы x′=Axx' = Ax равна V˙=−∥x∥2\dot V = -\|x\|^2.

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

Развернуть

Элементы etAe^{tA} — комбинации функций tjeλtt^je^{\lambda t} с Re⁡λ<0\operatorname{Re}\lambda < 0, поэтому ∥etA∥≤Me−γt\|e^{tA}\|\le Me^{-\gamma t} при некоторых MM, γ>0\gamma > 0. Интеграл сходится.

BB симметрична, и xTBx=∫0∞∥etAx∥2 dt>0x^TBx = \int_0^\infty\|e^{tA}x\|^2\,dt > 0 при x≠0x\ne 0, так как etAx≠0e^{tA}x\ne 0.

Далее,

ATB+BA=∫0∞ddt(etATetA) dt=[etATetA]0∞=−E.A^TB + BA = \int_0^\infty\frac{d}{dt}\bigl(e^{tA^T}e^{tA}\bigr)\,dt = \Bigl[e^{tA^T}e^{tA}\Bigr]_0^\infty = -E.

Единственность: если ATB0+B0A=0A^TB_0 + B_0A = 0, то ddt(etATB0etA)=0\frac{d}{dt}\bigl(e^{tA^T}B_0e^{tA}\bigr) = 0, так что etATB0etA≡B0e^{tA^T}B_0e^{tA}\equiv B_0. Но левая часть стремится к нулю, значит, B0=0B_0 = 0.

Наконец, V˙=(Ax)TBx+xTBAx=xT(ATB+BA)x=−xTx\dot V = (Ax)^TBx + x^TBAx = x^T(A^TB + BA)x = -x^Tx. ■\blacksquare

Источники:Эльсгольц, Дифференциальные уравнения и вариационное исчисление, гл. 4Филиппов, Введение в теорию дифференциальных уравнений, гл. 4
теорема

Устойчивость по первому приближению

Формулировка

Пусть x′=Ax+g(x)x' = Ax + g(x), g(x)=o(∥x∥)g(x) = o(\|x\|) при x→0x\to 0.

  1. Если все собственные значения AA имеют Re⁡λ<0\operatorname{Re}\lambda < 0, то нулевое решение асимптотически устойчиво.
  2. Если хотя бы одно собственное значение имеет Re⁡λ>0\operatorname{Re}\lambda > 0, то нулевое решение неустойчиво.

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

Развернуть
  1. Возьмём V=xTBxV = x^TBx из леммы о квадратичной функции Ляпунова. Для полной системы
V˙=xT(ATB+BA)x+2xTBg(x)=−∥x∥2+2xTBg(x).\dot V = x^T(A^TB + BA)x + 2x^TBg(x) = -\|x\|^2 + 2x^TBg(x).

Выберем ρ>0\rho > 0 так, чтобы ∥g(x)∥≤∥x∥4∥B∥\|g(x)\|\le\frac{\|x\|}{4\|B\|} при ∥x∥≤ρ\|x\|\le\rho. Тогда ∣2xTBg(x)∣≤12∥x∥2|2x^TBg(x)|\le\frac12\|x\|^2 и V˙≤−12∥x∥2\dot V\le -\frac12\|x\|^2 — отрицательно определена. По теореме Ляпунова об асимптотической устойчивости нуль асимптотически устойчив.

  1. В курсе приводится без доказательства. Идея: в базисе, разделяющем собственные значения с Re⁡λ>0\operatorname{Re}\lambda > 0 и остальные, строится квадратичная форма V=V+(x)−V−(x)V = V_+(x) - V_-(x), удовлетворяющая условиям теоремы Четаева с учётом малости gg. ■\blacksquare
Источники:Эльсгольц, Дифференциальные уравнения и вариационное исчисление, гл. 4Филиппов, Введение в теорию дифференциальных уравнений, гл. 4
лемма

Критический случай первого приближения

Формулировка

Если max⁡Re⁡λ=0\max\operatorname{Re}\lambda = 0, то первое приближение не решает вопроса об устойчивости. Например, для уравнения x′=ax3x' = ax^3 первое приближение x′=0x' = 0 одно и то же при всех aa, но нулевое решение асимптотически устойчиво при a<0a < 0, неустойчиво при a>0a > 0 и устойчиво (не асимптотически) при a=0a = 0.

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

Развернуть

При a<0a < 0: V=x2V = x^2, V˙=2ax4<0\dot V = 2ax^4 < 0 при x≠0x\ne 0, и по теореме Ляпунова нуль асимптотически устойчив.

При a>0a > 0: решение x=x01−2ax02(t−t0)x = \dfrac{x_0}{\sqrt{1 - 2ax_0^2(t - t_0)}} уходит на бесконечность за конечное время при любом x0≠0x_0\ne 0.

При a=0a = 0: решения постоянны. ■\blacksquare

Источники:Эльсгольц, Дифференциальные уравнения и вариационное исчисление, гл. 4Филиппов, Введение в теорию дифференциальных уравнений, гл. 4