Материал: Diff_lektsii__1

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

Теперь решаем задачу Коши

 

 

 

 

 

 

 

(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