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

Математический_анализ_2_семестр

.pdf
Скачиваний:
23
Добавлен:
24.05.2015
Размер:
2.05 Mб
Скачать

 

Пример.

Рассмотрим

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

уравнение второго

порядка

y

f x,y,y .

Область

существования и

единственности

решения

G R3 x,y,y заполнена непересекающимися интегральными кривыми.

Через

любую

точку

 

 

 

проходит единственная интегральная

кривая.

x0 , y0 , ,y0 G

 

Однако

через

«точку»

x

0

,y

0

R 2 x,y

проходит

бесконечно

 

много

интегральных кривых, все они различаются значениями

 

 

что в

y0 . Заметим,

R 3 x,y,y «точка» x0 ,y0 R 2 x,y представляет собой прямую x x0 ,

y y0 .

Понижение порядка дифференциальных уравнений.

Мы умеем аналитически решать всего пять типов дифференциальных уравнений первого порядка: с разделяющимися переменными, однородные, линейные, Бернулли, в полных дифференциалах. Причем однородные, линейные и Бернулли тоже сводятся к уравнениям с разделяющимися переменными.

Даже решить уравнение второго порядка, не говоря уж об уравнении n-го порядка – проблема. Поэтому стараются понизить порядок дифференциального уравнения, если это возможно, чтобы свести его к известным типам уравнений первого порядка.

Если правая часть дифференциального уравнения n-го порядка зависит только от x, то интегрируя его n раз, можно получить решение.

y n f x ,

y n 1 f x dx Cn 1 ,

y n 2 f x dxdx Cn 1 x Cn 2,...

y x ...

f x dx...dx Cn 1

x n 1

 

x n 2

 

Cn 2

 

... C1 x C0 .

 

 

 

 

n 1 !

 

n 2 !

Но это – очевидный случай. Рассмотрим менее очевидные случаи.

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

1)

Уравнение

не

содержит

явно

y

 

,

 

его

 

вид

 

 

 

 

 

0

или

 

 

 

 

F x,y , y

 

 

y f x,y .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь применяется подстановка

y

 

p x , y

 

 

 

 

-

 

вводится

новая

 

 

p x

 

функция y p x

старой переменной.

Уравнение сводится к уравнению

первого порядка p f x, p .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример. Найти

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

уравнения

 

 

 

 

 

 

и его частное

y x ln x y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

решение, удовлетворяющее начальным условиям y e 0, y

e 1.

 

 

 

 

 

 

dp

 

dx

, p C1 ln x,

y

 

C1 ln x,

y

C1 x ln x x C2

 

 

 

 

 

 

p x ln x p,

 

 

p

x ln x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-

общее

 

 

 

 

решение.

 

Найдем

 

 

 

 

частное

 

 

решение.

 

 

1,

y e e ln e e C2

 

e e C2

C2

0 .

 

 

Частное

y e C1 ln e C1

 

 

 

решение y x ln x x .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2)

Уравнение

не

содержит

явно

x

 

,

его

 

вид

 

 

 

 

 

0

или

 

 

 

F y,y , y

 

y f y,y .

51

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p y ,

 

 

 

 

 

dp

 

 

 

dp dy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

dy dx

Здесь применяется

 

подстановка

y

 

y

 

 

p y p y -

вводится новая функция

 

y p y

новой переменной.

Уравнение сводится к

уравнению первого порядка pp f y, p .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример.

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

yy y 2

0 и его частное

решение, удовлетворяющее начальным условиям

y 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

1 1 .

 

 

 

 

 

 

ypp p 2

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C1

 

Либо p 0 y C - решение,

либо yp

p 0,

ydp pdy,

p y ,

 

 

 

 

 

 

 

 

 

 

y 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ydy C1dx,

 

 

 

C1 x C2 - общее решение.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Найдем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

частное

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

решение.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

2

 

 

 

 

 

1

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y 1 C1

1,

 

y 1 1

 

 

1 C2

 

 

 

 

C2

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y 2

2x 1 - частное решение.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Однородное уравнение относительно y,y ,y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Уравнение называется однородным относительно

 

 

 

 

 

, если при замене

y,y ,y

 

y ky,

y

 

 

 

 

y

 

ky

 

уравнение не изменится.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ky ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь применяется подстановка y yz x .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример. Найти общее решение уравнения xyy x y 2

 

yy

 

 

 

 

 

y

 

yz,

y

 

 

 

 

 

 

 

 

yz

2

 

 

xy yz

2

yz

 

xy

2

z

2

y

2

z,

 

 

 

y z z y

 

z y,

 

 

 

 

 

 

xy 2 z y 2 z,

 

 

y 0 - решение.

xz z,

 

 

 

 

dz

 

 

dx

 

,

z C1 x,

 

y yC1 x ,

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy

 

 

 

 

 

 

 

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C x,

 

y C

eC1

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

- общее решение.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

1

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3)Уравнения, обе части которых являются полными производными каких-либо функций.

Пример.

yy y 2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln y ,

 

 

 

 

 

 

 

 

y

y ,

 

ln y C,

Запишем уравнение в виде

ln y

ln y

 

y C1 y,

 

dy

C1dx,

ln y C1 x C2 ,

y C3eC1x .

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

Существуют еще несколько случаев, которые встречаются реже и здесь не рассматриваются.

52

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

Линейное однородное дифференциальное уравнение n –ого порядка с

переменными коэффициентами может быть записано в виде a0 x y n a1 x y n 1 ... an 1 x y an x y 0

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

 

 

 

a0 x y n a1 x y n 1 ... an 1 x y an x y f x .

 

 

 

 

 

 

 

 

 

 

Если коэффициенты и правая часть – непрерывные функции и a0 x 0 ,

то условия теоремы Коши выполнены,

решения однородного и неоднородного

уравнений существуют и единственны.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Введем линейный дифференциальный оператор

 

 

 

 

 

 

 

L

 

p, x a

 

x

d n

a

x

d n 1

... a

 

x

d

a

 

x a

 

x p n

...a

 

x p a

 

x

n

0

 

 

n 1

 

n

0

n 1

n

 

 

 

 

dxn

1

 

dxn 1

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь

p обозначает оператор дифференцирования

 

d

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

Ln p, x y 0 ,

 

 

Тогда линейное однородное уравнение можно записать в виде

 

а

линейное неоднородное – в виде Ln p, x y f x .

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как Ln p, x линеен, то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ln p, x y1 y2 Ln p, x y1 Ln p, x y2 ,

Ln p, x y Ln p, x y .

 

 

Пользуясь линейностью оператора, легко доказать теоремы о свойствах решений однородного и неоднородного уравнений (ниже обозначено yo -

решение однородного уравнения, yн - решение неоднородного уравнения).

Теоремы о свойствах решений.

1)сумма или разность решений однородного уравнения есть решение однородного уравнения,

2)разность решений неоднородного уравнения есть решение однородного уравнения,

3)сумма решений однородного и неоднородного уравнений есть решение неоднородного уравнения.

Докажем эти теоремы.

1)L yo1 yo2 Lyo1 Ly

2)L yн1 yн2 Ly н1 Ly

3)L yo yн Lyo Ly н

o2

н2

0

f x f x 0 0 f x f x .

Теорема. Решения линейного однородного уравнения с переменными коэффициентами образуют линейное пространство.

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

Решения образуют аддитивную группу по сложению (абелев модуль). В самом деле, ассоциативность по сложению очевидна, y 0 (тривиальное решение) является решением однородного уравнения, для каждого решения

53

y1 , y2 ,... yn

y x противоположное решение y x тоже является решением.

Следовательно, решения однородного уравнения – группа по сложению. Аддитивность решений очевидна, поэтому эта группа аддитивна. Справедливость четырех аксиом из восьми показана. Существует число «1», такое что 1 y x - решение, справедлива ассоциативность по умножению на

число y y . Это – две аксиомы относительно операции умножения на число. Наконец, справедливы две аксиомы дистрибутивности, связывающие

операции

сложения

и

умножения

на

число

y1 y2 y1 y2 ,

y y y .

 

 

 

Итак, налицо полный набор из восьми аксиом. Продумайте их еще раз

подробнее дома.

 

 

 

 

 

 

Линейная зависимость и независимость.

 

 

Функции g1 x , g 2

x ,... g n x называются линейно независимыми, если

1 g1 x ... n gn x 0 1

0,... n 0

(допустима

только

тривиальная

линейная комбинация функций, тождественно равная нулю). В отличие от линейной независимости векторов здесь тождество линейной комбинации

нулю, а не равенство. Это и понятно, так как равенство линейной комбинации

нулю должно быть выполнено при любом значении аргумента.

Функции g1 x , g 2 x ,... g n x

называются линейно зависимыми, если

существует не нулевой набор констант (не все константы равны нулю) 1 ,... n ,

такой что g

x ... g

n

x 0

( 2

... 2

0) (существует нетривиальная

1 1

n

 

1

n

 

линейная комбинация функций, тождественно равная нулю).

Теорема. Для того чтобы функции были линейно зависимы, необходимо и достаточно, чтобы какая-либо из них линейно выражалась через остальные (представлялась в виде их линейной комбинации).

Докажите эту теорему самостоятельно, она доказывается так же, как аналогичная ей теорема о линейной зависимости векторов.

Определитель Вронского.

Определитель Вронского для функций вводится как

определитель, столбцами которого являются производные этих функций от нулевого (сами функции) до n-1 го порядка.

 

y1

y2 ...

yn

 

W x

y '

y

' ...

y

'

 

1

 

2

 

n

.

...

...

...

 

y n 1

y

n 1

y

n 1

 

 

1

 

2

 

n

 

Теорема. Если функции y1 x , y2 x ,... yn x линейно зависимы, то W x 0

Доказательство. Так как функции y1 x , y2 x ,... yn x линейно зависимы, то какая-либо из них линейно выражается через остальные, например,

y1 x 2 y2 x ... n yn x . Тождество можно дифференцировать, поэтому

54

y k x

2

y

k x ...

n

y

k x ,

k 1,2,3... n 1 . Тогда первый столбец

1

 

2

 

n

 

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

Теорема. Для того, чтобы решения линейного однородного дифференциального уравнения n-ого порядка были линейно зависимы, необходимо и достаточно, чтобы W x 0 .

Доказательство. Необходимость следует из предыдущей теоремы. Достаточность. Зафиксируем некоторую точку x0 . Так как W x0 0 , то

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

k, C1 ,...Ck

0,...Cn

, что выполнены соотношения

C1 y1 x0 ... Ck yk

x0 ... Cn yn x0 0

 

 

 

C y x

0

... C

k

y

x

0

... C

n

y

n

x

0

0

 

 

 

1 1

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

C y

n 1 x

... C

y

n 1 x

0

... C y

n 1 x

0

 

0.

1 1

 

 

0

 

 

 

k

 

k

 

 

 

 

 

 

n

 

n

 

 

Так как линейная комбинация решений линейного однородного уравнения

является его решением, то можно ввести решение вида

y x

C1 y1

x ... Ck yk x

... Cn yn x

 

-

линейную комбинацию

решений с теми же коэффициентами. ЗаметШм

55

следовательно, W x 0 .
Доказательство.
зависимы,

что при x x0 это решение удовлетворяет нулевым начальным условиям, это

следует из выписанной выше системы уравнений. Но тривиальное решение линейного однородного уравнения тоже удовлетворяет тем же нулевым

начальным условиям. Поэтому из теоремы Коши следует, что введенное решение тождественно равно тривиальному, следовательно,

y x C1 y1 x ... Ck yk x ... Cn yn x 0,

Ck 0 ,

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

 

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

Если W x0 0 , то решения линейно

Теорема. 1. Для линейной зависимости решений необходимо и достаточно W x 0 (или W x0 0 ).

2. Для линейной независимости решений необходимо и достаточно

W x0 0 .

Доказательство. Первое утверждение следует из доказанной выше

теоремы и следствия. Второе утверждение легко доказывается от противного.

Пусть решения линейно независимы. Если W x0 0 , то решения линейно

зависимы. Противоречие. Следовательно, W x0 0 x0 .

W x 0 ,

Пусть W x0 0 . Если

решения линейно зависимы, то

следовательно, W x0 0 ,

противоречие. Поэтому решения

линейно

независимы.

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

Отличие определителя Вронского от нуля является критерием линейной независимости решений линейного однородного уравнения.

Теорема. Размерность пространства решений линейного однородного уравнения n-ого порядка равна n.

Доказательство.

1. Покажем, что существуют n линейно независимых решений линейного однородного дифференциального уравнения n-го порядка. Рассмотрим решения y1 x , y2 x ,... yn x , удовлетворяющие следующим начальным

условиям:

y1 x0 1, y2 x0 0,... yn x0 0,

y1 x0 0, y2 x0 1,...yn x0 0,

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

y1 n 1 x0 0, y2 n 1 x0 0,...yn n 1 x0 1,

56

Такие решения существуют. В самом деле, по теореме Коши через точку x0 , y0, y0 , y0 ,...ynn 1 проходит единственная интегральная кривая –

решение. Через точку x0 , 1,

0, 0, ...0 проходит решение y1 x , через точку

x0 , 0, 1, 0, ...0 - решение

y2 x , через точку x0 , 0, 0, 0, ...1 - решение

yn x .

 

 

 

 

 

 

 

0...

0

 

 

 

1

 

Эти решения линейно независимы, так как W x

0

1...

0

1 0 .

 

 

0

0...

1

 

 

 

 

 

 

 

2. Покажем, что любое решение линейного однородного уравнения линейно выражается через эти решения (является их линейной

комбинацией).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y x с

Рассмотрим

два решения.

 

Одно

-

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

решение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

y

 

,

y

 

,... y

n 1

 

 

 

 

 

начальными условиями x

 

 

 

 

 

 

. Справедливо соотношение

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

0

 

 

0

 

 

 

0

 

 

 

 

 

 

y x0 C1 y1 x0

C2 y2 x0

... Cn yn x0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C1 y1 x0

C2 y2

x0 ... Cn yn

x0

 

 

 

 

 

y x0

 

 

 

 

 

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

 

 

 

 

 

y n 1 x

C y

n 1 x

C y

 

n 1 x

 

 

... C y

n

n 1 x

, где

 

 

 

0

 

1 1

 

0

 

 

 

2

 

2

 

 

 

 

0

 

 

n

0

 

 

C y

0

, C

2

y

0

...C

n

y

 

n 1 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Второе решение – это линейная комбинация решений y1 x ,... yn x

с теми

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x ... Cn yn x .

 

 

 

же коэффициентами y x C1 y1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x0

 

 

 

 

 

Вычисляя начальные условия в точке

для решения

y x , убеждаемся,

что они совпадают с начальными

условиями

для решения

y x .

Следовательно,

по

теореме

 

 

Коши,

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

решение

y x

представляется в виде линейной комбинации линейно независимых

решений

y1 x ,... yn x

y x y x .

 

 

 

Таким, образом, существует n линейно независимых решений линейного однородного дифференциального уравнения n-ого порядка, и произвольное решение линейно выражается через эти решения . Поэтому размерность пространства решений линейного однородного дифференциального уравнения n-ого порядка равна n. dim I n .

Любые n линейно независимых решений линейного однородного дифференциального уравнения n-ого порядка представляют собой базис пространства решений или фундаментальную систему решений.

Теорема о структуре общего решения однородного уравнения.

Общее решение линейного однородного уравнения есть линейная комбинация решений фундаментальной системы.

yoo x C1 y1 x ... Cn yn x .

Доказательство. Покажем, что линейная комбинация

57

yoo x C1 y1 x ... Cn yn x является общим решениям (удовлетворяет

пунктам определения общего решения)

1.yoo x - решение линейного однородного уравнения комбинация решений.

2.Зададим произвольные начальные условия y0 , y0 ,...y0

что можно подобрать константы C1 ,...Cn такие,

как линейная

n 1 , покажем, что yoo x

 

удовлетворяет этим начальным условиям.

 

yoo

x0

C1 y1

x0

 

...

Cn yn

x0

y0 .

 

 

 

 

 

y

oo

x

0

C y

x

0

...

C

n

y

n

x

0

 

y

.

 

 

 

 

 

1 1

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

y

oo

x

0

C y x

0

...

C

n

y

x

0

y

.

 

 

 

 

1

1

 

 

 

 

 

 

n

 

 

 

 

 

0

 

 

 

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

 

y

oo

n 1 x C y

n 1 x ... C y

 

n 1 x

0

y

n 1 .

 

 

 

 

0

1 1

 

 

 

0

 

 

 

 

n

 

 

 

n

 

 

 

0

Это – система линейных алгебраических уравнений относительно констант C1 ,...Cn . Определитель этой системы – определитель Вронского.

Он не равен нулю, так как решения y1 x ,... yn x линейно независимы. Поэтому константы C1 ,...Cn определяются из этой системы по начальным условиям – правым частям системы единственным образом.

Следовательно, yoo x C1 y1 x ... Cn yn x - общее решение.

Замечание. Определитель Вронского (как всякий определитель) представляет собой ориентированный n – мерный объем, натянутый на векторы решений фундаментальной системы решений.

Формула Остроградского – Лиувилля.

Рассмотрим линейное однородное уравнение a0 x y n a1 x y n 1 ... an 1 x y an x y 0 .

Определитель Вронского можно вычислить по формуле Остроградского – Лиувилля

 

 

 

 

 

 

 

a1

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W x Ce

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a0

x

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вывод формулы Остроградского – Лиувилля.

 

 

 

 

 

 

 

 

 

 

 

Известна формула для производной определителя

 

 

 

 

 

 

 

 

 

 

 

a11 x

a12 x

...

a1n x

 

a

 

a

...

a

 

 

 

 

a11

a12 ...

a1n

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

x

 

 

 

x

 

11

12

 

1n

 

 

 

 

 

 

 

 

 

 

 

d

a

21

a

22

...

a

2n

 

a21

a22 ...

a2n

 

...

 

a

21

a

22

a

2n

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

...

 

 

 

...

...

 

...

 

...

 

...

...

 

 

 

...

...

...

 

 

 

 

an1 x

an2 x

...

ann x

 

a

n1

a

n2

...

a

nn

 

 

 

an1

an2

ann

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

58

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

y

2

...

 

 

y

n

 

 

 

 

y '

 

y' ...

y'

 

 

 

 

 

 

 

 

 

 

 

 

 

dW x

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

n

 

 

 

 

Вычислим

 

 

d

 

 

y1'

 

 

 

 

y2'

 

 

yn'

 

 

 

 

y1'

 

 

y2'

yn'

 

...+

 

 

dx

 

dx

 

...

 

 

 

 

 

...

 

...

 

 

 

...

...

...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y n 1

 

y n 1

y n 1

 

 

y n 1

 

y n 1

y n 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

2

 

 

 

 

 

n

 

 

 

 

 

1

2

n

 

 

 

 

 

 

y1

 

y2 ...

 

 

 

yn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y1

 

 

 

 

 

 

 

...

 

 

yn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y '

 

 

 

 

 

 

 

...

 

 

y '

 

 

 

 

 

'

 

 

'

 

 

 

 

 

 

'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y1

 

y2 ...

 

 

 

yn

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

...

 

 

...

 

 

 

...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

...

 

 

 

 

 

 

 

...

...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a1

 

y1n 1

 

 

 

 

an

 

 

 

 

 

 

a1

ynn 1

 

 

an

 

 

 

 

n

 

 

n

...

 

 

n

 

 

 

 

...

 

y1

 

...

...

yn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y1

 

y2

 

yn

 

 

 

 

 

 

a0

 

 

 

 

 

 

 

 

 

 

a0

 

 

 

 

 

 

a0

 

 

a0

 

 

 

 

 

 

 

 

 

 

 

 

y1

 

 

 

...

 

 

 

 

 

yn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

'

 

 

 

 

 

...

 

 

 

 

 

y

'

 

 

 

 

 

 

 

 

a1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

W x .

 

 

 

 

 

 

 

 

0+...+0+

 

 

 

 

 

...

 

 

 

 

 

...

 

 

 

 

 

...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a1

 

y n 1

 

 

 

a1

 

y n 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a0

 

 

1

 

 

 

 

 

 

 

 

 

a0

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dW x

 

 

 

a1

x

 

 

W x Ce

a1

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

a0

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

0

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Замечание. В формуле Остроградского – Лиувилля участвуют только

коэффициенты при двух старших производных.

 

 

 

 

 

 

 

 

 

 

Рассмотрим частный случай уравнения второго порядка.

 

 

 

 

a0 x y a1

x y a2

x

y 0 . Здесь формулу Остроградского – Лиувилля

можно вывести проще. Рассмотрим y1 x ,

y2 x - два частных решения

 

 

 

a1

 

 

 

 

 

 

 

 

x y1

0 . ,

 

 

 

 

 

 

 

a1

 

 

a2 x y2

0 .

 

Умножим

a0 x y1

 

x y1 a2

 

 

a0 x y2

x y2

 

первое уравнение на y2 , а второе на y1

и вычтем первое уравнение из второго.

 

 

a

 

 

 

 

 

 

 

y

y

 

 

 

a

 

 

 

 

y

 

y

 

 

 

 

 

 

 

 

 

 

x y y

 

 

 

 

 

x y y

 

 

 

0 .

 

 

 

 

 

 

 

 

0

 

 

1

2

 

 

 

2 1

 

 

1

 

1

 

2

 

 

 

 

 

2 1

 

 

 

 

 

 

 

Так как W x

 

 

y1

 

y2

 

y y

'

y

 

 

 

 

'

,

 

 

 

'

'

''

' '

''

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y '

 

y '

 

2

 

 

y

 

то W x y1 y

2

y1 y2 y

2 y1 y

2 y1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

y y ''

y

2

y '' .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теперь уравнение можно переписать в виде a0 x W x a1 x W x 0 .

Решая

это

 

уравнение

 

с

разделяющимися

переменными,

 

получаем

формулу

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a1

x

 

 

 

 

 

 

 

Остроградского – Лиувилля W x Ce

 

 

 

dx

 

 

 

 

 

 

a0

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Формула для построения второго частного решения по известному

 

 

 

(построение фундаментальной системы).

 

 

 

 

 

 

 

 

 

 

 

 

 

y1

 

y2

 

 

 

 

 

 

 

 

 

y ' Ce

a1

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W x

 

 

y y ' y

a0

x

dx .

 

 

 

 

 

 

 

 

 

 

 

 

 

y '

 

 

y '

 

 

 

 

1

 

2

 

 

2 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Разделим обе части уравнения на

y 2

x 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

59

 

y y'

y

y'

 

 

 

y

 

 

 

 

 

 

 

1

 

 

a1 x

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a0 x

 

 

 

 

 

 

 

 

 

1 2

 

2 1

 

 

 

 

2

 

 

 

 

C

 

 

 

e

 

.

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

y1

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y1

 

 

 

 

 

 

 

 

 

 

 

 

y1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y2

 

 

 

 

 

 

 

 

 

e

a1

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

dx

 

 

 

 

Отсюда

 

 

 

 

 

C

 

a0

x

dx C

. Нам надо найти частное решение,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

y 2

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a1

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

a0

x

поэтому выберем С=1, C

1=0, получим

y

2

 

e

 

dx .

y 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

Теорема о структуре общего решения неоднородного уравнения.

Общее решение линейного неоднородного уравнения есть сумма частного решения линейного неоднородного уравнения и общего решения однородного уравнения.

yон x yчн x yоо x .

Доказательство. Покажем, что yон x yчн x yоо x - общее решение неоднородного уравнения.

1.yон x yчн x yоо x - решение неоднородного уравнения как сумма

решений однородного и неоднородного уравнений (теоремы о свойствах решений).

2.

Зададим

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

 

 

начальные

 

условия

 

 

 

x0 , y0 , y0 ,...y0

n 1 .

 

Вычислим

начальные

 

условия

 

для

выбранного

частного решения

 

неоднородного

 

 

уравнения

 

 

 

y

чн

x

0

,

y

 

x

0

,...y

чн

т 1 x

0

.

Получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

чн

 

 

 

 

 

 

 

 

 

 

систему линейных алгебраических уравнений для определения

 

констант:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yoo x0 C1 y1 x0 ... Cn yn x0 y0

 

yчн x0 .

 

 

 

 

 

 

 

 

 

 

 

 

y

x

0

C y

x

0

... C

n

y

n

x

0

y

 

y

x

0

.

 

 

 

 

 

 

 

 

oo

 

 

 

 

 

1 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

чн

 

 

 

 

 

 

 

 

 

 

 

 

y

 

x

0

 

C y x

0

... C

n

y

x

0

 

y

 

 

y

 

 

x

0

.

 

 

 

 

 

 

 

oo

 

 

 

 

 

1 1

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

0

 

 

 

чн

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

n 1 x

 

C y n 1 x

 

... C y

 

n 1 x

y

n 1

y

n 1 .

 

 

 

 

 

oo

 

 

 

 

0

 

1

1

 

 

0

 

 

 

 

 

 

n

 

n

 

 

 

 

 

 

0

 

 

 

 

0

 

 

 

 

 

чн

 

 

 

 

Определитель этой системы – определитель Вронского. Он не равен нулю,

так как

решения

y1

x ,... yn

x

 

 

 

линейно

 

 

независимы.

Поэтому

константы

C1 ,...Cn

определяются

из этой системы

 

по

начальным

условиям

– правым

частям системы единственным образом. Следовательно,

yон x yчн x yоо x -

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Метод вариации произвольной постоянной для линейного

неоднородного

 

 

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

 

 

 

 

 

 

 

уравнения

 

 

n-ого

 

порядка.

y n a1 x y n 1

... an 1 x y an x y f x . ( y n

Ly f x ).

 

 

 

 

Здесь обозначено

L a1 x y n 1

... an 1 x y an x y , заметим, если y -

решение однородного уравнения, то y n

 

Ly .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

60