САНКТ-ПЕТЕРБУРГСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ
КАФЕДРА МОДЕЛИРОВАНИЯ ЭКОНОМИЧЕСКИХ СИСТЕМ
Белоусова Мария Владимировна
Выпускная квалификационная работа бакалавра
Анализ функционирования региональной экономики
на основе модели межотраслевого баланса
Направление 010400
Прикладная математика и информатика
Научный руководитель,
доктор физ.-мат. наук,
профессор
Смирнов Н. В.
Санкт-Петербург
2016
Содержание
Введение
3
Обзор литературы
4
Постановка задачи
5
Глава 1. Построение статической модели
7
§ 1.1. База данных WIOD . . . . . . . . . . . . . . . . . . . . . . . . .
7
§ 1.2. Описание NIOT . . . . . . . . . . . . . . . . . . . . . . . . . . .
9
§ 1.3. Вывод статической модели . . . . . . . . . . . . . . . . . . . . . 10
§ 1.4. Агрегирование и редуцирование . . . . . . . . . . . . . . . . . . 12
Глава 2. Динамическая модель
14
§ 2.1. Классическая динамическая модель . . . . . . . . . . . . . . . . 15
§ 2.2. Расширенная динамическая модель МОБ . . . . . . . . . . . . . 16
§ 2.3. Переход от разностной системы к дифференциальной . . . . . 17
Глава 3. Прогноз основных экономических показателей
19
§ 3.1. Регрессионный подход к нахождению фондоемкостей . . . . . . 19
§ 3.2. Идентификация фондоемкостей . . . . . . . . . . . . . . . . . . 19
§ 3.3. Обзор регрессионных моделей . . . . . . . . . . . . . . . . . . . 20
§ 3.4. Вычисление прогноза. Оценка погрешности . . . . . . . . . . . 21
Выводы. Заключение
24
Список литературы
26
Приложение
28
2
Введение
Одной из существенных задач современной науки является совершенствование математических моделей, способных определять структуру национального хозяйства и прогнозировать основные макроэкономические тенденции.
Такие модели должны обеспечивать планирование и регулирование экономических показателей. Однако, со временем экономическая система становится
сложнее, что связано как с возрастанием разновидностей товаров и услуг, так
и с постоянным изменением самой структуры (эволюция технологии производства, трансформация производственных связей и общественных отношений). В связи с этим возникают трудности с выбором адекватных моделей и
инструментов для анализа и описания экономических процессов.
Наиболее известной является модель межотраслевого баланса (МОБ),
описанная В. В. Леонтьевым [1]. Она представляет собой систему взаимосвязей между различными секторами экономики либо экономическими системами. За ее разрабку Василий Леонтьев получил Нобелевскую премию в
1973 году. Модель отображает межотраслевые отношения, в которых произведенный товар или услуга одной индустрии является продуктом для другой.
Конечно, сбалансированность в процессе хозяйственной деятельности нечасто осуществляется, так как неизбежно влияние различных непредсказуемых
факторов (политических, стихийных, экологических и т.д.). Тем не менее, балансовый метод оказывается довольно эффективным.
Модель «затраты-выпуск» примечательна еще и тем, что почти целиком
описывает процесс воспроизводства. По этой причине данная модель представляет особый интерес для анализа конечного потребления, выпусков, промежуточного потребления, экспорта и импорта, добавленной стоимости, национального дохода, различных видов налогов и т.д.
Для использования модели межотраслевого баланса создаются специальные таблицы. Они составляются для многих стран мира, так как они
3
предоставляют довольно подробное отображение структуры промышленной
деятельности. Также таблицы «затраты-выпуск» дают возможность оценить
последствия от изменений в определенных отраслях экономики, влияния роста или падения цен на различные товары. Кроме того, важным достоинством
данных таблиц является то, что они позволяют выделить основополагающие
экономические отрасли.
Также эта модель - одна из тех, что позволяет находить специфические
состояния экономической системы (например, состояния равновесия).
В данной работе основное внимание уделено динамической модели межотраслевого баланса. Она обобщает статическую модель и строится на тех же
предположениях.
В анализе, проведенном в исследовании, используются таблицы
«затраты-выпуск» для экономики Германии, являющейся одной из самых
больших экономик мира и крупнейшей в Европе. Данный выбор обусловлен
тем, что экономическая деятельность этой страны примечательна тем, что
обладает высокоэффективными производствами разных направлений, основными из которых являются машиностроение, химическая промышленность,
электротехническая промышленность и многие другие.
Обзор литературы
Модель Леонтьева была изложена во многих источниках, представляется рациональным использовать информацию из первоисточника [1].
Было разработано несколько переходов к динамической модели. Один
из них осуществлен Н. И. Ведутой [2] методами экономической кибернетики. Переход к разностной и дифференциальной системам, суть которого заключается в использовании межотраслевых потоков капитальных вложений,
описан в монографии [3]. Алгоритм построения расширенной модели, когда
общественное потребление рассматривается как фазовая переменная, предло-
4
жен в [4, 5].
В настоящее время существует значительное количество ресурсов, публикующих информацию по экономике отдельных регионов в виде ежегодных
таблиц «затраты-выпуск». Например, «Федеральная служба государственной статистики» (Росстат) [6], статистическая служба Евросоюза (Евростат)
[7]. Особый интерес представляет «Мировая база данных «затраты-выпуск»
(World Input-Output Database, WIOD) [8], которая предоставляет доступ как
к всемирным таблицам, так и к таблицам по отдельным странам в различных модификациях, например, «National Input-Output Tables» (NIOT). Этот
ресурс является основным источником статистической информации в данной
дипломной работе.
Постановка задачи
Процесс исследования можно условно разделить на решение двух задач.
Первой из них является конкретизация основных балансовых соотношений и построение статической модели межотраслевого баланса для таблиц
NIOT экономики Германии для каждого года с 1995 по 2011. Иными словами,
составление матрицы относительных цен Rk [4], вектора выпуска Ik и вектора
потребления Yk , таких что:
Ik = Rk Ik + Yk ,
k = 1995, 2011,
и последующее получение системы
ek Iek
Iek = R
ek —
путем добавления ВВП в число фазовых переменных. Соответственно R
расширенная матрица модели.
Второй (и главной) задачей работы является получение динамических
модели МОБ в виде системы дифференциальных (разностных) уравнений от-
5
носительно выпусков и ВВП [4, 5]
˙ = DI(t) (Ik+1 = D1 Ik )
I(t)
с целью получения прогноза макроэкономических тенденций на будущий год.
За счет варьирования таких параметров, как фондоемкости отраслей, входящих в матрицу D (D1 ), необходимо выбрать наилучшую модель с точки
зрения функционала
b
min ||I(1) − I(1)||
(min ||I1 − Ib1 ||),
D
D1
b
где I(1) — реальный будущий выпуск, I(1)
— прогнозный выпуск,
b = eD I(0) (Ib1 = D1 I0 ).
I(1)
На последнем этапе провести анализ региональной экономики на основе
построенной модели.
6
Глава 1. Построение статической модели
В данной главе будет описана структура базы данных, на основе статистических данных которой, впоследствии проводится исследование. Также будут
охарактеризованы используемые таблицы, а затем произведен вывод основных балансовых соотношений для статической модели «затраты-выпуск».
§ 1.1. База данных WIOD
Ввиду ускоряющейся глобализации происходит переплетение экономик различных стран, перераспределение доходов между ними. С одной стороны,
этот процесс способствует возрастанию числа рабочих мест и экономическому росту. С другой стороны, изменение структуры международной торговли
может вызвать проблемы в распределении ресурсов. Для проведения качественного анализа этих процессов, необходимо работать с большим объемом
статистических данных, обновляемых каждый год. Именно таким источником информации является «Мировая база данных «затраты-выпуск»» (World
Input-Output Database, WIOD) [8].
Данный ресурс был впервые запущен в апреле 2012 года и финансируется Европейской комиссией. Он создавался в сотрудничестве с учеными из
Гронингенского университета, Нидерланды.
WIOD состоит из мировых таблиц различных модификаций. Эти таблицы содержатся в четырех крупных частях данного ресурса:
• мировые таблицы (World Tables),
• национальные таблицы (National Tables),
• социально-экономические счета (Socio-Economic Accounts),
• экологические счета (Environmental Accounts).
Первый раздел «World Tables» включает в себя:
7
1. Международные таблицы ресурсов и использования товаров и услуг с
потреблением внутренних и импортных товаров;
2. Мировые таблицы «затраты-выпуск»;
3. Межрегиональные таблицы «затраты-выпуск» для шести крупных регионов.
Особенностью мировых таблиц является тот факт, что по ним можно узнать
из какого зарубежного промышленного сектора производится импорт для
какой-либо страны, а также, как экспортная продукция данной страны используется в других странах.
Далее идет совокупность таблиц «National Tables», в которую входят:
1. Национальные таблицы использования товаров и услуг;
2. Национальные таблицы «затраты-выпуск».
Национальные таблицы позволяют провести анализ экономики отдельной страны в период с 1995 по 2011 годы. Они были построены для экономик 27 стран
Евросоюза и 13 экономик других стран, в сумме производящих около 90%
мировой продукции.
Следующий раздел «Socio-Economic Accounts» предоставляет детальные данные об использовании капитала и труда:
1. Промышленный выпуск, добавленная стоимость;
2. Основной капитал, инвестиции;
3. Заработная плата и занятость людей трех видов квалификации: низкой,
средней и высокой.
Последней составляющей базы WIOD является «Environmental Accounts».
Здесь содержатся такие разновидности информации:
1. Объем используемой секторами энергии;
2. Выбросы углекислого газа по секторам экономики;
8
3. Выбросы в атмосферу;
4. Использование земли и других природных материалов.
В данном исследовании основным источником статистической информации являются именно национальные таблицы — National Input-output tables
(NIOT), а вспомогательным — таблицы из раздела «Socio-Economic Accounts».
Рассмотрим структуру NIOT подробнее.
§ 1.2. Описание NIOT
В ежегодных таблицах NIOT представлены взаимоотношения 35-ти экономических секторов, все значения представлены в текущих ценах в американских
долларах (рис. 1).
Рис. 1. Схематическое представление таблицы NIOT
Таблицы «затраты-выпуск» являются симметричными в том смысле,
что в строках и столбцах представляются одинаковые названия экономических отраслей, только секторы, представленные в строках — производители,
а в столбцах — потребители продукции других секторов.
9
Элемент pij матрицы P размерности 35 × 35 представляет собой промежуточное потребление j-ым сектором услуг и товаров i-го сектора и задаются
соотношением pij = Pi aij Inj , где Pi — цена потребленной продукции, aij —
технологический коэффициент, а Inj — объем годового выпуска j-ого сектора
в натуральном выражении.
35
Матрица P m = {pm
ij }i,j=1 показывает промежуточное потребление им-
портных товаров. Вектор-cтрока V = (V1 , V2 , . . . , V35 ) характеризует добавленную
стоимость,
полученную
каждым
сектором.
Вектор-столбец
Y = (Y1 , Y2 , . . . , Y35 )T показывает конечное потребление продукции индустрий, он также включает в себя валовое накопление основного капитала
Cf и изменения товарно-материальных запасов. Далее идет столбец экспорта Ex = (Ex1 , Ex2 , . . . , Ex35 )T . Через T x обозначен налог на потребление
(сумма всех налогов, полученных не из добавленной стоимости секторов).
Элемент Ij 35-мерного вектора I обозначает выпуск продукции j-го сектора.
§ 1.3. Вывод статической модели
Согласно таблицам NIOT, выпуск продукции j-го сектора Ij может быть вычислен как
Ij =
n
X
(pij + pm
ij ) + Vj ,
∀j = 1, n,
(1)
∀i = 1, n.
(2)
i=1
а ежегодные продажи i-ой отрасли
Xi =
n
X
pij + Yi + Exi ,
j=1
Исходя из того, что в рассматриваемых таблицах выпуски совпадают с продажами Ii = Xi (ситуация динамического равновесия), просуммировав данные
равенства по i = 1, n, из (1) и (2) получим
n X
n
X
j=1
(pij +
pm
ij )
+ Vj =
i=1
n X
n
X
i=1
10
j=1
pij + Yi + Exi .
Отсюда
n
X
j=1
где Imi =
Vj =
n
X
(Yi + Exi − Imi ),
i=1
Pn
m
j=1 pij .
Далее введем понятие ВВП:
n
n
X
X
GDP =
Vj + T x =
(Yi + Exi − Imi ) + T x.
j=1
(3)
i=1
Введём теперь новые обозначения для дальнейшего вывода модели:
Iei = Ii − Imi ,
Yei = Yi + Exi − Imi ,
и перейдем к относительным ценам [5, 4]
pij
rij =
.
Iej
Тогда из (2) получим формулу для вычисления выпусков:
n
n
X
X
e
Ii = Ii − Imi =
pij + Yi + Exi − Imi =
rij Iej + Yei ,
j=1
i = 1, n,
(4)
j=1
или в матричном виде
(5)
Ie = RIe + Ye ,
где Ie = (Ie1 , Ie2 , . . . , Ien )T ,
Ye = (Ye1 , Ye2 , . . . , Yen )T ,
r r . . . r1n
11 12
r r . . . r
2n
21 22
R= .
.
.
.
..
. . ...
..
rn1 rn2 . . . rnn
Рассмотрим теперь модель, в которой ВВП является (n + 1)-ым компонентом фазового вектора:
Ie = (I1 − Im1 , I2 − Im2 , . . . , In − Imn , GDP )T ,
Будем использовать ряд относительных величин [5, 4]
Yei
e
, i = 1, n,
Y ri =
Ien+1
11
Ien+1 = GDP.
V rj =
Vj
, i = 1, n,
Iej
Tx
.
Ien+1
T xr =
Тогда (4) примет вид:
Iei =
n
X
rij Iej + Ye ri Ien+1 ,
i = 1, n,
j=1
Таким образом, имеем n уравнений с n + 1 неизвестной. Последнее уравнение
получим из формулы (3):
n
n
X
X
Ien+1 =
Vj + T x =
V rj Iej + T xrIen+1 .
j=1
j=1
В матричном виде полученная система из n + 1 уравнения относительно n
выпусков и ВВП будет выглядеть как
eI,
e
Ie = R
e=
R
R
Ye r
V r T xr
!
,
(6)
где V r = (V r1 , V r2 , . . . , V rn ), Y r = (Ye r1 , Ye r2 , . . . , Ye rn )T .
Модель (6), включающая ВВП как элемент фазового вектора, является
аналогом классической модели (5), не содержащей в себе ВВП.
e обладает нижеперечисленными свойствами:
Матрица R
1. Она является неотрицательной: reij > 0, i = 1, n + 1, j = 1, n + 1;
P
2. Суммы по столбцам дают единицу: n+1
eij = 1, j = 1, n + 1;
i=1 r
3. Всегда существует собственное число λ = 1. Оно же является максимальным по модулю.
e
Для удобства в дальнейшем переобозначим I = I.
§ 1.4. Агрегирование и редуцирование
В данном параграфе приведены два подхода к уменьшению размерности системы (6) [1]. Их использование является оправданным и уместным в зависимости от цели исследования.
12
Агрегирование. При анализе макроэкономических трендов может быть
нецелесообразно оперировать со значениями конкретных отраслей. Поэтому
логично объединять некоторые секторы в обобщенные агрегированные отрасли.
Рассмотрим построение агрегированной системы на примере с двумя
секторами:
s
X
Ia =
Ii ,
i=1
n+1
X
!T
Ii
.
i=s+1
Из (6) получим
Ps Ps
Ps Pn+1
r
e
I
eij Ij a
P
P
ij
j
j=1
j=s+1 r
i=1
i=1
n+1
s
a
I
+
I2 ,
I1a = i=1 j=1 reij Ij =
1
a
a
I
I
1
2
Pn+1 Ps
Pn+1 Pn+1
r
e
I
eij Ij a
P
P
ij
j
j=1
i=s+1
j=s+1 r
i=s+1
n+1
n+1
a
a
r
e
I
=
=
+
I
I
I2 .
2
1
j=1 ij j
i=s+1
I1a
I2a
В матричном виде
Ps
i=1
ea I a ,
Ia = R
Ps
eij Ij
j=1 r
I1a
a
e
R = Pn+1 P
i=s+1 sj=1 reij Ij
I1a
Ps
Pn+1
r
e
I
j=s+1 ij j
i=1
I2a
Pn+1 P
.
n+1
eij Ij
i=s+1
j=s+1 r
I2a
Редуцирование. Предположим теперь, что интерес представляют только n1 компонент вектора I, объединим их в I r1 , тогда I r2 — n2 -мерный вектор
выпусков оставшихся секторов (n1 + n2 = n + 1). Систему (6) можно переписать в виде
I r1 = R
e11 I r1 + R
e12 I r2 ,
I r1 = R
e21 I r1 + R
e22 I r2 .
e22 ), где En — единичная матрица размерности n2 × n2 ,
Если матрица (En2 − R
2
является
I r2
невырожденной, из второго уравнения сможем найти
e22 )−1 R
e21 I r1 . Подставив в первое уравнение, получим замкнутую
= (E − R
относительно выпусков I r1 систему
er I r1 ,
I r1 = R
er = R
e11 + R
e12 (E − R
e22 )−1 R
e21 .
R
13
e построенную по таблице NIOT
Пример. Рассмотрим систему I = RI,
для Германии за 2000 год.
Построим агрегированную систему, объединив в первую отрасль все секP
торы с 1 до 35 (всё производство): I1a = 35
i=1 Ii , во второй отрасли останется
сфера потребления I2a = I36 . Тогда получим
!
!
12
2,
876
·
10
0,
4146
0,
9200
ea =
Ia =
,R
.
1, 830 · 1012
0, 5854 0, 0800
Проиллюстрируем результат редуцирования для той же системы, оставив отрасль «Renting of M&Eq and Other Business Activities» и сферу потребления:
!
!
11
0,
1798
0,
1283
2,
863
·
10
er =
,R
.
I r1 =
1, 830 · 1012
0, 8202 0, 8717
ea и R
er обладают свойствами матрицы R,
e описанными
Отметим, что матрицы R
в предыдущем параграфе.
Глава 2. Динамическая модель
На основе модели межотраслевого баланса, описанной в главе 1, может быть
осуществлен переход к динамической модели. Её рассмотрение выглядит целесообразным ввиду того, что экономика представляет собой динамическую
структуру, показатели которой прямым образом зависят от действий в прошлом. Кроме того, использование такой модели позволяет прогнозировать
развитие экономики, что несомненно является актуальной задачей.
В § 2.1–2.2 представлены два подхода к построению динамической модели. В первом подходе изложен классический подход к построению разностной
и дифференциальной системы, во втором — метод построения расширенной
динамической модели.
В зависимости от задачи и методов её решения удобнее анализировать
дифференциальную, либо разностную систему. Они являются в некотором
14
смысле эквивалентными, от одной можно переходить к другой и обратно. В
§ 2.3 описан метод преобразования дифференциальной системы в разностную.
§ 2.1. Классическая динамическая модель
Данная модель основывается на статической модели (5), рассмотренной в прошлой главе:
I = RIe + Y,
Ключевой идеей является выделение коэффициентов капиталовложений i-го
сектора в j-ую отрасль ∆Φij из вектора потребления [3]:
Yi =
n
X
∆Φij + Yi0 .
j=1
Параметр ∆Φij характеризует количество продукции i-ой отрасли, направленной в j-ую в качестве производственных капитальных вложений в её основные фонды за рассматриваемый интервал времени.
Далее предполагается, что капиталовложения в j-ую отрасль пропорциональны приросту выпуска продукции предыдущего периода k − 1 до текущего периода k, ∆Ij = Ij,k − Ij,k−1 , с коэффициентом пропорциональности
f eij , называемым коэффициентом приростной фондоемкости:
∆Φij = f eij ∆Ij ,
i = 1, n,
j = 1, n.
Подставив это соотношение в (4), получим
Ii =
n
X
rij Ij +
j=1
n
X
f eij ∆Ij + Yi0 .
j=1
Иным образом это выражение можно переписать так:
n
n
X
X
0
Ii,k =
rij Ij,k +
f eij (Ij,k − Ij,k−1 ) + Yi,k
.
j=1
(7)
j=1
Введя обозначение F = {f e}ni,j=1 , можно переписать систему (7) в матричном
виде:
Ik = RIk + F (Ik − Ik−1 ) + Yk0 .
15
(8)
Система (8) — динамическая модель МОБ в виде разностных уравнений. Заменяя дискретные моменты времени k на непрерывное время t, а величину
˙
δIk на производную I(t),
получим динамическую систему МОБ в виде дифференциальных уравнений:
(9)
˙ + Y 0 (t),
I(t) = RI(t) + F I(t)
или
˙ = (E − RI)(t) − Y 0 (t).
F I(t)
§ 2.2. Расширенная динамическая модель МОБ
Согласно [4, 5], предположим, что изменение объема выпуска i-го сектора
I˙i (t) пропорционально инвестициям в эту отрасль Cpi (t), роль коэффициента
пропорциональности играет фондоемкость сектора F ei (t):
F ei (t)I˙i (t) = Cpi (t),
i = 1, n + 1.
Будем считать, что на инвестиции идет вся добавленная стоимость отрасли,
за исключением оплаты труда работников, т.е. Cpi = Vi −Wi , где Wi — фонды
зарплаты, или что то же Cpi = (V ri − W ri )Ii , Wi = Wi Ii . Пусть для сферы
потребления Cpn+1 = (T xr + Cf r)In+1 , где T xr и Cf r являются нормированными по ВВП значениями налога на потребление и валового накопления
основного капитала соответственно.
Тогда с учетом модели (6) получим систему
˙ = DI(t),
I(t)
где матрица
V r1 −W r1
F e1
M =
(10)
e
D = M R,
..
.
...
...
0
..
.
0
..
.
0
...
V rn −W rn
F en
0
0
...
0
T xr+Cf r
F en+1
16
.
Ввиду того, что таблицы NIOT ежегодно пересчитываются, представляется целесообразным ежегодно пересчитывать и матрицу D. Тогда получим
систему
˙ = Dk I(t), t = [k, k + 1),
I(t)
I(k) = Ik ,
k = 1995, 2011,
(11)
где матрица Dk и вектор Ik представляют собой матрицу D и вектор I, построенные по таблицам за k-ый год.
§ 2.3. Переход от разностной системы к дифференциальной
В § 2.1 был совершен переход от дискретной модели (8) к непрервыной системе (9). С математической точки зрения, данный переход является не совсем
корректным. Изложим в данном параграфе метод перехода от разностной
линейной системы уравнений к дифференциальной системе, при котором решения обоих систем будут совпадать при t = k [10, 11].
Рассмотрим разностное уравнение:
(12)
xk+1 = Axk + b.
Если матрица A не имеет собственных чисел, абсолютное значение которых
равняется единице, |λj | =
6 1, j = 1, n, то его решение в общем виде:
xk = Ak x0 + (A − E)−1 (Ak − E)b,
k = 0, 1, 2, . . .
Теперь введем в рассмотрение непрерывную систему уравнений:
ż(t) = Cz(t) + d.
(13)
При невырожденной матрице C её решением будет:
z(t) = eCt z0 + C −1 (eCt − E)d.
В случае одинаковых начальных данных x0 = z0 , решения систем (12) и (13)
будут совпадать при целочисленных значениях аргумента, если их коэффициенты связаны соотношениями:
17
1) eC = A,
2) C −1 d = (A − E)−1 b,
откуда
1) C = ln A,
2) d = (A − E)−1 ln A · b.
18
Глава 3. Прогноз основных экономических показателей
В данной главе вычисляется прогноз выпусков секторов и ВВП на следующий
год на основе динамической модели (11). Особое внимание уделено нахождению фондоемкостей F ei (t), i = 1, n + 1. Помимо способа вычисления этого
параметра по определению, рассмотрен ряд регрессионных моделей. Высчитываются оценки отклонения прогноза.
§ 3.1. Регрессионный подход к нахождению фондоемкостей
В расширенной модели МОБ, описанной в § 2.2, фондоемкости являются коэффициентами пропорциональности между ростом выпуска и объемом инвестиций и по определению равны [5]
F ei (k) =
Cpi (k)
,
Ii (k) − Ii (k − 1)
i = 1, n + 1,
k = 1995, 2011.
(14)
Однако в данной работе с целью достижения лучшей точности прогноза рассмотрим регрессионные способы определения этих параметров. Для поиска
взаимосвязей с различными показателями и составления уравнений зависимости от них, необходимо иметь значения фондоемкостей за прошлые временные отрезки. Поэтому в следующем параграфе приведен способ определения
их значений в прошлом по историческим данным.
§ 3.2. Идентификация фондоемкостей
Согласно (11) с учетом формулы Коши, прогноз расширенного вектора выпусков на год вперед может быть вычислен как
Ibk+1 = eDk Ik ,
k = 1995, 2011,
здесь Dk , Ik определены по информации из таблицы NIOT за k-ый год. В
идеальном случае данный прогноз должен совпадать со значением Ik+1 , полученным из таблиц за следующий год. Поэтому на исторических данных
19
может быть поставлена задача определения матрицы Dk таким образом, чтобы Ibk+1 и Ik+1 совпадали, либо были в какой-то мере близки. В данной работе
варьирование Dk осуществляется за счет вектора фондоемкостей. Тогда получаем задачу из (n + 1)-го уравнения с (n + 1)-ой неизвестной для каждого
доступного года:
k
k
Ik+1 = eD(F e1 ,...,F en+1 ) Ik ,
(15)
k = 1995, 2010.
Эта система была решена в среде Matlab с использованием функции fsolve.
Таким образом, был найден временной многомерный ряд {F ek }2010
k=1995 , который в дальнейшем используется в качестве учебной выборки для составления
регрессионных моделей.
В таблице 2 приведены идентифицированные значения фондоемкостей
для секторов № 3 «Food, Beverages and Tobacco», № 11 «Other Non-Metallic
Mineral», №18 «Construction», №27 «Post and Telecommunications».
XXX
год
XX
№ индустрии XXXX
X
1995
1996
1997
1998
1999
2000
2001
2002
3
-0,6472
-1,1523
0,0678
-0,7762
-0,6215
-0,3664
-0,1805
0,6145
11
-0,6538
-0,9889
-0,1606
-0,1417
-0,8977
-0,6102
0,0362
0,7272
18
-0,5175
-1,1232
-0,3250
-0,2370
-1,0917
-0,5319
-0,1349
0,7350
27
-0,2014
-0,2831
0,1211
0,0790
-0,2751
0,1053
0,3386
0,5629
год
XXX
XXX
№ индустрии
X
2003
2004
2005
2006
2007
2008
2009
2010
3
0,1809
-0,2591
0,0140
0,3360
0,0416
-1,1540
0,0315
0,2539
11
0,6110
-0,1398
0,4030
0,5597
0,2117
-1,2822
-0,0054
0,5957
18
0,1511
-0,2118
0,2775
0,4785
0,4691
-0,2088
-0,2166
0,3090
27
0,2350
-0,0688
-0,0381
0,0385
0,2379
-0,4549
-0,0273
0,1842
XXX
XXX
XXX
Таблица 1. Таблица фондоемкостей
§ 3.3. Обзор регрессионных моделей
Исследуем функциональную зависимость фондоемкостей от ряда параметров.
Будем использовать модель линейной регрессии, где параметры высчитываются методом наименьших квадратов. В качестве объясняемых переменных
1
рассмотрим обратные величины
, i = 1, n + 1, т.к. именно в таком виде
F ei
20
они входят в матрицу M системы (10).
Среди всех моделей выделим следующие три:
1
(1)
(1)
= ai k + b i ,
k
F ei
i = 1, n + 1,
1
(2)
(2)
= ai inf lk + bi , i = 1, n + 1,
k
F ei
где inf lk обозначает уровень инфляции в год k,
1
(3)
(3)
= ai (Ii (k) − Ii (k − 1)) + bi ,
k
F ei
i = 1, n + 1.
(16)
(17)
(18)
§ 3.4. Вычисление прогноза. Оценка погрешности
Предположим, что имеется информация в виде таблиц NIOT за 11 последних
лет: с (s − 10)-го по s-ый. Тогда прогноз вектора выпусков и ВВП на год
вперед Ibs+1 может быть вычислен по следующей схеме:
1. Идентификация вектора фондоемкостей за предыдущие 10 лет {F ek }s−1
k=s−10
(см. § 3.2).
ces из уравнения регрессии (см. § 3.3),
2. Нахождение текущего значения F
либо по определению по формуле (14).
s
c
3. Вычисление прогноза по формуле Ibs+1 = eDs (F e ) Is .
4. Оценка погрешности rs+1
||Is+1 − Ibs+1 ||
=
100%.
||Is+1 ||
Данный алгоритм был реализован в виде программного комплекса в
среде Matlab.
Замечание 1. Учитывая объем имеющейся информации, данный процесс
может быть проделан для s = 2005, 2010.
Замечание 2. Глубина учебной выборки, равная 10, была определена
эмпирически по историческим данным.
1
Обозначим за rs+1
относительное отклонение годового прогноза от исl
торических данных при вычислении фондоемкостей по определению, за rs+1
,
21
где l = 2, 3, 4 — погрешности при вычислении фондоемкостей по формулам
(16), (17), (18) соответственно.
Найдем относительные ошибки аппроксимации по формуле
2010
1 X l
Al =
r .
6 s=2005 s+1
Полученные результаты проиллюстрированы в таблице 2.
A1
A2
A3
A4
8, 25%
8, 31%
6, 21%
6, 46%
Таблица 2. Относительные ошибки аппроксимации
В таблице 3 продемонстрированы результаты прогноза на 2010 год для
некоторых секторов. I2009 , I2010 — реальные выпуски за 2009 и 2010 годы,
i
, i = 1, 4 — прогнозы при вычислении фондоемкостей по формулам (14),
Ib2010
(16)–(18) соответственно.
№ индустрии
I2009
I2010
1
Ib2010
2
Ib2010
3
Ib2010
4
Ib2010
1
4, 86 · 104
4, 81 · 104
4, 17 · 104
5, 03 · 104
4, 49 · 104
4, 47 · 104
6
1, 70 · 104
1, 74 · 104
1, 25 · 104
1, 95 · 104
1, 53 · 104
1, 43 · 104
11
3, 22 · 104
3, 22 · 104
2, 42 · 104
3, 52 · 104
2, 97 · 104
2, 66 · 104
16
3, 31 · 104
3, 28 · 104
2, 43 · 104
3, 71 · 104
3, 11 · 104
2, 47 · 104
21
1, 51 · 105
1, 50 · 105
1, 32 · 105
1, 63 · 105
1, 49 · 105
1, 47 · 105
26
9, 68 · 104
10, 11 · 104
8, 31 · 104
11, 1 · 104
9, 60 · 104
9, 07 · 104
31
2, 09 · 105
2, 01 · 105
2, 07 · 105
2, 33 · 105
2, 01 · 105
2, 05 · 105
36
2, 79 · 106
2, 75 · 106
2, 53 · 106
3, 09 · 106
2, 72 · 106
2, 72 · 106
Таблица 3. Результаты прогнозирования
На рис. 2 представлен годовой прогноз ВВП Германии с 2007 по 2011
года. Реальное значение в ценах 1995 года сравнивается прогнозом, полученным путем интегрирования системы (11), где фондоемкости высчитываются
по определению (модель 1), либо по формулам (16)–(18) (модели 2–4 соответственно).
22
Рис. 2. Прогноз ВВП
По полученным результатам можно сделать вывод, что среди рассмотренных регрессионных моделей на данной выборке лучший результат показала модель, предполагающая зависимость от уровня инфляции.
23
Выводы. Заключение
К результатам работы можно отнести следующее:
1. Анализ базы данных WIOD.
2. Модификация статической и динамической модели МОБ под таблицы
NIOT.
3. Исследование методов сведения линейной разностной системы к системе
дифференциальных уравнений.
4. Построение регрессионных моделей зависимости фондоемкостей от различных показателей.
5. Вычисление прогноза выпусков и ВВП, оценка погрешности.
6. Реализация приведенных алгоритмов в среде Matlab.
Результаты проделанной работы были изложены на нескольких научных конференциях с последующей публикацией:
• Белоусова М. В., Дорофеев Б. В. «Оценка корректности перехода от линейных разностных уравнений к дифференциальным». Конференция «Процессы управления и устойчивость — 2015».
• M. V. Belousova, B. V. Dorofeev. «The method of finding the coefficients of a
dynamic input-output model». «Stability and Control Processes» in Memory
of V.I. Zubov (SCP), 2015 International Conference.
• Белоусова М.В., Попков А.С. «Построение динамической модели МОБ
на основе WIOD». Конференция «Процессы управления и устойчивость
— 2016».
В работе был изложен способ построения динамической модели межотраслевого баланса. Это удобный инструмент для анализа функционирования
экономики региона и получения прогноза различных показателей. Метод был
24
проиллюстрирован на примере экономики Германии в 1995–2011 гг. С целью
уменьшения ошибки прогноза были рассмотрены различные варианты определения фондоемкостей.
Продолжение данной работы может быть осуществлено засчет поиска
экономических показателей, наиболее коррелирующих с фондоемкостями, и
построения уравнений, предполагающих зависимость от этих показателей.
25
Список литературы
1. Леонтьев В. В. Межотраслевая экономика / пер. с англ., автор предисл. и
науч. ред. А. Г. Гранберг. М.: Экономика, 1997. 479 c.
2. Ведута Н. И. Социально эффективная экономика. Под общей ред. докт.
экон. наук. Ведута Е. Н. — М.: РЭА, 1999. — 254 с.
3. Федосеев В. В., Гармаш Д. М. Экономико-математические методы и прикладные модели / под ред. В. В. Федосеева. М.: ЮНИТИ, 1999. 391 с.
4. Пересада В. П. Управление динамикой развития экономики на базе межотраслевого баланса. СПб.: Политехника-сервис. 2010. 169 c.
5. Пересада В. П., Смирнов Н. В., Смирнова Т. Е. Управление развитием
многопродуктовой экономики на основе динамической модели «затратывыпуск»// Вестник Санкт-Петербургского университета. Серия 10. Прикладная математика, информатика, процессы управления, 2014. — № 4. —
P. 121–134.
6. Росстат [Электронный ресурс]: URL:http://www.gks.ru/ (дата обращения: 18.03.16)
7. Евростат [Электронный ресурс]: URL:http://ec.europa.eu/eurostat/
(дата обращения: 18.03.16)
8. WIOD [Электронный ресурс]: URL:http://www.wiod.org/ (дата обращения: 18.03.16)
9. Попков А. С. Идентификация динамической модели межотраслевого баланса для экономики России и оптимальное распределение инвестиций
на ее основе // Процессы управления и устойчивость. 2015. Т. 2 (18). С.
696–701.
26
10. Белоусова М. В., Дорофеев Б. В. Оценка корректности перехода от линейных разностных уравнений к дифференциальным // Процессы управления и устойчивость. 2015. T. 2. № 1. С. 541–544.
11. Belousova M. V, Dorofeev B. V. The method of finding the coeffcients of a
dynamic input-output model // 2015 International Conference on "Stability
and Control Processes" in Memory of V. I. Zubov, SCP 2015 – Proceedings,
2015. P. 455–456.
27
Приложение
Здесь приведем коды скриптов и функций, написанных в процессе выполнения данной работы. Все они выполнены в среде Matlab на встроенном языке.
Скрипт mob.m
ek и вектор Ik , k = 1995, 2011
Программа записывающая матрицы Mk , R
в переменную tb типа cell. В процессе выполнения скрипта подгружаются
Excel-таблицы Deu.xlsx, W.xlsx, Exr.xlsx, содержащие таблицы NIOT, информацию по зарплатам в отрасли за год, курс евро/доллар к концу каждого
года (требуется для пересчета зарплат из евро в доллары) соответственно.
1
l s = xlsread ('Deu.xlsx' ) ;
2
Comp = x l s r e a d ( 'W.xlsx ' , 'DATA' , 'E8211 : U8245' ) ;
3
Exr = x l s r e a d ( 'EXR.xlsx' , 'EXR' , 'C14 : S14' ) ;
4
for i = 1: length ( l s )
5
tb1 = x l s r e a d ( ' D e u . x l s x ' , num2str ( l s ( i ) ) ) ;
6
tb { 1 , i } = tb1 ( 5 : e n d , 5 : e n d ) ;
7
tb { 2 , i } = Exr ( i ) ;
8
tb { 3 , i } = Comp ( : , i ) * Exr ( i ) ;
tb { 3 , i } ( e n d ) = tb1 ( 8 2 , 3 9 ) ;
9
10
end
11
for j = 1: length ( l s )
12
Tb = tb { 1 , j } ;
13
I 1 = Tb ( 7 8 , 1 : 3 5 ) ' ;
14
A = Tb ( 1 : 3 5 , 1 : 3 5 ) ;
15
Im1 = Tb ( 3 6 : 7 0 , 1 : 3 5 ) ;
16
Im = sum ( Im1 , 1 ) ' ;
17
V1 = Tb ( 7 2 : 7 7 , 1 : 3 5 ) ;
18
V = sum (V1 , 1 ) ;
19
Y1 = Tb ( 1 : 3 5 , 3 6 : 4 0 ) ;
20
Y = sum (Y1 , 2 ) ;
21
Ex = Tb ( 1 : 3 5 , 4 1 ) ;
22
Tx1 = Tb ( 7 2 : 7 7 , 3 6 : 4 0 ) ;
23
Tx = sum ( sum ( Tx1 ) ) ;
24
Gdp1 = sum (V)+Tx ;
25
Gdp1 = sum (Y+Ex - Im)+Tx ;
26
A0 = [ A,Y+Ex - Im ; V, Tx ] ;
27
I = I1 - Im ;
28
I 0 = [ I ; Gdp1 ] ;
29
n = length ( I ) ;
30
f o r i = 1 : n+1
R0 ( : , i ) = A0 ( : , i ) / I 0 ( i ) ;
31
32
end
33
tb { 4 , j } = R0 ;
34
Rs = sum ( R0 ( 1 : n , 1 : n ) ,1) '+ tb { 3 , j } . / I ;
35
tb { 5 , j } = 1 - Rs ;
28
36
tb { 5 , j } ( 3 6 , 1 ) = (Tx + sum (Tb ( 1 : 3 5 , 3 9 ) ) ) / Gdp1 ;
37
tb { 6 , j } = d i a g ( tb { 5 , j } ) ;
tb { 7 , j } = I 0 ;
38
39
end
40
s a v e ( ' tb .m at ' , ' tb ' ) ;
Функция reduction.m
ek и вектора Ik (см. § 1.4). Входным параметРедукция матриц Mk , R
ром является вектор, состоящий из номеров отраслей, которые должны быть
оставлены в процессе редукции.
1
f u n c t i o n tb1 = r e d u c t i o n ( n1 )
2
l o a d ( ' tb .m at ' )
3
N = l e n g t h ( tb ) ;
4
n = l e n g t h ( tb { 4 , 1 } ) ;
5
n2 = s e t d i f f ( 1 : n , n1 ) ;
6
f o r i = 1 :N
7
R = tb { 4 , i } ;
8
M = tb { 6 , i } ;
I = tb { 8 , i } ;
9
I v = I ( n1 , 1 ) ;
10
11
Rv = R( n1 , n1 ) + R( n1 , n2 ) * ( ( e y e ( l e n g t h ( n2 ) ) -R( n2 , n2 ) ) \R( n2 , n1 ) ) ;
12
Mv = M( n1 , n1 ) ;
13
tb1 { 1 , i } = Rv ;
14
tb1 { 2 , i } = Mv;
tb1 { 3 , i } = I v ;
15
16
end
Скрипт ident.m
Идентификация вектора фондоемкостей F ek путем решения нелинейной
системы (15) с помощью функции fsolve.
1
clear
2
load ( 'tb1.mat')
all ;
3
F1 = @( Fe ) [ Fe ( 1 ) 0 ; 0 Fe ( 2 ) ] ;
4
M1 = tb1 { 2 , 1 } ;
5
R1 = tb1 { 1 , 1 } ;
6
I 0 = tb1 { 3 , 1 } ;
7
I 1 = tb1 { 3 , 2 } ;
8
F = @( Fe ) I1 - expm ( F1 ( Fe ) *M1) * I 0 ;
9
10
F10 = o n e s ( l e n g t h ( n1 ) , 1 ) ;
Fe = f s o l v e (F , F10 ) ;
Функция mnk.m
Функция, реализующая метод наименьших квадратов. На вход поступа29
ют вектор наблюдений независимой переменных x, вектор наблюдений объясняемой y, вектор x1, для которого будет высчитываться соответствующие
значения зависимой переменной. Далее вычисляются коэффициенты a1 , a2 ,
после чего находится выходное значение y1 = a1 x1 + a2 .
1
f u n c t i o n y1 = mnk( x , y , x1 )
2
A = [ sum ( x . ^2) sum ( x ) ; sum ( x ) l e n g t h ( x ) ] ;
3
b = [ sum ( x . * y ) ; sum ( y ) ] ;
4
a = A\b ;
5
y1 = [ x1 o n e s ( l e n g t h ( x1 ) , 1 ) ] * a ;
Скрипт prognoz.m
Вычисление прогноза вектора I на год вперед. Здесь представлен код
нахождения прогноза и его отклонения от реальных цифр при определении
фондоемкостей по формуле (17).
1
clear
2
l o a d ( ' tb .m at ' ) ;
all ;
3
load (' i n f l . m a t ' ) ;
4
load ('Inf.mat' ) ;
5
l1 = 10;
6
f o r j = 1:17 - l1 -1
m1 = j ;
7
8
[ T, Y, tb1 ] = i d e n t 0 ( n1 , m1 , l 1 ) ;
9
f o r i = 1 : l e n g t h ( n1 )
Y1( i ) = mnk( i n f l (T) ,Y ( : , i ) , ( i n f l (T ( e n d ) + 1 ) ) ) ;
10
11
end
12
I 1 = expm ( d i a g (Y1) * tb1 { 2 ,m1+l 1 } * tb1 { 1 ,m1+l 1 }) * tb1 { 3 ,m1+l 1 } ;
r 3 ( j ) = norm ( I1 - tb1 { 3 ,m1+l 1 +1})/norm ( tb1 { 3 ,m1+l 1 +1});
13
I1_3 ( : , j )= I 1 ;
14
15
end
16
R3 = mean ( r 3 ) ;
30
Отзывы:
Авторизуйтесь, чтобы оставить отзыв