Характеристическое уравнение
В п. 13.1 было введено определение собственного значения и собственного вектора матрицы. Пусть — собственный вектор квадратной матрицы А порядка n. Тогда имеет место матричное уравнение
или
где λ — собственное значение матрицы А, а E и — соответственно единичная матрица и нулевой вектор-столбец. Уравнение (15.17) эквивалентно системе однородных уравнений
В уравнениях (15.18) aij — элементы матрицы А, xj — координаты собственного вектора х. Поскольку собственный вектор не является нулевым, то однородная система (15.18) должна иметь ненулевое решение, т.е. в силу следствия 2 (см. выше) определитель этой системы равен нулю:
Определитель системы однородных уравнений (15.18) называется характеристическим многочленом, а уравнение (15.19) — характеристическим уравнением матрицы А.
Уравнение (15.19) имеет степень n относительно неизвестной λ. Его корни являются собственными числами матрицы А. Определив набор этих чисел, для каждого из них можно найти соответствующий собственный вектор как решение однородной системы (15.18).
Пример 2. Найти собственные числа и собственные векторы матрицы
Решение. Характеристическое уравнение для этой матрицы имеет вид
откуда, раскрывая определитель, получаем
Корни этого уравнения суть λ1 = 2, λ2 = 5. Для нахождения собственных векторов подставим найденные собственные значения в систему однородных уравнений (15.18) при n = 2 с соответствующими элементами заданной матрицы А. Собственный вектор, соответствующий собственному значению λ1 = 2, является решением системы
Пo сути дела, это одно уравнение, поскольку определитель системы равен нулю. Полагая x2 = b свободной переменной, получаем первый собственный вектор 1 = (—2b, b) = b (-2, 1). Подстановка второго собственного значения λ2 = 5 приводит к системе уравнений
которая через свободную переменную x2 = с определяет второй собственный вектор матрицы А: 2 = (с, с) = с (1, 1).
Поскольку b и с — произвольные числа, то одному собственному значению может соответствовать несколько собственных векторов разной длины. Например, собственные векторы, соответствующие фундаментальным решениям однородных систем (в данном случае их будет по одному на каждое собственное значение), имеют вид 1 = (-2, 1), 2 = (1, 1).
ПРИМЕНЕНИЕ ЭЛЕМЕНТОВ ЛИНЕЙНОЙ АЛГЕБРЫ В ЭКОНОМИКЕ
Использование алгебры матриц
Использование элементов алгебры матриц является одним из основных методов решения многих экономических задач. Особенно этот вопрос стал актуальным при разработке и использовании баз данных: при работе с ними почти вся информация хранится и обрабатывается в матричной форме.
Матричные вычисления
Рассмотрим типичные задачи, использующие понятие вектора и его свойства.
1. Предприятие выпускает ежесуточно четыре вида изделий, основные производственно-экономические показатели которых приведены в табл. 16.1.
Требуется определить следующие ежесуточные показатели: расход сырья S, затраты рабочего времени Т и стоимость Р выпускаемой продуции предприятия.
Решение. По данным табл. 16.1 составим четыре вектора, характеризующие весь производственный цикл:
= (20, 50, 30,40) — вектор ассортимента,
= (5, 2, 7, 4) — вектор расхода сырья,
= (10, 5, 15, 8) — вектор затраты рабочего времени,
= (30, 15, 45, 20) — ценовой вектор.
Тогда искомые величины будут представлять собой соответствующие скалярные произведения вектора ассортимента на три других вектора, т.е.
2. Предприятие выпускает 4 вида изделий с использованием 4-х видов сырья. Нормы расхода сырья даны как элементы матрицы А:
Требуется найти затраты сырья на каждый вид изделия при заданном плане их выпуска: соответственно 60, 50, 35 и 40 ед.
Решение. Составим вектор-план выпуска продукции
Тогда решение задачи дается вектором затрат, координаты которого и являются величинами затрат сырья по каждому его виду; этот вектор затрат вычисляется как произведение вектора на матрицу А:
3. Пусть затраты 4-х видов сырья на выпуск 4-х видов продукции характеризуются матрицей А, приведенной в предыдущей задаче. Требуется найти: а) общие затраты на сырье для каждого вида продукции и его перевозку; б) общие затраты на сырье и его транспортировку при условии заданного вектора-плана предыдущей задачи, если известны себестоимости каждого вида сырья и его доставки (соответственно 4, 6, 5, 8 и 2, 1, 3, 2 ден. ед.).
Решение. Составим матрицу себестоимостей сырья и его доставки (соответственно 1-я и 2-я строки):
Тогда ответ на первый вопрос задачи дается в виде произведения матрицы А на транспонированную матрицу CT:
Суммарные затраты на сырье и его доставку (в денежных единицах) при векторе-плане выпуска продукции = (60, 50, 35, 40) определяются произведением вектора на матрицу АСT:
4. В табл. 16.2 приведены данные о дневной производительности 5 предприятий, выпускающих 4 вида продукции с потреблением 3-х видов сырья, а также продолжительность работы каждого предприятия в году и цена каждого вида сырья.
Требуется определить:
1) годовую производительность каждого предприятия по каждому виду изделий;
2) годовую потребность каждого предприятия по каждому виду сырья;
3) годовую сумму кредитования каждого предприятия для закупки сырья, необходимого для выпуска продукции указанных видов и количеств.
Решение. Нужно составить матрицы, характеризующие весь интересующий нас экономический спектр производства, а затем при помощи соответствующих операций над ними получить решение данной задачи. Прежде всего приведем матрицу производительности предприятий по всем видам продукции:
Каждый столбец этой матрицы соответствует дневной производительности отдельного предприятия по каждому виду продукции. Следовательно, годовая производительность j-гопредприятия по каждому виду продукции получается умножением j-гo столбца матрицы А на количество рабочих дней в году для этого предприятия (j = 1, 2, 3, 4, 5). Таким образом, годовая производительность каждого предприятия по каждому из изделий описывается матрицей
Матрица затрат сырья на единицу изделия (эти показатели по условию одинаковы для всех предприятий) имеет вид
Дневной расход по типам сырья на предприятиях описывается произведением матрицы В на матрицу А:
где i-я строка соответствует номеру типа сырья, а j-й столбец — номеру предприятия согласно табл. 16.2 (i = 1, 2, 3; j = 1, 2, 3, 4, 5). Ответ на второй вопрос задачи получим по аналогии с матрицей Агод умножением столбцов матрицы ВА на соответствующие количества рабочих дней в году для предприятий — это годовая потребность каждого предприятия в каждом виде сырья:
Введем вектор стоимости сырья
Тогда стоимость общего годового запаса сырья для каждого предприятия получается умножением вектора на матрицу ВAгод:
Следовательно, суммы кредитования предприятий для закупки сырья определяются соответствующими компонентами вектора .
5. Отрасль состоит из п предприятий, выпускающих по одному виду продуции каждое; обозначим объем продукции i-го предприятия через xi. Каждое из предприятий отрасли для обеспечения своего производства потребляет часть продукции, выпускаемой им самим и другими предприятиями. Например, в отрасли электротехнического оборудования часть продукции предприятий, выпускающих электродвигатели, силовые кабели, электрокары и т.д., употребляется практически всей отраслью. Пусть aij — доля продукции i-го предприятия, потребляемая j-м предприятием для обеспечения выпуска своей продукции объема xj. Возникает естественный вопрос о величине yi — количестве продукции i-го предприятия, предназначенной для реализации вне данной отрасли (объем конечного продукта). Эта величина легко может быть подсчитана по формуле
Введем в рассмотрение матрицу порядка n, описывающую внутреннее потребление отрасли:
Тогда вектор конечного продукта является решением матричного уравнения
или с использованием единичной матрицы Е получаем
Рассмотрим конкретный пример при п = 3. Пусть вектор выпуска продукции отрасли и матрица внутреннего потребления имеют соответственно вид
Используя формулу (16.1) и правило сложения матриц, получаем вектор объемов конечного продукта, предназначенного для реализации вне отрасли, состоящей из 3-х предприятий: