4Дифференциальные уравнения в многомерных фазовых пространствах

4.1Многомерные фазовые пространства

До сих пор мы рассматривали дифференциальные уравнения с одномерным фазовым пространством: искомая функция принимала значения во множестве вещественных чисел. На практике нам зачастую нужны более сложные дифференциальные уравнения с многомерными фазовыми пространствами. Для простоты обозначений здесь и ниже мы будем рассматривать преимущественно двумерные фазовые пространства. Большинство определений мгновенно обобщаются на случай произвольной размерности.

4.1.1Отступление: что такое фазовое пространство?

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

Пример 1. Зачастую просто использование понятия фазового пространства позволяет сильно продвинуться в понимании задачи. Приведем задачу Н.Н.Константинова из книги Арнольда:

Пусть есть два города — A и B — и между ними проведены две дороги. Из города А в какой-то момент вышли Ромео и Джульетта и направились в сторону города B. При этом они решили идти каждый по своей дороге, держась за веревочку длиной 5 метров — Ромео за один конец, Джульетта за другой. Веревочка может не быть натянутой (то есть расстояние между Ромео и Джульеттой может быть меньше 5 метров, но не может быть больше), они могут притормаживать, ускоряться, пропускать друг друга вперед, идти назад и т.д. — двигаться как угодно. Известно, что двигаясь таким образом им удалось попасть из города A в город B.

В другой день из города B в город A выехал воз, нагруженный товаром, а из города A в город B выехал другой воз, тоже нагруженные товаром, и ехали они по разным тропинкам. Возы круглые (если смотреть сверху), их радиус 3 метра.

Утверждение. Если два человека смогли пройти, держась за верёвочку, то возы не смогут разъехаться.
Доказательство. Для решения мы нарисуем фазовое пространство — в данном случае, фазовую плоскость. По горизонтальной оси отложим расстояние до объекта (будь то человек или воз) до города А по одной дороге, а по вертикальной — расстояние до другого объекта по другой дороге. Таким образом, расположение обоих объектов будет задано одной точкой на плоскости. Движению обоих объектов будет соответствовать некоторая кривая на фазовом пространстве — фазовая кривая (или траектория).
import matplotlib.pyplot as plt
import numpy as np
import qqmbr.odebook as ob
# see https://github.com/ischurov/qqmbr/blob/master/qqmbr/odebook.py

plt.axis([-0.1,1.1,-0.1,1.1])
ob.center_spines()
plt.plot([1, 1, 0, 0, 1],[0, 1, 1, 0, 0], linewidth=2)
plt.text(0.02, 1.02, "B",fontsize=20)
plt.text(1.02, 0.02, "B",fontsize=20)
plt.text(0.02, 0.02, "A",fontsize=20)
ob.mplot(np.linspace(0, 1), 
         lambda x: 1 - ((1 - x)**2 + np.sin(3 * np.pi * (1 - x)) * 0.3 * (x - x**2)))
ob.mplot(np.linspace(0, 1), lambda x: (1 - x**2) - np.cos(3 * np.pi * x) * 0.3 * (x - x**2))
Рис. 4.1: Фазовая плоскость
В начальный момент времени Ромео и Джульетта находились в городе A, что соответствует точке на картинке. В конечный момент времени они оказались в городе B, что соответствует точке на картинке. Таким образом, если смотреть на траектории, они «прошли» из левого нижнего угла в правый верхний. Наоборот, возы в фазовом пространстве двигались из левого верхнего угла в правый нижний. На картинке рыжим цветом изображено движения возов, зеленым — движение людей. Очевидно, что траектория может быть не только графиком функции, но и произвольной непрерывной кривой. Теперь заметим, что две кривые, соединяющие противоположные углы квадрата, в какой-то момент обязательно пересекутся. (Строгое доказательство этого утверждения нетривиально — оно требует сложной в доказательстве леммы Жордана — но интуитивно оно не вызывает сомнений.) Это означает, что как минимум в одном месте фазового пространства одновременно будет выполнены условия, что расстояние между возами не более пяти (люди держатся за веревочку), но и не менее шести (сумма радиусов возов). Противоречие.

Пример 2. В модели Солоу рассматривалось одномерное фазовое пространство — единственной переменной была капиталовооруженность , потребление считалось постоянным. В более сложной модели Рамсея учитывается, что потребление может меняться со временем, и оно вступает в игру как одна из неизвестных функций. Фазовое пространство становится двумерными, фазовыми переменными являются капиталовооруженность и потребление .

Пример 3. Когда мы обсуждали мальтузианский рост популяции, наше пространство было одномерным: нас интересовал только размер популяции в данный момент времени. Если бы мы рассмотрели более сложную модель, включающую, например, взаимодействующие популяции двух разных видов, нам потребовалось бы два числа для описания состояния системы — количество особей одного вида и другого вида. Рассмотрим, скажем, простейшую модель взаимодействия двух видов, один из которых является хищником, а другой — жертвой — например, взаимодействие кроликов и лис, или щук и карасей.

Пусть — число лис, — число кроликов. Если ни одной лисы нет, скорость роста числа кроликов пропорциональна числу самих кроликов (как в мальтузианской модели — какая-то часть популяции воспроизводится за единицу времени). Наоборот, если нет кроликов, то лисы вымирают от голода — за единицу времени какая-то доля популяции погибает. Каждая встреча лисы с кроликом (которые происходят с частотой, пропорциональной произведению ) вносит положительный вклад в динамику лис и отрицательный — в динамику кроликов. Запишем дифференциальное уравнение (систему из двух уравнений):

где , , , — некоторые положительные параметры. Фазовым пространством здесь является первая координатная четверть . Эта система называется моделью Лотки—Вольтерры. Мы к ней ещё вернёмся в будущем.

4.2Дифференциальные уравнения в многомерных пространствах

Нам потребуется напомнить некоторые определения из математического анализа.

Определение 2. Областью в пространстве мы будем называть связное подмножество в . Часто под областью подразумевают открытое подмножество, но мы не будем накладывать этого требования. Например, отрезок интервал, полуинтервал, луч — все они являются областями в . Первая четверть (открытая или замкнутая — не важно) — область в .Подробнее

Определение 3. Вектор-функцией называется отображение из области числовой прямой в многомерное векторное пространство. Вектор-функцию можно представлять себе как упорядоченный набор числовых функций (то есть отображений числовой прямой в себя):
где для всех .

Определение 4. Производной вектор-функции , называется вектор-функция, составленная из производных каждой компоненты функции :

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

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

Его решением является вектор-функция , такая, что
для всех из области определения .

Понятия «система дифференциальных уравнений» и «дифференциальное уравнение с многомерным фазовым пространством» эквивалентны. Действительно, рассмотрим, например, систему из двух дифференциальных уравнений:

Систему (4.3) можно представить как одно дифференциальное уравнение вида (4.1) на вектор-функцию , где .

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

4.2.1Теорема существования и единственности решения задачи Коши

Мы уже формулировали кратко теорему существования и единственности для многомерных дифференциальных уравнений, но сейчас не грех и повторить.

Теорема 1. Рассмотрим уравнение (4.1), в котором правая часть непрерывно дифференцируема (-гладкая, то есть дифференцируемая и все частные производные непрерывны) в окрестности точки расширенного фазового пространства.

Тогда найдётся такая окрестность , что на существует решение уравнения (4.1), удовлетворяющее условию , и при этом любое другое решение уравнения (4.1), удовлетворяющее этому же условию, совпадает с на некоторой окрестности точки .

4.3Автономные системы дифференциальных уравнений

4.3.1Геометрия вектор-функций

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

Определение 5. Траекторией вектор-функции называется кривая в , являющаяся образом для отображения , с отмеченным на ней направлением, соответствующим увеличению аргумента. Иными словами, траектория — это кривая
на которой нарисована стрелочка, соответствующая направлению увеличения параметра .

Пример 4. Рассмотрим вектор-функцию , . Её траекторией будет единичная окружность, ориентированная против часовой стрелки.
import matplotlib.pyplot as plt
import numpy as np
import qqmbr.odebook as ob
# see https://github.com/ischurov/qqmbr/blob/master/qqmbr/odebook.py

plt.figure(figsize=(3,3))

plt.axis([-1.1, 1.1, -1.1, 1.1])
ob.center_spines()
T = np.linspace(0, 2*np.pi, 200)
plt.plot(np.cos(T), np.sin(T))
ax = plt.gca()
ax.arrow(np.sqrt(1/2),np.sqrt(1/2),
         -0.00001, 0.00001, head_width=0.04, 
         head_length=0.07, fc='b', ec='b')

Вопрос 1. Что будет траекторией для вектор-функции ?
  Такая же окружность.

Не совсем. Что насчёт ориентации?

  Окружность радиуса 2.

Нет. Подставьте любое значение и посмотрите, что получится.

  Единичная окружность, ориентированная по часовой стрелке.

Верно!

Если кривая задана как траектория некоторой вектор-функции, говорят также, что эта кривая задана параметрически: параметр — это аргумент вектор-функции. Как следует из вопроса 1, одна и та же кривая может быть задана разными вектор-функциями.

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

Пример 5. Для вектор-функции производной является вектор-функция . Нетрудно видеть, что вектор перпендикулярен радиус-вектору точки . На картинке это будет выглядеть так.
import matplotlib.pyplot as plt
import numpy as np
import qqmbr.odebook as ob
# see https://github.com/ischurov/qqmbr/blob/master/qqmbr/odebook.py

plt.figure(figsize=(4,4))

plt.axis([-2, 2, -2, 2])
ob.center_spines()
T = np.linspace(0, 2*np.pi, 200)
plt.plot(np.cos(T), np.sin(T))
ax = plt.gca()

T = np.linspace(0, 2*np.pi, 4*4+1)
plt.quiver(np.cos(T), np.sin(T), -np.sin(T), np.cos(T), 
           scale=0.9, scale_units='x', color='Teal', headlength=7)

Утверждение 1. Вектор производной вектор-функции в точке касается траектории этой функции в точке .

Упражнение. Доказать утверждение 1.

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

Формальное определение такое:

Определение 6. Графиком вектор-функции называется подмножество декартова произведения , устроенное следующим образом:

Пример 6. Графиком вектор-функции будет спираль.
import plotly
import plotly.graph_objs as go
from plotly.offline import iplot as plot
from plotly.offline import init_notebook_mode

init_notebook_mode()

import numpy as np

t = np.linspace(0, 2*np.pi,50)

trace1 = go.Scatter3d(
    x = t,
    y = np.cos(t),
    z = np.sin(t),
    mode = 'lines+markers',
    marker = {'size': 3},
    hoverinfo = 'none',
    name = 'График',
    line = go.Line(width = 4),
    projection = {'x': {'show': True, 'opacity': 1}},
)
trace2 = go.Scatter3d(
    x = np.zeros_like(t),
    y = np.cos(t),
    z = np.sin(t),
    mode = 'lines',
    hoverinfo = 'none',
    name = 'Траектория',
    line = go.Line(width = 4)
)
trace3 = go.Scatter3d(
    x = t,
    y = np.zeros_like(t),
    z = np.zeros_like(t),
    mode = 'lines',
    hoverinfo = 'none',
    name = 'x=0, y=0'
)
layout = go.Layout(
                width = 800, 
                height = 600, 
                scene = go.Scene(
                    xaxis = {'title': 't'},
                    yaxis = {'title': 'x'},
                    zaxis = {'title': 'y'},
                    camera = {'eye': {'x': 1, 'y': -2, 'z': 1}}
                )
)
plot(dict(data=[trace1, trace3], layout=layout))
Рис. 4.4: Траектория и график вектор-функции . Картинку можно крутить!

Замечание 1. Траектория вектор-функции является проекцией её графика вдоль оси .

4.3.2Автономные дифференциальные уравнения на плоскости

Если правая часть не зависит от явно, система из двух дифференциальных уравнений становится автономной и принимает вид:
Напомним, что интегральной кривой дифференциального уравнения называется график его решения; интегральная кривая лежит в расширенном фазовом пространстве. Для автономных уравнений также оказываются содержательными геометрические объекты, определённые в фазовом пространстве.

Определение 7. Траекторией (или фазовой кривой) системы (4.4) называется траектория вектор-функции, являющейся решением этой системы.

Замечание 2. Фазовая кривая — это проекция интегральной кривой на фазовое пространство вдоль оси времени (ср. с замечанием 1).

Утверждение 2. Через любую точку фазовой плоскости уравнения (4.4) проходит фазовая кривая, причём ровно одна.

Доказательство. Возьмём произвольный момент времени . По теореме существования и единственности существует и единственно решение уравнения (4.4), определённое в некоторой окрестности точки , и такое, что , . Его траектория даёт фазовую кривую, проходящую через точку .

Остаётся доказать, что эта траектория единственна. Для конкретного единственность следует из единственности решения; остаётся доказать, что от выбора фазовая кривая не зависит. Действительно, уравнение (4.4) автономно, а, следовательно сдвиг по времени начального условия приведёт к такому же сдвигу по времени всего решения. Иными словами, для другого начального момента времени решением будет функция . Траектории у функций и совпадают. Следовательно, фазовая кривая единственна.

import plotly
import plotly.graph_objs as go
from plotly.offline import iplot as plot
from plotly.offline import init_notebook_mode

init_notebook_mode()

import numpy as np

t = np.linspace(0, 1.5, 20)
t0 = 0.5
data = []
def f(t):
    return (np.cos(2*t), np.sin(2*t))

data.append(go.Scatter3d(
    x = t,
    y = f(t)[0],
    z = f(t)[1],
    mode = 'lines+markers',
    hoverinfo = 'none',
    name = 'График',

    line = go.Line(width = 4),
    marker = {'size': 3},
    projection = {'x': {'show': True, 'opacity': 0.6}}
))
data.append(go.Scatter3d(
    x = [t0],
    y = [f(t0)[0]],
    z = [f(t0)[1]],
    hoverinfo = 'none',
    showlegend = False,
    projection = {'x': {'show': True}}
))
deltat = 0.5
data.append(go.Scatter3d(
    x = t + deltat,
    y = f(t)[0],
    z = f(t)[1],
    mode = 'lines+markers',
    hoverinfo = 'none',
    name = 'Траектория',
    line = go.Line(width = 4),
    marker = {'size': 3},
    projection = {'x': {'show': True, 'opacity': 0.7}}
))
data.append(go.Scatter3d(
    x = [t0 + deltat],
    y = [f(t0)[0]],
    z = [f(t0)[1]],
    hoverinfo = 'none',
    showlegend = False,
    projection = {'x': {'show': True}}
))
layout = go.Layout(
                width = 800, 
                height = 600, 
                scene = go.Scene(
                    xaxis = {'title': 't'},
                    yaxis = {'title': 'x', 'range': [-1.2, 1.2]},
                    zaxis = {'title': 'y', 'range': [-1.2, 1.2]},
                    camera = {'eye': {'x': 1, 'y': -2, 'z': 1}}
                ),
                showlegend = False,
)
plot(dict(data=data, layout=layout))
Рис. 4.5: Интегральные кривые, отличающиеся сдвигом по времени, проецируются в одну и ту же фазовую кривую

Определение 8. Если нарисовать на фазовом пространстве все возможные траектории, получится фазовый портрет.

С интегральными кривыми связаны поля направлений; для фазовых кривых есть похожий объект: векторное поле.

Определение 9. Векторным полем, заданным системой (4.4), называется следующий объект: в каждой точке фазового пространства нарисован вектор, выходящий из точки и имеющий координаты .

Замечание. Каждый вектор векторного поля показывает, с какой скоростью движется точка вдоль фазовой кривой, проходя определённое место в фазовом пространстве. Фазовые кривые системы касаются векторов соответствующего векторного поля в каждой своей точке. Это мгновенно следует из утверждения 1.

Рассмотрим несколько примеров.

Пример 7. Нулевое поле:
В каждой точке отложен нулевой вектор.
import matplotlib.pyplot as plt
import numpy as np
import qqmbr.odebook as ob
# see https://github.com/ischurov/qqmbr/blob/master/qqmbr/odebook.py

plt.figure(figsize=(6,6))
plt.axis([-4,4,-4,4])
X = np.arange(-4,4,0.7)
x, y = np.meshgrid(X, X)
plt.plot(x, y, 'o', color='Teal');
Рис. 4.6: Нулевое векторное поле: кружочки символизируют нулевые векторы

Вопрос 2. Как выглядят фазовые кривые этого уравнения?
  Произвольные кривые

Это неверно: у системы (4.5) и начального условия всегда существует тождественное решение , а по теореме существования и единственности других решений нет.

  Вся плоскость

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

  Каждая точка является отдельной фазовой кривой

Верно! У системы (4.5) и начального условия всегда существует тождественное решение , а по теореме существования и единственности других решений нет. Траектория постоянной вектор-функции — одна точка.

Пример 8. Постоянное поле:
В каждой точке отложен один и тот же вектор с координатами . Фазовыми кривыми являются параллельные прямые с угловым коэффициентом .
import matplotlib.pyplot as plt
import numpy as np
import qqmbr.odebook as ob
# see https://github.com/ischurov/qqmbr/blob/master/qqmbr/odebook.py

plt.figure(figsize=(4,4))
fs = lambda x, y: (1,2)
ob.vectorfield(np.linspace(-4, 4, 6),
               np.linspace(-4, 4, 6),
               fs, color='Teal', 
               headlength=7, scale=0.9,
               scale_units='x')
Рис. 4.7: Постоянное векторное поле

Пример 9. Эйлерово поле:

import matplotlib.pyplot as plt
import numpy as np
import qqmbr.odebook as ob
# see https://github.com/ischurov/qqmbr/blob/master/qqmbr/odebook.py

plt.figure(figsize=(4,4))
plt.xlim(-1, 1)
plt.ylim(-1, 1)
theta = np.linspace(0, 2*np.pi, 4*5+1)
r = np.linspace(0, 1, 5)**3

r, theta = np.meshgrid(r, theta)

x = np.cos(theta)*r
y = np.sin(theta)*r

plt.quiver(x, y, x, y, color='Teal', headlength=7, scale=1,
            scale_units='x')
Рис. 4.8: Эйлерово векторное поле

Вопрос 3. Каковы фазовые кривые этого уравнения?
  Прямые, проходящие через начало координат.

Неверно. Фазовая кривая, стартовавшая вне начала координат, не может попасть в начало координат: это противоречило бы теореме существования и единственности.

  Открытые лучи, выходящие из начала координат.

Не совсем верно. Есть ещё одна фазовая кривая.

  Открытые лучи, выходящие из начала координат, а также само начало координат.

Так и есть!

  Открытые лучи, выходящие из начала координат (кроме вертикального), и само начало координат

Нет, в данном случае вертикальный луч является полноправной фазовой кривой: ничто не запрещает его рассматривать.

4.3.3Особые и неособые точки векторного поля

Точка фазового пространства называется особой, если векторное поле в этой точке равно нулю (нулевому вектору). Иными словами, правая часть дифференциального уравнения в особой точке обнуляется. Все остальные точки называются неособыми. Особые точки также называются положениями равновесия (equilibrium point или steady state). Пусть — особая точка для некоторого дифференциального уравнения. Тогда вектор-функция , тождественно равная этой точке, очевидно, является решением соответствующего уравнения (ведь её производная равна нулевому вектору). Если уравнение удовлетворяет условию теоремы существования и единственности, эта функция является единственным решением уравнения с начальным условием . Иными словами, если решение начинается в особой точке, то оно вечно живёт в этой особой точке (и, более того, вечно жило в этой точке в прошлом). Никакое другое решение попасть в особую точку не может (но может долго-долго к ней приближаться). Особые точки являются отдельными траекториями и обычно изображаются на фазовом портрете маленькими кружочками.

4.4Связь между автономными и неавтономными уравнениями

Между автономными и неавтономными дифференциальными уравнениями есть важная связь, которую неформально можно сформулировать так: мир неавтономных систем ровно на одно измерение богаче мира автономных систем; чтобы перейти от неавтономной системы к автономной, нужно увеличить размерность фазового пространства на один, а для обратного перехода — уменьшить на один.

Сейчас мы обсудим одну из версий этого принципа.

Теорема 2. Рассмотрим систему
Для любой точки , такой, что , фазовая кривая системы (4.6), проходящая через , совпадает с интегральной кривой для уравнения
Доказательство. Построим векторное поле системы (4.6), а также поле направлений уравнения (4.7). Рассмотрим произвольную точку , близкую к точке . В ней поле направлений задается вектором . Угловой коэффициент прямой, на которой лежит этот вектор, равен . Угловой коэффициент прямой из поля направлений уравнения (4.7), проведенной в точке , будет равен тому же числу: . По теореме о существовании и единственности, кривая, касающаяся в каждой своей точке соответствующего направления, единственна. Значит, фазовая кривая автономного уравнения и интегральная кривая неавтономного совпадают.

Этот принцип работает во всех точках, когда . В точках, когда это не выполнено можно поменять ролями и : выражать не через , а через . Наконец, когда оба выражения обнуляются: — то эта точка является особой, и в ней не определено поле направлений.

Пример 10. Рассмотрим систему, задающую эйлерово векторное поле:
Уравнения, входящие в эту систему, не зависят друг от друга (в уравнении на фигурирует только , в уравнении на — только ). Решения этих уравнений нам уже известны (см. параграф 1.2.6):
где — начальные условия.

Уравнения на фазовые кривые можно получить, выразив через или через . Например, из первого уравнения следует, что и подставляя это во второе уравнение, имеем:

Иными словами, фазовые кривые лежат на прямых, проходящих через начало координат. Из выражения для решений (4.8) видно, что фазовые кривые являются открытыми лучами: решения стремятся к точке при , но никогда не достигают этой точки.

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

Это уравнение с разделяющимися переменными и оно легко решается: Поскольку и никогда не меняют знак, можно снять модули и получить уравнение , задающее то же семейство прямых, что и раньше.

При этом, однако, мы потеряли некоторую информацию: в частности, направление движения по фазовым кривым, а также тот факт, что на самом деле фазовые кривые являются лучами, а не прямыми. Мы также потеряли вертикальные фазовые кривые и фазовую кривую, соответствующую тождественно нулевому решению (она изображается одной точкой). Всё это связано с условием в теореме 2.

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

  1. Найти решения, то есть зависимость и .
  2. Найти фазовые кривые, то есть зависимость от или от .
Для решения второй из этих задач полезно использовать теорему 2, но при этом мы теряем некоторую информацию об исходной системе.

4.4.1Репараметризация фазовых кривых

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

Рассмотрим систему

где . Понятно, что уравнение (4.7) при этом не изменится, поскольку окажется в числителе и знаменателе и сократится. Из теоремы 2 теперь следует, что не изменятся и фазовые кривые. Однако решения уравнения при этом изменятся.

Можно сказать, что при этом происходит репараметризация фазовых кривых: мы посещаем те же точки фазового пространства, но в другие моменты времени. Сформулируем необходимое определение.

Определение 11. Рассмотрим гладкую вектор-функцию и пусть функция также гладкая, причём её производная всюду положительна и , . Тогда сложная функция (то есть ), , задаёт ту же траекторию, что и исходная вектор-функция . Переход от функции к функции называется репараметризацией кривой.

Если представить себе, что на каждой точке кривой написано значение параметра, при котором мы проходим эту точку, то репараметризация — просто изменение «номеров» точек кривой.

Утверждение 3. При репараметризации вектор производной умножается на число, но не меняет направления.

Доказательство. Это следует из теоремы о производной сложной функции:
При этом — это вектор, — число (причём мы потребовали, чтобы оно было положительным).

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

Вопрос 4. Что произойдёт с фазовыми кривыми при умножении векторного поля на отрицательное число?
  Ничего не произойдёт

Не вполне верно: на фазовой кривой есть стрелочка, показывающая направление движение. Что произойдёт с этими стрелочками?

  Фазовые кривые симметрично отразятся

Нет, это интегральные кривые симметрично отразятся относительно плоскости .

  Кривые останутся на месте, а вот стрелочки все поменяются на противоположные

Верно! Именно так всё и будет.