Дифференциальные уравнения · Фазовый портрет

Фазовый портрет

Фазовый портрет изображает каждое решение системы дифференциальных уравнений кривой на плоскости её переменных. Для линейной системы x′ = Ax след и определитель A задают собственные значения. Если они не совпадают, по ним определяется тип картины.

Фазовые портреты в этой визуализации

Фазовый портрет изображает решения системы дифференциальных уравнений кривыми на плоскости её переменных: через каждую начальную точку проходит одна кривая. Здесь показана линейная система x′=Axx' = Ax с матрицей AA размера 2×22 \times 2. Перемещение точки (τ,Δ)(\tau, \Delta) на карте следа и определителя меняет картину между седлом, узлом, фокусом и центром. На соседней плоскости собственных значений видна пара, которая определяет каждую картину.

Система, её след и определитель
(xy)′=(abcd)(xy),τ=a+d,Δ=ad−bc\begin{pmatrix} x \\ y \end{pmatrix}' = \begin{pmatrix} a & b \\ c & d \end{pmatrix} \begin{pmatrix} x \\ y \end{pmatrix}, \qquad \tau = a + d, \quad \Delta = ad - bc

Каждая кривая следует полю направлений. Все стрелки имеют одинаковую длину и показывают только направление. О скорости говорят метки через равные промежутки времени: там, где движение медленное, они расположены теснее. Некоторые решения идут по прямой через начало координат. Если Av=λvAv = \lambda v, то x(t)=eλtvx(t) = e^{\lambda t} v — решение: его производная λeλtv\lambda e^{\lambda t} v в точности равна результату применения к нему AA. Таким образом, эти прямые — собственные направления AA, а λ\lambda задаёт темп изменения вдоль каждого: при λ>0\lambda > 0 движение направлено наружу, при λ<0\lambda < 0 — внутрь, при λ=0\lambda = 0 движения нет.

Собственные значения удовлетворяют λ2−τλ+Δ=0\lambda^2 - \tau\lambda + \Delta = 0, поэтому τ\tau — их сумма, Δ\Delta — произведение, а τ2−4Δ\tau^2 - 4\Delta — дискриминант. Каждая граница на карте соответствует событию на плоскости собственных значений. Ниже Δ=0\Delta = 0 собственные значения имеют разные знаки, и возникает седло. На этой линии одно из них равно 0, поэтому целая прямая состоит из неподвижных точек. На параболе τ2=4Δ\tau^2 = 4\Delta собственные значения совпадают, а выше неё образуют комплексную пару. При Δ>0\Delta > 0 знак τ\tau отделяет затухание от роста, а на τ=0\tau = 0 пара лежит на мнимой оси. Итак, τ\tau и Δ\Delta определяют собственные значения и тип картины, кроме случая кратного собственного значения. Тогда тип зависит ещё и от того, является ли AA скалярным кратным единичной матрицы.

Для комплексных собственных значений α±βi\alpha \pm \beta i при β>0\beta > 0 любое ненулевое решение обходит начало координат за время 2π/β2\pi/\beta и за это время масштабируется с множителем e2πα/βe^{2\pi\alpha/\beta}. Поэтому метки с интервалом в один оборот лежат на одном луче из начала координат и связаны одним и тем же множителем масштаба. При α=0\alpha = 0 этот множитель равен 1, и любая ненулевая орбита замыкается в эллипс: это центр. При Δ=1\Delta = 1 и A=(01−1τ)A = \begin{pmatrix} 0 & 1 \\ -1 & \tau \end{pmatrix} система представляет собой затухающий осциллятор x′′+2ζx′+x=0x'' + 2\zeta x' + x = 0 с ζ=−τ/2\zeta = -\tau/2.

При пересечении параболы две прямые сходятся, как лезвия ножниц, в одну, а затем исчезают, когда картина становится спиральной. Вне вершины типичная матрица на параболе имеет только одну такую прямую. Все остальные решения входят в начало координат или выходят из него по касательной к ней. В точке (τ,Δ)=(−2,1)(\tau, \Delta) = (-2, 1) матрица −I-I имеет те же собственные значения, однако прямолинейные решения есть во всех направлениях. Замена базиса A↦PAP−1A \mapsto PAP^{-1} меняет форму картины, сохраняя τ\tau, Δ\Delta, собственные значения и тип. Скалярное кратное единичной матрицы остаётся неизменным.

На что обратить внимание

  • Две прямые узла обычно не перпендикулярны. Решение вне этих прямых входит в начало координат или выходит из него вдоль более медленной, а вдали поворачивает к более быстрой.
  • Синяя кривая — полное решение, проходящее через начальную точку. Жёлтый путь — его часть за время с момента добавления этой точки. Поэтому после изменения AA он показывает новое решение, а не историю старого.

Управление

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

Смотрите также

Дополнительные материалы: Википедия: Phase portrait (Английский).