книги из ГПНТБ / Тихомиров В.И. Линейное программирование в организации и планировании путевого хозяйства конспект лекций для студентов специальности Стр-во ж. д., путь и путевое хоз-во учеб. пособие
.pdfЕдиничная матрица. Если в диагональной матрице все эле менты являются равными единице, то мы получим единичную матрицу, которую обычно обозначают буквой Е.
Ленточная матрица. Если в квадратной матрице порядка т отличными от нуля являются лишь элементы, расположен ные на главной диагонали и на примыкающих к ней с каждой из сторон (сверху и снизу) к параллельных линиях, то такую матрицу называют ленточной (2 к + 1) -членной матрицей по рядка т, так при к=1 матрица будет являться трехчленной ленточной и иметь вид:
« и |
« 1 2 |
|
|
|
o 2t |
« 2 з |
« 2 з |
|
|
|
« 3 2 |
а з З |
« з 4 |
|
|
|
« 4 3 |
и и |
а 4 й |
(m - i ) d ' m m
В литературе такие матрицы встречаются под названием «якобиевы матрицы».
Верхняя и нижняя треугольные матрицы. Квадратная матрица, в которой все элементы, расположенные ниже глав ной диагонали, равны нулю, называется верхней треугольной матрицей.
« и |
« 1 2 |
«13 |
■ |
• |
■ |
« 1 т |
|
« 2 2 |
«2э |
|
• |
• |
« 2 т |
|
|
«S3 |
|
• |
' |
(19) |
|
|
|
«3 т |
|||
|
0 |
|
|
|
|
« т т |
Аналогично нижней треугольной называется матрица сле |
||||||
дующего вида: |
|
|
|
|
|
|
« , 1 |
|
|
|
|
|
0 |
|
|
|
|
|
|
|
«21 |
«2 3 |
|
|
|
|
(19 0 |
«31 |
« з а |
«33 |
|
|
|
|
|
|
|
|
|||
«m l |
|
« т з |
■ |
■ |
‘ |
«mm |
20
§ 4. Алгебраические операции над матрицами
1. Сложение и вычитание матриц
Складывать и вычитать можно лишь матрицы, имеющие одинаковое число строк и столбцов. Суммой матриц А и В на зывается матрица С, элементы которой C,j (i=l, 2,..., m; j= 1, 2,..., n) связаны с элементами a,j и ft,,- матриц А я В равенством:
Си = аи + Ьц.
Иными словами, |
|
|
а и |
^1 2 |
• |
аг, |
а 22 |
• • • |
А + В — |
|
|
а 1п |
|
612 |
. |
. ■ 6Jn * |
а 2п |
+ |
Ь » 1 62з |
. |
. • ь г „ i _ |
|
|
|
|
|
|
|
|
■ * ^ т п |
|
|
Ь т \ |
|
• ^ т п |
||
|
а,. + Ь ч |
Йи + ь п |
|
. • • ^1л "Ь |
|
||||||
= |
#2| |
+ |
b ix |
|
+ Ь 22 |
|
• |
• • |
^зл ”Ь Ь гп |
= |
|
|
|
|
|
|
|
|
|
|
|
||
|
1^ т п \ + b m , а тЧ + Ь т 2 . * * |
“Ъ^тл |
|
||||||||
|
|
|
|
С и |
C i2 |
• |
• |
• |
С1Л |
|
|
|
|
|
|
сг\ |
• |
• |
• |
<■,„ |
• |
(20) |
|
|
|
|
|
Cm i |
^ m 2 |
• |
• |
• |
Ст л |
|
|
Если А — -В = С , |
то |
С „ |
- - - - |
А,- |
|
|
|
|
2. Умножение матрицы на скаляр
Произведение матрицы А на скаляр а представляет собой матрицу В, каждый элемент которой равен произведению соответствующего элемента матрицы А на скаляр а, т. е. bij = % а.
21
|
au |
alt |
. . |
• |
am |
|
alta |
aia a . |
• alna |
4 |
atl |
ati |
. . |
• |
Q-ln |
• a = |
a2i a |
a„a . • |
■a2n a |
II « |
|
|
|
|
|
|
|
Q'mi am2 • • * ^mn
b u |
b i t |
. |
• |
ь » |
b 22 |
|
• • |
• |
• amna |
b l n |
|
b 2 „ |
(21) |
= B . |
Ь m l b m 2 ■ * • ^ m n
3.Умножение прямоугольных матриц
Пусть имеем матрицы типов ш Хп и pXQ-
|
а и |
^12 |
• |
• |
®in |
|
A = |
^21 |
|
|
■ |
|
(22) |
|
|
|
|
|
||
|
aml |
|
• |
' |
^mn |
|
|
bn |
Ь\2 |
■ |
• |
b l q |
|
В = |
b?i |
b-)0 |
|
■ b 2q |
(23) |
|
|
|
|
|
|
||
|
V |
Ьр2 |
■ ■ ■ ЬРЯ |
|
Если число столбцов матрицы А равно числу строк матри цы В, т. е. п = р, то для этих матриц получается матрица С типа m</, называемая их произведением:
|
С 11 С 12 |
-1<7 |
|
|
С |
С2Ч |
(24 ) |
|
|
||
|
сmi |
Стд |
|
где |
Си — ап blj + а[2 b2j + ... + ain bnj\ |
|
|
(i = |
1, 2.........m); (j = |
1, 2, . . ., q). |
|
22
Из определения вытекает следующее правило умножения матриц.
Чтобы получить элемент, стоящий в г'-той строке и в /-том столбце произведения двух матриц, нужно элементы t-той строки первой матрицы умножить на соответствующие эле менты /-того столбца второй и полученные произведения сложить.
Произведение АВ имеет смысл тогда и только тогда, ког да матрица А содержит в строках столько элементов, сколько элементов 'имеется в столбцах матрицы В.
Матрицы, которые можно перемножать, называются соот ветственными.
П р и м е р |
1 |
|
|
|
|
|
|
|
|
|
|
2 |
1 |
6 |
|
5 |
2 |
|
|
|
|
2 |
4 |
|
1 |
|
|
||
|
|
|
в = 1 |
7 |
|
||||
|
|
3 |
2 |
|
5 |
|
|||
|
|
|
|
4 |
2 |
|
|||
|
|
7 |
3 |
|
1 |
|
|
||
|
|
|
|
|
|
|
|||
|
2. 5 |
|
+ |
Ь 1 2+-62. +4 |
1-7 + |
6-2 |
35 23 |
||
_ |
_ 2 • 5 + 4 • 1 + 1-4 |
2 - 2 + 4 - 7 + 1 - 2 |
18 34 |
||||||
= |
~ 3-5 + 2-1 + 5 - -J |
3-2 + 2-7 + 5-2 |
37 30 |
||||||
|
7-5 + 3- 1 + 1-4 |
7.-2+ 3-7 + 1 • 2 |
42 37 |
||||||
П р и м е р |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 — 1 |
|
|
|
|
3 |
2 |
8 |
1 |
|
1 - 3 |
|
|
|
|
1 - 4 |
0 |
3 |
|
0 |
1 |
|
|
|
|
|
|
|
|
|
3 |
1 |
|
|
3-2 + 2 - 1 + 8 - 0 + 1-3 |
—3 -1 -2 -3 + 8-1 + 1-1 |
|||||||
АВ = С = |
|
|
|
|
|
—1-1 + 4-3 + 0-1 + 3-1 |
|||
|
1-2 — 4 - 1 + 0 0 + 3-3 |
11 0
7 14
Произведение двух матриц, как правило, не обладает пе реместительным свойством, т. е. АВф ВА .
23
П р и м е р 3 |
|
|
|
|
|
1 |
2 |
5 |
6 |
|
А = |
4 |
5 = |
8 |
|
3 |
7 |
||
Тогда |
19 |
22 I |
23 |
34 |
АВ = |
50 ! |
ВА = |
46 |
|
|
43 |
31 |
т. е. здесь АВф ВА .
В тех случаях, когда АВ = ВА, матрицы А и В называются перестановочными (коммутативными). Так, например, не трудно убедиться, что единичная матрица Е перестановочна с любой квадратной матрицей А того же порядка.
5. Транспонирование матриц
Положим, 'имеем матрицу
а , , |
Й ,2 |
• |
• |
• ^1Л |
|
a2i |
а 22 • • |
* а 2л |
(25) |
||
|
|
|
|
|
|
а т г |
Q ' m t |
■ |
• • ^ т п |
|
Если в этой матрице переменить местами строки и столбцы, то получим матрицу А \ которая называется матрицей, транс понированной по отношению к матрице А.
|
|
Ди |
а21 |
• |
• «ml |
|
А1 |
= |
alt |
£Zaa |
• |
* ат2 |
(26) |
|
|
|
|
|||
|
|
aUi |
а2п |
. * ^тп |
|
6. Обратная матрица
Обратной матрицей по отношению к данной называется матрица, которая, будучи умноженной как справа, так и слева на данную матрицу, дает единичную матрицу.
Для матрицы А обозначим обратную ей матрицу через А"1. Тогда в соответствии с указанным выше определением имеем:
АА~1=:А-' А = Е,
где Е — единичная матрица.
24
Всякая неособенная (невырожденная) матрица А (при Б[А]фО) имеет одну и только одну обратную матрицу Л-1.
Пусть, например, дана система
|
@nx\-\~@i2 x2 Jr @1зхз— Н\ |
|
|
||||||
|
®2l -^1 ~|~ ^22%2~t- ^23Х3== Я2 |
|
(27) |
||||||
|
@81•'-1 “Ь" |
|
-J- @33Л'з = |
Яз |
|
|
|||
или в матричной форме |
|
|
|
|
|
|
|
||
|
|
@ и |
а 1 2 |
@ 1 3 |
|
|
Н |
х |
|
Л-Х = Я; |
@ и |
@ 2 1 |
@ 2 3 |
X Х 2 |
— Я2 |
( 28) |
|||
|
|
@ 3 1 |
@ 3 1 |
@ 3 3 |
|
* 3 |
Я3 |
|
|
причем Det А ф 0. Здесь х — неизвестный вектор. |
|
||||||||
Решение такой системы, |
пользуясь |
правилом |
Крамера, |
||||||
будет |
|
|
|
|
|
|
|
|
|
я . |
@ - ч |
@ 1 3 |
|
|
|
@ 1 1 |
Я , |
@ 1 3 |
|
я 2 |
@ 2 2 |
@ 2 з |
|
|
|
@ 2 1 |
Я 2 |
@ 2 з |
|
Я 3 |
@ 3 1 |
@ 3 3 |
' |
|
2 |
@ 1 1 |
Я з |
@ 3 3 |
|
|
Det А |
|
|
Det А |
|
|
@ 1 1 |
@ 1 1 |
@ 2 1 |
@ 2 2 |
@ 3 1 |
@ 1 3 |
я,
я,
4 |
со |
; |
Det А
Здесь в числителе стоят определители, раскрывая которые по элементам Ht, получим:
@ 2 2 |
@ 2 8 |
- я |
@ 1 2 |
@ 1 1 |
|
@ 1 2 |
@ 1 1 |
|||
Я з |
|
|
|
2 |
|
|
+ Я , |
|
|
|
@ 3 1 |
@ 1 3 |
|
@ 3 2 |
@ 3 3 |
|
@ 2 2 |
@ 2 3 |
|||
|
|
|
|
|
Det А |
|
|
|
|
|
@12 |
@2з |
|
|
@ 1 2 |
@ 1 3 |
|
|
@11 |
@ 1 з |
|
@12 |
@3з |
|
И |
@31 |
@ 3 3 |
л _ |
и |
@21 |
@13 |
|
Det |
А |
|
|
Det А |
- т |
п , - |
Det А |
|||
|
— п 2 |
25
|
|
«21 |
'2 3 |
+ яа |
«п |
«,э |
|
|
«, |
a, |
||
- ч = |
Lf |
' « 3 1 ___ « S 3 |
«*, |
a , |
- |
|
|
a 23 |
||||
- Я , |
- d 7 F a - |
Det A |
|
Det A |
||||||||
|
|
|
|
|
|
|
|
|
Я 3Д 21 |
|
||
|
|
я21 |
|
«22| |
|
«и |
^1а 1 |
|
|
«... |
«12 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Я, |
а»1 |
|
«за |
— Н .- |
«31 |
fl“ - + |
|
я .- |
&21 |
&2а |
|
|
|
Пе/ |
A |
|
|
|
|
|
|
Det A |
||
■Коэффициенты при Hi образуют матрицу |
|
|
|
|
||||||||
|
|
(Х4 4 |
« 2 3 |
« 1 2 |
|
« 1 3 |
|
|
«13 |
« . , |
|
|
|
|
а „ , |
« 3 2 |
« 3 2 |
|
«з, |
|
|
|
«зз |
|
|
|
|
D e t А |
D e t А |
|
|
D e t A |
|
|||||
|
|
« „ |
С1 |
7 3 |
« и |
|
« 1 3 |
|
|
« а |
^13 |
|
|
|
« 3 1 |
«33 |
^3 1 |
|
«33 |
|
|
&2 \ |
«зз |
|
|
|
|
D e t А |
|
D e t А |
|
|
D e t A |
|
||||
|
|
«31 |
« 2 2 |
« 1 1 |
|
«13 |
|
|
« п |
«13 |
|
|
|
|
« Я 1 |
« 3 , |
«3. « 3 2 |
|
|
&ч\ |
Clод |
|
|||
|
|
D e t А |
|
D e t А |
|
|
D e t A |
|
||||
т. е. |
В = А~\ |
что |
легко |
проверить |
|
умножением Л"1А, |
в результате чего получим Е.
Сопоставляя матрицы А и В, нетрудно заметить, что каж дый элемент (а^-1) обратной матрицы В представляет собой алгебраическое дополнение элемента ам матрицы А, делен ное на детерминант этой матрицы, т. е.
А ы
Det А
где Akl—алгебраическое дополнение (адъюнкта) элемента матрицы А или определитель п—1 порядка, образованный из данного определителя зачеркиванием к-той строки, и /-того
•столбца, взятого со знаком, определяемым формулой (—1)*+/. К обратной матрице часто прибегают при решении системы
линейных уравнений.
П р и м е р . Пусть дана система уравнений:
X1+ Х%—3
Ах\ + х%—5х3 ——9
Х{ + 2хг+ 4хз= 17.
26
Матрица данной системы будет |
|
|||
|
1 |
1 |
0 |
|
А = |
4 |
1 —5 |
|
|
|
1 |
2 |
4 |
|
Детерминант матрицы А |
|
|
|
|
Det А =4—5—16+10 = —7. |
||||
Алгебраические дополнения |
элементов |
матрицы |
||
1 _с |
= 4 + 10 = 14 |
|||
Ли = (- 1 )г+1 |
4 |
|||
2 |
|
|
|
|
4 —5 |
|
|
|
|
Д а = ( - 1 Г 2 |
= — (16 + 5) = — |
|||
1 |
4 |
|
|
|
4 |
1 |
|
|
|
Д 3 = ( - 1 )1+3 |
= 8 — 1 = 7 |
|
||
I |
2 |
|
|
|
1 |
0 |
— 4 |
|
|
Л . = ( - 1)2+1 |
= |
|
||
2 |
4 |
|
|
|
1 |
0 |
4 |
|
|
Л22 = ( - 1 ) ^ |
= |
|
|
|
I |
4 |
|
|
|
1 |
1 |
|
|
|
Л2з= ( - 1 ) ^ |
------(2—1) = |
— 1 |
||
1 |
2 |
|
|
|
1 |
0 |
|
|
|
л„ = ( - 1 Г 1 |
= — 5 |
|
||
1 - 5 |
|
|
|
|
I |
С■ |
= |
5 |
|
л зг = ( ~ \ Г 2 |
- 5 |
|
||
4 |
|
|
|
|
1 |
1 |
|
|
|
Л зЗ = (—1)3+3 |
= |
1 — 4 = — 3. |
||
4 |
1 |
|
|
|
27
Обратная матрица |
|
|
|
|
|
|
|
|
14 |
4 |
5 |
|
—14 |
4 '} |
5 |
|
7 |
У |
7 |
|
|||
|
|
|
i |
|
|||
1 |
21 |
4 |
5 |
_ |
1 |
—5 |
|
II |
7 |
|
~ |
7 а [2» |
|||
|
7 |
7 |
|||||
|
|
|
|||||
|
7 |
1 |
3 |
|
— 7 |
1 |
3 |
|
~~ 7 |
7 |
7 |
|
|||
|
|
|
|
|
Решение системы уравнений:
х = A~l Н
или
х % _ 1 |
- 1 4 |
4 |
5 |
3 |
21 - 4 - 5 |
X - 9 |
|||
“ 7 |
—7 |
1 |
3 |
17 |
X, |
||||
(— 42—36 + 85) |
1 |
|||
|
(63 + 36—85) |
= 2 |
||
|
(—21— 9+51) |
3 |
||
т. е. *1 = 1; |
* 2 = 2 |
*з= 3. |
§5. Жордановы исключения
1.Обыкновенные Жордановы исключения
Для решения систем линейных уравнений и задач линей ного программирования используется метод последователь ных исключений, или как его называют, метод Жордана, пред ставляющий собой совокупность удобных вычислительных алгоритмов, построенных на последовательном применении, эквивалентных преобразований системы уравнений или сис темы векторов.
Пусть рассматривается система:
у1= ап х1+ аи*а + . . . + alnxn; (t — 1, . . . . т). (29>
Эта система состоит из т линейных форм с п независимым» переменными х,, *2, *„• Запишем эту систему в виде следующей таблицы (табл. 2):
28
|
|
|
Таблица 2 |
|
||
|
- у , . . • >xs, . |
* »х п |
|
|||
У\ = |
а и • ■ • а\s • . . Л| д |
|
||||
Уг в |
ап . |
• |
ars • |
. |
агп |
|
Ут — |
&т\ • |
■• |
am s•• • |
|
атп |
|
Шагом обыкновенного Жорданова |
|
исключения, произве |
||||
денного над табл. 2 |
с |
разрешающим элементом |
аг 5ф О |
|||
с r-той разрешающей строкой и s-тым |
разрешающим |
столб |
цом, называется операция перемены ролей между зависимой
переменной уг и |
независимой xs, т. е. операция решения |
уравнения |
|
У г = |
^21 Х2 -f- • . . -f- CLrs Л', ”1“ *• - "Ь ^гп Xп |
относительно xs подстановка его во все остальные уравнения
•системы (29) и запись полученной системы в виде новой таблицы, аналогичной табл. 2. Легко проверить, что получен ная таким образом новая таблица (табл. 3) будет иметь вид:
|
|
|
|
|
|
|
Таблица 3 |
|
|
|
Х1» *2, .. . . Уг. ■ • ■ хп |
|
|
||||||
зч - |
би |
613. |
• |
а\S * • |
• |
Ь\п |
|
|
|
у„ = |
ЬпЛ |
|
|
. |
a,s . ■ |
■ |
bin |
|
|
Xs = -- &г\ —1 |
• . |
. |
1 . . ■ — агп |
|
|
||||
Ут = |
|
• |
* |
ams • • • |
&тп |
|
|
||
где Ьц = au ars — a!sarJ {i ф г, |
j Ф s), |
причем |
все |
элементы |
|||||
таблицы следует разделить на |
|
ars. Следовательно, |
один шаг |
||||||
Жорданова исключения |
с |
разрешающим |
элементом ars |
||||||
переводит табл. 2 в табл. |
3 по схеме, состоящей из следующих |
||||||||
пяти правил: |
|
|
|
|
|
|
|
|
|
1.Разрешающий элемент заменяется единицей.
2.Остальные элементы разрешающего столбца (s-того) остаются без изменения.
3.Остальные элементы разрешающей строки (r-той меня ют лишь свои знаки.
4.Элементы, не принадлежащие к разрешающей строке или столбцу, вычисляются по формуле
Ь-^= аи Qrs uls cirj.
29