3447
.pdfИз курса высшей математики, изучаемого в 1-ом семестре, нам известно, что в алгебре большую роль играют теоремы о числе решений алгебраических уравнений. Так, в случае линейных систем уравнений, определитель которых не равен нулю, система имеет единственное решение. Известна также основная теоремы алгебры, утверждающая что многочлен n -й степени всегда имеет ровно n корней на множестве комплексных чисел (считая корни с их кратностями). В теории дифференциальных уравнений также важнейшим является вопрос о том, насколько много решений имеет дифференциальное уравнение. Выше мы убедились, что дифференциальное уравнение имеет бесконечное множество решений. Поэтому ставится вопрос не о числе решений, а о том, как можно описать совокупность всех решений данного уравнения, как можно выделить интересующее нас решение и единственно ли решение, описываемое интегральной кривой, проходящей через точку (x0 , y0 ) области.
Будем говорить, что задача Коши с начальными условиями (1.10) имеет единственное решение, если
существует |
такое число h |
0 , что |
в |
интервале |
x |
x0 |
|
h |
||
определено |
решение |
y |
(x) такое, |
что |
(x0 ) y0 |
и |
не |
|||
существует |
решения, |
определенное |
в |
этом |
интервале и |
не |
совпадающее с решением |
y |
(x) хотя бы в одной точке |
||
|
h , отличной от точки x x0 (рис. 3). |
|||
интервала |
x x0 |
|||
|
y |
|
|
|
|
|
|
y |
(x) |
|
y0 |
M 0 |
|
|
0 |
|
x0 |
x |
|
|
|
|
Рис. 3 |
|
10
В случае, если задача Коши с начальным условием (1.10) имеет не одно решение или совсем не имеет решений, говорят,
что в точке (x0 , y0 ) нарушается единственность решения
задачи Коши.
Вопрос о единственности решения задачи Коши исключительно интересен и для теории дифференциальных уравнений, и для многочисленных практических приложений. Зная, что решение задачи Коши единственно, мы, отыскав решение, удовлетворяющее заданным начальным условиям, будем уверены в том, что получили определенный, единственный закон, описывающий физическое явление, и других траекторий развития процесса при данных условиях нет.
2. В этом разделе будут сформулированы достаточные условия существования и единственности решения задачи Коши (теорема Пикара).
Ниже будут рассматриваться только непрерывно дифференцируемые задачи (1.10).
Итак, ставится вопрос: каким условиям достаточно подчинить правую часть дифференциального уравнения (1.7) в окрестности начальных данных x0 , y0 , чтобы через точку
P0 (x0 , y0 ) проходила одна и только одна интегральная кривая
этого уравнения?
Теорема. Пусть |
дано уравнение (1.7) |
dy |
f (x, y) , |
||||||||||
|
|
||||||||||||
dx |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||
и поставлено начальное условие (1.10) |
|
|
|
||||||||||
|
y |
y0 |
|
при |
x x0 . |
|
|
|
|||||
Предположим, |
что |
функция f (x, y) |
определена в |
||||||||||
некоторой замкнутой ограниченной области (рис. 4) |
|
||||||||||||
D : |
|
x x0 |
|
a, |
|
y y0 |
|
b |
|
|
|
||
|
|
|
|
|
|
|
сточкой (x0 , y0 ) внутри области ( a, b - заданные
положительные числа) и удовлетворяет в области D следующим условиям.
11
|
1) |
Функция |
f (x, y) непрерывна и, |
|
следовательно, |
|||||||
ограничена (поскольку D замкнутая область), т.е. |
|
|||||||||||
|
|
|
|
|
|
f (x, y) |
|
|
M , |
|
|
(1.11) |
|
|
|
|
|
|
|
|
|
|
|||
где |
M - постоянное положительное число, |
|
P(x, y) |
- любая |
||||||||
точка области D . |
|
|
|
|
|
|
|
|
|
|||
|
2) |
Функция |
f (x, y) имеет ограниченную |
частную |
||||||||
производную по аргументу y , т.е. |
|
|
|
|||||||||
|
|
|
|
|
|
f (x, y) |
|
K , |
|
|
(1.12) |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
y |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
||
где |
K - постоянное положительное число, |
P(x, y) |
- любая |
|||||||||
точка области D . |
|
|
|
|
|
|
|
|
|
|||
|
При этих условиях уравнение (1.7) имеет единственное |
|||||||||||
решение |
y |
(x) , |
|
удовлетворяющее начальному |
условию |
|||||||
(1.10). |
Это |
решение |
|
определено |
и |
непрерывно |
дифференцируемо в некоторой окрестности начального значения x0 независимой переменной x ; именно, оно заведомо определено в интервале
|
|
|
x x0 |
h, |
(1.13) |
|
где h - наименьшее из чисел a и b M |
|
|||||
|
|
h |
min a,b M . |
(1.14) |
||
|
y |
|
|
|
|
|
y0 |
b |
|
|
|
|
|
|
|
P0 |
|
|
||
|
y0 |
|
|
|
|
|
|
|
|
|
D |
|
|
y0 |
b |
|
|
|
|
|
|
|
|
|
|
||
|
0 x0 |
a x0 |
h x0 x0 h x0 a |
x |
||
|
Рис. 4. Иллюстрация |
|
||||
|
Области в теореме Пикара |
|
12
Доказательство теоремы Пикара мы не приводим. Из теоремы следует, что если правая часть уравнения (1.7) –
полином относительно переменных |
x и y или любая другая |
|
функция, определенная и непрерывная относительно x |
и y |
|
вместе с частной производной по y |
при всех значениях |
x и |
y , то через любую точку P0 (x0 , y0 ) |
проходит одна и только |
одна интегральная кривая, поскольку в любом прямоугольнике D с центром в точке P0 (x0 , y0 ) оба условия теоремы Пикара
будут выполняться. В этом случае, вся плоскость XOY будет заполнена не пересекающимися и не касающимися друг друга
гладкими интегральными кривыми.
2. ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ ПЕРВОГО ПОРЯДКА, ИНТЕГРИРУЕМЫЕ В КВАДРАТУРАХ
Введем терминологические дополнения.
Интегрированием дифференциального уравнения назовем совокупность операций выполняемых для нахождения общего решения.
Интегрируемость дифференциального уравнения в квадратурах означает, что в общий или частный интегралы дифференциального уравнения входят интегралы от функций, зависящих от y и x (по отдельности), которые не обязательно
вычисляются через элементарные функции. Например, для дифференциального уравнения
y xy2 2xy
общий интеграл представляется квадратурами
dy |
xdx C . |
|
|
||
y( y 2) |
||
|
В данном случае интегралы вычисляются в элементарных функциях.
13
2.1. Дифференциальные уравнения с разделяющимися переменными
1. К дифференциальным уравнениям, интегрируемым в квадратурах, в первую очередь принадлежат уравнения с
разделяющимися переменными. В этом случае функция f (x, y) ,
в уравнении представляет собой произведение функции
переменного x |
на функцию переменного y : |
|
|
|
|||||||
|
|
|
dy |
f (x) ( y) . |
|
|
|
(2.1) |
|||
|
|
|
|
|
|
|
|
||||
|
|
|
dx |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
||
Теорема |
1. |
Если в |
уравнении |
с |
разделяющимися |
||||||
переменными |
(2.1) |
функции |
f (x) |
и |
( y) |
непрерывны в |
|||||
интервалах (a,b) и (c, d) соответственно |
и |
(y) 0 , |
то |
||||||||
общий интеграл уравнения выражается в квадратурах: |
|
||||||||||
|
|
|
dy |
|
f (x)dx |
C , |
|
|
(2.2) |
||
|
|
|
|
|
|
|
|
||||
|
|
|
( y) |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
||
причем |
заданным |
|
начальным |
условием |
y(x0 ) |
y0 |
|||||
определяется |
единственное |
решение этого |
уравнения, |
где |
|||||||
M 0 (x0 , y0 ) - |
любая |
точка |
прямоугольника |
(a,b) (c, d) |
на |
плоскости XOY .
Частный интеграл уравнения (2.1), определяющий решение, удовлетворяющее начальному условию может быть записан в виде
y |
dy |
x |
|
|
|
|||
|
|
f (x)dx |
С . |
|
(2.3) |
|||
|
|
|
|
|
|
|||
|
|
( y) |
|
|||||
y0 |
x0 |
|
|
|
||||
|
|
|
|
|
|
|||
Доказательство. |
|
Поскольку по условию теоремы |
f (x) и |
|||||
( y) непрерывны при |
x (a,b), y |
(c, d) и |
(y) 0 при |
|||||
y (c, d) , то умножая обе части уравнения (2.1) |
на dx |
( y) , |
||||||
мы добиваемся разделения переменных |
|
|
|
|||||
|
|
dy |
f (x)dx . |
|
|
(2.4) |
||
|
|
|
|
|
|
|
||
|
|
|
( y) |
|
|
|||
|
|
|
|
|
|
|
14
Функция |
1 |
|
и f (x) |
непрерывна, и следовательно, |
||||
|
|
|
||||||
|
( y) |
|||||||
|
|
|
|
|
|
|
||
имеют первообразные |
|
|
|
|
|
|||
|
Ф( y) |
|
dy |
, |
F (x) |
f (x)dx. |
||
|
|
( y) |
||||||
|
|
|
|
|
|
|
|
Поэтому равенство (2.4) можно записать следующим образом:
dФ( y) dF(x). |
(2.5) |
В формуле (2.5) y рассматривается как функция от x ,
определяемая дифференциальным уравнением. Из равенства дифференциалов двух функций следует, что эти функции отличаются на постоянную:
Ф( y) F(x) C ,
или
dy |
f (x)dx C . |
(2.6) |
|
|
|||
( y) |
|||
|
|
Докажем, что соотношение (2.6) представляет собой общий интеграл уравнения (2.1). Запишем соотношение (2.6) в
виде |
|
|
|
|
|
|
|
|
|
|
|
|
|
G(x, y) 0, |
|
|
|
|
|
|
|
(2.7) |
|||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
G(x, y) F(x) Ф( y) |
C. |
|
|
|
|
(2.8) |
|||||||
Соотношение (2.7) удовлетворяет условиям теоремы о |
|||||||||||||
неявной функции: производные |
|
G |
f (x) , |
|
|
G |
1 |
- |
|||||
|
|
|
|
|
|
|
|
|
|||||
|
x |
|
|
y |
|
( y) |
|||||||
|
|
|
|
|
|
|
|
||||||
непрерывны в области a x |
b , c y |
d ; |
G |
0 |
. Поэтому |
||||||||
|
|
||||||||||||
y |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||
уравнение (2.7) определяет y |
как функцию от x , непрерывную |
и дифференцируемую. При этом имеет место равенство:
y |
G |
|
G |
f (x) ( y) . |
|
|
|
||
x |
|
y |
||
|
|
|
15
Таким образом, функции, определяемые уравнением (2.7), а следовательно и уравнением (2.6), являются решениями данного дифференциального уравнения (2.1). Эти функции получаются при различных значениях C , выбор которых определяется начальными условиями (x0 , y0 ) . При
фиксированном значении параметра C решение единственное. Справедливо и обратное утверждение. Всякая функция, являющаяся решением уравнения (2.1), т.е. тождественно ему удовлетворяющая, обязана удовлетворять и вытекающему из
него соотношению (2.6).
Следовательно, формула (2.6) действительно определяет общий интеграл уравнения.
Из (2.6) видно, что любые начальные условия (x0 , y0 ) из
прямоугольника a |
x |
b , |
c y d |
|
однозначно определяют |
||
надлежащие значения C . |
|
|
|
|
|
||
Действительно, из (2.7), |
|
(2.8) получаем |
|||||
|
|
C Ф( y0 ) F (x0 ), |
|||||
следовательно |
|
|
|
|
|
|
|
Ф( y) |
F (x) Ф( y0 ) |
|
F (x0 ). |
||||
Поскольку же |
Ф ( y) |
1 |
0 |
, |
то функция Ф(y) |
||
|
|
||||||
|
( y) |
||||||
|
|
|
|
|
|
|
допускает однозначное обращение. Обозначив обратную
функцию Ф 1 , |
получим |
|
|
искомое частное |
решение, |
|||
удовлетворяющее заданному начальному условию |
|
|||||||
y |
Ф 1(F (x) |
|
|
F (x0 ) Ф( y0 )) |
(2.9) |
|||
т.о., решение задачи Коши единственно и теорема доказана. |
||||||||
Пример. Найти общее решение уравнения: |
|
|||||||
|
|
dy |
1 |
y 2 |
|
|||
|
|
|
|
|
|
|
. |
|
|
|
dx |
1 |
x2 |
|
16
Решение. |
В представленном уравнении f (x) |
|
|
1 |
, |
||
|
|
|
|||||
1 |
x2 |
||||||
( y) 1 |
y2 . |
Эти функции непрерывны при |
|
x |
|
, |
|
y |
, т.е. на всей плоскости XOY , функция |
1 |
0 . |
||||
|
|
||||||
( y) |
|||||||
|
|
|
|
|
Таким образом, условия теоремы 1 выполнены. Разделяя переменные, получим уравнение с разделенными переменными
dy |
|
dx |
|
|
|
|
. |
1 y 2 |
1 x2 |
Отсюда получаем соотношение
d |
dy |
d |
|
dx |
. |
|
|
|
|||
1 y 2 |
|
x2 |
|||
|
1 |
|
Освобождаясь от дифференциалов, имеем общий интеграл:
|
|
dy |
|
|
dx |
C; |
arctgy arctgx |
C. |
|
||||
|
|
|
|
|
|
|
|
||||||
1 |
y 2 |
1 |
x2 |
|
|||||||||
|
|
|
|
|
|
|
|||||||
Разрешив общий |
интеграл |
|
относительно |
y , |
получим |
||||||||
общее решение уравнения: |
|
|
|
|
|
|
|
||||||
y tg(arctgx C) , |
|
|
|
arctgx C |
|
. |
|
||||||
|
|
2 |
2 |
|
|||||||||
Отметим, что в полученном общем решении значения |
|||||||||||||
параметра |
C |
(“произвольной |
|
постоянной”) |
|
не |
вполне |
произвольны. При выбранном интервале изменения независимой переменной x , параметр C может принимать бесконечное множество значений, но таких, чтобы не нарушались указанные неравенства.
Можно провести дальнейшее упрощение общего решения, если воспользоваться формулой для тангенса суммы дуг и обозначить C1 tgC . В результате получим
y |
|
x |
C1 |
. |
|
|
|
||
|
1 |
C1x |
17
Пример. Найти частное решение предыдущего дифференциального уравнения, удовлетворяющее начальному условию y(0) 1.
Решение. |
Подставляя в общее решение x |
0 , y |
1, |
||||||||||
получим C1 1. Искомое решение задачи Коши имеет вид |
|
||||||||||||
y |
1 |
x |
, |
или |
y |
1 |
|
2 |
. |
|
|
||
|
|
|
|
|
|
||||||||
|
|
1 |
x |
|
|
|
|
1 |
x |
|
|
||
Интегральной кривой является гипербола, смещенная и по |
|||||||||||||
оси OX , и по оси OY . |
|
|
|
|
|
|
|
|
|
|
|
||
2. Доказанные в теореме 1 утверждения установлены в |
|||||||||||||
предположении, что |
(y) |
0 |
ни при каком значении y |
из |
|||||||||
рассматриваемой |
области. |
Если |
( ) |
0, |
где |
(c, d) , |
то |
видно, что дифференциальное уравнение (2.1) имеет решение
y |
(в этом можно убедиться непосредственно подстановкой |
||||||||||||||
в уравнение y |
). |
|
|
|
|
|
|
|
|
|
|
|
|||
|
Однако |
при |
y |
интеграл |
|
dy |
не |
существует по |
|||||||
|
|
|
|
||||||||||||
|
|
( y) |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
крайне мере |
как собственный |
интеграл. |
Поэтому |
решение |
|||||||||||
y |
не входит в состав общего интеграла (2.6). |
|
|
|
|||||||||||
|
Если в дифференциальном уравнении |
y |
f (x) |
( y) |
при |
||||||||||
c |
d |
( |
) |
0, то уравнение, кроме общего интеграла, |
|||||||||||
имеет еще решение, не получающееся из общего. |
|
|
|
||||||||||||
|
Будет |
ли это решение |
y |
|
особым (т.е. |
будет |
ли в |
||||||||
каждой его точке |
нарушаться условие единственности) – этот |
||||||||||||||
вопрос требует специального рассмотрения. |
|
|
|
|
|||||||||||
|
Пример. Найти решение дифференциального уравнения |
||||||||||||||
|
|
|
|
|
y |
xy2 |
2xy. |
|
|
|
|
|
|||
|
Решение. |
Запишем |
данное |
|
уравнение |
|
в |
виде |
|||||||
y |
xy(2 |
y) . Отсюда |
видно, |
что |
функции |
y |
0 |
и y |
2 |
||||||
являются решениями уравнения. |
|
|
|
|
|
|
|
|
|||||||
|
Остальные |
решения |
уравнения |
найдем, |
|
применив |
процедуру разделения переменных и интегрируя уравнение с разделенными переменными:
18
dy |
xdx , |
( y 0, y |
2) . |
|
|
||||
y( y 2) |
||||
|
|
|
Для вычисления интеграла в левой части равенства разложим правильную рациональную дробь на простейшие
дроби |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
A |
|
|
B |
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y( y |
2) |
|
|
y |
|
|
y |
2 |
|
|
|
|
|
||||||||||||||
Отсюда, приравняв числители дробей в левой и правой |
||||||||||||||||||||||||||||||||||
частях равенства, получаем |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
A( y 2) By 1. |
|
|
|
|
|
|
|
|
|
|||||||||||||||||
Решив систему уравнений |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
A |
B |
0, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
2A |
1, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
имеем: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dy |
|
|
|
|
|
1 |
|
|
dy |
1 |
|
|
|
dy |
|
|
|
|
1 |
|
|
y |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ln |
|
|
|
. |
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
y( y 2) 2 y 2 y 2 2 |
y 2 |
|
|
|
|||||||||||||||||||||||||||||
Общий интеграл уравнения имеет вид |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
ln |
|
y |
|
|
|
|
x2 |
ln C , |
C |
0 , |
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
y |
2 |
|
|
|
|
|
|
|
1 |
|
|
|
1 |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
где обозначили 2C lnC1 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
Отсюда имеем |
|
|
|
|
y |
|
|
|
|
|
x 2 |
|
|
y |
|
|
|
|
~ x 2 |
~ |
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
C e |
, |
|
|
|
|
|
|
|
Ce |
, (C |
C |
0). |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
y 2 |
1 |
|
|
|
|
|
y |
2 |
|
|
|
|
|
|
|
|
|
1 |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
Решение y 0 |
можно |
|
|
получить |
из |
|
общего |
интеграла, |
||||||||||||||||||||||||||
расширив |
область |
|
|
|
|
|
|
|
|
изменения |
|
|
|
параметра |
||||||||||||||||||||
~ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
~ |
|
|
0 . Но решение y |
2 |
|||||||||||
C , включив в эту область значение C |
|
|
нельзя получить из общего интеграла ни при каком значении этого параметра.
Таким образом, вся совокупность решений рассматриваемого дифференциального уравнения включает общий интеграл и решение, не входящее в этот интеграл:
y Ce x 2 ( y 2); C ( , ); y 2.
19