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

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

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

Умножим заданное уравнение на |(х ):

у — х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

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