|
|
|
|
(1) |
|
(2) |
|
(1) |
0 |
|
|
|
|
|
|
1 |
1 |
2 |
|
||||
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
(2) |
|
|
(2) |
0 |
|
|
|
|
|
|
1 |
2 |
|
||||
|
|
|
|
(1) |
|
(1) |
|
|
(2) |
0 |
|
|
|
|
|
|
|
||||||
|
|
|
|
1 |
|
2 |
|
2 |
|
|
|
|
|
|
|
|
(2) |
|
|
(2) |
|
|
|
|
|
|
|
|
|
|
|
0. |
|
||
|
|
|
|
|
1 |
2 |
|
||||
|
Следовательно, из искомых постоянных две произвольны, |
например, |
|||||||||
2(2) |
C2 |
1(1) |
C1, поэтому из этой системы для остальных 1(2), |
2(1) будет |
|||||||
2(1) C1 C2. Отсюда общее решение исходной системы будет 1(2) C2. y C1ex C2xex
z( C1 C2)ex C2xex.
3.Если среди корней характеристического уравнения (5.63) есть
комплексное |
спряженные |
пары |
|
m m i m, |
m 1 m i m , то |
этим |
||||||||||||||||||
корням будут соответствовать частные решения |
|
|
||||||||||||||||||||||
|
|
m(i) |
|
|
ym(1) m(1) |
e( m i m), |
ym(2) m(2) e( m i m), |
|
||||||||||||||||
причем |
определяются |
из |
|
|
систем |
(5.63). |
Можно показать, |
что |
||||||||||||||||
действительные и мнимые части |
также являются решениями системы (5.57). |
|||||||||||||||||||||||
Таким образом, можно записать два частных решения: |
|
|||||||||||||||||||||||
|
|
|
|
|
|
y |
m(1) |
|
e m (m1)cos mx (m2)sin mx |
|
||||||||||||||
|
|
|
|
|
|
|
|
m(2) |
e m |
|
|
(m2)sin mx , |
|
|||||||||||
|
|
|
|
|
|
|
y |
m cos mx |
|
|||||||||||||||
где m(1), |
m(2), |
|
|
(m1), |
|
|
m(2) |
- действительные коэффициенты, связанные с |
||||||||||||||||
|
|
|||||||||||||||||||||||
m(1), |
m(2). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Следовательно, можно составить фундаментальную систему решений |
|||||||||||||||||||||||
только из действительных решений. |
|
|
|
|
|
|
|
|||||||||||||||||
|
Пример 3. |
|
|
Дана |
система |
|
дифференциальных уравнений |
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
dy1 |
|
|
7y |
y |
2 |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
dx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dy2 |
|
2y1 |
5y2. |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
dx |
|
|
|
|
|
|
|
|
|
|
|
|
Найдем ее общее решение.
Решение. Характеристическое уравнение для этой системы имеет вид
|
7 |
1 |
0 |
или |
|
2 |
12 |
|
37 |
|
0, |
|
2 |
5 |
|||||||||
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
||
корни этого уравнения 1 6 i, |
2 6 i. |
|
|
|
|
|
|
||||
Подставляя эти корни в (5.63), получим |
|
|
|
|
|
|
|||||
|
|
|
|
150 |
|
|
|
|
|
|
|
1(1) 1, 2(1) 1 i.
Соответствующее частное решение yij (x) будет
y1(1) e( 6 i)x; y2(1) (1 i)e( 6 i)x .
Для 2 6 i, получим соответственно
1(2) 1, 2(2) 1 i;
y1(2) e( 6 i)x; y2(2) (1 i)e( 6 i)x.
Поскольку из четырех получившихся функций в общее решение для1, 2 должны войти лишь два линейно независимых решения с произвольными постоянными C1, C2, преобразуем получившиеся решения следующим образом
y(1) |
e 6x (cosx isin x); |
1 |
|
y2(1) |
(1 i)e 6x(cosx isin x); или |
y(1) |
e 6x cosx ie 6x sin x; |
1 |
|
y(1) |
e 6x(cosx sin x) ie 6x(cosx sin x); |
2 |
|
для |
y1(2), y2(2) |
будет |
|
|
|
|
|
|
|
|
|||||||
y(2) |
e 6x cosx ie 6x sin x), |
|
|
|
|
|
|
|
|
||||||||
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y2(2) |
e 6x (cosx sin x) ie 6x (cosx sin x). |
|
|
|
|||||||||||||
Поэтому составляя общее решение из действительных и мнимых частей |
|||||||||||||||||
этих решений для yij (x).найдем |
|
|
|
|
|
|
|
|
|||||||||
|
|
|
y |
|
C e 6x |
cosx C |
2 |
e 6x sin x. |
|
|
|
||||||
|
1 |
|
1 |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
y |
2 |
C e 6x |
(cosx sin x) C |
2 |
e 6x (cosx sin x). |
|
||||||||
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
n-го |
||
Известно, что линейное однородное дифференциальное уравнение |
|||||||||||||||||
порядка |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dn y |
|
|
|
dn 1y |
|
|
dn 2 y |
|
dy |
|
|||||
|
|
|
an 1 |
(x) |
|
|
an 2 |
(x) |
|
|
... a1(x) |
|
a0y |
(5.65) |
|||
|
|
|
|
|
dxn 2 |
|
|||||||||||
|
|
dxn |
|
|
|
dxn 1 |
|
|
|
dx |
|
||||||
эквивалентно линейной, однородной системе
|
dy0 |
y |
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
||||
dx |
1 |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dy1 |
y2 |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
dx |
|
|
|
|
|
|
|
|
|
|
|
(5.66) |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
dyn 1 |
|
|
|
|
|
|
|
|
|
|
|
|||
|
a |
y |
0 |
a y |
... a |
n 1 |
y |
n 1 |
. |
|
|||||
|
|
||||||||||||||
|
dx |
|
|
0 |
|
1 1 |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
151 |
|
|
|
|
|
|
Будем считать все ak const.
Характеристическая матрица для (5.66) в этом случае будет иметь вид
1 |
0 |
|
0 |
|
|
||
|
0 |
|
1 |
|
0 |
|
|
|
|
|
|
||||
|
|
|
|
|
|
||
|
|
||||||
|
a |
0 |
a |
|
a |
n 1 |
|
|
|
1 |
|
|
|
||
или
det( ) n an 1 n 1 ... a1 a0 0.
Если же вычеркнуть из матрицы первый столбец и последнюю строку,
получим определитель равный |
+I или –I. Следовательно, характеристический |
||||
многочлен, если имеет корень k |
кратности |
Sk , |
то матрица |
имеет |
|
элементарный делитель ( k)Sk . Поэтому корню |
k будут соответствовать |
||||
Sk линейно независимых решений вида |
|
|
|
|
|
(C0 C1x ... CSk 1 xSk 1)e kx. |
|
|
(5.67) |
||
Таким образом, методы |
решения |
систем |
типа |
(5.58) могут |
быть |
использованы для решения линейных дифференциальных уравнений n-го порядка, разрешенных относительно производных.
Пример 4. Найти общее решение дифференциального уравнения y 3y 3y y 0,
Обозначим y y0, |
y y0 |
y1, |
|
y |
y1 |
y2, |
y y2. |
||||||
Тогда |
|
dy0 |
|
|
|
|
|
|
|
|
|
||
|
|
|
y , |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|||||
|
|
|
dx |
1 |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
||
|
dy |
|
|
|
|
|
|
|
|
|
|||
|
|
1 |
|
y2 |
, |
|
|
|
|
(5.68) |
|||
|
|
|
|
|
|
|
|
||||||
|
dx |
|
|
|
|
|
|
|
|
||||
|
|
dy2 |
|
y |
0 |
3y 3y |
2 |
. |
|
||||
|
|
|
|
||||||||||
|
|
|
dx |
|
|
|
1 |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
||
Характеристическое уравнение системы (5.68) имеет вид
|
|
1 |
0 |
|
|
|
|
0 |
|
1 |
|
0, |
т.е |
|
1 |
3 |
3 |
|
|
|
2(3 ) 3 1 0 |
или |
( 1)3 0. |
||||
Поскольку корень 1 1 |
имеет кратность S 3, то, учитывая (5.67), |
|||||
получим общее решение (5.68), а следовательно, и исходного уравнения третьего порядка в виде
(C0 C1x C2x2)ex y(x).
152
Для нахождения общего решения нормальных систем дифференциальных уравнений вида
|
dy1 |
|
|
f (x, y ,..., y |
|
|
), |
|
|||
|
|
|
|
||||||||
dx |
1 |
1 |
n |
|
|
|
|||||
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
||
dy2 |
|
|
f2(x, y1,..., yn), |
|
|||||||
|
|
|
|
|
|||||||
|
|
|
|||||||||
dx |
|
|
|
|
|
|
|
|
|
(5.69) |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
dyn |
|
|
f |
n |
(x, y ,..., y |
n |
) |
|
||
|
|
|
|
||||||||
|
|
|
|
1 |
|
|
|
||||
dx |
|
|
|
|
|
|
|
|
|
|
|
применяется также метод последовательного исключения неизвестных функций, то есть сведение системы (5.69) к уравнению вида (5.65). Для этого, дифференцируя первое уравнение из (5.69) по x, получим
|
d2 y |
|
|
f |
1 |
|
|
|
f |
|
|
|
|
dy |
|
|
|
|
f |
1 |
|
dy |
2 |
|
|
|
f |
1 |
|
dy |
n |
|
||||
|
1 |
|
|
|
|
|
|
1 |
|
|
1 |
|
|
|
|
... |
|
|
|
. |
||||||||||||||||
dx2 |
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
y1 dx |
|
|
|
|
y2 dx |
|
|
yn |
dx |
|||||||||||||||||||||
Используя тот факт, |
что |
|
dyk |
|
fk , |
|
|
найдем |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
dx |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
d2 y |
|
|
|
f |
1 |
|
f |
1 |
|
|
|
|
|
|
f |
1 |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
f1 ... |
|
|
|
fn. |
|
|
|
|
(5.70) |
|||||||||
|
|
|
|
|
|
dx2 |
|
|
x |
y1 |
|
yn |
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
Определим xn из первого уравнения системы (5.69) |
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
yn g(x, y1, y2,...,yn 1, y1) |
|
|
|
|
|
|
||||||||||||||||||||||
Подставим yn |
|
|
в (5.70) |
|
и, |
таким образом, |
исключим эту неизвестную |
|||||||||||||||||||||||||||||
функцию из (5.70). Продолжая эту процедуру, можно получить уравнение n-го
порядка для одной неизвестной |
функции. |
|
|
||||||||
|
Пример 5. Найти |
общее решение |
|
|
|||||||
|
|
|
|
|
|
dy |
z, |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
dx |
|
|
|
||
|
|
|
|
|
|
|
dz |
|
y. |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
dx |
|
|
|
||
|
Решение. |
d2 y |
|
dz |
, т.е. |
y |
y 0 . |
||||
|
dx2 |
dx |
|||||||||
|
|
|
|
|
|
|
|
|
|
||
Характеристическое уравнение имеет вид |
|
|
|||||||||
|
|
|
|
|
k2 1 0; |
k i, |
|||||
|
|
|
|
|
|
|
|
|
1,2 |
||
|
|
|
|
|
y(x) C eix C |
2 |
e ix |
||||
|
|
|
|
|
|
1 |
|
|
|
||
или |
y(x) C1(cosx isin x) C2(cosx isin x), тогда общее решение можно |
||||||||||
записать |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
153 |
|
|
|||
y(x) C1 cosx C2 sin x, z(x) C2 cosx C1 sin x.
5.6.Физические задачи, приводящие к уравнениям 2-го порядка
исистемам дифференциальных уравнений
1.Уравнение колебаний математического маятника.
Пусть шарик массы m закреплен на конце B невесомого стержня, длина которого l. Стержень шарнирно закреплен в точке O так, что колебания маятника происходят в одной плоскости. На шарик действуют две силы: сила тяжести, направленная вертикально вниз, F mg; направленная по радиусу BO сила реакции стержня T (рис. 13).
O
l
T
B
Рис. 13
Обозначим a - проекцию ускорения на направление касательной к траектории маятника. Возвращая в положение равновесия, сила – проекция силы тяжести на то же направление, равна - mgsin .
|
По второму закону Ньютона |
|
|
|
|
|||||
|
|
|
|
ma |
mgsin . |
(5.71) |
||||
|
При малых |
углах отклонения sin |
s |
|
(по теореме синусов для |
|||||
l |
||||||||||
AOB), |
где s - |
|
|
|
|
|
||||
смещение маятника, измеряемое |
длиной дуги. Поскольку |
|||||||||
a |
d2s |
, то уравнение (5.71) можно записать в виде |
|
|||||||
|
|
|||||||||
|
|
|
|
|
|
|
|
|
||
|
dt2 |
|
|
d2s |
|
|
|
|
|
|
|
|
|
|
2s 0, |
(5.72) |
|||||
|
|
|
|
dt2 |
||||||
|
|
|
|
|
|
|
|
|||
где 2 g . Уравнение (5.72) является уравнением движения математического l
154