Лекции ДиРУр
.pdfy 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