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

высшая математика

.pdf
Скачиваний:
73
Добавлен:
18.02.2016
Размер:
5.88 Mб
Скачать

Лекция 50

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

Рассмотрены основные свойства решений однородных и неоднородных линейных дифференциальных уравнений высших порядков.

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

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

 

y(n) + p

n–1

(x)y(n–1) + . . . + p

(x)y(1)

+ p

(x)y = f(x),

(1)

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

0

 

 

 

 

где функции pi(x), i =

 

 

 

 

 

 

 

0, n 1 и f (x) заданы и непрерывны на некотором

интервале I = (a; b).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для краткости записи уравнения (1) введем понятие линейного

дифференциального оператора (ЛДО) n-го порядка.

 

ЛДО n-го порядка называют выражение:

 

 

def

d n

 

 

 

 

 

 

 

d n

1

 

 

 

 

 

 

d

 

 

Ln =

 

 

 

 

 

+

 

pn1(x )

 

 

 

+

... + p1

(x )

 

 

 

+ p0 (x ) .

(2)

dx n

 

 

 

dx n

1

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

i

def

y (i)

 

 

 

 

 

 

 

 

 

 

 

 

Считаем, что

 

 

 

[y] =

 

,

 

i =

 

1, n

.

 

 

 

 

 

 

 

 

dx i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда ЛДУ (1) с использованием обозначения (2) примет вид

 

 

 

 

 

 

 

 

 

 

Ln[y] = f (x).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1)

Частным решением ЛДУ (1) называется такое решение этого

уравнения, которое удовлетворяет начальным условиям:

 

y(x

) = y

0

, y(1)(x

) = y (1), ..., y(n–1)(x

 

) = y (n–1),

(3)

0

 

 

 

 

 

 

0

 

0

 

 

 

 

 

 

 

0

 

0

 

где заданное x0 I, a y0 , y0(1), ..., y0(n–1) – любые заданные числа.

 

Задачей Коши (или начальной задачей) называют решение ЛДУ

(1) при начальных условиях (3).

 

Нетрудно убедиться, что ЛДО удовлетворяет условиям:

 

а) однородности: Ln[сy] = с Ln[y];

 

б) аддитивности: Ln[y1 + y2] = Ln[y1] + Ln[y2].

если

Уравнение (1) называется линейным неоднородным,

f (x) 0

и линейным однородным, соответствующим (1),

если

f (x) 0,

т.е. линейное однородное дифференциальное уравнение

(ЛОДУ) имеет вид

 

 

Ln[y] = 0.

(4)

3 11

Общим

решением ЛДУ (1) называется такая функция

y = ϕ (x, C1, ...,

Cn), зависящая от n произвольных постоянных, что

выполняются условия: 1) при любых постоянных C1, ..., Cn эта функ-

ция является решением ЛДУ (1); 2) каковы бы ни были начальные

условия (3), можно подобрать такие значения C10, ..., Cn0 постоянных

C

, ..., C ,

что решение y = ϕ (x, C

0, ..., C 0) ЛДУ (1) будет удовлетво-

1

n

1

n

 

 

 

 

 

 

рять заданным начальным условиям (3).

 

 

 

 

 

 

 

20. Свойства решений ЛОДУ (4). Пусть y

,

y

,

..., y

n

 

 

 

1

 

2

 

 

 

 

 

 

n

 

 

 

 

 

 

решения уравнения (4). Тогда y =

C k yk – также решение уравнения

 

 

 

k = 1

 

 

 

 

 

 

(4) при любых произвольных постоянных Ci ,

i =

1, n

.

 

Действительно, по условию,

Ln[yk] = 0.

Имеем:

 

 

n

 

=

n

n

 

 

n

= 0 .

Ln [y] = Ln

C k yk

Ln [C k yk ]= C k Ln [yk ] = Ck 0

k = 1

 

 

k = 1

k = 1

 

k = 1

 

Замечание 1. Систему функций y1(x), y2(x), ..., yn(x) называют линейно

зависимой (ЛЗ) на I, если существуют постоянные числа

α

,

α

, ... ,

α

 

, хотя

бы одно из которых отлично от нуля,

что выполняется:

1

 

2

 

n

 

α 1 y1(x) + α 2 y2(x) + ... + α

n yn(x) = 0 x (a; b).

 

 

 

 

 

(5)

Если равенство (5) имеет место только при α

1 = α

2

= ... = α

n

= 0, то

система функций y1(x), y2(x), ..., yn(x)

называется ЛНЗ (линейно независи-

мой) на I.

 

 

 

 

 

 

 

 

 

Пример 1. Функции 1, х, х2 ЛНЗ на любом интервале I. Дей-

ствительно, равенство (5) в данном случае имеет вид

 

 

 

 

 

α 1 + α 2x + α 3x2 = 0

 

x (a; b), что

 

 

 

 

 

 

 

эквивалентно α 1 = α 2 = α 3 = 0, т.к. многочлен обращается в тождественный нуль на любом интервале I тогда и только тогда, когда все его коэффициенты равны нулю.

 

Структура общего решения ЛОДУ (4) следующая: 1) уравнение

(4)

имеет ровно n ЛНЗ решений y1(x), y2(x),

..., yn(x), x I = (a; b);

2) общее решение этого уравнения имеет вид :

 

 

 

n

 

y(x ) = C1y1(x ) + C2 y2 (x ) + ... + Cn yn (x ) =

Ci yi (x ) ,

где

C1, ..., Cn – произвольные постоянные.

i= 1

 

Фундаментальной системой решений уравнения (4) называют систему n ЛНЗ решений этого уравнения.

Таким образом, уравнение (4) имеет фундаментальную систему решений.

3 12

Рассмотрим систему n функций y1(x), y2(x), ..., yn(x), непрерывных вместе со своими производными до (n–1)-го включительно порядка на I. Определителем Вронского, или вронскианом, W [y1, y2, ..., yn] этой системы функций называется определитель из функций и их производных, составленный следующим образом:

 

 

 

 

 

y1

y2

%

yn

 

 

 

 

 

 

 

y

y

%

y

 

 

 

W [y1, ..., yn ] =

 

1

2

 

n

.

 

 

 

&

&

%

&

 

 

 

 

 

 

 

 

 

 

 

 

y1( n1)

y2(n1)

% yn(n1)

 

 

Если указанные функции y1,

y2, ...,

yn

ЛЗ на I, то определитель

Вронского

W [y1, y2, ...,

yn]

равен тождественно нулю на I,

т.е. если

найдется

хотя

бы

одно

значение

x (a; b), при

котором

W [y1, y2, ..., yn]

0, то эти функции y1, y2,

..., yn ЛНЗ на I.

 

Пример 2. Является ли ЛНЗ система следующих функций: {x, cos x, sin x}.

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

 

 

 

 

 

 

 

 

 

 

 

x

cos x

sin x

 

 

 

x

cos x

sin x

 

 

 

 

 

 

W [x, cos x,sin x ] =

1

sin x

cos x

 

=

 

1

sin x

cos x

 

=

 

0

cos x

sin x

 

 

 

x

0

0

 

 

=

x (cos2 x +

sin2 x ) =

 

x

0

 

 

 

на любом интервале I.

Таким образом, эта система ЛНЗ на любом интервале I. Для уравнения (4) имеет место формула Остроградского-Лиу-

вилля:

 

 

x

 

 

 

 

pn1 (t )dt

 

 

W [y1, y2 , ..., yn ] = W [y1, y2 , ..., yn ]

 

x = x 0 e x0

.

(6)

 

 

 

Сформулированные выше свойства позволяют решить следующую задачу: по данной системе n ЛНЗ функций y1(x), y2(x), ..., yn(x), непрерывных вместе со своими производными до n-го порядка включительно на I = (а; b), таких, что W [y1, y2, ..., yn] 0 на I, построить ЛОДУ n-го порядка, решениями которого являются данные функции. Искомым дифференциальным уравнением является следующее равенство:

3 13

 

 

y

y1

y2

%

yn

 

W [y, y1

, ..., yn ] =

y

y

y

%

y

= 0 .

%

1

2

%

n

 

 

%

%

%

 

 

 

y (n)

y1( n)

y2(n)

% yn( n)

 

Пример 3. Получить ЛОДУ 3-го порядка, решениями которого являются функции х, sin x, cos x.

Решение. Составим определитель Вронского и равенство:

 

y

x

 

sin x

 

cos x

 

W [y, x, sin x, cos x ] =

y

1

 

cos x

sin x

= 0 .

y′′

0

sin x

cos x

 

y′′′

0

cos x

 

sin x

 

После вычисления определителя получаем уравнение:

 

 

y +

xy′ −

y′′ +

xy′′′ =

0 .

Ответ: y′′′ −

1

y′′ +

y′ −

1

y =

0 ,

x I

(0 I ).

x

x

30. Линейные неоднородные дифференциальные уравнения (1). Пусть y0(x) – частное решение ЛДУ (1), а {y1(x), y2(x), ... , yn(x)} образуют фундаментальную систему решений соответствующего уравнению (1) однородного уравнения (4), тогда

 

n

Ci yi (x ) , где С , ...

 

 

 

y(x ) =

y0 (x ) +

, C

n

– произвольные постоянные,

 

i = 1

1

 

 

есть общее решение ЛДУ (1).

Таким образом, общее решение линейного неоднородного дифференциального уравнения (1) есть сумма его частного решения y0(x) и общего решения соответствующего ему однородного уравнения (4)

n

y = Ci yi (x ) .

i= 1

Для отыскания частного решения линейного неоднородного уравнения (1) применяется метод вариации произвольных постоянных, который заключается в следующем. Пусть известна фундаментальная система {y1, y2, ... , yn} решений соответствующего однородного уравнения (4). Тогда общее решение неоднородного уравнения (1) будем искать в виде общего решения ЛОУ:

3 14

n

 

y = C k (x )yk ,

(7)

k = 1

но с переменными произвольными постоянными. Для определения Ck(x),

k =

1, n

, которых

n, а уравнение

одно,

требуется наложить

(n – 1) условие при вычислении производных функции у. Имеем:

 

 

y

 

n

C(x) y

 

 

n

 

 

(x) y

 

 

 

=

k

+

C

k

,

 

 

 

 

k

 

 

 

k

 

 

 

 

k = 1

 

 

 

k =

1

 

 

 

 

 

 

 

$!#!"

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

y′′

 

n

C(x) y

 

n

 

 

(x) y′′

 

 

 

=

+

C

k

,

 

 

 

 

k

 

k

 

 

k

 

 

 

 

k =

1

 

 

 

k =

1

 

 

 

 

 

 

 

$!#!"

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

. . . . . .

 

. . . . ,

 

 

 

 

 

$!#!"

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

n

 

 

 

 

 

y (n) = C k(x )yk( n

1) + Ck (x )yk( n) .

 

 

 

 

k = 1

 

 

 

 

k = 1

 

 

 

 

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

уравнение для нахождения C k(x ) :

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

Ln [y]

= Ln

[y ]+

Ck(x )yk(n1) = f (x ) .

 

 

 

 

$#"

k = 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

Таким образом, получаем следующую квадратную систему линейных алгебраических уравнений для функций C k(x ) :

C1y1 +

C2y2

+

... +

Cnyn

=

0,

 

C y′ +

C y

+

... +

C y

=

0,

 

 

1 1

 

2 2

 

 

n n

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

(n1)

 

 

(n1)

+ ... +

 

(n1)

= f (x ).

C1y1

 

+ C2y2

Cnyn

Решение системы (8) можно найти по формулам Крамера:

C k

(x ) =

k

 

 

 

 

 

,

k = 1, n ,

W [y, y1

, ..., yn ]

 

 

 

 

 

 

(8)

(9)

причем определителем этой системы является определитель Вронского.

3 15

Интегрируя (9), получаем функции Ck(x), k = 1, n .

Пример 4. В примерах 2 и 3 показано, что {x, cos x, sin x} – фундаментальная система решений ОДУ для уравнения:

y′′′ −

1

y′′ + y′ −

1

y = x .

x

x

Найти частное и общее решения этого уравнения.

Решение. y = C1(x) x + C2(x) cos x + C3(x) sin x.

Вычислим определители:

 

 

 

 

 

 

 

 

 

 

x

 

 

cos x

sin x

 

 

 

 

 

W [x, cos x, sin x ] =

 

1

 

sin x

cos x

=

x ,

 

 

 

 

 

 

 

 

 

 

0

 

cos x

sin x

 

 

 

0

 

cos x

 

sin x

 

= x ,

 

 

 

 

x

0

sin x

 

 

 

x(x cos x sin x ) ,

 

 

 

 

 

 

 

 

 

 

1 =

0

sin x

 

cos x

 

2

=

 

1

0

cos x

 

 

=

 

x

cos x

 

sin x

 

 

 

 

 

0

x

sin x

 

 

 

 

 

 

 

 

 

 

x

 

cos x

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

=

 

1

sin x

0

=

x(x sin x cos x ) .

 

 

 

 

 

 

0

cos x

 

x

 

 

 

 

 

 

 

 

 

Вычислим Ck(x), k = 1, 2, 3.

 

 

 

C1′ = 1 C1(x ) = x + C1 ,

 

 

C2′ =

x cos x +

sin x

 

C2 (x ) = − x cos x dx

cos x + C2

=

 

 

 

= − x sin x

2 cos x + C2 ,

 

 

C3′ =

x sin x

cos x

 

C3 (x ) = − x sin x dx

sin x + C3

=

 

 

 

=

x cos x

 

2 sin x + C3 .

 

 

Окончательно получим:

 

 

 

 

 

 

y = (x + C1 )x + (

x sin x

2 cos x + C2 ) cos x + (x cos x 2 sin x + C3 ) sin x =

 

=

x 2

2 +

C

1

x + C

2

cos x + C

3

sin x .

 

 

 

 

$#"

 

 

 

 

 

 

 

 

 

y0

$!!!!#!!!!"

 

 

 

 

 

 

 

 

 

y

 

 

 

 

3 16

! Задания для самостоятельной работы

 

 

 

 

 

 

 

2

 

 

1. Проинтегрировать уравнение

y′′ +

 

 

y′ +

y =

0 , имеющее ча-

x

 

 

 

sin x

 

 

 

 

 

 

стное решение y

 

=

.

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Указание. Сделать замену y = y1 z dz , где z – новая неизвестная функция.

2. Выписать решение системы (8) для уравнения (1) второго порядка:

y′′ + p1(x )y′ + p0 (x )y =

f (x ) .

 

 

 

 

3. Дано уравнение y′′′ −

y′ = 0 .

Составляют ли фундаменталь-

 

x

–x

 

 

ex

+ e

x

ную систему решений функции e , e

 

,

ch x =

 

 

 

, являющиеся,

 

 

2

 

 

 

 

 

 

 

 

 

как легко проверить, решениями этого уравнения.

 

 

4. Установить, будут ли ЛНЗ в промежутке своего существова-

ния следующие функции:

 

 

 

 

 

 

 

 

а) 2x2 + 1, x2 – 1, x + 2; б)

x , x + a , x + 2a ;

 

 

в) ln (2x), ln (3x), ln (4x).

 

 

 

 

 

 

 

 

3 17

Лекция 51

Линейные однородные дифференциальные уравнения n-го порядка с постоянными коэффициентами

Приведены формулы общего решения ЛОДУ n-го порядка как для случая простых, так и для случая кратных корней характеристического уравнения, рассмотрены линейные однородные уравнения Эйлера.

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

 

y(n) + p

n–1

y(n–1) + ... + p

1

y(1) + p

y = 0,

(1)

 

 

0

 

 

где pi , i = 0, n 1 – вещественные постоянные.

Для нахождения общего решения уравнения (1) поступим так. Составляем так называемое характеристическое уравнение для уравнения (1) заменой производных y( j ) на степени λj:

λn + p

n–1

λn–1 + ... + p

λ + p

0

= 0.

(2)

 

1

 

 

 

Пусть λ1, λ2, ... , λn – корни уравнения (2), причем среди них могут быть и кратные.

Возможны следующие случаи.

10. λ1, λ2, ..., λn – вещественные и различные. Тогда фундаментальная система решений уравнения (1) имеет вид

eλ1x , eλ 2 x , ..., eλ n x ,

(3)

и общим решением уравнения (1) будет

 

y = C1eλ1x + C2 eλ2 x + ... + Cneλ n x ,

 

где C1, C2, ... , Cn – произвольные постоянные.

 

Пример 1. Найти общее решение уравнения:

 

y′′′ − 2y′′ −

3y′ = 0 .

 

Решение. Составляем

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

уравнение

λ3 – 2λ2 – 3λ = 0. Находим его корни: λ1 = 0, λ2 = –1, λ3 = 3. Так как они

действительные и различные, то общее решение имеет вид: y = C1 + C2 ex + C3e3x .

3 18

20. Корни характеристического уравнения вещественные, но среди них есть кратные. Пусть, например, λ1 = λ2 = ... = λk, т.е. λ1 является k-крат- ным корнем уравнения (2), а все остальные (n – k) корней различные. Фундаментальная система решений уравнения (1) в этом случае имеет вид

eλ1x , xe λ1x , ..., x k 1eλ1x , eλ k + 1x , ..., eλ n x ,

(4)

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

y = C1eλ1x +

C2 xe λ1x +

... + Ck x k 1eλ1x + Ck + 1eλ k + 1x + ... + Cn eλ n x .

Пример 2.

Найти общее решение уравнения

 

 

y′′′ + 2y′′ + y′ = 0 .

Решение.

Характеристическое уравнение имеет вид

 

 

λ3 + 2λ2 + λ = 0.

Отсюда λ1

= λ2 = –1,

λ3 = 0. Корни действительные, причем один из

них, а именно

λ = –1 двукратный, поэтому общее решение имеет вид:

y = C1ex + C2 xe x + C3 .

30. Среди корней характеристического уравнения (2) есть комплексные. Пусть для определенности λ1 = α + i β , λ2 = α i β , λ3 = ν + i δ, λ4 = ν i δ, а остальные корни уравнения (2) вещественные и различные.

Отметим, что так как по предположению коэффициенты pi , i = 0, n 1 уравнения (2) вещественные, то комплексные корни этого уравнения попарно сопряженные.

Фундаментальная система решений в этом случае имеет вид

eα x cosβ x , eα x sin β x , eνx cos δx ,

eνx sin δx , eλ5 x , eλ 6 x , ..., eλ n x , (5)

а общее решение:

 

y = C1eα x cosβ x + C2 eα x sin β x +

C3eνx cosδx + C4 eνx sin δx +

+ C5 eλ5 x + ... + Cn eλ n x .

Пример 3. Найти общее решение уравнения y′′′ + 4y′′ + 13y′ = 0 .

Решение. Характеристическое уравнение

λ3 + 4λ2 + 13λ = 0 имеет корни λ1 = 0, λ2 = –2 – 3i, λ3 = –2 + 3i.

3 19

Общее решение:

 

 

 

 

y =

 

C1 + C2 e

2x

cos 3x + C3e2x sin 3x .

40. B

случае, если λ1 = α

+ i β является k-кратным корнем уравне-

ния (2) (k

n

), λ

 

= α i β ,

также будет k-кратным корнем, и фундамен-

2

2

 

 

 

 

 

 

 

 

 

 

 

тальная система решений уравнения (1) (считаем, что остальные корни уравнения (2) вещественные и различные) имеет вид:

eα x cosβ x , eα x sin β x , xe α x cosβ

x , xe α

x sin β

x , ...,

(6)

x k 1eα x cosβ

x , x k 1eα x sin β x ,

eλ 2 k + 1x

, ..., eλ n x .

 

Значит, общее решение ДУ (1):

 

 

 

 

 

y = C1eα x cosβ x + C2 eα

x sin β x + C3 xe α x cosβ x +

C4 xe α

x sin β x + ... +

 

+ C2k 1x k 1eα x cosβ x +

C2k x k 1eα x sin β x + C2k + 1eλ 2 k + 1x

+

... + Cn eλ n x .

 

Пример 4. Найти общее решение уравнения

 

 

 

yV 2y IV + 2y′′′ − 4y′′ + y′ − 2y =

0 .

 

 

 

Решение. Характеристическое уравнение имеет вид

λ5 – 2λ4 + 2λ3 – 4λ2 + λ – 2 = 0.

Запишем его так: (λ – 2)(λ2 + 1)2 = 0. Отсюда получаем λ = 2 – однократный и λ = ± i – пара двукратных мнимых корней.

Общее решение есть

y = C1e2x + (C2 + C3 x ) cos x + (C4 + C5 x ) sin x .

5o. Остальные возможные случаи рассматриваются аналогично 20-

40.

Уравнение Эйлера. Линейные дифференциальные уравнения

вида

xny(n) + pn–1xn–1y(n–1) + ... + p1xy(1) + p0 y = 0,

(7)

где все pi , i = 0, n 1 – постоянные, называются уравнениями Эйлера. Эти уравнения заменой независимого переменного x = et преобразуются в ЛОДУ с постоянными коэффициентами

y (n) + a

n1

y(n1) + ... +

a y′ +

a y = 0 ,

(8)

t

t

1 t

0

 

где yt(i ) i-тая производная функции у по t .

3 2 0