Материал: Обыкновенные дифференциальные уравнения с приложениями к задачам механики, физики, термодинамики и экологии. Ряжских В.И., Бырдин А.П

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

 

 

 

 

(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

Источник: https://studfile.net/preview/16568742/