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

Лекции ДиРУр

.pdf
Скачиваний:
19
Добавлен:
27.03.2015
Размер:
1.98 Mб
Скачать

y n

dy

a(x)y1 n f (x)

(8)

 

 

dx

 

и введем замену переменных z y1 n . Подставив в уравнение (8) выраже-

ние y n

dy

через производную

dz

, приходим к линейному уравнению

 

dx

 

 

 

 

 

 

 

 

dx

 

 

 

 

1

 

dz

a(x)z f (x).

(9)

 

 

 

 

 

 

 

1 n dx

 

С помощью формулы общего решения неоднородного линейного уравнения

(6) находим решение уравнения (9)

zCe(n 1) a(x)dx (1 n)e(n 1) a(x)dx f (x)e(1 n) a(x)dxdx

иобщее решение уравнения Бернулли (7)

 

(n 1) a(x)dx

 

 

 

 

 

 

 

 

 

 

 

 

(n 1) a(x)dx

(1 n) a(x)dx

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

1 n

 

y Ce

 

 

 

 

 

 

 

 

(1 n)e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (x)e

dx

 

. (10)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy

 

y

 

 

 

 

 

1

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

x

 

 

y x3 1

 

 

 

 

 

 

 

 

 

Для данного уравнения a(x)

1

,

 

 

f (x)

 

 

 

1

 

 

 

 

 

, n 1. Общее решение нахо-

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

x3 1

 

 

 

 

дится с помощью формулы (10)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

2

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

2

 

e

x

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

y

Ce

 

 

 

 

 

2e

 

 

 

 

 

 

 

 

 

 

 

dx

.

 

 

 

 

 

 

 

 

 

 

x

3

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из последнего равенства окончательно получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

C

 

 

 

 

ln

x

 

 

1

 

.

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

21

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

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

22

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

 

dy

f (x,y),

(11)

 

 

 

dx

 

 

с известным начальным условием

 

 

 

y(x0) y0 .

(12)

Теорема 1. Пусть в уравнении (11) функция

f (x,y) непрерывна в прямо-

угольнике G:

 

 

x0 a x x0 a,

y0 b y y0 b,

и удовлетворяет в G условию Липшица:

f (x,y1) f (x,y2) L y1 y2

для некоторой постоянной L 0, тогда существует единственное решение уравнения (11), удовлетворяющее начальному условию (12), на отрезке:

 

 

 

x0 h x x0 h,

 

b

, K max

 

f (x,y)

 

 

 

где h min a,

 

 

 

.

 

 

K

(x,y) M

 

 

 

Аналогично можно доказать теорему существования и единственности для системы дифференциальных уравнений

dyi

fi(x,y1, ,yn), yi(0) yi0

i 1,2, ,n .

(13)

 

dx

 

 

Теорема 2. Пусть правые части системы уравнений (13) в области G, определяемой неравенствами:

x0 a x x0 a, yi0 bi yi yi0 bi , i 1,2, ,n

удовлетворяют условиям:

1)функции fi (x,y1, ,yn) i 1,2, ,n непрерывны, а следовательно, огра-

ничены fi K;

23

2) функции fi (x,y1, ,yn) i 1,2, ,n удовлетворяют условию Липшица:

n

fi(x,y1, ,yn) fi(x,z1, ,zn) L yi zi .

i 1

Тогда существует единственное решение системы дифференциальных уравнений (13), состоящее из непрерывных функций y1(x), ,yn(x),определен-

ных на отрезке

 

 

b

b

 

x0 h x x0 h,

где h min a,

1

, ,

n

.

 

K

 

 

K

 

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

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

F x,y,y , ,y(n) 0.

(1)

24

Будем предполагать, что F - непрерывная функция всех своих аргументов; в

окрестности начальных значений x0,y0,y0, ,y(0)(n) выполняются условия:

F x0,y0,y0

(n)

0,

F

 

 

 

 

 

0 .

, ,y(0)

 

 

 

 

 

 

 

y

(n)

 

 

 

 

 

 

 

 

 

 

x x0

,y y0,y y0, ,y

(n)

y

(n)

 

 

 

 

 

 

 

(0)

Тогда согласно теореме о неявной функции можно утверждать, что в этой окрестности уравнение (1) однозначно разрешимо относительно y(n) , т.е.

существует функция f такая, что уравнение (1) эквивалентно уравнению

y

(n)

 

(n 1)

 

.

(2)

 

f x,y,y , ,y

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

n 1

 

y y1, y1 y2,

...........

yn 2 yn 1,

f x,y,y1, ,yn 1 .

Теорема 1. Существует единственное решение дифференциального уравнения n- го порядка (2),удовлетворяющее условиям

y(x0) y0

,

 

 

y

(n 1)

(n 1)

,

 

y

(x0) y0

 

(x0) y0

 

если в окрестности начальных значений (x0,y0,y0, ,y((0)n 1))функция f явля-

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

25

Общим решением дифференциального уравнения n- го порядка называет-

ся функция y (x,C1,C2, ,Cn), которая зависит от аргументов x и n неза-

висимых постоянных C1,C2, ,Cn , обращающая вместе со своими производ-

ными

 

(n)

это

уравнение

в

тождество.

y ,y , ,y

 

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

зывается решение, которое получается из общего решения, если придавать постоянным C1,C2, ,Cn определенные числовые значения.

5.1. Интегрирование уравнения

y(n) f (x).

(3)

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

x

y(n 1) f (x)dx C1 ,

 

 

x0

 

 

 

 

 

 

 

 

 

 

 

 

 

x

x

 

 

 

 

 

 

 

 

 

 

y(n 2) dx f (x)dx C1(x x0) C2 ,

 

 

 

 

 

x0

x0

 

 

 

 

 

 

 

 

 

 

 

(n 3)

x

x

x

 

 

C (x x )2

 

 

 

 

 

y

 

dx dx

f (x)dx

1

0

C2(x x0) C3 ,

 

 

 

2

 

 

 

x0

x0

x0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

……………………………………………………………………………………

 

 

 

 

x

x

x

 

 

C (x x )(n 1)

C

(x x )(n 2)

 

y dx dx

 

f (x)dx

1

0

 

 

2

0

 

 

 

(n 1)!

 

 

(n 2)!

 

x0

x0

x0

 

 

 

 

 

 

 

 

n

Cn 1(x x0) Cn . (4)

Итак, получено общее решение уравнения (2), удовлетворяющее

начальным данным:

 

y(x0) Cn , y (x0) Cn 1

y(n 1)(x0) C1 .

 

26

Если положить в формуле (4) значения констант Ci (i 1, ,n)равными ну-

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

x

x

x

 

y dx dx f (x)dx.

(5)

x0

x0

x0

 

Можно показать, что для интеграла в формуле (5) справедливо представление (формула Коши)

x

x

x

1

x

 

dx dx f (x)dx

(x z)n 1 f (z)dz.

(6)

 

x

0

x

x

(n 1)!x

 

 

0

0

 

0

 

Пример 1.

d3y

sinx.

Начальные значения x0 0,

y0,y0,y0 - некоторые за-

dx3

 

данные числа. Решение находим согласно формулам (4), (6)

 

 

1 x

 

2

1

 

2

 

 

y

 

0

(x z)

 

sinzdz

 

C1x

 

C2x C3

(7)

2

 

2

 

Берем интеграл в (7) методом интегрирования по частям

1

x

 

 

 

 

1

 

 

 

 

x

 

 

(x z)2 sinzdz

(x z)2 cosz

 

0x (x z)coszdz

 

 

 

 

 

 

 

2

2

 

 

0

 

 

0

 

 

x2

 

 

 

x

 

x2

 

 

 

 

 

 

(x z)sinz

 

0x sinzdz

cosx 1.

(8)

 

 

 

 

 

2

 

 

 

 

0

 

2

 

 

 

 

Примем в (7) значения констант равные начальным значениям функции y и

ее производных: C3 y0 , C2 y0 , C1 y0 . С учетом (8) окончательно получа-

ем

 

x

2

 

 

2

 

y

 

cosx 1

y0x

 

y0x y0 .

 

 

 

 

2

2

 

 

 

 

 

27

5.2. Линейные дифференциальные уравнения n-го порядка

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

нение линейное относительно неизвестной функции и ее производных, т.е. уравнение вида

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

(9)

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

a0(x)y(n) a1(x)y(n 1) an(x)y 0.

(10)

Будем предполагать, что решение уравнения (10) ищется на некотором от-

резке a x b и на этом отрезке коэффициент a0(x)

не равен нулю. Разде-

лив уравнение (10) на a0(x), приведем его к виду

 

y(n) p1(x)y(n 1)

pn(x)y 0.

(11)

Обозначим L[y] y(n) p1(x)y(n 1)

pn(x) - линейный дифференциаль-

ный оператор, тогда уравнение (11) можно записать в виде L[y] 0.

Свойства линейного оператора:

1) постоянный множитель выносится за знак оператора L:

L[cy] cL[y];

2)применение оператора L к сумме двух функций равно сумме результатов применения этого оператора к каждой функции в отдельности:

L[y1 y2] L[y1] L[y2] .

Из свойств 1), 2) оператора Lследует

n

n

 

L[ ci

yi ] ciL[yi].

(12)

i 1

i 1

 

28

Тождество (12) позволяет доказать следующее утверждение: линейная

n

комбинация с произвольными постоянными коэффициентами ci yi реше-

i 1

ний y1,y2, ,yn линейного однородного уравнения L[y] 0 является реше-

нием этого уравнения.

Функции y1(x),y2(x), ,yn(x)называются линейно зависимыми на некотором отрезке a x b, если существуют постоянные величины 1, 2, n , не все равные нулю, такие, что для любого x [a,b] выполняется равенство

1y1(x) 2 y2(x) n yn(x) 0.

Теорема 2. Если функции y1,y2, ,yn

линейно зависимы при x [a,b], то на

отрезке [a,b] определитель Вронского

 

 

 

 

 

 

y1

y2

 

yn

 

 

 

 

 

W[y1,y2, ,yn

]

y1

y2

 

yn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y1(n 1)

y2(n 1)

yn(n 1)

 

 

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

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

являются реше-

ниямилинейного однородного уравнения

 

y(n) p1(x)y(n 1) pn(x)y 0

 

с непрерывными на отрезке [a,b]коэффициентами pi(x)

(i 1, ,n), то

определитель Вронского W[y1,y2, ,yn ] не может обратиться в нуль ни в одной точке отрезка [a,b].

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

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

29

симых решений этого уравнения. Фундаментальную систему можно постро-

ить, задавая n2

произвольных чисел yi(k)(x0)

(k 0, ,n 1;i 1, ,n) для лю-

бой точки x0 [a,b] таких, что

 

 

 

 

 

y1(x0)

y2(x0)

yn(x0)

 

 

 

 

 

 

y1(x0)

y2(x0)

yn(x0)

0.

 

 

 

 

 

 

 

 

 

 

 

y1(n 1)(x0)

y2(n 1)(x0)

yn(n 1)(x0)

 

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

yi(x) (i 1, ,n), соответствующие начальным значениям

yi(k)(x0) линейно независимы, т.к. определитель Вронского в точке x0 не ра-

вен нулю. Следовательно, эти решения образуют фундаментальную систему.

Теорема 4. Если линейное однородное дифференциальное уравнение

y(n) p1(x)y(n 1) pn(x)y 0

с действительными коэффициентами pi(x) имеет комплексное решение y(x) u(x) iv(x), то действительная часть этого решения u(x) и его мнимая часть v(x) являются решениями этого уравнения.

Теорема 5. Общим решением на отрезке [a,b] линейного однородного уравнения

y(n) p1(x)y(n 1) pn(x)y 0

с непрерывными на [a,b] коэффициентами pi(x) (i 1, ,n) является ли-

нейная комбинация y n ci yi n линейно независимых на отрезке [a,b]

i 1

частных решений с произвольными постоянными коэффициентами.

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

30