Файл: Медич Д. Статистически оптимальные линейные оценки и управление.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
В квадратичной форме без ограничения общности можно считать матрицу А симметрической. В этом лег ко убедиться, замечая, что для іф\ члены в разложении формы имеют вид Если матрица А несим метрическая, ее можно заменить матрицей Ä, диаго нальные элементы которой совпадают с диагональными элементами А, а остальные имеют вид
При этом значение квадратичной формы не изменит ся, т. е. х' Ах = х'Ах.
Например, если
А : - 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