Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Системы линейных уравнений. Матрицы и определители

.pdf
Скачиваний:
67
Добавлен:
01.05.2014
Размер:
155.71 Кб
Скачать

Для матриц большего размера этот метод требует слишком больших усилий. Метод Гаусса гораздо эффективнее.

Пример. Найти матрицу, обратную к матрице

1

1

1

1

2

3

A = 1

3

6 :

Решение. Запишем составную матрицу (A|E) и пребразуем ее с помощью элементарных

преобразований строк в соответствии с методом Гаусса:

 

 

 

 

 

 

 

 

 

 

 

 

1

1

1

 

 

1

0

0

 

 

1

1

1

 

 

1

 

0

 

0

 

 

1

1

1

 

1

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(A

E) = 1 2 3

 

 

0 1 0 0 1 2

 

 

1 1 0

0 1 2

 

1

1 0

|

 

1 3 6

 

 

0 0 1 0 2 5

 

 

1 0 1

0 0

1

 

1 2 1

 

 

 

 

1

1

0

 

 

 

 

0

2

 

1

 

1

0

0

 

 

3

3

 

1

 

 

 

 

 

 

 

 

0 1 0

 

 

 

 

3

5

2

0

1 0

 

 

3

5

 

2

 

 

 

 

 

Ответ:

 

0 0 1

 

 

 

1 2

1

 

0

0 1

 

 

1 2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

3

 

1

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A1 =

 

3

 

5 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

1

 

 

 

 

 

 

 

 

 

 

 

Умножив полученную матрицу на A, нетрудно убедиться, что эти две матрицы действи-

тельно взаимно обратны.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Существует модификация этого метода:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 1 1

 

 

 

y1

 

1

1 1

 

 

y1

 

 

1 1 1

 

 

y1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(A

y) =

1 2 3

 

 

y2

0

1 2

−y1+y2

0 1 2

 

−y1 +y2

 

 

 

|

 

1 3 6

 

 

y3

 

0

2 5

−y1+y3

 

0 0 1

 

 

y12y2+y3

 

 

 

 

 

1 1 0

 

 

 

 

2y2 −y3

 

1 0 0

 

3y13y2 +y3

 

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 1 0

 

 

3y1+5y2

 

2y3

0 1 0

 

3y1+5y2

2y3

 

 

 

 

 

 

 

0 0 1

 

 

y12y2

+y3

 

0 0 1

 

y12y2+y3

 

 

 

 

 

Матрица коэффициентов при y1; y2; y3 есть A1. Результат совпадает с выписанным выше, разумеется.

Свойства операции обращения матриц. Отображение A 7→A1 называется операцией обращения. Она имеет следующие свойства:

1)(A1)1 = A;

2)(AT )1 = (A1)T ;

3)(AB)1 = B1A1.

Проверим для примера последнее свойство:

(AB)(AB)1 = (AB)(B1A1) = A(BB1)A1 = AEA1 = AA1 = E:

Здесь учтены свойство 3, ассоциативность умножения матриц и определение обратной матрицы.

21

5 Ранг матрицы

Рассмотрим произвольную матрицу

 

 

 

 

 

a11

a12

: : :

a1n

 

A =

a21

a22

: : :

a2n

:

 

 

 

 

 

 

 

a: : : : : :a: : : : ::::::

: : :a: : :

 

 

m1

m2

 

mn

 

Рангом матрицы называется максимальное число линейно независимых строк этой матрицы.

Введем обозначение для i-ой строки матрицы A:

Li = (ai1; ai2; : : : ; ain):

Напомним понятия линейной зависимости и линейной независимости, рассмотренные в §4.2. Составим уравнение

L1x1 + L2x2 + · · · + Lrxr = 0;

(5.1)

лежащее в основе определений. Если уравнение (5.1) может выполняться только при x1 = x2 = · · · = xr = 0, то строки L1; L2; : : : ; Lr называются линейно независимыми, в противном случае, когда некоторые из xi могут быть не равными нулю, – строки линейно зависимыми.

Теорема 1. Пусть строки L1; L2; : : : ; Lr линейно независимы, тогда пребразование строк

k

 

r

 

Li= SikLk; i = 1; 2; : : : ; r;

(5.2)

=1

 

с невырожденной матрицей S (det S ≠ 0) приводит к линейно независимой системе из

штрихованных строк.

 

 

 

L

x

 

+ Lx

 

+

 

+ L

x

 

= 0

 

 

 

 

 

Доказательство. Составляем уравнение

1

2

· · ·

r

или в краткой

 

1

 

2

 

 

r

 

 

 

форме

 

∑ ∑k

 

 

 

 

 

∑ ∑

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

r

r

 

 

 

 

 

r

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xiLi= xi

 

SikLk = Lk

 

 

xiSik = 0:

 

 

 

 

 

 

 

 

 

i=1

i=1

=1

 

 

 

 

 

k=1

i=1

 

 

 

 

 

 

 

 

 

 

 

 

Обозначив yk =

 

r

 

 

r

 

 

 

 

= 0. Так как строки Li линейно незави-

 

i=1 xiSik, получаем

 

k=1 Lkyk

 

теоремы, то все величины yk равны нуюлю, то есть

 

r

 

xiSik = 0;

k.

симы по условию

 

 

 

 

xS = 0

 

 

 

 

 

 

 

 

i=1

t

x

t

 

Переписав это уравнение в матричной форме

 

 

 

, или в более

 

 

 

 

 

 

 

 

 

заключаем, что строка x нулевая, так как det S ≠ 0. Таким образом, мы нашли, что из уравнения L1x1 +L2x2 +· · ·+Lrxr = 0 следует xi = 0; i. Следовательно, строки Li линейно независимы, что и требовалось.

Заметим, что исходную матрицу A можно рассматривать как столбец из ее строк, а формула (5.2) показывает, что операции над строками матрицы A – сложение, умножение на число и образование линейных комбинаций, – сводятся к умножению этой матрицы

слева на некоторую невырожденную матрицу:

 

 

 

 

 

 

 

 

c11

c12

: : : c1m

 

L1

 

 

 

c11L1 + c12L2 +

· · ·

+ c1mLm

 

 

SA =

 

c21

c22

: : : c2m

 

L2

=

 

c21L1 + c22L2 +

+ c2mLm

 

:

 

 

 

 

 

 

· · ·

 

 

 

 

 

 

 

 

 

 

· · ·

 

 

 

 

 

: : : : : : : : : : : : : : : : : : :

 

 

 

 

 

. . . . . . . . . . . . . . . . . . . . . . . . . . . . .

 

 

 

 

 

...

 

 

 

 

 

 

cm1

cm2 : : : cmm

Lm

 

cm1L1 + cm2L2 +

 

+ cmmLm

 

22

Отсюда нетрудно получить матрицы всех элементарных операций над строками. Например, диагональной матрице S соответствует умножение всех строк на числа. Если в диагональной матрице S все элементы, кроме i-го, равны 1, то этой матрице соответствует умножению i-ой строки на число. Если в матрице S главная диагональ состоит из единиц, а остальные элементы, кроме c21 равны нулю, то этой матрице соответствует прибавление ко второй строке первой строки, умноженной на число c21. Если матрица S имеет следующую блочную структуру

 

0

 

0

1

S = (0

E);

где

= (1

0);

а E – единичная матрица порядка m − 2, то умножение матрицы A слева на матрицу S приводит к перестановке первой и второй строк матрицы A. Нетрудно проверить, что определители всех матриц элементарных операций отличны от нуля.

Учитывая все изложенное, заключаем, что элементарные операции над независимыми строками матрицы не меняют ее ранга. Можно доказать большее: элементарные операции над всеми строками матрицы не меняют ее ранга.

Рассмотрение квадратных матриц приводит к следующему урверждению: определитель матрицы равен нулю в том и только том случае, когда строки этой матрицы линейно зависимы. Неизменность определителя при транспонировании позволяет утверждать, что определитель матрицы равен нулю в том и только том случае, когда столбцы этой матрицы линейно зависимы. А отсюда ясно, что строки квадратной матрицы линейно зависимы в том и только том случае, когда линейно зависимы ее столбцы.

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

Для практики самым важным является следующее утверждение: элементарные операции над строками и столбцами матрицы не меняют ее ранга. Оно позволяет вычислять ранг матрицы при помощи модифицированного алгоритма Гаусса. Суть метода в следующем: матрица приводится элементарными преобразованиями строк и столбцов к единичной матрице; число единиц на главной диагонали равно рангу матрицы.

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

 

 

 

 

 

 

 

 

1

2

3

2

 

A =

2 3

1

4

:

 

 

 

 

 

 

 

 

1

9

 

 

 

8

2

 

 

1

12

7

2

 

Решение. Первый шаг выполняем по методу Гаусса, а затем, путем прибавления первого столбца к последующим столбцам, уничтожаем элементы первой строки, начиная со второго элемента:

 

 

 

1

2

3

2

 

 

 

1

0

0

0

 

A

 

0 7

5

0

 

0 7

5

0

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 10

10

0

 

0 10

10

0

 

 

 

 

0

7

 

0

 

 

 

0

7

 

 

 

 

 

 

5

 

 

 

5

0

 

Далее вычеркнем вторую строку, которая пропорциональна третьей строке, и нулевой четвертый столбец. Затем прибавим ко второму столбцу третий столбец:

 

1

0

0

1

0

0

 

 

A=

0

7

5

0

12

5

:

 

0

10

10 0

0

10

 

23

На этом вычисления можно прекратить, так как мы получили квадратную матрицу 3-го порядка с ненулевым определителем. Ответ: ранг матрицы A равен 3.

Список литературы

[1]Смирнов В.И. Курс высшей математики. Т. 3, часть 1. М., Наука, 1974.

[2]Фаддеев Д.К. Лекции по алгебре. М., Наука, 1984.

[3]Мальцев А.И. Основы линейной алгебры. М., Наука, 1970.

[4]Шевцов Г.С. Линейная алгебра. Теория и прикладные аспекты. М., Финансы и статистика, 2003.

[5]Лунгу К.Н., Письменный Д.Т. и др. Сборник задач по высшей математике. 1 курс.

М., Айрис Пресс, 2007.

24