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

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

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

внешних возмущений; б — матрица (типа / х 1 ) управля­ ющих воздействий; т0, тх, г, /0, — квадратные матрицы порядка п, а а — матрица типа п 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).

 

 

Длина вектора * равна

 

 

 

 

 

 

 

|* |

=

*,

 

 

S 1]

а условие ортогональности векторов х и у принимает вид

у хх —X ху = 0.

Обозначая звездочкой переход к эрмитово сопряженной величине, отметим еще такое свойство:

(ухх)* — у хх = Xху.

Пусть R — л-мерное унитарное пространство и g =

*= {ех е2...еп) — какой-нибудь

базис в нем.

Обозначим

через х столбцовую матрицу,

составленную из

координат

хх, х2, ..., хп вектора х , а через у — столбцовую матрицу,

составленную из координат ух, уй,

...» уп вектора у в этом

базисе. Тогда

у =

%у.

(1.1)

х =

В силу свойств 1) — 4)

 

 

 

 

УХХ = { Ы хх = (Jje,y}j

х

= ^Jj№?jх ,

 

или

 

 

 

(1.2)

у хх — у*&хх,

 

где у* = (уху2 ...уп) — строчная

 

матрица, эрмитово

со­

пряженная столбцовой матрице у,

а

 

 

&хлг =

Учитывая еще первое соотношение (1.1), из (1.2) находим

у хх = у*Н'х.

(1.3)

Здесь

Н' = g x g =

е х е х

е х ег

е х еп

(«i. e j

(е2, е х)

(бл, <?i)

е х е х

е $ е й

е х еп

(fii* еа)

(^2*

(бл, <?а)

 

 

 

 

 

 

ех еа „

-(^1*

 

(®а» ^л)

•••

J $ e x е х е%

 

&п)

 

{fin*еа) _

 

 

 

так что скалярное произведение (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.

Отсюда

§ = # ( £ „ - Г)"'

(Е„ — единичная матрица порядка п).

От ортогонального базиса g можно перейти к ортонормированному базису путем умножения $ справа на диагональную матрицу

diaS ( V{ev ex) * Viez,ег) ’ ’ ‘

V(en,en) \ '

Таким образом, в унитарном

(евклидовом) простран­

стве существует и ортогональный

и

ортонормированный

базис.

 

 

§3. Линейные операторы в унитарном пространстве

3.1.Сопряженный оператор. Если для любых векторов

хи у из унитарного пространства между линейными опе­ раторами А и В в R выполняется соотношение

(Ах, у) = (xt By),

то оператор В называется сопряженным по отношению к оператору А.

Оператор, сопряженный оператору А, обозначается че­ рез А*.

Каждому линейному оператору А отвечает единствен­ ный сопряженный А*. Покажем это.

Выберем в R ортогональный базис $ = (е1 е2...еп). Если В — линейный оператор в R, то для произвольного

вектора у из R справедливо соотношение

П

 

By = 2 . (By, ek) ек.

(3.1)

ft=i

 

Всамом деле, имея в виду, что В — матрица оператора

Вв базисе g, Вь Вг, .... Вп— строки матрицы В, а у — столб­ цовая матрица координат уг, у2, ..., уп вектора у в том же базисе, получаем

By = Bgt/ = %Ву = s екВку = 2 В„уек (3.2)

fes=»l

fc = a l

(так как 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{,

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