книги из ГПНТБ / Михельсон В.С. Элементы вычислительной математики учебное пособие для электроприборостроительных техникумов
.pdfУмножим заданное уравнение на |(х ):
у — х2 , |
х 2 ц 2 + х |
(17) |
---- — dx — |
------ ------ dy = 0. |
Теперь оно стало уравнением в полных дифференциалах, в чем нетрудно убедиться, вычислив
|
У- |
х2 |
у2 |
+ X |
|
Х(х, у) — |
ду |
+' |
дх |
“ |
I2 +(~ *) = 0' |
Уравнение |
(17) |
запишем в |
виде |
|
|
|
-У— dx — dx — у2 dy — — dy — 0 |
||||
|
|
, „ , |
xdy — ydx |
Л |
|
|
— dx — y2 dy— ------------ = О, |
||||
поэтому
d\ — x — — — — ) = 0 .
Это значит, что
или
За-2 + ху3 +3у + ЗхС = 0
является общим интегралом уравнения (17).
§3. Линейные дифференциальные уравнения первого порядка
Клинейным дифференциальным уравнениям первого порядка относятся уравнения вида
f ( x ) ~ + v ( x) y - i q ( x) = 0> |
(18) |
—130
в которых неизвестная функция у и ее производная |
входят линейно, |
|||||||||
в первой степени. |
|
|
dy |
|
|
|
|
|
||
Например, (3x2+sinx) |
|
|
|
|
|
|||||
——+ (cos х )у + хъ = 0 является линейным |
||||||||||
уравнением. |
|
|
|
ах |
|
|
|
|
|
|
|
т(х) |
|
ё (х) |
|
|
|
|
|||
Пусть |
<pi (х) |
qi{x)= |
тогда |
уравнение запи- |
||||||
— ■ |
|
I \х) |
||||||||
|
|
/ (х) |
|
|
|
|
|
|||
шется в более простом виде |
|
|
|
|
|
|
||||
|
|
dii |
|
|
‘7iW = |
0. |
|
|
|
|
|
|
-Г- + Ф1 (х)У + |
|
|
|
|||||
|
|
ах |
|
|
|
|
|
|
|
|
Нахождение общего |
интеграла. |
Решим |
это уравнение |
методом |
||||||
подстановки. |
Для |
этого |
произведем |
замену |
переменных |
y — uv, где |
||||
и a v неизвестные пока функции от х. Тогда |
dy |
du |
" ' |
dv |
||||||
dx |
dx |
dx |
||||||||
и уравнение (18) молено записать как |
|
|
|
|
||||||
|
du |
г dv |
+ (fi(x) v) + q ^x) = 0 |
. |
|
(19) |
||||
|
v ~ |
+u \ — |
|
|||||||
Выберем o(x) |
таким образом, чтобы |
|
|
|
|
|||||
|
|
|
dv |
+ Ф: (x) v = 0 . |
|
|
|
|
||
|
|
|
— |
|
|
|
(20) |
|||
|
|
|
dx |
|
|
|
|
|
|
|
Из уравнения (20) |
получаем |
|
|
|
|
|
|
|||
|
|
|
dv |
|
(x) dx7 |
|
|
|
|
|
|
|
|
— = — |
|
|
|
|
|||
|
|
|
v |
|
|
|
|
|
|
|
или |
|
|
|
|
|
|
|
|
|
|
|
|
In V |
(x) = — f Ф1 (x) + In C. |
|
|
|
||||
(Для удобства дальнейших преобразований постоянную интегриро вания берем в виде In с.)
Поэтому
0 |
1 |
e JГtpi {x)dx |
(21) |
v = C e ~ h ‘ |
dx или |
|
|
|
^ 0 |
|
|
131
Учитывая равенство (20), уравнение (19) можно переписать в виде
|
du |
|
|
|
|
|
( 22) |
|
dx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
поэтому |
|
|
|
|
|
|
|
|
du ■ |
|
■qi (х) dx |
|
|
||
ГГ |
1 |
f <Pi ( X ) dx |
■q1 |
"I |
_ |
(23) |
|
u = \ |
— — |
eJ |
|
(x)\dx + C1. |
|||
Окончательно получаем общее решение |
|
|
|
||||
y = tw = - е~ Р‘ <*> dx [ J e.f |
{x)dx qi (x) dx] + |
C0 Cj e~ I |
<x)dx . |
||||
В общем решении остается один параметр C= C0Ci.
Если мы хотим найти частное решение уравнения (18), удовлет воряющее начальному условию у1 х^ Х!1 =Уо, то вместо равенств (21) и (23) определяем
|
X |
(x)dx |
. j 4>i |
||
: Сп е |
|
1 |
Х о |
----q1 (x)dx + yot |
|
|
_ |
— I V |
|
*0 |
|
Функция y(x)=vu является частным интегралом уравнения (19), |
||
так как при х = х 0 |
функция v тождественно равна Со и поэтому функ |
|
ция и равна i/o, т. е. у(х0) =у0. |
||
Рассмотренный |
способ подстановки позволяет свести задачу |
|
интегрирования одного линейного дифференциального уравнения (18) |
||
к нахождению решений двух уравнений с разделяющимися перемен ными (20) и (22).
Пр и м е р . Найти общий интеграл уравнения |
|
||||
|
|
du |
(хг + |
2) = 0. |
(24) |
|
|
х —q— + 3(/ + |
|||
Р е ш е н и е . |
|
dx |
|
ср(х) =3 |
и q{x) = (х2 +2), |
В этом уравнении {(х) |
|||||
поэтому фх (х) — |
3 |
хг |
2 |
|
|
— |
и <7х ( х )= -------- . |
|
|
||
132
Из (21) н (23) находим
— общий интеграл уравнения (24).
Однородные линейные дифференциальные уравнения первого
порядка. |
Если в уравнении (18) функция |
q(x) тождественно рав |
на нулю, |
то получаем однородное уравнение |
|
f {х) ~7~ +ф(*)2/ = 0 .
ах
В этом уравнении можно разделить переменные
dy ____у(х)
У~ Цх)
поэтому
. |
Г ФМ |
■dx + In С, |
|
т. с. |
— J / м |
||
|
|
|
|
|
- ( |
f(v) |
|
|
У = Со е |
J |
|
является общим решением однородного уравнения.
§ 4. Линейные дифференциальные уравнения высших порядков
Линейным дифференциальным уравнением «-го порядка называ ется уравнение вида
У{п) + a1 (x)yin~ l) + а 2 (х )у (п~ 2)-]------- |
1- а ^у '+ а „(х )у = Ь (х ), (25) |
—133
которое содержит искомую функцию и ее производные только в пер вой степени и не содержит их произведений.
Например, уравнение
d*y |
. |
d2 У |
dy |
|
у = 5 cos х |
||
|
— х2 sin |
х — —-f tg х |
—— + х3 |
||||
|
dx4 |
|
dx2 |
dx |
|
|
|
является |
линейным |
дифференциальным |
уравнением четвертого по |
||||
рядка, а уравнение |
|
|
|
|
|
|
|
не является линейным уравнением, |
так |
как |
содержит произведения |
||||
d3 у |
d2 у |
I dy |
v2 |
|
|
dy |
|
у на ~7Т |
и -Т 7 11а (” ;— |
и квадрат производного —— . |
|||||
и л |
и л |
\ |
и Л |
/ |
|
|
и л |
Общий интеграл |
линейного уравнения |
« го порядка содержит |
|||||
п параметров С,,. . Сп. Для нахождения этих параметров задают п начальных условий:
у U*. = У 0 ’ |
У' U *0= »!•■ ■ •• |
I*=*„ = Уп- 1 • • • |
(26) |
||||
Т е о р е м а . |
Пусть |
коэффициенты |
а,(х),.. , ап(х) линейного |
||||
дифференциального уравнения п-го порядка |
(25) |
непрерывны |
на |
||||
некотором отрезке [с, d) |
и (26)— произвольная |
система начальных |
|||||
условий. Если х0 |
принадлежит интервалу (с, |
d), |
то |
существует одна |
|||
и только одна функция у={(х), которая определена и непрерывна на всем интервале (с, d), является решением уравнения (25) и удов летворяет начальным условиям (26).
Если в уравнении (25) Ь(х) = 0, то оно называется однородным,
при Ь(х) Ф 0 — неоднородным.
Однородное линейное дифференциальное уравнение. Однород ное линейное дифференциальное уравнение n-го порядка имеет вид
У{п) + ах ( х )у {п~ Х) + |
а2 |
(х) у <'г~ 2) + ----- Ь «л_ i (■*) У' + |
|
|
+ |
ап (х) у=* 0, |
(27) |
Частные решения этого уравнения обладают рядом свойств: |
то и |
||
1. Если функция У \ ( х ) |
является решением уравнения (27), |
||
функция Су\(х) является решением этого уравнения.
134
2. Если функции у\(х) и у2 (х) являются решениями уравнения (27), то и функция y ( x ) - y i ( x ) + y 2 (x) является решением этого Уравнения.
3. Если iji{x), у2 (х) , ..., у п (х) — решения уравнения (27), то н Функция
у (х) = Ci yi (х) + С2 Уг (х) + • • • + Сп уп (х).
есть решение этого уравнения.
П р и м е р . |
Функция г/i (х) =sinx |
является решением линейного |
||||
однородного уравнения |
|
|
|
|
||
так как (sin jc)"+sin х = 0. |
У" + У = 0, |
|
(28) |
|||
также является решением этого |
уравне |
|||||
Функция у ( х ) = С sin х |
||||||
ния, так как (С sinx)"+(C |
sinx)s9. |
|
|
|||
Рассматриваемое уравнение имеет два частных решения: i/,= |
||||||
= sinjc, y2 =cosx. Составим |
их линейную комбинацию; |
|
||||
|
Уз (х) — Ci sin х + |
С2 cos х. |
|
|||
Убедимся, |
что У з ( х ) также |
является решением этого уравнения: |
||||
у3 -\-у3 |
= (Сг sin x-t- С2 |
cos х) |
-|- (Cj sin х + Cg cos х) |
= |
||
= — Ci sin x — C2 cos x + Ci sin x + C2 cos x = 0.
Рассмотренное свойство 3 показывает, что при помощи частных решений уравнения (27) можно составить решение, содержащее п параметров. Является ли такое решение общим решением уравнения (27)? Решение уравнения (27), содержащее п параметров, является общим решением только в том случае, если при помощи любых за данных п начальных условий можно однозначно определить значения
всех п параметров. |
|
|
является |
|
решением уравнения |
Например, у\=Схsin х+С2 sin х |
|
||||
второго порядка (28). Возьмем начальные условия |
|||||
У |
х= TZ |
4 и |
ТС |
= |
2. |
|
|
|
|||
|
4 |
|
4 |
|
|
Попытаемся определить значения Ci и С2, чтобы у\ удовлетворя ло этим начальным условиям:
у'х = Схcos х + С2 cos х,
135
поэтому
У1 |
/ 2 |
- СС, Т—2— - l - c . V |
|
у 1 |
У 2 |
=С, V + c s I72" = 2. |
Составленная система несовместна, т. е. при любых значениях Ci и С2 решение г/i не может одновременно удовлетворять обоим началь ным условиям. Это означает, что решение у\ не является общим ре шением.
Решение (/2=С, sin х+С 2 cos х уравнения (28) является общим решением этого уравнения. Чтобы в этом убедиться, рассмотрим на
чальные условия у2 (х0)=а и у±(х0 )=Ь. Составим систему уравне ний для нахождения С| и С2;
У2 (*о) = С, sin х„ + С2 cos .v0 ■-= а,
У'2 ( *(>) “ ^1 |
Ху Сд Sin Х у |
b . |
|
|
Из этой системы находим |
I sin хп а I |
|
||
I a |
cos х0 |
|
||
I Ь —sin х0 |
| C O S |
Хд Ь | |
|
|
sin х0 |
cos х0 |
I sin х0 |
cos |
I |
cos х0 |
—sin Хд |
| cos х0 |
—sin |
I |
При любых значениях a, b и Хд составленная система совместна, так как при любых значениях х0 знаменатель этих дробей равен —1.
sin х |
cos х |
|
|
|
|
= — sin2 x — cos2 x = — 1 |
|
cos х —sin x |
|
||
Ci = |
a sin x0 -|- b cos x0 и C2 = a cos x0 |
— b sin x0. |
|
Таким образом, |
решение i/2=C, sin x+C 2 cos x является общим |
||
решением уравнения (28). |
у2 =С{sin х+С2 cos х. |
||
Сравним |
функции У\ = С\ sin x-f-C2 sin х н |
||
Хотя обе эти функции содержат по два параметра и степень урав нения (28) равна двум, первая не является общим решением этого уравнения, а вторая является.
136
Если у\(х) , ..., у п(х) являются частными решениями уравнения
(27), то
у (х) = С-1 У! (х) + • • • + Сп уп (х)
будет общим решением этого уравнения тогда и только тогда, когда функции у i(д:), .. ,,ijn(x) линейно независимы.
Линейная зависимость системы функций определяется аналогич но линейной зависимости системы столбцов матрицы, рассмотренной в гл. II.
Пусть У\(х).......Уп(х) — система функций, |
определенных |
и Не |
|||||
прерывных на одном п том же отрезке [с, d\. |
такие п— 1 |
чисел |
|||||
О п р е д е л е н и е |
1. |
Если |
существуют |
||||
«1, • . . , а , _ 1, |
а |
1+1, . . . |
а „ , |
что |
|
|
|
yt (х) |
== |
ttj ух (х) н-------- 1- а |
ух_ х(х) + |
у 1+х (х) + |
|
||
+ ' " + о п уп {х).
для всех х, принадлежащих отрезку [с, d], то говорят, что функция У,■(х) линейно выражается через остальные функции системы на от резке [с, d].
О п р е д е л е н и е 2. Система функций У\.......у п называется ли нейной независимой на отрезке [с, d], если ни одна из функций этой
системы |
не выражается |
линейно через остальные на этом отрезке. |
|
О п р е д е л е н и е 3. |
Если хотя бы одна функция y t |
из системы |
|
У \ , - - , У п |
линейно выражается через остальные функции |
на отрезке |
|
[с, d], то система называется линейно зависимой на этом отрезке.
П р и м е р ы . 1) |
Функция |
у\ = tgx+21nx+3 линейно выража |
ется через функции |
J/2=4tgx+sin х и (/3= 0,25sinx—21пх—3 на ин |
|
тервале (0; + °о ), так как |
|
|
|
У\ — |
Уч + (—О Уз |
при х>0 Так как функция У\ линейно выражается через функции у2 и у3,
то |
система (/,, у2 |
и у3 линейно зависима |
на интервале (0; |
оо). |
|
[с, |
2) Система |
функций у\=\, у2=х и |
у3 = х2 на любом |
отрезке |
|
d] линейно независима. Действительно, если У\ линейно выража |
|||||
ется на отрезке [с, d] через у2 |
и у3, то |
|
|
||
|
|
У\ = |
ач -Уч + a.i Уз■ |
|
|
137
Это значит, что для всех х, принадлежащих отрезку [с, d], долж но выполняться равенство
1 = а2 х + «з х2. |
(29) |
Это выражение является квадратным уравнением, поэтому ра венство (29) возможно только при двух значениях х.
Можно также проверить, что функция у2 не выражается линейно через у\ и у2 ни на каком отрезке и функция у3 также не выража ется линейно через у\ и у2. Окончательно приходим к выводу, что функции системы у\, у2 и у2 на любом отрезке линейно независимы.
Приведем теперь правило, позволяющее |
определить, является |
||||
ли заданная система функции |
|
|
|
||
У\(х), |
|
у2 (х), |
у„(х). |
|
(30) |
линейно независимой. |
|
|
|
|
|
Для этого составим определитель п-го порядка |
|
||||
|
|
У\ |
У 2 • • . |
Уп |
|
|
|
h |
у'2 ... |
У п |
|
|
|
|
|
|
|
W{yu . • , Уп) |
= |
У1 |
у"ч • •' |
Уп |
(31) |
|
|
|
|||
у[п
вкотором каждый столбец состоит из одной функции системы (30)
ивсех ее^ производных до (п—1)-го порядка включительно.
Такой определитель называется определителем Вронского си
стемы функций (по имени польского математика Вронского, впервые изучавшего его свойства).
Т е о р е м а . Если функции системы (30) линейно зависимы на отрезке [с, d], то определитель Вронского (31) системы (30) тожде ственно равен нулю.
Т е о р е м а . Если функции системы (30) линейно независимы на отрезке [с, d] и являются частными решениями некоторого линей ного однородного дифференциального уравнения, то для всех х, принадлежащих этому отрезку, определитель Вронского (31) систе мы (30) не равен нулю.
Эти теоремы позволяют определить, является ли система функ ций (30) линейно независима или нет.
П р и м е р . Система функций
yi — sin х и уг = cos х.
является линейной независимой на всей числовой оси, так как оп ределитель Вронского этой системы (см. стр. 136) не равен нулю при
любых значениях х.
Если система функций (30) линейно независима на отрезке [с, dl а каждая функция является на этом отрезке частным реше нием некоторого линейного однородного дифференциального уравне ния, то такая система называется фундаментальной системой реше- •
ний этого уравнения.
Найдя фундаментальную систему решений уравнения (27), мож но сразу построить общее решение.
Т е о р е м а . Если система функций |
|
Уу(х) ..........Уп(х). |
(32) |
является фундаментальной системой решений линейного однородного дифференциального уравнения п-го порядка
y in) + а, (х) y (n~ u + • • • + а„ (х) у = 0, |
(33)' |
то линейная комбинация функций системы (32)
У — Ci Уг + • • • + Сп уп
является общим решением уравнения (33).
Таким образом, нахождение общего решения уравнения (33) сводится к нахождению фундаментальной системы решений, т. е. частных линейно независимых решений.
Метод нахождений таких решений мы рассмотрим на примере решения однородных линейных уравнений с постоянными коэффи
циентами.
Однородные линейные уравнения с постоянными коэффициентами. Если в уравнении (33) все коэффициенты а,(х) являются постоян ными числами, то оно называется однородным линейным дифферен циальным уравнением с постоянными коэффициентами. Такое урав нение можно записать в виде
У(п) + «1 У{п~ Х) Ч-------\-°п- 1 У + ап У =" °- |
(34) |
где а, — постоянные числа.
139
