Все решения данной системы будут асимптотически устойчивы тогда и только тогда, когда матрица
|
a |
a |
... |
a |
|
|
11 |
12 |
... |
1n |
|
A |
a21 |
a22 |
a2n |
||
... ... |
... |
... |
|
||
|
|
|
|
|
|
|
|
an2 |
... |
|
|
|
an1 |
ann |
|||
будет иметь только собственные значения с отрицательной действительной частью.
В случае n 3 целесообразно находить собственные значения матриц в среде MathCAD (или в другом математическом пакете).
Достаточно ввести матрицу коэффициентов системы
|
a |
a |
... |
a |
|
|
11 |
12 |
... |
1n |
|
A: |
a21 |
a22 |
a2n |
||
... ... |
... |
... |
|
||
|
|
|
|
|
|
|
|
an2 |
... |
|
|
|
an1 |
ann |
|||
и получить вектор собственных значений матрицы (включая комплексные собственные значения) с помощью функции
eigenvals(A)
3. Устойчивость по линейному приближению (линеаризация системы).
Пусть дана система нелинейных ОДУ 1-го порядка с постоянными ко-
эффициентами
y1 |
f1 y1, y2,...,yn , |
|
|
|
|
|
|
|
|||||||
|
|
|
|
y1, y2,...,yn , |
|
|
|
|
|
|
|
||||
y2 |
f2 |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|||||||||
.................................. |
или в векторной форме y F y , |
||||||||||||||
|
|
|
|
y , y |
|
|
|
. |
|
|
|
|
|
|
|
y |
|
f |
n |
2 |
,...y |
n |
|
|
|
|
|
|
|
||
|
n |
|
1 |
|
|
|
|
|
|
|
|
|
|||
для которой выполняется F 0 0.
Построим матрицу
16
|
f |
|
|
|
|
f |
|
|
|
|
|||||
0 |
|
0 |
|||||||||||||
|
1 |
|
|
|
... |
|
1 |
|
|
|
|
|
|
|
|
y1 |
yn |
|
|||||||||||||
|
... |
|
|
|
|
||||||||||
A ... |
... |
|
|
|
|
. |
|||||||||
|
fn 0 |
|
|
fn |
|
|
|||||||||
|
... |
0 |
|||||||||||||
|
|
|
|
|
|
||||||||||
|
y |
|
|
y |
n |
|
|
||||||||
|
1 |
|
|
|
|
|
|
|
|
||||||
Если:
1) действительные части всех собственных значений матрицы A отри-
цательны, то начало координат является асимптотически устойчивым поло-
жением равновесия;
2) хотя бы одно собственное значение матрицы A имеет положитель-
ную действительную часть, то начало координат является неустойчивым по-
ложением равновесия; 3) в остальных случаях исследовать нелинейную систему на устойчи-
вость с помощью линеаризации нельзя.
Замечание. Если fi представляют собой многочлены, то для линеари-
зации системы достаточно удалить все члены, степень которых превосходит
1.
2 Разностные уравнения
2.1 Линейные разностные уравнения второго порядка с постоян-
ными коэффициентами
Определение. Линейным разностным уравнением второго порядка с
постоянными коэффициентами называется уравнение вида
y n 2 py n 1 qy n f n q 0 .
17
Таблица 8 – Общее решение однородного разностного уравнения вто-
рого порядка с постоянными коэффициентами |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
Разностное |
|
|
|
y n 2 py n 1 qy n 0 |
|
|
||||
уравнение |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Характери- |
|
|
|
|
|
|
|
|
|
|
стическое |
|
|
|
|
2 p q 0 |
|
|
|
||
уравнение |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Корни |
Два простых |
действи- |
Один |
действитель- |
Два комплексно |
сопряжен- |
||||
|
тельных корня |
|
|
ный |
корень |
второй |
ных корня |
|
|
|
|
1, 2 |
|
кратности |
|
1,2 |
i |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Общее ре- |
y n C n |
C n |
y C n C n |
y n C cosn C sinn |
||||||
шение |
1 |
1 |
2 |
2 |
|
1 |
2 |
1 |
|
2 |
Таблица 9 - Общее решение линейного неоднородного разностного
уравнения второго порядка с постоянными коэффициентами и специальной
правой частью |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Линейное неоднородное раз- |
|
|
|
|
|
|
|
|
|
ностное уравнение второ- |
y n 2 py n 1 qy n f n |
q 0 |
|
||||||
го порядка с постоянными |
|
||||||||
коэффициентами и специаль- |
|
|
|
|
|
|
|
|
|
ной правой частью |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Характеристическое уравне- |
2 |
p q 0 |
|
|
|
|
|
||
ние |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
f n |
Корни характеристи- |
|
|
Вид частного решения |
|
||||
|
ческого уравнения |
|
|
|
|
|
|
|
|
|
1 не является корнем |
|
~ |
n Am n |
m |
|
Am 1 n |
m 1 |
|
|
характеристического |
|
Pm |
|
|
||||
|
|
... A1 n A0 |
|
|
|
|
|||
Pm n |
уравнения |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
1 является корнем |
|
|
|
|
|
|
|
|
|
|
характеристического |
|
|
ns P~ n |
|
|
|||
|
уравнения кратности |
|
|
|
|
m |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
s |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
18 |
|
|
|
|
|
|
|
|
Продолжение таблицы 9
Линейное |
неоднородное |
раз- |
|
|
|
|
|
|
|
|
|
|
|
ностное |
уравнение второго |
|
y n 2 py n 1 qy n f n |
q 0 |
|
||||||||
порядка с постоянными коэф- |
|
|
|||||||||||
фициентами и |
специальной |
|
|
|
|
|
|
|
|
|
|
||
правой частью |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Характеристическое уравнение |
|
|
2 |
p q 0 |
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|||||
|
f n |
|
Корни характеристи- |
|
|
|
Вид частного решения |
|
|||||
|
|
|
|
ческого уравнения |
|
|
|
|
|
|
|
||
|
|
|
|
|
не является кор- |
|
|
|
|
|
|
|
|
|
|
|
|
нем |
характеристиче- |
|
|
|
A n |
|
|
|
|
|
|
|
|
ского уравнения |
|
|
|
|
|
|
|
||
|
c n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
является |
корнем |
|
|
|
|
|
|
|
||
|
|
|
|
характеристического |
|
|
|
Ans n |
|
|
|||
|
|
|
|
уравнения кратности |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
s. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
cos isin |
|
|
|
|
|
|
|
||
|
|
|
|
не |
является |
корнем |
|
|
|
Acosn Bsinn |
|
||
|
|
|
|
характеристического |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|||
acosn , |
asinn , |
|
уравнения |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|||
acosn bsinn |
|
cos isin |
|
|
|
|
|
|
|
||||
|
является корнем ха- |
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
рактеристического |
|
|
|
n(Acosn Bsinn ) |
|
||||
|
|
|
|
уравнения |
первой |
|
|
|
|
|
|
|
|
|
|
|
|
кратности |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
(cos isin ) |
|
n |
|
~ |
~ |
|
|
||
|
|
|
|
не |
является |
корнем |
|
Ps n cosn Qs n sinn |
|
||||
|
|
|
|
характеристического |
|
|
|
, |
|
|
|
||
|
|
|
|
|
|
|
где s max m,r |
|
|||||
|
|
|
|
уравнения |
|
|
|
|
|
||||
n P n cosn Q n sinn |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
||||
(cos isin ) |
|
|
|
|
|
|
|
||||||
m |
|
r |
|
|
|
|
~ |
~ |
n sinn |
|
|||
|
|
|
|
является корнем ха- |
n |
n |
|||||||
|
|
|
|
|
P n cosn Q |
||||||||
|
|
|
|
|
|
|
|
|
|
s |
s |
|
|
|
|
|
|
рактеристического |
|
|
|
, |
|
|
|
||
|
|
|
|
уравнения |
первой |
|
|
|
где s max m,r |
|
|||
|
|
|
|
кратности |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
19
2.2 Элементы теории устойчивости для разностных уравнений
Пусть дано линейное разностное уравнение m-го порядка с постоянны-
ми коэффициентами
y n m am 1y n m 1 ... a1y n 1 a0y n f n .
Тогда:
1) если все корни характеристического уравнения удовлетворяют нера-
венству 1, то все решения данного уравнения будут асимптотически ус-
тойчивы; 2) если все корни характеристического уравнения удовлетворяют нера-
венству 1, причем все корни 1 являются простыми, то все решения данного уравнения будут устойчивы;
3) все решения данного уравнения будут неустойчивы во всех осталь-
ных случаях.
3 Уравнения в частных производных
3.1 Уравнения в частных производных первого порядка
Определение. Линейным уравнением в частных производных (УЧП)
первого порядка называется уравнение вида
A |
u |
B |
u |
C, |
(2) |
x |
|
||||
|
|
y |
|
||
где A, B и C - функции от x, y и u; u u x,y - неизвестная функ-
ция. |
|
|
|
|
|
|
|
Пусть |
система ОДУ |
dx |
|
dy |
|
du |
имеет общие интегралы |
|
|
|
|||||
|
|
A B C |
|||||
F1 x,y,u C1 |
и F2 x,y,u C2, тогда общий интеграл УЧП (2) имеет вид |
||||||
|
F1 x,y,u ,F2 x,y,u 0. |
||||||
|
20 |
|
|
||||