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

ИсследОпераций

.pdf
Скачиваний:
28
Добавлен:
08.05.2015
Размер:
4.21 Mб
Скачать

 

 

 

 

 

Часть I. ЭЛЕМЕНТЫ ЛИНЕЙНОЙ АЛГЕБРЫ

 

 

 

 

 

 

 

 

 

 

 

 

 

§ 1. Основные обозначения

 

 

 

 

 

 

 

 

En

n -мерное евклидово пространство;

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x x2 n -мерный вектор-столбец или точка пространства E

n

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x x , x , ..., x

– вектор-строка;

 

 

 

 

 

 

 

 

 

 

 

 

x;Q

1

2

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

– множество всех точек x En , обладающих свойством Q ;

 

 

 

 

 

 

O x

– окрестность точки x , т.е. любое открытое ограниченное множество то-

 

 

 

 

чек x E

n

, содержащее точку x ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D S означает, что множество D является подмножеством множества S ;

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x, y xi yi

– скалярное произведение векторов x

и y ;

 

 

 

 

 

 

 

 

x, y

i 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

отрезок

 

с концами x

и

y , т.е. множество {z : z x y ,

где

 

 

 

0, 0, 1};

 

 

 

 

 

 

 

 

 

 

 

 

AT

– матрица, транспонированная к матрице A ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i 1, n означает, что i пробегает все целые значения от 1 до n ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x y означает, что xi yi при i 1, n ;

 

 

 

 

 

 

 

 

 

 

f x f x1, x2,..., xn – функция n переменных;

 

 

 

 

 

 

 

 

 

min f x – запись задачи минимизации функции f x на множестве D ;

 

 

 

 

x D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

x

– градиент функции

f x

в точке x , т.е. n -мерный вектор с компонен-

 

 

 

 

 

f

x k

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

тами

1, n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H f

– гессиан функции

f x

в точке x , т.е. квадратная матрица n -го по-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 f

x

 

 

 

 

 

 

рядка,

элементы

 

которой

определены

равенством

h

 

 

i

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i, j

 

xi

x j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i1, n ; j 1, n ;

– означает «равносильно» или «эквивалентно»;

– символ, означающий конец доказательства.

3

§ 2. Основные понятия и определения

Определение. Матрицей размера m n называется прямоугольная таблица

чисел, состоящая из m строк и n столбцов.

 

 

a

a

...

a

 

 

 

 

 

11

12

 

1n

 

 

 

 

 

Am n a21

a22

...

a2n

,

A aij ,

i 1,..., m,

j 1,..., n .

 

...

...

...

 

 

 

 

 

...

 

 

 

 

 

am1

am2

...

amn

 

 

 

 

Если число строк в матрице равно числу столбцов, то матрица An n называ-

ется квадратной матрицей n -го порядка.

Определение. Квадратная матрица E , у которой на главной диагонали стоят единицы, а на остальных местах нули, называется единичной матрицей

 

1

0

...

0

 

 

0

1

...

0

 

E

.

 

 

...

...

...

 

...

 

 

0

0

...

1

 

Определение. Суммой матриц A и B одинакового размера называется матрица C такого же размера, каждый элемент которой равен сумме соответствующих элементов слагаемых матриц.

Cm n Am n Bm n , cij aij bij .

Определение. Произведением матрицы A на число называется матрица B того же размера, полученная умножением каждого элементы матрицы A на число .

Bm n Am n , bij aij .

Определение. Произведением матрицы Am e на матрицу Be n называется матрица Cm n , каждый элемент которой cij равен сумме попарных произведений элементов i -ой строки первой матрицы и элементов j -го столбца второй матрицы.

cij ai1 b1 j ai2 b2 j ... aie bej . Cm n Am e Be n .

Умножение матриц определено лишь для таких двух матриц, что число столбцов первого сомножителя равно числу строк второго. Заметим, что AB BA даже когда оба произведения существуют.

Определение. Матрица AT , полученная из матрицы A заменой ее строк столбцами с теми же номерами, называется транспонированной к данной.

4

a11

a

Am n 21

...

am1

2

Пример. A 4

Р е ш е н и е.

a

 

...

a

 

 

 

12

 

 

1n

 

 

a22

 

...

a2n

,

...

 

...

...

 

 

 

 

 

 

am2

 

...

amn

 

1

1

 

 

1

 

 

 

 

3

2

 

,

B

3

 

 

 

 

1

 

 

 

 

 

 

a

a

...

a

 

 

11

21

 

m1

 

AT

a12

a22

...

am2

.

n m

 

 

... ...

 

 

... ...

 

 

a1n

a2n

...

amn

 

2

1 . Найти A B .

0

 

 

 

 

 

 

 

2

1 1

 

1

2

 

 

6

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A B

 

3 2

 

 

3 1

 

 

 

.

 

 

 

 

 

 

 

 

 

4

 

 

1

0

 

 

7

11

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c11 2 1 1 3 1 1 6,

 

 

 

 

 

 

 

 

с12 2 2 1 1 1 0 3,

 

 

 

 

 

 

 

 

 

 

 

с21 4 1 3 3 2 1 11,

 

 

 

 

 

 

 

 

с22 4 2 3 1 2

0 11,

 

 

 

 

 

 

 

 

 

 

 

 

A2 3 B3 2 C2 2 .

 

 

 

 

 

Пример. A

1 2

1

 

3

 

1

. Найти: 1)

A B ;

2) B A ;

3) AT B .

 

 

 

 

,

B

 

 

 

 

 

 

 

3

 

2 1

 

 

1

 

0

 

 

 

 

 

 

 

 

Р е ш е н и е.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1) Произведение

A B найти невозможно, так как число столбцов матри-

цы A не равно числу строк матрицы B .

 

 

 

 

 

 

 

 

 

 

 

3

1 1 2

1

 

0 8 4

 

 

 

 

 

2) B

A

0

 

 

 

 

 

 

 

1 2 1

.

 

 

 

 

 

 

 

1

3

2 1

 

 

 

 

 

 

 

 

1

3

 

 

 

1

 

 

3

3

1

 

6 1

 

T

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

3) A

 

2

2

 

, A

B

 

2

 

2

 

 

1

0

 

 

4 2

.

 

 

 

1

1

 

 

 

 

 

1

 

 

1

 

 

 

 

2 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Определение. Пусть A – квадратная матрица 2-го или 3-го порядков. Определителем квадратной матрицы называется число, полученное по следующему правилу:

a

a

 

 

 

 

a

a

 

a11

a22

a12 a21 .

 

 

 

 

A 11

12

 

,

A

 

11

12

 

a21

a22

 

 

 

 

a21

a22

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

a11

 

a12

 

a13

 

 

 

 

 

 

a11

a12

a13

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A a

 

a

 

a

 

,

 

A

 

 

a

a

a

 

 

 

 

 

 

 

 

 

 

21

 

22

23

 

 

 

 

 

 

21

22

23

 

 

 

 

 

 

a32

 

a33

 

 

 

 

 

 

a31

a32

a33

 

 

 

 

a31

 

 

 

 

 

 

 

 

 

a11a22a33 a12a23a31

a13a21a32

a13a22a31 a12a21a33 a11a23a32 .

Обозначается определитель

 

A

 

A .

 

 

 

 

 

 

 

 

 

 

2

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример.

 

3

1

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Р е ш е н и е.

 

 

 

 

1

2

 

 

 

 

 

 

2

 

 

 

 

 

 

A

 

 

3

1

0

2 0 12 2 3 0 5 .

 

 

 

 

 

 

 

 

 

 

 

1

2

1

 

 

 

 

 

Определение. Если определитель квадратной матрицы A не равен нулю,

 

A

 

0 , то A называется невырожденной матрицей.

 

 

 

 

 

 

 

Определение. Пусть A

квадратная матрицы n -го порядка. Квадратная

матрица A 1 того же порядка называется обратной к A , если

A A 1 A 1 A E .

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

Определение. Линейной системой m уравнений с n неизвестными называется запись вида

a11x1 a12 x2

... a1n xn b1,

 

a22 x2

... a2n xn b2 ,

a21x1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a x

a

x

 

... a

x

b ,

m1 1

m2 2

 

mn n

m

где x1, x2,..., xn – неизвестные,

aij

– коэффициенты при неизвестных, b1, ..., bm

свободные члены.

 

 

 

 

 

 

Определение. Числа c1, c2,..., cn

называются решением системы, если, при

подстановке вместо неизвестных,

они обращают каждое уравнение системы в

верное равенство. Если система имеет хотя бы одно решение, она называется совместной, в противном случае – несовместной.

Определение. Матрица, составленная из коэффициентов при неизвестных, называется «матрица системы»

6

a

a

...

a

 

11

12

 

1n

 

A a21

a22

...

a2n .

 

...

...

...

 

...

 

am1

am2

...

amn

 

x

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

Матрица-столбец X x2

называется «столбец неизвестных», матрица-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xn

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

столбец B b2

называется «столбец свободных членов». Матрица, получаемая

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

bm

 

 

 

 

 

 

 

 

 

присоединением к матрице A столбца свободных членов, называется «расши-

ренная матрица системы»

 

 

 

 

 

 

 

 

 

 

 

a

a

 

...

a

 

b

 

 

 

 

11

12

 

1n

 

1

 

 

A

a21

a22

...

a2n

b2 .

 

 

 

 

... ... ...

 

...

 

 

 

...

 

 

 

 

am1

am2

...

amn

bm

Определение. Пусть дана линейная система m уравнений с n неизвестны-

ми

 

 

 

 

 

 

 

 

 

 

 

a11x1 a12 x2

... a1n xn b1,

 

 

 

 

 

 

... a2n xn b2 ,

 

a21x1 a22 x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

x

a

x

 

... a

x

b .

 

m1 1

 

m2 2

 

mn

n

 

m

По правилу умножения матриц

 

a

a

...

 

11

12

 

A X

a21

a22 ...

 

... ...

 

 

...

 

am1

am2 ...

a

 

 

x

 

 

a x

1n

 

 

1

 

 

11 1

a2n

 

x2

 

 

a21x1

...

 

 

 

 

 

 

 

 

...

 

 

...

amn

 

xn

 

am1x1

a

x

...

a

x

 

12

2

 

1n

n

 

a22 x2

...

a2n xn

...

...

...

 

am2 x2

 

 

 

 

...

amn xn

– матрица размера m 1.

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

A X B .

7

Определение. Евклидово пространство En – это множество n -мерных векторов x x1, x2, ..., xn , в котором определены операции:

сложения

xy x1, x2, ..., xn y1, y2, ..., yn

x1 y1, x2 y2, ..., xn yn ,

умножения на число

x x1, x2, ..., xn ,

скалярное произведение векторов

x, y x1y1 x2 y2 ... xn yn .

Длина вектора вычисляется по формуле

x x12 x22 ... xn2 .

§ 3. Матрицы

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

Определение. Пусть дана произвольная матрица A и aij – некоторый, отличный от нуля ее элемент. Жордановым исключением с ведущим элементом aij

называется переход от матрицы A к матрице B того же размера, осуществляемый по следующим правилам:

1) все элементы ведущей (содержащей ведущий элемент) строки делятся на

aij ;

2)все остальные элементы ведущего столбца заменяются нулями;

3)остальные элементы матрицы B bkl (k i , l j ) вычисляются по фор-

муле

bkl aijakl ail akj . aij

З а м е ч а н и е 1. Расчет по последней формуле удобно производить, пользуясь «правилом прямоугольника»: для вычисления элемента bkl матрицы B ,

стоящего на месте элемента akl исходной матрицы A , мысленно выделяем прямоугольник с вершинами на элементе akl и ведущем элементе aij со сторонами, идущими по строкам и столбцам исходной матрицы A ,

8

 

 

l

 

j

 

 

 

 

столбец

столбец

 

 

 

 

 

 

 

 

 

 

 

 

akl

 

akj

 

k -я строка

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ail

 

aij

 

 

i-я строка

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ведущий элемент

Для нахождения bkl нужно из произведения элементов, стоящих в вершинах вы-

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

Пример 1. Дана матрица

 

2

3

7

5

 

 

 

 

 

 

 

A

1

0

2

4

.

 

3

1

1

3

 

 

 

Выполнить жорданово исключение с ведущим элементом a23 2 .

Р е ш е н и е. В результате жорданова исключения получим матрицу B того же размера, т.е. содержащую 3 строки и 4 столбца. Вычисление элементов матрицы B удобно производить в том порядке, который указан в определении жорданова исключения.

1) Сначала заполняем вторую строку матрицы B (она соответствует ведущей строке матрицы A ). Для этого соответствующие элементы ведущей строки

делим на ведущий элемент. Получим

 

 

 

 

 

 

 

 

 

 

 

 

0 1 2

 

B 1 / 2

.

 

 

 

 

 

 

 

 

 

 

 

 

2) Затем заполняем нулями третий (соответствующий ведущему) столбец матрицы B . Получим

9

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

0 1 2

 

 

 

 

 

B 1 / 2

.

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

3) Наконец, вычисляем по правилу прямоугольника значения остальных

элементов матрицы B . Например,

 

 

 

 

 

 

 

 

 

 

b

 

2 2 1 7

 

 

3

;

b

2 3 0 7

3;

 

 

 

 

11

2

 

 

2

 

12

 

 

 

2

 

 

 

 

 

 

 

 

 

 

b

2 5 4 7

9

и т.д.

 

 

 

 

 

 

 

 

 

14

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В результате получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 / 2

3 0

9

 

 

И

 

1 / 2

0 1 2

 

 

 

 

.

 

 

 

 

 

7 / 2

1

0

1

 

 

 

 

 

 

 

При выполнении жорданова исключения ведущий элемент всегда будем выделять рамкой, а исходную матрицу A с полученной матрицей B часто будем соединять знаком . С учетом этих обозначений выполненное преобразование запишется в виде:

2

3 7

5

3 / 2

3 0

9

 

1 0

2

4

 

 

1 / 2

0

1

2

 

 

 

 

.

 

3

1

1

3

 

 

7 / 2

1

0

1

 

 

 

 

 

Жордановы исключения широко используются при решении многих математических задач (особенно в курсах линейной алгебры и математического программирования). Поэтому необходимо научиться выполнять эти преобразования быстро и безошибочно.

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

З а м е ч а н и е. Если в ведущей строке (ведущем столбце) содержится нулевой элемент, то при выполнении жорданова исключения все элементы соответствующего столбца (строки) не изменяются. Так в примере 1 второй столбец матрицы A после выполнения жорданова исключения не изменился.

Пример 2. Дана матрица

 

3

0

7

1

 

 

 

 

 

 

 

A

2

1

0

0

.

 

0

2

4

5

 

 

 

 

 

10

 

 

 

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

Р е ш е н и е. Если нет никаких дополнительных условий, при выборе ведущего элемента обычно руководствуются следующими соображениями:

1)ввиду того, что при выполнении жорданова исключения приходится делить на ведущий элемент, удобно (если это возможно) выбрать ведущий элемент, равный единице;

2)учитывая замечание, ведущий элемент следует выбрать так, чтобы в его строке и в его столбце было как можно больше нулей.

Учитывая оба указанных соображения, для данной матрицы A в качестве ведущего следует взять элемент a22 1.

Выполним жорданово исключение с выбранным ведущим элементом a22 .

При преобразовании матрицы A в матрицу B :

а) вторая (ведущая) строка не изменится, так как ведущий элемент равен единице;

б) второй (ведущий) столбец заполнится нулями:

в) не изменятся также первая строка, третий и четвертый столбцы матрицы A , ввиду замечания.

Поэтому имеем

 

 

 

 

 

 

3

0

7

1

 

 

2

1

0

0

 

B

 

 

 

0

4

5

 

 

 

и по правилу прямоугольника, за счет удачного выбора ведущего элемента, придется вычислить только элемент

 

 

 

b

1 0 2 2

4 .

 

 

 

 

 

 

 

 

 

 

 

13

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом,

 

 

 

 

 

 

 

 

 

 

 

3

0

7

1

3

0

7

1

 

2

1

0

0

 

 

2

1

0

0

 

 

 

 

.

 

0

2

4

5

 

 

4

0

4

5

 

 

 

 

 

Рассмотрим теперь применение жордановых исключений для нахождения ранга и обращения матриц.

Пусть дана произвольная матрица A . Требуется найти ее ранг r A .

Алгоритм нахождения ранга матрицы.

1) Если все элементы данной матрицы A равны нулю, то задача решена, r A 0 .

11

2)Если среди элементов матрицы A имеются ненулевые, то любой из таких элементов выбираем в качестве ведущего.

3)Выполняем жорданово исключение с выбранным ведущим элементом.

4)У полученной матрицы помечаем (например, символом V ) строку, соответствующую ведущей строке выполненного исключения, а также все строки, которые в предшествующей матрице имели ту же пометку (в начале процесса таких строк нет, но они появятся по ходу вычислений).

5)В полученной матрице вычеркиваем строки, состоящие из одних нулей (если такие строки существуют).

6)Рассматриваем наличие пометок у строк полученной матрицы. Возможно два случая:

а) помечены не все строки. Тогда любой ненулевой элемент, стоящий в непомеченной строке, выберем в качестве ведущего и переходим к шагу 3 данного алгоритма;

б) помечены все строки. Тогда вычислительный процесс закончен. Ранг исходной матрицы A равен числу строк матрицы, полученной в результате вычислений.

Пример 3. Найти ранг матрицы

 

3

5

4

2

 

 

1

2

1

3

 

A

.

 

1

3

2

0

 

 

 

 

2

1

3

3

 

Р е ш е н и е. Применяя описанный алгоритм, получаем

 

 

 

 

№ 1

 

 

 

 

№ 2

 

 

 

 

№ 3

 

 

3

5 4 2

 

 

0 4

2

2

 

0

6 0

8

 

 

1

2

1 3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

5

1 3

 

V

0

5 1

3

 

 

1

3 2 0

 

 

V

1

3

2 0

 

V

1

7 0

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

1 3 3

 

 

 

0

5

1 3

 

 

0

0 0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

№ 4

 

 

 

 

 

 

 

 

 

 

 

V

0

3 / 4 0

1

 

 

 

 

 

 

 

 

 

 

V

 

0

11 / 4 1

0

 

 

 

 

 

 

 

 

 

 

 

r A 3 .

 

 

 

 

 

V

 

1

5 / 2

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Прокомментируем предложенное решение, занумеровав матрицы, полученные в процессе вычислений.

В качестве ведущего элемента для первого жорданова исключения может быть выбран любой ненулевой элемент исходной матрицы №1. Для простоты вы-

12