Добавил:
хачю сдать сессию Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

лекции вышмат

.pdf
Скачиваний:
93
Добавлен:
01.06.2024
Размер:
5.41 Mб
Скачать

ЛЕКЦИЯ 10

ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ ПЕРВОГО ПОРЯДКА (продолжение)

Линейные уравнения первого порядка это уравнения вида

y'+P(x) y = Q(x) .

(1)

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

Первый способ (метод Бернулли). Будем искать решение уравнения (1) в виде

y

=

uv

, где

u = u(x) ,

v = v(x) – некоторые функции. Тогда

y'= u'v +uv' и (1) принимает

вид

u'v +uv'+P(x)uv = Q(x) . Перепишем последнюю формулу следующим образом:

 

 

u'v +u(v'+P(x)v) = Q(x)

(2)

Теперь выберем функцию v(x) такой, чтобы

 

 

 

v'+P(x)v = 0

,

(3)

а затем найдем все функции u(x) , при которых справедливо равенство (2) (т.е. при

нахождении решения в виде произведения двух сомножителей мы выбираем один из этих сомножителей как нам удобно, а затем находим второй сомножитель так, чтобы их произведение было решением).

Разделяя переменные, имеем:

dv

= −P(x)v ;

dv

= −P(x)dx ;

dv

= − P(x)dx + c ; так

 

 

 

 

 

dx

 

 

v

 

 

v

 

 

 

достаточно найти одно решение уравнения (3), то возьмем в последней формуле

c

Далее из (2) и (3) имеем: u'v = Q(x) ; u'=

Q(x)

 

(v(x) 0) ; u =

 

Q(x)dx + c .

 

v(x)

 

 

 

 

 

 

 

v(x)

 

 

 

 

 

 

 

 

 

как нам

= 0 .

Пример.

Решить дифференциальное уравнение y '

3

= 0 .

 

y + x

 

 

 

 

x

 

 

Решение. Ищем

решение уравнения в виде

y = uv .

Тогда

y '

= u ' v + uv ' , и уравнение

 

3

;

 

 

 

3

 

 

 

. Потребуем, чтобы

 

3

, тогда

dv

 

u ' v + uv '

uv + x = 0

u ' v + u

 

v '

 

v

 

+ x = 0

v '

v = 0

 

=

 

x

 

 

 

 

x

 

 

 

 

 

x

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

принимает вид

3

dv

 

3dx

v ;

 

=

;

x

v

 

x

 

dv

= 3

v

 

 

 

 

y =

 

1

+

 

 

 

 

x

 

dx ; ln |

x

 

3

=

c x

 

 

 

v |= 3 ln

cx

3

+ x

2

 

 

|

.

x

|

; v = x3 . Теперь u ' x3 + x = 0 ;

u ' = −

1

 

2

 

x

 

 

;отсюда

u = −

dx

+ c =

x

2

 

 

 

1 x

+ c

,

и

Второй способ (метод вариации произвольной постоянной). Рассмотрим так называемое линейное однородное уравнение, соответствующее данному уравнению (1):

 

 

y'+P(x) y = 0 .

 

 

(4)

 

 

 

 

 

 

dy

 

dy

= −P(x)dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решаем это уравнение (с разделяющимися переменными): dx = −P(x) y;

y

;

 

 

dy

= − P(x)dx + c ; ln | y |= − P(x)dx + c ; | y |= eP( x)dx+c = eP( x)dxec ;

y = ec eP ( x)dx , или

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y = ce

P( x)dx

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

(5)

 

 

 

 

 

 

 

 

 

Теперь будем искать решение уравнения (1) по той же формуле (5), считая, что в ней c = c(x) (отсюда и название метода). Тогда

y'= c'eP( x)dx + c(eP( x)dx )'= c 'eP( x)dx + ceP( x)dx (P(x)dx ) = c 'eP( x)dx ceP( x)dx P(x) .

Подставляя эту производную в (1), имеем:

160

 

 

 

c'e

P( x)dx

cP(x)e

P( x)dx

+ P(x)ce

P( x)dx

= Q(x) , или c'e

P( x)dx

= Q(x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(т.е. члены с

c

всегда сокращаются, остается только член с c' ). Отсюда c'= Q(x)e

 

P( x)dx

и,

 

 

 

интегрируя, мы находим с=с(х).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

) y

 

 

 

 

 

 

 

 

2

)

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1 + x

'2xy = (1 + x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример.

Решить задачу Коши

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(0) = 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Переписав это уравнение в виде

 

y '

2x

y = 1 + x

2

, мы видим, что оно действительно является линейным.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 + x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Соответствующее

 

однородное уравнение

имеет

вид

 

y '

 

2x

 

y

= 0 . Решаем это

уравнение:

dy

=

2x

 

y ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 + x

2

 

 

 

 

 

 

 

 

dx

 

 

1 + x

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy

 

2x

 

 

dy

 

 

 

2x

 

 

 

 

 

 

 

2

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

2 dx

;

=

 

2 dx + c ; ln |

y

|=

d (1 + x

+ ln

|

 

c |=

 

ln | 1 + x2 | + ln | c |= ln | c(1 + x2 ) | ;

 

 

 

 

 

 

y

1 + x

 

y

1

+ x

1

+ x

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

| y |=| c(1 + x2 ) | ;

y = c(1 + x2 ) , или y = c(1 + x2 ) .

Теперь ищем решение исходного уравнения по этой же формуле,

считая, что в ней c = c(x) . Подставляя в уравнение, имеем:

c '(1 + x2 ) + c 2x

2x

 

c(1 + x2 ) = = 1 + x

2 ; c '

= 1 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 + x2

 

 

 

 

 

 

 

 

 

c =

 

dx + c = x + c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

, и ( c заменяем на c )

y = (x + c)(1 + x

) . Подставляя сюда х = 0,

имеем: 2 = с

и y = (x + 2)(1 + x

) .

Уравнения Бернулли это уравнения вида

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y '+ P(x) y = Q(x) y

α

,

 

 

 

 

 

 

 

 

 

 

(6)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

α 0 ,

α 1

(при

α = 0 получаем линейное уравнение, а при

α =1

уравнение с

разделяющимися переменными).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

Пример. Решить дифференциальное уравнение

Решение.

1) y = uv

;

u ' v + uv '

4

uv ;

u ' v +

uv = x

 

 

 

 

 

 

 

 

x

 

 

dv

 

dx

 

4

 

 

 

 

 

 

 

= 4

;

v = x

.

 

 

 

 

 

 

 

 

 

 

 

 

v

 

x

 

 

 

 

 

 

 

 

y '

4

y =

 

 

x

 

 

 

4

u(v '

v)

 

 

x

x

y

= x

 

.

uv

. Потребуем, чтобы

v '

4 v

x

=

0

, тогда

Тогда

u ' x

4

 

2) Решаем

 

 

2

 

u ' x =

du

=

dx

= x

u x

;

u ;

 

 

 

 

 

 

 

 

 

u

 

x

соответствующее однородное

; 2

u = ln |

уравнение

x | + ln

y ' =

4

 

 

x

| c

y .

|= ln | cx |

Имеем:

; u

dy dx

=

=

1 ln

4

4 y

x

2 | cx | ; y = uv =

1

x4 ln 2 | cx | .

 

 

 

 

 

 

 

4

 

 

 

 

 

;

dy

=

4

dx ;

 

dy

= 4

 

dx

+ c ;

 

 

 

 

 

y

 

x

 

 

 

y

 

x

ln |

y |= 4 ln |

x | + ln |

c | ;

ln |

y |= ln | cx

4

| ; y = cx

4

;

 

4

. Теперь ищем решение исходного уравнения по

этой

 

 

y = cx

 

 

 

 

 

 

 

 

c = c(x) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

3

 

4

 

 

4

 

 

 

4

 

 

 

 

формуле,

считая,

что в

ней

 

Подставляя

в

уравнение,

имеем:

 

c ' x

+ c4x

cx

= x cx

; c ' x =

 

c ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

dc

=

dx

;

dc

=

dx

+ c ;

2 c = ln |

x | + ln | c |= ln

| cx | ;

c =

1

ln

2

| cx

|

;

y = c( x) x

4

 

=

1

x

4

ln

2

| cx | .

 

 

 

 

 

c

x

c

x

4

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Уравнения в полных дифференциалах это уравнение вида

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P(x, y)dx +Q(x, y)dy = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(7)

в котором левая часть является (полным) дифференциалом некоторой функции двух переменных, т.е. существует функция u(x, y) , такая, что du(x, y) = P(x, y)dx +Q(x, y)dy .

161

Пример. Решить дифференциальное уравнение

(x

3

2

x)dx +

 

y

( y

3

 

x

2

y)dy

 

=

0

.

Решение.

P(x,

y)

плоскости 0xy .

 

P( x, y)

= −2 yx

 

 

 

,

y

 

 

= x

3

y

2

x ,

Q(x, y) = y

3

x

2

y

непрерывны вместе со

своими частными

производными

на всей

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

я

 

Q( x, y)

= −2xy , т.е.

 

P( x,

y)

Q( x, y)

, значит,

уравнение является

уравнением в

полных

 

 

 

 

 

=

 

 

 

 

 

 

 

y

 

 

 

 

x

 

 

 

 

 

 

x

 

 

 

дифференциалах

u

= P = x

3

y

2

x

 

 

 

 

 

 

 

 

 

y

4

 

 

 

 

c( y) =

 

u =

 

 

 

 

4

 

 

 

 

(кстати, отметим,

x u =

(x

3

y

2

x

4

 

x

2

y

2

 

 

y

4

 

 

 

 

+

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

2

 

 

 

4

 

что оно также является однородным). Далее имеем:

x)dx = x

4

y

2

x

2

u

 

 

 

 

 

 

 

 

 

 

 

 

 

+ c( y)

= −x2 y + c ( y) = Q = y3

 

 

 

 

4

 

2

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

4

 

x

2

y

2

 

y

4

 

 

 

 

 

 

 

 

 

 

 

+

 

=

Общий интеграл уравнения имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

2

 

 

4

 

x2

C

yc ( y) = y3

x4 2x2 y2 + y4

=

C

.

Дифференциальные уравнения высших порядков

Дифференциальное уравнение n-го порядка имеет вид F(x, y, y', y'',..., y

(n)

)

 

этого уравнения можно выразить старшую производную

y

(n)

, то мы

 

 

 

называемое уравнение, разрешенное относительно старшей производной (

функция (n +1)

-й переменной).

 

 

 

 

 

y

(n)

= f (x, y, y',..., y

(n1)

) .

 

 

 

Определение 1. Задачей Коши для уравнения (9) называется задача

= 0

. Если из

получим так

f

некоторая

 

(9)

y

(n)

 

=

f (x, y, y',..., y

(n1)

)

 

 

 

 

y(x

 

)

= y

 

 

 

 

 

 

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y'(x0 ) = y'0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.....................................

 

 

 

(n1)

 

 

 

 

 

(n1)

 

 

 

y

(x

 

) = y

,

 

 

 

 

 

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(10)

где

x

 

, y

 

, y'

 

,..., y

(n1)

0

0

0

0

 

 

 

 

 

некоторые числа.

Теорема 1 (существования и единственности решения задачи Коши) (без

доказательства).

Пусть функция f и ее частные производные первого порядка по всем

аргументам, кроме x , непрерывны в некоторой области D

(n +1) - мерного пространства и

точка

(x

 

, y

 

, y'

 

,..., y

(n1)

) D

. Тогда задача Коши (10)

имеет единственное решение

0

0

0

0

 

 

 

 

 

(определенное в некоторой окрестности точки x0 ).

 

Определение 2. Пусть выполняются условия теоремы 1. Функция

 

 

 

 

 

 

 

 

 

 

y = φ(x, c1, c2 ,..., cn ) ,

(11)

162

где c1 ,c2

,...,cn постоянные, называется

общим решением

уравнения

(9) в некоторой

окрестности точки x0

U (x0 ) , если:

 

 

 

 

 

 

 

 

1.

При x U (x0 ) и

наборе

(c1

,c2 ,...,cn ) C , где

C некоторое множество (в

простых случаях c1 , c2 ,...,cn

, будут любыми числами) функция (11) является решением

уравнения (9).

 

 

 

 

 

 

 

 

 

 

 

 

Какие бы начальные условия y(x0 ) = y0 , y'(x0 ) = y0 ' ,…,

 

(n1)

(n1)

 

 

2.

y

 

(x0 ) = y0

,

где

 

(n1)

 

 

 

 

 

 

(c1 , c2 ,..., cn ) C

 

 

точка (x0 , y0 , y'0 ,..., y0

 

) D , мы не

задали, существует

набор

 

,

при

котором функция (11) удовлетворяет этим начальным условиям.

Определение 3. Равенство вида (x, y, c1, c2 ,..., cn ) = 0 , неявно задающее общее решение, называется общим интегралом уравнения (9).

Уравнения, допускающие понижение порядка

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

1. Уравнение не содержит явным образом искомую функцию y и, может быть, несколько ее первых производных, т.е. имеет вид

F(x, y

(k )

, y

(k +1)

,..., y

(n)

)

 

 

 

=

0

.

(12)

Сделаем в этом уравнении замену

y

(k )

 

=

z

, где

z =

z(x)

новая неизвестная функция (т.е. за

новую неизвестную функцию берется производная наименьшего порядка, входящая в это

уравнение).

Тогда

y

(k +1)

= z' ,

y

(k +2)

= z'' , …,

y

(n)

 

 

 

F(x, z, z',..., z

(nk )

) = 0

, т.е. порядок уравнения понизился.

 

=

z

(nk )

 

,

и (12) примет вид

Пример. Решить задачу Коши для дифференциального уравнения

x

Решение. Обозначая

y ' = z ,

z = z(x) , имеем

x

2

z ' = z

2

. Это уравнение

 

 

2

y '' = y '

2

; y(1)

= 1; y '(1)

=

1

 

 

 

 

 

 

 

 

 

2

первого порядка является

.

уравнением с

разделяющимися переменными:

1

=

1

+ c1

=

1 + c1 x

;

z =

x

,

z

x

 

x

1 + c x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

dz

 

2

dz

 

x

 

 

= z

;

 

 

 

 

 

 

2

 

 

 

dx

 

 

z

 

 

 

 

 

 

 

т.е.

y ' =

 

x

 

.

 

 

 

 

 

 

 

 

1 + c x

 

 

 

 

 

 

 

1

 

 

=

dx

. Интегрируя, получаем:

 

dz

=

x

2

z

2

 

 

 

 

 

 

Прежде чем интегрировать еще раз, найдем

dx

+ c

;

1

= −

1

+ c

;

 

 

 

 

x

2

1

 

 

z

 

x

1

 

 

 

 

 

 

 

 

c1

из второго начального

условия. При x =

1 из него

1

 

 

1 + c1

 

 

 

 

 

 

 

=

x + 1 1

dx + c2 =

dx

 

 

dx

+ c2 =

 

1

+ x

 

1 + x

 

 

 

=

=

1

2

x

;

1 + c

=

1

 

ln | 1 + x |

2

+

;

c 2

c1 = 1 .

. При

Значит

y ' =

x

. Отсюда y =

 

x

dx + c

 

 

 

2

 

 

 

=

 

 

1 + x

 

 

1 + x

 

 

x = 1 из

первого

начального условия

находим

постоянную

c 2

:

1 ln 2 + c 2

=

1

;

c 2

=

ln

2

. Таким образом, y = x ln | 1 + x | + ln 2 = x ln 1 + x .

2

2. Уравнение не содержит явным образом

F( y, y',

независимую переменную х, т.е. имеет вид

y'',..., y

(n)

) = 0 .

 

(13)

Сделаем в этом уравнении замену y'= z , где z = z( y) , т.е. за новую независимую переменную мы берем y, а за новую независимую функцию y'= z( y) . Покажем, что при такой замене порядок уравнения понижается на единицу:

y''= ( y')/x = zx/ = z y/ yx/ = z' z , т.е. y''= z' z ,

163

y'''= ( y'')

/

= (z' z)

/

= (z' z)

/

y

/

= (z'' z + z'

2

)z = z'' z

2

+ z'

2

z

x

x

y

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

Пример. Решить дифференциальное уравнение

2 yy '' = y '

2

+1 .

 

 

 

 

Решение. Обозначая y ' = z , z = z( y) , y '' = z ' z , имеем

2 yz ' z = z

2

+1

. Это уравнение первого порядка, опять-таки,

 

является уравнением с разделяющимися переменными:

2 yz

dz

= z

2

 

 

 

dy

 

 

+ 1

;

 

2z

dz =

dy

 

2

 

 

z

+ 1

 

y

 

 

. Интегрируя, получаем:

2z

dz z 2 + 1

Еще раз

=

 

dy

 

;

 

1

 

 

+ c

 

 

 

y

 

 

разделяем

 

d ( z

2

+ 1)

 

 

 

=

z

2

 

+ 1

 

 

 

 

переменные:

ln | y

dy dx

| +c

; ln( z

2

 

1

 

 

=

c y

 

1

 

+

1

1)

;

= ln | c y | 1

dy

c y 1 1

;

=

z

2

+ 1

dx .

= c y ; 1

Далее

z =

c

 

1

имеем:

y 1

1

 

c

 

1

 

, т.е.

d (c y 1

c y 1

y '

=

c y 1

 

 

1

1)

= dx + c2

1

 

 

.

;

2

c 1

c y 1 1

=

( x + c

)

2

 

. Возводим в квадрат:

4

c

2

 

1

(c y 1) 1

= ( x + c

)

 

2

2

 

;

 

c

2

c y 1 =

 

1

1

4

 

( x + c

)

2

 

 

2

 

;

y =

c

( x + c

)

2

+

1

1

 

 

 

 

 

 

 

4

 

2

 

 

c

 

 

 

 

 

 

 

 

 

 

 

1

.

164

ЛЕКЦИЯ 11

ЛИНЕЙНЫЕ ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ ВЫСШИХ ПОРЯДКОВ

Определение 1. Линейным дифференциальным оператором называется оператор вида

 

L( y) = y

(n)

+ p y

(n1)

+ p

 

y

(n2)

+... + p

 

y

,

(1)

 

 

 

2

 

n

 

 

 

1

 

 

 

 

 

 

где pi

= pi (x) , i =1,2,...,n определены и непрерывны на некотором интервале (a,b) .

 

Линейный дифференциальный оператор ставит в соответствие функции

функцию

L( y) , определенную по формуле (1).

Свойства линейного дифференциального оператора: 1. Ly) = αL( y) , где α постоянная.

y = y(x)

новую

L(αy) = (αy)(n) + p (αy)(n1) +... + p

(αy) = α( y(n) + p y(n1)

+... + p

n

y) =

αL( y) . ■

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

2.

L( y1 + y2 ) = L( y1 ) + L( y2 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L( y + y

2

) =

( y

+ y

2

)(n) + p ( y

+ y

2

)(n1)

+... + p

n

( y + y

2

) =

 

 

 

 

 

1

 

 

 

1

 

 

1 1

 

 

 

 

 

 

1

 

 

 

 

y

(n)

+ p

y

(n1)

+... + p

 

y + y

 

(n)

+ p y

 

(n1)

+... + p

 

y

2 = L( y1 ) + L( y2 ) . ■

 

 

 

 

 

2

 

2

n

 

= 1

 

 

1

1

 

 

 

 

n 1

 

 

1

 

 

 

 

 

 

Определение 2. Линейным однородным дифференциальным уравнением называется

уравнение вида L( y) = 0

, или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

(n)

+ p (x) y

(n1)

+ p

 

(x) y

(n2)

+...

+ p

 

(x) y = 0

.

 

 

 

(2)

 

 

 

 

 

2

 

 

n

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

Свойства решений линейного однородного дифференциального уравнения.

 

 

 

 

1. y = y(x) решение (2),

c

произвольная постоянная cy тоже решение (2).

▲Согласно свойству 1 линейного дифференциального оператора, L(cy) = cL( y) = 0

. ■

 

 

2. y1 = y1 (x)

и

y2

= y2 (x) решения (2)

y1

+ y2

тоже решение (2).

 

 

 

Согласно

свойству

2

 

линейного

дифференциального

оператора,

L( y1 + y2 ) = L( y1 ) + L( y2 ) = 0 + 0 = 0 . ■

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3. y1 , y2

,...yn решения

(2),

 

c1 ,c2

,...,cn произвольные

постоянные

 

 

y = c1 y1 + c2 y2

+... + cn yn тоже решение (2).

 

 

 

 

 

 

 

 

 

 

 

▲ Это свойство является очевидным следствием свойств 1 и 2. ■

 

 

 

 

Функция y = c1 y1 + c2 y2 +... + cn yn зависит от

x

и от n

произвольных постоянных.

При

любых значениях этих постоянных она является решением уравнения (2). Ниже будут изучаться условия, при которых эта функция является общим решением уравнения (2). т.е. условия, при которых она удовлетворяет определению 2 предыдущей лекции.

Линейная зависимость и независимость функций

Определение 3. Функции y1 (x), y2 (x),...yk (x) называются линейно зависимыми на

некотором множестве M , если существуют постоянные

c

, c

2

,...,c

k

1

 

 

, хотя бы одна из

которых отлична от нуля, такие, что

c1 y1 (x) + c2 y2 (x) +... + ck yk (x) 0

(здесь знак тождества « » означает выполнение равенства (3) для y1 (x), y2 (x),...yk (x) называются линейно независимыми на множестве

(3) выполняется только при всех коэффициентах ci = 0 .

(3)

x M ). Функции M , если тождество

165

Пример 1. Покажем, что функции

 

y0

=1

,

y1

 

= x

,

y2

= x

2

, …,

 

 

 

yk

= x

k

 

линейно независимы

на

любом

 

 

 

 

 

 

 

 

 

 

 

конечном или бесконечном промежутке M .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

▲ Из равенства

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c0 + c1x + c2 x

2

 

+... + ck x

k

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

следует, что любой

x M

является корнем уравнения (4). Но любой многочлен степени

k имеет только k

корней, откуда следует, что единственным возможным является случай c

 

 

= c

 

= c

 

= ... = c

 

= 0

. ■

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

1

 

 

2

 

 

 

 

k

 

 

 

 

Пример 2. Покажем,

что

 

функции y

 

= e

k x

 

,

 

y

 

 

= e

k

 

x

,

 

 

 

…, y

 

= e

k

 

x

, k

 

 

k

 

 

 

при

i j ,

линейно

 

 

 

1

 

 

 

 

 

2

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

независимы на любом конечном или бесконечном промежутке M .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

▲ Пусть

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c e

k x

+ c

 

e

k

 

x

+ c

e

k

 

x

+... + c

 

 

e

k

 

x

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

2

 

 

3

 

 

n

 

n

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5)

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Докажем, что в равенстве (5) все коэффициенты

 

ci

 

обязательно равны 0, например,

докажем,

что

cn = 0 .

k x

:

c + c

 

 

e

(k

 

k ) x

+ c

e

(k

 

k ) x

 

 

+

... + c

 

 

e

(k

k

 

) x

0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Разделим (5) на e 1

2

 

 

2

 

1

 

3

 

 

1

 

 

 

 

 

n

 

 

n

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Дифференцируем по x

:

c

 

 

(k

 

 

k )e

( k

 

k ) x

+ c

 

(k

 

k )e

( k

 

k

) x

 

+ ...

+ c

 

 

 

(k

 

 

k )e

( k

k ) x

0

.

 

 

 

 

 

2

 

 

2

 

1

 

 

3

 

 

 

3

 

1

 

 

 

n

 

n

 

 

n

 

1

 

 

 

 

 

2

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

Разделим последнее равенство на

 

e

(k2 k1 ) x

, учитывая, что

 

k k

(k

 

k ) = k k

 

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

(k

 

 

k ) + c

(k

 

 

k )e

( k

 

k

) x

+ ... + c

 

(k

 

 

 

k )e

( k

k

 

) x

 

0

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

2

 

 

 

 

 

 

 

n

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

 

 

 

1

 

 

3

 

 

 

3

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

n

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Дифференцируем по x

и т.д. В итоге получим:

c

 

(k

 

k )(k

 

 

k

 

 

) ... (k

 

 

k

 

 

 

)e

( k

k

n1

) x

0

 

.

 

 

 

 

n

n

2

n

n1

 

n

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В этом равенстве все сомножители, кроме первого, отличны от 0, отсюда cn

 

= 0 . ■

 

 

 

 

 

 

 

 

 

 

 

Линейно зависимые и линейно независимые функции обладают обычными свойствами линейно зависимых и линейно независимых векторов.

Кроме того, оказывается, что линейная зависимость и независимость функций тесно связана с так называемым определителем Вронского.

Определение 4. Пусть даны k 1 раз дифференцируемые функции

y1 (x)

,…,

y

k

(x)

 

 

.

Определителем Вронского этих функций называется следующий определитель порядка:

k

-го

 

 

y

 

y

2

...

 

y

k

 

 

 

 

 

 

1

 

 

 

 

 

 

 

y'

y'

2

...

y'

k

 

 

1

 

 

 

 

 

 

 

W = W (x) =

y''

y''

2

...

y''

k

 

 

1

 

 

 

 

 

 

 

 

 

 

 

...

 

...

...

 

...

 

y

(k 1)

y

(k 1)

...

y

(k 1)

 

 

2

 

 

 

k

 

 

 

 

1

 

 

 

 

 

 

 

 

 

(все функции в этом определителе берутся в некоторой точке x ).

Теорема 1 (первая теорема об определителе Вронского). Пусть функции

y1 (x)

yk (x) линейно зависимы на некотором множестве M . Тогда на этом определитель Вронского тождественно равен 0: W (x) 0, х M.

, y2 (x) , …,

множестве

▲ По условию теоремы существуют постоянные c1 , c2 ,...,ck , хотя отлична от 0, такие, что c1 y + c2 y2 +... + ck yk 0 . Пусть, например, ck

бы

0

одна из которых

. Тогда, разделив

это равенство на

c

 

и обозначив

ci

= α

, его можно переписать в виде

k

 

 

 

 

 

i

 

 

 

 

 

ck

 

y

= α y + α

y +... + α

y

. Отсюда

y(i) = α y(i) + α

2

y(i) +... + α

k 1

y(i)

, i =1,2,...,k 1, и

k

1 1

2 2

k 1 k 1

 

k

1 1

2

k 1

 

166

 

y

 

 

y

2

 

...

y

k 1

 

α y

 

+ α

2

y

2

 

+... + α

k 1

y

k 1

 

 

 

 

1

 

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

'

y

2

'

...

y

k 1

'

α

y

 

'+ α

2

y

2

 

'+... + α

k 1

y

k 1

'

 

 

W (x) =

1

 

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

...

 

...

 

... ...

 

.....................................

 

 

 

 

 

 

 

 

 

 

 

 

 

y

(k 1)

y

(k 1)

...

y

(k 1)

α y

(k 1)

+ α

 

 

y

(k 1)

+...

+ α

 

 

y

( k 1)

 

 

 

2

 

k 1

 

 

 

2

2

k

1

k 1

 

1

 

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

,

так как последний столбец этого определителя является линейной комбинацией остальных его столбцов. ■

Примеры показывают, что теорема, обратная данной, не верна.

Теорема 2 (вторая теорема об определителе Вронского). Пусть

n

решений линейного

однородного дифференциального уравнения n -го порядка L( y) = 0

 

y1 (x), y2 (x), …, yn (x)

линейно независимы на интервале (a,b) . Тогда определитель Вронского этих функций не обращается в 0 ни в одной точке интервала (a,b) : W (x) 0 , x (a,b) .

Замечание: Первая и вторая теоремы относятся к разным объектам: первая – к произвольным функциям, а вторая – к решениям линейных однородных дифференциальных уравнений; для последних определитель Вронского либо тождественно равен 0, либо он не обращается в 0 ни в одной точке.

▲ Пусть

x0 (a,b) :

W (x0 ) = 0

. Рассмотрим функцию

y(x) = c1 y1(x) +... + cn yn (x) .

Согласно свойству 3. решений линейного однородного дифференциального уравнения, при

любых значениях постоянных ci , i =1, 2,

..., n функция y(x) является решением уравнения

L( y) = 0. Теперь подберем постоянные ci

так, чтобы эта функция в точке

x0

удовлетворяла

нулевым начальным условиям y(x0 ) = 0 ,

y'(x0 ) = 0 , y''(x0 ) = 0 …

y

(n1)

(x

 

)

= 0

. Имеем:

 

0

 

 

 

 

 

y(x

0

) = c y

(x

0

) + c

2

y

2

(x

0

) +

...

 

+ c

n

y

n

(x

0

) = 0

 

 

 

 

 

 

 

 

 

 

 

 

1 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y'(x

 

 

 

) = c y' (x

 

) + c

 

y'

 

(x

 

) +

 

+ c

 

 

y'

 

 

(x

 

) = 0

 

 

 

 

0

0

2

2

0

...

n

n

0

 

 

 

 

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y''(x

 

 

) = c y''

(x

 

) + c

 

y''

 

 

(x

 

) +

 

 

+ c

 

y''

 

(x

 

) = 0

 

 

 

0

0

2

2

0

 

 

n

n

0

 

 

 

 

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

........................................................................

 

 

 

 

 

 

 

 

(n1)

 

 

 

 

 

 

(n1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(n1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(n1)

 

 

 

y

(x

 

) = c y

(x

 

) + c

 

 

y

(x

 

) +

 

 

 

+ c

 

y

(x

 

) = 0

 

 

 

 

 

0

 

 

 

 

 

0

2

2

 

 

0

...

 

 

n

n

0

 

 

 

 

 

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6)

(6) является системой

n линейных однородных уравнений с

n

неизвестными

c1 , c2 ,...,cn .

Определитель этой системы – это определитель, составленный из коэффициентов при

неизвестных

ci

, т.е. определитель Вронского W (x0 ) = 0

. Но тогда система (6), как система

n линейных однородных уравнений с n неизвестными и определителем, равным 0, имеет ненулевые решения.

Пусть

c

*

, c

*

,...,c

*

одно

из

таких решений.

Тогда функция

y

= c

*

y

+ c

*

y

 

+... + c

*

y

 

 

2

n

 

2

2

n

n

1

 

 

 

 

 

 

1

1

 

 

 

 

удовлетворяет на интервале

(a,b)

уравнению

L( y) = 0 и нулевым начальным условиям

y(x ) = 0 ,

y'(x

0

) = 0 , …,

y (n1) (x

0

) = 0 . Но этому же уравнению и этим же начальным

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

условиям удовлетворяет и функция y 0 . В силу единственности решения задачи Коши, эти два решения должны совпадать:

y = c

*

y

+ c

*

y

 

+... + c

*

y

 

0

,

x (a,b) .

(7)

 

2

2

n

n

1

1

 

 

 

 

 

Так как в тождестве (7) хотя бы один из коэффициентов ci* отличен от 0, то это тождество означает линейную зависимость на интервале (a,b) функций y1 , y2 ,..., yn , что противоречит условию теоремы. Значит, наше предположение W (x0 ) = 0 не верно. ■

167

Структура общего решения линейного однородного уравнения

Определение 5.

Фундаментальной системой решений линейного однородного

дифференциального

уравнения n -го порядка

L( y) = 0

называется любая система

n

линейно независимых решений этого уравнения.

 

 

 

Теорема

3 (структура общего решения линейного однородного уравнения). Пусть

функции

y1, y2 ,..., yn образуют фундаментальную систему решений линейного однородного

дифференциального уравнения n -го порядка

L( y) = 0 . Тогда функция

 

y = c1 y1 + c2 y2

+... + cn yn ,

(8)

где c1 ,c2 ,...,cn произвольные постоянные, является общим решением этого уравнения. ▲ Проверим, что (8) удовлетворяет определению общего решения. При любых значениях

постоянных ci , i = 0,

1, 2, ..., n

эта функция,

согласно свойству 3 решений линейного

однородного дифференциального уравнения, является решением уравнения

L( y) = 0 .

Теперь проверим второе условие определения 7: зададим в точке x0

(a,b)

произвольные

начальные условия y(x0 ) = y0 , y'(x0 ) = y'0 , …,

 

(n1)

 

 

(n1)

и подберем постоянные ci

y

 

(x0 ) = y0

так, чтобы функция (8) удовлетворяла этим начальным условиям. Имеем:

 

y(x0 ) = c1 y1 (x0 ) + c2 y2 (x0 ) +... + cn yn (x0 ) = y0

 

 

 

 

 

 

 

 

 

 

 

 

'

 

 

 

y '(x0 ) = c1 y '1 (x0 ) + c2 y '2 (x0 ) +

... + cn y 'n (x0 ) = y0

 

 

 

 

 

 

 

 

 

 

 

 

 

"

.

(9)

y ''(x0 ) = c1 y ''1 (x0 ) + c2 y ''2 (x0 ) +... + cn y ''n (x0 ) = y0

........................................................................

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(n1) (x

) = c y(n1) (x

) + c y(n1) (x )

+... + c

n

y(n1)

(x )

= y(n1)

 

 

0

1 1

0

2 2

 

0

 

n

0

0

 

(9) является системой

n линейных уравнений с

n неизвестными c1 , c2 ,...,cn . Определитель

этой системы – это определитель, составленный из коэффициентов при неизвестных ci , т.е. определитель Вронского W (x0 ) , который отличен от 0 в силу линейной независимости

функций

y

,

1

 

y

2

,...,

 

 

yn

. Но тогда система (9), как система

n

линейных уравнений с n

неизвестными и определителем, не равным 0, имеет единственное решение. ■

Таким образом, для нахождения общего решения уравнения (1) нужно знать фундаментальную систему его решений. Однако в общем случае методов нахождения такой фундаментальной системы не существует.

Ниже будет описан способ нахождения фундаментальной системы решений для одного класса уравнений вида L( y) = 0.

168

ЛЕКЦИЯ 12

ЛИНЕЙНЫЕ ОДНОРОДНЫЕ УРАВНЕНИЯ С ПОСТОЯННЫМИ КОЭФФИЦИЕНТАМИ

В этом параграфе будет рассматриваться оператор вида

 

 

 

L( y) = a y

(n)

+ a y

(n1)

+

... + a

 

y'+a

 

y

,

 

 

 

 

 

 

 

 

 

 

 

(1)

 

 

 

 

 

n1

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x R

a0

0 , a1, a2 ,..., an постоянные действительные числа, и уравнение вида

L( y) = 0 :

 

 

 

a0 y

(n)

+ a1 y

(n1)

+... + an1 y'+an y = 0 .

 

 

 

 

 

 

 

 

 

 

 

 

(2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Уравнение (2) можно разделить на

a0 , поэтому к нему

применимы

 

предыдущие

рассуждения, и общее решение (2) ищется по формуле:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y = c1 y1

+ c2 y2

+... + cn yn ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3)

где

c1, c2 ,...,cn произвольные

постоянные,

а

y1, y2

,..., yn

фундаментальная

система

решений.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для нахождения последней будем искать решения (2) в виде

y = e

kx

, где

k некоторое

 

 

 

число.

Подставляя

y = e

kx

 

в

(2)

 

и

учитывая, что

(e

kx

)

(m)

= k

m

e

kx

,

имеем:

 

 

 

 

 

 

 

 

 

a0k nekx + a1k n1ekx +... + an1kekx + anekx

= 0 , или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т.е.

y = ekx

a k

n

+ a k

n1

+... + a

 

k + a

 

=

 

 

n1

n

0

 

1

 

 

 

 

решение (2) тогда и только тогда, когда число

0

k

.

корень (4).

(4)

Определение 1. Уравнение (4) называется характеристическим уравнением для дифференциального уравнения (2).

Характеристическое уравнение (4) получается из дифференциального уравнения (2)

заменой производной

y

i

 

на

k

i

 

,

i = 0, 1, 2, ..., n

(под нулевой производной функции

понимается сама эта функция).

Определение 2. Левую часть характеристического уравнения (4)

характеристическим многочленом и обозначим Ф(k) .

Уравнение (4), как и всякое алгебраическое уравнение степени

n , имеет ровно n

учетом их кратности. Рассмотрим следующие 4 случая:

 

1. Все корни характеристического уравнения действительные и различные

назовем

корней с

Пусть k1, k2 ,...,kn эти корни. Этим корням соответствуют

n

решений уравнения (2):

y1

=

ek1x

,

y2

=

k

x

e

 

2

 

, …, yn = ek n x . Эти функции линейно независимы (см. пример выше),

значит, они образуют фундаментальную систему решений уравнения (2), и общее решение этого уравнения задается формулой (3).

 

 

Пример.

Решить дифференциальное уравнение y''' 2 y"

8y'

 

= 0 .

 

 

 

 

 

Решение.

Характеристическое

уравнение

имеет

вид:

k

3

2k

2

8k = 0 .

Решаем

это

уравнение:

 

 

k (k

2

2k 8) = 0; k

= 0, k

= 4, k

= −2.

Этим трем действительным различным корням соответствуют три решения из

 

 

 

 

 

1

2

3

 

фундаментальной

системы

решений:

y

= 1, y

= e4 x , y

= e2 x ,

и общее решение исходного уравнения

имеет вид:

 

 

 

 

 

 

 

 

1

2

3

 

 

 

 

 

 

 

 

 

y = c 1 + c e

4 x

+ c e

2 x

 

 

1

2

 

3

 

= c

+ c e

 

4 x

1

2

+ c e

2 x

 

3

 

.

2. Корни характеристического уравнения разные, но среди них есть комплексные

Так как в предыдущих рассуждениях нигде не использовалось, что корни

169