книги / Матричные и асимптотические методы в теории линейных систем
..pdfФормула (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» • • •» |
^<т)» |
|
ное решение представляется равенствами
*т(*,е) = 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).
Приведенные оценки свидетельствуют об асимптотиче ском характере построенных приближенных решений.
Г л а в а 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.