Линейная алгебра для экономистов. Учебное пособие
.pdf4. С к а л я р н о е п р о и з в е д е н и е в е к т о р о в одинаковой размерности, обозначаемое круглыми скобками:
X
(x; y) = xiyi:
Например,
00 1 0 11
11
@@ 1 A; @ 2 AA = 1 ( 1) + ( 1) 2 + 3 1 = ( 1) + ( 2) + 3 = 0:
31
5. П р о и з в е д е н и е |
д в у х |
м а т р и ц согласованной |
|||
размерности: |
|
|
|
1 |
|
0 (a2; B1) |
(a2; B2) |
: : : (a2 |
; Bn) |
|
|
(a1; B1) |
(a1; B2) |
: : : (a1 |
; Bn) |
C |
|
C = A B = B ... |
... |
... ... |
|
: |
|
B |
|
|
|
C |
|
B (am; B1) |
(am; B2) : : : (am; Bn) |
C |
|
||
@ |
|
|
|
A |
|
Результат представляет собой таблицу из всевозможных скалярных произведений векторов-строк первого множителя на векторы-столбцы второго множителя. Чтобы эти скалярные произведения имели смысл, первый матричный множитель должен иметь столько же столбцов, сколько второй строк. Напротив, число строк результирующей матрицы совпадает с числом строк первого множителя, а число ее столбцов с числом столбцов второго множителя.
Например,
3 0 1 |
0 |
1 |
4 |
1 |
= |
3 |
1 + 0 |
2 + 1 |
1 |
3 |
4 + 0 |
0 + 1 |
1 |
= |
|
|
|
||||||||||||||
|
1 |
1 |
|
||||||||||||
1 1 2 |
|
@ |
2 |
0 |
A |
|
1 |
1 + 1 |
2 + 2 |
1 |
1 |
4 + 1 |
0 + 2 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
||||||||
5 6
=:
4 13
Замечание 1.1. В силу требования согласованности размерности менять местами матричные множители в общем случае нельзя.
Замечание 1.2. Перечисленные операции неявно вводят в ы ч и - т а н и е матриц одинаковой размерности по естественному правилу
Cm n = Am n Bm n = Am n + ( 1)Bm n = (aij bij)m n:
11
Например, |
|
|
3 1 |
5 2 |
1 1 |
|
|
|
|
3 5 1 |
|
1 2 1 |
2 3 0 |
||||||
4 0 2 |
|
2 0 1 |
= |
4 2 |
0 0 |
2 1 |
= |
2 0 1 |
: |
Перечисленные выше арифметические операции над матрицами обладают следующими легко проверяемыми свойствами2:
1: A + B = B + A; |
9: |
(AB)C = A(BC); |
|
2: (A + B) + C = A + (B + C); |
10: |
AB + AC = A(B + C); |
|
3: |
(A + B)> = A> + B>; |
11: |
AB + CB = (A + C)B; |
4: ( A)( B) = ( )(AB); |
12: |
(AB)> = B>A>; |
|
5: ( A) = ( )A; |
13: |
A + 0 = 0 + A = A; |
|
6: |
( A)> = A>; |
14: |
A 0 = 0 A = 0; |
7: A + B = (A + B); |
15: |
A E = E A = A; |
|
8: |
A + A = ( + )A; |
16: |
1 A = A: |
Убедитесь в справедливости этих свойств самостоятельно. Для примера приведем доказательство свойства 9. Предположим, что размерности всех сомножителей согласованы. Обозначим F = (AB)C, H = A(BC). Совпадение размерности этих матриц очевидно. Элементы первой из них равны
01
X X XX
fij = @ aisbskAckj = aisbskckj:
(k) (s) (k) (s)
Элементы второй вычисляются как
01
X X XX
hij = @ais bskckjA = aisbskckj:
(s) (k) (s) (k)
Поскольку порядок суммирования чисел можно менять, то
fij = X(k) |
X(s) |
aisbskckj = X(s) |
X(k) |
aisbskckj = hij; |
что и требовалось. |
|
|
|
|
2На самом деле это привычные всем переместительный, сочетательный и распределительный законы, а также ряд других. Исключение составляет умножение матриц, где в общем случае AB 6= BA.
12
1.4. Элементарные преобразования матриц
Применение перечисленных выше операций к специальным матрицам обладает рядом важных особенностей.
Так, д и а г о н а л ь н ы е матрицы D = diag( 1; : : : ; s) подходящей размерности, используемые в качестве левых или правых множителей, реализуют соответственно масштабирование строк или столбцов других матриц, а именно:
D Am n = 0 2a2 |
1; Am n D = 1A1 2A2 : : : nAn : |
||||||||||
|
B |
|
1a1 |
|
C |
|
|
|
|
||
|
|
|
|
|
j j j |
||||||
|
: : : |
|
|||||||||
|
|
|
a |
m |
|
|
|
|
|
||
|
B |
m |
|
C |
|
|
|
|
|||
|
@ |
|
|
|
A |
|
|
|
|
||
Умножение произвольной матрицы на п е р е с т а н о в о ч н у ю слева или справа приводит соответственно к перестановке строк или столбцов исходной матрицы в том порядке, в каком следуют единицы перестановочной матрицы. Например,
0 1 0 0 1A3 n = |
0 a1 |
1; |
Bm 3 |
0 1 0 0 1 |
= (B2 B1 B3) : |
||
0 1 0 |
|
a2 |
|
|
|
0 1 0 |
|
@ 0 0 1 A |
@ |
|
|
A |
|
@ 0 0 1 A |
j j |
a3 |
|
|
|||||
Н и ж н е - и в е р х н е т р е у г о л ь н ы е матрицы интересны тем, что при перемножении снова дают нижне- и верхнетреугольные матрицы. Кроме того, они раскладываются в произведения элементарных матриц, например:
L = |
0 a 1 |
0 0 |
1 = |
0 a 1 0 0 10 0 1 0 0 10 0 1 |
0 0 1 = |
|||||||||||||||||||||||
|
|
1 |
0 |
0 |
0 |
|
|
1 |
0 |
0 |
0 |
CB |
1 |
0 |
0 |
0 |
CB |
1 |
0 |
0 |
0 |
C |
|
|||||
|
B d e f 1 |
C B d 0 0 1 |
0 |
e |
0 1 |
0 0 |
f 1 |
|
||||||||||||||||||||
|
B |
b c |
1 0 |
C |
B |
b 0 1 0 |
CB |
0 |
c |
1 0 |
CB |
0 0 |
1 0 |
C |
|
|||||||||||||
|
@ |
|
|
|
|
A |
@ |
|
|
|
|
A@ |
|
|
|
|
|
A@ |
|
|
|
|
|
A |
|
|||
|
|
|
|
|
|
= 0 a 1 0 0 |
10 0 1 0 0 |
10 0 |
1 0 |
0 1 |
: |
|||||||||||||||||
|
|
|
|
|
|
|
B |
1 |
0 |
0 |
0 |
CB |
1 |
0 0 |
0 |
|
|
1 |
0 |
0 |
0 |
C |
|
|||||
|
|
|
|
|
|
|
0 0 0 1 |
0 0 0 1 |
CB d e f 1 |
|
||||||||||||||||||
|
|
|
|
|
|
|
B |
0 0 1 0 |
CB |
b c 1 0 |
CB |
0 |
0 1 |
0 |
C |
|
||||||||||||
|
|
|
|
|
|
|
@ |
|
|
|
|
|
A@ |
|
|
|
|
|
A@ |
|
|
|
|
|
A |
|
||
Умножение произвольной матрицы на в е к т о р - с т р о к у слева
или на |
в е к т о р - с т о л б е ц справа дает соответственно л и н е й - |
н у ю |
к о м б и н а ц и ю ее строк или столбцов, в силу чего имеем |
13
равенства
m
X x>A = xiai;
i=1
где
|
0 |
a1 |
1 |
= 0 A1 |
||
|
|
|
|
|||
A = |
|
|
a2 |
|
||
B |
: : : |
C |
||||
|
a |
m |
|
|||
|
B |
|
C @ |
|||
|
@ |
|
|
A |
|
|
В частности, для
x>y = (x; y) =
Однако
0
u1v1 u1v2
uv> = B u2v1 u2v2
B
@ : : : : : :
umv1 umv2
n
X
yiAi;
i=1
1
: : : An A = (aij)m n:
xiyi:
Отсюда для в е к т о р о в - о р т о в выполняются равенства
|
0 |
0 : : : 0 |
w1 |
0 : : : 0 |
1 |
|
|
Aek = Ak; ek>A = ak; wek> = |
0 : : : 0 |
w2 |
0 : : : 0 |
k |
; |
||
B |
: : : : : : : : : : : : : : : : |
C |
; (ej; ek) = j |
||||
|
B |
0 : : : 0 |
wm 0 : : : 0 |
C |
|
|
|
|
@ |
|
" |
|
A |
|
|
k-й столбец
где jk символ Кронекера (равен 1 при k = j и 0 в противном случае).
Для элементарной матрицы имеем представление
00 1 : : : w2 |
: : : 01 |
|
|
1 0 : : : w1 |
: : : 0 |
C |
|
T = B ... ... ... |
... |
= E + wek>; |
|
B0 0 : : : wn : : : 1C |
|
||
B |
|
C |
|
@ |
|
A |
|
14
где w вектор-столбец, отличающий матрицу T от единичной и расположенный в ее k-й позиции (при этом tkk = 1 + wk).
Теорема 1.1. Умножение матрицы Am n слева на элементарную матрицу Tm m = E + we>k приводит к тому, что k-я строка исходной матрицы последовательно добавляется ко всем другим ее строкам в промасштабированном виде, а в случае wk 6= 0 и сама к себе.
Доказательство. В самом деле, |
|
|
0 |
|
1 |
|
||
|
|
|
|
w2 ak |
|
|||
|
|
|
|
= A + B |
|
w1 ak |
C |
|
TA = (E + wek>)A = A + wek>A = A + wak |
|
: : : |
: |
|||||
|
|
|
|
B |
|
C |
|
|
|
|
|
|
B |
|
C |
|
|
|
|
|
|
B |
wm ak |
C |
|
|
0 |
|
1 |
@ |
|
A |
|
||
|
|
|
|
|
|
|||
В итоге |
a1 + w1 ak |
|
|
|
|
|
|
|
TA = B |
|
|
C |
: |
|
|
|
|
a2 + w2 ak |
|
|
|
|
|
|||
B |
|
C |
|
|
|
|
|
|
B |
|
C |
|
|
|
|
|
|
B |
: : : |
C |
|
|
|
|
|
|
@ |
|
A |
|
|
|
|
|
|
am + wm ak
Определение 1.5. Поэлементное добавление одной из строк матрицы ко всем другим ее строкам (быть может, в промасштабированном виде) называется э л е м е н т а р н ы м преобразованием этой матрицы.
В качестве упражнения проверьте, что умножение слева (3 3)– матрицы A на матрицу вида
T = |
0 |
2 |
1 |
0 |
1 |
|
|
1 |
0 |
0 |
A |
|
@ 3 |
0 |
1 |
||
приведет к тому, что у исходной матрицы 1-я строка дважды прибавится ко 2-й и трижды вычтется из 3-й.
Аналогично умножение слева (3 3)–матрицы A на матрицу вида
01
T = @ |
1 |
0 |
2 |
0 |
1 |
3 |
|
0 |
0 |
1 A |
приведет к тому, что у исходной матрицы 3-я строка дважды прибавится к 1-й и трижды вычтется из 2-й.
15
1.5. Блочные матрицы и действия над ними
Определение 1.6. Блочной называется матрица, разделенная на ряд горизонтальных и/или вертикальных полос.
Блочную матрицу можно рассматривать как матрицу, чьи элементы сами являются матрицами и/или векторами. Например, следующая
матрица может быть представлена как блочная: |
|
|
|
|
|||||||||
|
M = 0 |
1 |
2 |
3 |
|
1 = |
M11 |
M12 |
|
|
|
|
|
|
4 |
|
|
|
|||||||||
|
|
|
; |
|
|
||||||||
|
2 3 4 |
5 |
|
|
|||||||||
|
B |
3 4 5 |
1 |
C |
M21 |
M22 |
|
|
|
||||
где |
@ |
|
|
|
|
A |
2 |
3 |
4 |
|
M22 = |
5 |
: |
|
|
|
|
|
|||||||||
M11 = ( 1 |
2 3 ); M12 = (4); |
|
M21 = 3 4 5 ; |
1 |
|||||||||
Большинство операций над блочными матрицами, размеры блоков которых с о г л а с о в а н ы, можно производить так, словно их блоки обычные числа. Необходимо только помнить, что матричные множители в общем случае нельзя менять местами.
Проиллюстрируем основные операции над блочными матрицами на простом примере. Пусть
A11 A12 B11 B12
A = ; B = :
Если размеры блоков согласованы, то верны следующие формулы:
умножение матрицы на число
|
A11 |
A12 |
; |
|
|
|
A = |
A22 |
|
||
|
A21 |
|
|
|
|
сложение матриц |
|
|
|
|
|
|
A11 + B11 |
A12 + B12 |
|
||
A + B = A21 + B21 |
A22 + B22 ; |
|
|||
перемножение матриц |
|
|
|
!; |
|
AB = |
A11B11 + A12B21 |
A11B12 |
+ A12B22 |
||
A21B11 + A22B21 |
A21B12 |
+ A22B22 |
|||
16
транспонирование
A11> |
A21> |
: |
A> = A12> |
A22> |
В частности, при умножении матриц
0 1
a1B
BC
Ba2B C
AB = B |
am...B |
C = (AB1jAB2j : : : jABn); |
B |
|
C |
@ |
|
A |
где ai строки матрицы A; а Bj столбцы матрицы B:
Блочные матрицы полезны при агрегированном описании и изучении больших систем, состоящих из большого числа взаимосвязанных подсистем, каждая из которых может иметь свои структурные и иные особенности.
Вопросы для самопроверки
1.Что такое матрицы? Из чего они состоят?
2.Что такое размерность матрицы?
3.Какие матрицы признаются равными?
4.Какие специальные классы матриц рассматриваются в линейной алгебре?
5.Перечислите основные арифметические операции над матрицами.
6.Каким законам подчиняются арифметические матричные операции?
7.В чем состоит специфика операций над диагональными матрицами?
8.В чем состоит специфика операций над треугольными матрицами?
9.В чем состоит специфика операций над матрицами, состоящими из одного столбца или из одной строки?
10.Чем интересны элементарные матрицы? Что такое элементарное преобразование матрицы?
11.Что такое блочная матрица? Сколько горизонтальных или вертикальных полос она может содержать?
17
12.Как складывать блочные матрицы?
13.Как умножать блочные матрицы на число?
14.Как умножать блочные матрицы одна на другую?
15.Как транспонировать блочные матрицы?
16.В табл. 1 приведены нормы расхода различных видов сырья на производство единицы той или иной продукции некоторого пищевого производства:
|
|
|
|
|
|
|
|
|
|
|
Таблица 1 |
||
|
Нормы расходы сырья на единицу готовой продукции |
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
В и д |
|
|
|
В и д п р о д у к ц и и |
|
|
|
|
||||
|
с ы р ь я |
1 |
2 |
|
3 |
4 |
5 |
|
6 |
|
7 |
|
|
1. |
Мука, |
кг |
7 |
3 |
|
12 |
1 |
3 |
|
9 |
|
3 |
|
2. Масло, |
кг |
3 |
9 |
|
5 |
8 |
1 |
|
7 |
|
5 |
|
|
3. Молоко, л |
11 |
1 |
|
8 |
4 |
5 |
|
2 |
|
9 |
|
||
4. Сахар, |
гр |
9 |
15 |
|
2 |
10 |
1 |
|
8 |
|
1 |
|
|
5. |
Цена, руб. |
7 |
5 |
|
8 |
2 |
9 |
|
4 |
|
3 |
|
|
В табл. 2 приведены данные о предполагаемых объемах производства перечисленной выше продукции:
|
|
|
|
|
|
Таблица 2 |
|
|
Предполагаемые объемы производства |
|
|
||||
|
|
|
|
|
|
|
|
|
|
В и д п р о д у к ц и и |
|
|
|
||
1 |
2 |
3 |
4 |
5 |
6 |
7 |
|
100 |
101 |
900 |
120 |
190 |
400 |
300 |
|
Как рассчитать суммарную потребность в сырье для производства перечисленных видов продукции в заданных объемах? Как рассчитать суммарную выручку от продажи этой продукции?
17. Обозначим матрицу норм расхода из предыдущей задачи как |
|
|||||||||||||
|
A = 0 3 |
9 |
5 |
8 |
1 |
7 |
5 1 |
; |
|
|||||
|
B |
|
7 |
3 |
12 |
1 |
3 |
9 |
3 |
C |
|
|
||
|
|
9 |
15 |
2 |
10 |
1 |
8 |
1 |
|
|
||||
|
B |
11 |
1 |
8 |
4 |
5 |
2 |
9 |
C |
|
|
|||
|
@ |
|
|
|
|
|
|
|
|
|
A |
|
|
|
а векторы цен и объемов производства соответственно через |
: |
|||||||||||||
b = |
7 5 8 2 9 4 3 |
|
и |
c = |
100 101 900 120 190 400 300 |
|||||||||
Как выразить суммарную потребность в сырье и суммарную выручку из предыдущей задачи на языке векторов и матриц?
2.ОПРЕДЕЛИТЕЛИ
2.1.Что такое определитель матрицы
Каждой квадратной матрице A ставится в соответствие некоторое число, называемое ее определителем (или детерминантом) и обозначаемое jAj или det(A).
Чтобы понять, как это делается, нам потребуется
Определение 2.1. Подстановкой порядка n называется произвольное взаимно-однозначное отображение множества Nn = f1; 2; : : : ; ng в себя.
Такого рода отображения задаются (2 n)-матрицами соответствия, первая строка которых содержит числа от 1 до n, а вторая образы этих чисел при рассматриваемом отображении.
Ниже приведены примеры подстановочных матриц 3-го порядка:
P = |
1 |
2 |
3 |
; P0 |
= |
1 |
2 |
3 |
; P00 |
= |
1 |
2 |
3 |
: |
2 3 1 |
3 2 1 |
1 3 2 |
||||||||||||
Обозначим через Pn совокупность всех подстановок порядка n.
r |
|
p |
подстановочной мат- |
Определение 2.2. Два столбца s |
и q |
||
рицы P образуют инверсию, если r < p, |
но s > q. |
||
Обозначим через J(P) число инверсий в подстановке P: В частности, для приведенных выше примеров
J(P) = 2; J(P0) = 3; J(P00) = 1:
Заметим, что перестановка строк или столбцов подстановочной матрицы не меняет числа ее инверсий.
Определение 2.3. Детерминантом (определителем) квадратной матрицы называется величина
|
|
X: : : |
|
|
J(P) |
jAn nj = |
1 2 |
n |
( 1) |
a1j1 a2j2 : : : anjn : |
|
P= |
j1 j2 |
: : : |
jn |
!2Pn |
|
Здесь суммируются всевозможные произведения ее элементов, взятых по одному из каждой строки и каждого столбца, причем каждое такое произведение берется со своим или с противоположным знаком в
19
зависимости от четности или нечетности числа инверсий в подстановке, образуемой номерами строк и столбцов, в которые эти элементы входят.
Проиллюстрируем приведенное выше общее правило вычисления определителя на матрицах небольшого порядка3.
Для матрицы 1-го порядка |
jaj = a: |
|
|
||
Для матрицы 2-го порядка |
|
|
|
|
|
|
|
|
|
|
|
a11 |
a12 |
= a11a22 |
a12a21: |
||
jA2 2j = a21 |
a22 |
||||
|
|
|
|
|
|
|
|
|
|
|
|
Для матрицы 3-го порядка |
|
|
|
|
|
|
a11 |
a12 a13 |
|
||
jA3 3j = |
|
|
|
|
|
|
a21 a22 a23 |
|
= |
||
|
|
|
|
|
|
a31 a32 a33
=a11a22a33 + a12a23a31 + a21a32a13 a31a22a13 a21a12a33 a11a32a23:
Перегруппировав слагаемые в последней формуле, приходим к соотношению
|
a11 |
a12 |
a13 |
= a11 |
a32 |
a33 |
a12 |
a31 |
a33 |
+ a13 |
a31 |
a32 |
; |
|||
a21 |
a22 |
a23 |
||||||||||||||
|
31 |
32 |
33 |
|
a22 |
a23 |
|
|
a21 |
a23 |
|
|
a21 |
a22 |
|
|
|
a |
a |
a |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где определитель матрицы большего размера вычисляется через определители ее подматриц размера, на единицу меньшего. В дальнейшем этот подход будет перенесен на определители n-го порядка.
2.2. Простейшие свойства определителей
Следующие утверждения являются достаточно очевидными следствиями общего определения либо последовательно вытекают друг из друга.
Свойство 2.1. jA>j = jAj: Свойство 2.2. j An nj = njAn nj.
3Отметьте на схеме матрицы те ее элементы, произведения которых входят в определитель со своим знаком, и те элементы, которые входят в него с обратным знаком. Какие фигуры они составляют?
20
