ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 07.09.2024
Просмотров: 184
Скачиваний: 0
тором этапе преобразований на месте матрицы A получится единичная мат рица, и тогда каждый столбец пристроенной матрицы будет представлять решение соответствующей подсистемы уравнений, т. е. на месте приписанной единичной матрицы появится обратная матрица. Схема обращения невырож денной матрицы A кратко может быть записана в виде:
(A | E)→(E | A−1 ) |
(3.4.5) |
Если же на некотором этапе процесса преобразований вспомогательной матрицы (3.4.4) на месте одной из строк матрицы A появится строка нулей, то это означает необратимость матрицы A.
Вычисление обратной матрицы в Microsoft Excel производится с помощью функции
A−1 = МОБР(матрица A),
где «матрица A» — ссылка на ячейки рабочего листа, содержащие данную матрицу. Данная формула должна быть введена в рабочий лист как ф о р м у л а м а с с и в а Microsoft Excel.
3.4.1. Найти (если это возможно) матрицу, обратную к матрице А из примера 1.3.1.
Решение. Припишем справа к матрице A единичную матрицу, и с помощью эле ментарных преобразований строк приведем матрицу (A | E) к такому виду, чтобы на месте матрицы А оказалась единичная матрица, тогда на месте единичной матрицы будет искомая матрица A−1 . Процесс элементарных преобразований иллюстрирует ся табл. 3.4.1.
|
|
|
|
|
|
|
Т а б л и ц а 3.4.1 |
|
|
|
|
(A | E) |
|
|
Примечания |
|
|
|
|
|
|
|
|
|
|
|
0 |
1 |
2 |
|
1 |
0 |
0 |
(1/4) II → I |
|
4 |
0 |
1 |
|
0 |
1 |
0 |
I → II |
|
3 |
–1 |
1 |
|
0 |
0 |
1 |
III → III |
|
1 |
0 |
1/4 |
|
0 |
1/4 |
0 |
I → I |
|
0 |
1 |
2 |
|
1 |
0 |
0 |
II → II |
|
3 |
–1 |
1 |
|
0 |
0 |
1 |
III – 3 I → III |
|
1 |
0 |
1/4 |
|
0 |
1/4 |
0 |
I → I |
|
0 |
1 |
2 |
|
1 |
0 |
0 |
II → II |
|
0 |
–1 |
1/4 |
|
0 |
–3/4 |
1 |
III + II → III |
|
1 |
0 |
1/4 |
|
0 |
1/4 |
0 |
I → I |
|
0 |
1 |
2 |
|
1 |
0 |
0 |
II → II |
|
0 |
0 |
9/4 |
|
1 |
–3/4 |
1 |
(1/3) III → III |
|
1 |
0 |
1/4 |
|
0 |
1/4 |
0 |
I – (1/4) III → I |
|
0 |
1 |
2 |
|
1 |
0 |
0 |
II – 2 III → II |
|
0 |
0 |
1 |
|
4/9 |
–1/3 |
4/9 |
III → III |
|
1 |
0 |
0 |
|
–1/9 |
1/3 |
–1/9 |
I – (1/4) III → I |
|
0 |
1 |
0 |
|
1/9 |
2/3 |
–8/9 |
II – 2 III → II |
|
0 |
0 |
1 |
|
4/9 |
–1/3 |
4/9 |
III → III |
|
В результате элементарных преобразований матрица
63
|
|
|
0 |
1 |
|
2 |
|
1 |
0 |
0 |
|
|||
|
|
|
|
|
|
|||||||||
|
|
|
|
4 0 |
1 |
|
0 |
1 |
0 |
|
|
|||
|
(A | E) = |
|
|
|
||||||||||
|
|
|
|
|
−1 |
|
|
|
|
|
|
|
||
приведена к виду |
|
|
3 |
1 |
|
0 |
0 |
1 |
|
|||||
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
0 |
|
0 |
|
−1/9 |
|
1/3 |
|
−1/9 |
|||
|
|
|
|
|
|
|||||||||
(E | A |
−1 |
|
1 |
|
0 |
|
1/9 |
|
|
2/3 |
|
−8/9 |
|
|
|
) = 0 |
|
|
|
|
|
. |
|||||||
|
|
|
|
|
|
|
|
|
|
|
−1/3 |
|
|
|
|
|
0 |
0 |
|
1 |
|
4/9 |
|
|
|
4/9 |
|||
|
|
|
|
|
|
|
Таким образом, получена искомая обратная матрица
|
|
−1/9 |
1/3 |
−1/9 |
||
A |
−1 |
|
1/9 |
2/3 |
−8/9 |
|
|
= |
. |
||||
|
|
|
|
−1/3 |
|
|
|
|
4/9 |
4/9 |
Сделаем проверку: по определению обратной матрицы должно выполняться ра
венство A−1A = AA−1 = E . В нашем случае |
|
|
|
|
|
|||
|
|
−1/9 |
3/9 |
−1/9 0 |
1 |
2 |
|
|
A |
−1 |
|
6/9 |
|
0 |
1 |
|
= |
|
A = 1/9 |
−8/9 4 |
|
|||||
|
|
|
−3/9 |
|
−1 |
1 |
|
|
|
|
4/9 |
4/9 3 |
|
|
(−1 0 +3 4 −1 3)/9 |
[−1 1+3 0 −1 (−1)]/9 |
(−1 2 +3 1−1 1)/9 |
1 |
0 |
0 |
|
|
||||
|
(1 0 |
+6 4 −8 3)/9 |
[1 1 +6 0 −8 (−1)]/9 |
(1 2 +6 1−8 1)/9 |
|
|
0 |
1 |
0 |
|
= E — |
= |
|
= |
|
||||||||
|
|
−3 4 +4 3)/9 |
[4 1 −3 0 +4 (−1)]/9 |
|
|
|
|
0 |
1 |
|
|
(4 0 |
(4 2 −3 1+4 1)/9 |
0 |
|
|
верно. Аналогично можно проверить, что AA−1 = E .
Теперь поясним, как получить тот же результат в пакете Microsoft Excel. Введем матрицу A в ячейки A2:C4 рабочего листа Microsoft Excel, как показано на рис. 3.4.1, а.
|
A |
B |
C |
D |
E |
F |
G |
H |
1 |
A |
|
|
|
A–1 |
|
|
|
2 |
0 |
1 |
2 |
|
=МОБР(A2:C4) |
|
|
|
3 |
4 |
0 |
1 |
|
|
|
|
|
4 |
3 |
–1 |
1 |
|
|
|
|
|
|
|
а) формула Microsoft Excel |
|
|
||||
|
A |
B |
C |
D |
E |
F |
G |
H |
1 |
A |
|
|
|
A–1 |
|
|
|
2 |
0 |
1 |
2 |
|
–0,11 |
0,33 |
–0,11 |
|
3 |
4 |
0 |
1 |
|
0,11 |
0,67 |
–0,89 |
|
4 |
3 |
–1 |
1 |
|
0,44 |
–0,33 |
0,44 |
|
б) результаты расчета
Рис. 3.4.1. Вычисление обратной матрицы в Microsoft Excel
64
Матрица A имеет две строки и три столбца, значит, матрица A–1 будет иметь три стро ки и два столбца. Отведем под результат ячейки E2:G4 (они как раз занимают три строки и два столбца). В ячейку E2 введем формулу «=МОБР(A2:C3)», причем эту фор мулу необходимо ввести как ф о р м у л у м а с с и в а. Для этого нужно мышью выделить диапазон E2:G4, начиная с ячейки E2, содержащей формулу, затем на жать клавишу <F2>, а затем — комбинацию клавиш <Ctrl> + <Shift> + <Enter>.
Результат представлен на рис. 3.4.1, б (в ячейках E2:G4). Замечаем, что результаты руч ных и компьютерных вычислений совпали (с точностью до ошибок округления). Если формулу ввести не как формулу массива, то будет рассчитан только левый верх ний элемент результата — число (–0,11).
В о п р о с ы д л я с а м о п р о в е р к и
1.Какая матрица называется обратной к данной квадратной матрице A?
2.Ко всем ли матрицам существуют обратные матрицы?
3.Может ли существовать несколько различных матриц, обратных к дан ной матрице?
4.Чему равен определитель матрицы, обратной к данной матрице A?
5.Всегда ли верно, что (AT)–1 = (A–1)T?
6.Всегда ли верно, что (A–1)–1 = A?
7.Всегда ли верно, что (AB)–1 = B–1A–1?
8.Как найти матрицу, обратную к данной, с помощью метода Жордана — Гаусса?
З а д а ч и д л я с а м о с т о я т е л ь н о г о р е ш е н и я
1. Найдите матрицу, обратную к матрице
1 0 A = .
1 1
2. Для данных матриц найдите обратные или докажите необратимость
2 |
0 |
|
|
|
12 4 7 |
|
5 −1 11 |
|||||
, |
|
−7 −2 |
|
|
, |
|
|
|
|
|||
|
2 |
|
|
4 |
|
|
4 |
1 7 |
. |
|||
0 |
|
|
|
1 −1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
1 |
−1 3 |
§ 3.5. ПРЯМЫЕ, ПЛОСКОСТИ И ЛИНЕЙНЫЕ НЕРАВЕНСТВА
В декартовой прямоугольной системе координат на плоскости о б щ е е
у р а в н е н и е п р я м о й имеет вид |
|
Ax1 +Bx2 +C =0 , |
(3.5.1) |
при этом вектор
A n = B
перпендикулярен данной прямой и носит название нормального вектора.
65
Если B ≠0 , то из общего уравнения прямой (3.5.1) c помощью замены
k = − |
A |
, |
b = − |
C |
|
B |
|||
|
B |
|
получается у р а в н е н и е п р я м о й с у г л о в ы м к о э ф ф и ц и е н т о м: y =kx +b .
Здесь угловой коэффициент k =tg α равен тангенсу угла наклона прямой с положительным направлением оси абсцисс. Уравнение с угловым коэффици ентом не описывает прямые, параллельные оси ординат.
Уравнение прямой, проходящей через данную точку P(x0 , y0 ) и имеющую угловой коэффициент k , таково:
y −y0 =k(x −x0 )
Г е о м е т р и ч е с к и й с м ы с л л и н е й н ы х н е р а в е н с т в определя ется следующими теоремами.
ТЕОРЕМА 3.5.1. Пусть прямая задана общим уравнением Ax1 +Bx2 +C =0 . Если отложить нормальный вектор от любой точки этой прямой, то ко нец отложенного вектора будет находиться в положительной полуплоско сти от данной прямой, т. е. A(x1 + A) +B(x2 +B) +C >0 .
Доказательство. Доказательство очень короткое:
A(x1 + A) +B(x2 +B) +C = Ax1 +Bx2 +C + A2 +B2 =0 + A2 +B2 = A2 +B2 >0 , так как A и B не могут быть равны нулю одновременно.
ТЕОРЕМА 3.5.2. Всякая прямая Ax1 +Bx2 +C =0 разбивает плоскость на две полуплоскости. Для координат всех точек, лежащих в одной полуплоско сти, выполняется неравенство Ax1 +Bx2 +C >0 , для координат всех точек
другой полуплоскости справедливо |
противоположное |
неравенство |
Ax1 +Bx2 +C <0 . |
|
|
Доказательство этой теоремы мы оставляем читателю. |
|
|
В трехмерном пространстве о б щ е е |
у р а в н е н и е п л о с к о с т и имеет |
|
вид: |
|
|
Ax1 +Bx2 +Cx3 +D =0 , |
(3.5.2) |
|
при этом вектор |
|
|
A |
|
|
|
|
|
n = B |
|
|
|
|
|
С |
|
перпендикулярен данной плоскости и носит название нормального вектора. Плоскость (3.5.2) делит трехмерное пространство на два полупространства,
в одном из которых выполняется неравенство Ax1 +Bx2 +Cx3 +D >0 , а в дру гом — неравенство Ax1 +Bx2 +Cx3 +D <0 .
Система m линейных алгебраических неравенств с n неизвестными мо жет быть записана в виде
66
n |
-bi , i =1,2,…p, |
∑aijxj |
|
j=1 |
(3.5.3) |
n |
∑aijxj .bi , i =p +1,p +2,…m,
j=1
Совокупность n чисел α1,α2,…,αn , взятых в определенном порядке, называ ется решением системы неравенств (3.5.3), если при подстановке этих чисел на место соответствующих неизвестных неравенства системы не нарушаются. Решение (α1,α2,…,αn ) системы неравенств называется неотрицательным, если все αj >0 . Определения совместности, несовместности, определенно
сти, неопределенности и эквивалентности систем линейных алгебраиче ских неравенств формулируются точно так же, как и соответствующие опре деления для систем линейных алгебраических уравнений.
Для решения системы линейных алгебраических неравенств вида (3.5.3) эту
систему путем |
введения дополнительных неотрицательных |
неизвестных |
xn+1, xn+2,…, xn+m |
преобразовывают в систему линейных алгебраических урав |
|
нений: |
|
|
|
n |
|
|
∑aijxj +xn+i =bi , i =1,2,…p, |
|
|
j=1 |
(3.5.4) |
|
n |
|
|
∑aljxj −xn+l =bl , l =p +1,p +2,…m. |
|
|
l=1 |
|
Каждому решению (α1,α2,…,αn ,αn+1,αn+2,…,αn+m ) системы линейных урав нений (3.5.4), где последние m компонент αn+1,αn+2,…,αn+m
можно поставить в соответствие вполне определенное решение (α1,α2,…,αn ) системы линейных неравенств (3.5.3), и наоборот. В этом смысле говорят, что система уравнений (3.5.4) заменяет систему неравенств (3.5.3). Исследование и решение системы m линейных неравенств с n неизвестными сводится к ис следованию и решению соответствующей системы m линейных уравнений с (n +m) неизвестными. В частности, вопрос о нахождении неотрицательных решений системы линейных неравенств (3.5.3) сводится к вопросу о нахожде нии .неотрицательных решений системы линейных уравнений (3.5.4).
Решение системы линейных неравенств удобно иллюстрировать гео метрически в случае двух и трех неизвестных. Например, каждое неравенство системы с двумя неизвестными определяет полуплоскость вместе с ограничи вающей прямой, а вся система — общую часть таких полуплоскостей, кото рая, очевидно, будет выпуклым многоугольником, называемым многоугольни ком решений данной системы неравенств. Важнейшее свойство выпуклого многоугольника состоит в том, что он вместе с любыми двумя своими точками содержит весь соединяющий их отрезок. По аналогии, каждое неравенство системы с тремя неизвестными определяет полупространство вместе с гра ничной плоскостью, а множество решений такой системы, в случае ее совме стности, образует выпуклый многогранник.
67