35) Метод Эйлера интегрирования линейной однородной системы с постоянной матрицей. Случай различных корней. Для решения системы (где ) или,в векторной записи, ,где -вектор, a-матрица:

надо найти корни характеристического уравнения

36) Метод Эйлера интегрирования линейной однородной системы с постоянной матрицей. Случай кратных корней. Для решения системы (где ) или,в векторной записи, ,где -вектор, a-матрица:

надо найти корни характеристического уравнения

37) Автономные системы. Свойства.

Автономной системой дифференциальных уравнений n –го порядка называется система, которая в нормальной форме записывается в виде

В векторной форме автономная система имеет вид x' = F(x) (не зависит от t), где

Название автономная система связано с тем, что поскольку производная  x' зависит только от x и не зависит от t, то решение само управляет своим изменением. Автономные системы называют также динамическими системами.

Любую систему дифференциальных уравнений, записанную в нормальной форме, можно свести к автономной системе, увеличив число неизвестных функций на единицу:

Будем полагать, что для рассматриваемых автономных систем выполнены условия теоремы существования и единственности решения задачи Коши.

Пусть x = φ(t) — решение автономной системы, определенное на отрезке [a, b]. Множество точек  x = φ(t) ,   t ∈ [a, b] — кривая в пространствеR_xⁿ . Эту кривую называют фазовой траекторией или просто траекторией системы, а пространство R_xⁿ , в котором расположены фазовые траектории, называют фазовым пространством автономной системы.

Точка a называется положением равновесия (точкой покоя) автономной системы, если F(a) = 0 .

Равенство x = φ(t) ,   t ∈ [a, b] — параметрические уравнения фазовой траектории.

Интегральная кривая системы изображается в (n + 1) –мерном пространстве R_x, tⁿ⁺¹ и может быть определена уравнениями

Ясно, что соответствующая фазовая траектория — проекция интегральной кривой на пространство R_x.

Свойства: Если   - решение автономной системы дифференциальных уравнений (в векторном виде), то эта функция остаётся решением и при сдвиге аргумента. Автономная система моделирует автономные процессы, т.е. процесс, не подверженные внешним влияниям, и стационарные процессы, т.е. процессы, установившиеся во времени. Все эти процессы полностью определяются начальными значениями переменных состояния, т.е.  , и не зависят от выбора начального значения аргумента  .

38) Положения равновесия. Циклы.

39) Особые точки. Узлы, центр, седло.

40) Основные понятия устойчивости по Ляпунову.

41) Устойчивость линейных систем. Для линейной системы

x′ = A(t)x + b(t),

(ЛС)

aij, bi ∈ C([t0, +∞), R),

и любого ее решения x = φ(t) приведенная система совпадает с соответствующей (ЛОС).

Д о к а з а т е л ь с т в о.  Произведем замену y = x – φ(t):

y′ = A(t)x + b(t) – A(t)φ(t) – b(t) = A(t)(x – φ(t)) = A(t)y.
	Критерии устойчивости (ЛС). Пусть Φt0(t) — фундаментальная матрица (ЛОС), нормальная в t0. Утверждается, что

(а) (ЛС) устойчива ⇔ Φt0(t) ограничена на [t0, +∞);

(б) (ЛС) асимптотически устойчива ⇔ Φt0(t) → 0 при t → +∞ ⇔ (ЛС) асимптотически устойчива в целом;

(в) (ЛС) экспоненциально устойчива ⇔ (M > 0, γ > 0) ∀ (t ≥ t0) [||Φt0(t)|| ≤ Me–γ(t–t0)] ⇔ (ЛС) экспоненциально устойчива в целом.

Д о к а з а т е л ь с т в о.  (а) Пусть (ЛС) устойчива, т. е. устойчиво нулевое решение (ЛОС). Положив в определении устойчивости ε = 1, найдем δ > 0 такое, что

||x0|| < δ ⇒ ||gt0t(x0)||= ||Φt0(t)x0|| < 1    (t ≥ t0).

Следовательно, если ||x|| = 1, то ||δx/2|| < δ и

||Φt0(t)x0|| =

2 δ

||Φt0(t)(δx/2)|| <

2 δ

.

Поэтому ||Φ_t0(t)|| < 2/δ, т. е. Φ_t0(t) ограничена.

Если, наоборот, известно, что ||Φ_t0(t)|| ≤ H   (t ≥ t₀),

то ||g_t0^t(x₀)||≤ H||x₀||,

так что для любого ε > 0 в определении устойчивости нулевого решения (ЛОС) можно взять δ = ε/H.

(б) Пусть (ЛС) асимптотически устойчива. Тогда ||x0|| < Δ ⇒ ||Φt0(t)x0|| → 0 при t → +∞.

В частности для орта e_k

||Φt0(t)ek|| = 2||ek|| Δ · || Φt0(t) ( ek· Δ 2||ek|| ) || → 0 при t → +∞

(мы рассматриваем произвольную норму в Rn, поэтому, возможно, ||ek|| ≠ 1). Это означает, что все столбцы матрицы Φt0(t) стремятся к нулю при t → +∞; но тогда и сама матрица стремится к нулю.

Пусть дано, что Φt0(t) → 0 при t → +∞. Тогда для любого x0 ∈ Rn

gt0t(x0)= Φt0(t)x0 → 0 при t → +∞,

т. е. (ЛС) асимптотически устойчива в целом.

Наконец, из асимптотической устойчивости в целом следует асимптотическая устойчивость.

(в) Если (ЛС) экспоненциально устойчива, то существуют Δ₁ > 0, M > 0 и γ > 0 такие, что

||x0|| < Δ1 ⇒ ||Φt0(t)x0|| ≤ Me–γ(t–t0)||x0||    (t ≥ t0).

Поэтому для любого x, удовлетворяющего условию ||x|| = 1, будем иметь:

||Φt0(t)x|| = 2 Δ1 || Φt0(t) ( x· Δ1 2 ) || ≤

≤ 2 Δ1 Me–γ(t–t0) || x0· Δ1 2 || = Me–γ(t–t0)||x||.

Следовательно,

||Φt0(t)|| ≤ Me–γ(t–t0)    (t ≥ t0).

Наоборот, если выполнено последнее неравенство, то для любого x₀

||Φt0(t)x0|| ≤ ||Φt0(t)||·||x0|| ≤ Me–γ(t–t0)||x0||   (t ≥ t0),

т. е. (ЛС) экспоненциально устойчива в целом.

Остается заметить, что экспоненциальная устойчивость в целом влечет экспоненциальную устойчивость.