Теперь решаем задачу Коши |
|
|
|
|
|
|
|
(y(0)0 |
= 1. |
|
|
|
|
||
y = 1 + x, |
|
|
|
||||
Ее решение |
|
|
|
|
|
|
|
x |
|
|
|
|
|
x2 |
|
y3(x) = 1 + Z0 |
|
|
|
|
|
||
(1 + t) dt = 1 + x + |
|
. |
|||||
2 |
|||||||
После n-го шага получим |
|
|
|
|
|
|
|
|
|
x2 |
xn |
n−→ ex. |
|||
yn+1(x) = 1 + x + |
|
+ ... + |
|
||||
2 |
n! |
||||||
|
|
|
|
|
→∞ |
||
1.1. Уравнение с разделяющимися переменными
Уравнением с разделяющимися переменными называют дифференциальное уравнение первого порядка вида
y0 = p(x)q(y), |
(1.6) |
в котором функция p(x) предполагается непрерывной на интервале (a, b), а q(y) – непрерывно дифференцируемой на интервале (c, d). В прямоугольнике
D = {(x, y) | a < x < b, c < y < d}
для уравнения (1.6) выполнены все условия теоремы 1.1 и, следовательно, через каждую точку (x0, y0) D проходит единственная интегральная кривая этого уравнения.
Нетрудно видеть, что если q(y0) = 0, то решением уравнения (1.6) является функция
y(x) ≡ y0, x (a, b). |
(1.7) |
Поэтому при поиске решений уравнения (1.6) прежде всего находят нули функции q(y) и соответствующие им решения данного уравнения вида (1.7).
Например, такими решениями уравнения
y0 = x cos(y)
являются функции y(x) ≡ π2 + πk, k Z, x R.
Предположим теперь, что q(y0) 6= 0, и рассмотрим задачу Коши
(
y0 = p(x)q(y),
(1.8)
y(x0) = y0
10
(разумеется, x0 (a, b), y0 (c, d)). Пусть y = y(x), x (α, β) – решение этой задачи (его существование и единственность гарантируются теоремой 1.1 ), т. е.
y0(x) = p(x)q(y(x)). |
(1.9) |
Ясно, что при всех x (α, β) q(y(x)) 6= 0, ибо в противном случае рассматриваемое решение в некоторой точке совпало бы с решением вида (1.7), что невозможно в силу единственности. Поэтому (1.9) можно разделить на q(y(x)):
y0(x)
q(y(x))
= p(x).
Из равенства функций следует равенство интегралов от них (x (α, β)):
x
Z y0(x) dx
q(y(x))
x0
x
Z
=p(x) dx.
x0
Производя в первом интеграле замену y = y(x) с учетом y(x0) = y0, получаем
y(x) |
|
|
x |
|
|
|
|
Z |
dy |
= Z |
p(x) dx. |
|
(1.10) |
||
|
|
|
|||||
q(y) |
|
||||||
y0 |
|
|
x0 |
|
|
|
|
Пусть Q(y) – некоторая первообразная функции |
1 |
, а P (x) – функции |
|||||
q(y) |
|||||||
|
|
|
|
|
|
||
p(x). Тогда (1.10) с учетом формулы Ньютона–Лейбница можно переписать в виде
Q(y(x)) − Q(y0) = P (x) − P (x0). |
(1.11) |
Поскольку Q(y0) и P (x0) – постоянные, из (1.11) следует, что решение любой задачи Коши (1.8) при q(y0) 6= 0 должно удовлетворять уравнению
Q(y) = P (x) + C, |
(1.12) |
где C – некоторая постоянная. Уравнение (1.12) задает решения уравнения (1.6) неявно. Желательно искомую переменную y выразить через x, что в общем случае не всегда возможно. Вспоминая, что Q и P – это первооб-
разные функций 1q и p, можем понять, почему иногда (1.12) записывают в виде
Z |
q(y) = Z |
p(x) dx. |
(1.13) |
|
dy |
|
|
Заметим, что иногда и решения y(x) ≡ y0, q(y0) = 0 также содержатся в формуле (1.12) при некоторых значениях C.
11
Пример. |
|
|
xy2 |
|
|||
y0 |
|
|
|
||||
= |
|
|
. |
|
|||
x2 + 1 |
|
||||||
В принятых обозначениях p(x) = |
|
|
x |
|
, q(y) = y2 |
, D = R2. Функция q(y) |
|
|
x2 + 1 |
||||||
|
|
|
|
|
|||
имеет единственный нуль y = 0, и, следовательно, данное уравнение имеет единственное решение вида (1.7), а именно
y(x) ≡ 0, x R.
Теперь найдем остальные решения этого уравнения. Записывая для него (1.13):
Z |
y2 = Z |
x2 + 1 |
|
dy |
x dx |
и находя первообразные подынтегральных функций, получаем уравнение вида (1.12):
1 |
1 |
|
|
|
|
|
|
|
|||||
− |
|
= |
|
|
|
ln(x2 + 1) + C, |
|
|
|||||
y |
2 |
|
|
||||||||||
или |
1 |
|
|
|
|||||||||
|
|
|
|
|
|||||||||
y = − |
|
|
|
|
|
|
. |
|
|
||||
C + |
1 |
|
ln(x2 + 1) |
|
|
||||||||
|
|
|
|
|
|
|
|||||||
|
|
2 |
|
|
|
|
|
|
|||||
Таким образом, решениями данного уравнения являются функции |
|||||||||||||
и |
y(x) ≡ 0, x R, |
|
|
||||||||||
1 |
|
|
|
||||||||||
|
|
|
|
(1.14) |
|||||||||
y(x) = − |
|
|
|
, |
|||||||||
C + |
1 |
ln(x2 + 1) |
|||||||||||
|
|
|
|
|
|
|
|
||||||
|
|
2 |
|
|
|
|
|||||||
где C – произвольная постоянная.
Заметим, что интервал (α, β) задания любого решения из найденного семейства зависит от C и может быть определен только по начальным условиям для этого решения.
Рассмотрим, например, решение данного уравнения, удовлетворяющее начальному условию y(0) = 1. Подставляя это условие в (1.14), находим:
y(0) = −C1 = 1,
т. е. C = −1. Следовательно, искомое решение определяется функцией
1
y(x) = ,
1 − 12 ln(x2 + 1)
12
которая задана в трех интервалах (из R исключаются нули знаменателя):
I2 |
= −∞√e2− |
|
|
|
|
|
|
|
|
1 ; |
||||||
|
|
1, |
√−e2 |
|||||||||||||
I1 |
= |
√e2 |
, |
|
√e2 |
1 |
; |
|
|
|||||||
I3 |
= |
|
|
|
|
|
|
|
||||||||
|
1, + . |
− |
|
|||||||||||||
|
|
− |
|
− |
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
∞ |
|
|
|
||
|
|
|
− |
|
|
|
|
|
|
|
|
|
||||
Поскольку любое решение дифференциального уравнения по определению задано на некотором интервале числовой оси, из трех интервалов I1, I2, I3 выбирается тот, которому принадлежит x0 = 0 из начального условия, т. е. I2. Итак, установлено, что решением задачи Коши
|
|
|
|
y0 |
|
|
|
xy2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= |
|
, |
|
|
|
|
|
|
|
|
|
|||
|
|
|
x2 + 1 |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
y(0) = 1 |
|
|
|
|
|
|
|
|
|
|||||
является функция |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
, x − |
|
|
|
|
|
|
|
|
. |
|||
y(x) = |
|
|
|
|
|
e2 |
− 1, |
e2 |
− 1 |
|||||||||
|
|
|
|
|
|
|||||||||||||
1 |
|
1 |
ln(x2 + 1) |
|||||||||||||||
|
− |
2 |
p |
|
|
p |
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Рассмотрим еще один пример, относящийся к электротехнике, а именно процесс исчезновения тока при размыкании цепи, изображенной на рис. 1.4 (ключ K при этом переводится из положения 1 в положение 2).
...
... .
.
RL
..•.
2 |
• |
. . |
|
• |
K |
. . |
|
U 1 |
|
Рис. 1.4
Размыкание цепи приводит к уменьшению силы тока и, следовательно, к возникновению ЭДС самоиндукции UL в катушке индуктивности L. Будем считать, что мгновенные значения силы тока на всех участках цепи одинаковы и подчиняются законам постоянного тока. В этом случае сила тока J может быть найдена из закона Ома:
dJ JR = UL = −L dt .
13
Считая, что в момент размыкания сила тока равнялась J0, приходим к следующей задаче Коши:
−LdJdt = JR,
J(0) = J0.
Для ее решения воспользуемся формулой (1.10):
J(t) |
J |
= − Z0 |
t |
||
Z |
|
L dt. |
|||
|
dJ |
|
|
R |
|
J0 |
|
|
|
|
|
Отсюда
R ln(J(t)) − ln(J0) = −L t
и, следовательно,
−R t
J(t) = J0e L .
Это так называемый экстраток размыкания (ток, проходящий в цепи при отсутствии в ней напряжения, т. е. под действием одной лишь электродвижущей силы самоиндукции). С возрастанием t он стремится к нулю (рис. 1.5).
J . .
J0
|
. . |
..... |
|
0 |
.t |
Рис. 1.5
Замечание. При решении дифференциальных уравнений бывает удобно записывать произвольную постоянную, принимающую все действительные значения, в виде некоторой функции g(C), C (a1, a2), множество значений которой на интервале (a1, a2) совпадает с множеством R. В качестве такой функции можно взять, например, ln(C), C (0, +∞) или tg(C),
C −π2 , π2 , либо другую с указанием соответствующего интервала.
Пример. Найти кривые y = y(x), для которых отрезок касательной,
14