14-es_1 / Высшая математика (РТФ) / умк_Вакульчик_Элементы линейной алгебры

Последнюю формулу называют разложением определителя по первой строке. Легко проверить, что число, вычисленное по этой форму- ле, совпадает с числом, найденным по определению 1.3.2.

Из общих соображений ясно, что первая строка матрицы А не «луч- ше» других ее строк или столбцов. Поэтому такое же разложение можно записать для любого ряда матрицы и получить тот же определитель А.

Вычислим detА разложением по третьему столбцу, элементы которо-

Тогда detА = а13А13 + а23А23 + а33А33 = 3×1 + 1×1 + 2×1 = 6.

Применим теперь для вычисления того же определителя разложение по 2-й строке: а21 = 0, а22 = 1, а23 = 1, А21 = 2, А22 = 5, А23 = 1. Тогда

detA = 0×2 + 1×5 + 1×1 = 6.

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

Сформулируем основные свойства определителей, присущие опреде- лителям всех порядков. Некоторые из этих свойств поясним на определите-

лях 3-го порядка. (Доказательство свойств изучить самостоятельно).

1.При транспонировании значение определителя не изменяется.

2.При перестановке двух параллельных рядов определитель меняет знак на противоположный.

3.Если определитель содержит нулевой ряд, то его значение равно нулю.

a11 0 a13

a21 0 a23 = 0.

a31 0 a33

4.Если в определителе имеется два одинаковых ряда, то его значе- ние равно нулю.

5.Если какой-либо ряд определителя содержит общий множитель, то его можно вынести за знак определителя.

31

6. Если у определителя пропорциональны элементы каких-либо па- раллельных рядов, то его значение равно 0.

7. Если элементы какого-либо ряда определителя представляют со- бой суммы двух слагаемых, то определитель может быть разложен на сум- му двух соответствующих определителей:

8. Если к элементам какого-либо ряда определителя прибавить соот- ветствующие элементы параллельного ряда, умноженные на одно и то же число, то его значение не изменится:

9. Сумма произведений элементов какого-либо ряда определителя на алгебраические дополнения соответствующих элементов параллельного ряда равна нулю:

a11 a12 a13

a21 a22 a23 = а12А11 + а22А21 + а32А31 = 0.

a31 a32 a33

Рассмотрим несколько примеров.

Пример 10. Вычислим определитель верхнетреугольной матрицы

32

Решение. Будем вычислять его, раскладывая по1-му столбцу. Так как только один элемент этого столбца отличен от нуля a11 = 3 ¹ 0, а ос- тальные элементы равны нулю, то

Замечание 1.4.1. Легко видеть, что определитель верхнетреуголь- ной матрицы равен произведению диагональных элементов.

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

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

33

Пример 40. Вычислить определитель, пользуясь его свойствами

Решение. Заметим, что можно вынести за знак определителя об- щий множитель 2 из 1-го столбца и еще множитель 2 из 1-й строки, а из 2- го столбца выносится множитель 3

В последнем определителе 1-й и 2-й столбцы совпадают и определи- тель равен нулю (свойство 4).

1.5. Произведение матриц и его свойства

Определение 1.5.1. Матрицы А и В называются согласованны- ми, если число столбцов (ширина) первой матрицы равно числу строк (вы- соте) второй матрицы.

34

Свойства произведения матриц

6. (А + В) С = АС + ВС; 7. (ab) А = a (bА);

8. (aА) В = a (АВ);

9. $ А ¹ 0, В ¹ 0, такие, что А×В = 0.

35

Эта последняя матрица и есть обратная для А. Ее еще можно запи-

сать, вынеся за знак матрицы 1 4

36

Как правило, такой способ не применяют для вычисления обратной матрицы более высокого порядка, чем 2 или 3. Способ этот слишком гро- моздкий, требует много операций при вычислении определителей, даже на современных компьютерах потребуется много машинного времени, и на- копятся большие ошибки в вычислениях.

Рассмотрим способ, как вычислить обратную матрицу более эффек- тивным методом. Будем применять метод приведения матрицы с помощью элементарных преобразований к каноническому виду. Идея метода состоит в следующем.

Напишем исходную невырожденную матрицу А.

Припишем к ней за вертикальной чертой единичную матрицу Е того же порядка (А|Е).

Будем приводить матрицу А сначала к каноническому виду, причем преобразования будем совершать с «длинной» строкой, т.е. по тем же пра- вилам будем преобразовывать соответствующие строки матрицы Е. Про- делав конечное число шагов, приведем А к виду Р. Матрица Р в нашем случае будет треугольной с угловыми элементами, отличными от нуля, так как исходная матрица А – невырожденная.

Получим промежуточный результат (Р|В), где В – матрица, которая

37

1.7. Системы линейных уравнений. Основные понятия и определения

Определение 1.7.1. Системой линейных алгебраических уравне-

ний, содержащей m линейных уравнений и n неизвестных, называется система вида

тельно совпадает с числом неизвестных n.

В качестве такой системы уравнений можно предложить следующую задачу: поставщик продавал в течение 4-х кварталов года заказчику 4 вида товаров (вентиляторы или насосы различной мощности); аkl – это извест- ное количество (в условных единицах) товара l-ного вида (l = 1, 2, 3, 4), проданного в k-том квартале (k = 1, 2, 3, 4). Цена условной единицы каж- дого из 4-х видов товара x1, x2, x3, x4, неизвестна. Требуется найти зна- чение этих цен по известным значениям сумм квартальных платежей b1, b2, b3, b4. Математической моделью такой задачи будет система 4-х урав- нений с четырьмя неизвестными.

К системам линейных уравнений приводят, например, некоторые за- дачи теории цепей. Пусть схема электрической цепи с источником сину- соидального напряжения и тока одинаковой частоты содержит (q – 1) неза-

38

висимых узлов и l независимых контуров. Система уравнений, составлен- ных для данной схемы методом контурных токов, запишется следующим образом

z11ik1 + z12ik 2 + + z1likl = Ek1,

zl1ik1 + zl 2ik 2 + + zllikl = Ekl.

рядок m × n, состоит из коэффициентов при неизвестных в системе (1.7.1)

и называется матрицей системы уравнений. Вектор x = (х1, х2, ..., хn), ко-

ординатами которого являются неизвестные нашей системы, называется

вектором неизвестных, а вектор b = (b1, b2, ..., bm) – это вектор-столбец из свободных членов правых частей уравнений. Обычно оба эти вектора x и b записываются как векторы-столбцы

39

называется расширенной матрицей системы.

Ее последний столбец играет особую роль, поэтому будем отделять

нулю, то система (1.7.1) называется неоднородной системой. Если же в пра- вой части системы уравнения (1.7.1) стоят только нули (b1 = b2 = … = bm = 0), то систему называют однородной.

Определение 1.7.5. Вектор-столбец x , т.е. набор чисел х1, х2, ..., хn, называется решением системы уравнений, если при подстановке чисел

х1, х2, ..., хn в уравнения системы получаются верные равенства.

Решить систему уравнений (1.7.1) значит найти или описать все та- кие векторы-решения.

Определение 1.7.6. Две системы называются эквивалентными, если каждое решение первой является решением второй и каждое решение второй – решением первой. Системы, не имеющие решений, являются эк- вивалентными.

Определение 1.7.7. Элементарными преобразованиями системы линейных уравнений называют:

1.Перестановку уравнений в системе местами.

2.Умножение обоих частей уравнения на число отличное от нуля.

3.Перестановку местами соответствующих слагаемых в уравнении.

4.Прибавление к одному уравнению другого, умноженного на число отличное от нуля.

Замечание. Очевидно, элементарные преобразования решения сис- темы не меняют, т.е. полученная новая система остается эквивалентной прежней. Заметим, что эти операции над системой уравнений сводятся к

элементарным преобразованиям над расширенной матрицей A . Таким об-

разом, элементарные преобразования над системой (1.7.1) не меняют совокупности ее решений.

40