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

книги из ГПНТБ / Абгарян К.А. Матричные и асимптотические методы в теории линейных систем

.pdf
Скачиваний:
97
Добавлен:
25.10.2023
Размер:
17.28 Mб
Скачать

3 4 0

К В А Д Р А Т И Ч Н Ы Е И Э Р М И Т О В Ы Ф ОРМ Ы

[ГЛ . X I I I

Покажем сначала, что

 

 

 

 

 

сц = сц

(і,/ =

1,2,

. . . , /г).

(4.4)

Согласно

(4.1)

 

 

 

 

A gj = сііё 1 +

c2jg 2 +

•••

+ c kig k,

 

A gi cuSi +

c2ig 2-f •••

+ c kig k.

 

Умножим первое равенство на вектор g lt а второе на вектор gj. Получим, учитывая, что векторы g lt g 2, ..., g k нор­ мированы и попарно ортогональны:

Ctl = (gt>Agj),

cu = {g ,,A g l).

Но, так как А — симметрический оператор,

(gp Agt) = (Ag,, gi) = (gi} Agj)

и

cll Cij,

что и доказывает соотношение (4.4).

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

Покажем, что алгебраическое уравнение (4.3) имеет толь­ ко вещественные корни.

Допустим противное, а именно, пусть уравнение (4.3) имеет комплексный корень ^0. Тогда среди корней этого

уравнения имеется и корень Ä.0, комплексно сопряженный корню Х0.

Систему (4.2) для удобства последующего изложения

представим в матричной записи:

 

 

Сх = Хх

(С = (сц)).

(4.5)

(

х°Л

 

X — Д,0)

Пусть х0 =

. . . — решение системы (4.5) при

V х і )

 

 

так что

 

и —

(4.6)

 

 

$ 4] Л И Н Е Й Н Ы Е О П Е Р А Т О Р Ы В Е В К Л И Д О В О М П Р О С Т Р А Н С Т В Е 341

Тогда решение системы (4.5) при X = Х0 представится столбцовой матрицей

Умножим (4.6) слева на транспонированную матрицу

Хо:

XQ С х 0 = ^ оХо-'-О-

Здесь х0х0 — вещественное число. Вещественным является также скаляр х0Схо, так как

х0Сх0= (х0Сх0У — хоСх0 = ХоСхй

(С — С). Но тогда Х0, как отношение двух вещественных

чисел

х0Сх0 и х0х0, не может быть комплексным числом.

Это

противоречие показывает, что корни алгебраиче­

ского уравнения (4.3) вещественны. Пусть Х0— корень этого уравнения. При X = Х0 система (4.2) имеет ненулевое реше­

ние х°, х°, ..., х°к, которое к тому же, будучи решением алгебраической системы с вещественными коэффициентами, вещественно. Лемма доказана.

Л е м м а 4.2. Пусть А симметрический оператор в п-мерном евклидовом пространстве R, а е его собствен­ ный вектор. Тогда совокупность Rx векторов х, ортогональ­ ных к е, есть (п — 1 )-мерное подпространство, инвариант­ ное относительно оператора А.

Д о к а з а т е л ь с т в о леммы 4.2 подобно доказатель­ ству леммы 3.3.

С помощью приведенных лемм легко устанавливаются и следующие теоремы (аналогичные теоремам 3.1 и 3.2).

Т е о р е м а 4.1. В п-мерном евклидовом пространстве R существует п попарно ортогональных собственных векто­ ров симметрического оператора А . Соответствующие собст­ венные значения вещественны.

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

342

К В А Д Р А Т И Ч Н Ы Е И Э Р М И Т О В Ы Ф ОРМ Ы

[ГЛ . X II I

Т е о р е м а 4.2. Для того чтобы оператор А в евкли­ довом пространстве R был симметрическим, необходимо и достаточно, чтобы существовал ортогональный базис, в котором матрица оператора диагональна и вещественна.

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

4.2.

Ортогональный

оператор. Линейный

оператор О

называется

ортогональным, если

 

 

0 0 ’ =

О'О = Е.

(4.7)

Из (4.7)

следует: ] О |2 =

1, т. е.

 

 

| 0

| = ± 1.

 

Если I О I = 1, то О называется ортогональным операто­ ром первого рода; если же | 0 | = —1, — ортогональным оператором второго рода.

Свойства ортогональных операторов в евклидовом про­ странстве аналогичны свойствам унитарных операторов в

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

(см. §

3). Приведем некоторые

из них.

 

 

Ортогональный оператор сохраняет скалярное произ­

ведение векторов:

 

 

(Ох, Оу) =

( X , у)

( X , у е /?)•

Ортогональный оператор сохраняет длину векторов:

(Ох,

Ох) =

(х, х).

Если О — матрица ортогонального оператора в ортонор- • мированном базисе, то

00' = О’О = Е.

(4.8)

Матрица О, обладающая свойством (4.8), называется ортогональной матрицей. Значит, ортогональному операто­ ру в ортонормированном базисе отвечает ортогональная матрица.

Для того чтобы оператор О был ортогональным, необхо­ димо и достаточно, чтобы он переводил какой-нибудь орто­ нормированный базис в R снова в ортонормированный базис.

§ 5] К В А Д Р А Т И Ч Н Ы Е Ф О Р М Ы 3 4 3

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

4.3.

Преобразование симметрической матрицы к диаго­

нальному виду с помощью ортогональной матрицы. Пусть

А — симметрическая матрица порядка я. Будем рассматри­

вать А как матрицу симметрического оператора в ортонор-

мированном базисе

$ =

(g1 g 2 ... g n)

я-мерного евклидо­

ва пространства R,

так

что

 

 

 

А& = &А.

(4.9)

Согласно теореме 4.2 в пространстве R существует орто­ нормированный базис g = (ех е2 ...еп), в котором матрица оператора А диагональна и вещественна. Обозначим эту матрицу через Л. Тогда

Л § = §Л.

(4.10)

Далее, существует ортогональная матрица О, которая преобразует ортонормированный базис g в ортонормирован­ ный базис

# = gO.

(4.11)

Подставляя (4.11) в (4.9), получаем А% = %ОАО~1. Сравнивая последнее соотношение с (4.10), находим

Л= ОАО~1= ОАО'.

Всоответствии с этим имеем также

А= СГ'ЛО = О'ЛО.

Из подобия матриц Л и Л следует, что диагональными элементами матрицы Л служат собственные значения матри­ цы А.

§ 5. Квадратичные формы

Пусть А — симметрический оператор в я-мерном евк­ лидовом пространстве R, а А — матрица этого оператора в некотором ортонормированном базисе g = {еу е2 ... еп). Матрица симметрического оператора в ортонормированном базисе симметрична = А ’).

3 4 4

К В А Д Р А Т И Ч Н Ы Е И Э Р М И Т О В Ы

Ф ОРМ Ы

[ГЛ.

X I I I

 

Скалярное произведение векторов А х

и у , где х ,у

£ /?>

удовлетворяет тождеству

[Ах, у) = (х, Ау).

Учитывая, что

= gx, у -- gу,

где хм у — столбцовые матрицы, составленные из коорди­ нат векторов X и у соответственно, получаем

СА х , 31) = (Л gx, gy) = ( 21 Лв,х,,

 

21 од Л

=

 

 

\ /=1

 

 

ft=»l

/

 

 

 

 

 

 

 

П

 

 

 

 

 

=

S

 

(Ле,, **)*/&■

Обозначим

 

 

 

 

/,*=і

 

 

 

 

 

 

, n).

(Л ßi,

£fe) = а,-*

(t, А =

1, 2, . . .

Тогда

 

П

 

 

 

 

 

 

 

 

 

 

 

 

 

(Ах, у) =

21

а/*дг,%

 

 

(5.1)

и, в частности,

 

i,k=\

 

 

 

 

 

 

л

 

 

 

 

 

 

 

 

 

 

 

 

 

(Ах, X)

21

aikX[Xk.

 

 

(5.2)

 

 

I.Ä=1

 

 

 

 

 

Так как А — симметрический

оператор

 

в евклидовом

пространстве,

то

 

 

 

 

 

 

аік = (Ае(, ek) = (е{, Aek) = (Aek, et) = aki.

Квадратичной формой называется однородный много­ член второй степени над полем вещественных чисел отно­ сительно переменных xlt х2, хп. Любую квадратичную форму от п переменных хі можно представить в виде

П

21 aikXcxk, i,k=1

где aik = akl (i, k — 1, 2, n), a в соответствии с (5.2) эту квадратичную форму можно рассматривать как скаляр­ ное произведение векторов А х и х , где Л — некоторый сим­ метрический оператор в n-мерном евклидовом пространстве.

П

Билинейную форму 21 aikxilJk можно трактовать как скалярное произведение векторов Л л: и у.

§ 5]

К В А Д Р А Т И Ч Н Ы Е Ф ОРМ Ы

3 4 5

А

Формы (5.1) и (5.2), которые принято обозначать через

(X, у) и А

{х, х), коротко записываются в виде

 

 

А(х, у) = х'Ау = у'Ах

(А = (а^)),

(5.3)

 

 

А (х, х) = х'Ах.

(5.4)

 

Матрица

А является

матрицей

оператора А

в базисе

g

= (ег е2... еп), ибо, как

нетрудно проверить,

 

АЪ = ѢА.

Определитель матрицы А — det А — называется диск­ риминантом квадратичной формы А (х, х). Если дискри­ минант равен нулю, то форма называется сингулярной.

5.1. Замена переменных. В формах (5.3) и (5.4) произ­ ведем замену переменных:

 

 

х = ТЬ,

у=Тѵ\,

 

 

(5.5)

где

£ и

т) — столбцовые матрицы, составленные

из ко­

ординат

..., £п и гх, г2,

...,

г)асоответственно. Получим

где

 

А (х, у) = 1'Ах\,

 

А (х, х) =

£'

 

 

 

 

 

 

 

 

 

 

 

А =

TA T .

 

 

(5.6)

 

Формула (5.6) связывает

матрицу

А =

(aik)

формы

Л

 

 

 

 

п

 

 

2

ackltlk с матрицей первоначальной формы

2 aikxtxk-

i,k=1

 

 

 

i.k= 1

 

. Матрица А преобразованной формы также является симметрической матрицей.

Если Т — невырожденная матрица, то ранг матрицы формы при замене переменных (5.5), как это видно из (5.6), не меняется.

Две симметрические матрицы А и А, связанные друг с другом равенством (5.6), в котором det Т ф 0, называются конгруэнтными. С каждой квадратичной формой связан целый класс конгруэнтных матриц. Все матрицы данного класса имеют один и тот же ранг, которому равен по опре­ делению ранг соответствующей формы.

5.2. Закон инерции. Допустим, что квадратичная форма

А {х, х) = х'Ах

(5.7)

3 4 6

К В А Д Р А Т И Ч Н Ы Е И Э Р М И Т О В Ы

ФОРМЫ

[ГЛ. X II I

каким-нибудь способом приведена к виду

 

 

А(х, х ) =

/■=!

 

(5.8)

где %і Ф

0 = 1, 2......г) и

 

 

 

 

6і = £ <*«**

(£=1,2,

.. •. г)

(5.9)

 

А=I

 

 

 

— независимые линейные формы от переменных хѵ х2, .... хп

(элементов столбцовой матрицы х). Обозначим

dt = {dndi2... din).

Тогда замену переменных (5.9) молено записать так:

I = Dx.

Перейдем к матричной записи и

в соотношении

(5.8):

Г

 

 

 

А(х, х) = 2 ^ix'didix =

x'D'ADx.

(5.10)

i=i

 

 

 

Вычитая из (5.10) равенство (5.7), получаем

 

х' (D'AD — А) X =

0.

 

Но последнее равенство при любых значениях хѵ

х2, ...

..., х„ справедливо лишь тогда, когда

 

А =

D'AD.

 

(5.11)

Ранг матрицы D типа г х

п равен г (матрица D набрана

из г линейно независимых строк dlt d2, ..., dr). Ранг диаго­ нальной матрицы А типа г х г также равен г (так как все диагональные элементы отличны от нуля). В этих усло­ виях ранг матрицы D'AD равен числу г. Отсюда в силу ра­ венства (5.11) может быть сформулирован следующий вывод: число квадратов в представлении (5.8) всегда равно рангу формы. Более того, имеет место

Т е о р е м а 5.1

( з а к о н и н е р ц и и

к в а д р а ­

т и ч н ы х ф о р м ) .

При представлении

вещественной

§ 5]

КВАДРАТИ ЧНЫ Е ФОРМЫ

3 4 7

квадратичной формы А (х, х) в виде суммы квадратов:

 

Г

 

 

А(х, х)= 2 ^

-

 

і=і

 

где Х[ Ф 0 (г = 1, 2,

г), а £х, 12,

линейно неза­

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

Д о к а з а т е л ь с т в о .

 

Пусть имеет место и

другое

представление формы А

 

{х, х)

в виде суммы квадратов:

 

 

 

 

А (х, х) =

 

Г

 

 

 

 

 

 

 

 

 

 

2

ИД1?.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

і=Д

 

п

 

 

 

 

где р,( ф О

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,

2,

 

г),

 

 

а

TJ,- =

2

С;А — независи-

 

 

 

 

 

 

 

 

 

 

 

 

 

*=і

 

 

 

мые линейные формы от переменных хѵ х2, ..., хп.

 

Пусть

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Хг >

О,

Х2^>0, . . . ,

ХЛ> 0 ,

 

< 0 ...........\г <с0,

| Д > 0 ,

р2> 0 ,

. . . ,

(Лг> 0 ,

 

р,/+ 1<

О, . . .

, рг < 0.

Предположим,

что h Ф

I,

 

например,

h <

/.

удовлетво­

Переменным л'х, х2, ..., хп дадим

значения,

ряющие системе г — (/ — /г) уравнений

 

 

 

 

h — 0

(г = 1,

2, . . . .

 

/г),

 

 

ту = 0

(/ = /

+

1.......... г)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.12)

и не обращающие в нуль хотя бы

одну

из

форм

lh+2 i

Ъ-

Такие

значения

 

существуют, так

как

в про­

тивном

случае

из

 

 

=

0,

 

..., Ъ,г =

0

следовало

бы, что

все г уравнений

£(- = 0

=

 

1,

..., г)

являются следствием

г — (I h)

уравнений

 

(5.12),

 

но

это

невозможно

в силу

линейной независимости форм

 

£v

g2,

Ir- При таким обра­

зом выбранных значениях хѵ х2, .... хп в тождестве

 

 

 

 

 

У

i t ?

=

 

у

 

 

 

 

 

 

 

 

 

 

 

IЛ;ы

 

~ і

 

 

 

 

 

левая часть

равна

i=i

 

 

 

 

i=i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 K é < o ,

А-Н

3 4 8

К В А Д Р А Т И Ч Н Ы Е

И Э Р М И Т О В Ы Ф ОРМ Ы

[ГЛ. X I I I

а правая равна

 

 

 

I

 

 

 

2

іѵі* > о.

 

 

*=і

 

 

Предположение h Ф I привело к противоречию. Теоре­ ма доказана.

5.3. Приведение квадратичной формы к главным осям. Рассмотрим произвольную вещественную квадратичную форму

П

 

А (X, х) = 2 aikXiXy = х'Ах

= А').

/.*=1

 

Так как матрица квадратичной формы А является ве­ щественной симметрической, то существует такая вещест­ венная ортогональная матрица О, что

0~ ‘Л0 = А —- diag(^1, Хй..........

Хп).

Здесь Хѵ Хг, ..., Хп — собственные значения матрицы А. Учитывая, что О-1 = О', при замене переменных

X 0^

(5.13)

получаем следующее выражение для квадратичной формы:

А(х, *) = Г А £ = S U ? -

(5-14)

і=і

 

Таким образом, справедлива

форма

Т е о р е м а 5.2. Вещественная квадратичная

А (х, х) ~ х’Ах всегда может быть приведена посредством ортогонального преобразования координат к канонической форме (5.14), где Хх, Хг, ..., Хп собственные значения мат­ рицы А.

Рассмотрим гиперповерхность второго порядка, задан­ ную с помощью квадратичной формы:

П

 

х'Ах = 2 üikXiXb = с

(с — const Ф 0). (5.15)

i.k= 1

 

При ортогональном преобразовании координат (5.13) уравнение (5.15) принимает вид

1,

§ 5]

КВ АДРАТИ Ч Н Ы Е

ФОРМЫ

349

где

 

 

 

~~Г = ~Т ,

8і = sign

(&1= 0 при

Xt = 0).

Будем рассматривать матрицу А как матрицу некоторо­ го симметрического оператора А в некотором ортонорми­

рованием базисе $ = (ffi g 2 ••• gn) евклидова

пространства

R; при этом хѵ х2, ..., хп представляют собой координаты

вектора х в базисе $ . Тогда Л = diag (Xlt ...,

Хп) есть мат­

рица

оператора

А в новом ортонормированном базисе g =

= (ег

е2 ■■■еп),

а £lt £2,

\п — координаты

вектора х в

базисе g. Векторы базиса определяют направления осей координат в пространстве R. Поворот осей координат опре­ деляется ортогональным преобразованием

g = $0-

Новые оси являются осями симметрии центральной по­ верхности (5.15). Оси симметрии поверхности обычно назы­ ваются главными осями.

В связи с этим приведение квадратичной формы А (х,х) посредством ортогонального преобразования к канониче­ ской форме (5.14) называется приведением к главным осям.

Из (5.14) следует, что ранг г формы А (х, х) равен числу не равных нулю собственных значений матрицы А.

Вещественная квадратичная форма А (х, х) называется

неотрицательной (положительно определенной), если при любых значениях переменных А (х, х) >• 0 (А (х, х) > 0, X =5^ 0). Аналогичным образом определяются неположитель­ ные (отрицательно определенные) квадратичные формы.

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

Наконец, соотношение (5.14) позволяет получить сле­ дующие важные неравенства. Обозначим через А,тах и Я,тіп

соответственно

максимальное и минимальное собственные

значения матрицы А квадратичной формы А (х, х).

Тогда

п

п

 

 

2 М < 2 Х а * Й = Я т . * 2 &

 

f = l

/ = I

I

 

n

n 2

Ä-mln 2 I £/•

/2- 1

> 2/-1 ?ЧтіІп£?=

Соседние файлы в папке книги из ГПНТБ