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

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

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

х €

R в базисе (е,, е2,

е

 

п),

а у,, yv

ут — координаты век­

тора у €

S в базисе (g,, g2,

g

m). Тогда

 

 

 

 

 

 

 

 

 

х = &х,

 

 

(2.6.3)

 

 

 

 

 

 

 

У = &У,

 

 

(2.6.4)

где

=

(ер е2, ...) вп),

^ =

(вр 8 2 » • ••»

 

 

 

 

 

 

 

 

hхгi

 

 

(У,)

 

 

 

 

 

X =

.

У =

Уг

 

 

 

 

 

• ««

9 9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

х п

 

 

*■4

 

 

 

 

 

 

 

 

п!

 

 

 

 

Подставим (2.6.3)

и (2.6.4) в (2.6.1):

 

 

 

 

&у = А8?х = (Ае,, Ае2, ..., Ае„)х.

(2.6.5)

Векторы

Aej, Ае2, ...» Ае„

 

принадлежат пространству

S, поэтому

их можно представить через базисные векторы gp g2,

..., gm:

 

Ае* = au gj + alkg2+ ... + amk%m (&=1,2, ...,/*),

или

 

 

 

 

 

 

 

 

 

 

 

 

Aek = &ak (k = 1, 2,... ,«),

 

(2.6.6)

где

 

 

 

 

 

 

 

ulк

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ak ~~ a2k

 

 

 

 

 

 

 

 

 

 

 

lmk

 

 

 

Подставляя (2.6.6) в (2.6.5), получим

 

 

 

или

 

&у = (?ах

З а г

 

 

^ а п)х = ^(щ аг

ап)х ,

 

 

 

 

 

9 у — &Ах,

 

 

(2.6.7)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/ «и

 

 

 

 

 

 

 

 

 

 

 

«12

«1л

 

 

 

А = ( а ,

а2

 

 

««)

«21

«22

«2л

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

«ml

«m2

а тп

 

Строчная матрица 9 набрана из линейно независимых векторов gj, g2, ...» gm. Поэтому из (2.6.7) вытекает равенство

у = Ах.

(2.6.8)

Таким образом, линейному оператору А при выбранных базисах в R и S отвечает некоторая матрица А, которая является матрицей

линейного преобразования координат исходного вектора х в коор­ динаты преобразованного вектора у = Ах. И обратно, т х л-матри- ца А при выбранных базисах Ши & соответственно в «-мерном про­ странстве R и /«-мерном пространстве S представляет некоторый оператор А, который каждому вектору х = е R соотносит неко­ торый вектор у = Ах = S y G S. Связь оператора А с соответствую­ щей матрицей А при выбранных базисах £ и & в пространствах R и S представляется равенствами (2.6.6), которые можно более ком­ пактно записать в виде

А& = SA .

(2 .6 .9 )

Пусть оператору А в базисах & и & в пространствах R и S от­ вечает матрица А. Выясним, как изменится матрица оператора А при изменении базисов. Наряду с % и & рассмотрим новые базисы и в R и S соответственно, связанные со старыми базисами

соотношениями

(2.6. 10)

Здесь Т и N — невырожденные матрицы порядков соответственно п и т . Пусть в базисах 8? и & оператору А отвечает матрица А , а в базисах %х и & х — матрица Av так что

АГ =

АШх — S XAV

(2.6.11)

Используя (2.6.10), из второго равенства (2.6.11) находим

A S = S N A {T~X.

 

Сравнивая полученное

соотношение с первым

равенством

(2 .6 .1 1 ), находим

 

 

 

A = N A {T -1.

(2 .6 .1 2 )

Таким образом, один и тот же линейный оператор А, отобража­ ющий R в S, в зависимости от выбора базисов в R и S представля­ ется разными матрицами, общий вид которых дается формулой

(2.6.12).

§ 2.7. Матрица как линейный оператор в численных пространствах

Пусть А — оператор, который каждому вектору х из л-мернош векторного пространства R соотносит вектор у из m-мерного век­ торного пространства S:

у = Ах,

(2.7.1)

и пусть & и 3? — соответственно базисы в R и S. Если х — столб­

цовая матрица координат х 1г х2, х

п вектора х в базисе S’, а у

столбцовая матрица координат у{, у2, ..., ут вектора у в базисе

то (см. § 2.6)

 

у = А х ,

(2.7.2)

где А — матрица оператора А при выбранных базисах в R и S. Введем теперь в рассмотрение n-мерное численное простран­

ство R, изоморфное пространству R, и m-мерное численное

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

S, изоморфное пространству

S. Каждому вектору

х из R с координатами хр х2,

хп в

базисе Ш поставим в

соответствие

вектор

 

 

 

(х,\

 

 

 

х =

R

 

 

х„.

 

 

 

п)

 

 

и каждому вектору у из S с координатами yv у2>..., ут в базисе & поставим в соответствие вектор

' л '

s.

У = Уг

гт

 

В силу введенного взаимно однозначного соответствия между векторными пространствами R и S, с одной стороны, и векторными численными пространствами R и S, с другой, оператору А, отобра­ жающему пространство R в S, соответствует матрица А линейного преобразования (2.7.2), которое каждому векторуиз R соотносит вектор у из S.

Таким образом, т X «-матрица А выступает как линейный опе­ ратор, отображающий «-мерное численное пространство R в т- мерное численное пространство S.

Матрицу-оператор А можно рассматривать как упорядоченную систему m-мерных векторов — столбцовых матриц

а

II

( а

\

 

 

 

а , \

“ н

 

 

 

f a ‘2 l

in

а 21

 

 

 

а 2п

»

/ )

fl22

 

а 2

»

» а п ~

mlj

 

 

Л 12;

тп)

 

 

 

 

из m-мерного численного пространства S: А = (а{ а2

ап). Мно­

жество всевозможных линейных комбинаций линейно независимых столбцов матрицы А образует подпространство 5, пространства S.

Преобразование (2.7.2) соотносит каждому вектору х е R век­ тор у подпространства 5, пространства S. Действительно,

П

у —,Ах = ^ x iai G Sv

/=1

Сдругой стороны, каждый вектор у подпространства 5,, являясь

линейной комбинацией столбцов матрицы А, представляется произ­ ведением Ах, где А — столбцовая п X 1 матрица, х — вектор п- мерного пространства R.

Итак, совокупность векторов Ах, где х — любой вектор из R , является подпространством 51m-мерного пространства S.

Выясним, какова размерность подпространства Sr Покажем

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

ранг матрицы А = (а{

а2

ап) равен г и не равный нулю минор

порядка г

находится

на пересечении столбцов а*, а , , ..., а,

(1 < ij < i2 <

< ir < п)

*i

г

и некоторых строк матрицы.

 

Допустим, что эти столбцы линейно зависимы, т.е. имеются

числа а 1, а2, ..., аг €: Ж,

не все равные нулю и такие, что

 

 

а 1а1 + а 2at - f ... + arai = 0 .

(2.7.3)

Равенство (2.7.3) эквивалентно следующей системе алгебраических уравнений:

аи ъ + аи а2 + ... + аи аг = 0,

1

2

г

а2/,а 1 + a2ia2+ ... + a2iar = 0,

Эта система однородных уравнений относительно а 1? а2, а г имеет только нулевое решение, так как ранг матрицы коэффициентов равен

числу неизвестных. Но это значит, что столбцы

Ч

, а} ,

а.г линей-

 

*2

Ч

но независимы.

Число линейно независимых столбцов матрицы А, с другой сто­ роны, не может быть больше, чем г. В самом деле, предположим, что имеются / (/ > г) линейно независимых столбцов

а-,, а},..., а, . Но тогда равенство

Ч*2 Ч

а.аf + а2а.- + ... + а,я* = 0

1

в2

1 Ч

может выполняться только тогда, когда = се2= ... = 0^ —0. Одна­ ко это не так, ибо эквивалентная система алгебраических уравнений

л1/«1 + аи«2 + -

+ «u«/ = °>

а2г«1 + % а2 + -

+ в а в, = 0,

ami,a1+ ami3a2 + -

+

= °>

как система однородных уравнений относительно ор а2, ..., а;, в

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

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

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

Пусть ранг матрицы-оператора А в преобразовании (2.7.2) ра­

вен г и линейно независимыми являются столбцы а; , а-,,..., а,

*1 *2

V

этой матрицы. Каждыйстолбец матрицыА есть линейная комби­

 

нация г ее линейно независимых столбцов а-., а; ,..., af. Значит, и

*1 *2

V

каждый вектор подпространства 5, есть линейная комбинация этих г столбцов, т.е. Sj есть подпространство, порожденное г линейно независимыми векторами а( , а, , ..., а-., и потому размерность мат-

рицы равна г, т.е. равна рангу матрицы А.

удовлетворяю­

Рассмотрим совокупностьвсех векторов

щих уравнению

 

Ах = 0.

(2.7.4)

Эти векторы образуют в R некоторое подпространство НА.

Размерность этого подпространства равна п — г.

В самом деле,

так как ранг матрицы А равен г, то система алгебраических уравнений

a U x l +

«12*2 +

••• +

а 1пХп -

°>

«21*1 +

а 22х 2 +

+

а 2пХ п =

° ’

a m lX l +

а т2х 2 +

+

а тпХ п = °>

эквивалентная соотношению (2.7.4), имеет ровно п — г линейно независимых решений.

Число измерений d пространства RA, состоящего из всех векто­ ров х Е R, удовлетворяющих условию (2.7.4), называется дефек­ том матрицы-оператора А.

На основании вышеизложенного

 

d — n — r.

(2.7.5)

§ 2.8. Неравенства Сильвестра

 

Пусть даны численные пространства: m-мерное

R, я-мерное

5, ^-мерное Т и линейные операторы — прямоугольные матрицы: А с размерами q х я и В с размерами пХ т . Пусть В отображает R в S, а оператор А отображает 5 в Г, так что

у = Вх (у Е 5, х Е R), z — Ау (г Е у Е 5).

Тогда оператор С =

Л2? — матрица с размерами q x m

— отобра­

жает R в Т:

 

 

z = A B x = C x (х Е R, z Е Г).

(2.8.1)

Обозначим через гЛ

гс ранги операторов (матриц)

А, В и С.

Множество всех векторов Ау(у Е S) образует подпространство AS, размерность которого равна рангу матрицы А, т.е. гА. Множество всех векторов Вх(х Е S) образует подпространство BR, размерность которого равна рангу матрицы В, т.е. гв . Наконец, множество всех векторов Ау, где у = Вх, а л: Е R, образует подпространство A(BR), размерность которого равна рангу матрицы А В = С, т.е. гс.

Так как BR С S,то A(BR) С AS, т.е. ABR — подпространство размерности гс — есть часть подпространства AS, имеющего раз­

мерность гА. Значит,

 

гс < гА.

(2.8.2)

Число линейно независимых решений уравнений

 

Ay = 0 ( у Е 5 )

(2.8.3)

равно дефекту d матрицы А. Имеем (см. (2.7.5))

 

d = л — гА.

(2.8.4)

Через d{ обозначим число линейно независимых решений уравнения (2.8.3), принадлежащих подпространству BR. Так как BR С S, то

dx d. (2.8.5)

Совокупность векторов z, удовлетворяющих равенству (2.8.1), можно представить так:

z = Ау (у Е BR).

(2.8.6)

Число линейно независимых векторов z, определенных равенст­ вом (2.8.1) или (2.8.6), равно, как указывалось выше, гс.

В подпространстве BR, размерность которого равна гв , имеются dx линейно независимых решений уравнения (2.8.3). Значит, в этом подпространстве имеются гв <1Хлинейно независимых век­

торов, которые уже не являются решениями уравнения (2.8.3). Эти векторы образуют некоторое подпространство S{ размерности

гв d {. Учитывая это, всю совокупность векторов z, удовлетворя­

ющих равенству

(2.8.1), можно представить и так:

 

 

 

Z = Ау

(у G

5, С BR).

 

 

Покажем, что размерность / подпространства S, равна гс, т.е.

 

 

/ = r B

 

dx= гQ.

 

В

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

Sl выберем

I

линейно независимых

векторов

y,

, у2, .... у,. Этим векторам соответствуют векторы

 

 

 

2, = Ayt

 

( /= 1 ,2 , ... ,/)

(2.8.7)

пространства Т. Векторы (2.8.7) линейно независимы. В самом де­ ле, допуская их линейную зависимость, будем иметь

Ct,Z,

+

(*2z 2 "Ь ••• “Ь а 12 1 = 0 .

 

Это ведет к равенству

 

 

 

Л(а,у, +

а2у2 + ... + а, у,) = 0 .

(2.8.8)

Вектор а ,^ + а2у2 4-

+ а1у1 не равен нулю (равенство нулю оз­

начало бы, что векторы у,, у2, ..., у{ линейно зависимы), но тогда

равенство (2.8.8) противоречит тому условию, что векторы подпро­ странства £, не являются решениями уравнения (2.8.3). Значит,

Z,, z2, ..., Z[ линейно

независимы. Отсюда можно

сделать вывод,

что I = гв dx $ гс.

 

 

 

 

 

 

Покажем,что гс не может быть больше, чем /. Допустим, что

rc > I. Пусть линейно

независимыми

являются

векторы

z, , z2, ..., zr . Имеем

 

 

 

 

 

2,

i

(з*/ ^

/ 1» 2,...,

г^.)

 

 

и

+ arczrc = ^(a,y, + *гУг +

 

 

a,z, + a2z2 +

+

 

По предположению

векторы

у,, у2, ..., уг

линейно

зависимы

(гс > /). Поэтому имеются такие числа а,, а2, .... аг^ , не все рав-

ные нулю, что

<*!?! +«2^2 +

+ « ГсУГс = 0.

Но тогда

alzl 4* a2z2 + ... + arczrc =

что означает линейную зависимость векторов zp z2, ..., zr . Остает-

*г>

С

ся одна возможность, а именно:

гс ” Гв dV

(2.8.9)

Отсюда, в частности, следует

 

гс <

(2.8.10)

Из (2.8.4), (2.8.5) и (2.8.9) получаем

 

rc > rB - d = rA + r B - n'

(2.8.11)

Объединяя неравенства (2.8.2), (2.8.10) и (2.8.11), получаем нера­ венства Сильвестра

( 2.8. 12)

определяющие соотношение между рангами матриц А и В с разме­ рами <7 х и и л х т соответственно и рангом их произведения А В.

С л е д с т в и е . При умножении матрицы ранга г в любом порядке на неособенную матрицу ранг произведения остается равным г.

Пусть А — неособенная матрица порядка л, а В — матрица размеров п х т и ранга г (г ^ т> п). Применяя неравенства Силь­ вестра к произведению С = АВ, получим

Г + П — П < Гс «а Г , Л .

Отсюда гс — г.

Точно так же, если А — неособенная матрица порядка т , то для произведения С = ВА будем иметь

т + г — т ^ гс ^ г, т.

Отсюда снова гс = г.

§ 2.9. Разложение матрицы на прямоугольные множители

Пусть дана прямоугольная матрица с размерами т х л, т.е. мат­ рица вида С = (с, с2 сп), где ct — столбцовые матрицы-векто­

ры /л-мерного численного пространства R. Пусть ранг матрицы С равен г.

Столбцы матрицы С порождают г-мерное подпространство R{про­ странства R. Выберем произвольную систему г линейно независимых

векторов а1Уаъ ..., аг е и составим матрицу А — (ах а2 яг), имеющую размеры тХг .

Для любой, таким образом построенной матрицы А, существует такая прямоугольная матрица В размеров г х п и ранга г, что

С = АВ.

(2.9.1)

Разложение(2.9.1) принято называть скелетнымразложением прямоугольной матрицы. Для существования соотношения (2.9.1) необходимо, чтобы столбцы матрицы В = (At Ъг Ьп) удовлетво­ ряли уравнениям

 

 

Abi — ci

( /= 1, 2

, п).

(2.9.2)

Ранг матрицы коэффициентов уравнений А равен рангу расши­

ренной матрицы (а, аг

аг с*), ибо А составлена из г линейно

независимых

векторов

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

Ry и Cj Е Rv

Поэтому

уравнения (2.9.2) разрешимы относительно bt.

 

Остается показать, что ранг матрицы В также равен г. Учиты­

вая, что

в

данном случае

гс = гЛ= г,

неравенства Сильвестра

(2 .8 .1 2 )

можно записать в виде

 

 

г + гв - г §**г^г, гв.

Отсюда гв — г.

§2.10. Задачи и упражнения

1.Найти координаты многочлена

/(дс)= а0+ а 1д:+ а2х2 + ... + а лхя

а)

в базисе (1, х, х 1, ..., хп); б)

в базисе

(1, (х — а), (х — а)2, ....

Ос— a)n).

 

2. Найти матрицу перехода от

базиса (1, х, х 2,...,

х")

к базису

(1, (х — а), (х —а)2, .... Ос—а)”)

пространства многочленов степени: а)

меньше п;

б)

равной п.

 

 

 

 

3.Как изменится матрица перехода от одного базиса к другому, если: а) поме­ нять местами два вектора первого базиса? б) поменять местами два вектора второго базиса? в) записать векторы обоих базисов в обратном порядке?

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

5.Доказать, что все n-мерные векторы, у которых первая и последняя коорди­

наты одинаковы, образуют линейное подпространств^; найти размерность и базис этого подпространства.

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

7. Доказать, что следующие множества являются линейными пространствами по указанным операциям сложения и умножения.

1) Поле: вещественные числа; множество: вещественные числа. С л о ж е н и е : сложение вещественных чисел. У м н о ж е н и е на число: умножение веществен­ ного числа на вещественное число.

2) Поле: рациональные числа; множество: вещественные числа. С л о ж е н и е : сложение вещественных чисел. У м н о ж е н и е на число: умножение веществен­ ного числа на рациональное число.

3) Поле: любое число; множество: один вектор. С л о ж е н и е : определяется ра­ венством а + а —а. У м н о ж е н и е на число: умножение вектора а на любое чис­ ло а определяется правилом а а = а .

4) Поле: вещественные числа; множество: многочлены с действительными ко­ эффициентами от одной переменной, в том числе константы. С л о ж е н и е : сложе­ ние многочленов. У м н о ж е н и е на число: умножение многочлена на вещест­

енное число.

 

8. Доказать, что все многочлены степени

от одного неизвестного с вещест­

венными коэффициентами (или с коэффициентами от любого поля Р) образуют век­ торное пространство, если за операции взять обычные правила сложения многочленов и умножения многочлена на число. Определить базис и размерность этого про­ странства.

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

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

подпространства. (Матрица А называется кососимметрической, если АТ— —А.)

11. Найти какой-нибудь базис и размерность линейного подпространства L про­ странства Rп, если L задано уравнением х х + х 2+ • • • + х п=0 .

12.Доказать, что любая система попарно ортогональных ненулевых векторов линейно независима.

13.Доказать, что если размерность суммы двух линейных подпространств про­

странства Rn на единицу больше размерности их пересечения, то сумма совпадает с одним из этих подпространств, а пересечение — с другим.

14. Бели произведение матриц АВ определено, то доказать

rank {АВ) <min {rank А , rank В}.

15. Пусть А — неособенная матрица, произведения АВ и СА определены. До­ казать, что

rank (АВ) = rank В; rank (СД) =rank С.

16. Пусть сумма матриц А + В определена. Доказать, что

rank (А+ В) * rank А + rank В.

17. Показать, что дифференцирование является линейным преобразованием пространства всех многочленов степени <п от одного неизвестного с вещественными коэффициентами. 'Найти матрицу этого преобразования в базисе

а)

1, х, х 2, ..., хп;

 

 

б)

1, (х — о). (х-а)2

С*-а)

 

 

2!

 

 

о — вещественное число.

 

 

18.

Как изменится

матрица

линейного преобразования, если в базисе

е,, е

е„ поменять местами е, и е .?

1*&2

п

I

)

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

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