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

книги из ГПНТБ / Михельсон В.С. Элементы вычислительной математики учебное пособие для электроприборостроительных техникумов

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

Решением дифференциального уравнения (1) называется вся­ кая функция y=f(x), при подстановке которой в это уравнение полу­ чается тождество, т. е.

F( x, f ( x) , Г ( x ) , . . . , f w (*)) = 0 .

Например:

1) (jr(x) = 2sin х — 7х+5 является решением дифференциального уравнения у"—2 y'+ y + 4cosx+7x— 19=0, так как q"(x)—2q'(x) +

+<7(х)-|-4 cosx + 7x—19= 0.

2)Решением уравнения

у ” — 6 =

0 . . .

( 2)

является функция

-f- С2х Сз,

(3)

q (х) — х3 4" Ci х2

так как q"'(x)—6 s 0; Сь С2, и С3 в этом решении являются постоян­ ными числовыми коэффициентами. Выбирая для них различные чис­ ловые значения, можно получить различные функции q\(x), q2 (x),...

и т. д., являющиеся решениями рассматриваемого уравнения. Так, при

С ,= О, С2= 4 и С3

= —1 получаем qi=x3 + 4x—1, а при С ,= —8 и С2=

= С3= 6 получаем

q2 = x3 —8х2+6х + 6, при С! = С2= С3 = 0 получаем

q3 (x) = х3.

 

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

Решение дифференциального уравнения может содержать про­ извольные постоянные, как в примере, рассмотренном выше. У раз­ личных уравнений решения могут содержать различное число таких постоянных.

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

П р и м е р . Рассмотрим уравнение

 

 

(Ру

 

 

- sin х.

 

( 4 )

dx2

 

 

dy

d2y

dz

Обозначим функцию —— через г, тогда

— ■= ——, поэтому за-

(IX

(IX

и х

данное уравнение примет вид

 

 

dz

— = sin х dx

120

или

 

 

 

йг — sin x-dx.

 

 

Интегрируя обе части, получаем

 

 

J dz = J sin xdx

 

или

 

 

 

г =

cos х -f- С:

 

 

11

 

поэтому

 

 

 

ах

 

 

 

Из полученного уравнения находим

 

 

j dy = J (— cos х +

Ci) dx

 

или

 

 

(5)

у = — sin X -f- Cl X -j- C2.

Таким образом, решение рассматриваемого уравнения содержит две произвольные постоянные. В рассмотренном примере 2 реше­ ние содержало три постоянные.

Функция (3) представляет собой семейство решения уравнения (2), функция (5)— семейство решения уравнения (4).

Решение дифференциального уравнения, содержащее столько произвольных постоянных (параметров), каков его порядок, называ­ ется общим решением или общим интегралом этого уравнения.

Как видно из рассмотренных примеров, связь между функцией у и параметрами, входящими в общее решение, выражается в уравне­ нии семейства функции.

Если в таком уравнении выразить у через аргумент х и парамет­ ры С \ , , Сп. то функция y=g(x, Си . . ., С») называется общим ре­ шением уравнения (1), в противном случае оно называется общим интегралом этого уравнения.

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

шение.

 

3

Например,

гл=—sin х,

у2 =ъ\пх--- —лт+17, у3= —sinx+5 явля­

ются частными

решениями

уравнения (4), которые получаются из

121

общего решения (5) при Ci=0 и С2 = 0, при С1= -^ -и С2= 17 и при

С, = 0 и С2=5.

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

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

 

 

 

У — 4>(х, C

i , , Сп)

 

(6)

содержит

п

произвольных

постоянных. Пусть

заданы числа

у0, . . ., у п-\

и требуется найти такое частное решение,

ф(х, сь ..., с„),

что при х = х0

 

 

 

 

 

 

 

 

ф(х. С1........С„) = у 0,

 

 

 

 

 

dy(x,

Ci , . . . ,

Сп)

 

 

 

 

 

 

 

 

 

(7)

 

 

 

d<n~ 1 \ ( x , C 1 }. . . , C n)

 

 

 

 

 

 

dxn~ l

— Уп—\ •

 

 

 

 

 

 

 

 

_

Для_этого достаточно из системы уравнений

(7)

найти значения

Си

, Сп и подставить их

в общее решение

(6).

 

 

Уравнения

(7), при помощи которых определяются значения па­

раметров С.........Сп, называются начальными условиями дифферен­ циального уравнения.

§ 2. Основные методы интегрирования дифференциальных уравнений первого порядка

Для решения таких уравнений применяются различные методы. Выбор метода зависит от того, к какому классу относится решаемое уравнение. Рассмотрим некоторые из них.

122

Уравнение с разделяющимися переменными.

Уравнение первого

порядка

 

 

можно записать в виде

Р(х, У, У') — О

 

 

 

f ( х ,

у) dx = ф ( х , у) dy,

(8)

где левая и правая части

представляют собой

произведение диф­

ференциала переменного на функцию от х и у. Например, уравнение

(х + У) ~ + ху sin х = О

ах

может быть записано как

+ у) dy — ху sin xdx.

В некоторых случаях дифференциальное уравнение (8) можно записать в виде

f(x)dx==(f(y)dy,

(9)

где каждая часть (9) представляет собой произведение функции од­ ной переменной на дифференциал этой переменной.

В этом случае говорят, что в исходном уравнении можно разде­ лить переменные, записав его в виде уравнения (9) с разделенными переменными.

П р и м е р . Найти частный интеграл дифференциального урав­ нения

х3 (1 + У3) dx = (I -(- х4) уЧу,

удовлетворяющий начальному условию у(0 ) 0 .

Р е ш е н и е . Разделим обе части уравнения на (1+ </3) (1+х4), получим

x3dx уЧу

1 + х4 = 1 + у3’

поэтому

Г X3dx ^ Г у2 dy

J l + x4 J l + ( / 3 +

123

или

 

 

- £ - ! n ( l+ * « ) = - y lnd

+*/3) + C.

Преобразуем полученное равенство, получим

(!+**)»

12 С.

(1 + 2/3)4

 

 

Определим значение С, используя

начальное условие 1/(0) =0,

при х=0 и у = 0, получаем 12С= InI=0,

 

поэтому С=0.

Окончательно получаем частный интеграл

( 1 + * * ) * = ( ! +

г/3 ) 4

заданного уравнения, удовлетворяющий начальному условию.

В общем виде уравнения, в которых можно разделить перемен­

ные,

имеют вид

 

 

 

 

 

 

K(x)L(y)dx +

M(x)N(y)dy = 0.

(10)

 

Разделим обе части уравнения на L(y)M(x), получим уравнение

вида

(9) при f (х)= -

К (х)

 

N(y)

Чтобы найти об-

~ и ср ((/) =

.

 

 

взять

^ (у/

от обеих частей

щий интеграл уравнения, нужно

интеграл

уравнения (9).

 

 

 

 

 

 

Рассмотрим уравнение

 

 

 

 

 

f(ax-\- by +

с) dx + dy =

0.

(11)

Оно не является уравнением вида (9), но при помощи замены переменных получается уравнение, в котором можно разделить пе­ ременные. Возьмем z —ax+by+с, тогда dz=adx+bdy, поэтому урав­ нение (11) запишется в виде

/ (г) dx +

dz adx

= 0

 

b

или

(bf(z) a) dx + dz = 0.

Таким образом, получили уравнение вида (10), в котором можно раз­ делить переменные:

dz

dx —

bf (г) —а

124 —

Из этого уравнения находим

Сdz

* = -

J

ГТТТ------

+ с = F (г) + С.

 

bf (г) — а

Поэтому, подставив вместо г сумму ax+by+с, получаем интеграл уравнения (11)

х — F (ах + by + с) + С — 0.

П р и м е р . Найти общий интеграл уравнения

 

 

(4х +

у) 2 dx + dy = 0.

 

 

Р е ш е н и е .

 

В этом

уравнении а = 4, 6=1,

с=0 и г= 4х+ у и

{(г) = г 2; поэтому

 

 

 

 

 

 

 

С

dz

 

С

dz

1

г—2

^

 

bf(z)J — a

J

г2 — '

In ----- .

 

4

г+2

Значит,

 

 

 

 

 

 

 

 

 

1

 

4х 4- у — 2

0

 

 

 

х +

In -----—

-------+ С =

 

'4 4х + у + 2 ^

общий интеграл решаемого уравнения.

Однородное уравнение.

Если в

уравнении (8) функции

1 (х,у) и (р(д:,у) являются однородными

функциями двух переменных

одинаковой степени однородности, то

оно называется однородным

дифференциальным уравнением.

Чтобы определить, является ли дан­

ная функция ф(х, у) однородной, и найти ее степень, необходимо вместо переменных х н у подставить соответственно произведения tx и ty. Если функцию ф(<х, ty) можно преобразовать в виде произ­ ведения первоначальной функции и какой-либо степени числа t, т. е. t ",то она является однородной, а степень однородности равна п. Таким образом, для однородной функции должно выполняться тож­ дество

4?(tx, ty) = tn $(x, у)

для любого числа t.

Например, функция ф(х, У ) — х 2+ х у является однородной вто­ рой степени, так как ф(/л', ty) = (<*)2+ (tx) (ty)=t2 (x2 +xy) = г2ф((х, у ) .

125

Если уравнение является однородным, то в нем можно разде­ лить переменные. Для этого вводят новую переменную

У

и — — .

X

Рассмотрим подробнее метод решения однородных дифферен­ циальных уравнений разделением переменных.

Пусть уравнение (8) является однородным, т. е. функции f{x,y)

и ф (х, у)

являются однородными функциями одинаковой степени од­

нородности. Подставив

вместо х и у в эти функции произведения

1

1

 

 

X•— и у -----, получим

 

 

X

X

 

 

f(x,

у) =

и <р(*, у) =

т Г

 

 

 

где п —ту поэтому рассматриваемое уравнение можно записать как

/ ( 1, — ] dx = Ф

у

Так как и= — , у=их, то dy=udx+xdu, поэтому

X

 

f (1, и) dx = ф (1, и) (udx +

xdu)

или

 

[/(1, и) — (1, и) u]dx = ф (1,

и) xdu.

Получаем уравнение вида (10), в котором можно разделить пе­ ременные:

dx

 

ф (1, и) du

X

f ( I , и) — ф (1, и) и '

Поэтому

 

 

 

Ф (1,

u)du

/(1,

и)

-C = F(u) + C,

ф( 1, и) и

Окончательно получаем интеграл решаемого уравнения в виде

In х — F

У_ = С,

 

х

126

Чтобы составить функции f(\,u) и <р(1,н), нужно в функциях Г(х, У) и qi(x,y) вместо х подставить 1, а вместо у — переменную и. Например, в уравнении

(7х4 + ху3) dx = 4 (х4 + у4) dy

f (х, у) = 7х*+ ху3 и ф (х, у) = 4 (х4 + г/4)

являются однородными функциями со степенью однородности, рав­ ной четырем; из сказанного выше следует

/(1, ц) = 7 + «з и ф (1, и) = Ф(1 + « 4).

Уравнение в полных дифференциалах. Полным дифференциалом du функции и(х, у) от двух переменных х и у называется сумма

ди ди

du = ——dx + ——d y .

дх ду

Например, если «(х, у) = x2+sin(2x+i/), то

du = [2х + 2 cos (2х + у )\dx -f cos (2х + у) dy.

Рассмотрим дифференциальное уравнение f (х, y)dx-Ьф(х, y)dy=0 .

Если функции f(x,y) и ф(х, у) такие, что

df(x, у) _

д(р (х, У) _ Q

(12)

ду

дх

 

то существует такая функция F(x,y), для которой полный дифферен­ циал dF(x,y) равен:

dF (х, у) = / (х, y)dx + <p(x, y)dy.

(13)

Всвязи с этим рассматриваемое уравнение называется уравне­ нием в полных дифференциалах.

Вэтом случае из уравнения (13) определяется функция F(x,y);

х

у

 

Е(х,г/) = |

f ( l , y 0 ) d l + \

Ф(*.Е)<*Б.

Xо

Уо

где х0 и Уо произвольно взятые числа.

Выражение F(x,y)=C и является интегралом рассматриваемой

Функции.

Например, решим уравнение (уех +2x+y2 )dx+ (ех +2xy)dy = 0 .

127

Так как f(x, у) = уе * j-2x+y2 и ф(х, у) =ех -г2 ух, то

df(x, у)

 

дф (х, у)

= е 1' +

2у.

 

ду

 

 

дх

 

 

 

Последнее равенство

показывает,

что заданное уравнение яв­

ляется уравнением в полных дифференциалах.

 

 

Найдем функцию F(x,y),

взяв х0 0 =0:

 

 

X

 

 

 

У

 

=

F (х , у)= J ( У0

+ 2£ + у,)2 d% + J (ех + 2 x1 )

Х о

еХа

 

 

У *

 

 

= у0 ех + х2 + yl х — у0

 

 

 

1

 

 

 

— 4 — Уохо + УеХ + ХУ — Уо е*ху\.

Окончательно находим F(х,

у )= х 2 +уех +ху2,

поэтому

уех + х 2+

+ у2 х = С — общий интеграл заданного

уравнения.

 

Может оказаться, что для заданного уравнения условие (12) не выполняется и оно не является уравнением в полных дифференциа­ лах. Однако иногда уравнение можно превратить в уравнение в полных дифференциалах, если умножить на подходящую функцию.

Интегрирующий множитель уравнения (8 )— такая функция Ъ,(х,у), что от умножения на нее оно обращается в уравнение в пол­ ных дифференциалах. Интегрирующий множитель существует для любого дифференциального уравнения, но часто его трудно опреде­ лить. Задача нахождения интегрирующего множителя облегчается, если он зависит только от одного переменного х и у.

Чтобы это определить, можно

поступить следующим образом.

Определяем функцию

дц>(х, у)

df(x, у)

X (х, у)=

дх

ду

Если %(х,у) = 0, то выполняется условие (12) и заданное урав­ нение является уравнением в полных дифференциалах. Когда функ­ ция Х(х,у) не равна тождественно нулю, то находим

X (х , у)

— v ■■ (14) Ф(*. У)

Если это выражение не зависит от у, то £(х) определяется по формуле

f' х <*• У)

d x

 

U x ) = e J 9 (х, у)

.

• (15)

128

Если выражение (14) зависит

от у, то составляем выражение

М *'

У)

(16)

f(x,

у )

 

и определяем, зависит ли оно от х. Если выражение (16) не зависит от х, то 1 (у) определяется по формуле

_ с U x , у) йу

1(») = е * П х -у) .

Если выражение (14) зависит от х и у, то интегрирующий мно­ житель является функцией переменных х и у.

П р и м е р. Найти общий интеграл уравнения

х2) dx — (х2 у2 + х) dy = 0.

Р е ш е н и е . В рассматриваемом уравнении f(x,y)= y —х2 и

ф(*. 1/ ) = - ( х!!/Ч х).

Проверим, является ли оно уравнением в полных дифференциа-

лах, для этого вычисляем

df (х, у)

 

д у ( х ,у )

------------

и ------------ :

 

ду

 

дх

 

df

dq>

- 2хуг -

I .

— = 1

и —

 

ду

дх

 

 

 

 

Так как

 

 

 

 

 

Х(дг, у) =

2лгг/2 +

2 ф 0,

 

то заданное уравнение не является

уравнением в полных дифферен-

циалах. Составим выражение

^ (х, у)

 

 

—;---- ::

 

 

 

ф (х ,у )

 

 

1 ( х , у ) _

2 (хуг -)- 1)

= _ _ А (

Ф (х, у)

— х(хуг +

1)

х '

Так как составленное выражение не зависит от у, то интегриру­ ющий множитель находим по формуле (15)

_2

6(*) =

X = е - 2 1 п д г

X2 *

 

9 — 4 4 0

129

Соседние файлы в папке книги из ГПНТБ