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

книги / Матричные и асимптотические методы в теории линейных систем

..pdf
Скачиваний:
1
Добавлен:
19.11.2023
Размер:
26.59 Mб
Скачать

Формула (4.17) определяет все субматрицы блочной мат­ рицы Q$ 1, кроме одной — QaaИз вышеизложенного ясно,

что в качестве

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

рица типа ka X 1,

имеющая

производные по т всех по­

рядков. В частности, можно принять

 

= 0.

 

 

Таким образом,

матрица QE^

определяется

полностью.

Через эту матрицу последний столбец матрицы

/С?]

выра­

жается так:

 

 

 

 

К(&кК

 

 

 

 

 

 

 

£ $ =

 

 

 

 

Остальные

столбцы матрицы

Яо'

<т.......Йо)

опре-

деляются соотношениями

(4.13).

 

 

 

 

 

Остается указать способ построения членов разложений

матриц Ма (т, е) и R (т, е),

обращающих

равенство (4.9)

в тождество.

 

 

 

 

 

 

 

(4.10) выполня­

Как показано в гл. VIII, § 2, равенство

ется тождественно

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

е,

если квадратную мат­

рицу

 

 

 

 

 

 

 

 

 

 

 

 

 

/И =

V BkMw (т)

 

 

 

 

 

 

 

А=0

 

 

 

 

 

 

и члены разложения

матрицы

R (т, е) определить

форму­

лами

 

 

 

 

 

 

 

 

 

 

М™= М,

 

 

= 2

KlnMlk~l\

 

 

 

 

 

 

 

1

 

 

 

(4.18)

 

 

 

Rk=

- V

i ; а ,Rt- t.

 

 

= V

,

 

 

Здесь

 

 

 

 

 

 

1=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

(К!'1

 

KP),

 

 

 

Mw = I

Полученные соотношения позволяют последовательно определить члены разложения (4.7), посредством которых представляется формальное решение системы (4.3) в форме (4.5) — (4.6). Тем самым теорема доказана.

 

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

 

 

 

 

Т е о р е м а

4.2. Пусть на сегменте 10, L] а) матрицы

A k (т),

Bk (т) (k — 0, 1 ,2 ,...) имеют производные по т всех

порядков, а А0 (т),

кроме того, является невырожденной

матрицей;

б)

собственные значения

матрицы

U (т)

=

=

АГ 1 (т) В0(т)

разбиты на р групп

^ 0), ...,

=

 

 

 

р

 

п) при условии г (4.4); в) соответствую­

=

1,

р;

2&<х =

 

 

 

щих

 

подпространства R lt /?2, ...,

R p явля­

щие этим группам

ются инвариантными и циклическими подпространствами п-мерного пространства R относительно линейного опе­ ратора U, которому внекотором базисеотвечаетматрица U. Тогда формальное решение системы (4.3) может быть представлено равенствами

X =

2

(Т, е) уа,

(4.19)

 

0=1

 

 

-^ Г “ =

л а (т,

е) уа + М а (т, е) R (т, е) f (/, т, е),

(4.20)

где Ко, Ла, Ма, R — матрицы типа соответственно п хka,

ka х ka, ka x n,

n x n, представленные формальными ря­

дами

 

 

 

 

 

со

 

 

Ко (т, е ) =

 

 

 

 

 

 

 

2

 

 

] (т),

Л„ (т, е) = 2

Л \М

(т),

 

А=0

 

 

 

 

 

 

(4.21)

 

 

 

 

 

 

 

 

Мо {Ъ е) =

S еАЛ ^ ] (т),

Я (т,

е) =

(*)»

 

fc=0

 

 

 

 

 

 

 

причем Л™ есть матрица типа (4.2), а

 

 

 

«la-

— 0&2а^

 

^AfT—lo

1

 

 

 

 

'[£

 

 

ASP =

 

О

 

о

 

о

 

о

 

 

 

 

 

 

 

 

 

 

о

 

о

 

о

 

о

Члены рядов (4.21) в данном случае

определяются сле­

дующими

рекуррентными соотношениями:

 

 

 

 

 

a /а

(/

1» • • •»

^<т)»

 

где

+

Т" U&ka—la “Г

,

= (Ло” 'а 0 Ло° 2Д0

#ст)?

причем

 

/■—Otlo

Л . = Л™ =

1

 

О

а

 

 

Qa

'&2а

0

О

• • • |

СС/гст_1а

a *o<A

0

0

 

i

0

)

p)t

о

 

 

1

 

(Заметим попутно, что здесь det

= (—I)2

•)

Столбцы матрицы

определяются, так.

 

Первый столбец

 

 

 

где

Й? = KQln] = s2= l к,<№,

Й

? = Фа ' (A s) M s4 ft~ 13

(S * О ),

a Qo3 — произвольная матрица типа ka х I, имеющая про­ изводные по т всех порядков. Остальные столбцы матрицы

К1а^ определяются формулами

 

 

Й? = и& Ъ + «

М

+

 

а =

2, з , ..., ka).

Что касается членов

разложений

матриц Мо и /?, то

для них остаются в силе соотношения

(4.18), имея в виду,

что /(М построены по соответствующим формулам.

§ 5. Приближенное решение системы

 

Вектор (столбцовую

матрицу) хт(/, е), определенный

равенствами (4.5) и (4.6) (или(4.19)и (4.20)) в предположении, что в разложениях (4.7) (соответственно (4.21)) оставлены лишь члены порядка не выше т относительно е, назовем приближенным решением системы (4.3). Итак, приближен-

ное решение представляется равенствами

*т(*,е) = 2 /С Г(*. г ) ^ т),

где

т

т

З а м е ч а н и е . Для построения приближенного реше­ ния условия дифференцируемости матриц Av>BVyсформули­ рованные в теоремах 4.1 и 4.2, могут быть ослаблены: для формального построения приближенного решения хт до­ статочно существования лишь первых т — v производных матриц Л v и Вч (v С т).

При условии, что матрицы Av и fiv (v < т) имеют на ГО, L] производные по т до — v -f- 1)-го порядка вклю­ чительно, a f (/, т, е) — непрерывная вектор-функция, ре­ гулярная относительно е в окрестности точки е = 0, име­ ют место следующие оценки для приближенного решения хт(см. Приложение).

Если х (0) =

хт (0),

то существуют

такие ех >- 0 и по­

стоянные ст>

0

е2 >

0 (е2 £ (0, ej),

что для

всех t £

Uv h] cz [0,

L/e2] ||л —

( e <e

2).

Если, помимо сделанных выше предположений, все соб­ ственные значения эрмитовой матрицы 1/2 (А + А*) не­ положительны, то

В случае однородной системы (fs O ) имеет место оценка

1 % | ‘С

Последняя оценка остается в силе и для приближенного решения однородной системы, представленного в форме (3.30), (3.31).

Приведенные оценки свидетельствуют об асимптотиче­ ском характере построенных приближенных решений.

Рис. 10.1.

Г л а в а X

ДИНАМИЧЕСКИЕ ХАРАКТЕРИСТИКИ ЛИНЕЙНЫХ СИСТЕМ

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

Настоящая глава посвя­ щается такого рода характеристикам многомерной линей­

ной системы (рис.

10.1), процессы в которой описываются

уравнением

 

 

Л $ - § -

=B(()x+H(t)u, det Л =£ 0,

(0.1)

где х — матрица выходных сигналов xlt х2, ..., хп (столб­ цовая матрица с размерами п х I); и — матрица входных сигналов и1г и21 •••» Щ (столбцовая матрица с размерами / X 1); А, В, Н — матрицы динамических коэффициентов системы с размерами п х п , п х п н п х 1 соответственно.

§1. Единичная ступенчатая функция

идельта-функция

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

является симметричная единичная функция

 

'

0

при

г <

О,

 

 

l(z)

=

1/2 при

г =

0,

 

(1. 1)

 

L 1

при

г > 0 .

 

 

График этой функции

изображен на рис. 10.2.

Дельта-функция Дирака, или

симметричная единичная

импульсная функция

б (z)

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

переменной z,

определяется условиями

 

 

 

 

 

 

 

 

 

0

 

 

при

t £ [a, b]t

 

 

- r f i t + O)

 

при

t = а,

 

 

- Y H t - o )

 

при t = b,

- r i f v - w + f v + w

при

t £ (а, Ь),

 

 

 

 

 

 

 

( 1.2)

где a <^b, a f(z) — произвольная функция, являющаяся в окрестности точки z — t функцией ограниченной вариации.

Для произвольных функций f (г), непрерывных в точке z — t, в частности, имеем

 

0

при

t £

[а, b],

 

а

-^-/(/)

при

t =а и t =b,

(1.3)

f(t)

при

/ £ (д, Ь).

 

 

 

Полагая / (г) в

1, из условий (1.3)

получаем

следую­

щие свойства дельта-функции:

 

 

 

 

6 <г) =

0 (г^ О ),

J 8 ( i ) d | = l ,

(1.4)

которые часто принимаются в качестве определения дель­ та-функции.

Среди функций, понимаемых в обычном смысле, нет функций со свойствами (1.4), б (г) является символической (обобщенной) функцией, с помощью которой функциональ­ ное преобразование / © - * - / (/) формально можно пред­

ставить как интегральное преобразование. Формальное применение дельта-функции приводит к удобным построе­ ниям, позволяющим получить обобщения многих матема­ тических соотношений, ко­ торые, однако, вообще го­ воря, нуждаются в стро­ гих обоснованиях. Мы здесь не будем касаться вопросов обоснования при­ меняемых далее операций с использованием дельта­ функции, отсылая читателя к монографиям, в которых более детально рассмотре­

ны применения дельта-функции Дирака в теории линейных

систем (см.,

например, [17, 37, 46]).

 

 

Производные б' (.г),

6" (г),

б ю (z),...

дельта-функ­

ции определяются условиями

 

 

 

ь

 

 

 

 

 

] / © 8 w (!-/)< *£ =

 

 

 

 

 

0

 

 

при

16 [а, »].

<

- i)r 4 - f(o ('- и ))

 

при

t = a,

 

 

 

 

при

(1.5)

 

 

 

 

t — b,

( - О' 4 -

(< -

0) +

(/ + 0)1

при

16 (а, 5).

где а ■< Ь, а

/ (г) — произвольная

функция, производ­

ная /w (г) которой имеет односторонние пределы /<Г) (t — 0)

и /<Г>(/ + 0). Соотношения (1.5) могут быть

получены пу­

тем r-кратного формального

интегрирования

по

частям

с учетом (1.2) и (1.4).

 

 

производные

/<г> (z)

Для произвольных функций / (z),

которых в точке г =

t непрерывны, в частности, имеем

 

 

0

при

t \

[а, Ь],

1 /© S W< 6 - * ) 4 =

(— 1)"

/ (Л) (0

при

Ь

а и

] = б,

 

(— 1)7И (it)

при

t € (а, Ь).

 

 

 

 

 

 

( 1.6)

Согласно приведенным соотношениям в случае совпаде­ ния одного из пределов интегрирования с моментом дейст­ вия дельта-функции t перед функцией f<r) (t) (г = 0, 1 ,2 ,...) в правых частях соотношений появляется коэффициент 1/2. Это является следствием того, что как дельта-функция, так и ее производные обладают симметрией относительно мо­ мента t. Однако в большинстве случаев, связанных с прак­ тическим применением, коэффициент 1/2 опускают, пред­ полагая такое расширение пределов интегрирования, при котором импульс целиком оказывается в интервале интегри­

рования. Следуя

этому, мы также будем считать, что

ь

 

 

ь

 

 

f f (I)6м (l- 0

®

=

]' m 8м ( I - f ) d i =

( -

1)7"’ (t)

t

 

 

t—0

 

 

И

 

 

H-o

 

 

 

 

 

 

 

j f ( i ) (l - 1) dl

=

i / (I) s'0 (I - 0 dl -

( -

1)7"’ (t).

Связь между дельта-функцией 6 (z) и единичной функ­ цией 1 (г) представляется символическим соотношением

6 (г) =

ЛЦг)

 

dz

которое легко устанавливается, например, с помощью пре­

образования Лапласа:

оо

« ( < - ! ) = L~'L [f>(t - £)1 = L~' j 8 (t — I) e~p,dt =

= L~' [e_s<l] = l r ' { p L { 1 (t - 0]) = L~'L | Г (/ - 0] =

П р и м е ч а н и е . Все соотношения настоящего пара­ графа сохраняют свой внешний вид и в том случае, когда вместо скалярной функции f (г) стоит прямоугольная мат­ рица, элементами которой служат скалярные функции с соответствующими свойствами.

§ 2. Реакция системы на входной сигнал в виде дельта-функции. Импульсная переходная функция

Реакцию предварительно невозбужденной системы на входной сигнал в виде дельта-функции принято называть

импульсной переходной функцией (иногда эту реакцию называют весовой функцией). Допустим, что на /-й вход пред­

варительно невозбужденной системы подается сигнал в ви­ де дельта-функции 5 (t — |). На г-м выходе появится выход­ ной сигнал — импульсная переходная функция, которую обозначим через gtj (t, £). Сигнал в виде дельта-функции,

поданный

на ;-й вход, вызовет

на разных выходах систе­

мы, вообще говоря, разные сигналы g^

(t, £), g2j (t, £), .

..., S n t

(t>I)- С другой

стороны,

сигналы в

виде дельта­

функции,

поданные

на разные

входы,

вызовут на одном

и том же

(например,

£-м) выходе

разные, выходные

сиг­

налы

gtl

(t, l), gi2

(t,

£) , ...,gtl (t, l).

Многомерная

ли­

нейная система с

/

входами и я выходами характеризует­

ся nl

импульсными переходными функциями

 

 

 

giliUl)

(i = 1,2,

я;

/ =

1,2,

.. ,/) •

 

Эти функции удобно собрать

в

одну

матрицу G (I,

Q с

размерами п х I:

 

 

 

 

 

 

 

 

 

 

 

 

£ll(*. S) § L 2 & l)

gu V, l)

 

 

 

 

 

§21(^» i) §22(^» i)

§21(С I)

 

Столбцы матрицы G (t, £) обозначим через gj (/, |):

Каждый такой столбец представляет собой набор им­ пульсных переходных функций по всем выходам системы, отвечающий какому-нибудь ее входу.

Выходные сигналы не могут появиться раньше, чем бу­ дет приложен входной сигнал, поэтому в реальных систе­ мах

gu(ty| ) S B O при t < t

Это свойство реальных систем принято называть усло­ вием физической осуществимости или физической реализуе­ мости системы.

Для определения матрицы G (t, £) уравнение (0.1) за­ меним эквивалентным соотношением

t

x = X ( t) c + ^ X ( t) Х~' (?) А~' (?) Н (?) и d?,

где X (t) — фундаментальная

матрица однородного

урав­

нения

 

 

л (< )-£ -=

*<*)*•

(2.1)

Если система до подачи входного сигнала находилась в покое, так что х (t0) = 0, то с = 0 и в случае предваритель­ но невозбужденной системы имеем

г

 

x = $ X ( t) X ~ l (?) А~‘ (?) Н (?) и d?

(2.2)

*0

 

Пусть на /-й вход предварительно невозбужденной си­ стемы подан сигнал в виде дельта-функции, т. е.

д5 = 0

(s + j),

uf =

b ( t - l )

&Z[t0ft]).

Тогда, имея в виду, что

 

 

 

Н = (h-Jiz

hj)t

 

и используя

(2.2), получаем

 

 

 

t

 

 

 

g i ( t . I) = U (0

Х ~ ' ( П

( ? ) h , ( ? ) S ( ? - D d ? ,

или

io

 

 

 

g, (t, |) =

X (t) X~' © A -1(Q h, (l).

 

В соответствии с этим

G(t, l) = X (0 X-' (i) A-' (l) H © .

(2.3)

Согласно вышеизложенному. G (/, £) представляет собой решение матричного дифференциального уравнения

A(t)-%- = B(t)G + H (t)6 (t - Q

(2.4)

при условии

с (£ — °> £) = 0.

Соседние файлы в папке книги