![](/user_photo/_userpic.png)
книги / Матричные и асимптотические методы в теории линейных систем
..pdfвнешних возмущений; б — матрица (типа / х 1 ) управля ющих воздействий; т0, тх, г, /0, 1г — квадратные матрицы порядка п, а а — матрица типа п X /.
Результаты, полученные в § 2, непосредственно приме нимы к уравнению (3.1). Для этого нужно лишь в формулах (2.26) векторную функцию ф0 заменить функцией аЬ -{- q>, Ф2 положить равной 0 и в выражениях коэффициентов при нять е = 1 .
Этим путем, удерживая в уравнениях члены, содержа щие е в степени не выше первой, получим преобразованную
систему уравнений регулируемого объекта в виде |
|
|||||
|
+ |
« У - ^ Г |
+ (“ » + |
“ 2 » ) = I1»(6t0) + |
(3.2) |
|
|
|
|
(о = 1 , 2 , |
. . . , |
р), |
|
ГДе |
|
|
|
|
|
|
а |
~ |
Ра |
-(- 2 ^ |
j 4" |
аа —^1аа&а> |
|
СС^а |
—- — |АоЯ!о |
(ttl^Ha^iG— |
-j- Ota^aa ^2ao®a» |
|
bm — trio' (a6 -f- ф),
= — (x^V + т 1щ) rn ' (аб -h <p) — *i‘W iT 1 (a 4 f + 4 ?) •
В случае, |
когда векторная функция аб мала, полагая |
Ф0 = Ф и = |
аб, будем иметь |
+ а»1 |
+ («о + “W ) га = |
= Ра [Ь0— (4 ПР + mV'mL) 1ЩХф] — xS^p/no-14 г (а = 1 , 2 , . . . . р).
Преобразованные уравнения объекта управления имеют более простой вид, поскольку левая часть системы оказы вается расщепленной и связь между подсистемами полной системы сохраняется только через матрицу б в правой час ти системы. Такая система более удобна при всевозможных исследованиях, в частности при моделировании на анало говых машинах.
П р и м е ч а н и е . При замене переменных желатель но, чтобы норма вектора za совпадала бы или же мало
отличалась от нормы вектора q0. Этого можно достичь так. Пусть
<7о = ех[с - J - -f («а га (3.4)
— вектор, который при е = 1 представляет приближенное решение системы (3.1), соответствующее решению га систе мы (3.2) или (3.3). Квадрат нормы этого вектора равен
tJoCfa — ZoX(jX(jZo |
8 |
|
|
X<j2c + ZoXaxfi1 |
dza |
|
|
|
~dt + |
||||
|
|
“f" |
* |
• дч |
■}" ZQXJO K(jZ(j 1 -j- 8' |
|
Матрица x0 всегда может быть выбрана так, что |
||||||
|
|
ХаИо = |
Ека. |
|
(3.5) |
|
Используя произвол, имеющийся в выборе матриц |
||||||
можно добиться |
выполнения |
равенств |
|
|||
|
нах\а = 0 |
|
(t = |
1,2). |
(3.6) |
|
Действительно, |
|
|
|
|
|
|
|
|
|
|
$ 3 |
|
|
x jx [i3 = Ха (XiX2 |
Хр) |
|
—Ха2 |
Xs<7/Sa, |
||
|
|
|
|
|
5 = 1 |
|
_ ^1'ро _
или, учитывая (3.5), получим
XcxSi1= Й + Ха 2 Х ^З.
S*<J
Отсюда ясно, что равенства (3.6) будут иметь место, если принять
------ Ха 2 |
Xs^[so- |
s^o |
|
Тогда |
|
ЯоЦо = ZQZC |
82 . . . , |
т. е. с точностью до членов, содержащих е2, норма векторов при преобразовании (3.4) сохраняется.
Г л а в а XIII
КВАДРАТИЧНЫЕ И ЭРМИТОВЫ ФОРМЫ
§ 1. Метризация векторного пространства
Рассмотрим векторное пространство R над полем комп лексных чисел ЗС, в котором введена бинарная операция, ставящая каждой паре *, у векторов из R в соответствие скаляр (*, у) — скалярное произведение* и у , причем для любых векторов * , y t z из R и любого комплексного числа а из di
О |
(лг, У) = (У* х) (эрмитова |
симметрия), |
|
2) |
(а*, у) = а (*, у) |
(ассоциативный закон), |
|
3) |
{х + у уz) = (*, z) + |
(у , z) |
(дистрибутивный закон). |
Если для векторов пространства R определена бинар ная операция со свойствами 1) — 3), то говорят, что в R
внесена эрмитова метрит.
Из 2) и 3) следует еще:
2') |
(х,ау) = |
а (х ,у), |
3') |
(х, У + Z) = |
(*, у) + (*, Z). |
Если для любого вектора * из R |
||
4) |
(*, *) > О, |
то эрмитова метрика называется неотрицательной. При выполнении же дополнительного условия
5) |
(* ,* ) > 0 |
(* Ф 0) |
эрмитова метрика называется положительно определенной. Векторное пространство с положительно определенной эрмитовой метрикой называется унитарным пространством.
Под длиной вектора * подразумевают арифметическое (неотрицательное) значение корня квадратного из скаляр ного произведения (х , х):
\х\ = У (.х , *).
Из 2) и 5) следует, что каждый вектор унитарного про странства, отличный от нуля, имеет положительную длину, а длина нулевого вектора равна нулю. Вектор х называется нормированным, если | х \ = 1. Любой вектор х можно про нормировать, умножив его на какое-нибудь комплексное
число, модуль которого равен -Д -р
Два вектора х н у называются ортогональными (*_|_ J_y), если (*, у ) = 0 .
В случае численного пространства, элементами которого служат столбцовые матрицы, скалярное произведение опре деляется так:
(о, Ь) = Ь*а.
По аналогии с этим, наряду с общепринятым обозна чением скалярного произведения векторов пространства R — {х, у), мы будем пользоваться и другим условным обо
значением—^ * , полагая, что
(х,у) = у хх.
|
Здесь у х |
выступает в роли вектора, «сопряженного» |
||||||
вектору у. |
|
|
|
|
|
|
||
|
Свойства 1)—5) при этом представляются так: |
|||||||
1) |
|
у хх = |
(*,з>) = |
(у, *) = х ху , |
|
|||
2) |
з>х |
( а * ) |
= |
( а * , у) = |
а ( * , у) = |
ау х * , |
|
|
2') |
{ау)хх = |
(*, ау) = а (*, у) = |
аухх, |
|
||||
3) |
г х (х |
у) = (* + у у z) — |
(*, г) + |
(Уу z) = |
z xx -f z xy, |
|||
30 |
(У + |
z)xx |
= (x,y-\-z) = (x, у) + (*, z) = |
31х* -f z xx t |
||||
4) |
|
x xx = |
(*, *) > |
0, |
|
|
|
|
5) |
|
x xx — (*, *) > |
0 |
(* Ф 0). |
|
|||
|
Длина вектора * равна |
|
|
|
||||
|
|
|
|
|* | |
= |
Xх*, |
|
|
так что скалярное произведение (1.3) можно записать и в форме
П_
у хх |
= 2 hikxiyk |
(hik = {ес, ek) = |
exet). |
(1.4) |
||
|
i,k=i |
|
|
|
|
|
В частности, имеем |
|
П |
|
|
|
|
|
х хх = |
hikXtxk. |
|
|
||
|
2 |
|
(1.5) |
|||
|
|
|
l,k=1 |
|
|
|
На основании свойства 1) эрмитовой метрики |
|
|||||
|
Hik ~ |
@k |
== |
Gfc=s hki. |
|
(1 .6) |
П |
|
hik = hki (i,k = 1 , |
|
..., n), |
||
Форма 2 |
hikXiXk, |
где |
2, |
|||
l |
,6=1 |
|
|
|
|
|
называется эрмитовой формой; выражение, стоящее в правой части равенства (1.4), называется билинейной эр митовой формой. Таким образом, квадрат длины вектора унитарного пространства представляется в виде эрмитовой формы его координат.
Матрица Я эрмитовой (билинейной эрмитовой) формы в силу (1 .6) является эрмитово сопряженной (эрмитовой), т. е. Я* = Я.
Если базис %в R составлен из ортонормированных век
торов |
ev eit ..., еп>т. е. таких, что |
|
||
|
0 |
при |
i Ф к, |
|
|
1 |
при |
i — к |
(t, к = 1 , 2 , . . . , л), |
то в |
этом случае Я — диагональная |
(единичная) матрица |
||
и формы (1.4), (1.5) принимают вид |
|
|||
|
|
|
П |
__ |
|
у хх = {хУу) = 2 |
|
||
|
|
|
/=1X |
|
|
X х X = |
(Ху х) = 2 |
I xtI2. |
|
|
|
|
( = \ |
|
В следующем параграфе будет показано, что ортонормированный базис существует в каждом унитарном простран стве.
Если в качестве числового поля ЭСпринято поле вещест венных чисел, то метрика, удовлетворяющая условиям 1) — 5), называется евклидовой.
Векторное пространство /? над полем вещественных, чичел с положительно определенной евклидовой метрикой на
зывается евклидовым пространством.
В евклидовом пространстве скалярное произведение век
торов хм^у с |
координатами соответственно |
х( и у( (i = 1 , |
|||||
2 , |
/у представляется |
равенством |
п |
|
|
||
|
Ч |
|
|
y 'g x %x = |
|
|
|
|
^(х>у) = у хх = |
23 |
|
|
|||
|
|
|
|
|
i,k=1 |
|
|
Здесь у' — строчная матрица, полученная из у |
транспониро |
||||||
ванием, S(k = skl (t, k = |
1 , |
2 , .... n) — вещественные |
числа. |
||||
В |
частности, |
|
|
П |
|
|
|
|
|
|
|
|
|
|
|
|
(X, лг) = X хX = |
= |
23 SikXiXk. |
|
|||
|
|
п |
|
|
l,k=1 |
|
|
|
|
|
|
|
|
|
|
Выражение |
23 sikxixk называется квадратичной формой |
||||||
|
|
i,k=1 |
|
|
|
П |
|
|
|
|
|
|
|
|
|
относительно хг, х2, |
хп. Квадратичной форме 23 |
slhxtxk |
|||||
|
|
|
|
|
|
ltk=l |
|
П
отвечает билинейная форма 23 si i,k=1
Квадрат длины вектора евклидова пространства пред ставляется в виде квадратичной формы его координат.
З а м е ч а н и е . Из положительной определенности мет рик, введенных в унитарном и евклидовом пространствах, вытекает положительная определенность соответствующих эрмитовых и квадратичных форм, т. е.
x*gxg x > 0 |
и x 'g xg ;t> 0 |
(знаки равенств имеют место только при х = 0).
§2 . Ортонормированные базисы в унитарном
иевклидовом пространствах
Здесь мы покажем, что как в унитарном, так и в евкли довом пространстве имеется ортогональный базис, т. е. та
кая система векторов |
е2, |
.... еп (п — размерность прост |
|||
ранства), что |
I = |
0 |
при |
i Ф kt |
|
|
(2. 1) |
||||
(е{, ek) = \ |
Л |
при |
. . |
||
' |
1Ф 0 |
i = k. |
|
Заметим, что векторы, обладающие свойством (2.1), линей но независимы. В самом деле, равенство
CtiBi -}- |
+ |
"f" |
~ |
® |
|
возможно только тогда, когда |
а± — а2 = |
|
= ос„ = |
0 , так |
|
как по умножении этого |
равенства скалярно, например, |
||||
на ек получаем в силу (2 .1) ak(ek, ek) = |
0 , |
откуда а л = 0 . |
|||
Доказательством существования ортогонального |
базиса |
в рассматриваемых пространствах может служить приво димый ниже процесс ортогонализации*), который позволяет из данных п линейно независимых векторов построить п
попарно ортогональных векторов. |
|
|
|
||
Пусть & = |
(gv g 2, ..., gn) — какой-нибудь базис в про |
||||
странстве (унитарном или евклидовом). |
eit ...» еп, |
||||
Построим систему |
ненулевых |
векторов ev |
|||
удовлетворяющих условию (2.1). Положим |
|
||||
- |
gv |
|
|
|
|
et = g2 + Yia^i. |
|
|
(2 .2) |
||
|
• |
• |
• |
• |
|
е п = |
ё п + |
Y l/ie l + |
~Ь |
Уп—\пвп—\- |
, |
Условие (2.1) приводит к следующим выражениям для чи
сел |
уц: |
|
|
|
|
|
УН |
(gh *() |
(i= 1 , 2 , |
, п — 1; j = |
2, 3, . . . , |
ti). |
|
(et, et) |
||||||
|
|
|
|
|
|
(2.3) |
|
Соотношения (2.2), (2.3) позволяют последовательно |
|||||
строить взаимно ортогональные векторыelt е2, .... |
что |
|||||
и подтверждает существование в R ортогонального базиса. |
||||||
|
Равенства (2 .2 ) можно записать в виде матричного соот |
|||||
ношения |
g = & + |
g r, |
|
|
|
|
где |
|
|
|
|
||
|
|
Уп |
YlS • • • |
Ут |
|
|
|
|
0 |
|
|||
|
g — (^1^8 • |
0 |
0 |
Таз |
У2п |
|
|
* |
|
* |
• |
|
|
|
|
0 |
0 |
0 |
Уп—In |
|
|
|
0 |
0 |
0 |
0 |
|
*) Этот процесс ортогонализации в иных обозначениях приведен
и в § 5 гл. VIII.
(так как BLy — скаляр). Но, с другой стороны, так как g ортонормированный базис,
П
ВкУ= ( 2 в# е0 ek) = (By, ek). i—I
Учитывая последнее соотношение, из (3.2) непосредст венно приходим к (3.1).
Теперь легко установить, что оператор В, определенный соотношением
B y = 'L (y ,A e k)e„ |
(3.3) |
k=[ |
|
является сопряженным по отношению к оператору Л . Дей ствительно, если А * — оператор, сопряженный оператору А, то для произвольных х и у
(х, By) = (х, 2 (у, А ek) ед = |
(л:, 2 |
(А*у, ek) ek) = |
|
к=я1 |
ft=l |
|
|
откуда В = А*. |
|
- |
(*, А*у), |
|
|
|
|
Равенство (3.3) определяет оператор, сопряженный опе |
|||
ратору А , единственным образом. |
|
|
|
Найдем матрицу сопряженного |
оператора Л* |
в орто |
|
нормированием базисе g. |
|
А в базисе g. Через В |
|
Пусть А — матрица оператора |
|||
обозначим матрицу сопряженного оператора Л* в |
том же |
||
базисе. Имеем |
|
|
|
Л g = gЛ, |
Л*g = gЯ. |
|
Отсюда, так как g — ортонормированный базис и потому g x g = Eni находим
л « 8 хл jg. в = 8 хл *8-
Далее,
( |
е Л |
(е?А*е1 |
е?А*ея |
В = 1 |
|
\(А*е1А*ег ...А*е„)= \ |
|
\ |
ея ) |
\etA * el |
е*А*еа |
Так как Л и Л* — взаимно сопряженные операторы, e*A*ej = (A*eJt е{) — (е/г Aet) = (Aet, et) = efA e{,