Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
1235.pdf
Скачиваний:
7
Добавлен:
15.11.2022
Размер:
12.71 Mб
Скачать

Поскольку p = sliw,

получим ch u=t V P1Jr 1-

Учитывая выражения

для sh и и ch и через

показательные функции,

находим

 

откуда

 

е» = К > + 1 +/>,

e-« = V p * + l - p ,

 

так что

 

 

 

г = С, {V p -+ 1 + р)п+

С.Л V p 1+

 

р)п.

 

 

 

1 -

 

Для X (р) получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

( V Р 2 +

1 +

р ) ' г +

С г ( У Р 2+ 1 ~~ р)п

 

 

 

 

 

 

 

r >

+ 1

 

 

 

 

 

 

 

 

 

 

 

 

При

п = 0 из

(12)

найдем

 

 

 

 

 

 

Выбирая

Cj + С2 =

1,

получим решение х (/) = J0(/) — бесселеву функ­

цию первого рода порядка нуль.

С ^ О ,

^ 2=

1,

получим решение

Полагая

л = 1

и

выбирая

x(t) — J l (t)

уравнения

( 10).

 

 

 

 

 

 

Найти ре1иения уравнений:

742.lx" + ( 2 t - \ ) x ' + ( t - l ) x = 0.

743.*х" + 2.г' = 0.

744.

x"

(i -\-1) x9-\- tx —0,

v (0) =

1, *'(0)=

—1.

 

745.

x" + (t + b)x' = 0, A(0) = —1,xf (0)—

0 (b —любое

действительное

число).

 

 

 

 

 

 

746.. хГ+

tx9- (/ -I-1) A = 0,

v (0) -

 

x9(0) =

1.

 

 

747.

а" — lx9-^nx — Q, n —целое,

/i> 0 (уравнение

Че­

бышева — Эр*мита):

 

 

 

 

 

 

 

а)

x (0) =

1,

x9(0) - 0, / I -

2k:

 

 

 

 

 

б)

x (0) = 0,

х* (0) —11 /i = 2k+l.

 

 

 

 

 

§

16. Интеграл

Дюамеля

 

 

 

 

 

 

Если

функция /(0

непрерывна на [0,

°о),

а

функция

ф(0

непрерывно дифференцируема на 10, + со)

и

 

 

 

 

 

 

 

 

F (р) Vs*/ (0.

ф (Р)

Ф(0.

 

 

 

то

Отсюда по теореме о дифференцировании оригинала

 

 

pF (р) Ф (Р) =4 / (О ф (0) +

t

f (т) ф' (t

г) dr.

 

 

 

 

6

 

(1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Это —так

называемая формула Дюамеля.

 

 

уравнение

Пусть

требуется

решить

линейное дифференциальное

с постоянными коэффициентами л-го порядка

 

 

 

 

 

L [х] =

 

(0 + ахх'п-" (t) + ... + аях (/) = /(/), «о Ф 0,

(2)

при

нулевых

начальных

условиях

 

 

 

 

 

 

 

 

 

 

 

 

* (0) = х' (0) = ... = х<*~1>(0) = 0.

 

 

(3)

Допустим,

что известно решение

уравнения

 

 

 

 

 

 

 

 

 

 

 

Ц х )

=

\

 

 

 

 

 

(4)

с той же левой

частью и правой

частью, равндй единице,

при

усло­

виях

(3). Переходя к операторным уравнениям, будем иметь {А (р)

известный

многочлен

от р)

 

 

 

 

'

 

 

 

 

 

для

(2) и

 

 

 

 

Л (Р) X{p) = F(p)

 

 

 

 

(5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А (р) X, (р) =

i

 

 

 

 

(6)

для

(4).

Из

(5) находим

X (р) =

^ ^

, а

из

(6)

А(р)

РХ, (Р) ’

откуда X (р) = рХ х (р) F (р)

Согласно формуле

(I)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

(7)

 

 

рХх(Р) F (р) =

/ (0 *i (0) +

/ (т) х[ (I -

т) dr.

 

 

 

 

 

 

 

 

 

 

 

8

 

 

 

 

 

Учитывая,

что

(0) = 0,

получаем

 

 

 

 

 

 

 

 

 

 

 

X (р) = рХ, (р) F (р) =

/ (т)

(/ —т) dr.

 

(8)

 

 

 

 

 

 

 

 

 

б

 

 

 

 

 

 

Отсюда решение *(/) уравнения (2) при нулевых начальных

усло­

виях

(3) буде!

 

иметь

вид

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x(t) = \f(r) x[(l—r)dr,

 

 

 

(9)

 

 

 

 

 

 

 

б

 

 

 

 

 

 

 

 

где хх (t) — решение задачи

(4) —(3).

 

 

 

 

 

уравнение

П р и м е р

 

1. Используя формулу Дюамеля, решить

при заданных

начальных

условиях

 

 

 

 

 

 

 

 

 

 

 

дг" (0 —дг (0 =

- ^

,

* (0) = * '

(0) = 9-

 

 

Р е ш е н и е .

Рассмотрим

вспомогательную

задачу

 

 

 

 

%

 

*; (0 -

(О=

1,

 

(9) = Х[ (о)= о.

 

 

Применяя

операционный метод,

находим

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

откуда

 

 

 

 

 

 

(Р) = Р (Р2— 1)

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(/) — ^ sh т dr = ch t — 1.

 

По формуле (9)

 

 

 

 

 

 

 

 

 

 

х

 

1

 

sh (t-т ) dT=

~l. (e< -te< -l)-f sh t ■In ]+ el

 

 

1 +с*

 

 

 

 

 

 

 

 

2 '

Требование, чтобы начальные условия были нулевыми,

является

несущественным:

простой

заменой

искомой

функции задачу с нену­

левыми

начальными

условиями

можно

свести, к задаче с нулевыми

условиями.

Покажем

это

па примере дифференциального уравнения

2-го порядка.

решение уравнения

 

 

 

Пусть

ищется

 

 

 

 

 

 

 

о0*»(0 +

а д '( 0 + а д ( 0 = ) ( 0 .

(10)

удовлетворяющее ненулевым

начальным

условиям

 

Положим

 

 

 

A (0) = .V0,

х' (0) = х,.

 

 

 

 

1/(0 = -'"(0

V -

 

(П)

Тогда

 

 

 

 

 

 

 

у’ (/) =

* '( 0 - х „

у" (t) = x"(t),

 

 

 

 

 

и уравнение (10)

преобразуется

к виду

 

 

 

где

 

 

 

аау" (0 +

О|‘/' ( 0 + а д

(0 =

/i (0.

 

 

 

 

h (О= / (О -

а д - а д -

ад<-

 

Далее,'

в силу (1 1 )

 

 

 

 

 

 

 

 

 

 

 

 

у (0) = * (0) — *„ = 0,

 

у' (0) = х* (0) — я, = 0.

 

Таким образом, приходим к следующей задаче Коши:

( 0 + « , ; / ' ( / ) + а д ( ' ) = = / , (0 ,

1/( 0) = 0, у' (0) = 0

снулевыми начальными условиями.

Пр и м е р 2. С помощью формулы Дюамеля решить следующую

задачу Коши:

х" + 2 х '+ х =

-И .

( 12)

х(0) = —2,

х' (0) =

1.

(13)

Р е ш е н и е . Сводим задачу

(12) —(13)

к задаче с нулевыми

на­

чальными условиями. Для этого полагаем

1/(0 = * (0 + 2 - / .

Тогда

*'(/) = * '( / ) - 1 , у" (t) = x" (I),

и уравнение (12) преобразуется в следующее:

V* (0 + 2у' (0 + У (0 = YJTpTjF*

где

У(0) = 0, у '( 0) = 0.

Решая последнюю задачу с помощью формулы Дюамеля, найдем

(/=*-' [/-In (1+ОЬ

Решение исходной задачи (12) -(13)

х(0 = в - ' [ / - 1п ( ! + 0 ] - 2 + /.

Спомощью формулы Дюамеля найти решения урав­ нений, удовлетворяющие заданным начальным условиям:

748.

х" -

х' =

—— -,

 

х (0) -

дг' (0) = 0.

 

 

 

(Н -с'Г-

 

 

 

749.

х" +

2х' + х

^ - ^

,

х (0) -

х' (0) = 0.

750.

х " - х ’

е"-<

 

х (0) а-' (0) = 0.

2 -1- с '’

 

 

 

 

 

 

 

751.

х " - х ' =

1

 

A:( 0 ) = A-' (0 ) = 0 .

 

 

 

752.

** + * = ■

1

 

x (0) = x' (0) = 0.

 

 

2 +

cos / ’

 

 

 

75 3 - * ' + J‘ = ? + W

 

A' (0) = x' (0) = 0.

 

 

 

754.

x" + ,v=

 

1

 

x (0) = x' (0) = 0.

 

 

 

 

■COS2 t 9

 

 

755.

*•" + * =

 

1

 

x (0) = x' (0) = 0.

1 +

sin21

9

756.

.v"—* = th/,

 

x (0) = x' (0) = 0.

757.

x'“+ x' =

2 +W r*

A- (0) = x' (0) = x" (0) = 0.

§ 17. Решение систем линейных дифференциальных уравнений операционным методом

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

Пусть, например, нужно решить систему дифференциальных

уравнений

второго

порядка

 

 

'

 

 

 

 

 

п

 

 

 

 

 

 

 

 

 

 

2

(а<* I

F

+

b<xi

r +

-

л (/)

(1 )

 

 

/< =

1

( / =

1. 2,

...,

л),

 

 

 

 

 

 

 

 

где aib> bikt

ciff — const, при

начальных

условиях

 

 

 

 

хк (0)= а * .

4 (0) = Ра.

 

(2)

Обозначая

с

через

Xk (p)

и

Z7/ (р)

изображения

л> (/)’ и /,•(/),

от си­

стемы (I)

учетом (2) перейдем к операторной

системе

 

п

 

 

 

 

 

 

п

 

 

У! (в |* Р 3 +

& /* Р + с/*)

(p) = Fj (р ) -j- 2

[ ( я /* р + &,•*) а *

(3)

* = 1

 

 

 

 

 

 

* = |

 

 

 

 

 

 

(< =

1. 2.........я).

 

 

Решая систему (3) как линейную алгебраическую систему уравнений

относительно

Х к (р),

найдем

X к (р),

а

затем

их

оригиналы xk (/)

(Л =1, 2, ....

я). Эти

последние

будут

решениями

задачи Коши для

системы ( 1).

Решить систему

уравнений

 

 

П р и м е р .

 

 

 

 

(

= 3 (у-х-\-г),

 

 

 

 

{

tj" = x - tj,

 

 

 

 

 

 

I г" — — г,

 

 

 

 

 

 

 

А*(0) ~

х ' (0) — 0,

 

 

 

 

У(0) = 0,

//'(0) = -1 ,

 

 

 

 

2 (0) = I.

г '( 0) = 0.

 

 

Р е ш е н и е . Переходя

к операторной

системе,

получим

 

 

(

р * Х - 3 ( К — X + Z),

 

 

 

 

| р"-К+ 1 = Х - К , '

 

 

 

 

 

{ p ! Z - p = - Z ,

 

 

 

 

Л

(Р) == А- (0 .

и

(Р) - - У (0 ;

2 (Р )

; • г (/).

Решая последнюю систему относительно X (р), Y (р) и Z (р), по­ лучим

* (Р )

3 (р — 1)

Y (р) =

3 ( Р - 1 )

1

Р2 (Рг + 4) ’

P*(P2+ l ) ( P s+ 4 )

Р2+ 1 ’

 

 

2(Р)

Р

Р -+ 1 -

Находя оригиналы для X (р),

У (р), Z (р),

получаем

 

x(t) = ~

( I - 0 - - 4

COS2/+ |

Sin 2/,

\

V (/)■=

(1 0 +

 

cos 2/ — -g- sin 2t — cost,

 

z (/) — cos t.

Решить системы уравнений:

7 5 8 - { , , ' + L O ,

x -4-я? = у -j~ ,

{

f x' — у' — 2x + 2y = 1 —2£,

m { * - + v + * = o .

j

x" - 3x' + 2x + y' - y = 0,

761* \

— x' + x + i/”- 5y' + 4y = 0,

 

x (0) = x' (0) = y' (0) = 0, y {0) = 1.

Ш-{Т=*+2У.

 

( 2x”— x' -f 9x y"— 1/ 3y = 0,

 

763,

\

2x" + x' + 7x -

</" + y' - by = 0,

 

 

 

 

 

x (0) = x' (0) = 1,

y(0)= y '(0 ) = 0.

 

f x' + y' — у = e7,

 

 

" М ь ' + г ч - * - * » » .

 

 

 

 

x' = — X -fy + z + e',

 

 

765.

{

y^

x - y + z + e*,

x (0) = у (0) = г (0) == 0.

 

 

2' =

X + < / + 2 +

4,

 

 

 

 

 

< / - 2 ,

 

 

 

766.

 

y' = — x — z,

x(0) = —1, y(0) = 0, 2(0) == 1.

 

 

2

x

y,

 

 

 

 

 

x'=*y + z,

 

0, у (0) = 1,

z (0) = 1,

 

 

i/'= 3 x +

2 , x (0) =

г>— 3ft -j- y,

х — 3у х,

768.

У' =

7/-М +е°',

* ( 0 ) = 1 ,

У (0) = 1 .

769.

а' =

2х - у

+ г,

 

 

у' —х-\-г,

 

х(0) =

1, т/(0) = 1, г (0) = 0.

г' = Зх - \ - у - 2г,

х' — —2х —2у — 4г,

770.У' = 2x + y — 2z, х (0)= у (0) = г (0) = I.

г' — Ьх+ 2у -\- 7г,

tx' = — л'

У z -\-t,

771. ty' = x - y

+ ? + t3, z (1) = 7/(1) = г(1) = 0.

te' — х - \ - у z -\- 4,

*0-------СЛ'0|

772.Xi = схх-\- СХПу

 

Хцсхп

схп-ч%

 

 

 

 

 

х„(0) = 1,

л:, (0) = х, (0) = . . . = х„ (0) = 0.

773.

За:' +

+ у’ = 1,

х (0) = у (0) = 0.

х

, , ,

 

, 7 *

 

+4у

+3у —0,

'

J ' '

 

3tx' = 2х-\-у — г,

 

 

774.

2ly' = x + 3y + z,

,v(l) = //(l) = 2 (l) = l.

 

Ыг' х

7У ~\r5г,

 

Гх' —x — 2y = t,

776.Электрон вылетает из начала координат с началь­ ной скоростью ©0, направленной по оси Ох. Найти закон

движения электрона, предполагая, что напряжение маг­ нитного поля Н постоянно и направлено перпендикулярно

кплоскости хОу.

777.Снаряд вылетает из орудия со скоростью ©0 м/сек под углом 45° к горизонту. Найти, пренебрегая сопротив­ лением воздуха, наибольшую высоту, на которую подни­ мается снаряд, и место его падения.

778.Электрон движется в магнитном поле постоянного напряжения Н. Найти траекторию, если начальная ско­ рость ©о образует угол а с направлением магнитного поля.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]