Файл: Медич Д. Статистически оптимальные линейные оценки и управление.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 11.04.2024

Просмотров: 228

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

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

Собственные

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

векторы

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

нонической

жордановой

форме. Эти результаты приме­

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

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

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

 

 

(А—ХІ)х

= 0

 

 

имеет

нетривиальное решение. Здесь

А — матрица раз­

мера

пХп, а X — п-вектор. Эта задача

называется

зада­

чей отыскания собственных

значений;

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

значения

X называются

собственными

значениями

мат­

рицы А, а

векторы, для которых справедливо соотноше­

ние Ах = Хх, называются

собственными

векторами

матри­

цы А.

 

 

 

 

 

 

Необходимое и достаточное условие того, что задача отыскания собственных значений имеет нетривиальное решение, выражается уравнением

\А—ХІ\ =0 .

 

 

Это алгебраическое уравнение п-го

порядка

назы­

вается характеристическим

уравнением

матрицы А. Его

решениями являются п (не обязательно

различных)

соб­

ственных значений матрицы A: Хі,...,Хп

(они также

обо­

значаются через Хі (А), і =

1,.. .,п].

 

 

30


Можно показать, что

[A\=f]b{A)

і=.і

 

SpA 4 S « « =

2 ^ ( i 4 ) .

1=1

i=\

Последняя величина, SpA,

называемая следом мат­

рицы А, является суммой элементов, расположенных вдоль главной диагонали матрицы А (знаком Л обозна­ чается равенство по определению). Легко доказываются два свойства следа матрицы:

Sp (А+В)

=Sp Л + S p ß ;

Sp (ABC) =Sp

{ВСA) =Sip

(CAB).

Если матрица А—симметрическая,

то можно пока­

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

Наконец,

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

матрицы

А (не

обязательно

симметрической)

различны,

a M — матрица,

составленная

из собственных

векторов

матрицы А,

то

 

 

X,

о

.

.

. 0

 

 

 

 

М-1АЖ=

0

 

 

• . о

 

 

 

 

 

о

о

.

.

. х п

,

 

 

где собственные векторы матрицы А в матрице M упо­ рядочены так же, как соответствующие им собственные значения %и . . . , %п. Если матрица А, кроме того, симме­ трическая, то M является ортогональной матрицей и

 

 

X,

о

.

. .

о

М'АМ

= О

Х2

.

.

. (Г

 

 

О

О

.

.

. X.

Квадратичные

формы

 

 

 

 

Если А — матрица

размера

 

пХп, а х—«-вектор, то

легко показать, что

пп

х'Ах — S S OijXiXj.

1=1 /=1

Это произведение, являющееся скаляром, называется

квадратичной формой.

31


(ctij-\-aji)xtXj.

В квадратичной форме без ограничения общности можно считать матрицу А симметрической. В этом лег­ ко убедиться, замечая, что для іф\ члены в разложении формы имеют вид Если матрица А несим­ метрическая, ее можно заменить матрицей Ä, диаго­ нальные элементы которой совпадают с диагональными элементами А, а остальные имеют вид

При этом значение квадратичной формы не изменит­ ся, т. е. х' Ах = х'Ах.

Например, если

А : - I о И

— 2 3

ТО

хАх, —• х^ 2x1^Cg " I - '^-^'2

Сдругой стороны, для

Ä-

ап

 

2

« 1 2 +

i7j

Д 2 2

 

 

 

 

имеем:

 

x

' ä x

х\

— 2.г1

л:2 -\- Ъх\-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

Если х'Ах>0

для всех хфО,

то говорят,

что квадра­

тичная

форма положительно

определена.

В этом

случае

обычно

матрицу

А называют положительно определен­

ной. Если

х'Ах^О

для всех х^О, то квадратичная фор­

ма и матрица А

называются неотрицательно

определен­

ными

или

положительно

полуопределеннщми.

Если

х'Ах<0

для всех

хфО, то

используют

термин

отрица­

тельно

определенная

квадратичная форма

(матрица).

Можно

показать,

что

необходимым

и

достаточным

условием положительной определенности матрицы А является условие, чтобы все главные миноры матрицы

«12

А\

«22

 

были положительны. Отсюда видно, что если матрица А положительно определена, то И | > 0 и матрица А не­ сингулярна.

32


 

Теперь

пусть

В— несингулярная

матрица

размера

пХп,

X и у — я-векторы и у = Вх. Ясно, что

у'у=х'В'Вх.

Но

у'у^О,

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

выполняется

тогда

и толь­

ко

тогда,

когда

г/ = 0. Однако

из равенства

г/ = 0

следует

х = 0,

так

как матрица А несингулярна. Следовательно,

матрица В'В положительно определена.

 

 

 

С другой стороны, предположим,

что А является сим­

метрической положительно определенной матрицей. Тог­ да существует ортогональная матрица М, для которой матрица

 

 

А*

= М'АМ =

\\а*и\\

 

 

 

является

диагональной

при

а*ц>0

для всех

і.

Полагая

 

 

 

В* = | | у ^ | | ,

 

 

 

заметим,

что матрица

В* — симметрическая

и

несингу­

лярная,

причем

В*В*=А*.

Обозначая В = В*М\

имеем:

А = М(М'АМ)М'

 

= МА*М' = МВ*В*М'

=

В'В.

Иными словами, для каждой положительно опреде­

ленной

матрицы

А

существует

несингулярная

матри­

ца В, для которой А = В'В.

Матрицу В можно

рассма­

тривать как «квадратный корень» из матрицы А.

(Обыч­

но под

квадратным

корнем

из матрицы

А

понимается

матрица В такая, что В2=А, т. е., В = МВ*М'. В даль­

нейшем автор сам пользуется именно таким определе­ нием (Прим. ред.).

Наконец, так как матрица В несингулярна,.

Л- і = 5 - і ( 5 ' ) - * ,

аэто означает, что если матрица А положительно опре­ делена, то матрица, обратная ей, также положительно

определена. • Заметим, что когда в дальнейшем будут упоминаться

положительно, неотрицательно или отрицательно опре­

деленные матрицы, везде будут подразумеваться

матри­

цы, которые также являются симметрическими.

 

Градиент

 

 

 

Пусть f(x)

определяет

скалярную функцию

«-век­

т о р а X, т. е. f(x)

=f(xi, . . . ,

хп). Определим градиент V *

как оператор, имеющий вид

векгора-строки

 

3 - 85


Тогда

Vxf(x)

является

д-мерным

вектором-строкой

вида

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i L

 

" дхп

 

 

 

в предположении,

что

указанные

частные производные

существуют.

 

 

 

 

 

 

 

 

 

 

Если

f(x)

билинейная

форма

f(x)=y'Bx,

где

В —

матрица

размера

пХп,

не обязательно

симметрическая,

а X я у — «-векторы, то

 

 

 

 

 

 

 

w w = v

* ^

2

2 f c « < J

=

і=і

і=і

=

У'В.

 

 

 

 

 

 

 

 

 

Аналогично, если /(х) —х'Ву,

то

 

 

 

 

 

 

 

Ѵ ж /(х)=г/'Л' .

 

 

 

Наконец,

если

f(x)=x'Ax,

причем матрица

А не обя­

зательно

симметрическая, то

 

 

 

 

 

 

Vxf (X)

= V * ( 2

2

 

)

=

 

2 ß-iiXi

•••2

Яіп-Хі

+

 

 

 

 

 

 

 

 

1=1

i=l

 

 

 

+

2 a i Л

• • • 2 anjx5

=

x'A +

x'A'.

 

 

Если же матрица

А симметрическая, то, очевидно,

Ѵх(х'Ах)=2х'А.

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

2-2. НЕПРЕРЫВНЫЕ ЛИНЕЙНЫЕ СИСТЕМЫ

Уравнения системы

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

£ = F(t) x + G(t)w(t)+C(t)u

(t)

(2-1 )

14