Линейные операции над матрицами
Суммой матриц А и В одинаковой размерности m n называется матрица С той же размерности, каждый элемент которой равен сумме элементов матриц А и В, стоящих на тех же местах:
Пример 1. Найти сумму матриц и .
Решение.
Вычислим элементы матрицы С = А + В, складывая элементы исходных матриц, стоящие на одинаковых местах:
Следовательно,
Произведением матрицы на число называется матрица той же размерности, что и исходная, все элементы которой равны элементам исходной матрицы, умноженным на данное число.
Пример 2. Найти матрицу 5А – 2В, если
.
Решение.
.
Итак, 5А – 2В .
Перемножение матриц
Произведением матрицы А размерности m p и матрицы В размерности называется матрица С размерности , каждый элемент которой определяется формулой: Таким образом, элемент представляет собой сумму произведений элементов i-й cтроки матрицы А на соответствующие элементы j-го столбца матрицы В.
Операция перемножения матриц некоммутативна, т.е. Действи-тельно, если существует произведение АВ, то ВА может вообще не существовать из-за несовпадения размерностей. Если существуют и АВ, и ВА, то они могут иметь разные размерности (если ).
Для квадратных матриц одного порядка произведения АВ и ВА существуют и имеют одинаковую размерность, но их соответствующие элементы в общем случае не равны.
Пример 3. Выяснить, можно ли умножить друг на друга матрицы
и .
Если произведение существует, вычислить его.
Решение.
Сравним размерности матриц А и В: A[3×2], B[2×2]. Следовательно, поэтому произведение АВ[3×2] существует, а произведение ВА – нет.
Найдем элементы АВ:
(ab)11 = 0 · 5 + 3 · 7 = 21; (ab)12 = 0 · 6 + 3 · 8 = 24; (ab)21 = 4 · 5 – 2 · 7 = 6;
(ab)22 = 4 · 6 – 2 · 8 = 8; (ab)31 = 1 · 5 – 1 · 7 = -2; (ab)32 = 1 · 6 – 1 · 8 = -2.
Таким образом, , ВА не существует.
Пример 4. Найти АВ и ВА, если
.
Решение.
Проверим возможность перемножения матриц, определив их размерность.
A[2×4], B[4×2]. Следовательно, n = l = 4, m = k = 2, поэтому матрицы АВ и ВА существуют, причем АВ[2×2], BA[4×4].
Для вычисления элементов матрицы С = АВ элементы строк матрицы А умножаются на соответствующие элементы столбцов матрицы В:
с11 = 2 · 2 + (-2)(-1) + 1 · 1 + 0 · 2 = 9
(сумма произведений элементов первой строки А на элементы первого столбца В; первый индекс вычисляемого элемента задает номер строки А, второй индекс – номер столбца В);
с12 = 2 · 2 + (-2) · 0 + 1 · 1 + 0 · 4 = 5;
с21 = -3 · 3 + 1 · (-1) + (-1) · 1 + 1 · 2 = -9;
с22 = -3 · 2 + 1 · 0 + (-1_ · 1 + 1 · 4 = -3.
Следовательно,
.
При вычислении элементов матрицы D = BA элементы строк В умножаются на элементы столбцов А:
d11 = 3 · 2 + 2 · (-3) = 0; d12 = 3 · (-2) + 2 · 1 = -4; d13 = 3 · 1 + 2 · (-1) = 1;
d14 = 3 · 0 + 2 · 1 = 2; d21 = -1 · 2 + 0 · (-3) = -2; d22 = -1 · (-2) + 0 · 1 = 2;
d23 = -1 · 1 + 0 · (-1) = -1; d24 = -1 · 0 + 0 · 1 = 0; d31 = 1 · 2 + 1 · (-3) = -1;
d32 = 1 · (-2) + 1 · 1 = -1; d33 = 1 · 1 + 1 · (-1) = 0; d34 = 1 · 0 + 1 · 1 = 1;
d41 = 2 · 2 + 4 · (-3) = -8; d42 = 2 · (-2) + 4 · 1 = 0; d43 = 2 · 1 + 4 · (-1) = -2;
d44 = 2 · 0 + 4 · 1 = 4.
Таким образом,
.
Определители
Определителем второго порядка называется число, полученное с помощью элементов квадратной матрицы 2-го порядка следующим образом:
.
При этом из произведения элементов, стоящих на так называемой главной диагонали матрицы (идущей из левого верхнего в правый нижний угол) вычитается произведение элементов, находящихся на второй, или побочной, диагонали.
Пример 5.
Определителем третьего порядка называется число, определяемое с помощью элементов квадратной матрицы 3-го порядка следующим образом:
Для того чтобы легче запомнить эту формулу, можно использовать так называемое правило треугольников. Оно заключается в следующем: элементы, произведения которых входят в определитель со знаком «+», располагаются так:
образуя два треугольника, симметричных относительно главной диагонали. Элементы, произведения которых входят в определитель со знаком «-», располагаются аналогичным образом относительно побочной диагонали:
Пример 6.
Вычислить определитель
.
Решение.
Вычислим определитель 3-го порядка, используя его определение:
Δ = 2·0·(-1) + (-3)·(-4)·2 + 5·1·1 - 2·0·5 -1·(-4)·2 – (-1)·1·(-3) =
= 0 + 24 + 5 – 0 + 8 – 3 = 34.
Пред тем, как перечислить основные свойства определителей, приведем определение понятия транспонирования матрицы.
Транспонированием матрицы называется операция, в результате которой меняются местами строки и столбцы с сохранением порядка их следования. В результате получается матрица А′, называемая транспонированнойпо отношению к матрице А, элементы которой связаны с элементами А соотношением a′ij = aji .
Основные свойства определителей
1. Определитель не изменяется при транспонировании, т.е.
2. При умножении элементов строки определителя на некоторое число весь определитель умножается на это число, т.е.
3. Определитель, имеющий нулевую строку, равен нулю:
4. Определитель, имеющий две равные строки, равен нулю:
5. Определитель, две строки которого пропорциональны, равен нулю:
6. При перестановке двух строк определителя он умножается на –1:
7.
8. Величина определителя не изменится, если к элементам одной строки прибавить соответствующие элементы другой строки, умноженные на одно и то же число:
Разложение определителя по строке
Минором элемента определителя называется определитель, полученный из данного путем вычеркивания строки и столбца, в которых стоит выбранный элемент.
Обозначение: выбранный элемент определителя, его минор.
Пример 7.
Для
Алгебраическим дополнением элемента определителя называется его минор, если сумма индексов данного элемента i+j есть число четное, или число, противоположное минору, если i+j нечетно, т.е.
При этом справедливо следующее утверждение: определитель равен сумме произведений элементов любой его строки или столбца на их алгебраические дополнения, т.е.
где i=1,2,3.
Таким образом, для вычисления определителя достаточно найти алгебраические дополнения к элементам какой-либо строки или столбца и вычислить сумму их произведений на соответствующие элементы определителя.
Пример 8.
Вычислим определитель из примера 6 с помощью разложения по строке. Для удобства вычисления выберем 2-ю строку, содержащую нулевой элемент (а22 = 0), поскольку при этом нет необходимости находить А22, так как произведение а22 А22 = 0. Итак,
(напомним, что определитель второго порядка, входящий в алгебраическое дополнение Aij, получается вычеркиванием из исходного определителя i-й строки и j-го столбца).
Тогда Δ = а21 А21 + а23 А23 = 1·2 + (-4)(-8) = 34.
Определители более высоких порядков
Определитель n-го порядка
есть сумма n! членов каждый из которых соответствует одному из n! упорядоченных множеств полученных r попар-ными перестановками элементов из множества 1,2,…,n.
Свойства определителей 3-го порядка справедливы и для определителей n-го порядка.
На практике определители высоких порядков вычисляют с помощью разложения по строке или столбцу. Это позволяет понизить порядок вычисляемых определителей и в конечном счете свести задачу к нахождению определителей 3-го порядка.
Пример 9.
Вычислить определитель 4-го порядка
.
Решение.
Преобразуем определитель так, чтобы три из четырех элементов какой-либо строки или столбца стали равными нулю. Для этого воспользуемся свойством 8. Его особенно удобно применять, если в определителе существует элемент, равный ±1. Выберем в качестве такого элемента а13 = 1 и с его помощью обратим все остальные элементы 3-го столбца в нуль. С этой целью:
а) к элементам 2-й строки прибавим соответствующие элементы 1-й строки;
б) из элементов 3-й строки вычтем элементы 1-й строки, умноженные на 2;
в) из элементов 4-й строки вычтем элементы 1-й строки
(напомним, что при этом величина определителя не изменится). Тогда
.
Разложим полученный определитель по 3-му столбцу:
.
Вычтем из элементов 1-й строки нового определителя удвоенные элементы 2-й строки:
и разложим этот определитель по 1-й строке:
Обратная матрица
Квадратная матрица А называется вырожденной, если , и невырожденной, если .
Квадратная матрица В называется обратной к квадратной матрице А того же порядка, если АВ = ВА = Е. При этом В обозначается .
Для существования обратной матрицы необходимо и достаточно, чтобы исходная матрица была невырожденной. Тогда
,
то есть ее элементами являются алгебраические дополнения к элементам транспонированной матрицы А, деленные на ее определитель.
Пример 10.
Найти обратную матрицу для матрицы
.
Решение.
Вычислим определитель матрицы А разложением по первому столбцу:
.
Следовательно, обратная матрица для матрицы А существует.
Найдем алгебраические дополнения а элементам матрицы А:
Значит,
.
Ранг матрицы
Минором порядка k матрицы А называется определитель, составленный из элементов, стоящих на пересечении любых k строк и k столбцов данной матрицы. Таким образом, каждый элемент матрицы является ее минором 1-го порядка.
Ранг матрицы – это порядок ее наибольшего ненулевого минора (обозначения: r(A), R(A), Rang A).
Пример 11.
Определить ранг матрицы
.
Решение.
Единственным минором максимального (3-го) порядка для матрицы А является ее определитель. Если ΔА ≠ 0, r(A) = 3; если ΔА = 0, r(A) < 3.
Найдем ΔА разложением по первой строке:
Следовательно, r(A) < 3. Поскольку матрица А содержит ненулевые элементы, r(A) > 0. Значит, r(A) = 1 или r(A) = 2. Если найдется минор 2-го порядка, не равный нулю, то r(A) = 2.
Вычислим минор из элементов, стоящих на пересечении двух первых строк и двух первых столбцов:
Для матриц большой размерности непосредственное вычисление всех миноров затруднительно. Поэтому в этом случае можно преобразовать матрицу к так называемому треугольному виду (когда элементы, стоящие ниже равны 0), воспользовавшись операциями, не изменяющими ранг матрицы (эквивалентными преобразованиями). К ним относятся:
1) транспонирование;
2) умножение строки на ненулевое число;
3) перестановка строк;
4) прибавление к элементам данной строки элементов любой другой строки, умноженных на ненулевое число;
5) вычеркивание нулевой строки.
Действительно, любая из этих операций переводит нулевые миноры в нулевые, а ненулевые – в ненулевые. Матрица, полученная в результате, не равна исходной, но имеет тот же ранг.
Пример 12.
Определить ранг матрицы
.
Решение.
У матрицы А существуют миноры до 4-го порядка включительно, поэтому
r(A) ≤ 4. Разумеется, непосредственное вычисление всех миноров 4-го, 3-го и т.д. порядка потребовало бы слишком много времени. Поэтому, используя элементарные преобразования, приведем матрицу А к треугольному виду. Поменяем местами 1-ю и 2-ю строки, чтобы элемент а11 стал равным 1:
А ~ .
Прибавим к третьей строке первую, ко второй – удвоенную первую, к четвертой – первую, умноженную на 3. Тогда все элементы 1-го столбца, кроме а11, окажутся равными нулю:
А ~ .
Вычтем вторую строку полученной матрицы из третьей и четвертой строк:
А ~
и вычеркнем нулевые строки:
А ~ .
Итак, ранг матрицы А равен рангу полученной матрицы размера 2×6, т.е.
r(A) ≤ 2. Минор
следовательно, r(A) = 2.