Дифференциальные уравнения. Учебное пособие-3
.pdfтворяющая однородным начальным условиям и, значит, равная нулю в некоторой окрестности точки x0 . Следовательно, решение СЛДУ –
тривиальное решение, т. е. Yi линейно зависимы и определитель Вронского равен нулю в любой точке отрезка [a, b].
Определение. Совокупность решений Y1, Y2 , , Yn системы (9)
называется фундаментальнойсистемойрешенийна [a, b], еслиопре-
делитель Вронского не обращается в нуль ни в одной точке [a, b].
Определение. Если для однородной системы ДУ (9) начальные условия определены следующим образом: Y (x0 ) E, то система реше-
ний Y1, Y2 , , Yn называется нормальной фундаментальной систе-
мой решений.
Замечание. Если Y1, Y2 , , Yn – фундаментальная система или нормальная фундаментальная система, то линейная комбинация
n
∑CiYi (x) – общее решение (9).
i 1
n
Теорема 2. Линейная комбинация ∑CiYi (x) линейно независимых
i 1
решений Yi , i 1, ..., n, однороднойсистемыДУ(9) снепрерывнымина отрезке [a, b] коэффициентами aij (x) будет общим решением (9) на
этом же отрезке.
Доказательство. Так как коэффициенты aij (x) непрерывны на
[a, b], то система удовлетворяет условиям теоремы о существовании и
единственности. Следовательно, для доказательства теоремы достаточно показать, что подбором постоянных Сi можно удовлетворить не-
которому произвольно выбранному начальному условию (7), т. е.
n
можно удовлетворить векторному уравнению∑ CiYi (x0 ) Y0 . Так как
i 1
n
∑CiYi (x) – общее решение (9), то система разрешима относительно Сi ,
i 1
посколькувсе Yi линейнонезависимыи W (x0 ) 0. Однозначноопределяем Сi ; так как Сi линейно независимы, то W (x) 0 x [a, b] .
21
Теорема 3. Если Y – решение неоднородной системы ДУ (8), а Y –
решение однородной системы ДУ (9), то Y + Y будет решением неоднородной системы (8).
Доказательство. По свойствам линейного оператора
L[Y Y ] L[Y ] L[Y ] F.
Теорема 4. Общее решение неоднородной системы ДУ (8) на отрезке [a, b] с непрерывными на этом отрезке коэффициентами aij (x)
и правыми частями fi (x) равно сумме общего решения соответствую-
щей однородной системы (9) и частного решения Y неоднородной системы (8).
Доказательство. Так как условия теоремы о существовании и един-
n
ственности выполнены, то остается доказать, что Y (x) ∑CiYi Y бу-
i 1
дет удовлетворять произвольно заданным значениям (7), т. е.
n |
|
Y (x0 ) Y0 ∑ CiYi (x0 ) Y (x0 ) . |
(11) |
i 1
Для системы (11) всегда можно определитьзначения Сi . Это можно сделать, так как Yi , i = 1, …, n, – фундаментальная система решений.
m
Теорема 5 (принцип суперпозиции). Решение СЛДУ L[Y ] ∑ Fi
|
|
|
m |
|
i 1 |
|
|
|
|
|
|
может быть представлено в виде ∑Yi , где каждое Yi |
удовлетворяет |
||||
|
|
|
i 1 |
|
|
|
|
|
m |
m |
m |
уравнению L[Yi ] Fi |
i |
1,m |
, тогда L[∑Yi ] ∑ L[Yi ] |
∑ Fi . |
|
|
|
|
i 1 |
i 1 |
i 1 |
Замечание. Принцип суперпозиций можно распространить на случай m при условии, что ряд, составленный из Yi , равномерно схо-
дится.
22
ОДНОРОДНЫЕ СИСТЕМЫ ЛИНЕЙНЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ С ПОСТОЯННЫМИ КОЭФФИЦИЕНТАМИ
Рассмотрим однородную систему дифференциальных уравнений (ОСДУ)
dY |
AY |
(12) |
|
dx |
|||
|
|
где A – матрица с постоянными коэффициентами (n n); Y – неизвест-
ная вектор-функция (n). Общее решение системы (12) имеет вид Y (x) exp(Ax)C, где C – вектор произвольных постоянных [3].
Согласно теореме о приведении матрицы к квазидиагональному
виду [3] матрица A может быть представлена в виде A SJS 1 , где J – жордановаформаматрицыA; S – матрица, составленнаяизсобственных и присоединенных векторов матрицы A. Тогда систему (12) можно за-
писать в виде |
dY |
SJS 1Y, или |
d(S 1Y ) |
JS 1Y. Сделаем замену пе- |
dx |
dx |
ременных Z S 1Y и запишем ОЛСДУ покомпонентно для случая, когда J – одна жорданова клетка размерностью (n n) :
dzdx1 z1 z2 ; dzdx2 z2 z3 ;
.... ....
dzn 1 zn 1 zn ; dx
|
dzn |
zn . |
|
|
|
dx |
|
|
|
|
|
|
|
|
Решая |
эту систему, |
получаем: zn (x) Cn exp( x) , |
zn 1(x) |
|
(Cn x Cn 1)exp ( x) , ..., |
z1(x) Cn xn 1 Cn 1xn 2 C1 |
exp ( x). |
||
Тогда общее решение ОЛСДУ: Y (x) SZ(x) .
23
Если все собственные значения матрицы A действительны и раз-
n |
|
личны, то общее решение ОЛСДУ Y (x) ∑ Ciei exp( i x) , где Ci |
– про- |
i 1 |
|
извольные постоянные; ei – собственные векторы матрицы A; i |
– соб- |
ственные значения матрицы A. |
|
ЕслиматрицаA вещественнаяисреди еесобственных значенийесть комплексно-сопряженные, то частными решениями ОЛСДУ будут не только решения Yj (x) ej exp(( p iq)x) , но и решения вида
Y1 j (x) Re{ej exp(( p iq)x)} и Y2 j (x) Im{ej exp(( p iq)x)} [1].
НЕОДНОРОДНЫЕ СИСТЕМЫ ЛИНЕЙНЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
МЕТОД ВАРИАЦИИ ПОСТОЯННЫХ
Рассмотрим неоднородную линейную систему ДУ (НЛСДУ):
|
|
|
dY |
AY F(x) . |
(13) |
|
|
|
dx |
||
|
|
|
|
|
|
Пусть Y (x) [Y |
Y |
Y ] |
(Y – частные решения ОЛСДУ) – фун- |
||
1 |
2 |
|
n |
i |
|
даментальнаясистемарешенийоднороднойсистемыДУ(12). Будемискать решение неоднородной системы (13) в виде Y (x) Y (x)C(x) , где C(x) – неизвестная вектор-функция. Подставляя решение в си-
стему (13), получаем соотношения для C(x) |
dC(x) |
Y |
1F(x) , тогда |
|
dx |
||||
C(x) ∫Y 1(x)F(x)dx C; |
C – вектор произвольных |
постоянных. |
||
Тогда общее решение неоднородной системы ДУ будет иметь вид
Y (x) Y (x) ∫Y 1(x)F (x)dx C .
24
МЕТОД КОШИ
Рассмотрим однородную систему ДУ (12) и сопряженную к ней систему:
|
|
dZ |
|
|
|
|
|
A Z. |
(14) |
|
|
dx |
||
Пусть 1(x), , n (x) |
– нормальная фундаментальная система ре- |
|||
шений СЛДУ (12), 1(x), |
, n (x) , – нормальная фундаментальная |
|||
системарешенийСЛДУ(14). Дляэтихсистемрешенийвыполненыусловия взаимной ортогональности: ( i (x), j (x)) 0, i j [1]. Это условие
выполняется для всех x [a, b] . Будем искать решение неоднородной
|
|
|
n |
|
|
|
|
|
|
|
системы ДУ (13) в виде Y (x) ∑ui (x) i (x) , где ui (x) |
– неизвестные |
|||||||||
|
|
i 1 |
|
|
|
|
|
|
||
функции. |
Подставляя решение |
|
в систему |
(13) |
и |
учитывая, что |
||||
i (x), i 1, |
|
|
|
|
|
|
n |
|
|
|
, n, – решения(12), получаем ∑ ui (x) i (x) F(x) . Умно- |
||||||||||
|
|
|
|
|
|
|
i 1 |
|
|
|
жая последнее соотношение скалярно на |
j (x) , получаем уравнения |
|||||||||
для определения функций ui (x) : |
|
|
|
|
|
|
|
|
||
|
u j (x) ∫ |
(F, j ) |
dx C j , |
j 1, , n . |
|
|||||
|
( j , j ) |
|
||||||||
|
|
|
|
|
|
|
|
|
||
Тогда общее решение НЛСДУ (13) будет иметь вид |
|
|
||||||||
|
|
n |
|
(F, j ) |
|
|
|
|
||
|
Y (x) ∑ j (x) ∫ |
dx |
C j |
. |
|
|||||
|
|
|
||||||||
|
j 1 |
|
( j , j ) |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
Замечание. Метод Коши применим для СЛНДУ с матрицей А, приводящейся к диагональному виду.
25
МАТРИЧНЫЕ ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ
Пусть дано матричное дифференциальное уравнение
dY (x)
dx
где
y11(x)
Y (x) y21(x)
yn1(x)
A(x)Y (x) ,
y12 (x) y1n (x) |
||
y22 (x) |
|
|
y2n (x) |
||
|
|
|
yn2 (x) |
|
|
ynn (x) |
||
(15)
(16)
с начальными значениями |
|
Y (x0 ) Y0 . |
(17) |
Теорема. Если матрица A(x) непрерывна на [a, b], а определитель матрицы Y0 0 , то на отрезке [a, b] существует единственное решение Y (x) уравнения (15) и определитель Вронского этого решения не обращается в нуль ни в одной точке [a, b].
Будем рассматривать одновременно с системой (15) систему вида
dZ(x) |
Z(x)A(x) . |
(18) |
|
dx |
|||
|
|
Определение. СЛДУ (18) называется сопряженной к СЛНДУ (15).
Теорема. Пусть Y (x) – решение (15), а матрица A(x) непрерывна
на [a, b], тогда Y 1(x) существует на [a, b] и является решением си-
стемы (16).
Доказательство: по предположению теоремы решение Y (x) на отрезке [a, b] существует и определитель этого решения не равен нулю в любой точке [a, b], следовательно, существует обратная матри-
26
ца Y 1(x) . А так как Y (x)Y 1(x) E, то, продифференцировав это соотношение, получаем:
dxd Y (x)Y 1(x) 0,
Y (x) dY 1(x) dY (x) Y 1(x) 0, dx dx
где Y (x) – решение уравнения (15), следовательно,
dY |
A(x)Y (x) Y (x) |
dY 1(x) |
A(x)Y (x)Y 1(x) 0, |
|
dx |
dx |
|||
|
|
dY 1(x) Y 1(x)A(x). dx
Рассмотрим систему n-линейно независимых неоднородных дифференциальных уравнений:
dY (x) |
A(x)Y (x) F(x) . |
(19) |
|
dx |
|||
|
|
Будем искать решение этой системы методом вариации постоянной. Обозначим Y (x) – решение (15). Будем искатьрешение системы (19)
в виде
|
|
|
|
Y (x) Y (x)C(x) , |
(20) |
||||
где C(x) – неизвестная вектор-функция. |
|
|
|||||||
Подставим (20) в (19): |
|
|
|
|
|||||
|
|
dY |
|
|
dC(x) |
|
|
|
|
|
|
dx C(x) Y (x) |
dx |
A(x)Y (x)C(x) F(x); |
|||||
|
dC(x) |
|
dC |
|
|
|
|
x |
|
Y |
F(x), |
Y 1(x)F(x), |
C(x) |
∫ Y 1( )F( )d C. |
|||||
dx |
|
||||||||
|
|
dx |
|
|
|
x0 |
|||
|
|
|
|
|
|
|
|
||
27
Тогда общее решение неоднородной системы будет иметь вид
x
Y (x) Y (x) ∫ Y 1
x0
( )F( )d C . (21)
Теорема. Общее решение неоднородной системы (19) можно представить в виде
x |
Y (x)Y 1 |
|
|
|
Y (x) (x) Y (x) ∫ |
( )F( )d |
, |
(22) |
|
|
|
|
|
|
x0 |
|
|
|
|
где (x) – общее решение соответствующей однородной системы.
Замечание. ПустьматрицаАпостояннаиначальныеусловияимеют
вид
dY |
AY (x), |
|
(23) |
dx |
Y (0) E.
Покажем, что решение CЛНДУ (23) – Y (x) удовлетворяет функци-
ональному уравнению |
|
|
|
Y (x ) Y (x)Y ( ) . |
|
(24) |
|
При любом фиксированном |
матрицы Y (x ) |
и Y (x)Y ( ) |
будут |
решением (23). При x 0 эти матрицы будут совпадать в силу начальных условий и, следовательно, по теореме о существовании и единственности они будут совпадать для любых x . Если воспользоваться (24), то можно записать(сделавзамену x x иумноживсоот-
ношение (24) справа на Y 1( )) :
Y (x)Y 1( ) Y (x )Y ( )Y 1( ) Y (x ) . |
(25) |
Подставляем (25) в(21) иполучаемвидобщегорешения НСДУспостоянной матрицей А:
x |
|
Y (x) (x) ∫ Y (x )F( )d . |
(26) |
x0 |
|
28
ПРИМЕР РЕШЕНИЯ СНЛДУ
1. Решить неоднородную систему дифференциальных уравнений методом вариации постоянной:
|
|
dY |
|
|
|
0 |
3 |
3 |
|
0 |
|
||
|
|
AY F(x) , где |
A |
1 |
8 |
6 |
, |
F (x) 0 |
. |
||||
|
|
|
|||||||||||
|
|
dx |
|
|
|
2 |
14 |
|
|
|
|
||
|
|
|
|
|
10 |
x |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Собственные значения матрицы A: 1 |
0, 2,3 |
1. Собственный век- |
|||||||||||
тор для |
– |
e (2, 1, |
1)T |
, собственный и присоединенный векторы |
|||||||||
|
|
1 |
|
|
1 |
4)T , |
e (3, 0, 2)T . Тогда |
|
|
||||
для |
2 |
– |
|
e |
(3, 3, |
матрица перехода |
|||||||
|
|
2 |
|
3 |
|
|
|
|
|
|
|||
S
J
dY
dx
2 |
3 |
3 |
|
|
|
|
|
1 |
3 |
|
и жорданова |
форма |
матрицы A имеет вид |
|
0 |
|||||
|
1 |
4 |
2 |
|
|
|
|
|
|
|
|
|
|
0 |
0 |
0 |
|
|
1 |
|
|
|
1 |
|
Сделаем замену |
Z S |
|
0 |
1 . |
Y в однородной СЛДУ |
||||
0 |
0 |
1 |
|
|
|
|
|
|
|
|
|
|
|
AY и получим систему линейных дифференциальных уравнений:
dz1 |
0, |
dz2 |
z |
2 |
z , |
dz3 |
z . |
|
|
|
|||||
dx |
dx |
3 |
dx |
3 |
|||
|
|
|
|||||
Общее решение однородной системы:
z1(x) C1, |
z2 (x) (C3x C2 )exp( x), |
z3 (x) C3 exp( x) . |
||||||||||
|
2 |
3 |
3 |
|
|
|
C1 |
|
|
|||
Y (x) SZ (x) |
1 |
3 |
0 |
(C |
x C |
2 |
)exp ( x) |
|
|
|||
|
|
1 |
4 |
2 |
|
3 |
|
|
|
|
|
|
|
|
|
|
C exp ( x) |
|
|
||||||
|
|
|
|
|
|
|
3 |
|
|
|
|
|
29
|
2 |
|
3e x |
|
3(x 1)e x |
C1 |
|
|||||
|
|
|
3e x |
|
|
3xe x |
|
|
||||
|
1 |
|
|
|
|
|
C |
2 |
. |
|||
|
|
|
|
x |
|
|
x |
|
|
|
||
|
1 |
4e |
( 4x 2)e |
|
|
|
|
|||||
|
|
|
C3 |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Фундаментальная система решений |
|
|
|
|
|
|
||||||
|
|
|
2 |
|
3e x |
3(x 1)e x |
||||||
|
|
|
|
|
3e x |
3xe x |
|
|
|
|||
Y (x) |
1 |
|
|
|
. |
|||||||
|
|
|
|
|
4e x |
( 4x 2)e x |
|
|||||
|
|
1 |
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Тогда решение неоднородной |
системы |
|
ДУ: Y (x) Y (x)C(x) , |
|||||||||
C(x) ∫Y 1(x)F(x)dx . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
9 |
|
|
|
||||
|
|
|
|
|
|
|
|
|
||||||
|
|
|
x 2 |
|
|
5x |
7 |
|
|
|||||
1 |
(x)F(x) |
|
|
x |
|
x |
||||||||
Y |
|
3 |
|
|
e |
|
3 |
|
|
e |
|
|||
|
|
|
1 |
|
x |
|
5 |
|
x |
|
|
|||
|
|
|
e |
|
e |
|
|
|||||||
|
|
|
|
3 |
|
|
|
|
3 |
|
|
|
||
Тогда
3 |
|
|
|
|
|
0 |
|
|
|
|
|
|
||
(1 x)ex |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
ex |
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3x |
|
|
|
|
(1 |
x)ex . |
|
|
xex |
|
|
|
|
C1(x) 3∫ xdx 32x2 C1 ,
C2 (x) ∫(x x2 )exdx (x2 2x 2)ex C2 , C3(x) ∫ xexdx (x 1)ex C3 .
30
