Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Тихомиров В.И. Линейное программирование в организации и планировании путевого хозяйства конспект лекций для студентов специальности Стр-во ж. д., путь и путевое хоз-во учеб. пособие

.pdf
Скачиваний:
11
Добавлен:
23.10.2023
Размер:
3.99 Mб
Скачать

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

Ленточная матрица. Если в квадратной матрице порядка т отличными от нуля являются лишь элементы, расположен­ ные на главной диагонали и на примыкающих к ней с каждой из сторон (сверху и снизу) к параллельных линиях, то такую матрицу называют ленточной (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

Соседние файлы в папке книги из ГПНТБ