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

eP+i

Решение. Изображение п е п равно (еР—е)2. Поэтому согласно теореме об изображении суммы

Л- 1

ер+1

 

 

 

т е т

 

 

 

(iгР — е )2 (еР— 1)

 

 

 

т = О

 

 

 

е Г еР

еР

е — 1

 

 

еР+1

 

^ Т у 2[ е Р - [

~ ё Р ^ е

е

(еР — е)

Ь—6я - (б — 1) /le*]

•’ ( e - U 2[

Следовательно,

л — 1

т е т = е (1 — gn) + (е— 1) пеп

(6- 1)2

ш =о

Найти следующие суммы:

 

л— 1

 

 

 

 

 

л — 1

 

 

 

870.

2

т2 -

 

 

871.

2

tficosmoc.

 

 

ш = 0

 

 

 

 

 

т = 0

 

 

 

 

л — 1

 

 

 

 

 

Л — 1

 

 

 

872.

2

т (т

1 ).

 

873.

2

emcosma (/гэ=2).

 

т=0

 

 

 

 

 

m e 0

 

 

 

Формула ^э браще ния .

Пусть

решетчатая

функция / (/г)

имеет своим изображением функцию Z7* (р) комплексного переменного

p= a + iT,' где

F* (р), в

силу своей периодичности:

Е*(р + 2ш) =

= Z7* (р), рассматривается в основной полосе

— я <

Im р л.

Если

извесшо

изображение

F* (р), то

оригинал /(п) можно

найти по формуле обращения

 

 

 

 

 

 

 

 

 

 

с-f in

 

 

 

 

 

/(«) =

2Ш

J

F*(P)enpdP>

(24>

 

 

 

 

 

с — |’л

 

 

 

где с—любое действительное число, большее, чем

абсцисса сходи­

мости s*.

 

когда

F* (р)

есть

правильная рациональная дробь

В случае,

относительно еР, имеем

 

 

 

 

 

 

 

 

 

/(П) =

2

ге« [f* (р) e'n~v Р],

(25)

 

 

 

 

V

Р \

 

 

 

 

 

где сумма берется по всем полюсам функции F* (р),

расположенным

в полосе — л < 1 тр ^ л ;

и на

ее границе 1шр = д. Если pv—про­

стой полюс, то

 

 

 

 

 

 

 

 

 

res [f * (р) еР |л"11] =

11m

[ f * (р) ( е Р - Л ) е Р 1'1' 1'];

 

pv

 

 

 

р —Ру

 

 

 

если

же pvполюс

порядка r.v,

то

 

res [F*(p) еР<п-^] =

 

 

 

 

PV

 

 

 

 

 

 

=

77— Ц тг

lim

_ n

[f* (p) (eP — ePvY v eP,n~1'].

(27)

 

(rv ~ H1

p-+pv dePVv l)

 

 

 

 

П р и м е р

14.

Используя

формулу (25), найти оригинал,

соот­

ветствующий изображению

 

еР

 

 

 

 

F* (о)=

 

 

 

 

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

 

 

 

 

w

 

е2Р — ЪгР+ 2

 

Р е ш е н и е . Находим нули знаменателя. Имеем е2Р — ЗеР + 2 = 0, откуда еР= 1 и еР= 2. Значит, рх = 0, р2 = In 2 — простые нули зна­

менателя, а следовательно, они являются простыми полюсами функ­ ции F* (р) в основной полосе. Находим вычеты функции F* (р)еР1п~1> относительно этих полюсов. Имеем

 

 

 

- l ;

 

fiP (гР_pin 2\ eP(n-l)

 

 

[F*

(eP — i)'(eP—2)

=

 

 

^ lln2

е1п2я

 

= p !!7n2? T i = >

2_ j

= 2= r = 2я.

Следовательно, по формуле (25) получаем f(n) = - 1+ 2*.

П р и м е р 15. С помощью формулы обращения найти оригинал для функции

Р е ш е н и е . Функция

еР

F* (Р) = :(сР— 1) (еР-\- 1)

имеет два простых

полюса в точках /4 = 0, ра= /я

основной полосы

— л < Im р ^

л. Находим

,.

еР (еР— 1) еРт-1)

 

res

[Z7* (р) еР{п~1)] =

1

р= о

 

р ™о

(е Р -\ )(е Р + \ )

2 1

 

 

 

=

еР*

g/ЛЯ

 

 

 

lim

е'я — 1

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

(25)

 

1

 

 

 

 

/ д о = 4 + ^ ( - 1)л- 1.

У

Найти оригиналы для следующих изображений:

874.

F*(p) = {^-

w .

875.

F* (р)

 

еР

 

 

 

 

 

 

 

е*Р — 1еР-\- ю*

876.

F* (р)=*

 

 

 

877.

F*(p) = ^

т .

878.

 

 

 

еР

 

 

 

 

 

F* (р)=е2Р + 2аеР + Ы> (а > 0)-

 

 

 

879.

F* {р)=

 

 

 

 

 

 

 

 

Р е ш е н и е р а з н о с т н ы х у р а в н е н и й с п о м о щ ь ю

д и с к р е т н о г о п р е о б р а з о в а н и я

Л а п л а с а .

Уравнение

вида

F(n,

f(n),

/( я + 1 ) ........./(я +

А)) = 0

 

 

(28)

или

 

 

F(n,

f(n), Д/ (я),

A4/(«))*=0,

 

 

(29)

 

 

 

где /(л) — искомая

решетчатая

функция,

называется

разностным

уравнением k-eo порядка.

сводится

к уравнению вида (29) при

помощи

Уравнение вида

(28)

формулы.

 

 

 

 

 

 

 

 

 

 

/ ( « + А)= /

(л)+

q д/ (л)+ q дv (л)+ . . . + Д*/ («).

(30)

а при помощи формулы

_

 

 

 

 

 

Д */(« )= /(Л + Л) — q / (п + А— 1)Н-С|/ (rt + A —2) —... + (—I)*/(rt)

 

 

(А =

0, 1,

2, ...),

(31)

где

^ • --JJ1 — т +

биномиальные коэффициенты,

урав-

нение (29) сводится к уравнению (28).

[ (п)

и ее разностей,

Если

уравнение (29) линейно

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

то оно Называется линейным разностным уравнением. Линейное раз­

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

 

порядка

k с постоянными

коэффициентами

имеет

вид

Ь0 Д*/ (л) +

Ьх Д*-1/ (Л) + ... +

buf (Л) = Ф (Л),

 

(32)

 

 

где ф (п) —заданная

решетчатая функция, /(л) —искомая

решетчатая

функция, Ь0, Ьъ ...,

bk— постоянные/.причем

Ь0 Ф 0,

Ьк Ф 0.

k) по

Заменяя в уравнении (32) разности

Дmf (р)

(т= 1, 2,

формуле (31), получим другую форму

разностного

линейного

уравнения:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а0/ (л +

А) +

aJ (я + А — 1)

. . . + akf (л) = Ф (л).

 

(33)

Если

ф (л) =

0,

то

уравнения (32)

и (33)

называются однород-

ними, е^ли же ф (л)=£0, то эти уравнения

называются

неоднородными.

Разностное

уравнение,

содержащее

/(л)

и

 

 

называется

разностным уравнением

k-vo

порядка

(k > 0).

Таким

образом,

урав­

нение (Зз)

при

а0 Ф 0 и ак Ф 0 есть

неоднородное линейное разност­

ное уравнение

&-го

порядка.

 

 

 

 

 

 

 

 

Порядок разностного уравнения может не совпадать с порядк0м наивысшей разности, входящей в него, если разностное уравнец**е записано в виде (32).

П р и м е р 16. Разностное уравнение Д3/ (п) + 4Д2/ (п) + 5Д/ (п) + 2/ (я) = О

после замены в нем разностей по формуле (31) приводится к виру

/(« + 3) + /(« + 2) = 0

или

Пп+ 1) + Нп) = О,

т.е. является разностным уравнением первого порядка.

Определить порядки следующих разностных уравнений:

880.Д4/ (п) + 4Д3/ (п) + 4Д2/ (л) - / (п) = 0.

881.Д3/ (я) + ЗД2/ (л) гЬ ЗД/ (л) + / (л) = л3 + 1 .

882.Д3/ (я) + 2Д2/ (л) + Д/(/г) = 2".

883.Д4/ (л) + 4Д3/ (п) + 6Д2/ (л)Ч- 5Д/(л) + 2/ (n)=sin

Начальные

условия

для

разностного уравнения &-го

порядка

задаются в виде

значений решетчатой

функции f (п) и ее

разностей-

до (/г—- 1)-го

порядка включительно

при

п — 0, если уравнение имеет

форму

(32),

или

в виде

значений'

решетчатой

функции

f (п) при

п = 0,

1, ...,

k — 1, если

уравнение имеет форму (33).

 

Решение

линейных разностных

уравнений с постоянными коэф­

фициентами

методом дискретного

преобразования

Лапласа

произво­

дится

по той

же

схеме,

что

и в

случае

классического преобразова­

ния Лапласа. Применяя

дискретное

преобразование Лапласа к таким

уравнениям и используя теоремы линейности, опережения или изобра­

жения разностей, мы приходим к

простому линейному алгебраиче­

скому

уравнению

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

изображения

F* (р) искомой

функции

/ (п)

(уже с учетом начальных условий). Разрешая это

алгебраическое

уравнение

относительно F* (р),

получим операторное

решение

разностного

уравнения, оригинал для

которого будет иско­

мым решением

исходного

разностного

уравнения,

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

поставленным начальным

условиям.

 

 

 

 

П р и м е р

17.

Найти решение уравнения

 

 

 

 

/ ( " +

1 ) - * / » = ! ,

/ (0) =

0.

(34)

Р е ш е н и е .

Пусть / (/г) -r1-F* (р).

По теореме

опережения

f(n+\)-F-ePF* (р).

Применяя к обеим частям (34) дискретное преобразование Лапласа, получим операторное.уравнение

ePF* (p)-eF* ( р ) = ~ ,

откуда

еР

F * № ( е Р - e ) ( е Р — 1 ) ‘

Функция F* (р)

имеет два

простых полюса р = 0,

р =

1:

res

[Z7* (р) е1п~1)Р\ = Пт

eln~DP = j ^— ,

р= о

 

р-+оеР

е

1

е

res

ff * (р) e,n- l>Р] ——

еп~1 — — г.

 

 

Р= 1

 

е — 1

 

е — 1

 

 

Следовательно, по формуле

(25)

 

 

 

 

Решить следующие линейные однородные разностные

уравнения:

 

 

 

 

 

884.

/ (я +

1) — 2/ (л) = 0,

/ (0) = 1.

 

 

885.

/(л +

2) + 2 /(л + 1 )-Ь /(л )= 0 ,

/(0) = 1,

/(1) = 0.

886.

/ (л +

2) — 2/ (л + 1 ) + 2 /(л) = 0,

/( 0)= 0,

/(1) = 1.

887.

/ (л 4* 3) — 3/(л +

2) + 4/(л-|- 1) — 2/ (л)=0, f (0)=0,

888.

/ ( л +

4 )+ /(«) = 0, /(0) = 0,

/ ( 1) = 0, /( 2) = 1.

 

 

 

 

 

 

/ ( 1 ) = 1 ,

/ ( 2) = 0 ,

/(3 ) == 0 .

889.

f (л + 3) +. 3/ (л +

2) + 3/ (л+ 1) + / (л) = 0,

 

 

 

 

 

/ ( 0) =?0 , / ( 1 ) = 0 , / ( 2 ) = 1 .

Решить следующие линейные неоднородные разностные

уравнения:

 

 

 

 

 

890.

/(л + 1 ) + 2/(л) = л,

/(0) = 0.

 

 

891.

/ ( л

2) — 4 /(л) == 4я,

/(0 )= /(1 ) = 1.

 

892.

/(л +

2 )+ /(л ) =

1- ( —1)», /(0) = 0 , f(l) = l.

893.

/ (л +

2) - 6/ (л +

1) +

9/ (л) = л • 3",

 

 

 

 

 

 

/ (0) = о, / (1 ) = о.

894.

/ (л + 3) + 3/ (л +

2) + 3/ (л + 1) + / (л) = cos ля,

 

 

 

 

/ ( 0) = 0 , / ( 1 ) = 0 , / ( 2) = 0 .

895.

/ (л + 3) - 3/ (л + 2) + 3/ (л +1) -

/ (л) = л2,

/ ( 0) = 0 , / ( 1 ) = 0 , / ( 2) = 0 .

ГЛАВА III

ТЕОРИЯ УСТОЙЧИВОСТИ

§ 22. Понятие об устойчивости решения системы дифференциальных уравнений. Простейшие типы точек покоя

Пусть имеем систему обыкновенных дифференциальных уравнений

 

 

 

 

 

 

 

Уъ У2>•••» Уп)

 

 

(*=1. 2» •••»

л),

(.0

где

(i, А = 1,

2.........л)

существуют

и

непрерывны,

и пусть

 

 

fyk

....

п)

 

есть

решение

этой

системы,

удовлетворяющее при

( i = l , 2,

 

t = t0 условиям

 

Ф;('о) = Ф<

 

 

(‘=

1 .2........«)•

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение

ср/(/)

( i = l , 2,

...,

л) системы

(1) называется устойчи­

вым по Ляпунову

при

что

для

 

если для

любого е >

0 можно

подо­

брать б (е) >

0 такое,

всякого

решения y-L(t) (i = 1 , 2, ..., п)

той

же системы

 

(1),

начальные

значения

 

которого

удовлетворяют

неравенствам

IУ -1 (У -

 

Ф; | <

 

6 (в>

 

( * = 1 . 2 .........п),

 

 

 

 

 

 

 

 

 

 

для

всех

r ^ , t 0 справедливы'неравенства

 

 

 

 

 

 

 

 

 

 

|у«(0 - ф

| ( 0 | < 8

 

 

( / = 1.

2 ,

.... л).

 

(2)

т. е. близкие по начальным

значениям

решения

остаются близкими

для

всех

t

/0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

п) устойчиво,

 

 

Иными словами, решение

ф /(/)

(/= 1 ,

2,

...,

если

Достаточно близкое

к

нему

в

начальный

момент t = t0 решение

</; (*)

(i= 1, 2,

...,

п)

 

для

всех

t ^ t 0

содержится

в сколь

угодно

узкой

е-трубке,

построенной

 

вокруг

решения

ф /(/)

( / = 1 , 2,

... , п).

реше­

ния

Если

при сколь

угодно

малом

б > 0

хотя бы для одного

yi(t)

(< = 1,

2,

...,

п)

неравенства

(2) не выполняются, то

решег

ние

(pi (0

( i = l ,

 

2,

...,

л)

называется неустойчивым.

устойчиво, но.

 

Если решение ф,- (t)

(i =

l, 2, ...,

п)

не

только

кроме того,

удовлетворяет

условиям

 

 

 

 

 

 

 

 

 

 

 

 

lim

|* //(0 -< P i(0 i =

0

 

(i =

l,

2,

...,

л),

(3)

 

 

 

<.— оо

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

если | У; (^о) ~ ф? | <

 

то

решение

ф; (/)

называется асимптотически

устойчивым.

 

устойчивости

 

решения

ф /(<)

системы

(1) может

быт*,

 

Вопрос об

 

вреден к

вопросу

об устойчивости

нулевого

решения

(t) = 0 некфч

торой новой системы уравнений, получающейся из (1) линейной заменой искомых функций

Xi(t)=yi(t) — q>i(t)

(1 = 1, 2,*..., п),

(4)

где Xi(t) — новые неизвестные функции, равные отклонениям прежних неизвестных функций ^ (/) от функций ф /'/), определяющих иссле­ дуемое решение. Поэтому'в дальнейшем будем считать, что на устой­ чивость исследуется именно нулевое решение xt (0 = 0 или, что то же, расположенная в начале координат точка покоя системы уравнений

 

■^- =

'ф/г(0 *i, х2,

,

хп)

( t = l , 2, ..., п),

(5)

Вместо термина

«нулевое

решение» будем употреблять

термин три­

виальное решение.

 

 

 

*/(0 = 0 (/= 1 , 2,

п)

условие

В применении к точке покоя

устойчивости выглядит так:

 

 

 

п) системы (5)

устойчива по

Точка

покоя xt ( 0 = 0

(i =

1,

2,

...,

Ляпунову,

если для каждого е > 0

можно подобрать б (е) >

0 такое,

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

|*/(/0) 1< б(в)

(* = 1,

2, ..., п) следует |* /(/)|< е

(/ = 1, 2,

..., п)

при всех i ^

t0.

 

 

 

 

 

П р и м е р

1.

Каждое решение уравнения

 

 

устойчиво.

Действительно, решени? хх (/) этого уравнения, удовлетворяющее начальному условию x1(tQ) = x{)l, есть хх (/) = х\ const.

Рассмотрим другое решение х2(0 уравнения (6), удовлетворяю­ щее начальному условию

 

 

 

 

x2(io) = 4-

 

 

(7)

Для этих

решений

имеем

| х2( 0 ~ * i (0 1= 1*§—*? I Для

всех

дователыю, для

всякого

е > 0 существует

б > 0, например, 6 = е

такое,

что

как только |хЗ —* ; |< б , то

для

решений

х2(0 и * i (0

будет

выполняться неравенство

 

 

 

 

 

|* 2 (0 —* i (0 1=

1*2 —*! I < е

при

всех

 

Следовательно, любое решение уравнения (6)f устойчиво. Однако асимптотической устойчивости нет:

1*2 ( О ”

*1

( 0 | =

| *§ —

| - т ^ О -

Пр и / - * - + 00 .

П р и м е р

2.

Каждое

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

 

 

 

 

dx .

л

(8)

 

 

 

 

л + л г = 0

 

 

 

 

 

асимптотически

устойчиво.

 

уравнения имеет вид

В самом деле,

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

 

 

 

 

х (/) =

Се~(.

(9)

Решения хх (/),

х2(/)

уравнения

(8), удовлетворяющие начальным

условиям *I (/O) = *P

*2 № ))=*3» сУть

 

 

( 0 =

4 е~

~ 'о).

*2 (0 =

4 е~ и ~ и)‘

220

Отсюда

I Х 2 (0 “~ * 1 (О I = ! Iе~~ /о) ->0 при /-V-J-00,

что означает асимптотическую устойчивость любого решения урав­ нения (8).

П р и м е р

3.

Решение

 

 

 

dx

 

х (t) = — 1 уравнения — = \ — x2 (t) не­

устойчиво,

так

как при

/- > + о о все решения

уравнения

 

 

 

 

 

(1 + х 0) е2 (<

<о) — (1 х0)

 

 

 

 

 

 

(1 +*а)«2(' -

<,) +

(1 -

^

 

стремятся

к + 1 .

Решение

*(<) = 1 этого

уравнения

согласно опре­

делению

асимптотически

устойчиво.

 

 

 

 

Пользуясь определением, исследовать на устойчивость

решения следующих уравнений и систем:

 

896.

 

dx .

 

,

'

x (0) = 1.

 

 

Ш+ х = 1>

 

 

897.

 

я —

< ( * - » ■

лг (0) = 1.

 

 

898.

~ — 2x = t

 

 

 

 

 

 

 

dt

хх

 

 

 

 

 

 

 

899.

dt

zxl>

 

X (0 ) = 0 .

 

 

900.

dx

 

.

 

* ( 0) = 1.

 

 

ш =

cos/>

 

 

 

901.

 

( * - » .

 

x (0) = «/ (0) = 0 .

 

 

 

( а —

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

902.

 

Ш=У>

 

x(0) = tj (0) = 0 .

 

 

k^

=

+ 3x,

 

 

 

 

 

 

 

 

П р о с т е й ш и е т и п ы

т о ч е к

п о к о я .

Пусть

имеем систему

дифференциальных

уравнений

 

 

 

 

 

 

 

 

 

 

d x

 

 

 

 

 

 

 

 

 

 

Ш- Р ( х ,

У),

 

 

М)

 

 

 

 

 

 

dtdy = Q(х, У)-

 

 

 

 

 

 

 

 

 

 

 

Точка (дг0, у0) называется точкой покоя или особой точкой системы (А),

если Я (х0, (/0)==0, Q (х0. у0) = 0.

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

 

2 f ~ an x -\~ai2У’

(1)

du

^ а21х а22У>

 

где Я;/ (/, / = V 2) —постоянные. Точка (0, 0) является точкой

покоя

системы (1). Исследуем расположение траекторий системы (1) в окрест­ ности этой точки. Ищем решение в виде

х = а 1е*‘,

y = a 2ekt.

(2)в

Для определения k получаем характеристическое уравнение

 

ап k

а12

(3)

O21

&22—^

 

Рассмотрим возможные случаи.

I. К о р н и х а р а к т е р и с т и ч е с к о г о у р а в н е н и я д е й ­

с т в и т е л ь н ы

и

р а з л и ч н ы .

Подслучаи:

1)

^ < 0 ,

k2< 0 .

Точка

покоя асимптотически устойчива (устойчивый узел).

2) ^ > 0 ,

&2 > 0 .

Точка покоя неустойчива (неустойчивый узел).

3) А:1 > 0 ,

k2< 0.

Точка

покоя

неустойчива

(седло). 4) ^ = 0,

/?2 > 0 . Точка

покоя

неустойчива.

5)

ki = 0,

62 <

0.

Точка

покоя

устойчива, но

не асимптотически.

 

 

 

 

 

 

 

у р а в н е н и я к о м ­

II.

К о р н и х а р а к т е р и с т и ч е с к о г о

п л е к с н ы е :

ki = p-\-qi,

k2 = p —qi.

Подслучаи:

1)

p < 0 , q ^ Q .

Точка

покоя

асимптотически устойчива (устойчивый фокус).

2) р > 0,

<7=5^ 0.

Точка покоя неустойчива (неустойчивый фокус).

3) р = 0,

q Ф 0.

Точка покоя устойчива (центр). Асимптотической

устойчиво­

сти нет.

 

 

к р а т н ы е :

kx= k2.

Подслучаи:

 

1)

ki = k2< 0 .

III. К о р н и

 

Точка

покоя

асимптотически устойчива

(устойчивый

узел). 2) kx =

= k2>

0.

Точка

покоя

неустойчива

(неустойчивый

 

узел).

3) k\ =

= k2 —0.

Точка покоя

неустойчива.

Возможен

исключительный слу­

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

коэффициентами

dx

=

п

(< =

1, 2........ п)

(4)

- ^

2 ‘ ''А '

 

 

/=»

 

 

 

 

характеристически^

уравнением

будет

 

 

 

an — k

а \2

а13

а\п

 

 

°21

 

Q22 k

°23

а 2П

= 0.

(5)

ап\

 

°П 2

апЗ

йпп — Ь

 

 

I) Если действительные части всех корней характеристического уравнения (5) системы (4) отрицательны, то точка шпокоя Xj(t) = Q (i'= l, 2, п) асимптотически устойчива.

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