Правила определения типа точки покоя

Можно определить тип точки покоя и характер ее устойчивости, не находя собственных значений матрицы системы (1), а зная только ее след tr A и определитель det A.

Правила определения типа точки покоя - student2.ru

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

Правила определения типа точки покоя - student2.ru

Алгоритм построения фазового портрета ЛДС

1.Определить положения равновесия, решив систему уравнений:

Правила определения типа точки покоя - student2.ru

2. Найти собственные значения матрицы системы, решив характеристическое уравнение:

Правила определения типа точки покоя - student2.ru

3. Определить тип точки покоя и сделать вывод об устойчивости.

4. Найти уравнения главных изоклин горизонтальной и вертикальной, и построить их на фазовой плоскости.

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

6. Нарисовать фазовые траектории.

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

Главные изоклины

Правила определения типа точки покоя - student2.ru

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

Фазовые траектории

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

Уравнения таких прямых можно искать в виде y = kx. Подставляя y = k x в уравнение:

Правила определения типа точки покоя - student2.ru

для определения k получим:

  Правила определения типа точки покоя - student2.ru (4)

(Уравнения прямых, содержащих фазовые траектории, можно искать и в виде x = ky . Тогда для нахождения коэффициентов следует решить уравнение

Правила определения типа точки покоя - student2.ru

Дадим описание фазовых траекторий в зависимости от количества и кратности корней уравнения (4).

Фазовые траектории

Правила определения типа точки покоя - student2.ru

Правила определения типа точки покоя - student2.ru

* Если уравнения прямых ищутся в виде x = ky, тогда это будут прямые x= k1 y и y = 0.

Правила определения типа точки покоя - student2.ru

Если положение равновесия является центром, то фазовые траектории являются эллипсами.

Если положение равновесия является фокусом, то фазовые траектории являются спиралями.

В случае, когда ЛДС имеет прямую точек покоя, то можно найти уравнения всех фазовых траекторий, решив уравнение:

Правила определения типа точки покоя - student2.ru

Его первый интеграл ax + by = C и определяет семейство фазовых прямых.

Направление движения

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

Правда, в случае фокуса требуется установить еще и направление закручивания (раскручивания), спирали по часовой или против часовой стрелки. Это можно сделать, например, так. Определить знак производной y′(t) в точках оси x.

Когда Правила определения типа точки покоя - student2.ru если x > 0, то ордината движущейся точки по фазовой траектории при пересечении «положительного луча оси » возрастает.

Значит «закручивание (раскручивание)» траекторий происходит против часовой стрелки.

Когда Правила определения типа точки покоя - student2.ru если x > 0, то «закручивание (раскручивание)» траекторий происходит по часовой стрелки.

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

Правила определения типа точки покоя - student2.ru

Следовательно, если положение равновесия седло, то достаточно установить направление движения по какой-нибудь траектории. И далее можно однозначно установить направление движения по всем остальным траекториям.

Направление движения (седло)

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

1 способ

Определить, какая из двух сепаратрис соответствует отрицательному собственному значению. Движение по ней происходит к точке покоя.

2 способ

Определить, как изменяется абсцисса движущейся точки по любой из сепаратрис. Например, для y = k1x имеем:

Правила определения типа точки покоя - student2.ru

Если x(t) → 0 при t → +∞, то движение по сепаратрисе y = k1x происходит к точке покоя.

Если x(t) → ±∞ при t→+∞, то движение происходит от точки покоя.

3 способ

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

Когда Правила определения типа точки покоя - student2.ru если x > 0,то ордината точки возрастает и, значит, движение по фазовым траекториям, пересекающим положительную часть оси x, происходит снизу вверх. Если же ордината убывает, то движение будет происходить сверху вниз.

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

4 способ

Правила определения типа точки покоя - student2.ru

Правила определения типа точки покоя - student2.ru

Пример 1.

Правила определения типа точки покоя - student2.ru

1. Система имеет единственное нулевое положение равновесия, так как

det A = 6 ≠ 0.

2. Построив соответствующее характеристическое уравнение l2 – 6 = 0, найдем его корни l1,2 = ± Ö6. Корни вещественные и разного знака. Следовательно, положение равновесия седло.

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

>> a=[-2 -2;-1 2];

>> eig(a)

ans =

-2.4495

2.4495

3. Сепаратрисы седла ищем в виде y = kx.

Правила определения типа точки покоя - student2.ru

4. Вертикальная изоклина: x + y = 0.

Горизонтальная изоклина: x + 2y = 0.

Правила определения типа точки покоя - student2.ru

Правила определения типа точки покоя - student2.ru

Правила определения типа точки покоя - student2.ru

Пример 2.

Правила определения типа точки покоя - student2.ru

1. Система имеет единственное нулевое положение равновесия, так как det A = 10 ≠ 0.

2. Построив соответствующее характеристическое уравнение l2 + 7l + 10 = 0,

найдем его корни l1 = 2, l2 = 5. Следовательно, положение равновесия неустойчивый узел.

>> a=[4 2;1 3];

>> eig(a)

ans =

3. Прямые: y = kx.

Правила определения типа точки покоя - student2.ru

4. Вертикальная изоклина: 2x + y = 0.

Горизонтальная изоклина: x + 3y = 0.

Правила определения типа точки покоя - student2.ru

Найдем собственные вектора матрицы А:

>> A=[4 2;1 3];

>> [U,L]=eig(A)

U =

0.8944 -0.7071

0.4472 0.7071

L =

5 0

0 2

Т.е., собственные значения l1 = 5 и l2 = 2.

Так как l2 = 2 является меньшим по абсолютной величине собственным значением, найдем соответствующий ему собственный вектор:

>> U(:,2)./0.7071

ans =

-1.0000

1.0000

Следовательно, остальные фазовые траектории, образующие параболы, касаются в начале координат прямой y = –x.

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

Правила определения типа точки покоя - student2.ru

Пример 3.

Правила определения типа точки покоя - student2.ru

1. Система имеет единственное нулевое положение равновесия, так как det A = 18 ≠ 0.

2. Построив соответствующее характеристическое уравнение l2 + 3l + 18 = 0,

найдем его дискриминант D = 63. Так как D < 0, то корни уравнения комплексные, причем Re l1,2 = 3/2. Следовательно, положение равновесия устойчивый фокус.

>> a=[-1 -4;4 -2];

>> eig(a)

ans =

-1.5000 + 3.9686i

-1.5000 - 3.9686i

3. Вертикальная изоклина: x + 4y = 0.

Горизонтальная изоклина: 2x – y = 0.

Фазовые траектории являются спиралями, движение по которым происходит к началу координат. Направления «закручивания траекторий» можно определить следующим образом.

Определим знак производной y′(t) в точках оси x. При y = 0 имеем:

Правила определения типа точки покоя - student2.ru если x > 0.

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

Значит, «закручивание» траектории происходит против часовой стрелки.

Правила определения типа точки покоя - student2.ru

Пример 4.

Правила определения типа точки покоя - student2.ru

1. Система имеет единственное нулевое положение равновесия, так как det A = 3 ≠ 0.

2. Построив соответствующее характеристическое уравнение l2+3 = 0, найдем его корни l1,2 = ±iÖ3. Следовательно, положение равновесия центр.

>> a=[1 -4;1 -1];

>> eig(a)

ans =

-0.0000 + 1.7321i

-0.0000 – 1.7321i

3. Вертикальная изоклина: x – 4y = 0.

Горизонтальная изоклина: x – y = 0.

Фазовые траектории системы эллипсы.

Направление движения по ним можно установить, например, так.

Определим знак производной y′(t) в точках оси x. При y = 0 имеем:

Правила определения типа точки покоя - student2.ru если x > 0.

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

Значит, движение по эллипсам происходит против часовой стрелки.

Правила определения типа точки покоя - student2.ru

Пример 5.

Правила определения типа точки покоя - student2.ru

Так как определитель матрицы системы det A = 0, то система имеет бесконечно много положений равновесия. Все они лежат на прямой y = 2x.

Построив соответствующее характеристическое уравнение l2 + 5l = 0,

найдем его корни l1=0, l2= –5.

>> a=[-4 2;2 -1];

>> eig(a)

ans =

-5

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

Построим уравнения остальных фазовых траекторий:

Правила определения типа точки покоя - student2.ru

Таким образом, фазовые траектории лежат на прямых Правила определения типа точки покоя - student2.ru

Правила определения типа точки покоя - student2.ru

Пример 6.

Правила определения типа точки покоя - student2.ru

Так как определитель матрицы системы det A = 0, то система имеет бесконечно много положений равновесия. Все они лежат на прямой y = 2x.

Так как и след матрицы системы trA=0, то корни характеристического уравнения l1 = l2 = 0.

>> a=[2 -1;4 -2];

>> eig(a)

ans =

1.0e-015 *

0.0650 + 0.3140i

0.0650 - 0.3140i

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

Построим уравнения остальных фазовых траекторий:

Правила определения типа точки покоя - student2.ru

Таким образом, фазовые траектории лежат на прямых y = 2x + C, C – const, и параллельны прямой точек покоя.

Установим направление движения по траекториям следующим образом.

Определим знак производной y′(t) в точках оси x. при y = 0 имеем:

Правила определения типа точки покоя - student2.ru

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

Значит, движение по фазовым траекториям правее прямой точек покоя будет снизу вверх, а левее – сверху вниз.

Правила определения типа точки покоя - student2.ru

Правила определения типа точки покоя - student2.ru


Наши рекомендации