
книги из ГПНТБ / Михельсон В.С. Элементы вычислительной математики учебное пособие для электроприборостроительных техникумов
.pdfРешением дифференциального уравнения (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 |