СЕМИНАРЫ 5 И 6
Система двух автономных обыкновенных линейных
дифференциальных уравнений. Фазовая плоскость. Изоклины. Построение фазовых портретов. Кинетические
кривые. Знакомство с программой TRAX.
Фазовой плоскостью называется плоскость с осями
координат, на которых отложены значения переменных x
и y, каждая точка плоскости соответствует определенному состоянию системы. Совокупность точек на фазовой
плоскости, положение которых соответствует состояниям системы в процессе изменения во времени переменных x(t ), y (t ) согласно заданным уравнениям исследуемой системы, называется фазовой траекторией. Совокупность фазовых траекторий при различных начальных значениях переменных дает портрет системы. Построение фазового портрета позволяет сделать выводы о характере изменений переменных x и y без знания аналитических решений исходной системы уравнений.
Для построения фазового портрета на фазовую плоскость наносят изоклины (метод изоклин). Изоклина — линия на плоскости, в каждой точке которой, касательные к фазовым траекториям исследуемой системы уравнений имеют один угол наклона.
Учебное пособие «Математические модели в биологии»
Пусть уравнения имеют вид:
dx dt = P( x, y ), (5.1) dy = Q( x, y ).
dt Тогда уравнение изоклины запишется как:
dy Q( x, y ) = = A = const. (5.2) dx P( x, y ) В уравнении (5.2) константа A есть тангенс угла наклона A = tg касательной к фазовой траектории. Через главные изоклины (нуль-изоклины) фазовые траектории проходят под углом = 0 (изоклина горизонтальных касательных) и = 90 (изоклина вертикальных касательных). Для изоклины горизонтальных касательных уравнение (5.2.) принимает вид:
dy Q( x, y ) = = tg 0 = 0 или Q( x, y ) = 0 ;
dx P( x, y ) для изоклины вертикальных касательных:
dy Q( x, y ) = = tg 90 = или P( x, y ) = 0.
dx P( x, y ) Все изоклины пересекаются в особой точке ( x, y ), для которой Q( x, y ) = P( x, y ) = 0.
Рассмотрим систему линейных дифференциальных уравнений:
dx dt = ax + by = P ( x, y ), (5.3) dy = cx + dy = Q( x, y ).
dt Семинар 5-6. Система двух линейных ОДУ. Фазовый и кинетический портрет Построение фазового портрета начинаем с построения главных изоклин. Для системы двух линейных уравнений — это всегда прямые, проходящие через начало координат. Уравнение изоклины горизонтальных касаc тельных: Q ( x, y ) = 0 y = x. Уравнение изоклины верd a тикальных касательных: P ( x, y ) = 0 y = x. Для дальb нейшего построения фазового портрета полезно построить изоклину касательных, проходящих под углом = ±45.
Для нахождения соответствующего уравнения изоклины cx + dy dy = tg = ±1. Можнеобходимо решить уравнение ax + by dx но находить и изоклины касательных других углов, пользуясь приблизительными значениями тангенсов углов. В построении фазового портрета также может помочь ответ на вопрос, под каким углом фазовые траектории должны пересекать координатные оси. Для этого в dy Q( x, y ) cx + dy = tg подставляем = = уравнение изоклины dx P ( x, y ) ax + by соответствующие равенства x = 0 (для определения угла пересечения с осью OY) и y = 0 (для определения угла пересечения с осью OХ).
Рассмотрим примеры построения фазового и кинетического портрета поведения траекторий системы вблизи особой точки. Построить кинетический портрет системы — означает построить графики зависимости величин переменных x, y от времени. По фазовому портрету можно построить кинетический, и наоборот. Одной фазовой траектории соответствует одна пара кинетических кривых x(t ), y (t ).
Учебное пособие «Математические модели в биологии»
ПРИМЕР 5.1. Рассмотрим систему линейных уравнений:
Рис. 5.1. Фазовый и кинетический портреты системы, пример 5.1.
Координаты особой точки — (0,0). Коэффициенты линейных уравнений равны: a = 1, b = 4, c = 1, d = 1. Определим тип стационарного состояния (см. Семинар 4, стр. 51):
Таким образом, характеристические корни являются мнимыми: 1,2 = ±i 3, следовательно, особая точка рассматриваемой линейной системы имеет тип «центр».
Семинар 5-6. Система двух линейных ОДУ. Фазовый и кинетический портрет Уравнение изоклины горизонтальных касательных:
y = x, уравнение изоклины вертикальных касательных:
y = x. Найдем уравнение изоклины, которую траектории системы пересекают под углом в 45 к оси абсцисс:
После построения фазового портрета необходимо определить направление движения по найденным траекториям.
Это можно сделать следующим образом. Возьмем произвольную точку на любой траектории. Например, на изоклине горизонтальных касательных (1,1). Подставим координаты этой точки в систему уравнений. Получим выражения для скоростей изменения переменных x, y в этой точке:
Получившиеся значения нам показывают, что скорость изменения переменной x — отрицательная, то есть ее значение должно уменьшаться, а переменная y — не изменяется. Отмечаем полученное направление стрелкой (против часовой стрелки). Таким образом, в рассматриваемом примере движение по фазовым траекториям направлено против часовой стрелки.
Подставляя в систему уравнений координаты разных точек, можно получить «карту» векторов, задающих направление и скорость движения изображающей точки по фазовой кривой, так называемое векторное поле.
Отметим, что на изоклине горизонтальных касательных переменная y достигает своего максимального или минимального значения на данной траектории. Наоборот, на изоклине вертикальных касательных, своего максимального по модулю значения для выбранной траектории достигает переменная x.
Учебное пособие «Математические модели в биологии»
ПРИМЕР 5.1. Окончание.
Перейдем к построению кинетического портрета системы. Выберем на фазовом портрете произвольную точку на произвольной фазовой траектории. Это начальная точка, соответствующая моменту времени t = 0. В зависимости от направления движения в рассматриваемой системе значения переменных x, y либо уменьшаются, либо увеличиваются.
Пусть координаты начальной точки — (1,1). Стартуя из этой точки, мы должны двигаться против часовой стрелки, координаты x и y уменьшаются. Координата x проходит через 0, значение y при этом положительно. Далее координаты x и y продолжают уменьшаться, координата y проходит через (значение x при этом отрицательно). Величина x достигает минимального значения на изоклине вертикальных касательных, затем начинает увеличиваться. Величина y своего минимального значения достигает на изоклине горизонтальных касательных (значение x в этот момент времени отрицательно). Далее и величина x, и величина x увеличиваются, возвращаясь к начальным значениям.
Семинар 5-6. Система двух линейных ОДУ. Фазовый и кинетический портрет ПРИМЕР 5.2. Рассмотрим систему линейных уравнений Рис. 5.2. Фазовый и кинетический портреты системы, пример 5.2.
Координаты особой точки — (0,0). Тип особой точки — «неустойчивый узел». Уравнение изоклины горизонтальных касательных: y = x, уравнение изоклины вертикальных касательных: x = 0. Найдем уравнение изоклины, которую траектории системы пересекают под углом в 45 к оси абсцисс:
Направление движения по траекториям можно определять аналитически по знаку (a + d ) (устойчивая или неустойчивая точка) или методом векторного поля.
Кинетический портрет строится с помощью рассуждений, аналогичных предыдущему случаю. В качестве начальной точки выбрана точка A с координатами ( Ax, Ay ).
Учебное пособие «Математические модели в биологии»
ПРИМЕР. 5.3. Рассмотрим систему линейных уравнений Рис. 5.3. Фазовый и кинетический портреты системы, пример 5.3.
Координаты особой точки — (0,0). Тип особой точки — «устойчивый фокус». Уравнение изоклины горизонтальных касательных: y = 3x, уравнение изоклины вертикальных касательных: y = x. Найдем уравнение изоклины, которую траектории системы пересекают под углом в 45 к оси абсцисс:
Аналогично случаю с узлом направление движения по траекториям можно определять аналитически по знаку (a + d ) (устойчивая или неустойчивая точка) или методом векторного поля.
Кинетический портрет строится с помощью рассуждений, аналогичных предыдущему случаю. В качестве начальной точки выбрана точка A с координатами ( Ax, Ay ).
Семинар 5-6. Система двух линейных ОДУ. Фазовый и кинетический портрет ПРИМЕР 5.4. Рассмотрим систему линейных уравнений Рис. 5.4. Фазовый и кинетический портреты системы, пример 5.4.
Координаты особой точки — (0,0). Тип особой точки — «седло». Уравнение изоклины горизонтальных касательных:
y = x, уравнение изоклины вертикальных касательных:
Определим, под каким углом фазовые траектории пересекают оси координат.
траектории должны пересекать под углом примерно (рис. 5.4), должны пересекать под углом примерно 75.
Учебное пособие «Математические модели в биологии»
Важной характеристикой фазового портрета особой точки типа «седло» являются две сепаратрисы. Эти прямые всегда направлены вдоль собственных векторов матрицы коэффициентов линейных уравнений систеa b мы:. Напомним, что собственным вектором матриc d цы M, соответствующим собственному числу, называется любой отличный от нуля вектор x, который удовлетворяет уравнению Mx = x. Итак, уравнения прямыхсепаратрис задаются уравнениями где 1,2 — характеристические числа матрицы коэффициентов системы. Одному значению соответствует одна прямая (выражения 1 и 2 задают совпадающие прямые).
Сепаратрисы могут совпадать с главными изоклинами.
Кроме того, в роли сепаратрис могут выступать оси координат: например, если коэффициент b = 0, то из уравнения ( a 1,2 ) x + 0 y = 0 получаем уравнение сепаратрисы x = 0 (ось OY); если характеристическое число совпадает, например, с коэффициентом a =, то получаем уравнение 0 x + b y = 0, из которого следует, что прямая y = (ось OX) является сепаратрисой.
Семинар 5-6. Система двух линейных ОДУ. Фазовый и кинетический портрет ПРИМЕР 5.4. Окончание.
В разбираемом примере характеристические числа равны 1,2 = (1 5) ± (1 5)2 4 (1 (5) 4 4 ) =.
ся, что вторая пара выражений даст те же уравнения сепаx Кинетический портрет строится с помощью рассуждений, аналогичных предыдущему случаю. В качестве начальной точки выбрана точка A с координатами ( Ax, Ay ).
На что необходимо обратить внимание: особая точка типа «седло» — всегда неустойчива. На больших временах фазовые траектории уходят в бесконечность или, при соответствующем расположении сепаратрис, асимптотически стремятся к одной из осей координат.
На практических занятиях для компьютерного построения фазовых и кинетических портретов мы используем программу TRAX. В качестве рабочей программы можно также использовать MathCad, MatLab (коммерческие программы), XPP/WinPP, Octave, Scilab (свободно распространяемые программы) и другие программные средства.
Учебное пособие «Математические модели в биологии»
ПАМЯТКА ДЛЯ РАБОТЫ В ПРОГРАММЕ TRAX
Создание модели 1) Выбрать пункт меню «Differential equations» (стрелкой «вниз»), Enter.2) Нажать клавишу Insert и ввести имя новой модели, Enter.
3) Ввести уравнения. Параметры модели обозначаются только символами p0, p1…p9.
Работа с моделями 1) Перемещение по окну параметров осуществляется с помощью клавиш Page Up, Page Down.
2) Для того, чтобы задать значение параметра, необходимо установить курсор в нужную строчку, ввести численное значение и нажать клавишу Enter.
3) Изменить начальное положение курсора на фазовой плоскости или в окне кинетической кривой можно с помощью стрелок «вверх», «вниз», «вправо», «влево».
Для увеличения шага после выбора направления нажать клавишу Tab.
4) Для просмотра уравнений нажмите Shift+F1.
Построение графиков F10 строить фазовую или кинетическую кривую в прямом времени F9 строить фазовую или кинетическую кривую в обратном времени F8 построить оси координат F7 очистить окно Esc остановить счет (при отжатой клавише Num Lock) Для переключения между окном «фазовая плоскость» (плоскость Y ( X ) ) и окном «кинетика» ( X (t ), Y (t ) ) в окне параметров переключаться между значениями window 1 и window 2.
Завершение работы F3 сохранение построенного графика F2 выход из программы Семинар 5-6. Система двух линейных ОДУ. Фазовый и кинетический портрет
ЗАДАЧИ К СЕМИНАРАМ 5 И
5.1. Постройте фазовый портрет для каждой из систем задачи 4.1 в окрестности стационарного состояния:а) отметьте стационарную точку на фазовой плоскости;
б) постройте главные изоклины систем, изоклины ±45 ;
в) определите, под каким углом фазовые траектории должны пересекать оси координат фазовой плоскости;
г) по изоклинам постройте эскиз фазового портрета, стрелкой укажите направление движения изображающей точки вдоль интегральных кривых при t.
5.2. Для произвольной начальной точки на фазовом портрете задачи 5.1. постройте кинетический портрет.