Студопедия.Орг Главная | Случайная страница | Контакты | Мы поможем в написании вашей работы!  
 

Краткие сведения, необходимые для выполнения и защиты контрольной работы. 1.1.1. Модель Леонтьева многоотраслевой экономики [1]



1.1.1. Модель Леонтьева многоотраслевой экономики [ 1 ]

Эффективное ведение народного хозяйства предполагает наличие баланса между отдельными отраслями. Каждая отрасль при этом выступает двояко: с одной стороны, как производитель некоторой продукции, а с другой - как потребитель продуктов, вырабатываемых другими отраслями. Для наглядного выражения взаимной связи между отраслями пользуются определенного вида таблицами, называемыми таблицами межотраслевого баланса.

Предполагаем, что имеется 3 различных отрасли O 1, O 2, O 3

каждая из которых производит свой продукт. В процессе производства своего продукта каждая отрасль нуждается в продукции других отраслей (производственное потребление). Пусть промежуток времени служит плановый год. Введем следующие обозначения: xi – общий объем продукции отрасли Oi за год (валовой выпуск отрасли); xij –объем продукции отрасли Oi, расходуемый отраслью Oj в процессе производства; yi — объем продукции отрасли Oi, предназначенный к потреблению в непроизводственной сфере (объем конечного потребления). В yi входят запасы, личное потребление граждан, обеспечение общественных потребностей, поставки на экспорт. Указанные величины можно свести в таблицу.

Производственное потребление Конечное потребление Валовой выпуск
x 11, x 12, x 13 y 1 x 1
x 21, x 22, x 23 y 2 x 2
x 31, x 32, x 33 y 3 x 3

Балансовый характер этой та6лицы выражается в том, что должно выполняться соотношение баланса

xi = xi 1 + xi 2 + xi 3 + yi

Экономист Леонтьев сделал допущение (гипотеза линейности): для выпуска любого объема xj продукции Oj необходимо затратить продукцию О i в количестве aij × xj, где aij - постоянный коэффициент (коэффициент прямых затрат). Вектор = (x 1, x 2, x 3) называется вектором валового выпуска, вектор = (y 1, y 2, y 3) – вектором конечного потребления, а матрица – матрицей прямых затрат. Соотношение

(1)

называется уравнением линейного межотраслевого баланса. Вместе с изложенной интерпретацией это соотношение (1) называют также моделью Леонтьева. Решение уравнения (1) в матричном виде:

(2)

Матрица A ≥ 0 (принимается, что все элементы матрицы неотрицательны) называется продуктивной, если для любого вектора y ≥ 0 существует решение x ≥ 0 уравнения (1). Известен второй критерий продуктивности: матрица A продуктивна тогда и только тогда, когда матрица (EA)–1 существует и неотрицательна. В этом случае и модель Леонтьева, определяемая матрицей A, тоже называется продуктивной.

Введем понятие собственных векторов и собственных значений матрицы A. Собственным значением матрицы называется число λ, удовлетворяющее уравнению

Û . (3)

Число собственных значений матрицы, при условии, что нет совпадающих значений, равно размерности матрицы (в данном случае - трем). Собственным вектором, отвечающим собственному значению (решение уравнения (3)), называется вектор , удовлетворяющий уравнению

Û (4)

где E - единичная матрица.

Числом Фробениуса F матрицы A называется максимальное из взятых по модулю собственных значений матрицы . Соответствующий этому собственному значению собственный вектор называется вектором Фробениуса

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

1) если сумма элементов любого столбца матрицы меньше единицы;

2) если число Фробениуса данной матрицы меньше единицы.

Поскольку для определения числа Фробениуса нужно только максимальное собственное значение положительной матрицы, эффективный способ вычисления F и дает алгоритм Ланцоша [ 2 ]. Из уравнения (4) следует, что

Û . (5)

Умножив последнее уравнение скалярно на вектор слева, получим , где круглыми скобками обозначено скалярное произведение. Окончательно,

(6)

Формулы (5) - (6) позволяют организовать итерационную процедуру вычисления собственного вектора и собственного значения. Пусть - некоторое начальное приближение для вектора Фробениуса. Положим для простоты модуль этого вектора равным 1: . Тогда соответствующее приближенное число Фробениуса вычисляется по формуле (6): . Следующее приближение для вектора запишем в виде (см. (5)) . Делить правую часть на не обязательно, поскольку вектор, как было условлено ранее, нормируется на 1. Продолжая процесс, получим последовательные приближения к числу и вектору Фробениуса. Заметим, что уравнения (5) и (6) справедливы для любого собственного значения и соответствующего вектора. Однако, в итерационном процессе последовательность значений сходится к F.

Зная число F, можно вычислить запас продуктивности матрицы. Запасом продуктивности матрицы A называется такое число a > 0, что все матрицы β A продуктивны, если β ≤ 1 + a и непродуктивны, если β > 1 + a. Матрица β A продуктивна, если матрица и обратная к ней имеют только положительные собственные значения. Поэтому особой точке отвечает такое значение β, при котором матрица B становится вырожденной, то есть имеет собственное значение равное нулю. Соответствующее уравнение Þ . Поскольку число Фробениуса есть максимальное из всех , постольку при возрастании β число есть наименьшее, при котором матрица B становится вырожденной. Следовательно, Þ .

Пример 1. Задана матрица межотраслевого баланса

и одно из собственных значений этой матрицы. Требуется найти: 1) собственные значения и собственные векторы матрицы, 2) число и вектор Фробениуса матрицы A, используя алгоритм Ланцоша, 3) запас продуктивности матрицы A.

Решение.

Рассмотрим матрицу межотраслевого баланса. Уравнение на собственные значения для этой матрицы имеет вид

(7)

После разложения определителя получим

(8)

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

Вычислив корни квадратного уравнения, получим , .

Для того, чтобы найти собственные векторы, выпишем уравнение (4) вычислив коэффициенты с точностью до 3 знака после запятой. Получим

(9)

Подставим по очереди собственные значения в систему ()

I.
(10)
Поскольку система (10) линейно зависима, одну из переменных можно выбрать произвольно. Выберем . Остальные переменные найдем из уравнений (10). Тогда собственный вектор, отвечающий выбранному собственному значению, есть .

II.
(11)
Решая систему (11) также, как и в предшествующем случае, выберем . Собственный вектор, отвечающий собственному значению , есть .

III.
(12)
Решая систему (12) найдем собственный вектор .

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

. (13)

Таким образом, межотраслевой баланс отраслей O 1, O 2, O 3 будет при соотношении валовых выпусков x 1: x 2: x 3 = 1.25: 1.50: 1. Доля конечного потребления в этом случае определяется формулой (13).

Таким образом, найдено число Фробениуса F = 0.923 и вектор Фробениуса матрицы A.

Найдем указанные величины, пользуясь алгоритмом Ланцоша. Выберем начальный вектор . Тогда , после нормировки получим вектор в первом приближении . Можно (но не обязательно) также вычислить нулевое приближение для числа Фробениуса . Для получения вектора следует вычислить вектор и нормировать его. Вычисления удобно свести в таблицу

Номер итерации i          
Вектор 0.570 0.610 0.549 0.568 0.632 0.526 0.568 0.647 0.507 0.569 0.658 0.494 0.569 0.665 0.483
Число 0.915 0.920 0.920 0.922 0.922

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

Зная число F, вычислим запас продуктивности матрицы:

Варианты заданий

Задание. Задана матрица межотраслевого баланса и одно из собственных значений этой матрицы. Требуется найти: 1) собственные значения и собственные векторы матрицы, 2) число и вектор Фробениуса матрицы A, используя алгоритм Ланцоша, 3) запас продуктивности матрицы A.





Дата публикования: 2015-11-01; Прочитано: 236 | Нарушение авторского права страницы | Мы поможем в написании вашей работы!



studopedia.org - Студопедия.Орг - 2014-2024 год. Студопедия не является автором материалов, которые размещены. Но предоставляет возможность бесплатного использования (0.011 с)...