книги / Матричное исчисление с приложениями в теории динамических систем
..pdfх € |
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г |
г |
|
и некоторых строк матрицы. |
|
|||
Допустим, что эти столбцы линейно зависимы, т.е. имеются |
||||
числа а 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. Доказать, что матрицы одного и того же линейного преобразования в двух базисах тогда и только тогда совпадают, когда матрица перехода от одного из этих базисов к другому перестановочна с матрицей этого линейного преобразования в од ном из этих базисов.