I. Обращение верхней треугольной матрицы

Практическое занятие 3. Операции с треугольными матрицами

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

Пусть А – верхняя треугольная матрица четвертого порядка

А = I. Обращение верхней треугольной матрицы - student2.ru .

а искомая обратная ей матрица, элементы которой I. Обращение верхней треугольной матрицы - student2.ru подлежат определению, запишется в виде

I. Обращение верхней треугольной матрицы - student2.ru = I. Обращение верхней треугольной матрицы - student2.ru .

По определению обратной матрицы должно выполняться равенство

I. Обращение верхней треугольной матрицы - student2.ru

I. Обращение верхней треугольной матрицы - student2.ru · I. Обращение верхней треугольной матрицы - student2.ru = I. Обращение верхней треугольной матрицы - student2.ru .

Используя правило умножения матриц, перемножим матрицы в левой части равенства

I. Обращение верхней треугольной матрицы - student2.ru =

= I. Обращение верхней треугольной матрицы - student2.ru .

Исходя из определения равенства двух матриц, для отыскивания неизвестных величин получаем уравнения, сравнив соответствующие элементы первых строк.

I. Обращение верхней треугольной матрицы - student2.ru = 1; (1)
I. Обращение верхней треугольной матрицы - student2.ru = 0; (2)
I. Обращение верхней треугольной матрицы - student2.ru = 0; (3)
I. Обращение верхней треугольной матрицы - student2.ru = 0. (4)

Из этих уравнений следует:

из (1) I. Обращение верхней треугольной матрицы - student2.ru ; (5)
из (2) I. Обращение верхней треугольной матрицы - student2.ru ; (6)
из (3) I. Обращение верхней треугольной матрицы - student2.ru (7)
из (4) I. Обращение верхней треугольной матрицы - student2.ru (8)

Проделав аналогичную работу для второй, третьей и четвертой строк, получим для второй строки

I. Обращение верхней треугольной матрицы - student2.ru = 0; (9)
I. Обращение верхней треугольной матрицы - student2.ru = 1; (10)
I. Обращение верхней треугольной матрицы - student2.ru = 0; (11)
I. Обращение верхней треугольной матрицы - student2.ru = 0 (12)
из (9) I. Обращение верхней треугольной матрицы - student2.ru = 0. (13)

С учетом, что I. Обращение верхней треугольной матрицы - student2.ru = 0, получаем:

из (10) I. Обращение верхней треугольной матрицы - student2.ru ; (14)
из (11) I. Обращение верхней треугольной матрицы - student2.ru ; (15)
из (12) I. Обращение верхней треугольной матрицы - student2.ru (16)

Для третьей строки

I. Обращение верхней треугольной матрицы - student2.ru = 0; (17)
I. Обращение верхней треугольной матрицы - student2.ru = 0; (18)
I. Обращение верхней треугольной матрицы - student2.ru = 1; (19)
I. Обращение верхней треугольной матрицы - student2.ru = 0. (20)

Из этих уравнений следует:

из (17) I. Обращение верхней треугольной матрицы - student2.ru = 0. (21)

Учитывая, что I. Обращение верхней треугольной матрицы - student2.ru = 0, из (18) получаем

I. Обращение верхней треугольной матрицы - student2.ru = 0, (22)

а учитывая, что и I. Обращение верхней треугольной матрицы - student2.ru = 0, и I. Обращение верхней треугольной матрицы - student2.ru = 0, находим:

из (19) I. Обращение верхней треугольной матрицы - student2.ru ; (23)
из (20) I. Обращение верхней треугольной матрицы - student2.ru ; (24)

И наконец, для четвертой строки ищем

I. Обращение верхней треугольной матрицы - student2.ru = 0; (25)
I. Обращение верхней треугольной матрицы - student2.ru = 0; (26)
I. Обращение верхней треугольной матрицы - student2.ru = 0; (27)
I. Обращение верхней треугольной матрицы - student2.ru = 1. (28)

Из (25), (26) и (27) следует, что

I. Обращение верхней треугольной матрицы - student2.ru = 0; (29)
I. Обращение верхней треугольной матрицы - student2.ru = 0; (30)
I. Обращение верхней треугольной матрицы - student2.ru = 0, (31)
а с учетом этого и (28) получаем  
I. Обращение верхней треугольной матрицы - student2.ru . (32)

Равенства (9), (21), (22), (29), (30) и (31) показывают, что равны нулю те элементы обратной матрицы I. Обращение верхней треугольной матрицы - student2.ru , у которых первый индекс i больше второго индекса j, т.е. если i > j, то

I. Обращение верхней треугольной матрицы - student2.ru . (3.1)

Из равенства (5), (14), (23), (32) диагональные элементы обратной матрицы I. Обращение верхней треугольной матрицы - student2.ru , у которой первый и второй индексы равны (i=j), определяются так:

I. Обращение верхней треугольной матрицы - student2.ru ; I. Обращение верхней треугольной матрицы - student2.ru ; I. Обращение верхней треугольной матрицы - student2.ru ; I. Обращение верхней треугольной матрицы - student2.ru ,

что можно объединить одной записью: если i = j, то

I. Обращение верхней треугольной матрицы - student2.ru . (3.2)

По формулам (6), (7), (8), (15), (16), (24) и (32) определяются те элементы обратной матрицы для верхней треугольной, у которых первый индекс i меньше второго индекса j (i<j), т.е. элементы, стоящие над главной диагональю. Полученная по этим формулам обратная матрица будет также верхней треугольной.

Когда верхняя треугольная матрица имеет порядок n, элементы обратной ей матрицы находятся по аналогичным формулам, которые имеют следующий вид:

если i = j, то I. Обращение верхней треугольной матрицы - student2.ru ; (3.3)
если i > j, то I. Обращение верхней треугольной матрицы - student2.ru ; (3.4)
если i < j, то I. Обращение верхней треугольной матрицы - student2.ru ; (3.5)

Например, по формуле (3.5) элемент обратной матрицы пятого порядка

I. Обращение верхней треугольной матрицы - student2.ru . (3.6)

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

1. Определить диагональные элементы обращенной матрицы по формуле (3.3)

2. После этого подписать матрицы одну под другой.

3. На места элементов, стоящих ниже главной, вписать нули.

4. Чтобы определить элемент, стоящий над главной диагональю обратной матрицы, надо составить алгебраическую сумму произведений элементов, стоящих в обратной матрице I. Обращение верхней треугольной матрицы - student2.ru , левее определяемого, на соответствующие элементы того столбца матрицы А (т.е. той матрицы, для которой ищется обратная), в котором стоит определяемый элемент. Эту алгебраическую сумму надо разделить на диагональный элемент матрицы А, стоящий в том же столбце, что и определяемый элемент. Определяемый элемент равен этому частному, взятому с обратным знаком.

По этому правилу выражение в скобках в (3.6) получается так: искомый элемент I. Обращение верхней треугольной матрицы - student2.ru матрицы I. Обращение верхней треугольной матрицы - student2.ru находится в первой строке и пятом столбце. Перед ним в первой строке обратной матрицы I. Обращение верхней треугольной матрицы - student2.ru стоят элементы I. Обращение верхней треугольной матрицы - student2.ru , I. Обращение верхней треугольной матрицы - student2.ru , I. Обращение верхней треугольной матрицы - student2.ru и I. Обращение верхней треугольной матрицы - student2.ru , а в пятом столбце матрицы А – элементы I. Обращение верхней треугольной матрицы - student2.ru , I. Обращение верхней треугольной матрицы - student2.ru , I. Обращение верхней треугольной матрицы - student2.ru и I. Обращение верхней треугольной матрицы - student2.ru . Составляется алгебраическая сумма произведений первого элемента в первой строке матрицы I. Обращение верхней треугольной матрицы - student2.ru на второй элемент в пятом столбце матрицы А и т.д. Для уяснения этого правила решим несколько задач.

Задача 3.1 Найти обратную матрицу для матрицы

А = I. Обращение верхней треугольной матрицы - student2.ru .

Решение.п. 1.Находим диагональные элементы I. Обращение верхней треугольной матрицы - student2.ru , I. Обращение верхней треугольной матрицы - student2.ru , I. Обращение верхней треугольной матрицы - student2.ru и I. Обращение верхней треугольной матрицы - student2.ru обращенной матрицы по формуле (3.3)

I. Обращение верхней треугольной матрицы - student2.ru ; I. Обращение верхней треугольной матрицы - student2.ru ; I. Обращение верхней треугольной матрицы - student2.ru ; I. Обращение верхней треугольной матрицы - student2.ru .

пп. 2 и 3. Подписываем матрицы одну под другой и вписываем нули на места элементов, стоящих под главной диагональю

А = I. Обращение верхней треугольной матрицы - student2.ru

I. Обращение верхней треугольной матрицы - student2.ru

п. 4. Определяем элементы I. Обращение верхней треугольной матрицы - student2.ru , стоящие над главной диагональю (i < j), по формуле (3.5), пользуясь указанным в п.4 правилом.

I. Обращение верхней треугольной матрицы - student2.ru ;

I. Обращение верхней треугольной матрицы - student2.ru = I. Обращение верхней треугольной матрицы - student2.ru

Для проверки надо перемножить матрицы А и I. Обращение верхней треугольной матрицы - student2.ru и убедится, что получится единичная матрица

I. Обращение верхней треугольной матрицы - student2.ru = I. Обращение верхней треугольной матрицы - student2.ru .

При решении этой задачи все элементы обратной матрицы были представлены в виде простых дробей для облегчения контроля. На практике же все вычисления ведутся в десятичных дробях.

Наши рекомендации