Чтобы изучить линейную алгебру, вы можете прочесть и вникнуть в книгу И. В. Белоусова "Матрицы и определители". Однако она написана строгим и сухим математическим языком, который людям со средним умом воспринимать тяжело. Поэтому я сделал пересказ наиболее трудных для понимания мест этой книги, стараясь изложить материал как можно понятнее, максимально используя для этого рисунки. Доказательства теорем я опустил. Признаться, я и сам не стал в них вникать. Верю г-ну Белоусову! Судя по его работе, он грамотный и толковый математик. Скачать его книгу можно по адресу http://eqworld.ipmnet.ru/ru/library/books/Belousov2006ru.pdf Если собираетесь вникать в мою работу, это нужно сделать, потому что я буду на Белоусова часто ссылаться.
Начнём с определений. Что такое матрица? Это прямоугольная таблица чисел, функций или алгебраических выражений. Зачем нужны матрицы? Они сильно облегчают сложные математические расчёты. У матрицы можно выделить строки и столбцы (рис. 1).
Строки и столбцы нумеруются, начиная слева сверху (рис. 1-1). Когда говорят: матрица размером m n (или m на n), подразумевают под m количество строк, а под n количество столбцов. Например, матрица на рисунке 1-1 имеет размер "4 на 3", а не "3 на 4".
Смотрите на рис. 1-3, какие бывают матрицы. Если матрица состоит из одной строки, она называется матрицей–строкой, а если из одного столбца, то матрицей–столбцом. Матрица называется квадратной n–го порядка, если число строк у неё равно числу столбцов и равно n. Если все элементы матрицы равны нулю, то это нулевая матрица. Квадратная матрица называется диагональной, если равны нулю все её элементы, кроме расположенных на главной диагонали.
Сразу объясняю, что такое главная диагональ. На ней номера строк и столбцов одинаковые. Идёт она слева направо сверху вниз. (рис. 3) Элементы называются диагональными, если они расположены на главной диагонали. Если все диагональные элементы равны единице (а остальные нулю), матрица называется единичной. Две матрицы A и B одинакового размера называются равными, если все их элементы одинаковые.
Умножение двух матриц определено лишь тогда (в переводе на русский: матрицы можно умножать лишь тогда),
когда число столбцов первой матрицы в произведении равно числу строк второй (рис. 7 , наверху, синие скобки). Чтобы лучше запомнить: цифра 1 больше похожа на столбец. В результате умножения получается матрица размером (смотри рисунок 6).
Чтобы было проще запомнить, что на что надо умножать, предлагаю следующий алгоритм: смотрим рисунок 7. Умножаем матрицу A на матрицу B. У матрицы A два столбца, у матрицы B две строки - умножать можно.
1) Займёмся первым столбиком матрицы B (он у неё один только и есть). Записываем этот столбик в строку (транспонируем столбик, о транспонировании чуть ниже).
2) Копируем эту строку, чтобы у нас получилась матрица размером с матрицу A.
3) Умножаем элементы этой матрицы на соответствующие элементы матрицы A.
4) Складываем получившиеся произведения в каждой строчке и получаем матрицу-произведение из двух строк и одного столбца.
На рисунке 7-1 даны примеры умножения матриц, которые размером поболее.
1) Здесь у первой матрицы три столбца, значит у второй должно быть три строчки. Алгоритм ровно тот же, что в предыдушем примере, только тут в каждой строчке три слагаемых, а не два.
2) Здесь у второй матрицы два столбца. Сначала проделываем алгоритм с первым столбцом, затем со вторым, и получаем матрицу "два на два".
3) Тут у второй матрицы столбец состоит из одного элемента, от транспонирования столбец не изменится. И складывать ничего не надо, так как в первой матрице всего один столбец. Проделываем алгоритм три раза и получаем матрицу "три на три".
Имеют место следующие свойства:
1. Если сумма
B + C и произведение AB существуют, то A (B + C ) = AB + AC
2. Если произведение AB существует,
то x (AB) = (xA) B = = A (xB).
3. Если произведения
AB и BC существуют, то A (BC) = (AB) C .
Если произведение матриц AB существует, то произведение BA может
не существовать.
Если даже произведения AB и BA существуют, то они могут оказаться
матрицами разных размеров.
Оба произведения AB и BA существуют и являются
матрицами одинакового размера лишь в случае квадратных матриц A и B
одного и того же порядка. Однако, даже в этом случае AB
может не равняться BA.
Транспонирование -это преобразование матрицы A в матрицу AT , при котором строки матрицы A записываются в столбцы AT с сохранением порядка. (рис. 8). Можно сказать по другому: столбцы матрицы A записываются в строки матрицы AT с сохранением порядка. Обратите внимание, как при транспонировании меняется размер
матрицы, то есть количество строк и столбцов. Также обратите внимание, что элементы на первой строке, первом столбце, и последней строке, последнем столбце остаются на месте.
На рисунке 9 вверху слева изображена симметричная матрица. Её элементы, симметричные относительно главной диагонали, равны.
А теперь определение: Квадратная матрица A называется симметричной, если AT=A .
То есть симметричная матрица при транспонировании не меняется.
В частности, симметричной является любая диагональная матрица.
(Такая матрица изображена на рис. 2).
Что такое определитель матрицы? Это такое число, которое вычисляется при помощи специальных операций с (квадратной!) матрицей. При расчётах в зависимости от значения определителя можно делать важные выводы, например, перпендикулярны или нет прямые.
Определителем |A| матрицы первого порядка A , или
определителем первого порядка, называется число , равное матричному элементу
Определителем |A| матрицы второго порядка , или
определителем второго порядка, называется число , определяемое формулой
(смотрите рис. 12) .
То есть, чтобы найти этот определитель, перемножаем элементы главной диагонали и вычитаем из них произведение
оставшихся. Первая цифра индекса элемента на рисунке (выделена голубым) означает номер его строки,
а вторая цифра, выделенная жёлтым, - номер столбца. Скобками выделены члены определителя.
Найти определитель матрицы третьего порядка (определитель третьего порядка)
гораздо сложней.
Далее изложен один из способов вычисления определителя третьего порядка. Другие способы изожены в разделах 5 и 6.
Сначала объясню некоторые понятия, которые понадобятся при вычислении.
Представьте себе последовательность чисел (именно вот с такими элементами, в таком порядке):
1 2 3 4 5 6(рис. 10).
a1 a2 a3 a4 a5 a6Чему равны эти числа из второй последовательности? Смотрим рис. 10. Каждое из чисел второй последовательности равно какому - то из чисел первой последовательности, причём одинаковых чисел во второй последовательности нет. Я выражусь может быть не очень правильно с точки зрения строгой математики, но более понятно: вторая последовательность - это та же первая последовательность, с теми же самыми элементами, но расположенными в ином порядке. Вот эта вторая последовательность называется перестановкой степени 6. Из чисел второй последовательности, то бишь перестановки можно образовать всевозможные пары, например: (a1=4 , a3=3); (a5=6 , a4=2); (a2=1 , a6=5); (a5=6 , a1=4). Теперь смотрите на рисунок 10. Если один элемент больше другого, а индекс его, наоборот, меньше, чем у другого, то такая пара называется инверсией.
Белоусов описал простой способ нахождения числа инверсий в перестановке. Смотрим рисунок 11.
Находим в перестановке число 1 и считаем, сколько чисел слева от него. Запоминаем или записываем.
Затем число 1 зачёркиваем. Зачёркнутые числа в дальнейших подсчётах не учитываются. Далее
находим число 2, считаем, сколько незачёркнутых чисел слева от него. Запоминаем или записываем.
Число 2 зачёркиваем. Затем находим число 3 и проделываем те же операции, что с числами 1 и 2.
И так проходим все числа перестановки. Затем те числа, которые мы запомнили или записали, суммируем.
Получаем число инверсий в перестановке.
А зачем нужно знать число инверсий в перестановке? У перестановки есть такая характеристика, как знак перестановки, или, не по русски, сигнатура перестановки (обозначается как sign). Так вот, если число инверсий в перестановке чётное, sign=1, и сама перестановка называется чётной. А если число инверсий в перестановке нечётное, sign=-1, и перестановка называется, соответственно, нечётной. Для тех, кому особо интересно, излагаю, как вычисляется этот знак перестановки. Он равен (-1)(число инверсий) Если число инверсий чётное, знак равен 1, если нечётное, то -1 .
Смотрим рис. 13, позиция 1 . На нём матрица и её определитель. Внимательно рассмотрим определитель. Рассмотрим какой-нибудь его член. Он состоит из трёх элементов. Посмотрите на первые цифры их индексов (голубые, означают номера строк). Теперь так же посмотрите на другие члены определителя. Вы видите закономерность ? В каждом члене определителя первые цифры образуют последовательность 1, 2, 3. То есть первый элемент члена определителя находится на первой строке, второй на второй, третий на третьей. Теперь посмотрите на жёлтые индексы, то бишь номера столбцов. А здесь какая закономерность, и чем она отличается от предыдушей? В каждом члене определителя одинаковых жёлтых номеров нет. Это значит, в каждом члене определителя элементы находятся на разных столбцах. В каждом члене определителя есть и 1, и 2, и 3, но они "перемешаны", то есть образуют перестановку. Вот каждому возможному варианту этой перестановки и соответствует член определителя. Я это попытался графически изобразить на рис. 13 . позиция 2. Смотрим снова на определитель. Некоторые его члены прибавляются к общей сумме, то есть имеют положительный знак, другие, наоборот, отнимаются, (знак их отрицательный). Как определить, прибавлять член или отнимать? А вот как - смотрим на жёлтые индексы. В каждом члене они образуют перестановку. То, что я вам толковал про инверсии и знак перестановки, ещё не вылетело у вас из головы? Если вылетело, снова изучите рис.11, перечитайте соответствующий абзац. Вычисляем знак перестановки. Если он положительный, член определителя прибавляется, если отрицательный - отнимается. Белоусов предлагает более удобный способ определить знак члена определителя. Смотрим рис. 13, позиция 3. Один член выделен синим, определяем его знак. Соединяем элементы отрезками (всеми возможными способами). Смотрим на рисунок. Два отрезка имеют так называемый положительный наклон, один - отрицательный. Что такое положительный наклон отрезка? Если в нём направление перехода на большую строку совпадает с направлением перехода на больший столбец, наклон отрезка положительный, если не совпадает - отрицательный. Отрезок с отрицательным наклоном соответствует инверсии в перестановке. Не поняли? Помозгуйте. Так вот, если число отрезков с отрицательным наклоном чётное, член определителя прибавляется, если нечётное, то вычитается.
Можете посмотреть, что такое определитель, определение его то бишь, у Белоусова на стр. 32 . После определения у Белоусова следуют примеры, внимательно их изучите, разжёвывать их мне неохота. Мы можем сделать некоторые важные выводы. Если какая–либо строка (или столбец) квадратной матрицы состоит из одних нулей, этот нуль входит множителем во все члены определителя, следовательно определитель равен нулю. После определения у Белоусова следуют примеры, внимательно их изучите, разжёвывать их мне неохота.
Посчитаем определитель какой - нибудь диагональной матрицы. Огого! Да у нас только один член определителя, который состоит только из диагональных элементов, не равен нулю! Делаем важный вывод: Определитель диагональной матрицы равен произведению её диагональных элементов. А следовательно определитель единичной матрицы равен 1.
Привожу без доказательств теорему, которая доказана в работе Белоусова, и следствия из неё.
Теорема: Определитель произведения двух (а также нескольких) квадратных матриц одного и того же порядка равен произведению их определителей.
Следствие: Определитель целой положительной степени квадратной матрицы равен определителю этой матрицы, возведённому в ту же степень.
Обратите внимание: если при умножении матрицы переставить местами, в результате получатся разные матрицы. Однако согласно теореме, определитель у них будет одинаковый.
Что такое минор? Возьмём какой нибудь элемент квадратной матрицы, например, элемент A22 на рисунке 15-1, позиция 1. Если у матрицы убрать строку, на которой расположен этот элемент, а также столбец, на котором расположен этот элемент, мы получим матрицу меньшего размера. Определитель этой матрицы и называется минором элемента (обозначается греческой буквой "мю"). Обратите внимание, что минор элемента вычислить гораздо легче, чем определитель матрицы. Если матрица второго порядка (рис. 15-1, позиция 2), то минор элемента и вовсе равен одному из других элементов.
Введём ещё понятие - алгебраическое дополнение элемента. Величина алгебраического дополнения зависит от суммы номеров столбца и строки, на которых расположен элемент. Если эта сумма чётная, алгебраическое дополнение равно минору элемента, если нечётная - то минору, взятому с отрицательным знаком. Обозначается алгебраическое дополнение греческой буквой "альфа".
Замечание 5.1 Его формулировку на математическом языке и доказательство смотрите у Белоусова. Я же изложу его "простым" языком. Если матрицу транспонировать, алгебраические дополнения её элементов (переместившихся на другие "места") останутся прежними. Смотрите рисунок 15-2.
Теорема 5.1 Определитель квадратной матрицы равен сумме произведений элементов любой строки на их алгебраические дополнения. Важнейшая теорема! На ней основан эффективный способ нахождения определителей.
Эта теорема подтверждает следствие 1 теоремы 4.3. Смотрим рис. 15-3. Две матрицы различаются только одной строкой, причём соответствующие элементы этой строки у матрицы C в два раза больше, чем у матрицы B. Если вычислить определители матриц через алгебраические дополнения этих строк, определитель матрицы C окажется в два раза больше матрицы B. Вывод: Общий множитель всех элементов строки матрицы можно вынести за знак определителя.
Теорема 5.2 Сумма произведений элементов какой-либо строки
определителя n–го порядка на алгебраические дополнения элементов
другой его строки равна нулю.
Теоремы 5.1 и 5.2 в равной степени применимы как строкам, так и к столбцам матрицы.
Смотрим рис. 16 позиция 1. На нём квадратная матрица четвёртого порядка. Нам надо эту матрицу привести к треугольному виду (рис. 16, позиция 2), потому что определитель треугольной матрицы равен простому произведению диагональных элементов, и его легко вычислить. Товарищ Белоусов описывал такие матрицы в примерах 4.2 и 4.3 и там же доказал это (утверждение). Я надеюсь, что всем
понятно, что такое матрица в треугольном виде. У неё все элементы, расположенные ниже диагональных, равны нулю.
Сначала вспомним два правила. Первое: если переставить местами две строки (или столбца) определителя, то он сохранит абсолютное значение, но поменяет знак на противоположный. Второе: если к элементам строки определителя прибавить соответствующие элементы другой его строки, умноженные на произвольное число, то определитель не изменится. Итак, начнём.
Смотрим на первый столбец первоначальной матрицы. Если бы все его элементы были равны нулю, то всю последущую работу нам не нужно было бы проводить, потому что определитель матрицы с таким столбиком (с одними нулями) равен нулю. Мы это с вами уже проходили. В нашей же матрице нам надо добиться, чтобы первый (верхний) элемент столбика не был равен нулю, а все остальные были равны нулю. Сначала делаем перестановку строк (результат на рис. 16, позиция 3). Я сделал перестановку строк два раза, а не один, чтобы знак определителя не изменился. Далее нам надо добиться, чтобы вместо пятерки на нижней
зелёной строке появился нуль. Используем второе правило. (Напоминаю, смотрим позицию 3, первый столбец). Вопрос: на сколько нужно умножить 2, чтобы получившееся произведение прибавить к 5, и в итоге получился нуль? Ответ: -5/2 .
Умножаем, прибавляем, и первый столбик приобретает нужный нам вид. Но нам придётся все остальные элементы первой строки тоже умножить на -5/2 и прибавить к соответствующим элементам последней строки. Делаем это и получаем матрицу на рис 16, поз. 4 . Элементы первой строки в дальнейших преобразованиях уже не участвуют. Первый же столбик от дальнейших преобразований не изменится. Вы в этом убедитесь. Далее преобразуем второй столбик. Если бы все его элементы, расположенные ниже первой строки, были равны нулю, то тогда и определитель был бы равен нулю. (Объяснение для крутых математиков: в этом случае два первых столбика оказались бы пропорциональными, и, согласно свойствам определителей, определитель был бы равен нулю.) Нам же придётся добиться, чтобы все элементы второго столбика, расположенные ниже главной диагонали, были равны нулю. Итак , позиция 4. Умножаем элементы второй (второй, не первой) строки на -4/1, то есть на -4, и прибавляем произведения к соответствующим элементам третьей строки. Результат на позиции 5. Далее умножаем элементы второй строки на 7/1, то есть на 7 и прибавляем к элементам четвёртой строки. Второй столбик принимает нужный нам вид (рис. 16, поз. 6). Вторая строка в дальнейших преобразованиях не участвует, а второй столбик от них не изменится. Нам осталось добиться, чтобы нижний элемент третьего столбика был равен нулю. Позиция 6. Умножаем элементы третьей строки на -16/-7= 16/7 и прибавляем к элементам четвёртой строки, и получаем треугольную матрицу в окончательном виде. (поз. 2) Теперь можно перемножить диагональные элементы и получить определитель. У меня получилось -29.
Сначала повторите теорему 5.1 и её следствие, а также следствия 1 и 3 теоремы 4.3 Будем вычислять определитель матрицы пятого порядка (рис. 16-0). В четвёртом столбике у неё имеется два нуля, что весьма кстати. Определитель можно вычислить через алгебраические дополнения элементов этого столбика. Чтобы уменьшить объём вычислений, сначала добъёмся, чтобы все элементы этого столбика, кроме одного, были равны нулю.
2) Умножили вторую строку на -1 и прибавили к первой.
3) Умножили пятую строку на -3 и прибавили ко второй.
4) Понижаем порядок определителя. Сумма номеров столбца (четвёртый) и строки (пятая) является нечётной, поэтому определитель пятого порядка равен -2 определителя четвёртого порядка.
5) Работаем со вторым столбиком определителя четвёртого порядка. Четвёртую строку прибавили ко второй.
6) Третью строку умножили на -7 и прибавили к четвёртой.
7) Третью строку умножили на 2 и прибавили к первой.
8) Понижаем порядок определителя. Сумма номеров третьей строки и второго столбика является нечётной, поэтому определитель пятого порядка равен 2 определителя третьего порядка.
9) Умножили первую строку на 2 и прибавили к третьей.
10) Вынесли общий множитель второй строки за знак определителя.
11) Выносим общий множитель третьей строки за знак определителя.
12) Теперь вспомним, что строки и столбцы у нас равноправны. Делаем "разложение определителя третьего порядка по элементам третьей строки". В третьей строке два элемента не равны нулю, и нам придётся вычислить два определителя второго порядка. В результате у меня получилось 1128.
Посмотрите ещё примеры у Белоусова на стр. 56 - 58.
Чтобы найти обратную матрицу, можно проделать следущее:
1) Найти определитель исходной матрицы. Если он равен нулю, матрица вырожденная, и обратной к ней матрицы не существует.
2)Транспонировать исходную матрицу.
3)Заменить в получившейся матрице все элементы их алгебраическими дополнениями.
4)Умножить получившуюся матрицу на число 1/A , где A - определитель исходной матрицы.
Однако есть более простые способы.
На рис. 16-1, позиция 1) изображена так называемая расширенная матрица. Она состоит из двух подматриц. Слева - исходная, к ней мы будем находить обратную матрицу. Справа - единичная матрица. С ней мы будем проделывать ровно те же операции, что и с исходной матрицей, и в результате исходная матрица превратится в единичную, а
единичная - в обратную к исходной.
Операции будем проводить следущие:
1) Умножение любой строки на число X, не равное нулю.
2) Прибавление любой строки, умноженной на число X, к другой строке. При этом остальные строки не меняются.
3) Перестановка двух строк между собой. При этом остальные строки не меняются
Почему именно эти операции, и почему в результате получится обратная матрица? Желающие докопаться до ответа на эти вопросы могут обратиться к работе Белоусова (глава 7).
Будем добиваться. чтобы в нули превратились все недиагональные элементы сначала в первом столбике, затем во втором, и так далее (впрочем, это не принципиально).
2) Умножили вторую строку на (-3) и прибавили к третьей строке.
3) Поменяли местами первую и вторую строки.
4)Умножили первую строку на (-2) и прибавили ко второй строке.
5) Умножили вторую строку на (-1) и прибавили к третьей строке.
6) Прибавили третью строку к первой строке.
7) Умножили третью строку на -1/2 .
8) Умножили третью строку на 3 и прибавили ко второй. Слева получаем единичную матрицу, справа - обратную к исходной
В работе Белоусова имеется несколько примеров нахождения обратных матриц таким способом (раздел 7.2). Рекомендую их посмотреть.
В единичной матрице все единички располагаются на разных строках и столбцах, и в то же время принадлежат главной диагонали (рис. 17, поз. 1). Строки и столбцы, образующие единичную матрицу, также назовём единичными.
Eсли бы единичные столбики находились в любом другом порядке, то, переставляя их, можно добиться, чтобы все единички расположились на главной диагонали.
Нам дана матрица А (рис. 17 позиция 2, левая часть , нам нужно найти обратную к ней матрицу. Будем преобразовывать матрицу A в единичную и одновременно будем аналогично преобразовывать единичную матрицу Е (средняя часть), чтобы она превратилась в матрицу, обратную к A. Будем пока рассматривать матрицы А и Е не как две разные матрицы, а как одну "расширенную" матрицу. Следует заметить, что в каждом столбце (как и в каждой строке) матрицы A хотя бы один
элемент должен быть отличен от нуля, потому что если это не так, то матрица является вырожденной (её определитель равен нулю, если забыли) и, следовательно, обратная к ней матрица не существует. И если у нас в левой части в результате какого-то шага преобразований получится строчка или столбец с одними нулями, дальнейшую работу можно не проводить.
Для контроля вычислений введём третий, самый левый столбец. В нём в каждой строчке находится сумма элементов соответствующих строчек матриц А и Е. Зачем она нужна, об этом чуть позже. Возьмём какой-либо элемент матрицы А, отличный от нуля (!) и назовём его разрешающим (рис. 17, поз. 2, первый столбик, средняя строка). Строку и столбец, на которых он расположен. так же
обзовём разрешающими. Обратите внимание на рисунке, что разрешающая строка проходит через обе матрицы, через А и через Е.
Будем превращать разрешающий столбец в единичный . Сначала превратим в единичку сам разрешающий элемент (если он не равен единице). Для этого разрешающую строчку разделим на ... правильно, разрешающий элемент (делим все элементы разрешающей строки, и те, что относятся к матрице А, и относящиеся к матрице В). Далее представим себе, что остальные элементы в разрешающем столбце мы преобразовали в нули известными нам способами. Тогда элементы матрицы А ,не лежащие на разрешающей строке и разрешающем столбце, вычисляются по формуле (рис. 17 поз. 5 первая формула) Элементы матрицы Е, не лежащие на разрешающей строке, вычисляются по формуле (вторая формула) Чтобы упростить вычисления, существует простое графическое правило
(рис. 17 поз. 6). Проведём
воображаемый прямоугольник через вычисляемый элемент, разрешающий элемент, и соответствующие элементы на разрешающей строке и разрешающем столбце. Вычисляемый элемент равен алгебраической сумме произведений элементов, лежащих на диагоналях этого прямоугольника, разделённой на разрешающий элемент, причём произведениё, содержащее разрешающий элемент, входит в сумму со знаком +, а другое произведение - со знаком - (минус). Однако гораздо удобнее использовать для вычислений матрицу с преобразованной разрешающей строчкой (рис. 17 поз. 3). Тогда сумму произведений элементов не надо делить на разрешающий элемент, поскольку он равен единице. После вычисления всех необходимых элементов в левой и средней части расширенной матрицы окончательно
преобразуем разрешающий столбец в единичный - заменяем все элементы в нём, кроме разрешающего, на нули (рис. 17 поз. 4 ). Мы сделали один полный шаг преобразований. Сейчас мы можем проверить правильность своих вычислений. Подсчитываем сумму элементов в каждой строчке расширенной матрицы, записываем её в правый столбец. Если мы всё правильно подсчитали, новую сумму в каждой строчке можно получить из первоначальной по формуле (рис. 17 поз. 5 третья формула). Далее преобразуем в единичный другой столбец матрицы А. Находим в нём элемент, отличный от нуля, не находящийся на разрешающей строчке , делаем его разрешающим и проводим вышеописанные преобразования - превращаем его в единичный.1
Товарищ Белоусов доказал, что предыдущий разрешающий столбец, который мы сделали единичным, при этом уже не изменится. И так делаем, пока все столбцы матрицы А не станут единичными. Затем, если нужно, переставляем столбцы и строки (аналогично и у матрицы А и у матрицы Е), чтобы единички расположились на главной диагонали матрицы А. Тогда в средней части получим искомую обратную матрицу. Проработайте ещё примеры у Белоусова на стр. 73 -75.
Белоусов в своей работе доказывает несколько свойств:
1) Если определитель матрицы равен x (x не равно нулю), определитель обратной к ней матрицы равен 1/x.
Это следует из того, что определитель единичной матрицы равен 1.
2) (A-1)-1=A Матрица, обратная к обратной - та же самая матрица.
3) (Am)-1=(A-1)m Если найти матрицу, обратную матрице A, а потом эту обратную матрицу (A-1) возвести в степень m, в итоге получится такая же матрица, которая получится, если проделать эти операции в другом порядке: сначала матрицу A возвести в степень m, а затем найти матрицу, обратную к Am.
4) (AB)-1=B-1A-1 Если две матрицы перемножить, а затем найти обратную к получившейся матрице, в итоге получится та же матрица, которая получится после умножения двух матриц, обратных к первоначальным.
5) (A-1)T=(AT)-1 Если найти матрицу, обратную матрице A, а потом эту обратную матрицу (A-1) транспонировать, в итоге получится такая же матрица, которая получится, если проделать эти операции в другом порядке: сначала транспонировать матрицу A, а затем найти обратную к транспониорванной.
Представим себе произвольную матрицу размерами m на n, причём n - наименьший размер, и натуральное число k,
k меньше либо равно n. Например, у матрицы на рисунке 18 m=4 , n=3 , k меньше либо равно 3 .
Вычеркивая из нашей матрицы произвольные строки и столбцы, можно получить квадратную
подматрицу размером k. А если из первоначальной матрицы вычеркнуть другие столбики и строки, можно получить
другую квадратную подматрицу, тоже размером k. (Для крутых математиков: сколько таких различных подматриц можно "выкроить" из матрицы? На странице 77
Белоусов приводит формулу вычисления количества таких подматриц. Восклицательный знак в формуле означает факториал.)
Определители этих подматриц
называют минорами k порядка данной матрицы.
Произвольный минор k порядка может быть равен или не равен нулю.
Белоусов привёл доказательство, что если все миноры k порядка матрицы равны нулю, то равны нулю и все её миноры более высокого порядка.
Если матрица А не нулевая, то всегда можно указать натуральное число r, обладающее следующими свойствами:
Вы помните, как выполняются арифметические действия со строками матрицы? Чтобы сложить две строки, складываем соотвествующие их элементы и получаем строку с элементами, равными сумме слагаемых элементов. Чтобы умножить строку на число, умножаем каждый элемент на это число и получаем строку с произведениями.
Смотрим матрицу из двух столбиков на рис. 19, позиция 1. Строку C можно получить, если строку A умножить на 2 и прибавить к строке B(19, поз. 1, 1). Говорят, что строка C является линейной комбинацией строк A и B. Частным случаем линейной комбинации является сумма строк и умножение строки на число (рис. 19, поз. 2). Вернёмся к строке C. Она получается из строк A и B. Следует ли из этого, что строку A можно получить из строк B и C? Следует (19, поз. 1-2). То есть строка А является линейной комбинацией строк B и C. Также и строка B является линейной комбинацией строк A и C (19, поз. 1-3). А вот ещё линейная комбинация из этих строк (19, поз. 1-4). Строка N состоит из одних нулей, то есть это нулевая строка. Строки A, B, C являются линейно зависимыми. Что это означает? Согласно определению линейной зависимости строк, найдутся такие числа (a=2, b=1,c=-1,), не равные одновременно нулю что линейная комбинация строк aA+bB+cC будет равна нулевой строке. Обратите внимание на требование, чтобы числа (a, b, c) не были равны одновременно нулю. Без него понятие линейной зависимости теряет смысл. Можно взять какие угодно строки. Если все коэффициенты при них будут равны нулю (19, поз. 1-5), их линейная комбинация, естественно, будет равна нулевой строке. Однако часть коэффициентов может быть равна нулю. Добавим к линейной комбинации (19, поз. 1-4) строку D, умноженную на 0 (19, поз. 1-6). Равенство осталось верным. Поэтому можно отнести к линейно зависимым строкам и строку D, какая бы она ни была. Мы можем сделать важный вывод: если часть строк линейно зависимы, то и все они линейно зависимы.
В книге Белоусова доказываются следущие утверждения:
Строки единичной матрицы (19, поз. 3) линейно независимы.
Замечание: Характер линейных зависимостей матрицы не меняется при произвольной перестановке её строк или столбцов.
Теорема 9.1: Строки матрицы зависимы тогда и только тогда, когда одна из них является линейной комбинацией остальных строк (19, поз. 1-4) и (19, поз. 1-1, 1-2, 1-3).
Теорема 9.2: Если часть строк линейно зависимы, то и все они линейно зависимы (19, поз. 1-1) и (19, поз. 1-6).
Следствие1: Если среди строк матрицы имеется нулевая строка, то эти строки линейно зависимы. В самом деле, если взять нулевую строку с каким угодно коэффициентом, не равным нулю, а остальные с нулевыми коэффициентами, получится нулевая строка (19, поз. 1-7).
Следствие 2 Если среди строк матрицы имеются пропорциональные, то все
они линейно зависимы (19, поз. 2).
Замечание 9.2 Понятие линейной зависимости (а стало быть, и приведённые выше теоремы и следствия)
применимо не только к строкам, но и к столбцам матрицы.
В последней главе Белоусов предлагает способ нахождения ранга матрицы. Матрица приводится к так называемому ступенчатому
виду при помощи элементарных преобразований. Вот эти преобразования:
