Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
М / Линейная алгебра,ПК,ПКЛК,по направлению Экономика, Менеджмент 2012.pdf
Скачиваний:
122
Добавлен:
15.02.2016
Размер:
1.44 Mб
Скачать

Решение. Данный определитель может быть вычислен непосредственно, используя теорему Лапласа. Однако, как нетрудно заметить, первая и третья строки матрицы – пропорциональны (элементы третьей строки в два раза больше соответствующих элементов первой строки). Поэтому, используя третье свойство определителей, заключаем, что значение его равно нулю:

A = 0

Тема 3. Обратная матрица.

Основные вопросы темы

1.Понятие обратной матрицы.

2.Необходимое и достаточное условие существования обратной матрицы.

3.Алгоритм вычисления обратной матрицы.

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

ОПРЕДЕЛЕНИЕ. Обратной, по отношению к квадратной матрице А, называется матрица, обозначаемая как A1 , удовлетворяющая соотношению:

A A1 = A1 A = E

3.2 Необходимое и достаточное условие существования обратной матрицы.

Не всякая матрица имеет обратную. Условия существования обратной матрицы определены следующей теоремой.

Теорема (необходимое и достаточное условие существования обратной матрицы). Обратная матрица существует и единственна тогда и только тогда, когда исходная матрица является невырожденной.

3.3 Алгоритм вычисления обратной матрицы

Алгоритм вычисления обратной матрицы включает 5 шагов.

Шаг 1 (проверка существования обратной матрицы). Вычисление определителя исходной матрицы А. Если матрица А – невырожденная

21

(т.е. ее определитель отличен от нуля), то обратная матрица существует.

Шаг 2. Вычисление транспонированной матрицы АТ.

Шаг 3. Вычисление алгебраических дополнений - АТij всех элемен-

тов транспонированной матрицы АТ и запись их в в матрицу ~ = (~ ),

A aij

называемую присоединенной ( aij = Aij ).

~

T

Шаг 4. Вычисление обратной матрицы А-1 по формуле:

 

1

 

 

1

~

A

 

=

 

 

 

 

 

A .

 

 

 

A

 

 

 

 

 

Шаг 5. Проверка правильности вычисления обратной матрицы, используя соотношение:

A A1 = E .

 

 

 

 

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

3

2

 

À =

 

 

.

 

 

5

4

 

 

 

 

Решение.

Шаг 1. (проверка существования обратной матрицы). Вычисляем определитель исходной матрицы А:

À =

3

2

= 3 4 2 5 = 2 0 .

 

5

4

 

Определитель матрицы отличен от нуля. Следовательно, обратная матрица существует.

Шаг 2. Транспонируем исходную матрицу А (вычисляем транспонированную матрицу АТ):

À

Ò

3

5

 

=

 

 

 

 

 

 

2

4

 

 

 

 

 

Шаг 3. Находим алгебраические дополнения - АТij всех элементов

транспонированной матрицы АТ и запись их в в матрицу ~ = (~ ), назы-

A aij

ваемую присоединенной ( aij = Aij ):

 

 

 

 

 

 

 

 

 

 

 

~

 

T

 

 

 

 

 

 

 

 

 

 

~

 

Ò

 

 

 

1+1

Ì

 

Ò

= 1 4 = 4

à11

= À11

= (1)

 

 

 

11

~

 

Ò

 

 

 

1+ 2

Ì

Ò

 

 

 

à12

= À12

= (1)

 

 

12= (1) 2 = −2

~

=

Ò

=

(1)

2+1

 

Ì

Ò

 

= (1) 5 = −5

à21

À21

 

 

 

21

 

~

 

 

Ò

= (1)

2+

2

Ì

 

Ò

= 1 3 = 3

à22

= À22

 

 

 

 

22

 

 

 

 

~

 

4

 

2

 

 

 

 

 

À

=

 

 

 

 

 

 

.

 

 

 

 

 

 

5

 

3

 

 

 

 

 

 

 

 

 

 

 

Шаг 4. Вычисляем обратную матрицу А-1 по формуле:

22

 

1

 

 

1 ~

1

4

2

2

1

A

 

=

 

 

A =

 

 

 

 

 

=

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

A

 

2

 

5

3

 

 

2,5

1,5

 

 

 

 

 

 

 

 

 

 

Шаг 5. Осуществляем проверку правильности вычисления обрат-

ной матрицы, используя соотношение A A1

= E :

 

 

 

 

 

 

A A

1

 

3

2

 

 

2

1

 

 

3

2

+ 2 (2,5)

3

(1)+ 2 1,5

 

 

1

0

 

= E .

 

=

 

 

 

 

 

 

 

=

 

 

 

 

 

 

=

 

 

 

 

 

 

5

4

 

 

2,5

1,5

 

 

5

2

+ 4 (2,5)

5

(1)+ 4 1,5

 

 

0

1

 

 

 

 

 

 

 

 

 

 

 

 

 

Условие A A1

= E - выполнено, следовательно, обратная матрица

найдена верно.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3.4 Ранг матрицы.

В теории матричных вычислений широко используется поня-

тие ранга матрицы. Определим это понятие.

Рассмотрим произвольную матрицу размером m×n - À . В общем

m×n

случае эта матрица не является квадратной, поэтому, понятие определителя для нее не существует. Однако, вычеркиванием отдельных ее строк или столбцов можно вычленить квадратные подматрицы k–ого порядка (k ≤ min(m,n)), для которых значение определителя существует. Эти определители получили название минора k–ого порядка мат-

рицы А. Например, для матрицы À существуют миноры третьего,

4×3

второго и первого порядка. Среди множества определителей (миноров) различных квадратных подматриц исходной матрицы А могут встречаться как равные, так и неравные нулю.

ОПРЕДЕЛЕНИЕ. Рангом матрицы А называется наивысший порядок отличных от нуля миноров этой матрицы и обозначается как rang(A) или r(A).

Свойства ранга матрицы.

1. Значение ранга матрицы À не превосходит меньшего из ее

m×n

размеров: rang( À ) ≤ min(m,n).

m×n

2. Значение ранга матрицы равно нулю тогда и только тогда, когда исходная матрица является нулевой: rang(A)=0 Ù A=O.

3. Для квадратной матрицы À n-го порядка значение ее ранга

n×n

равно n (rang(A)= n) тогда и только тогда, матрица А не является вырожденной.

23

При операциях сложения и умножения матриц, ранг результирующей матрицы удовлетворяет следующим соотношениям:

1)rang(A+B) ≤ rang(A) + rang (B)

2)rang(A+B) ≥ rang(A) - rang (B)

3)rang(A·B) ≤ min [rang(A); rang (B)]

4)rang(A·B)= rang(A), если В – квадратная невырожденная матрица.

5)rang(A·B) ≥ rang(A)+ rang (B) – n, где n - число столбцов мат-

рицы А (строк матрицы В)

6)rang(AT·A) = rang(A)

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

 

2

3

6

2

 

 

0

0

0

0

 

 

 

B =

3

1

4

2

.

 

 

 

8

2

3

3

 

 

 

Решение. Порядок матрицы В равен четырем. В силу указанного выше свойства 1, ее ранг не превосходит это значение: rang(В) ≤ 4. Однако ее ранг не может быть равным четырем, поскольку (как не трудно заметить) матрица В является вырожденной (ее определитель равен нулю, т.к. содержит нулевую строку). Поэтому: rang(В) ≤ 3. Вычеркнем эту нулевую строку и, например, четвертый столбец. Получим подматрицу B* третьего порядка:

 

2

3

6

 

 

3

1

4

 

B* =

.

 

8

2

 

 

 

3

Вычислим ее определитель: Â = 47 0 . Поскольку определитель отли-

чен от нуля, а определитель более высокого порядка равен нулю, делаем вывод, что ранг исходной матрицы В равен трем: rang(В) = 3.

В случае матриц больших размеров, число возможных подматриц – достаточно велико. Поэтому, непосредственное вычисление ранга матрицы может оказаться весьма трудоемким процессом. С целью минимизации вычислительных затрат использую специальные преобразования исходной матрицы, не изменяющие ранг исходной матрицы, но облегчающие его нахождение. Такие преобразования получили название элементарных преобразований. К ним относятся:

1. Отбрасывание нулевых строк и столбцов.

24

2.Умножение всех элементов строки (столбца) на число, не равное нулю.

3.Перестановка строк и столбцов.

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

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

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

Линейная зависимость строк и столбцов матрицы.

Рассмотрим матрицу общего вида:

A =

m×n

a

11a21

a31

Lai1

Lam1

a12 a13 a21 a21 a32 a33

LL

ai2 ai3

L L am2 am3

Kaij

La2 j

La3 j

LL

Laij

LL

Lamj

K a1n

 

 

L a2n

L a

 

 

3n

L L .

 

 

L ain

L L

L amn

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

e1 = (a11 a12 a13 … a1n) e2 = (a21 a22 a23 … a2n)

em = (am1 am2 am3 … amn)

Каждая из этих строк является матрицей – строкой. Поэтому к ним применимы все рассмотренные ранее свойства матриц. В частности:

а) строки равны, если их элементы, находящиеся в одинаковых позициях - совпадают;

б) для умножения строки на некоторое число λ необходимо умножить все элементы строки на это число:

λek = (λak1 λak2 λak3 … λakn) ;

в) для сложения (вычитания) двух строк необходимо сложить (вычесть) их элементы, находящиеся в одинаковых позициях; и др.

ОПРЕДЕЛЕНИЕ. Строка е называется линейной комбинацией строк e1, e2 , …, em , если она может быть представлена в виде

25