Элементы матричной алгебры

Матрица. Действия с матрицами

Матрицей размеров mn будем называть совокупность mn чисел расположенных в виде таблицы состоящей из m строк и n столбцов и записывать в виде:

.

Элементы матрицы – это числа a_ij ( ) составляющие её, где i – номер строки, j – номер столбца на пересечении которых находится элемент матрицы.

Основные операции над матрицами

Сложение и вычитание матриц

Определяется для матриц одинакового размера. Суммой (разностью) матриц A и B, обозначаемой A+B (A-B), называется матрица C, элементы которой определяются по формуле: c_ij=a_ij+b_ij, где a_ij и b_ij – соответственно элементы матриц A и B.

Умножение матрицы на число

Произведением матрицы A и числа , обозначаемым A, называется матрица B той же размерности, элементы которой b_ij=la_ij, где a_ij элементы матрицы A, т.е. при умножении матрицы на число надо все элементы матрицы умножить на это число.

Свойства

Пусть A, B, C – матрицы одного размера, a,  любые действительные числа, тогда:

1. A+ B= B+ A

2. (A+ B)+ C= A+(B+ C)

3. a(A+B)= aA+aB

4. (a+)A=aA+dA

5. (ad)A=a(dA)

Матрица, все элементы которой равны нулю называется нулевой. Пусть O – нулевая матрица, тогда:

6. A+ O= A (-1)A – противоположная к A и обозначается – A.

7. A+(- A)= O.

Транспонирование матриц

Матрица A^T, полученная из данной матрицы A заменой её строк столбцами с теми же номерами называется транспонированной:

Умножение матриц

Произведением матриц A_{m n} и B_{n p} называется матрица C_{m p}= A B (или проще AB), элементы которой , где , - элементы матриц A и B. Произведение AB существует только в том случае, когда первый множитель A имеет число столбцов, равное числу строк второго множителя B.

Свойства умножения

1. AB BA даже если оба произведения определены, но существуют матрицы A,B, такие что AB= BA, тогда они называются перестановочными.

Матрица E вида:

называется единичной матрицей. E – перестановочная с любой квадратной матрицей того же размера, т.е. AE=EA=A.

2. Умножение матриц ассоциативно, т.е. если определены произведения AB и (AB)C, то определены BC и A(BC) и выполняется равенство:

(AB)C=A(BC).

3. Умножение матриц дистрибутивно по отношению к сложению, т.е.:

A(B+C)=AB+AC

(A+B)C=AC+BC.

4. Для любого числа :

(AB)=(lA)B=A(lB).

5. 5. Если существует AB, то определено (существует) B^ТA^Т и выполняется равенство:

(AB)^Т= B^ТA^Т.

Обратная матрица

Матрица X, удовлетворяющая вместе с заданной матрицей A равенствам XA=AX=E,

называется обратной к A и обозначается A^-1.

Определители

Определители (детерминанты) рассматриваются только для квадратных матриц.

Определитель n-го порядка это число, записываемое в виде таблицы

и может быть вычислено по элементам этой таблицы в соответствии с указанными ниже правилами.

Минором M_ij элемента a_ij называется определитель (n-1)-го порядка, полученный из определителя n-го порядка вычеркиванием i-й строки и j-го столбца.

Алгебраическое дополнение A_ij элемента a_ij определяется равенством:

A_ij=(-1)^i+j M_ij.

Определитель _n (det) находится по правилу:

,

а миноры M_1k – являются определителем (n-1)-го порядка, полученным из _n вычеркиванием 1-й строки и k-го столбца. Эта формула называется разложением по строке. Можно раскладывать по столбцу:

.

Определители первого, второго и третьего порядков

Определитель первого порядка .

Определитель второго порядка:

₂= .

Определитель третьего порядка, вычисленный разложением по первой строке:

₃= .

Правило вычисления определителя 3-го порядка равносильно правилу треугольников (правилу Саррюса):

₃=a₁₁a₂₂a₃₃+ a₁₂a₂₃a₃₁+ a₂₁a₃₂a₁₃-(a₁₃a₂₂a₃₁+ a₁₂a₂₁a₃₃+ a₂₃a₃₂a₁₁).

Основные свойства определителей

1. Определитель не изменится, если его разложить по любому столбцу (строке).

2. Значение определителя не изменится после замены всех его строк соответствующими столбцами и наоборот: detA=detA^T.

3. Если поменять местами две каких-нибудь строки (столбца), то определитель изменит свой знак.

4. Определитель с двумя одинаковыми строками (столбцами) равен нулю.

5. Если все элементы некоторой строки (столбца) имеют общий множитель, то его можно вынести за знак определителя. Следовательно, если элементы какой-либо строки (столбца) умножить на число l, то определитель умножится на это число l.

6. Если все элементы какой-либо строки (столбца) равны нулю, то определитель тоже равен нулю.

7. Определитель у которого элементы двух строк (столбцов) соответственно пропорциональны равен нулю.

8. Если каждый элемент какой-либо строки (столбца) определителя представляет собой сумму 2-х слагаемых, то такой определитель равен сумме 2-х определителей. В первом из которых соответствующая строка (столбец) состоит из 1-х слагаемых, а во втором – из вторых.

9. Определитель не изменится, если ко всем элементам какой-либо строки (столбца) прибавить соответствующие элементы другой строки (столбца), умноженные на одно и тоже произвольное число l.

Обратная матрица (продолжение)

Поскольку AA^-1=A^-1A, то они перестановочные и обе должны быть квадратными одинакового порядка n.

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

Утверждение. Каждая невырожденная квадратная матрица имеет обратную, и притом только одну.

.

Матрица A^* называется присоединенной, её элементы являются алгебраические дополнения A_ij транспонированной матрицы A^T.

Пример 1. Дана матрица A. Найти обратную ей A^-1. .

Вычислим определитель det A.

det A=2(-5)1+(-4)31+1(-1)1-1(-5)1-2(-1)3-1(-4)1=-10-12-1+5+6+4=-13+5=-8.

det A=-8  0. Вычислим алгебраические дополнения.

A₁₁ = -2, A₁₂ = 2, A₁₃ = 4, A₂₁ = 3, A₂₂ = 1, A₂₃ = -2, A₃₁ = -7, A₃₂ = -5, A₃₃ = -6.

Запишем результат: .

Ранг матрицы

Рангом матрицы A, обозначается rangA (RgA), называется наибольший порядок порожденных ею определителей, отличных от нуля.

Утверждение. Ранг матрицы не изменится, если:

1. поменять местами любые две строки (столбца);

2. умножить каждый элемент строки (столбца) на одно и тоже число l0;

3. прибавить к элементам строки (столбца) соответствующие элементы любой другой строки (столбца), умноженные на один и тот же множитель.

Преобразования 1 – 3 называются элементарными.

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

Базисным минором матрицы называется всякий отличный от нуля минор, порядок которого равен рангу данной матрицы.

Основные методы нахождения ранга матрицы

1. Метод нулей и единиц. C помощью элементарных преобразований любую матрицу можно привести к виду, когда каждый её ряд будет состоять только из нулей или из нулей и одной единицы. Тогда число оставшихся единиц и определит ранг исходной матрицы.

Пример 1. Найти ранг матрицы

.

Умножим 3-й столбец на1/2, затем 1-ю строку умножим на 2 и вычтем её из четвертой строки. Теперь 3-й столбец содержит три нуля и единицу в первой строке. Легко делаем нули в первой строке на первой, второй, четвертой и пятой позициях. Имеем

A .

Теперь 4-ю строку складываем со 2-й и 3-ей, после чего делаем нули в 4-й строке всюду, кроме единицы на пересечении четвертой строки и второго столбца. В результате этих преобразований имеем:

A   .

Получили три единицы. Следовательно, rangA=3.

2. Метод окаймляющих миноров. Минор M_k+1 порядка k+1, содержащий в себе минор M_k порядка k, называется окаймляющим минор M_k.. Если у матрицы A существует минор M_k 0, а все окаймляющие его миноры M_k+1 = 0, то rangA=k.

Пример 2. Найти ранг матрицы

.

Имеем . Для M₂ окаймляющими будут только

,

каждый из которых равен нулю. Поэтому rangA=2.

Вопросы для самоконтроля

1. Что такое называют mn матрицей?

2. Как определяется произведение матрицы на число? Сумма матриц? Всякие ли две матрицы можно сложить?

3. Как определяется произведение двух матриц?

4. Обладает ли операция умножения матриц коммутативностью?

5. Что такое единичная матрица?

6. Что называют определителем квадратной матрицы?

7. Что такое обратная матрица?

8. Что такое ранг матрицы?

Решение систем линейных уравнений

Системой линейных уравнений с тремя неизвестными называется система вида:

Решить систему (1) – значит найти такие значения (x, y, z), при подстановке которых в систему уравнения становятся верными тождествами. Для решения таких систем существует много способов. Рассмотрим три основных: метод Гаусса, метод Крамера, метод обратной матрицы.

Метод Гаусса

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

• умножение строки на любое число;

• прибавление к одной строке другой строки или суммы нескольких строк;

• умножение строки на любое число и прибавление к другой строке, также умноженной на любое число.

После приведения системы к "треугольному виду" появляется возможность найти значение переменной z из последнего уравнения, затем оно подставляется во второе уравнение и находится значение y. Полученные значения z и y подставляются в первое уравнение, откуда находится x.

Вычисления лучше производить с коэффициентами системы, записанными в форме таблицы, которая называется расширенной матрицей системы:

В результате преобразований должна получиться матрица вида:

Возвращаясь к интерпретации каждой строки матрицы как коэффициентов уравнения, можно записать:

C₃z=D₃, откуда z=C₃/D₃;

B₂y+C₂z=D₂, откуда y=(D₂-C₂z)/B₂;

a₁x+b₁y+c₁z=d₁, откуда x=(d₁-b₁y-c₁z)/a₁.

Таким образом, получили решение системы (1): (x,y,z).

Пример 1. Решить систему уравнений методом Гаусса:

Составим матрицу системы:

Умножим первую строку на (-3) и прибавим ко второй строке, записывая сумму на месте второй строки:

Умножим первую строку на (-9) и прибавим к третьей строке, записывая сумму на месте третьей строки:

Умножим вторую строку на (-4), а третью строку на 7 и сложим их, записывая сумму на месте третьей строки:

Из третьей строки имеем:

11z=22, откуда z=2.

Из второй строки имеем:

7y-15z=-16, откуда y=(-16+15z)/7, подставим z=2: y=(-16+30)/7=14/7=2.

Из первой строки находим x:

x-2y+3z=3, откуда x=3+2y-3z, подставляем z=2, y=2: x=3+4-6=1.

Ответ x=1, y=2, z=2.

Метод Крамера

Этот метод позволяет находить единственное решение системы с помощью определителей.

Определителем второго порядка называется число, вычисляемое по формуле:

Определителем третьего порядка называется число, вычисляемое по формуле:

Эта формула называется "разложение по первой строке".

Пример 2.

Пусть требуется решить систему уравнений:

1. Вычисляем определитель системы: если Δ ≠ 0, тогда система имеет единственное решение, в противном случае решений либо бесконечно много, либо не существует.

2. Вычисляем вспомогательные определители:

3. Находим x, y, z:

.

Матричный метод

Пусть для системы (1) основная матрица A вида - невырожденная, т.е. detA0. Тогда для A существует единственная обратная матрица A^-1. Введем в рассмотрение матрицы столбцы для неизвестных и свободных членов:

.

Тогда систему (1) можно переписать в матричной форме:

AX=D.

Умножив это матричное уравнение слева на A^-1, получим A^-1 AX= A^-1D, откуда EX=X= A^-1D. Следовательно, матрица-решение X находится как произведение A^-1 и D.

Пример 3. Решить систему уравнений матричным методом:

Имеем:

Обратная матрица

(см. пример). Находим:

,

т.е. x=2, y=0, z=-1 – решение данной системы.

Теорема (Кронекера – Капелли). Для того чтобы система m линейных алгебраических уравнений относительно n неизвестных x₁, x₂, … x_n

была совместна (имела решение), необходимо и достаточно, чтобы ранг основной матрицы

системы и ранг так называемой расширенной матрицы

системы были равны, т.е. rangA= rangB=r. Далее, если rangA= rangB и r=n, то система имеет единственное решение; если r<n, то система имеет бесконечное множество решений, зависящих от n – r произвольных параметров.

Вопросы для самоконтроля

1. Что называют системой линейных уравнений?

2. Что значит решить систему уравнений?

3. Что называют определителем системы уравнений?

4. В чем состоит метод Гаусса?

5. Каким образом при помощи определителей можно решить систему уравнений?

6. Что можно сказать, если определитель системы равен нулю?

7. Как записать в матричном виде систему уравнений?

8. О чем говорится в теореме Кронекера – Капелли?