элементы линейной алгебры - Кыргызско

advertisement
КЫРГЫЗСКО- РОССИЙСКИЙ СЛАВЯНСКИЙ УНИВЕРСИТЕТ
ЕСТЕСТВЕННО-ТЕХНИЧЕСКИЙ ФАКУЛЬТЕТ
КАФЕДРА ВЫСШЕЙ МАТЕМАТИКИ
И. А. УСЕНОВ, Р.К. УСЕНОВА
ЭЛЕМЕНТЫ ЛИНЕЙНОЙ АЛГЕБРЫ
ЛЕКЦИИ ПО КУРСУ МАТЕМАТИКИ
Бишкек-2011
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
УДК 512.643(575.2) (075)
Рецензенты:
канд. пед. наук, доц. Ж.Р. Жаналиева,
канд. ф.-м. наук, доц. А.К. Курманбаева
Рекомендовано к изданию УС ЕТФ КРСУ
Усенов И.А., Усенова Р.К.
У 74 ЭЛЕМЕНТЫ ЛИНЕЙНОЙ АЛГЕБРЫ. Лекции по курсу математика:
Бишкек: КРСУ, 2011. 86 с.
Содержатся теоретические сведения по линейной алгебре. Основные внимание
уделяется матрицам, определителям и системам линейных уравнений, а также
практическому освоению студентами материала. Для достижения этой цели приводиться
большое число упражнений.
Пособие предназначено для студентов 1-2 курсов экономических специальностей
очной, заочной и дистантной форм обучения.
© КРСУ, 2011
2
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
ОГЛАВЛЕНИЕ
Предисловие .................................................................................................. 3
Глава 1. МАТРИЦЫ И ОПРЕДЕЛИТЕЛИ
Матрицы.........................................................................................................
4
Определители................................................................................................. 10
Вектор……………………………………………………………………….
27
Глава 2. СИСТЕМЫ ЛИНЕЙНЫХ АЛГЕБРАИЧЕСКИХ
УРАВНЕНИЙ
Системы линейных уравнений..................................................................... 29
Использование систем линейных уравнений при решении
экономических задач……………………………………………………… 46
Глава 3. ЭЛЕМЕНТЫ МАТРИЧНОГО АНАЛИЗА
Векторные (линейные) пространства……………………………………... 52
Размерность и базис векторного пространства…………………………... 56
Евклидовы пространства…………………………………………………... 62
Линейные операторы
67
Литература………………………………………………………………….. 77
3
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
ПРЕДИСЛОВИЕ
Учебное пособие “Высшая математика для экономистов” включает
такие разделы высшей математики, изучение которых дает математический
аппарат,
наиболее
активно
применяемый
для
решения
прикладных
экономических и управленческих задач. Это аналитическая геометрия,
линейная алгебра и математический анализ.
В разделе “Линейная алгебра” основное внимание уделяется матрицам,
определителям и системам линейных уравнений, поскольку в экономических
исследованиях широко используются различные матричные модели межотраслевого баланса, в плановых расчетах, при расчетах фонда
заработной платы и т.д. Линейные модели, сводящиеся к системам
алгебраических линейных уравнений или неравенств, c достаточно высокой
точностью соответствуют описываемым ими явлениям; с их помощью
решаются многие управленческие задачи.
В начале каждого параграфа приводятся краткие сведения из теории.
Основное
внимание
уделяется
практическому
освоению
студентами
изучаемого материала. Для достижения этой цели приводится большое число
примеров. Их выполнение будет способствовать выработке навыков
рационального решения типовых примеров и задач, а также задач
экономического и производственного содержания, развивающих навыки
применения изученного математического инструментария.
4
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
Глава 1. МАТРИЦЫ И ОПРЕДЕЛИТЕЛИ
§ 1.1 Матрицы
1. Понятие матрицы. Понятия матрицы и основанный на нем раздел
математики- матричная алгебра – имеют важное значение для экономистов.
Так, как значительная часть
математических моделей экономических
объектов и процессов записывается в достаточно простой и компактной
матричной форме.
Прямоугольная таблица чисел
⎛ a11
⎜
⎜ a 21
⎜ ...
⎜
⎜a
⎝ m1
a12
a 22
...
am 2
a1n ⎞
⎟
... a 2 n ⎟
,
... ... ⎟
⎟
... a mn ⎟⎠
...
(1)
содержащая m строк и n столбцов, называется матрицей размеров m × n .
Числа aij называются элементами матрицы. Каждый элемент матрицы
снабжен двумя индексами: первый индекс указывает номер строки, второй –
номер столбца, в которых расположен этот элемент. В дальнейшем будем
обозначать матрицы большими буквами латинского алфавита: A, B и т.д.
Часто вместо подробной записи (1) употребляют сокращенную:
A = ( aij ) im, j,n=1 или даже A = ( aij ).
Две матрицы A = ( aij ) и B = ( bij ) считаются равными, если совпадают их
размеры и aij = bij при любых i и j .
С помощью матриц удобно записывать некоторые экономические
зависимости. Например, таблица распределения ресурсов по отдельным
отраслям экономики
Ресурсы
Отрасли экономики
промышленность
сельское хозяйство
Электроэнергия
5.6
6.8
Трудовые ресурсы
2.9
3.4
Водные ресурсы
6.7
4.8
5
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
может быть записана в компактной форме в виде матрицы распределения
ресурсов по отраслям:
⎛ 5.6 6.8 ⎞
⎟
⎜
A = ⎜ 2.9 3.4 ⎟ .
3x 2
⎜ 6.7 4.8 ⎟
⎠
⎝
В этой записи, например, матричный элемент, a11 = 5.6 показывает,
сколько электроэнергии потребляет промышленность, а элемент a22 = 3,4 сколько трудовых ресурсов потребляет сельское хозяйство.
2. Квадратные матрицы. Если число строк матрицы равно числу ее
столбцов, то матрица называется квадратной, а число n ее строк (равное
числу столбцов) – порядком квадратной матрицы:
⎛ a11
⎜
⎜ a 21
⎜ ...
⎜
⎜a
⎝ n1
Диагональ
a12
...
a 22
...
...
an 2
...
...
a1n ⎞
⎟
a2n ⎟
.
... ⎟
⎟
a nn ⎟⎠
квадратной
a11 , a22 , ..., ann
матрицы
называется
главной
диагональю, а диагональ a1n , a2,n−1 , ..., an1 – побочной диагональю.
Квадратная матрица называется треугольной, если все ее элементы,
которые находятся ниже (или выше) главной диагонали, равны нулю, т.е.
треугольная матрица имеет вид
⎛ a11
⎜
⎜ 0
A=⎜
...
⎜
⎜ 0
⎝
a12
...
a 22
...
...
...
0
...
a1 n ⎞
⎟
a2n ⎟
... ⎟
⎟
a nn ⎟⎠
или
⎛ b11
⎜
⎜b
B = ⎜ 21
...
⎜
⎜b
⎝ n1
0
...
b22
...
...
...
bn 2
...
0 ⎞
⎟
0 ⎟.
... ⎟
⎟
bnn ⎟⎠
При этом матрицу A называют верхнетреугольной, а матрицу B –
нижнетреугольной.
3. Действия с матрицами.
1) Умножение матрицы на число. Для того чтобы умножить матрицу
A = ( aij ) на число λ , нужно каждый элемент матрицы A умножить на это
число: λA = ( λaij ) .
6
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
2) Сложение матриц(Вычитание матриц). Складывать (Вычитать) можно
только матрицы с одинаковым числом строк и столбцов, т.е. матрицы
одинаковых размеров. Суммой (Разностью) матриц
A = ( aij ) и
B = ( bij )
называется матрица C = ( cij ) , элементы которой равны суммам (разностью)
соответствующих элементов матриц A и B , т.е. cij = aij ± bij для любых
индексов i , j .
3) Умножение
матриц.
Произведение
матрицы
A
на
матрицу
B
(обозначается AB ) определено только в том случае, когда число столбцов
матрицы A равно числу строк матрицы B . В результате умножения получим
матрицу C = AB , у которой столько же строк, сколько их в матрице A , и
столько же столбцов, сколько их в матрице B . Для удобства запоминания
запишем это кратко:
A⋅B =C
( m× n )( n× k ) = m× k
A = ( aij ),
Если
B = ( bij )
и C = AB = ( cij ) , то элементы cij
определяются
следующим образом:
n
cij = ai1b1 j + ai 2b2 j + ... + ainbnj = ∑ aip b pj ,
(2)
p =1
где
i = 1,2,..., m; j = 1,2,..., k .
Это правило можно сформулировать и словесно: элемент cij , стоящий
на пересечении i -й строки и j -го столбца матрицы C = AB , равен сумме
попарных произведений соответствующих элементов i -й строки матрицы A
и
j -го столбца матрицы
B . Другими словами, элемент
cij
является
результатом скалярного произведения i -й вектор-строки и j -го векторстолбца.
В качестве примера применения указанного правила приведем
формулу перемножения квадратных матриц 2-го порядка:
⎛ a11 a12 ⎞⎛ b11 b12 ⎞ ⎛ a11b11 + a12b21 a11b12 + a12b22 ⎞
⎟⎟ .
⎟⎟ = ⎜⎜
⎟⎟⎜⎜
⎜⎜
⎝ a21 a22 ⎠⎝ b21 b22 ⎠ ⎝ a21b11 + a22b21 a21b12 + a22b22 ⎠
7
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
Заметим, что оба произведения AB и BA можно определить лишь в том
случае, когда число столбцов матрицы A совпадает с числом строк матрицы
B , а число строк матрицы A совпадает с числом столбцов матрицы
B. А
именно, матрица A имеет размеры m × n , а B – размеры n × m . При этом,
вообще говоря, AB ≠ BA , например: пусть даны матрицы
⎛1 2
А= ⎜⎜
⎝3 1
⎛1 2
⎜
В= ⎜ 3 1
⎜4 2
⎝
1⎞
⎟,
0 ⎟⎠
⎞
⎟
⎟.
⎟
⎠
Их произведение
⎛11 6 ⎞
⎟⎟,
AB = ⎜⎜
⎝ 6 7⎠
Следовательно
⎛ 7 4 1⎞
⎜
⎟
BA = ⎜ 6 7 3 ⎟ .
⎜10 10 4 ⎟
⎝
⎠
AB ≠ BA .
Из формулы (2) вытекают следующие свойства умножения матриц:
1) ( AB )C = A( BC ) (ассоциативность умножения);
2) ( A + B )C = AC + BC или
A( B + C ) = AB + AC (дистрибутивность умножения относительно сложения).
4) Транспонирование матрицы. Наряду с матрицей
⎛ a11
⎜
⎜a
A = ⎜ 21
...
⎜
⎜a
⎝ m1
a12
a 22
...
...
...
am 2
...
...
a1n ⎞
⎟
a2n ⎟
... ⎟
⎟
a mn ⎟⎠
часто приходится рассматривать матрицу, столбцами которой являются
строки матрицы A (т.е. столбцы и строки меняются местами). Эта матрица
называется транспонированной к A и обозначается через AT :
⎛ a11
⎜
⎜a
T
A = ⎜ 12
...
⎜
⎜a
⎝ 1n
Например. Если
⎛1 2
A = ⎜⎜
⎝3 1
1⎞
⎟ то
0 ⎟⎠
a 21
...
a 22
...
...
...
a2n
...
a m1 ⎞
⎟
am 2 ⎟ .
... ⎟
⎟
a mn ⎟⎠
⎛ 1 3⎞
⎜
⎟
AT = ⎜ 2 1 ⎟ .
⎜ 1 0⎟
⎝
⎠
8
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
5) Возведение в степень. Целой положительной степенью A n (n > 1)
квадратной матрицы A называется произведение n матриц, равных A т.е.
A n = A A ... A .
n − раз
По определению имеем A 0 = E , A1 = A . Нетрудно показать, что
A n A k = A n+k ,
(A )
n k
= A nk .
⎛1 2⎞
⎟⎟ .
Например. Найти A 2 , где A = ⎜⎜
⎝3 4⎠
⎛ 1 2 ⎞⎛ 1 2 ⎞ ⎛ 7 10 ⎞
⎟⎟ .
⎟⎟⎜⎜
⎟⎟ = ⎜⎜
⎝ 3 4 ⎠⎝ 3 4 ⎠ ⎝15 22 ⎠
Решение. A 2 = AA = ⎜⎜
Пример1. Найти произведение матриц
⎛1 2
А= ⎜⎜
⎝3 1
1⎞
⎟ и В=
0 ⎟⎠
⎛1 2 3 ⎞
⎜
⎟
⎜ 2 0 1⎟ .
⎜
⎟
⎝ 3 5 4⎠
Решение. Имеем: матрица А размера 2×3, матрица В размера 3×3, тогда
произведение АВ = С существует и элементы матрицы С равны
с11 = 1⋅1 +2⋅2 + 1⋅3 = 8, с21 = 3⋅1 + 1⋅2 + 0⋅3 = 5,с12 = 1⋅2 + 2⋅0 + 1⋅5 = 7,
с22 =3⋅2 + 1⋅0 + 0⋅5 = 6, с13 = 1⋅3 + 2⋅1 + 1⋅4 = 9,с23 = 3⋅3 + 1⋅1 + 0⋅4 = 10.
⎛8 7 9 ⎞
⎟⎟ ,
⎝ 5 6 10 ⎠
AB = ⎜⎜
а произведение BA не существует.
Пример 2. В таблице указано количество единиц продукции,
отгружаемой ежедневно на молокозаводах 1 и 2 в магазины М1, М2 и М3,
причем
доставка единицы продукции с каждого молокозавода в магазин М1 стоит 50
ден. ед., в магазин М2 - 70, а в М3 - 130 ден. ед. Подсчитать ежедневные
транспортные расходы каждого завода.
Молокозавод
Магазин
М1
М2
М3
1
20
35
10
2
15
27
8
9
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
Решение. Обозначим через А матрицу, данную нам в условии, а через В
- матрицу, характеризующую стоимость доставки единицы продукции в
магазины, т.е.
⎛ 20 35 10⎞
⎟ , В = (50, 70, 130).
⎝ 15 27 8 ⎠
А =⎜
Тогда матрица затрат на перевозки будет иметь вид:
⎛ 20 35 10⎞
АВ = ⎜
⎟
⎝ 15 27 8 ⎠
T
⎛ 50 ⎞
⎟ ⎛ 20 ⋅ 50 + 35 ⋅ 70 + 10 ⋅ 130 ⎞ ⎛ 4750 ⎞
⎜
⎟⎟ .
⎟⎟ = ⎜⎜
⎜ 70 ⎟ = ⎜⎜
⋅
+
⋅
+
⋅
15
50
27
70
8
130
3680
⎝
⎠
⎝
⎠
⎜130 ⎟
⎠
⎝
Итак, первый завод ежедневно тратит на перевозки 4750 ден. ед., второй 3680 ден.ед.
Пример
3.
Швейное
предприятие
производит
зимние
пальто,
демисезонные пальто и плащи. Плановый выпуск за декаду характеризуется
вектором X = (10, 15, 23). Используются ткани четырех типов Т1, Т2, Т3, Т4. В
таблице приведены нормы расхода ткани (в метрах) на каждое изделие.
Вектор С = (40, 35, 24, 16) задает стоимость метра ткани каждого типа, а
вектор P = (5, 3, 2, 2) - стоимость перевозки метра ткани каждого вида.
Изделие
Расход ткани
Т1
Т2
Т3
Т4
Зимнее пальто
5
1
0
3
Демисезонное пальто
3
2
0
2
Плащ
0
0
4
3
1. Сколько метров ткани каждого типа потребуется для выполнения плана?
2. Найти стоимость ткани, расходуемой на пошив изделия каждого вида.
3. Определить стоимость всей ткани, необходимой для выполнения плана.
4. Подсчитать стоимость всей ткани с учетом ее транспортировки.
Решение. Обозначим через А матрицу, данную нам в условии, т. е.
10
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
⎛ 5 1 0 3⎞
⎜
⎟
A = ⎜ 3 2 0 2⎟ ,
⎜
⎟
⎝ 0 0 4 3⎠
тогда для нахождения количества метров ткани, необходимой для
выполнения плана, нужно вектор X умножить на матрицу А:
⎛ 5 1 0 3⎞
⎟
⎜
XА = (10,15, 23) ⎜ 3 2 0 2⎟ = (10 ⋅ 5 + 15 ⋅ 3, 10 ⋅ 1 + 15 ⋅ 2, 23 ⋅ 4, 10 ⋅ 3 + 15 ⋅ 2 + 23 ⋅ 3) =
⎟
⎜
⎝ 0 0 4 3⎠
= (95, 40, 92, 129).
Стоимость ткани, расходуемой на пошив изделия каждого вида, найдем,
перемножив матрицу А и вектор CT:
⎛ 40⎞
⎛ 5 1 0 3 ⎞ ⎜ ⎟ ⎛ 5 ⋅ 40 + 1 ⋅ 35 + 0 ⋅ 24 + 3 ⋅ 16 ⎞ ⎛ 283⎞
⎟
⎜
⎟ 35 ⎜
⎟ ⎜
А CT = ⎜ 3 2 0 2⎟ ⎜⎜ ⎟⎟ = ⎜ 3 ⋅ 40 + 2 ⋅ 35 + 0 ⋅ 24 + 2 ⋅ 16⎟ = ⎜ 222⎟ .
⎟
⎜
⎟ 24 ⎜
⎟ ⎜
⎝ 0 0 4 3⎠ ⎜ ⎟ ⎝ 0 ⋅ 40 + 0 ⋅ 35 + 4 ⋅ 24 + 3 ⋅ 16⎠ ⎝ 144 ⎠
⎝ 16 ⎠
Стоимость всей ткани, необходимой для выполнения плана, определится
по формуле:
⎛ 283⎞
⎟
⎜
X А C = (10,15,23) ⎜ 222⎟ = 10 ⋅ 283 + 15 ⋅ 222 + 23 ⋅ 144 = 9472 .
⎟
⎜
⎝144 ⎠
T
Наконец, с учетом транспортных расходов вся сумма будет равна
стоимости ткани, т. е. 9472 ден. ед., плюс величина
⎛ 5⎞
⎜ ⎟
3
T
X А P = (95, 40, 92, 129) ⎜⎜ ⎟⎟ = 95 ⋅ 5 + 40 ⋅ 3 + 92 ⋅ 2 + 129 ⋅ 2 = 1037 .
2
⎜ ⎟
⎝ 2⎠
Итак,
X А C T + X А P T = 9472 + 1037 = 10509 (ден. ед).
11
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
§ 1.2 Определители
Каждой квадратной матрице A может быть поставлено в соответствие
некоторое число, вычисляемое по определенному правилу с помощью
элементов
матрицы.
Такое
число
называют
определителем
(или
детерминантом) матрицы A и обозначают символом | A | или det A . При этом
порядком определителя называют порядок соответствующей матрицы.
Правила вычисления определителей 2-го и 3-го порядков легко
выписать:
a11 a12
= a11a22 − a12 a21 ,
a21 a22
a11
a12
a21 a22
a31 a32
a13
a23 = a11a22 a33 + a12 a23a31 + a13a21a32 − a11a23a32 − a12 a21a33 − a13a22 a31 .
a33
Последнюю формулу, несмотря на внешнюю сложность записи, нетрудно
запомнить. Если соединить линией каждые три элемента определителя,
произведение которых входит в правую часть последней формулы со знаком
« + », то получим легко запоминающуюся схему 1. Аналогично для
произведений, входящих со знаком «–», имеем схему 2.
Схема 1
Схема 2
Это правило вычисления определителей 3-го порядка называется
правилом треугольников.
Прежде, чем сформулировать определение определителя n -го порядка,
рассмотрим одно вспомогательное понятие.
12
© КРСУ Усенов И.А., Усенова Р.К.
1.
Перестановки.
Кафедра высшей математики
Перестановкой
чисел
(1,2,..., n)
называют
расположение этих чисел в каком-либо определенном порядке (не
обязательно в порядке возрастания). Например, (3,2,1,4) – одна из возможных
перестановок чисел (1,2,3,4) .
Число различных перестановок, которые можно составить из чисел
(1, 2,..., n) , равно произведению 1 ⋅ 2 ⋅ 3 ⋅ ... ⋅ n = n ! (читается: «n факториал»).
Пусть дана какая-то перестановка
J = ( j1 , j2 ,..., jn )
чисел (1, 2,..., n) .
Назовем инверсией (или беспорядком) в перестановке J любую пару чисел в
этой перестановке, из которых большее число расположено левее меньшего.
Пример 1. В перестановке (3,2,1,4) имеются 3 инверсии: их образуют
пары (3,2) , (3,1) , (2,1) .
Условимся обозначать общее число инверсий в перестановке J
символом σ ( J ) . Перестановка J называется четной, если число σ ( J ) –
четное, и нечетной, если число σ ( J ) – нечетное.
Так в рассмотренном выше примере перестановка содержала 3
инверсии и, следовательно, является нечетной. Заметим, что перестановка
(1, 2,..., n) не содержит ни одной инверсии, иначе говоря, содержит 0
инверсий. Следовательно, эта перестановка является четной.
2. Определитель n -го порядка. Определителем n -го порядка (или
определителем матрицы n -го порядка) называется число, равное
a11
a12
... a1n
a21
a22
... a2 n
...
...
...
an1
an 2 ... ann
...
= ∑ (−1)σ ( J ) a1 j1 a2 j2 ... an jn ,
(4)
J
где суммирование распространяется на все перестановки J = ( j1 , j2 ,..., jn ) ,
которые можно составить из чисел (1, 2,..., n) . Количество слагаемых в правой
части равенства (4) равно n ! , так как количество всех перестановок
множества из n элементов равно 1 ⋅ 2 ⋅ 3 ⋅... ⋅ n = n !
13
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
Упражнение 1. В качестве примеров применения общей формулы (4)
рассмотрите формулы для вычисления определителей 2-го и 3-го порядков.
3. Свойства определителей. Сформулируем без доказательства ряд
свойств, которыми обладает произвольный определитель n -го порядка.
1)
Свойство
равноправности
строк
и
столбцов.
При
транспонировании величина определителя сохраняется, т.е. AT =| A | .
Это свойство означает полную равноправность строк и столбцов и
позволяет нам все последующие свойства формулировать лишь для строк и
быть уверенными в справедливости их и для столбцов.
2) При перестановке местами 2-х строк определитель сохраняет свою
абсолютную величину, но меняет знак на противоположный.
3) Линейное свойство определителя. Если все элементы i -ой строки
определителя n -го порядка представлены в виде суммы двух слагаемых
aij = λ aij′ + µ aij′′
( j = 1,2,..., n) ,
то исходный определитель можно представить в виде суммы двух
определителей, у которых элементами i -й строки являются соответственно
aij′ и aij′′ ( j = 1,2,..., n) , а все остальные строки – такие же, как у исходного
определителя. При этом определители умножаются на λ и µ соответственно:
...
a11
...
a11 ... a1n
a1n
...
...
...
...
... λain′ + µain′′ = λ ai′1
...
...
...
a n1
...
λai′1 + µai′′1
...
...
...
... ain′ + µ ai′′1
...
...
a n1 ... a nn
a nn
a11 ... a1n
...
...
...
... ain′′ .
... ...
a n1 ... a nn
Конечно линейное свойство справедливо и для случая, когда i -я строка
является линейной комбинацией не двух, а нескольких строк. На этот случай
(любого
конечного
числа
слагаемых)
сформулированное
свойство
обобщается индукцией по числу слагаемых.
Приведенные
три
свойства
являются
основными
свойствами
определителя. Все следующие свойства являются логическими следствиями
трех основных свойств.
14
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
4) Определитель с двумя одинаковыми строками равен нулю.
В самом деле, при перестановке двух одинаковых строк, с одной
стороны, определитель ∆ не изменится, а с другой стороны, в силу свойства
2 изменит свой знак на противоположный. Таким образом, ∆ = − ∆ , т.е. 2∆ = 0
или ∆ = 0 .
5) Умножение всех элементов некоторой строки определителя на число
λ равносильно умножению определителя на это число λ .
Иными словами, общий множитель всех элементов некоторой строки
определителя можно вынести за знак этого определителя. Это свойство
вытекает из свойства 3 при µ = 0 .
6) Если все элементы некоторой строки определителя равны нулю, то и
сам определитель равен нулю.
Это свойство вытекает из свойства 5 при λ = 0 .
7) Если элементы двух строк определителя пропорциональны, то
определитель равен нулю.
В самом деле, в силу свойства 5 множитель пропорциональности
можно вынести за знак определителя, после чего останется определитель с
двумя одинаковыми строками, который равен нулю согласно свойству 4.
8) Если к элементам некоторой строки определителя прибавить
соответствующие элементы другой строки, умноженные на произвольный
множитель λ , то величина определителя не изменится.
В самом деле, полученный в результате указанного прибавления
определитель можно в силу свойства 3 разбить на сумму двух определителей,
первый из которых совпадает с исходным, а второй равен нулю в силу
пропорциональности двух строк и свойства 7.
9) Определитель произведения матриц. Если C = AB , где A и B –
квадратные матрицы (одинакового порядка), то | C |=| A | ⋅ | B | .
4. Вычисление определителей n -го порядка. Пусть дана матрица A
n -го порядка. Минором любого элемента aij называют определитель порядка
n − 1 , соответствующий той матрице, которая получается из матрицы A в
15
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
результате вычеркивания i -й строки и j -го столбца (т.е. той строки и того
столбца, на пересечении которых стоит элемент aij ). Минор элемента aij
будем обозначать символом M ij .
Алгебраическим дополнением Aij элемента aij матрицы A называют
минор M ij этого элемента, умноженный на (−1)i+ j , т.е.
Aij = (−1)i + j M ij .
Теорема. Определитель матрицы A
n -го порядка равен сумме
произведений всех элементов какой-нибудь одной фиксированной строки на
их алгебраические дополнения, т.е. для любого i = 1,2,..., n имеет место
равенство
| A |= ai1 Ai1 + ai 2 Ai 2 + ... + ain Ain ,
называемое разложением определителя | A | по элементам i -й строки.
Аналогично для k = 1,2,..., n имеет место разложение определителя | A |
по элементам k -го столбца:
| A |= a1k A1k + a2 k A2 k + ... + ank Ank .
Пример 2. Не вычисляя определителя
1 2 3
4 5 6
7 8 9
показать, что он равен нулю.
Решение. Вычтем из второй строки первую, получим определитель
1
2
3
3
3
7
8
3,
9
равный исходному. Если из третьей строки также вычесть первую, то
получится определитель
1 2 3
3 3 3 ,
6 6 6
в котором две строки пропорциональны. Такой определитель равен нулю.
16
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
Пример 3. Вычислить определитель
1
D= 3
4
-2
5
1
3
-1 ,
2
разложив его по элементам второго столбца.
Решение. Разложим определитель по элементам второго столбца:
D = a12A12 + a22A22+a32A32=
= ( −2) ⋅ ( −1)
1+ 2
3 -1
4
2
+ 5 ⋅ ( −1)
2+2
1
3
4
2
+ 1 ⋅ ( −1)
3+ 2
1
3
3 -1
= −20 .
Пример 4. Вычислить определитель
a 11
0
0 ... 0
a 21
a 22
0 ... 0
A = a 31
a 32
a 33 ... 0
,
----------------a n1
a n 2 a n 3 ... a n n
в котором все элементы по одну сторону от главной диагонали равны нулю.
Решение. Разложим определитель А по первой строке:
A = a11 A11 = a 11
a 22
0 ... 0
a 32
a 33 ... 0
.
----------a n 2 a n 3 ...a n n
Определитель, стоящий справа, можно снова разложить по первой строке,
тогда получим:
A = a 11 ⋅ a 22
a 33
0 ... 0
a 43
a 44 ... 0
----------a n 3 a n 4 ...a n n
.
И так далее. После n шагов придем к равенству A = а11 а22... ann.
Пример 5. Вычислить определитель
17
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
1
2
3 ... n
-1 0 3 ... n
-1 - 2 0 ... n
--------------
.
-1 - 2 - 3 ... 0
Решение. Если к каждой строке определителя, начиная со второй,
прибавить первую строку, то получится определитель, в котором все
элементы, находящиеся, ниже главной диагонали, будут равны нулю. А
именно, получим определитель:
1
0
2
2
3 ... n
6 ... 2n
0 0 3 ... 2n ,
--------------0
0
0 ... n
равный исходному.
Рассуждая, как в предыдущем примере найдем, что он равен
произведению элементов главной диагонали, т.е. n!. Способ, с помощью
которого вычислен данный определитель, называется способом приведения к
треугольному виду.
Упражнение 2. Вычислите каждый из следующих определителей двумя
способами (с помощью правила треугольников и с помощью разложения по
элементам строки или столбца):
2 3 4
1 17 − 7
а) 5 − 2 1 , б) − 1 13 1 ,
1 2 3
1 7
1
2 0 5
в) 1 3 16 .
0 − 1 10
5. Обратная матрица. Среди квадратных матриц одного и того же
порядка (например, порядка n , т.е. размеров n × n ) важную роль играет
матрица вида
⎛1
⎜
⎜0
E =⎜
...
⎜
⎜0
⎝
0
1
...
0
...
...
...
...
0⎞
⎟
0⎟
,
...⎟
⎟
1 ⎟⎠
(3)
18
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
которую называют единичной матрицей. Легко проверить, что для любой
матрицы A n -го порядка имеют место равенства
A⋅ E = E ⋅ A = A.
Эти равенства показывают особую роль единичной матрицы E , аналогичную
той роли, которую играет число 1 при перемножении вещественных чисел.
Как известно, для каждого числа a ≠ 0 существует такое число b , что
ab = 1 . Число
b
называется обратным для a . Если мы зафиксируем
натуральное число n и будем рассматривать квадратные матрицы n -го
порядка, то в этом множестве матриц единичная матрица E будет играть
роль единицы. Естественно поставить вопрос о существовании обратной
матрицы, т.е. такой матрицы, которая в произведении с данной матрицей
дает единичную матрицу.
Определение. Пусть A – квадратная матрица n -го порядка. Квадратная
матрица B (того же порядка n ) называется обратной для A , если AB = BA = E .
Матрицу, обратную к матрице A , принято обозначать символом A−1 .
Если для квадратной матрицы A существует обратная матрица A−1 , то
справедливо равенство AA−1 = E , где E – единичная матрица. Переходя в этом
равенстве к определителям (и учитывая свойство 9 определителей), имеем
| A | ⋅ | A−1 |=| E | , или | A | ⋅ | A−1 |= 1 . Отсюда заключаем, что | A |≠ 0 (в противном
случае левая часть последнего равенства равнялась бы нулю). Этим доказано,
что если | A |= 0 , то для матрицы A не существует обратной. Другими
словами, условие | A |≠ 0 является необходимым условием существования
обратной матрицы. Оказывается, это условие является и достаточным.
Определение. Квадратная матрица A называется невырожденной, если
ее определитель не равен нулю ( | A |≠ 0 ). В противном случае матрица A
называется вырожденной ( | A |= 0 ).
Пусть матрица A имеет вид
19
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
⎛ a11
⎜
⎜a
A = ⎜ 21
...
⎜
⎜a
⎝ n1
a12
a 22
...
...
...
an 2
...
...
a1n ⎞
⎟
a2n ⎟
.
... ⎟
⎟
a nn ⎟⎠
Теорема. Если A – невырожденная матрица, то для нее существует
обратная матрица A−1 , которая вычисляется по формуле
A −1
⎛ A11
⎜
1
⎜A
=
= ⎜ 12
...
| A|
⎜
⎜A
⎝ 1n
A21
...
A22
...
...
A2 n
...
...
An1 ⎞
⎟
An 2 ⎟ ,
... ⎟
⎟
Ann ⎟⎠
(5)
где Aij – алгебраическое дополнение для элемента aij матрицы A .
Замечание. Обратим внимание на расположение чисел Aij в правой
части формулы (5): число Aij расположено не в i -й строке и j -м столбце, а
наоборот, в j -й строке и i -м столбце. Таким образом, матрица, стоящая в
правой части (5), является транспонированной матрицей алгебраических
дополнений
элементов
алгебраических
матрицы
дополнений,
A.
Матрица,
называется
составленная
присоединенной
из
матрицей.
Следовательно, нахождение обратной матрицы с помощью формулы (5)
называется методом присоединенной матрицы.
Если матрица немалого размера, то нахождение обратной матрицы
методом присоединенной матрицы составляет практическое неудобство, т.е.
требует больше вычислений. В этом случае используется другой метод,
который называется методом элементарных преобразований или методом
Гаусса.
Метод элементарных преобразований (метод Гаусса) вычисления
обратной матрицы к невырожденной матрице А состоит в следующем.
Приписывая справа к матрице А размера n x n единичную матрицу размера n
x n, получим прямоугольную матрицу Γ = ( A/E ) размера n x2 n. С помощью
элементарных преобразований над строками матрицы Γ сначала приведем ее
к ступенчатому виду Γ1 = ( A1 /B) , где матрица A1 - треугольная, а затем к виду
Γ2 = (E/A -1 ) .
20
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
Пример 6. Для матрицы
⎛2
⎜
А = ⎜2
⎜
⎝1
- 2 1⎞
⎟
1 - 2⎟
⎟
2 2⎠
найти обратную матрицу.
Решение. Находим сначала детерминант матрицы А
2 -2 1
det A = 2
1 - 2 = 27 ≠ 0 ,
1
2 2
значит, обратная матрица существует, и мы ее можем найти по формуле:
⎛ A11 A 21 A 31 ⎞
⎟
1⎜
A = ⎜ A12 A 22 A 32 ⎟ ,
∆⎜
⎟
⎝ A13 A 23 A 33 ⎠
−1
где Аi j (i,j=1,2,3) - алгебраические дополнения элементов аi j исходной
матрицы. Имеем:
A 11 = ( −1) 1+1
1 -2
= 2 + 4 = 6,
A 12 = ( −1) 1+2
= 4 − 1 = 3,
A 21 = ( −1) 2 +1
= 4 − 1 = 3,
A 23 = ( −1) 2 +3
A 31 = ( −1) 3+1
−2 1
= 4 − 1 = 3,
1 -2
A 32 = ( −1) 3+2
A 33 = (−1) 3+3
2 -2
= 2 + 4 = 6.
2 1
A 13 = ( −1) 1+3
A 22 = ( −1) 2 +2
2
2
2
1
1
2
2
1
1
2
2
-2
1
2
= − ( 4 + 2) = −6,
−2
1
2
2
2
-2
1
2
= − ( −4 − 2) = 6,
= − (4 + 2) = −6,
2 1
= − ( −4 − 2) = 6,
2 -2
Откуда
6
⎛6
1 ⎜
А = ⎜- 6 3
27 ⎜
⎝3 - 6
−1
3⎞
⎟
6⎟ =
6 ⎟⎠
⎛2 2
1⎜
⎜- 2 1
9⎜
⎝1 - 2
1⎞
⎟
2⎟ .
2 ⎟⎠
Пример 7. Методом элементарных преобразований найти обратную
матрицу для матрицы:
21
© КРСУ Усенов И.А., Усенова Р.К.
⎛2
⎜
А= ⎜ 5
⎜7
⎝
Кафедра высшей математики
1 -1 ⎞
⎟
2
4⎟ .
3
2 ⎟⎠
Решение. Приписываем к исходной матрице справа единичную
матрицу того же порядка:
⎛ 2 1 -1 1 0 0⎞
⎟
⎜
⎜ 5 2 4 0 1 0⎟ .
⎜ 7 3 2 0 0 1⎟
⎠
⎝
С помощью элементарных преобразований столбцов приведем левую
“половину”
к
единичной,
совершая
одновременно
точно
такие
преобразования над правой матрицей. Для этого поменяем местами первый
и второй столбцы:
⎛ 2 1 -1 1 0 0⎞ ⎛ 1 2 -1 0 1 0⎞
⎟ ⎜
⎟
⎜
⎜ 5 2 4 0 1 0⎟ ∼ ⎜ 2 5 4 1 0 0⎟ .
⎜ 7 3 2 0 0 1⎟ ⎜ 3 7 2 0 0 1⎟
⎠ ⎝
⎠
⎝
К третьему столбцу прибавим первый, а ко второму - первый, умноженный
на -2:
⎛1 0 0 0 1 0⎞
⎜
⎟
⎜2 1 6 1 - 2 1⎟ .
⎜ 3 1 5 0 0 1⎟
⎝
⎠
Из первого столбца вычтем удвоенный второй, а из третьего - умноженный
на 6 второй;
⎛1 0 0 - 2 1 - 6⎞
⎟
⎜
⎜ 0 1 0 5 - 2 13 ⎟ .
⎜1 1 -1 0 0 1 ⎟
⎠
⎝
Прибавим третий столбец к первому и второму:
⎛1 0 0 - 8 - 5 - 6⎞
⎟
⎜
⎜ 0 1 0 18 11 13 ⎟ .
⎜ 0 0 -1 1 1 1 ⎟
⎠
⎝
Умножим последний столбец на -1:
22
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
⎛1 0 0 - 8 - 5 6 ⎞
⎜
⎟
⎜ 0 1 0 18 11 - 13 ⎟ .
⎜ 0 0 -1 1 1 -1 ⎟
⎝
⎠
Полученная справа от вертикальной черты квадратная матрица
является обратной к данной матрице А. Итак,
−1
А
⎛ − 8 - 5 6⎞
⎟
⎜
= ⎜18 11 - 13 ⎟ .
⎜ 1 1 -1 ⎟
⎠
⎝
Проверка:
⎛2
⎜
AA −1 = ⎜ 5
⎜7
⎝
1 - 1 ⎞⎛ − 8 − 5 6 ⎞ ⎛ 1 0 0 ⎞
⎟
⎟⎜
⎟ ⎜
2
4 ⎟⎜ 18 11 − 13 ⎟ = ⎜ 0 1 0 ⎟ = E .
3
2 ⎟⎠⎜⎝ 1
1
− 1 ⎟⎠ ⎜⎝ 0 0 1 ⎟⎠
или
⎛ − 8 − 5 6 ⎞⎛ 2
⎜
⎟⎜
A −1 A = ⎜ 18 11 − 13 ⎟⎜ 5
⎜ 1
− 1 ⎟⎠⎜⎝ 7
1
⎝
1 -1 ⎞ ⎛1 0 0⎞
⎟
⎟ ⎜
2
4⎟ = ⎜ 0 1 0⎟ = E .
3
2 ⎟⎠ ⎜⎝ 0 0 1 ⎟⎠
Упражнение 3. Найти обратные матрицы для следующих матриц и
проверить результат:
7 ⎞
⎛2 5
⎛3 − 4 5 ⎞
⎟
⎟
⎜
⎜
⎛1 2⎞
⎟⎟ , б) ⎜ 6 3
4 ⎟ , в) ⎜ 2 − 3 1 ⎟ .
а) ⎜⎜
⎝3 4⎠
⎜ 5 − 2 − 3⎟
⎜ 3 − 5 − 1⎟
⎠
⎠
⎝
⎝
6. Ранг матрицы. Ранее для квадратной матрицы A = (aij ) n -го порядка
было введено понятие минора M ij элемента aij . Напомним, что так был
назван определитель порядка n − 1 , полученный из определителя | A |
вычеркиванием i -й строки и j -го столбца.
Введем теперь общее понятие минора. Рассмотрим некоторую, не
обязательно квадратную матрицу A . Выберем какие-нибудь s номеров
строк i1 , i2 , ..., is и s номеров столбцов j1 , j2 ,..., js .
Определение. Минором порядка s матрицы A (соответствующим
выбранным строкам и столбцам) называется определитель порядка s ,
23
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
образованный элементами, стоящими на пересечении выбранных строк и
столбцов, т.е. число
M
i1 ...is
j1 ... js
=
ai1 j1
ai1 j2
ai2 j1
ai2 j2
...
ais j1
...
ais j2
... ai1 js
... ai2 js
.
... ...
... ais js
Каждая матрица имеет столько миноров данного порядка s , сколькими
способами можно выбрать номера строк i1 , i2 , ..., is и столбцов j1 , j2 ,..., js .
Определение. В матрице
A размеров
m× n
минор порядка
s
называется базисным, если он отличен от нуля, а все миноры порядка s + 1
равны нулю или миноров порядка s + 1 вообще нет, т.е. s совпадает с
меньшим из чисел m или n .
Ясно, что в матрице может быть несколько разных базисных миноров.
Все базисные миноры имеют один и тот же порядок. Действительно, если все
миноры порядка s + 1 равны нулю, то равны нулю и все миноры порядка
s + 2 , а, следовательно, и всех больших порядков. Это становится
очевидным, если разложить минор порядка s + 2 по элементам какой-либо
строки
(столбца):
все
миноры
элементов
этой
строки
являются
определителями порядка s + 1 , а поэтому равны нулю.
Определение. Рангом матрицы называется порядок базисного минора,
или, иначе, самый большой порядок, для которого существуют отличные от
нуля миноры. Если все элементы матрицы равны нулю, то ранг такой
матрицы, по определению, считают нулем.
Ранг матрицы A будем обозначать символом r ( A) . Из определения
ранга следует, что для матрицы A размеров m × n справедливо соотношение
0 ≤ r ( A) ≤ min( m, n) .
Два способа вычисления ранга матрицы.
а) Метод окаймляющих миноров
Пусть в матрице найден минор M k k -го порядка, отличный от нуля.
Рассмотрим лишь те миноры (k + 1) -го порядка, которые содержат в себе
24
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
(окаймляют) минор M : если все они равны нулю, то ранг матрицы равен k .
В противном случае среди окаймляющих миноров найдется ненулевой минор
( k + 1) -го порядка, и вся процедура повторяется. Ясно, что перебирать таким
способом миноры в поисках базисного – задача, связанная с большими
вычислениями, если размеры матрицы не очень малы. Существует, однако,
более простой способ нахождения ранга матрицы – при помощи
элементарных преобразований.
б) Метод элементарных преобразований
Определение. Элементарными преобразованиями матрицы называют
следующие преобразования:
1) умножение строки на число, отличное от нуля;
2) прибавление к одной строке другой строки;
3) перестановку строк;
4) такие преобразования столбцов.
Преобразования 1 и 2 выполняются поэлементно.
Комбинируя преобразования первого и второго вида, мы можем к
любой строке прибавить линейную комбинацию остальных строк.
Теорема. Элементарные преобразования не меняют ранга матрицы.
Идея практического метода вычисления ранга матрицы
⎛ a11
⎜
⎜a
A = ⎜ 21
...
⎜
⎜a
⎝ m1
a12
a 22
...
...
...
am 2
...
...
a1n ⎞
⎟
a2n ⎟
... ⎟
⎟
a mn ⎟⎠
заключается в том, что с помощью элементарных преобразований
данную матрицу A приводят к виду
⎛ b11
⎜
⎜ 0
⎜ ...
⎜
B=⎜ 0
⎜ 0
⎜
⎜ ...
⎜
⎝ 0
b12
...
b1 r
b1, r +1
...
b 22
...
b2 r
b 2 , r +1
...
...
...
...
...
...
0
...
b rr
b r , r +1
...
0
...
0
0
...
...
0
...
...
...
0
...
0
...
...
b1 n ⎞
⎟
b2 n ⎟
... ⎟
⎟,
br n ⎟
0 ⎟⎟
... ⎟
⎟
0 ⎠
25
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
в котором «диагональные» элементы b11 , b22 , ..., brr отличны от нуля, а
элементы, расположенные ниже «диагональных», равны нулю. Условимся
называть матрицу
B такого вида треугольной (иначе, ее называют
диагональной, трапециевидной или лестничной). После приведения матрицы
A к треугольному виду можно сразу записать, что
r ( A) = r .
В самом деле, r ( A) = r ( B ) (т.к. элементарные преобразования не меняют
ранга). Но у матрицы B существует отличный от нуля минор порядка r :
b11 b12
0 b22
... ...
0
0
а любой минор порядка
r +1
... b1r
... b2 r
= b11 ⋅ b22 ⋅...⋅ brr ≠ 0 ,
... ...
... brr
содержит нулевую строку и поэтому равен нулю.
Сформулируем теперь практическое правило вычисления ранга матрицы
A с помощью элементарных преобразований: для нахождения ранга
матрицы A следует с помощью элементарных преобразований привести ее к
треугольному виду B . Тогда ранг матрицы A будет равен числу ненулевых
строк в полученной матрице B .
Пример 8. Найти методом окаймления миноров ранг матрицы
⎛1 2 - 1 - 2 ⎞
⎜
⎟
0⎟ .
⎜2 4 3
⎜-1 - 2 6 6 ⎟
⎝
⎠
Решение. Начинаем с миноров 1-го порядка, т.е. с элементов матрицы
А. Выберем, например, минор (элемент) М1 = 1, расположенный в первой
строке и первом столбце. Окаймляя при помощи второй строки и третьего
столбца, получаем минор
M2 =
1
-1
2
3
,
отличный от нуля. Переходим теперь к минорам 3-го порядка, окаймляющим
М2. Их всего два (можно добавить второй столбец или четвертый).
Вычисляем их:
26
© КРСУ Усенов И.А., Усенова Р.К.
1
2 -1
2
4 3 = 0,
-1 - 2 6
Кафедра высшей математики
1 -1 - 2
2 3 0 =0 .
-1 6 6
Таким образом, все окаймляющие миноры третьего порядка оказались
равными нулю. Ранг матрицы А равен r ( A) = 2 .
Пример 9. Найти методом элементарных преобразований ранг
матрицы
⎛2
⎜
А= ⎜ 3
⎜5
⎝
3 5 − 3 − 2⎞
⎟
4 3 -1 - 3 ⎟
6 - 1 3 - 5 ⎟⎠
и привести ее к каноническому виду.
Решение. Из второй строки вычтем первую и переставим эти строки:
⎛1 1 - 2 2 − 1 ⎞
⎜
⎟
⎜ 2 3 5 - 3 - 2⎟ .
⎜ 5 6 -1 3 - 5 ⎟
⎝
⎠
Теперь из второй и третьей строк вычтем первую, умноженную
соответственно на 2 и 5:
⎛1 1 - 2 2 − 1 ⎞
⎟
⎜
⎜0 1 9 - 7 0 ⎟ ;
⎜0 1 9 - 7 0 ⎟
⎠
⎝
из третьей строки вычтем первую; получим матрицу
⎛1 1 - 2 2 − 1 ⎞
⎜
⎟
B = ⎜0 1 9 - 7 0 ⎟ ,
⎜0 0 0 0 0⎟
⎝
⎠
которая эквивалентна матрице А, так как получена из нее с помощью
конечного множества элементарных преобразований. Очевидно, что ранг
матрицы В равен r ( B ) = 2 , а следовательно, и r ( A) = 2 .
Матрицу В легко привести к канонической вид. Вычитая первый
столбец, умноженный на подходящие числа, из всех последующих, обратим
в нуль все элементы первой строки, кроме первого, причем элементы
остальных
строк
не
изменяются.
Затем,
вычитая
второй
столбец,
27
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
умноженный на подходящие числа, из всех последующих, обратим в нуль все
элементы второй строки, кроме второго, и получим каноническую матрицу:
⎛1
⎜
⎜0
⎜
⎝0
0
1
0
0
0
0
0
0
0
0⎞
⎟
0⎟.
⎟
0⎠
Упражнение 4. Найти методом элементарных преобразований и
методом окаймления миноров ранг матрицы и сравнит результат
⎛ − 2 3 0 − 1⎞
⎟
⎜
A = ⎜ − 3 −1 3 2 ⎟ .
⎜−8 1 6 3 ⎟
⎠
⎝
§ 1.3 Понятие вектора. Векторы в экономике
Упорядоченный набор чисел, называется вектором. Например,
следующие наборы чисел (-2, 1), (3, 0, 2), (5, 3, -1, 3) есть векторы. Векторы
вкратце записывают одной буквой. Например, А=(-2, 1), В=(3, 0, 2), С=(5, 3,
-1, 3). Число задающие вектор называют компонентами или координатами
вектора. Например, в векторе В=(3, 0, 2) число 3 есть 1-я компонента, число 0
– 2-я, число 2 – 3-я компонента. Число компонент вектора называется его
размерностью. Например, выше записанный вектор А – двумерный, В –
трехмерный, С– четырехмерный векторы.
n-мерный вектор – это упорядоченный набор из n чисел. В общем виде
n-мерный вектор записывается в виде
А=(а1, а2, …, аn)
Из одних и тех же чисел можно образовать различные векторы.
Например, из чисел 2, 1, 6 можно образовать векторы А1=(2, 1, 6), А2=(2, 6,
1), А3=(1, 2, 6), А4=(6, 1, 2). Поэтому в данном векторе менять, местами
компоненты (координаты) нельзя. Иногда в обозначении вектора используют
символы
и т.д.
Векторы широко используются в экономике. Пусть некая фирма
выпускает три вида продукции. За год продукция 1-го вида выпущена 5 тыс.
28
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
штук, 2-го – 3 тыс. штук, а 3-го вида на 7 тыс. штук. Тогда валовой выпуск за
год можно записать в виде вектора А=(5, 3, 7).
Пусть произведенная продукция 1-го вида реализуется по 100 сомов за
штуку,
2-го по 150 сомов, а 3-го по 80 сомов за штуку. Тогда набор Р=(100,
150, 80) есть вектор цен.
Вектор, все координаты которого равны нулю называется нулевым и
обозначается 0 = (0, 0, …, 0). Таким образом векторы можно записать в виде
строки А=(а1, а2, …, аn) или в виде столбца
B=
Строчная запись вектора можно перевести в столбцовую запись и наоборот.
Такая операция называется транспонированием вектора и обозначения
символом АТ.
Если
А=( а1, а2, …, аn)
то
АТ=
.
Если
В=
,
то
ВТ=(b1, b2, …, bn).
Легко видеть, что (АТ)Т=А
29
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
Глава 2. СИСТЕМЫ ЛИНЕЙНЫХ АЛГЕБРАИЧЕСКИХ УРАВНЕНИЙ
§ 2.1 Системы линейных уравнений
1. Общие понятия. Линейным (относительно неизвестных x1 , x2 , ..., xn )
называют алгебраическое уравнение первой степени, т.е. уравнение вида
a1 x1 + a2 x2 + ... + an xn = b , где a1 , ..., an , b – числа. Причина такого названия в том,
что уравнение первой степени с двумя переменными ax + by = c определяет на
плоскости в декартовой прямоугольной системе координат прямую линию.
Система m линейных уравнений с n неизвестными имеет вид
⎧a11 x1 + a12 x 2 + ... + a1n x n = b1 ,
⎪a x + a x + ... + a x = b ,
⎪ 21 1
22 2
2n n
2
⎨
..........
..........
..........
..........
.......
⎪
⎪⎩a m1 x1 + a m 2 x 2 + ... + a mn x n = bm .
(6)
В общем случае число уравнений в системе не обязательно совпадает с
числом неизвестных: m может быть меньше, равно или больше числа n .
Числа
(вещественные
aij
или
комплексные)
называются
коэффициентами системы (6); bi – свободными членами; x1 , x2 , ..., xn –
неизвестными.
Систему (6) можно записать в матричной форме:
AX =B,
где
⎛ a11
⎜
⎜a
A = ⎜ 21
...
⎜
⎜a
⎝ m1
a12
a 22
...
...
...
am 2
...
...
a1n ⎞
⎟
a2n ⎟
,
... ⎟
⎟
a mn ⎟⎠
⎛ x1 ⎞
⎜ ⎟
⎜x ⎟
X =⎜ 2⎟,
#
⎜ ⎟
⎜x ⎟
⎝ n⎠
⎛ b1 ⎞
⎜ ⎟
⎜b ⎟
B = ⎜ 2 ⎟.
#
⎜ ⎟
⎜b ⎟
⎝ m⎠
Если b1 = b2 = ... = bm = 0 , то система называется однородной, в противном
случае она называется неоднородной.
Определение. Совокупность n чисел α1 , α 2 ,..., α n называется решением
системы (6), если после замены неизвестных x1 , x2 , ..., xn числами α1 , α 2 ,..., α n
соответственно каждое из уравнений системы превращается в верное
равенство.
30
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
Пример 1. Рассмотрим систему линейных уравнений
⎧ x + y + z = 1,
⎨
⎩2 x + 2 y + 2 z = 3 .
Эта система 2-х уравнений с тремя неизвестными решений не имеет, так как
любая тройка чисел, удовлетворяющая первому уравнению, не может
удовлетворять второму.
Пример 2. Система
⎧ x + y = 1,
⎨
⎩2 x + 7 y = −3.
имеет единственное решение x = 2 , y = −1 .
Пример 3. Рассмотрим систему линейных уравнений
⎧ x + y = 1,
⎪
⎨2 x + 2 y = 2 ,
⎪3x + 3 y = 3 .
⎩
Пара чисел (1; 0) есть одно из решений этой системы трех уравнений с
двумя неизвестными, (−1; 2) – другое решение. Эта система имеет бесконечно
много решений: значения
x=C,
y =1− C
при любом действительном
значении C удовлетворяют данной системе.
Рассмотренные примеры систем линейных уравнений показывают, что,
вообще говоря, система может либо вовсе не иметь решений, либо иметь
единственное решение, либо иметь их несколько (в последнем случае,
оказывается, система всегда имеет бесконечное множество решений).
Определение. Система линейных уравнений, не имеющая ни одного
решения, называется несовместной. Система, обладающая хотя бы одним
решением, называется совместной.
Относительно каждой системы линейных уравнений могут быть
поставлены следующие вопросы:
1) Совместна заданная система или нет?
2) В случае, если система совместна, как определить, сколько она имеет
решений – одно или несколько?
31
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
3) Как найти все решения системы?
Ответ на все эти вопросы дает теория систем линейных уравнений.
2. Правило Крамера. Ограничимся сначала рассмотрением систем, у
которых число уравнений равно числу неизвестных (такие системы называют
квадратными).
Пусть дана система n линейных уравнений с n неизвестными:
⎧a11 x1 + a12 x 2 + ... + a1n x n = b1 ,
⎪a x + a x + ... + a x = b ,
⎪ 21 1
22 2
2n n
2
⎨
⎪...............................................
⎪⎩a n1 x1 + a n 2 x 2 + ... + a nn x n = bn .
(7)
Определитель
a11 a12
a
a22
∆ = 21
... ...
an1 an 2
составленный
из
коэффициентов
... a1n
... a2 n
,
... ...
... ann
при
неизвестных,
называется
определителем системы (7).
Теорема. Если определитель ∆ квадратной системы (6) отличен от
нуля, то эта система имеет единственное решение. Это решение может
быть найдено по формулам
x1 =
∆1
,
∆
x2 =
∆
∆2
, ..., xn = n ,
∆
∆
где ∆ k – определитель, получаемый из определителя ∆ заменой k -го
столбца на столбец свободных членов.
Формулы для неизвестных носят название формул Крамера.
Пример 4. Решить систему уравнений по формуле Крамера
⎧ x1 + 2 x 2 + 3x3 = 6 ,
⎪
⎨4 x1 + 5 x 2 + 6 x3 = 9 ,
⎪7 x + 8 x
= −6 .
2
⎩ 1
Вычислим основной определитель системы
32
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
1 2 3
∆ = 4 5 6 = 27.
7 8 0
Так, как ∆ ≠ 0 то решение системы существует и единственно. Найдем
дополнительные определители системы
6 2 3
∆ 1 = 9 5 6 = −54,
−6 8 0
1 6 3
∆ 2 = 4 9 6 = 27,
7 −6 0
1 2 6
∆ 3 = 4 5 9 = 54.
7 8 −6
Отсюда по формуле Крамера получаем решение системы уравнений
x1 =
∆1 − 54
=
= −2,
∆
27
x2 =
∆ 2 27
=
= 1,
∆ 27
x3 =
∆ 3 54
=
= 2.
∆ 27
3. Матричный метод. Если матрица А системы линейных уравнений
невырожденная, т.е. det A ≠ 0 , то матрица А имеет обратную матрицу, и
решение системы (7) совпадает с вектором C = A −1 B . Иначе говоря, данная
система имеет единственное решение. Отыскание решения системы по
формуле X = C , C = A −1 B называют матричным способом решения системы,
или решением по методу обратной матрицы.
Пример 5. Решить матричным способом систему уравнений
⎧ x1 − x 2 + x3 = 6 ,
⎪
⎨2 x1 + x 2 + x3 = 3 ,
⎪x + x + 2 x = 5.
2
3
⎩ 1
Решение. Обозначим
⎛1 - 1 1 ⎞
⎟
⎜
T
A = ⎜ 2 1 1⎟ , X = (x1 , x2 , x3 ) , B = (6, 3, 5)T .
⎜1 1 2 ⎟
⎠
⎝
Тогда данная система уравнений запишется матричным уравнением
AX =B.
⎛1 - 1 1 ⎞
⎟
⎜
Поскольку ∆ = det A = ⎜ 2 1 1⎟ = 5 ≠ 0 , то матрица A невырождена и поэтому
⎜1 1 2 ⎟
⎠
⎝
имеет обратную матрицу:
33
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
⎛ A11 A 21 A 31 ⎞
⎟
1⎜
A = ⎜ A 12 A 22 A 32 ⎟ .
∆⎜
⎟
⎝ A 13 A 23 A 33 ⎠
−1
Для получения решения X мы должны умножить вектор-столбец B слева
на матрицу A-1: X = A−1B. В данном случае
−1
A
⎛ 1 3 - 2⎞
⎟
1⎜
= ⎜ -3 1 1⎟
5⎜
⎟
⎝ 1 - 2 3⎠
и, следовательно,
⎛ x1 ⎞
⎛ 1 3 - 2 ⎞ ⎛6 ⎞
⎟⎜ ⎟
⎜ ⎟ 1⎜
x
3
1
1
=
⎜
⎟ ⎜3 ⎟ .
2
⎜ ⎟
⎜ ⎟ 5 ⎜ 1 - 2 3 ⎟ ⎜5 ⎟
⎝ x3 ⎠
⎝
⎠⎝ ⎠
Выполняя действия над матрицами, получим:
⎛ x1 ⎞
⎛ 5 ⎞
⎟
⎜ ⎟ 1⎜
⎜ x 2 ⎟ = ⎜ − 10 ⎟ или
⎜ ⎟ 5 ⎜ 15 ⎟
⎝ x3 ⎠
⎝
⎠
x1 = 1
x 2 = −2
x3 = 3
Итак, C = (1,−2,3)T .
4.
Критерий
совместности
системы
линейных
уравнений.
Рассмотрим снова произвольную систему m линейных уравнений с n
неизвестными, которую запишем, как и раньше, в матричной форме:
AX =B,
где
⎛ a11
⎜
⎜a
A = ⎜ 21
...
⎜
⎜a
⎝ m1
a12
a 22
...
...
...
am 2
...
...
a1n ⎞
⎟
a2n ⎟
,
... ⎟
⎟
a mn ⎟⎠
⎛ x1 ⎞
⎜ ⎟
⎜x ⎟
X =⎜ 2⎟,
#
⎜ ⎟
⎜x ⎟
⎝ n⎠
(8)
⎛ b1 ⎞
⎜ ⎟
⎜b ⎟
B = ⎜ 2 ⎟.
#
⎜ ⎟
⎜b ⎟
⎝ m⎠
Матрицу A называют матрицей системы (8), а матрицу, полученную из
матрицы A добавлением столбца свободных членов B , – расширенной
матрицей системы (8). Обозначим расширенную матрицу системы (8)
символом ( A | B ) :
34
© КРСУ Усенов И.А., Усенова Р.К.
⎛ a11
⎜
⎜a
( A | B) = ⎜ 21
...
⎜
⎜a
⎝ m1
Кафедра высшей математики
a12
a 22
...
am2
... a1n
... a 2 n
... ...
... a mn
b1 ⎞
⎟
b2 ⎟
.
... ⎟
⎟
bm ⎟⎠
Очевидно, что ранги матриц A и ( A | B ) связаны неравенством
r ( A | B ) ≥ r ( A) .
Упражнение. Ранг матрицы
( A | B ) может быть лишь на единицу
больше ранга матрицы A . Объясните, почему.
Вопрос о совместности системы (8) полностью решается следующей
теоремой.
Теорема Кронекера-Капелли. Для того чтобы система линейных
уравнений была совместна, необходимо и достаточно, чтобы ранг матрицы
этой системы был равен рангу ее расширенной матрицы, т.е. чтобы
r ( A | B ) = r ( A) .(Без доказательства)
Пример 6. Исследовать систему уравнений и решить ее, если она
совместна:
⎧5 x1 − x 2 + 2 x3 + x 4 = 7 ,
⎪
⎨2 x1 + x 2 + 4 x3 − 2 x 4 = 1,
⎪ x − 3x − 6 x + 5 x = 0 .
2
3
4
⎩ 1
Решение. Выписываем основную и расширенную матрицу системы:
⎛5 -1 2 1
⎜
A = ⎜2 1
4 -2
⎜1 - 3 - 6 5
⎝
⎞
⎟
⎟,
⎟
⎠
⎛5 -1 2 1
⎜
⎯ (A / B) = ⎜ 2 1 4 - 2
⎜1 - 3 - 6 5
⎝
7⎞
⎟
1⎟ .
0 ⎟⎠
Вычислим ранг основной матрицы системы. Очевидно, что, например, минор
второго порядка в левом верхнем углу
5
2
-1
= 7 ≠ 0; содержащие его
1
миноры третьего порядка равны нулю:
5 -1 2
M′3 = 2 1 4 = 0,
1 -3 -6
5 -1 2
M″3 = 2 1 4 = 0.
1 -3 -6
35
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
Следовательно, ранг основной матрицы системы равен r ( A) = 2 . Для
вычисления ранга расширенной матрицы ( A / B ) рассмотрим окаймляющий
минор
5 -1
2 1
1 -3
7
5 14
1 = 2 7
0 1 0
7
1 = -35 ≠ 0,
0
значит, ранг расширенной матрицы r ( A / B ) = 3 . Поскольку r ( A / B ) ≠ r ( A) , то
система несовместна.
5.
Метод
Гаусса
(метод
последовательного
исключения
неизвестных) решения систем линейных уравнений.
Под элементарными преобразованиями системы линейных уравнений
понимаются следующие операции:
1) умножение какого-либо уравнения системы на число, отличное от нуля;
2) прибавление к одному уравнению другого уравнения;
3) перемена местами уравнений в системе.
Комбинируя элементарные преобразования первого и второго типов,
мы можем к любому уравнению прибавить другое уравнение, умноженное на
произвольное число.
Производя элементарные преобразования в системе, мы получаем
новую систему. Очевидно, что каждому элементарному преобразованию
системы
соответствуют
аналогичные
преобразования
над
строками
расширенной матрицы этой системы, и наоборот, каждому элементарному
преобразованию строк расширенной матрицы соответствует некоторое
элементарное преобразование в системе. Таким образом, элементарные
преобразования в системе сводятся к соответствующим преобразованиям над
строками ее расширенной матрицы.
Определение. Две системы линейных уравнений от одних и тех же
неизвестных называются равносильными, если каждое решение одной из них
является решением другой, и наоборот (или если обе системы несовместны).
36
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
Заметим, что число уравнений в равносильных системах может быть
различным.
Теорема. При элементарных преобразованиях система линейных
уравнений переходит в равносильную систему.
Сущность метода Гаусса заключается в том, что с помощью
элементарных преобразований система уравнений приводится к такому виду,
чтобы матрица системы оказалась треугольной. Для упрощения изложения
мы будем иметь дело не с самой системой (6), а с расширенной матрицей
этой системы (производя при этом элементарные преобразования только над
строками матрицы).
Рассмотрим алгоритм применения метода Гаусса на простых примерах.
Пример 7. Решить систему
⎧ x + y = 1,
⎨
⎩− x − y = 2.
Решение. Преобразуем расширенную матрицу системы:
⎛ 1 1 1⎞
⎜⎜
⎟⎟ ~
⎝ −1 −1 2⎠
⎛ 1 1 1⎞
⎜⎜
⎟⎟ .
⎝ 0 0 3⎠
Отсюда следует, что r ( A) = 1 , r ( A | B ) = 2 , т.е. исходная система несовместна.
Заметим, что, применяя метод Гаусса (т.е. исключая неизвестные), мы
одновременно проводим исследование системы на совместность (т.е.
отыскиваем ранги матрицы системы и расширенной матрицы).
Пример 8. Решить систему
⎧ x + y = 1,
⎨
⎩− x − y = −1.
Решение. Исследуем систему на совместность:
⎛1 1 1⎞
⎟⎟ ~
⎜⎜
⎝ − 1 − 1 − 1⎠
⎛1 1 1⎞
⎟⎟ .
⎜⎜
⎝ 0 0 0⎠
Отсюда следует, что r ( A) = r ( A | B ) = 1 – система совместна.
Итак, полученная система, равносильная исходной, содержит одно
уравнение с двумя неизвестными. Решение этой системы может быть
найдено
только
в
том
случае,
если
мы
придадим
произвольное
37
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
действительное значение одному из неизвестных. Тогда другое неизвестное
можно выразить через первое.
Заметим, что неизвестные, значения которых можно выбирать
произвольно,
называют
свободными.
Число
свободных
неизвестных
определяется по формуле n − r , где n – число неизвестных в исходной
системе, r – ранг матрицы системы (совпадающий с рангом расширенной
матрицы в силу совместности системы).
В данном случае n − r = 1 . Положим y = C ; тогда x = 1 − C . В итоге
получаем общее решение системы:
⎛1 − C ⎞
⎟⎟ ,
x = ⎜⎜
⎝ C ⎠
где C – произвольная постоянная.
Придавая
постоянной
C
различные
действительные
значения,
получаем бесконечное множество решений исходной системы.
При желании можно произвести проверку:
⎛ 1 1 ⎞⎛1 − C ⎞ ⎛ 1 − C + C ⎞ ⎛ 1 ⎞
⎜⎜
⎟⎟⎜⎜
⎟⎟ = ⎜⎜
⎟⎟ = ⎜⎜ ⎟⎟ .
⎝ − 1 − 1⎠⎝ C ⎠ ⎝ − 1 + C − C ⎠ ⎝ − 1⎠
Пример 9. Решить систему
⎧ x − 2 y + 3z = 0,
⎨
⎩2 x + y − z = 0.
Решение. Данная система apriori является совместной, т.к. она
однородна (все свободные члены равны нулю). Однородная система всегда
имеет нулевое (или тривиальное) решение: x = y = z = 0 . Для однородных
систем особый интерес представляет вопрос о существовании ненулевых (или
нетривиальных) решений. Так называют всякое решение системы, у которого
значение хоть одного неизвестного отлично от нуля.
Преобразуем расширенную матрицу системы:
⎛1 − 2 3 0⎞
⎜⎜
⎟⎟
⎝ 2 1 −1 0⎠
C 2 − 2 C1
⎛1 − 2 3 0⎞
⎟⎟ .
⎝0 5 − 7 0⎠
~ ⎜⎜
38
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
Имеем r ( A) = r ( A | B ) = 2 – система совместна. Тогда n − r = 3 − 2 = 1 – количество
свободных неизвестных. Полагая z = C (где C – произвольная постоянная),
получим
⎧ x − 2 y + 3C = 0,
⎨
⎩ 5 y − 7C = 0.
Отсюда x = −
C
7C
, y=
. Таким образом, общее решение системы имеет вид
5
5
⎛ x ⎞ ⎛ − C 5⎞
⎟
⎜ ⎟ ⎜
X = ⎜ y ⎟ = ⎜ 7C 5 ⎟ ,
⎜z⎟ ⎜ C ⎟
⎠
⎝ ⎠ ⎝
где C – произвольная постоянная.
Как мы отмечали ранее, система линейных уравнений может либо
вовсе не иметь решений, либо иметь единственное решение, либо иметь их
несколько (в последнем случае, оказывается, система всегда имеет
бесконечное множество решений). Это утверждение исчерпывает все
возможные ситуации. При изложении (на примерах) метода Гаусса мы
получили возможность эвристически дать ответ на вопрос о числе решений
(в случае совместности системы). Строгий ответ на этот вопрос дает
следующая теорема.
Теорема (о числе решений). Пусть для системы m линейных уравнений
с n неизвестными выполнено условие совместности, т.е. ранг r матрицы
системы равен рангу ее расширенной матрицы. Тогда, если ранг матрицы
системы равен числу неизвестных ( r = n ), то система имеет единственное
решение. Если же ранг матрицы системы меньше числа неизвестных ( r < n ),
то система имеет бесконечное множество решений, а именно: некоторым
n−r
неизвестным можно придавать произвольные значения, тогда
оставшиеся r неизвестных определятся уже единственным образом.
6. Системы линейных уравнений общего вида. Если система (6)
оказалась совместной, т. е. матрицы A и ( A | B ) имеют один и тот же ранг, то
могут представиться две возможности.
39
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
а) Если r = n , то имеем n независимых уравнений с n неизвестными,
причем определитель ∆ этой системы отличен от нуля. Такая система имеет
единственное решение, получаемое по формулам Крамера.
б) Если
r < n,
то число независимых уравнений меньше числа
неизвестных.
Перенесем
лишние
неизвестные,
которые
x r +1 , x r + 2 ,..., x n
принято
называть свободными, в правые части; наша система линейных уравнений
примет вид:
⎧a11 x1 + a12 x 2 + ... + a1r x r = b1 − a1r +1 x r +1 − a1n x n ,
⎪a x + a x + ... + a x = b − a x − a x ,
⎪ 21 1
22 2
2r r
2
2 r +1 r +1
2n n
⎨
⎪...............................................
⎪⎩a r1 x1 + a r 2 x 2 + ... + a rr x r = bm − a mr +1 x r +1 − a mn x n .
Ее можно решить относительно x1, x2,..., xr, так как определитель этой
системы (r-го порядка) отличен от нуля. Придавая свободным неизвестным
произвольные
числовые
значения,
получим
по
формулам
Крамера
соответствующие числовые значения для x1, x2,..., xr. Таким образом, при r < n
имеем бесчисленное множество решений.
Система (6) называется однородной, если все bi = 0, т. е. она имеет вид:
⎧a11 x1 + a12 x 2 + ... + a1n x n = 0,
⎪a x + a x + ... + a x = 0,
⎪ 21 1
22 2
2n n
⎨
⎪...............................................
⎪⎩a m1 x1 + a m 2 x 2 + ... + a mn x n = 0.
(9)
Из теоремы Кронекера-Капелли следует, что она всегда совместна, так
как добавление столбца из нулей не может повысить ранга матрицы. Это,
впрочем, видно и непосредственно - система (9) заведомо обладает нулевым,
или тривиальным, решением x1 = x2 = ... = xn = 0 . Пусть матрица А системы (9)
имеет ранг r.
Если r = n , то нулевое решение будет единственным решением системы
(9); при r < n система обладает решениями, отличными от нулевого, и для их
разыскания применяют тот же прием, как и в случае произвольной системы
уравнений.
40
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
Пример 10. Исследовать систему уравнений и решить ее, если она
совместна.
⎧ x1 + x 2 − 2 x3 − x 4 + x5 = 1,
⎪
⎨3x1 − x 2 + x3 + 4 x 4 + 3x5 = 4 ,
⎪x + 5x − 9 x − 8x + x = 0 .
2
3
4
5
⎩ 1
Решение. Будем находить ранги матриц A и
элементарных
преобразований,
приводя
( A / B ) методом
одновременно
систему
к
ступенчатому виду:
⎛1 1 - 2 - 1 1 1 ⎞
⎜
⎟
⎜ 3 -1 1 4 3 4 ⎟ ∼
⎜1 5 - 9 - 8 1 0 ⎟
⎝
⎠
Очевидно,
что
⎛ 1 1 - 2 - 1 1 1 ⎞ ⎛1 1 - 2 - 1 1 1 ⎞
⎜
⎟ ⎜
⎟
⎜ 0 - 4 7 7 0 1⎟ ∼ ⎜ 0 - 4 7 7 0 1 ⎟ .
⎜ 0 4 - 7 - 7 0 -1 ⎟ ⎜ 0 0 0 0 0 0 ⎟
⎝
⎠ ⎝
⎠
r ( A) = r ( A / B ) = 2 .
Исходная
система
равносильна
следующей системе, приведенной к ступенчатому виду:
⎧ x1 + x 2 − 2 x3 − x 4 + x5 = 1,
⎨
⎩− 4 x 2 + 7 x3 + 7 x 4 = 1.
Поскольку определитель при неизвестных x1 и x2 отличен от нуля, то их
можно принять в качестве главных и переписать систему в виде:
⎧ x1 + x 2 = 1 + 2 x3 + x 4 − x5 ,
⎨
⎩− 4 x 2 = −7 x3 − 7 x 4 + 1,
откуда x2 = 7/4 x3 + 7/4 x4 -1/4, x1 = 1/4 x3 -3/4 x4 - x5 + 5/4 - общее решение
системы, имеющей бесчисленное множество решений. Придавая свободным
неизвестным x3, x4, x5 конкретные числовые значения, будем получать
частные решения. Например, при x3 = x4 = x5 = 0 x1= 5/4, x2 = - 1/4. Вектор
C(5/4, - 1/4, 0, 0, 0) является частным решением данной системы.
Пример 11. Исследовать систему уравнений и найти общее решение в
зависимости от значения параметра а.
⎧2 x1 − x 2 + x3 + x 4 = 1,
⎪
⎨ x1 + 2 x 2 − x3 + 4 x 4 = 2 ,
⎪ x + 7 x − 4 x + 11x = a .
2
3
4
⎩ 1
Решение. Данной системе соответствует матрица⎯
41
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
⎛ 2 -1 1 1 1 ⎞
⎟
⎜
( A / B ) = ⎜1 2 - 1 4 2 ⎟ .
⎜1 7 - 4 11 a ⎟
⎠
⎝
Имеем ⎯
2 ⎞
⎛1
⎛1 2 - 1 4
⎜
⎟
( A / B ) ≈ ⎜ 0 5 - 3 7 a - 2 ⎟ ~ ⎜⎜ 0
⎜
⎜0 - 5 3 - 7 - 3 ⎟
⎝0
⎝
⎠
2 -1
5 -3
0 0
4
7
0
2 ⎞
⎟
a - 2⎟ ,
⎟
a - 5⎠
следовательно, исходная система равносильна такой системе:
⎧ x1 + 2 x 2 − x3 + 4 x 4 = 2 ,
⎪
⎨ 5 x 2 − 3 x3 + 7 x 4 = a − 2 ,
⎪
0 = a − 5.
⎩
Отсюда видно, что система совместна только при a=5. Общее решение
в этом случае имеет вид:
x2 = 3/5 + 3/5x3 - 7/5x4, x1 = 4/5 - 1/5x3 - 6/5x4.
Пример 12. Выяснить, будет ли линейно зависимой система векторов:
a1 = (1, 1, 4, 2),
a2 = (1, -1, -2, 4),
a3 = (0, 2, 6, -2),
a4 = (-3, -1, 3, 4),
a5 = (-1, 0, - 4, -7).
Решение. Система векторов является линейно зависимой, если
найдутся такие числа x1, x2, x3, x4, x5, из которых хотя бы одно отлично от
нуля, что выполняется векторное равенство:
x1 a1 + x2 a2 + x3 a3 + x4 a4 + x5 a5 = 0.
В координатной записи оно равносильно системе уравнений:
− 3 x 4 − x5 = 0 ,
⎧ x1 + x 2
⎪x − x + 2x − x
= 4,
⎪ 1
2
3
4
⎨
⎪4 x1 − 2 x 2 + 6 x3 + 3 x 4 − 4 x5 = 0,
⎪⎩2 x1 + 4 x 2 − 2 x3 + 4 x 4 − 7 x5 = 0 .
Итак, получили систему линейных однородных уравнений. Решаем ее
методом исключения неизвестных:
42
© КРСУ Усенов И.А., Усенова Р.К.
1
⎛1
⎜
⎜1 - 1
⎜4 - 2
⎜
4
⎝2
⎛1 1
⎜
0 -2
∼ ⎜⎜
0
0
⎜
0
⎝0
0
2
-3
-1
6
3
-2
4
0
2
0
0
Кафедра высшей математики
- 1⎞ ⎛ 1 1
0 - 3 - 1⎞
⎟ ⎜
⎟
0⎟ ⎜ 0 - 2
2
2 1⎟
∼
∼
- 4 ⎟ ⎜ 0 - 10 10 - 5 10⎟
⎟ ⎜
⎟
- 7⎠ ⎝ 0
2 - 2 10 - 5⎠
- 3 - 1⎞
⎟
2 1⎟
∼
- 3 1⎟
⎟
12 - 4⎠
⎛1 1
⎜
⎜0 - 2
⎜0
0
⎜
0
⎝0
0
2
0
0
0
⎛1 1
⎜
2
⎜0 - 2
⎜0 - 2
2
⎜
2 -2
⎝0
- 3 - 1⎞
⎟
2 1⎟
∼
- 3 1⎟
⎟
0 0⎠
⎛1 1 - 3
⎜
2
⎜0 - 2
⎜0
0 -3
⎜
0
0
⎝0
- 3 -1⎞
⎟
2 1⎟
∼
-1
2⎟
⎟
10 - 5⎠
0 - 1⎞
⎟
2 1⎟
.
0 1⎟
⎟
0 0⎠
Система приведена к ступенчатому виду, ранг матрицы равен 3, значит,
однородная система уравнений имеет решения, отличные от нулевого (r < n).
Определитель при неизвестных x1, x2, x4 отличен от нуля, поэтому их можно
выбрать в качестве главных и переписать систему в виде:
⎧ x1 + x 2 − 3x 4 = x5 ,
⎪
⎨− 2 x 2 + 2 x 4 = −2 x3 − x5 ,
⎪− 3 x = − x .
4
5
⎩
Имеем:
x4 = 1/3 x5, x2 = 5/6x5+x3, x1 = 7/6 x5 -x3.
Система имеет бесчисленное множество решений; если свободные
неизвестные x3 и x5 не равны нулю одновременно, то и главные неизвестные
отличны от нуля. Следовательно, векторное уравнение
x1 a1 + x2 a2 + x3 a3 + x4 a4 + x5 a5 = 0
имеет коэффициенты, не равные нулю одновременно; пусть например, x5 = 6,
x3 = 1. Тогда x4=2, x2 = 6, x1=6 и мы получим соотношение
6a1 + 6a2 + a3 + 2a4 + 6a5 = 0,
т.е. данная система векторов линейно независима.
43
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
§ 2.2 Использование систем линейных уравнений при решении
экономических задач
Пример 13. Из некоторого листового материала необходимо выкроить
360 заготовок типа А, 300 заготовок типа Б и 675 заготовок типа В. При этом
можно применять три способа раскроя. Количество заготовок, получаемых из
каждого листа при каждом способе раскроя, указано в таблице:
Тип
Способ раскроя
заготовки
1
2
3
А
3
2
1
Б
1
6
2
В
4
1
5
Записать в математической форме условия выполнения задания.
Решение. Обозначим через x, y, z количество листов материала,
раскраиваемых соответственно первым, вторым и третьим способами. Тогда
при первом способе раскроя x листов будет получено 3x заготовок типа А,
при втором - 2y, при третьем - z.
заготовкам типа А сумма
Для полного выполнения задания по
3x +2y +z должна равняться 360, т.е.
3x +2y + z =360.
Аналогично получаем уравнения
x + 6y +2z = 300
4x + y + 5z = 675,
которым должны удовлетворять неизвестные x, y, z для того, чтобы
выполнить задание по заготовкам Б и В. Полученная система линейных
уравнений
⎧3 x + 2 y + z = 360 ,
⎪
⎨ x + 6 y + 2 z = 300 ,
⎪4 x + y + 5 z = 675 .
⎩
44
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
и выражает в математической форме условия выполнения всего задания по
заготовкам А, Б и В. Решим систему методом исключения неизвестных.
Запишем расширенную матрицу системы и приведем ее с помощью
элементарных преобразований к треугольному виду.
⎛1
⎜
⎜3
⎜
⎝4
⎛1
⎜
∼ ⎜0
⎜
⎝0
6
2
1
2
1
5
6 2
16 5
2 9
300 ⎞
⎟
360⎟ ∼
⎟
675⎠
⎛ 1 6 2 300 ⎞
⎜
⎟
⎜ 0 - 16 - 5 - 540⎟ ∼
⎜
⎟
⎝ 0 - 7 2 15 ⎠
300 ⎞
⎛1 6
⎟
⎜
540⎟ ∼ ⎜ 0 2
⎟
⎜
570⎠
⎝ 0 16
2
9
5
⎛ 1 6 2 300 ⎞
⎜
⎟
⎜ 0 16 5 540⎟ ∼
⎜
⎟
⎝ 0 - 14 4 30 ⎠
300 ⎞
⎟
570⎟ ∼
⎟
540⎠
⎛1
⎜
⎜0
⎜
⎝0
6 2 300 ⎞
⎟
2 9 570 ⎟ .
⎟
0 - 67 - 4020⎠
Следовательно, исходная система равносильна следующей системе:
⎧ x + 6 y + 2 z = 300 ,
⎪
⎨ 2 y + 9 z = 570 ,
⎪
− 67 z = −4020 .
⎩
Из последнего уравнения находим z = 60; подставляя найденное значение z
во второе уравнение, получим y = 15 и, наконец, из первого имеем x = 90.
Итак, вектор C (90, 15, 60) есть решение системы.
Пример 14. Три судна доставили в порт 6000 т чугуна, 4000 т железной
руды и 3000 т апатитов. Разгрузку можно производить как непосредственно в
железнодорожные вагоны для последующей доставки потребителям, так и на
портовые склады. В вагоны можно разгрузить 8000 т, а остаток груза
придется направить на склады. Необходимо учесть, что поданные в порт
вагоны не приспособлены для перевозки апатитов. Стоимость выгрузки 1 т в
вагоны составляет соответственно 4,30, 5,25 и 2,20 ден. ед.
Записать в математической форме условия полной разгрузки судов,
если затраты на нее должны составить 58850 ден. ед.
Решение. По условию задачи, доставленные в порт чугун, железную
руду
и
апатиты
можно
разгрузить
двумя
способами:
либо
железнодорожные вагоны, либо в портовые склады. Обозначим через x
в
i j
количество груза (в тоннах) i-го вида (i= 1,2,3), которое предполагается
45
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
разгрузить j-м способом (j = 1, 2). Таким образом, задача содержит шесть
неизвестных. Условие полной разгрузки чугуна можно записать в виде
x11 + x12 = 6000,
где x11, x12 - части чугуна, разгружаемого соответственно в вагоны и на
склады. Аналогичное условие должно выполняться и для железной руды:
x2 1 + x22 = 4000.
Что же касается апатитов, то их можно разгружать только на склады, а
поэтому неизвестное x31 = 0, и условие полной разгрузки апатитов принимает
вид
x32 =3000.
Условие полной загрузки всех поданных в порт вагонов запишется так:
x11 + x21 = 8000.
Затраты на разгрузку, по условию, определены в 58850 ден. ед., что
можно выразить записью:
4,3x 11 + 7,8 x 12 + 5,25 x 21 + 6,4x 22 + 3,25x 32 = 58850.
Итак, с учетом сложившейся в порту ситуации условия полной
разгрузки судов выражаются в математической форме системой линейных
уравнений (1) - (2). С учетом (3) уравнение (5) перепишется в виде:
4,3x 11 + 7,8x 12 +5,25x 21 +6,4x 22 = 49100,
и теперь мы имеем систему из четырех уравнений с четырьмя неизвестными
x 11, x 12, x 21, x 22, расширенная матрица которой имеет вид:
1
⎛1
⎜
0
0
⎯ ( A / B ) = ⎜⎜
1
0
⎜
⎜ 4,3 7,8
⎝
0
0
1
1
1
0
5,25 6,4
6000 ⎞
⎟
4000 ⎟
.
8000 ⎟
⎟
49100 ⎟⎠
Преобразуем ее к треугольному виду:
46
© КРСУ Усенов И.А., Усенова Р.К.
1
0
0
⎛1
⎜
1
0
1
0
( A / B ) ∼ ⎜⎜
0
0
1
1
⎜
⎝ 4,3 7,8 5,25 6,4
⎛1 1
⎜
0 -1
∼ ⎜⎜
0 0
⎜
⎝0 0
Кафедра высшей математики
6000 ⎞
0
⎛1 1
⎜
⎟
8000 ⎟
1
⎜0 - 1
∼
⎜0 0
4000⎟
1
⎜
⎟
49100⎠
⎝ 0 3,5 5,25
6000 ⎞
⎛1 1
⎜
⎟
1
0
2000 ⎟
⎜0 - 1
∼
⎜0
1
1
4000⎟
0
⎜
⎟
8,75 6,4 30300⎠
0
⎝0
0
6000 ⎞
⎟
2000 ⎟
∼
1
4000 ⎟
⎟
6,4 23300⎠
0
0
6000 ⎞
⎟
1
0
2000 ⎟
.
1
1
4000 ⎟
⎟
0 - 2,35 - 4700⎠
0
0
0
Наша система равносильна следующей системе:
⎧ x11 + x12 = 6000 ,
⎪
− x12 + x 21 = 2000 ,
⎪
⎨
x 21 + x 22 = 4000
⎪
⎪⎩
− 2,35 x 22 = . − 4700,
откуда x 22 = 2000, x 21 = 2000, x 12 = 0, x 11 = 6000.
Пример 15. На предприятии имеется четыре технологических способа
изготовления изделий А и Б из некоторого сырья. В таблице указано
количество изделий, которое может быть произведено из единицы сырья
каждым из технологических способов.
Записать в математической форме условия выбора технологий при
производстве из 94 ед. сырья 574 изделий А и 328 изделий Б.
Изделие
Выход из единицы сырья
I
II
III
IV
А
2
1
7
4
Б
6
12
2
3
Решение. Обозначим через x1, x2, x3, x4 количество сырья, которое
следует переработать по каждой технологии, чтобы выполнить плановое
задание.
Получим
систему
трех
линейных
уравнений
с
четырьмя
неизвестными:
47
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
⎧ x1 + x 2 + x3 + x 4 = 94,
⎪
⎨2 x1 + x 2 + 7 x3 + 4 x 4 = 574 ,
⎪6 x + 12 x + 2 x + 3x = 328 .
2
3
4
⎩ 1
Решаем ее методом Гаусса:
⎛1 1
⎜
⎜2 1
⎜
⎝ 6 12
1 94 ⎞ ⎛ 1 1 1 1 94 ⎞
⎛1 1 1
⎟ ⎜
⎟
⎜
4 574 ⎟ ∼ ⎜ 0 - 1 5 2 386 ⎟ ∼ ⎜ 0 - 1 5
⎟ ⎜
⎟
⎜
3 328 ⎠ ⎝ 0 6 - 4 - 3 - 236⎠
⎝ 0 0 26
1
7
2
1 94 ⎞
⎟
2 386⎟ .
⎟
9 2080⎠
Имеем: r (А) = r (А) = 3, следовательно, число главных неизвестных
равно трем, одно неизвестное x4 - свободное. Исходная система равносильна
следующей:
⎧ x1 + x 2 + x3 = 94 − x 4 ,
⎪
⎨ − x 2 + 5 x3 = 386 − 2 x 4 ,
⎪
26 x3 = 2080 − 9 x 4 .
⎩
Из последнего уравнения находим x3 = 80 - 9/26 x4, подставляя x3 во второе
уравнение, будем иметь: x2 = 14 + 7/26x4 и, наконец, из первого уравнения
получим: x1 = - 12/13 x4. С математической точки зрения система имеет
бесчисленное множество решений, т. е. неопределенна. С учетом реального
экономического
содержания
величины
x1
и
x4
не
могут
быть
отрицательными, тогда из соотношения x1 = - 12/13 x4 получим: x1 = x4 = 0.
Тогда вектор (0, 14, 80, 0) является решением данной системы.
Математическая
модель
межотраслевого
баланса.
Модель
межотраслевого баланса, разработанная профессором В. Леонтьевым
(Гарвардский университет, США), имеет вид:
n
∑a
ij
x j + yi = x i
(i = 1, n) ,
(10)
j=1
или, в матричной форме,
AX + Y = X,
(11)
где А = (a i j) - матрица коэффициентов прямых затрат, Х - вектор валовых
выпусков, Y - вектор конечного продукта.
Перепишем систему (11) в виде
48
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
(E - A) X = Y,
(12)
где E - единичная матрица n-го порядка, тогда решение системы (12)
относительно неизвестных значений объемов производства продукции при
заданном векторе конечного продукта находится по формуле
X = (E - A) −1 Y.
(13)
Здесь (E - A) −1 - матрица коэффициентов полных затрат. Элемент b i j
матрицы (E - A) −1 характеризует потребность в валовом выпуске отрасли i,
который необходим для получения в процессе материального производства
единицы
конечного
продукта
отрасли
j.
Благодаря
возможность рассматривать валовые выпуски x
i
этому
имеется
в виде функций
планируемых значений y j конечных продуктов отраслей:
n
xi = ∑ bi jy j
.
j=1
Пример 16. Пусть дана Леонтьевская балансовая модель “затраты выпуск” X = AX +Y. Найти вектор конечной продукции Y при заданном X,
где
⎛ 100⎞
⎛ 0,1 0 0,6⎞
⎟
⎜
⎟
⎜
0 0,7 0,2⎟ ; X = ⎜ 200⎟
⎜
A= ⎜
⎟;
⎜
⎟
⎝ 150⎠
⎝ 0,7 0,1 0,1⎠
Решение. Имеем: Y = (E - A) X, где E - единичная матрица третьего
порядка.
⎛ 0,9 0 - 0,6 ⎞
⎜
⎟
E - A = ⎜ 0 0,3 - 0,2 ⎟ ,
⎜ - 0,7 - 0,1 0,9 ⎟
⎝
⎠
0 - 0,6 ⎞ ⎛ 100⎞ ⎛ 0,9 ⋅ 100 + 0 − 0,6 ⋅ 150⎞ ⎛ 0 ⎞
⎛ 0,9
⎟ ⎜
⎜
⎟ ⎜ ⎟
⎟ ⎜
0,3 - 0,2 ⎟ ⎜ 200⎟ = ⎜ 0 + 0,3 ⋅ 200 - 0,2 ⋅ 150 ⎟ = ⎜ 30⎟ .
значит, Y= ⎜ 0
⎟ ⎜
⎜
⎟ ⎜ ⎟
⎟ ⎜
⎝ -0,7 - 0,1 0,9⎠ ⎝ 150⎠ ⎝ -0,7 ⋅ 100 - 0,1 ⋅ 200 + 0,9 ⋅ 150 ⎠ ⎝ 45⎠
49
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
Пример 17. Пусть дана Леонтьевская балансовая модель “затратывыпуск”. Определить, будет ли продуктивной матрица технологических
коэффициентов A. Найти вектор валовой продукции X при заданном Y, где
⎛ 0,125 0,125 ⎞
⎛ 300 ⎞
⎟⎟ ; Y = ⎜⎜
⎟⎟ .
⎝1,125 0,125 ⎠
⎝ 400 ⎠
A= ⎜⎜
Решение. Для решения вопроса о продуктивности матрицы A следует
найти собственные значения этой матрицы. Составим характеристическое
уравнение:
0,125 - λ
0,125
0,125 - λ
1,125
=0,
Или (0,125 -λ)2 - 0,140625 = 0 ⇒ 0,125 - λ = ± 0,375.
Следовательно, λ1 = 0,5; λ2 = - 0,25. Оба корня по модулю меньше
единицы, значит, матрица технологических коэффициентов A продуктивная.
Для
определения
вектора
валовой
продукции
X
имеем
формулу
X = (E - A) −1 Y.
Найдем обратную матрицу для матрицы
⎛ 0,875 - 0,125⎞
⎟.
⎝ -1,125 0,875 ⎠
E - A= ⎜
Обозначим B = E-A, тогда B -1 =
1 ⎛ A 11 A 21 ⎞
1 ⎛ 0,875
⎟=
⎜
⎜
det B ⎝ A 12 A 22 ⎠ 5,4 ⎝ 1,125
0,125⎞
⎟.
0,875⎠
Следовательно,
X =
1 ⎛ 0,875 0,125⎞ ⎛ 300⎞
1 ⎛ 0,875 ⋅ 300 + 0,125 ⋅ 400⎞
1 ⎛ 312,5⎞ ⎛ 57,9 ⎞
⎟ ⋅⎜
⎜
⎟.
⎜
⎟ =⎜
⎜
⎟=
⎟=
5,4 ⎝ 1,125 0,875⎠ ⎝ 400⎠ 5,4 ⎝ 1125
,
⋅ 300 + 0,125 ⋅ 400 ⎠ 5,4 ⎝ 687,5⎠ ⎝ 127,3⎠
50
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
Глава 3. ЭЛЕМЕНТЫ МАТРИЧНОГО АНАЛИЗА
§ 3.1 Векторные (Линейные) пространства
1. Понятие векторного пространства.
Определение. Множество L элементов любой природы называется
векторным пространством, если выполнены следующие три требования:
I.
Имеется правило, посредством которого любым двум элементам x, y ∈ L
ставится в соответствие третий элемент этого множества, называемый
суммой элементов x и y и обозначаемый символом x + y .
II.
Имеется правило, посредством которого любому элементу x ∈ L и
любому числу λ
ставится в соответствие элемент этого множества,
называемый произведением элемента x на число λ и обозначаемый
символом λx .
III.
Для любых элементов x, y, z ∈ L и любых чисел λ и µ выполнены
следующие аксиомы:
1)
x+ y = y+ x;
2)
( x + y) + z = x + ( y + z) ;
3)
∃ 0 ∈ L : x + 0 = x ∀x ∈ L . Этот элемент 0 пространства L называют
нулевым (не путать с числом 0 !);
4)
∀ x ∈ L ∃ (− x) ∈ L : x + (− x) = 0 .
Такой
элемент
(− x)
называют
противоположным для x ;
5)
λ ( x + y ) = λx + λy ;
6)
(λ + µ ) x = λx + µx ;
7)
λ ( µx) = (λµ ) x ;
8)
1⋅ x = x .
Замечание 1. Если в пункте II мы ограничиваемся вещественными
числами, то L называется вещественным векторным пространством; если же
определено умножение на любое комплексное число, то векторное
пространство называется комплексным.
51
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
Замечание 2. Элементы произвольного векторного пространства
принято называть векторами. То обстоятельство, что часто термин «вектор»
употребляется
в
более
узком
смысле,
при
этом
не
приводит
к
недоразумениям, а, напротив, обращаясь к сложившимся геометрическим
представлениям, можно уяснить, а зачастую и предвидеть ряд результатов,
справедливых для векторных пространств произвольной природы.
При введении понятия векторного пространства мы абстрагируемся не
только от природы изучаемых объектов, но и от конкретного вида правил
образования суммы элементов и произведения элемента на число. Важно
лишь, чтобы эти правила удовлетворяли аксиомам векторного пространства.
Примеры векторных пространств
1. Пространство V2 геометрических векторов на плоскости.
2. Пространство V3 геометрических векторов в трехмерном пространстве.
3. Множество
Pn (R)
многочленов с вещественными коэффициентами
степени не выше n образует векторное пространство относительно операций
сложения многочленов и умножения многочлена на вещественное число.
4. Множество M m×n (R) матриц одинаковых размеров образуют векторное
пространство относительно операций сложения матриц и умножения
матрицы на число. В частности, часто встречается и используется векторное
пространство матриц-строк M 1×n ( R) ( M 1×n (C ) ). Для него принято другое
обозначение – R n ( C n ). Элементами этого векторного пространства служат
упорядоченные совокупности n произвольных вещественных (комплексных)
чисел ( x1 , x2 ,..., xn ) .
5. (Нестандартный пример). Рассмотрим множество R+ всех положительных
вещественных чисел. Определим «сумму» двух элементов x, y ∈ R+ как
произведение вещественных чисел x и y (понимаемое в обычном смысле):
x + y 6 xy . «Произведение» элемента
x ∈ R+
на вещественное число λ
определим как возведение числа x в степень λ : λx 6 x λ . Нулевым элементом
пространства будет служить вещественное число 1, а противоположным
52
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
элементом (для данного элемента x ) будет число 1 x . Проверьте выполнение
аксиом векторного пространства (которые в обычной записи принимают
другой вид: вместо λ ( µx) = (λµ ) x мы имеем (x µ ) = x λµ и т.д.). В этом примере,
λ
быть может, для обозначения суммы элементов пространства и для
произведения элемента пространства на число предпочтительнее другие
обозначения (например, ⊕ и ⊗ ).
2. Некоторые свойства произвольных векторных пространств. Из
определения
векторного
пространства
следует
ряд
утверждений,
справедливых для произвольных векторных пространств.
1) В векторном пространстве существует единственный нулевой
элемент.
2) Для каждого элемента векторного пространства существует
единственный противоположный элемент.
3) 0 ⋅ x = 0
∀ x∈ L.
4) Для любого элемента x ∈ L противоположный ему элемент (− x)
равен произведению элемента x на число (−1) , т.е. (−1) x = − x .
Отметим также, что из определения векторного пространства следует
существование
и
единственность
разности
любых
двух
элементов
векторного пространства x и y , которая определяется как элемент z ,
удовлетворяющий условию
x= z+ y.
Этим элементом служит сумма
z = x + ( −1) y .
53
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
§ 3.2 Размерность и базис векторного пространства
1. Линейно зависимые и независимые системы векторов. Пусть L –
произвольное вещественное векторное пространство, {xi } = {x1 , x2 ,..., xn } ⊂ L ,
{α i } = {α1 , α 2 ,..., α n } ⊂ R . Линейной комбинацией векторов {xi } с коэффициентами
{α i } называется вектор x ∈ L , получаемый по правилу
n
x = α1 x1 + α 2 x2 + ... + α n xn = ∑ α i xi .
(14)
i =1
Линейная комбинация α1 x1 + α 2 x2 + ... + α n xn называется нетривиальной,
если в ней, хотя бы один из коэффициентов α i отличен от нуля. Линейная
комбинация вида 0 ⋅ x1 + 0 ⋅ x2 + ... + 0 ⋅ xn называется тривиальной; она, очевидно,
равна нулевому вектору 0 .
Определение. Система векторов {xi } называется линейно зависимой,
если существует, хотя бы одна нетривиальная линейная комбинация этих
векторов, равная нулевому вектору. В противном случае, т.е. если только
тривиальная линейная комбинация данных векторов равна нулевому вектору,
векторы называются линейно независимыми.
Теорема (критерий линейной зависимости). Для того чтобы система
векторов
{xi }
линейного
пространства
была
линейно
зависимой,
необходимо и достаточно, чтобы, по крайней мере, один из этих векторов
являлся линейной комбинацией остальных.
Имеет место два простых, но важных предложения о линейной
зависимости.
1) Если среди векторов x1 , x2 ,..., xn имеется хотя бы один нулевой
вектор, то вся система векторов линейно зависима.
2) Если среди векторов x1 , x2 ,..., xn некоторые образуют линейно
зависимую систему, то и вся система x1 , x2 ,..., xn линейно зависима.
Примером
линейно
независимых
векторов
является
два
неколлинеарных, вектора a1 , a2 на плоскости. Действительно, условие
54
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
λ1a1 + λ2 a 2 = 0 будет выполняться лишь в случае, когда λ1 = λ2 = 0 , ибо если,
например λ2 ≠ 0 , то a 2 = −
Пример
1.
λ1
a1 и векторы a1 , a 2 коллинеарны.
λ2
Показать,
что
векторы
a1 = (1, 3,1, 3) ,
a 2 = (2,1,1, 2) ,
a3 = (3, − 1,1,1) линейно зависимы.
Решение. Составим векторное равенство λ1 a1 + λ2 a 2 + λ3 a3 = 0 . Записывая
a1 , a 2 , a3 в виде вектор-столбцов, получим
⎛1⎞
⎛ 2⎞
⎛ 3 ⎞ ⎛0⎞
⎜ ⎟
⎜ ⎟
⎜ ⎟ ⎜ ⎟
⎜ 3⎟
⎜1⎟
⎜ − 1⎟ ⎜ 0 ⎟
λ1 ⎜ ⎟ + λ 2 ⎜ ⎟ + λ3 ⎜ ⎟ = ⎜ ⎟.
1
1
1
0
⎜ ⎟
⎜ ⎟
⎜ ⎟ ⎜ ⎟
⎜ 3⎟
⎜ 2⎟
⎜ 1 ⎟ ⎜0⎟
⎝ ⎠
⎝ ⎠
⎝ ⎠ ⎝ ⎠
Задача свелась, таким образом, к решению системы
⎧λ1 + 2λ 2 + 3λ3 = 0,
⎪ 3λ + λ − λ = 0,
⎪ 1
2
3
⎨
⎪ λ1 + λ2 + λ3 = 0,
⎪⎩3λ1 + 2λ 2 + λ3 = 0.
Решив систему методом Гаусса, приведем ее к виду
⎧ λ1 + 2λ 2 + 3λ3 = 0,
⎪
λ2 + 2λ3 = 0,
⎪
.
⎨
0
0
,
=
⎪
⎪⎩
0 = 0.
Таким образом, данная система имеет бесконечное множество решений
λ1 = c, λ2 = −2c, λ3 = c , где с- произвольная константа.
Итак, для данных векторов a1 , a 2 , a3 условие λ1 a1 + λ2 a 2 + λ3 a3 = 0
выполняется не только при λ1 = λ2 = λ3 = 0 , а например при
λ1 = 1, λ2 = −2, λ3 = 1 ,
следовательно, эти векторы линейно зависимы.
2. Базис и размерность. Фундаментальным вопросом теории
векторных пространств является вопрос о том, можно ли, а если можно, то
как, произвольный вектор пространства представить в виде линейной
55
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
комбинации фиксированного набора векторов из этого пространства. Далее
мы получим ответ на этот вопрос.
Определение. Система линейно независимых векторов {e1 , e2 ,..., en }
векторного пространства L называется базисом этого пространства, если
любой вектор из L может быть представлен в виде линейной комбинации
векторов этой системы, т.е. для каждого вектора
x∈L
существуют
вещественные числа α1 , α 2 ,..., α n такие, что имеет место равенство
x = α 1e1 + α 2 e2 + ... + α n en .
Это равенство называется разложением вектора x по базису e1 , e2 ,..., en , а
числа α1 , α 2 ,..., α n называются координатами вектора x относительно базиса
(или в базисе) e1 , e2 ,..., en .
Очевидно, что нулевой вектор имеет все нулевые координаты, а вектор,
противоположный данному,- противоположные по знаку координаты.
Справедливы следующие теоремы:
Теорема (о единственности разложения по базису). Каждый вектор x
пространства L может быть разложен по базису e1 , e2 ,..., en единственным
образом, т.е. координаты каждого вектора
x
в базисе
e1 , e2 ,..., en
определяются однозначно.
Теорема. Если e1 , e2 ,..., en - система линейно независимых векторов
пространства L и любой вектор а линейно выражается через e1 , e2 ,..., en , то
пространства L является n-мерным, а векторы e1 , e2 ,..., en - его базисом.
Главное значение базиса заключается в том, что операции сложения
векторов и умножения их на числа при задании базиса превращаются в
соответствующие операции над числами – координатами этих векторов. А
именно, справедлива следующая
Теорема. При сложении двух любых векторов линейного пространства
L
их
координаты
(относительно
любого
базиса
пространства)
складываются; при умножении произвольного вектора на любое число λ все
координаты этого вектора умножаются на λ .
56
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
n
Пример 2. . Исследуем вопрос о базисе пространства R , введенного
ранее при рассмотрении примеров векторных пространств. Покажем, что n
элементов
указанного
e1 = (1,0,0,...,0,0), e2 = (0,1,0,...,0,0), ... , en = (0,0,0,...,0,1)
пространства образуют базис.
Решение. Во-первых, эти векторы линейно независимы. Проверка
линейной независимости набора {ei } состоит в определении значений
α1 , α 2 ,..., α n , при которых возможно равенство
n
∑α e
i =1
i i
= 0 = (0,0,...,0) .
Но в силу преведущий теоремы
n
∑α e
i =1
i i
= (α1 , α 2 ,..., α n ) ,
(15)
а последний вектор является нулевым лишь при условии α1 = α 2 = ... = α n = 0 .
Во-вторых, всякий вектор y = ( β1 , β 2 ,..., β n ) заведомо представим в виде
линейной комбинации векторов {ei } : y = β1e1 + β 2e2 + ... + β n en и, значит, набор
{ei } образует базис.
Определение. Векторное пространство L называется n -мерным, если в
нем существуют n линейно независимых векторов, а любые n + 1 векторов
уже являются линейно зависимыми. При этом число
n
называется
размерностью пространства L .
Размерность векторного пространства, состоящего из одного нулевого
вектора, принимается равной нулю.
Размерность пространства L обычно обозначают символом dim L .
Определение.
бесконечномерным,
Векторное
если
в
нем
пространство
существует
любое
L
называется
число
линейно
независимых векторов. В этом случае пишут dim L = ∞ .
Выясним связь между понятиями базиса и размерности пространства.
Теорема. Если L – векторное пространство размерности n , то любые
n линейно независимых векторов этого пространства образуют его базис.
57
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
Теорема. Если векторное пространство L имеет базис, состоящий из
n векторов, то dim L = n .
Пример 3. В базисе ee , e2 , e3 заданы векторы a1 = (1,1, 0) , a 2 = (1, − 1,1) ,
a3 = (−3, 5, − 6) . Показать, что векторы a e , a 2 , a3 образуют базис.
Решение. Составим векторное равенство λ1 a1 + λ2 a 2 + λ3 a3 = 0 . Записывая
a1 , a 2 , a3 в виде вектор-столбцов, получим
⎛ − 3⎞ ⎛ 0⎞
⎛1⎞
⎛1⎞
⎜ ⎟ ⎜ ⎟
⎜ ⎟
⎜ ⎟
λ1 ⎜ 1 ⎟ + λ2 ⎜ − 1⎟ + λ3 ⎜ 5 ⎟ = ⎜ 0 ⎟.
⎜ − 6⎟ ⎜ 0⎟
⎜1⎟
⎜0⎟
⎝ ⎠ ⎝ ⎠
⎝ ⎠
⎝ ⎠
Задача свелась, таким образом, к решению системы
⎧λ1 + λ2 − 3λ3 = 0,
⎪
⎨λ1 − λ2 + 5λ3 = 0,
⎪
λ2 − 6λ3 = 0.
⎩
Решив систему методом Гаусса, получим
⎧λ1 + 2λ 2 + 3λ3 = 0 ⎧λ1 + 2λ 2 + 3λ3 = 0 ⎧ λ1 = 0,
⎪
⎪
⎪
⎨ − 2λ 2 + 8λ3 = 0 ⇒⎨ − 2λ 2 + 8λ3 = 0 ⇒ ⎨λ 2 = 0,
⎪ λ − 6λ = 0
⎪ λ = 0.
⎪
− 2λ 3 = 0
2
3
⎩
⎩ 3
⎩
Таким образом, векторы a1 , a 2 , a3 составляют систему линейно независимых
векторов и, следовательно, образуют базис. Система трех линейно
независимых векторов ae , a 2 , a3 образуют базис то, размерность dim L = 3 .
3. Переход к новому базису. Пусть в пространстве R имеются два
базиса: старый e1 , e2 ,..., en и новый e€1 , e€2 ,..., e€n . Каждый из векторов нового
базиса можно выразить в виде линейной комбинации векторов старого
базиса:
⎧ e€1 = a11e1 + a12 e2 + ... + a1n en ,
⎪e€ = a e + a e + ... + a e ,
⎪ 2
21 1
22 2
21n n
.
⎨
⎪ ... ... ... ... ... ... ... ... ... ...
⎪⎩ e€n = a n1e1 + a n 2 e2 + ... + a nn en .
(16)
Полученная система означает, что переход от старого базиса e1 , e2 ,..., en к
новому e€1 , e€2 ,..., e€n задается матрицей перехода
58
© КРСУ Усенов И.А., Усенова Р.К.
⎛ a11
⎜
⎜a
A = ⎜ 21
...
⎜
⎜a
⎝ n1
Кафедра высшей математики
a12
a 22
...
...
...
an2
...
...
a1n ⎞
⎟
a2n ⎟
,
... ⎟
⎟
a nn ⎟⎠
причем коэффициенты размножения новых базисных векторов по старому
базису образуют столбцы этой матрицы.
Матрица А невырожденная, так, как столбцы базисные векторы.
Обратный переход от нового базиса e€1 , e€2 ,..., e€n к старому базису e1 , e2 ,..., en
осуществляется с помощью обратной матрицы A −1 .
Найдем зависимость между координатами вектора в разных
базисах. Пусть рассматриваемый вектор х имеет координаты (x1 , x2 ,..., x n )
относительно старого базиса и координаты (x€1 , x€2 ,..., x€n ) относительно нового
базиса т.е.
x = x€1e€1 + x€2 e€2 + ... + x€n e€n = x1e1 + x 2 e2 + ... + x n en . (17)
Подставляя значения e€1 , e€2 ,..., e€n из системы (16) в левую часть равенства
(17), получим после преобразований:
⎧ x1 = a11 x€1 + a12 x€2 + ... + a1n x€n ,
⎪ x = a x€ + a x€ + ... + a x€ ,
⎪ 2
21 1
22 2
21n n
,
⎨
⎪ ... ... ... ... ... ... ... ... ... ...
⎪⎩ x n = a n1 x€1 + a n 2 x€2 + ... + a nn x€n .
или в матричной форме
⎛ x1 ⎞
⎜ ⎟
⎜ x2 ⎟
⎜ . ⎟
⎜ ⎟=
⎜ . ⎟
⎜ . ⎟
⎜x ⎟
⎝ n⎠
⎛ x€1 ⎞
⎜ ⎟
⎜ x€2 ⎟
⎜ . ⎟
A⎜ ⎟
⎜ . ⎟
⎜ . ⎟
⎜ x€ ⎟
⎝ n⎠
или
⎛ x€1 ⎞
⎛ x1 ⎞
⎜ ⎟
⎜ ⎟
⎜ x€2 ⎟
⎜ x2 ⎟
⎜ . ⎟
⎜ . ⎟
−1
⎜ ⎟ = A ⎜ ⎟.
⎜ . ⎟
⎜ . ⎟
⎜ . ⎟
⎜ . ⎟
⎜ x€ ⎟
⎜x ⎟
⎝ n⎠
⎝ n⎠
(18)
Пример 4. Вектор b = (4, − 4, 5) , заданный в базисе e1 , e2 , e3 , выразить в
базисе a1 , a 2 , a3 , где a1 = (1,1, 0) , a 2 = (1, − 1,1) , a3 = (−3, 5, − 6) .
Решение. Выразим связь между базисами:
59
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
a1 = e1 + e2 ,
⎧
⎪
⎨ a 2 = e1 − e2 + e3 ,
⎪a = −3e + 5e − 6e .
1
2
3
⎩ 3
Матрица перехода от базиса e1 , e2 , e3 к базису a1 , a 2 , a3 имеет вид
⎛ 1 1 − 3⎞
⎜
⎟
A = ⎜1 −1 5 ⎟ .
⎜0 1 − 6⎟
⎝
⎠
Вычисляем A −1
2 ⎞
⎛1 3
⎟
1⎜
A = ⎜6 − 6 − 8⎟ .
4⎜
⎟
⎝1 −1 − 2⎠
−1
Тогда по формуле (18), имеем
2 ⎞⎛ 4 ⎞ ⎛ 0,5 ⎞
⎛1 3
⎛ x€1 ⎞
⎟
⎟⎜ ⎟ ⎜
⎜ ⎟ 1⎜
⎜ x€2 ⎟ = ⎜ 6 − 6 − 8 ⎟⎜ − 4 ⎟ = ⎜ 2 ⎟ ,
⎜ x€ ⎟ 4 ⎜ 1 − 1 − 2 ⎟⎜ 5 ⎟ ⎜ − 0,5 ⎟
⎠
⎠⎝ ⎠ ⎝
⎝
⎝ 3⎠
т.е. новые координаты вектора b в базисе a1 , a 2 , a3 есть 0,5; 2;-0,5 и вектор b
может быть представлен в виде:
b = 0,5a1 + 2a 2 − 0,5a3 .
60
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
§ 3.3 Евклидовы пространства
Выше мы определили линейное (векторное) пространства в котором
можно складывать векторы и умножить их на числа, ввели понятие
размерности и базиса, а теперь в данном пространстве введем метрику т.е.
способ измерять длины и углы. Это возможно если ввести понятие
скалярного произведения
1. Скалярное произведение.
Определение. Скалярным произведением двух векторов x = (x1 , x2 ,..., x n ) ,
y = ( y1 , y 2 ,..., y n ) называется числа
n
(x, y ) = x1 y1 + x2 y 2 + ... + xn y n = ∑ xi yi .
i =1
Скалярное произведение имеет экономический смысл. Если
x = ( x1 , x 2 ,..., x n ) есть вектор объемов различных товаров, а y = ( y1 , y 2 ,..., y n ) -
вектор их цен, то скалярное произведение (x, y ) выражает суммарную
стоимость этих товаров.
Скалярное произведение имеет следующие свойства:
1.
(x, y ) = ( y, x ) - коммутативное свойство;
2.
(x, y + z ) = ( y, x ) + (x, z ) - дистрибутивное свойство;
3.
(ax, y ) = a(x, y ) - для любого действительного числа а;
4.
(x, x ) > 0 , если х ненулевой вектор, (x, x ) = 0 если х нулевой вектор.
Мы рассмотрим векторные пространства произвольной природы, для
элементов которых каким-либо способом (причем, безразлично каким)
определено правило, ставящее в соответствие любым двум элементам число,
называемое скалярным произведением этих элементов. При этом важно
только, чтобы это правило обладало теми же четырьмя свойствами, что и
правило составления скалярного произведения для геометрических векторов
на плоскости или в пространстве. Векторные пространства, в которых
определено указанное правило, называются евклидовыми пространствами.
Дадим более строгое
61
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
Определение. Вещественное векторное пространство E называется
вещественным
евклидовым
пространством
(или
просто
евклидовым
пространством), если выполнены следующие два требования:
I.
Имеется правило, посредством которого любым двум элементам этого
пространства x и y ставится в соответствие вещественное число, называемое
скалярным произведением этих элементов (и обозначаемое символом ( x, y ) ).
II.
Указанное правило подчинено следующим четырем аксиомам:
1)
( x, y ) = ( y , x ) (коммутативность или симметрия);
2)
( x1 + x2 , y ) = ( x1 , y ) + ( x2 , y ) (дистрибутивность скалярного произведения
относительно сложения);
3)
( λ x, y ) = λ ( x, y ) ∀ λ ∈ R ;
4)
( x, x ) > 0 , если x ≠ 0 ; ( x, x ) = 0 , если x = 0 .
Пример
5.
Рассмотрим
векторные
пространства
V2
или
V3
геометрических векторов на плоскости или в пространстве. Скалярное
произведение любых двух векторов определим стандартным образом (как
произведение длин этих векторов на косинус угла между ними). В курсе
аналитической геометрии доказывается, что в этом случае выполняются
аксиомы 1) – 4). Следовательно, пространства V2 и V3 являются евклидовыми
пространствами.
Пример 6. Рассмотрим n -мерное координатное пространство R n . Если
x = ( x1 ,..., xn ) ∈ R n , y = ( y1 ,..., yn ) ∈ R n , то скалярное произведение ( x, y ) определим
def
равенством ( x, y ) = x1 y1 + ... + xn yn . В этом случае выполнение аксиом 1) – 4)
легко проверяется.
Длиной вектора х в евклидовом пространстве называется корень
квадратный из его скалярного квадрата
x =
( x, x ) =
x12 + x 22 + ... + x n2 .
Имеют место следующие свойства длины векторов:
1.
x = 0 тогда и только тогда, когда x = 0 ;
62
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
2.
λx = λ x , где λ -действительное числа;
3.
( x, y ) ≤
4.
x + y ≤ x + y - неравенство треугольника.
x y - неравенство Коши-Буняковского;
Угол ϕ между двумя векторами x и y определяется равенством
cos ϕ =
( x, y )
( x, y )
=
,
x ⋅ y
( x, x ) ⋅ ( y , y )
где 0 ≤ ϕ < π .
Данное нами определение угла корректно, ибо в силу неравенства
Коши-Буняковского дробь, стоящая в правой части последнего равенства, по
модулю не превосходит единицы т.е. cos ϕ ≤ 1 .
Будем называть два произвольных элемента x и
y
евклидова
∧
пространства E ортогональными, если ( x, y ) = 0 (в этом случае cos ( x, y) = 0 ).
x+ y
x
Так выглядит сумма
геометрических векторов
y
Проводя аналогию с геометрическими объектами, назовем сумму
x + y двух ортогональных элементов x и y гипотенузой прямоугольного
треугольника, построенного на элементах x и y .
Во всяком евклидовом пространстве справедлива теорема Пифагора:
квадрат гипотенузы равен сумме квадратов катетов. В самом деле, поскольку
x и
y
ортогональны, а, следовательно, ( x, y ) = 0 , то в силу аксиом и
определения длины
2
x + y = ( x + y , x + y ) = ( x , x ) + 2( x , y ) + ( y , y ) =
2
= ( x, x ) + ( y , y ) = x + y 2 .
Ранее
было
введено
понятие
базиса
n -мерного
векторного
пространства. Все базисы в произвольном векторном пространстве являлись
равноправными, и у нас не было оснований предпочитать один базис
другому. В евклидовом пространстве существуют специальные, особо
63
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
удобные базисы, называемые ортонормированными базисами. Эти базисы
играют ту же роль, что и декартов прямоугольный базис в аналитической
геометрии.
Определение. Говорят, что n элементов e1 , e2 ,..., en n -мерного евклидова
пространства E образуют ортонормированный базис этого пространства,
если эти элементы попарно ортогональны и норма каждого из них равна
единице, т.е. если
⎧1 при i = k ,
(ei , ek ) = ⎨
⎩0 при i ≠ k .
Ценность понятия ортонормированного базиса была бы невелика, если
бы не следующая теорема.
Теорема. Во всяком n -мерном евклидовом пространстве существует
ортонормированный базис.
Произвольный
ортонормированный
базис
любого
евклидова
пространства обладает свойствами, аналогичными свойствам декартова
прямоугольного базиса (на плоскости и в пространстве геометрических
векторов).
Так,
произведение
например,
двух
любых
в
ортонормированном
векторов
равно
базисе
сумме
скалярное
произведений
соответствующих координат этих векторов.
Еще раз подчеркнем, что существенным отличием произвольных
векторных пространств от их частного случая, евклидовых пространств,
является то, что в векторном пространстве не определены метрические
соотношения между его элементами.
Примером ортонормированного базиса является система n - единичных
векторов e1 = (1,0,0,...,0,0), e2 = (0,1,0,...,0,0), ... , en = (0,0,0,...,0,1) .
64
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
§ 3.4 Линейные операторы
1. Оператор. Одно из фундаментальных понятий матричной алгебры –
понятие линейного оператора.
Рассмотрим два линейных пространства: R n размерности n и R m
размерности m .
Определение. Если задан закон, по которому каждому вектору x
пространства R n ставится в соответствие единственный вектор y
~
пространства R m , то говорят, что задан оператор A( x) , действующий из R n в
~
R m и записываются y = A ( x ) .
Оператор называется линейным, если для любых векторов x и y
пространства R n и любого числа λ выполняются соотношения:
1.
~
~
~
A( x + y ) = A( x ) + A ( y ) - свойство аддитивности операторов;
2.
~
~
A (λx ) = λA( x ) - свойство однородности операторов.
~
Вектор y = A ( x ) называется образом вектора x , а сам вектор x прообразом вектора y .
Если пространства R n и R m совпадают, то оператор А отображает R n в
себя. Именно такие операторы будет изучаться в дальнейшем.
Выбираем в пространстве R n базис e1 , e2 ,..., en , то разложение вектора x по
данному базису имеет вид
x = x1e1 + x 2 e2 + ... + x n en
~
В силу линейности оператора A( x) получаем
~
~
~
~
A( x ) = x1 A(e1 ) + x 2 A(e2 ) + ... + x n A(en ) .
~
Поскольку A(ei ) (i = 1,2,..., n ) - также вектор из R n , то его можно
разложить по базису e1 , e2 ,..., en т.е.
~
A(ei ) = a1i e1 + a 2i e2 + ... + a ni en
(i = 1,2,..., n ) .
Тогда
65
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
~
A(ei ) = x1 (a11e1 + a 21e2 + ... + a n1en ) + x 2 (a12 e1 + a 22 e2 + ... + a n 2 en ) +
+ ... + x n (a1n e1 + a 2 n e2 + ... + a nn en ) = e1 (a11 x1 + a12 x 2 + ... + a1n x n ) +
(19)
+ e2 (a 21 x1 + a 22 x 2 + ... + a 2 n x n ) + ... + en (a n1 x1 + a n 2 x 2 + ... + a nn x n ).
~
С другой стороны, вектор y = A( x) , имеющий в том же базисе e1 , e2 ,..., en
координаты ( y1 , y 2 ,..., y n ) можно записать так:
~
A( x ) = y1e1 + y 2 e2 + ... + y n en .
(20)
Ввиду единственности разложение векторов по базису равны правые
части равенства (19) и (20), откуда
⎧ y1 = a11 x1 + a12 x 2 + ... + a1n x n ,
⎪ y = a x + a x + ... + a x ,
⎪ 2
21 1
22 2
21n n
⎨
⎪ ... ... ... ... ... ... ... ... ... ...
⎪⎩ y n = a n1 x1 + a n 2 x 2 + ... + a nn x n .
(21)
Матрица
⎛ a11
⎜
⎜a
A = ⎜ 21
...
⎜
⎜a
⎝ n1
a12
a 22
...
...
...
an2
...
...
a1n ⎞
⎟
a2n ⎟
... ⎟
⎟
a nn ⎟⎠
~
Называется матрицей оператора A в базисе e1 , e2 ,..., en , а ранг r матрицы
~
A - рангом оператора A .
Таким образом, каждому линейному оператору соответствует матрица
в данном базисе. Справедливо и обратное.
~
Связь между вектором x и его образом y = A ( x ) можно выразить в
матричной форме уравнением
Y = AX ,
где А- матрица линейного оператора, X = ( x1 ,..., xn ) , Y = ( y1 , y 2 ,..., y n ) матрицы столбцы из координат векторов x и y .
~
Пример 6. Пусть в пространстве R n линейный оператор A в базисе
e1 , e2 , e3 задан матрицей
66
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
⎛ 3 2 4⎞
⎟
⎜
A = ⎜ −1 5 6⎟ .
⎜ 1 8 2⎟
⎠
⎝
~
Найти образ y = A ( x ) вектора x = 4e1 ,−3e2 + e3 .
Решение. По формуле (21) имеем
⎛ y1 ⎞ ⎛ 3 2 4 ⎞⎛ 4 ⎞ ⎛ 10 ⎞
⎜ ⎟ ⎜
⎟⎜ ⎟ ⎜
⎟
⎜ y 2 ⎟ = ⎜ − 1 5 6 ⎟⎜ − 3 ⎟ = ⎜ − 13 ⎟ .
⎜ y ⎟ ⎜ 1 8 2 ⎟⎜ 1 ⎟ ⎜ − 18 ⎟
⎝ 3⎠ ⎝
⎠⎝ ⎠ ⎝
⎠
Следовательно y = 10e1 ,−13e2 − 18e3 .
Зависимости между матрицами одного и того же оператора в разных
базисах выражается формулой
A€ = C −1 AC ,
~
где A, A€ матрицы оператора A в базисах e1 , e2 ,..., en и e€1 , e€2 ,..., e€n , С- матрица
перехода от старого базиса к новому.
~
Пример 7. В базисе e1 , e2 оператор A имеет матрицу
⎛17 6 ⎞
⎟⎟ .
A = ⎜⎜
⎝ 6 8⎠
~
Найти матрицу оператора A в базисе e1∗ = e1 ,−2e2 , e2∗ = 2e1 ,+e2 .
Решение. Определяем матрицу перехода
⎛ 1 2⎞
⎟⎟ .
C = ⎜⎜
⎝− 2 1⎠
Находим обратную матрицу к матрицу перехода
1 ⎛1 − 2⎞
⎟.
C −1 = ⎜⎜
5 ⎝ 2 1 ⎟⎠
Следовательно
1 ⎛ 1 − 2 ⎞⎛17 6 ⎞⎛ 1 2 ⎞ ⎛ 5 0 ⎞
⎟⎜
⎟⎜
⎟=⎜
⎟.
A€ = C −1 AC = ⎜⎜
5 ⎝ 2 1 ⎟⎠⎜⎝ 6 8 ⎟⎠⎜⎝ − 2 1 ⎟⎠ ⎜⎝ 0 20 ⎟⎠
2.Собственные векторы и собственные значения оператора.
Определение. Всякий ненулевой вектор x ≠ 0 называется собственным
~
вектором линейного оператора A (квадратной матрицы A), если найдется
такое число λ, что будет выполняться равенство
67
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
~
A ( x ) = λx
~
Число λ называется собственным значением линейного оператора A
(матрицы A), соответствующим вектору x . Матрица A имеет порядок n.
В математической экономике большую роль играют так называемые
продуктивные матрицы. Доказано, что матрица A является продуктивной
тогда и только тогда, когда все собственные значения матрицы A по модулю
меньше единицы.
Для нахождения собственных значений матрицы A перепишем
~
равенство A(x ) = λx в векторном виде
( A − λE ) X
= 0,
(22)
где E- единичная матрица n-го порядка, X = (x1 , x2 ,..., xn )T -вектор-столбец
или в координатной форме:
⎧a11 − λx1 + a12 x 2 + ... + a1n x n = 0,
⎪a x + (a − λ )x + ... + a x = 0,
⎪ 21 1
22
2
2n n
⎨
⎪...............................................
⎪⎩a nm1 x1 + a n 2 x 2 + ... + (a nn − λ )x n = 0.
(23)
Получили систему линейных однородных уравнений, которая имеет
ненулевые решения тогда и только тогда, когда определитель этой системы
равен нулю, т.е.
a 11 − λ
A- λE =
a 21
a1 2
... a 1 n
a 2 2 − λ ... a 2 n
.... ... ... ... ... ... ...
a n1
an2
... a n n − λ
=0.
Получили уравнение n-ой степени относительно неизвестной λ,
которое называется характеристическим уравнением матрицы A, многочлен
A-λ E
называется характеристическим многочленом матрицы A, а его
корни - характеристическими числами, или собственными значениями,
матрицы A.
68
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
Для нахождения собственных векторов матрицы A в векторное
уравнение (22) или в соответствующую систему однородных уравнений (23)
нужно подставить найденные значения λ и решать обычным образом.
Пример 8. Найти собственные значения и собственные векторы
матрицы
⎛ 3 -1 0 0 ⎞
⎜
⎟
⎜1 1 0 0 ⎟
A= ⎜
.
3 0 -5 -3 ⎟
⎜
⎟
⎜ 4 -1 3 1 ⎟
⎝
⎠
Решение. Вычислим определитель матрицы A - λE:
⎛3 − λ
⎜
⎜ 1
A- λE = ⎜
3
⎜
⎜ 4
⎝
-1
0
0⎞
⎟
1- λ
0
0 ⎟
=
0 -5-λ -3 ⎟
⎟
-1
3 1 - λ ⎟⎠
λ 2 − 4λ + 4
=
3
1+ λ
0
- 5- λ
3
0
⎛
⎜ 2
⎜ λ - 4λ + 4
⎜
3
⎜
⎝ 1+ λ
-1
0
1- λ
0
0 - 5- λ
-1
3
0⎞
⎟
0⎟
=
-3 ⎟
⎟
1- λ ⎠
0
3
5+λ
- 3 = − ( λ 2 − 4λ + 4 )
.
3 1- λ
1- λ
Итак, A - λ E = (λ - 2)2 ⋅ (λ+2)2. Корни характеристического уравнения
A - λ E =0 - это числа
λ1 = 2 и λ2 = -2. Другими словами, мы нашли
собственные значения матрицы A. Для нахождения собственных векторов
матрицы A подставим найденные значения λ в систему (23): при λ = 2 имеем
систему линейных однородных уравнений
⎧ x1 − x 2 = 0 ,
⎪ x − x = 0,
⎪ 1
2
⇒
⎨
⎪3x1 − 7 x3 − 3x 4 = 0,
⎪⎩4 x1 − x 2 + 3x3 − x 4 = 0.
x1 = x 2 ,
⎧
⎪
⎨3x 2 = 7 x3 + 3x 4 ,
⎪ 5 x + x = 0.
3
4
⎩
Следовательно, собственному значению λ = 2 отвечают собственные векторы
вида α (8, 8, -3, 15), где α - любое отличное от нуля действительное число.
При λ = -2 имеем:
69
© КРСУ Усенов И.А., Усенова Р.К.
⎛5 - 1
⎜
1
3
A - λE = A +2E = ⎜⎜
3
0
⎜
⎝4 - 1
Кафедра высшей математики
0⎞
⎟
0⎟
∼
-3 -3⎟
⎟
3
3⎠
0
0
⎛1 3
⎜
⎜0 1
⎜
⎝0 0
0 0⎞
⎟
0 0⎟ ,
⎟
1 1⎠
и поэтому координаты собственных векторов должны удовлетворять системе
уравнений
⎧ x1 + 3 x 2 = 0 ,
⎪
⎨ x 2 = 0,
⎪ x + x = 0.
4
⎩ 3
Поэтому собственному значению λ = -2 отвечают собственные векторы вида
β (0, 0,-1, 1), где β - любое отличное от нуля действительное число.
~
Наиболее простой вид принимает матрица A линейного оператора A ,
имеющего
n
линейно
независимых
собственных
векторов
e1 , e2 ,..., en
соответственно соответствующих собственным значениям λ1 , λ2 ,..., λn .
Векторы e1 , e2 ,..., en
~
примем за базисные. Тогда A(ei ) = λi ei (i = 1,2,..., n ) или
~
A(ei ) = a1i e1 + a 2i e2 + ... + a ni en = λi ei ,
~
откуда a ij = 0 если i ≠ j , и a ij = λi i = j . Таким образом, матрица оператора A
в базисе, состоящих из его собственных векторов, является диагональной и
имеет вид:
⎛ λ1
⎜
⎜0
A=⎜
...
⎜
⎜0
⎝
0
λ2
...
...
...
0
...
...
0 ⎞
⎟
0 ⎟
.
... ⎟
⎟
λ n ⎟⎠
⎛1 4⎞
~
⎟⎟ линейного оператора A к
Пример 9. Привести матрицу A = ⎜⎜
⎝9 1⎠
диагональному виду.
Решение. Составим характеристическое уравнение
A − λE =
1− λ
4
= 0 или λ 2 − 2λ − 35 = 0 ,
9
1− λ
~
откуда собственные значения линейного оператора A λ1 = −5, λ2 = 7 .
70
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
Для нахождения собственных векторов матрицы A подставим найденные
значения λ в систему
⎧(1 − λ )x1 + 4 x 2 = 0,
⎨
⎩9 x1 + (1 − λ )x 2 = 0.
При λ1 = −5, имеем
⎧ x = c,
⎧6 x1 + 4 x 2 = 0 ⎧3 x1 + 2 x 2 = 0 ⎧
⇒⎨
⇒ ⎨3 x1 + 2 x 2 = 0 ⇒ ⎨ 1
⎨
⎩9 x1 + 6 x 2 = 0 ⎩3 x1 + 2 x 2 = 0 ⎩
⎩ x 2 = −1,5c.
При λ1 = 7, имеем
⎧ x = c1 ,
⎧− 6 x1 + 4 x 2 = 0 ⎧− 3 x1 + 2 x 2 = 0 ⎧
⇒⎨
⇒ ⎨3 x1 − 2 x 2 = 0 ⇒ ⎨ 1
⎨
⎩ 9 x1 − 6 x 2 = 0
⎩ 3 x1 − 2 x 2 = 0
⎩ x 2 = 1,5c1 .
⎩
Таким образом, найдены соответствующие собственные векторы
x (1) = (c,−1,5c ),
x (2 ) = (c1 ,1,5c1 ) .
Так, как координаты векторов
x (1) , x (2 )
не
пропорциональны, то векторы x (1) и x (2 ) линейно независимы. Поэтому
матрица A в базисе x (1) , x (2 ) имеет диагональный вид:
⎛ − 5 0⎞
⎟⎟.
A = ⎜⎜
⎝ 0 7⎠
3. Квадратичные формы. При решении различных прикладных задач
часто приходится исследовать квадратичные формы.
Определение. Квадратичной формой L(x1 , x 2 ,.., xn ) от n переменных
называется сумма, каждый член которой является либо квадратом одной из
переменных, либо произведением двух различных переменных, взятых с
некоторыми коэффициентом:
n
n
L( x1 , x 2 ,.., x n ) = ∑∑ aij xi x j .
i =1 j =1
Предполагается, что коэффициенты квадратичной формы aij действительные числа, причем a ij = a ji . Матрица, A = (aij ), (i, j = 1,2,...n )
составленная из коэффициентов, называется матрицей квадратичной формы.
Матричный запись квадратичный формы имеет вид:
L = X ′AX ,
где X = (x1 , x2 ,.., xn ) - матрица- столбец переменных.
71
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
В самом деле
⎛ a11
⎜
⎜a
L = ( x1 , x 2 ,.., x n )⎜ 21
...
⎜
⎜a
⎝ n1
a12
a 22
...
an2
⎛ ∑ a1 j x j ⎞
... a1n ⎞⎛ x1 ⎞
⎟
⎜
⎟⎜ ⎟
... a 2 n ⎟⎜ x 2 ⎟
⎜ ∑ a2 j x j ⎟
= ( x1 , x 2 ,.., x n )⎜
⎟=
...
... ... ⎟⎜ ... ⎟
⎟
⎜
⎟⎜ ⎟
⎜ ∑a x ⎟
... a nn ⎟⎠⎜⎝ x n ⎟⎠
nj
j
⎠
⎝
n
n
n
j =1
j =1
j =1
n
n
= ∑ a1 j x1 x j + ∑ a 2 j x 2 x j + ... + ∑ a nj x n x j = ∑∑ aij xi x j .
i =1 j =1
Пример 10. Дана квадратичная форма
L( x1 , x 2 , x3 ) = 4 x12 − 12 x1 x 2 − 10 x1 x3 + x 22 − 3 x32 .
Записать матричной форме.
Решение. Найдем матрицу квадратичной формы. Ее диагональные
элементы равны коэффициентам при квадратах переменных, а другиеполовинам соответствующих коэффициентов квадратичной формы. Поэтому
⎛ 4 − 6 − 5 ⎞⎛ x1 ⎞
⎜
⎟⎜ ⎟
L = ( x1 , x 2 , x3 )⎜ − 6 1
0 ⎟⎜ x 2 ⎟ .
⎜ − 5 0 − 3 ⎟⎜ x ⎟
⎝
⎠⎝ 3 ⎠
Выясним, как изменяется квадратичная форма при невырожденном
линейном преобразовании переменных.
Пусть матрицы-столбцы переменных X = (x1 , x2 ,.., xn ) и Y = ( y1 , y 2 ,.., y n )
связаны линейным соотношением X = CY , где C = (cij ), (i, j = 1,2,...n ) есть
некоторая невырожденная матрица n - го порядка. Тогда квадратичная форма
′
L = X ′AX = (CY ) A(CY ) = (Y ′C ′)A(CY ) = Y ′(C ′AC )Y .
Итак, при невырожденном линейном преобразовании X = CY матрица
квадратичной формы принимает вид:
A ∗ = C ′AC .
Пример 11. Дана квадратичная форма
L( x1 , x 2 ) = 2 x12 + 4 x1 x 2 − 3x 22 .
Найти квадратичную форму L( y1 , y 2 ) , полученную из данной линейным
преобразованием x1 = 2 y1 − 3 y 2 , x2 = y1 + y 2 .
72
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
Решение. Матрица данной квадратичной формы
⎛2 2 ⎞
⎟⎟ .
A = ⎜⎜
⎝ 2 − 3⎠
Матрица линейного преобразования
⎛ 2 − 3⎞
⎟⎟.
C = ⎜⎜
⎝1 1 ⎠
Следовательно, матрица искомой квадратичной формы:
⎛ 2 1⎞⎛ 2 2 ⎞⎛ 2 − 3 ⎞ ⎛ 13 − 17 ⎞
⎟⎟ = ⎜⎜
⎟⎟⎜⎜
⎟⎟⎜⎜
⎟.
A∗ = C ′AC = ⎜⎜
3 ⎟⎠
⎝ − 3 1⎠⎝ 2 − 3 ⎠⎝ 1 1 ⎠ ⎝ − 17
Тогда квадратичная форма L( y1 , y 2 ) имеет вид:
L( y1 , y 2 ) = 13 y12 − 34 y1 y 2 + 3 y 22 .
n
n
Определение. Квадратичная форма L(x1 , x 2 ,.., x n ) = ∑∑ aij xi x j называется
i =1 j =1
канонической, если все ее коэффициенты a ij = 0 при i ≠ j .
Из данной определение получим канонический вид квадратичной формы
n
L = a11 x12 + a 22 x 22 + ... _ + a nn x n2 = ∑ aii xi2 .
i =1
Матрица канонической квадратичной формы является диагональной
⎛ a11
⎜
⎜ 0
A=⎜
...
⎜
⎜ 0
⎝
0
a 22
...
...
...
0
...
...
0 ⎞
⎟
0 ⎟
.
... ⎟
⎟
a nn ⎟⎠
Пример 12. Привести к каноническому виду квадратичную форму
L( x1 , x 2 , x3 ) = x12 − 3 x1 x 2 + 4 x1 x3 + 2 x 2 x3 + x32 .
Решение. Вначале выделим полный квадрат при x1 , коэффициент при
квадрате которой отличен от нуля:
2
2
⎡ 2
⎛1
⎞ ⎛1
⎞ ⎤ ⎛1
⎞
L = ⎢ x1 − 2 x1 ⎜ (3x 2 − 4 x3 )⎟ + ⎜ (3x 2 − 4 x3 )⎟ ⎥ − ⎜ (3x 2 − 4 x3 )⎟ +
⎝2
⎠ ⎝2
⎠ ⎦⎥ ⎝ 2
⎠
⎣⎢
2
9
3
⎛
⎞
+ 2 x 2 x3 + x = ⎜ x1 − x 2 + 2 x3 ⎟ − x 22 + 8 x 2 x3 − 3x32 .
4
2
⎝
⎠
2
3
Теперь выделим полный квадрат при x 2 , коэффициент при квадрате
которой отличен от нуля
73
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
2
9⎛
32
256 2 ⎞ 9 256 2
3
⎛
⎞
L = ⎜ x1 − x 2 + 2 x3 ⎟ − ⎜ x 22 −
x 2 x3 +
x3 ⎟ +
x3 − 3 x32 =
4⎝
9
81
2
⎝
⎠
⎠ 4 81
2
2
3
9⎛
16 ⎞
37 2
⎛
⎞
= ⎜ x1 − x 2 + 2 x3 ⎟ − ⎜ x 2 − x3 ⎟ +
x3 .
2
4⎝
9 ⎠
9
⎝
⎠
Итак, невырожденное линейное преобразование
y1 = x1 −
3
16
x 2 + 2 x3 , y 2 = x 2 − x3 , y 3 = x3 .
2
9
Приводить данную квадратичную форму к каноническому виду
9
37 2
L1 ( y1 , y 2 , y 3 ) = y12 − y 22 +
y3 .
4
9
Канонический вид квадратичной формы не является однозначно
определенным, так, как одна и та же квадратичная форма может быть
приведена к каноническому виду многими способами. Однако полученные
различными способами канонические формы обладают рядом общих
свойств.
Например, квадратичную форму L в примере 12 можно было привести
к виду
L2 ( y1 , y 2 , y 3 ) =
37 2
y1 + y 22 − y 32 ,
4
применив невырожденное линейное преобразование
y1 = x1 , y 2 = 2 x1 + x 2 + x3 , y 3 =
7
x1 + x 2 .
2
Как видим, число положительных и отрицательных коэффициентов
сохранилось.
Следует отметить, что ранг матрицы квадратичной формы, называется
рангом квадратичной формы и равен числу отличных от нуля коэффициентов
канонической формы и не меняется при линейных преобразованиях.
Определение. Квадратичная форма L(x1 , x 2 ,.., xn ) называется
положительно (отрицательно) определенной, если при всех значениях
переменных, из которых хотя бы одно отлично от нуля
L( x1 , x 2 ,.., x n ) > 0
( L(x1 , x2 ,.., xn ) < 0 ).
Например, квадратичная форма
74
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
L1 ( x1 , x 2 , x3 ) = 3 x12 + 4 x 22 + 9 x32
является положительно определенной, а форма
L2 ( x1 , x 2 , x3 ) = − x12 + 2 x1 x 2 − x32
является отрицательно определенной.
Теорема. Для того чтобы квадратичная форма L = X ′AX , было
положительно (отрицательно) определенной необходимо и достаточно,
чтобы все собственные значений матрицы A были положительны
(отрицательны).
Пример 13. Показать, что квадратичная форма
L( x1 , x 2 ) = 13 x12 − 6 x1 x 2 + 5 x 22
является положительно определенной.
Решение. Матрица квадратичной формы имеет вид
⎛ 13 − 3 ⎞
⎟⎟ .
A = ⎜⎜
⎝− 3 5 ⎠
Составим характеристическое уравнение
A − λE =
13 − λ − 3
= 0 или λ 2 − 18λ + 56 = 0 ,
−3 5−λ
откуда собственные значения матрицы
λ1 = 14, λ2 = 4 . Так, как корни
характеристического уравнения матрицы A положительны, то на основании
приведенной теоремы квадратичная форма
L( x1 , x 2 ) = 13 x12 − 6 x1 x 2 + 5 x 22
положительно определенная.
4.Линейная модель обмена. В качестве примера математической
модели экономического процесса, приводящейся к понятию собственного
вектора и собственного значения матрицы, рассмотрим линейную модель
обмена.
Пусть имеется n стран S1 , S 2 ,..., S n национальный доход каждой из
которых равен соответственно x1 , x 2 ,.., xn . Обозначим коэффициентами aij
долю национального дохода, которую страна S J тратит на покупку товаров у
75
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
стран S i . Будем считать, что весь национальный доход тратится на закупку
товаров либо внутри стран, либо на импорт из других стран, т.е.
n
∑a
i =1
ij
= 1,
( j = 1,2,.., n ) .
(24)
Рассмотрим матрицу
⎛ a11
⎜
⎜a
A = ⎜ 21
...
⎜
⎜a
⎝ n1
a12
a 22
...
...
...
an2
...
...
a1n ⎞
⎟
a2n ⎟
,
... ⎟
⎟
a nn ⎟⎠
которая получила название структурной матрицы торговли. В соответствии с
(24) сумма элементов любого столбца матрицы A равна 1.
Для любой страны S J выручка от внутренней и внешней торговли
составит:
pi = ai1 x1 + ai 2 x 2 + ... + ain x n .
Для сбалансированной торговли необходима бездефицитность
торговли каждой страны S i , т.е. выручка от торговли каждой страны должна
быть не меньше ее национального дохода:
pi ≥ xi (i = 1,2,..., n ) .
Если считать, что pi > xi (i = 1,2,..., n ) , то получаем систему неравенств
⎧a11 x1 + a12 x 2 + ... + a1n x n > x1 ,
⎪a x + a x + ... + a x > x ,
⎪ 21 1
22 2
2n n
2
⎨
⎪...............................................
⎪⎩a n1 x1 + a n 2 x 2 + ... + a nn x n > x n .
(25)
Сложив все неравенства системы (25), получим после группировки
x1 (a11 + a 21 + ... + a n1 + ) + x 2 (a12 + a 22 + ... + a n 2 + ) + ... +
+ x n (a1n + a 2 n + ... + a nn + ) > x1 + x 2 + ... + x n .
Учитывая (24), выражения в скобках равны единице, и мы приходим к
противоречивому неравенству
x1 + x 2 + ... + x n . > x1 + x 2 + ... + x n .
Таким образом, неравенство pi > xi (i = 1,2,..., n ) невозможно, и условие
pi ≥ xi (i = 1,2,..., n ) принимает вид pi = xi (i = 1,2,..., n ) . С экономической точки
76
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
зрения это понятно, так как все страны не могут одновременно получать
прибыль.
Вводя вектор x = (x1 , x2 ,.., xn ) национальных доходов страны, получим
матричное уравнение
AX = X ,
где X - матрица- столбец из координат вектора x т.е. задача свелась к
отысканию собственного вектора матрицы A , отвечающего собственному
значению λ = 1 .
Пример 14. Структурная матрица торговли трех стран S1 , S 2 , S 3 имеет
вид
⎛1
⎜ 3
A=⎜1
⎜ 3
⎜1
⎝ 3
1
1
1
4
2
4
1 ⎞
2⎟
1 ⎟.
2⎟
0 ⎟
⎠
Найти национальные доходы стран для сбалансированной торговли.
Решение. Находим собственный вектор x , отвечающий собственному
значению λ = 1 . Для этого решим систему
⎛1
⎜ 3
⎜1
⎜ 3
⎜1
⎝ 3
1
1
1
4
2
4
1 ⎞ x
2 ⎟⎛⎜ 1 ⎞⎟ ⎛⎜ 0 ⎞⎟
1 ⎟⎜ x 2 ⎟ = ⎜ 0 ⎟
2⎟
⎜ ⎟ ⎜ ⎟
0 ⎟⎝ x3 ⎠ ⎝ 0 ⎠
⎠
Методом Гаусса. Найдем x1 = c, x 2 = 2c, x3 = c т.е. x = ⎛⎜ c, 2c, c ⎞⎟ .
3
2
3
⎝2
⎠
Полученный результат означает, что сбалансированность торговли
трех стран достигается при векторе национальных доходов x = ⎛⎜ c, 2c, c ⎞⎟ ,т.е.
3
⎝2
⎠
при соотношении национальных доходов стран 3 2 : 2 : 1 или 3 : 4 : 2.
77
© КРСУ Усенов И.А., Усенова Р.К.
Кафедра высшей математики
ЛИТЕРАТУРЫ
1.Тышкевич Р. И., Феденко А. С. Линейная алгебра и аналитическая
геометрия. Минск: Вышейш. школа, 1988.
2. Высшая математика: Общий курс / Под ред. А. И. Яблонского,
Минск: Вышейш. школа, 1993.
3. Баврин И. И., Матросов В. Л. Общий курс высшей математики.
М.: Просвещение, 1995.
4. Рублев А. Н. Линейная алгебра. М.: Высшая школа, 1998.
5.Е.С. Федорова, Т.А. Шемякина Линейная алгебра. Бишкек-2002г.
6. Ланкастер К. Математическая экономика. М.: 2002.
7. Высшая математика для экономистов // под ред. Н.Ш. Кремера.М.:ЮНИТИ, 2002г.
8.Общий курс высшей математики для экономистов // под ред. В.И.
Ермакова.-М.:ИНФРА-М, 2008г.
9. Минорский В. П. Сборник задач по высшей математике. М.: Наука, 87
10. Сборник задач по высшей математики для экономистов // под ред. В.И.
Ермакова.-М.: ИНФРА-М, 2008г.
78
Download