О развитии теории параметрического регулирования эволюции национальной экономики на базе

advertisement
О развитии теории параметрического регулирования эволюции национальной экономики на базе
неавтономных динамических систем
Ашимов А.А., Султанов Б.Т., Боровский Ю.В., Боровский Н.Ю., Айсакова В.А., Алшанов Р.А.
В работе изложены некоторые результаты развития теории параметрического регулирования на классы
математических моделей, представляемых неавтономными непрерывными и дискретными динамическими
системами и, в частности, на подкласс вычислимых моделей общего равновесия (CGE модели).
Представлены теоремы о существовании решений некоторых задач вариационного исчисления для
указанных классов моделей и о непрерывной зависимости от неуправляемых функций оптимальных
значений критериев указанных задач. На примере CGE модели отраслей экономики проиллюстрирована
эффективность применения предложенного метода параметрической идентификации большеразмерных
математических моделей. С помощью вычислительных экспериментов проверены некоторые положения
макроэкономической теории о причинах возникновения циклических процессов в экономике. На базе CGE
модели отраслей экономики продемонстрирована эффективность подхода теории параметрического
регулирования для проведения ациклической государственной экономической политики.
1.Введение
Как известно, существует широкое согласие среди макроэкономистов по применению
математических моделей для макроэкономического анализа [1].
В [2], [3] на уровне автономных динамических систем предложена теория параметрического
регулирования для оценки рациональных значений инструментов в сфере государственной экономической
политики.
На практике открытая национальная экономика эволюционирует в условиях влияния внешних и
внутренних факторов, представляемых с помощью векторных функций времени, что обуславливает
необходимость развития теории параметрического регулирования на классы математических моделей,
задаваемых неавтономными динамическими системами.
В данной работе содержатся некоторые результаты развития теории параметрического
регулирования на классы неавтономных непрерывных, дискретных динамических систем и иллюстрации
развитых положений на примере дискретной неавтономной вычислимой модели общего равновесия
отраслей экономики.
2. Элементы теории параметрического регулирования на базе неавтономных динамических
систем
В данном разделе приводятся результаты по развитию теории параметрического регулирования в
рамках двух ее компонентов: методов оценки оптимальных значений экономических инструментов и
исследование влияний неуправляемых факторов на результаты решения задачи выбора оптимальных
значений экономических инструментов.
Описываемое ниже развитие теории параметрического регулирования является обобщением ранее
полученных нами в [3] положений для автономных динамических систем на случай неавтономных
динамических систем.
2.1 Оптимальные законы параметрического регулирования непрерывной динамической системы
Рассматривается непрерывная управляемая система в следующем виде:
𝑥̇ (𝑡) = 𝑓(𝑥(𝑡), 𝑢(𝑡), 𝑎(𝑡)), 𝑡 ∈ [0, 𝑇],
𝑥(0) = 𝑥0
(1)
(2)
где t
– время,
𝑥 = 𝑥(𝑡) = (𝑥 1 (𝑡), … , 𝑥 𝑚 (𝑡))– – вектор-функция состояния системы,
𝑢 = 𝑢(𝑡) = (𝑢1 (𝑡), … , 𝑢𝑞 (𝑡)) – вектор-функция управления,
𝑎 = 𝑎(𝑡) = (𝑎1 (𝑡), … , 𝑎 𝑠 (𝑡)) – известная вектор-функция,
𝑥0 = (𝑥01 , … , 𝑥0𝑚 ) – начальное состояние системы, известный вектор,
f
– известная вектор-функция своих аргументов.
Метод выбора оптимальных значений экономических инструментов связан со следующей моделью,
представляемой
1
- максимизацией критерия оптимальности
𝑇
max ∫0 𝐹(𝑡, 𝑥(𝑡)) 𝑑𝑡,
(3)
𝑢
где F – известная функция.
- фазовыми ограничения на решения системы (1)-(2) вида
𝑥(𝑡) ∈ 𝑋(𝑡), 𝑡 ∈ [0, 𝑇],
(4)
где 𝑋(𝑡)– заданное множество и
- явными ограничениями на управление:
𝑢(𝑡) ∈ 𝑈(𝑡), 𝑡 ∈ [0, 𝑇],
(5)
где 𝑈(𝑡) – заданное множество.
На базе соотношений (1)-(5) получаем следующую задачу, называемую задачей вариационного
исчисления по синтезу оптимальных законов параметрического регулирования для непрерывной системы.
Задача 1. При известной функции а найти измеримое управление u, удовлетворяющее условию (5),
чтобы соответствующее ему решение динамической системы (1), (2) удовлетворяло условию (4) и
доставляло максимум(минимум) функционалу (3).
Определим для фиксированных 𝑡 ∈ [0, 𝑇] и 𝑥 ∈ 𝑋(𝑡) множество Γ𝑡,𝑥 = {𝑓(𝑥, 𝑤, 𝑎(𝑡))|𝑤 ∈ 𝑈(𝑡)} в 𝑅𝑚 .
Обозначим через 𝑉𝑎 множество допустимых пар «состояние – управление» рассматриваемой системы при
заданной известной функции а, т.е. таких пар вектор-функций (𝑥, 𝑢), которые удовлетворяют соотношениям
(1), (2), (4), (5). Пусть В есть единичный шар с центром в начале координат в 𝑅𝑚 , X есть замыкание
множества ⋃𝑡∈[0,𝑇] 𝑋(𝑡), U - замыкание множества ⋃𝑡∈[0,𝑇] 𝑈(𝑡). Пусть множество возможных значений
функции a принадлежит некоторому множеству A.
Следующие две теоремы следуют из предложения 4.2 главы 8 монографии [4]. Доказательство
предложения 4.2 основано на теореме 2.1 о полунепрерывности снизу из главы 8 [4].
Теорема 1. Пусть функция а непрерывна на отрезке [0, 𝑇], U есть компакт в 𝑅𝑞 , функция f непрерывна
в 𝑋 × 𝑈 × 𝐴, и для любого 𝜌 ≥ 0 существует такое 𝜎 ≥ 0, что справедливо неравенство
|𝑓(𝑥, 𝑢, 𝑎(𝑡)) − 𝑓(𝑥′, 𝑢, 𝑎(𝑡))| ≤ 𝜎|𝑥 − 𝑥 ′ | ∀𝑡 ∈ [0, 𝑇] ,
𝑥, 𝑥 ′ ∈ 𝜌𝐵,
𝑢∈𝑈
и существует такая константа 𝜂 ≥ 0, что справедливо неравенство
|𝑥𝑓(𝑥, 𝑢, 𝑎(𝑡))| ≤ 𝜂(1 + |𝑥|2 ) ∀𝑡 ∈ [0, 𝑇],
𝑥 ∈ 𝑅𝑚 ,
𝑢∈𝑈
Пусть X – компакт, и функция F непрерывна на [0, 𝑇] × 𝑋. Пусть, кроме того, отображения 𝑡 → 𝑋(𝑡), 𝑡 →
𝑈(𝑡) являются непрерывными для 𝑡 ∈ [0, 𝑇] в следующем смысле: если справедливы включения 𝑥𝑘 ∈ 𝑋(𝑡𝑘 ),
𝑢𝑘 ∈ 𝑈(𝑡𝑘 ), где 𝑡𝑘 ∈ [0, 𝑇], 𝑘 = 1, 2, … и имеются сходимости 𝑡𝑘 → 𝑡, 𝑥𝑘 → 𝑥, 𝑢𝑘 → 𝑢, то справедливы
включения 𝑥 ∈ 𝑋(𝑡), 𝑢 ∈ 𝑈(𝑡). Тогда в случае непустоты множества 𝑉𝑎 и выпуклости множества 𝛤𝑡,𝑥 для
всех 𝑡 ∈ [0, 𝑇] , 𝑥 ∈ 𝑋(𝑡) задача 1 имеет решение.
Будем рассматривать неуправляемые функции а в (1) в виде элементов пространства непрерывных
вектор-функций Θ = (𝐶[0, 𝑇])𝑠 . Следующая теорема устанавливает достаточные условия непрерывной
зависимости оптимальных значений критерия задачи 1 от неуправляемых функций.
Теорема 2. Предположим, что при выполнении условий теоремы 1 в некоторой окрестности точки 𝑎0
в Θ, функция f непрерывна по второму аргументу и удовлетворяет условию Липшица по первому и
третьему аргументам на 𝑋 × 𝐴 равномерно по второму аргументу. Тогда оптимальное значение критерия
для задачи 1 непрерывно зависит от неуправляемой функции в точке 𝑎0 .
2.2 Оптимальные законы параметрического регулирования дискретной динамической системы
Рассматривается дискретная управляемая система в следующем виде:
𝑥(𝑡 + 1) = 𝑓(𝑥(𝑡), 𝑢(𝑡), 𝑎(𝑡)), 𝑡 = 0, 1, … 𝑛 − 1;
𝑥(0) = 𝑥0 .
2
(6)
(7)
Здесь t – время, принимающее неотрицательные целочисленные значения;
𝑥 = 𝑥(𝑡) = (𝑥 1 (𝑡), … , 𝑥 𝑚 (𝑡))– вектор-функция состояния системы дискретного аргумента;
𝑢 = 𝑢(𝑡) = (𝑢1 (𝑡), … , 𝑢𝑞 (𝑡)) – управление, вектор-функция дискретного аргумента;
𝑎 = 𝑎(𝑡) = (𝑎1 (𝑡), … , 𝑎 𝑠 (𝑡)) – известная вектор-функция дискретного аргумента;
𝑥0 = (𝑥01 , … , 𝑥0𝑚 ) – начальное состояние системы, известный вектор;
𝑓 – известная вектор-функция своих аргументов.
Метод выбора оптимальных значений экономических инструментов связан со следующей моделью,
представляемой
- максимизацией критерия оптимальности
max ∑𝑛𝑡=1 𝐹[𝑡, 𝑥(𝑡)]
(8)
𝑥(𝑡) ∈ 𝑋(𝑡), 𝑡 = 1, … 𝑛,
(9)
𝑢
где F – известная функция;
- фазовыми ограничения на решения системы вида
где 𝑋(𝑡) заданное множество;
- явными ограничениями на управление:
𝑢(𝑡) ∈ 𝑈(𝑡), 𝑡 = 0, … 𝑛 − 1,
(10)
где 𝑈(𝑡) – заданное множество.
На базе соотношений (6)-(10) получаем следующую вариационную задачу, называемую задачей
вариационного исчисления по синтезу оптимальных законов параметрического регулирования для
дискретной системы.
Задача 2. При известной функции а найти управление u, удовлетворяющее условию (10), чтобы
соответствующее ему решение динамической системы (6), (7) удовлетворяло условию (9) и доставляло
максимум (минимум) функционалу (8).
Обозначим через 𝑉𝑎 множество допустимых пар «состояние – управление» рассматриваемой системы
при заданной известной функции а, т.е. таких пар вектор-функций (𝑥, 𝑢), которые удовлетворяют
соотношениям (6), (7), (9), (10). Введем обозначения: 𝑋 = ⋃𝑛𝑡=1 𝑋(𝑡), 𝑈 = ⋃𝑛−1
𝑡=0 𝑈(𝑡).
Справедливы следующие две теоремы, доказательства которых основаны на использовании свойств
непрерывных функций и, в частности, на использовании свойств функций непрерывных на компакте.
Теорема 3. Предположим, что при известной функции а множество 𝑉𝑎 не пусто, множества 𝑋(𝑡) и
𝑈(𝑡 − 1) замкнуты и ограничены для всех 𝑡 = 1, … 𝑛, функция f непрерывна по первым двум аргументам на
множестве 𝑋 × 𝑈, а функция F непрерывна по второму аргументу на множестве Х. Тогда задача 2 имеет
решение.
Будем рассматривать неуправляемые функции а в (6) в виде элементов евклидова пространства 𝑅 𝑠𝑚 .
Теорема 4. Пусть при выполнении условий теоремы 3 для любых значениях 𝑎 ∈ 𝐴 (где 𝐴 – некоторое
открытое множество в евклидовом пространстве 𝑅 𝑠 ) функция f непрерывна по третьему аргументу в A и
удовлетворяет условию Липшица по первому аргументу в X равномерно по второму и третьему
аргументам в 𝑈 × 𝐴. Тогда оптимальное значение критерия для задачи 2 непрерывно зависит от
неуправляемой функции 𝑎 принимающей значения в A.
3. Вычислительные эксперименты с вычислимой моделью общего равновесия
Эффективность развитых положений теории параметрического регулирования иллюстрируются
ниже на подклассе вычислимых моделей общего равновесия (CGE моделях). Описание в общем виде
автономных CGE моделей и подход к параметрической идентификации таких моделей, изложенный нами в
[3] в этом разделе развивается на случай неавтономных CGE моделей.
3.1 Представление CGE модели
Неавтономная CGE модель в общем виде представляется с помощью следующей системы
соотношений.
1) Подсистема разностных уравнений, связывающая значения эндогенных переменных для двух
последовательных лет:
𝑥1 (𝑡 + 1) = 𝑓1 (𝑥1 (𝑡), 𝑥2 (𝑡), 𝑥3 (𝑡), 𝑢(𝑡), 𝑎(𝑡)),
3
(11)
Здесь 𝑡 = 0, 1, … , 𝑛 − 1 – номер года, дискретное время; 𝑥(𝑡) = (𝑥1 (𝑡), 𝑥2 (𝑡), 𝑥3 (𝑡)) ∈ 𝑅𝑚 – вектор
эндогенных переменных системы;
𝑥𝑖 (𝑡) ∈ 𝑋𝑖 (𝑡) ⊂ 𝑅𝑚𝑖 , 𝑖 = 1, 2, 3.
(12)
Здесь переменные
𝑥1 (𝑡) включают в себя значения основных фондов секторов-производителей, остатки средств
агентов на счетах в банках и др.;
𝑥2 (𝑡) включают в себя значения спроса и предложения агентов на различных рынках и др.,
𝑥3 (𝑡) – различные виды рыночных цен и доли бюджета на рынках с государственными ценами для
различных экономических агентов; 𝑚1 + 𝑚2 + 𝑚3 = 𝑚;
𝑢(𝑡) ∈ 𝑈(𝑡) ⊂ 𝑅𝑞 – вектор-функция управляемых (регулируемых) параметров. Значения координат
этого вектора соответствует различным инструментам государственной экономической политики, например,
таким как доли государственного бюджета и бюджетов экономических агентов, различные налоговые
ставки, ставки по гос. облигациям и др.;
𝑎(𝑡) ∈ 𝐴 ⊂ 𝑅 𝑠 - вектор-функция неуправляемых параметров (факторов). Значения координат этого
вектора характеризуют различные зависящие от времени внешние и внутренние социально-экономические
факторы: цены экспортных и импортных товаров, численность населения страны, параметры
производственных функций и др.;
𝑋1 (𝑡), 𝑋2 (𝑡), 𝑋3 (𝑡), 𝑈(𝑡), – компактные множества с непустыми внутренностями; 𝑋𝑖 = ⋃𝑛𝑡=1 𝑋𝑖 (𝑡), 𝑖 =
1, 2, 3; 𝑋 = ⋃3𝑖=1 𝑋𝑖 ; 𝑈 = ⋃𝑛−1
𝑡=0 𝑈(𝑡), 𝐴 - открытое связное множество;
𝑓1 : 𝑋 × 𝑈 × 𝐴 → 𝑅𝑚1 – непрерывное отображение.
2) Подсистема алгебраических уравнений, описывающих поведение и взаимодействие агентов на
различных рынках в течение выбранного года, эти уравнения допускают выражение переменных 𝑥2 (𝑡) через
экзогенные параметры и остальные эндогенные переменные:
𝑥2 (𝑡) = 𝑓2 (𝑥1 (𝑡), 𝑥3 (𝑡), 𝑢(𝑡), 𝑎(𝑡)),
(13)
Здесь 𝑓2 : 𝑋1 × 𝑋3 × 𝑈 × 𝐴 → 𝑅 𝑚2 - непрерывное отображение.
3) Подсистема рекуррентных соотношений для итеративных вычислений равновесных значений
рыночных цен на различных рынках и долей бюджета на рынках с государственными ценами для различных
экономических агентов:
𝑥3 (𝑡)[𝑄 + 1] = 𝑓3 (𝑥2 (𝑡)[𝑄], 𝑥3 (𝑡)[𝑄], 𝐿, 𝑢(𝑡), 𝑎(𝑡))
(14)
Здесь 𝑄 = 0, 1, … – номер итерации; 𝐿 – набор из положительных чисел (настраиваемые константы
итераций, при уменьшении их значений экономическая система быстрее приходит в состояние равновесия,
однако при этом увеличивается опасность ухода цен в отрицательную область; 𝑓3 : 𝑋2 × 𝑋3 × (0, +∞)𝑚3 ×
𝑈 × 𝐴 → 𝑅 𝑚2 – непрерывное отображение (являющееся сжимающим при фиксированных 𝑡; 𝑥1 (𝑡) ∈
𝑋1 (𝑡); 𝑢(𝑡) ∈ 𝑈(𝑡); 𝑎(𝑡) ∈ 𝐴 и некоторых фиксированных 𝐿. В этом случае отображение 𝑓3 имеет
единственную неподвижную точку, к которой сходится итерационный процесс (13), (14).
Вычислимые модели (11), (13), (14) при фиксированных значениях функций 𝑢(𝑡) и 𝑎(𝑡) для каждого
момента времени t определяет значения эндогенных переменных 𝑥(𝑡), соответствующие равновесию цен
спроса и предложения на рынках товаров и услуг агентов в рамках следующего алгоритма.
1) На первом шаге полагается 𝑡 = 0 и задаются начальные значения переменных 𝑥1 (0).
2) На втором шаге для текущего 𝑡 задаются начальные значения переменных 𝑥3 (0)[0] на различных
рынках и для различных агентов; с помощью (13), вычисляются значения 𝑥2 (𝑡)[0] =
𝑓2 (𝑥1 (𝑡), 𝑥3 (𝑡)[0], 𝑢(𝑡), 𝑎(𝑡)), (начальные значения спроса и предложения агентов на рынках товаров и
услуг).
3) На третьем шаге для текущего 𝑡 запускается итерационный процесс (14). При этом для каждого
значения 𝑄 текущие значения спросов и предложений находятся из (13): 𝑥2 (𝑡)[𝑄] =
𝑓2 (𝑥1 (𝑡), 𝑥3 (𝑡)[𝑄], 𝑢(𝑡), 𝑎(𝑡)), через уточнения рыночных цен и долей бюджетов экономических агентов.
Условием остановки итерационного процесса является равенство значений спросов и предложений на
различных рынках с точностью до 0.01%. В результате определяются равновесные значения рыночных цен
на каждом рынке и долей бюджета на рынках с государственными ценами для различных экономических
агентов. Индекс 𝑄 для таких равновесных значений эндогенных переменных мы опускаем.
4) На следующем шаге по полученному равновесному решению для момента времени 𝑡 с помощью
разностных уравнений (11) находятся значения переменных 𝑥1 (𝑡 + 1). Значение 𝑡 увеличивается на
единицу. Переход на шаг 2.
Количество повторений шагов 2, 3, 4 определяются в соответствии с задачами параметрической
идентификации, прогноза и регулирования на заранее выбранных интервалах времени.
4
Рассматриваемая CGE модель может быть представлена в виде непрерывного отображение
𝑓: 𝑋 × 𝑈 × 𝐴 → 𝑅𝑚 , задающего преобразование значений эндогенных переменных системы для нулевого
года в соответствующие значения следующего года согласно приведенному выше алгоритму. Здесь
компакты 𝑋(𝑡) = 𝑋1 (𝑡) × 𝑋3 (𝑡) × 𝑋3 (𝑡), задающие компакт 𝑋 в пространстве эндогенных переменных
определяется множеством возможных значений переменных 𝑥1 и соответствующими равновесными
значениями переменных 𝑥2 и 𝑥3 рассчитываемых с помощью соотношений (13) и (14).
Будем предполагать, что при для выбранной точки 𝑥1 (0) ∈ Int(𝑋1 ) и соответствующей, рассчитанной
с помощью (13), (14) точки 𝑥(0) = (𝑥1 (0), 𝑥2 (0), 𝑥3 (0)) верно включение 𝑥(𝑡) = 𝑓 𝑡 (𝑥(0)) ∈ Int(𝑋(𝑡)) при
некоторых фиксированных 𝑢(𝑡) ∈ Int(𝑈(𝑡)), 𝑎(𝑡) ∈ 𝐴 для 𝑡 = 0, … , 𝑛. (𝑛 – фиксированное натуральное
число). Это отображение 𝑓 определяет дискретную динамическую систему вида (6), (7) в множестве 𝑋, на
траектории которого наложено соответствующее начальное условие:
{𝑓 𝑡 , 𝑡 = 0,1, … }, 𝑥|𝑡=0 = 𝑥0 .
(15)
На базе данного представления ниже рассматривается конкретная CGE модель отраслей экономики.
3.2. Параметрической идентификации CGE модели отраслей экономики
The considered model по статистическим данным экономики Республики Казахстан is presented by the
following nineteen economic agents (sectors):
Рассматриваемая модель по статистическим данным республики Казахстан представлена с помощью
следующих восемнадцати экономических агентов.
Экономический агент № 1. Сельское хозяйство, охота и лесоводство;
Экономический агент № 2. Рыболовство, рыбоводство;
Экономический агент № 3. Горнодобывающая промышленность;
Экономический агент № 4. Обрабатывающая промышленность;
Экономический агент № 5. Производство и распределение электроэнергии, газа и воды;
Экономический агент № 6. Строительство;
Экономический агент № 7. Торговля; ремонт автомобилей и изделий домашнего пользования;
Экономический агент № 8. Гостиницы и рестораны;
Экономический агент № 9. Транспорт и связь;
Экономический агент № 10. Финансовая деятельность;
Экономический агент № 11. Операции с недвижимым имуществом, аренда и услуги предприятиям;
Экономический агент № 12. Государственное управление;
Экономический агент № 13. Образование;
Экономический агент № 14. Здравоохранение и социальные услуги;
Экономический агент № 15. Прочие коммунальные, социальные и персональные услуги;
Экономический агент № 16. Услуги по ведению домашнего хозяйства;
Часть выпущенного продукта экономических агентов - производителей товаров и услуг № № 1–16
используется в производстве, другая часть уходит на инвестиции, а третья продается домашним хозяйствам.
Агенты–производители торгуют между собой промежуточными и инвестиционными товарами.
Экономический агент № 17. Совокупный потребитель, объединяющий в себя домашние хозяйства;
Совокупный потребитель покупает потребительские товары, производимые агентами–
производителями. Кроме того, он покупает импортные товары, предлагаемые внешним миром.
Экономический агент № 18. Правительство, представленное совокупностью центрального,
региональных и местных правительств, а также внебюджетными фондами. Правительство устанавливает
налоговые ставки и определяет сумму субсидий агентам-производителям и размеры социальных
трансфертов домашним хозяйствам. Кроме того, в этот сектор входят некоммерческие организации,
обслуживающие домашние хозяйства (политические партии, профсоюзы, общественные объединения и т.
д.);
Экономический агент № 19. Банковский сектор, включающий в себя Национальный банк и
коммерческие банки.
Здесь экономические сектора № 1-16 являются агентами производителями.
Рассматриваемая модель представляется в рамках общих выражений соотношений (11), (13), (14)
соответственно 𝑚1 = 67, 𝑚2 = 597, 𝑚3 = 34 выражениями, с помощью которых рассчитываются значения
ее 698 эндогенных переменных. Эта модель также содержит 2045 оцениваемых экзогенных параметров.
Задача параметрической идентификации исследуемой макроэкономической математической модели
состоит в нахождении оценок неизвестных значений ее параметров, при которых достигается минимальное
значение целевой функции, характеризующей отклонения значений выходных переменных модели от
соответствующих наблюдаемых значений (известных статистических данных для промежутка времени 𝑡 =
𝑡1 , 𝑡1 + 1, … , 𝑡2 ). Эта задача сводится к нахождению минимального значения функции нескольких
5
переменных (параметров) в некоторой замкнутой области D евклидова пространства с ограничениями вида
(12), накладываемыми на значения эндогенных переменных. В случае большой размерности области
возможных значений искомых параметров, стандартные методы нахождения экстремумов функции часто
бывают неэффективными в связи наличием нескольких локальных минимумов целевой функции. Ниже
предлагается алгоритм, учитывающий особенности задачи параметрической идентификации
макроэкономических моделей и позволяющий обойти указанную проблему «локальных экстремумов»
t2
В качестве области 𝐷 ⊂ ∏t=t
[𝑈(𝑡) × 𝐴(𝑡)] × 𝑋1 (𝑡1 ) для оценки возможных значений экзогенных
1
параметров (значений экзогенных функций 𝑢(𝑡), 𝑎(𝑡) и начальных условий динамических уравнений (11))
(𝑞+𝑠)(𝑡 −𝑡 +1)+𝑚1 𝑖 𝑖
рассматривалась область вида 𝐷 = ∏𝑖=1 2 1
[𝑎 , 𝑏 ], где [𝑎𝑖 , 𝑏 𝑖 ] - промежуток возможных значений
𝑖
параметра 𝑝 ; 𝑖 = 1, … , (𝑞 + 𝑠)(𝑡2 − 𝑡1 + 1) + 𝑚1 . При этом оценки параметров, для которых имелись
наблюдаемые значения, искались в промежутках [𝑎𝑖 , 𝑏 𝑖 ] с центрами в соответствующих наблюдаемых
значениях (в случае одного такого значения) или в некоторых промежутках, покрывающих наблюдаемые
значения (в случае нескольких таких значений). Прочие промежутки [𝑎𝑖 , 𝑏 𝑖 ] для поиска параметров
выбирались с помощью косвенных оценок их возможных значений. Для нахождения минимальных
значений непрерывной функции нескольких переменных 𝐾: 𝐷 → 𝑅 с дополнительными ограничениями на
эндогенные переменные вида (12) в вычислительных экспериментах использовался алгоритм направленного
поиска Нелдера-Мида [6]. Применение этого алгоритма для начальной точки 𝑝1 ∈ 𝐷 можно
интерпретировать в виде сходящейся к локальному минимуму 𝑝0 = argmin 𝐾 функции 𝐾
𝐷,(7)
последовательности {𝑝1 , 𝑝2 , 𝑝3 , … }, где 𝐾(𝑝𝑗+1 ) ≤ 𝐾(𝑝𝑗 ), 𝑝𝑗 ∈ 𝐷, 𝑗 = 1,2, … В описании следующего
алгоритма мы будем считать, что точка 𝑝0 может быть найдена достаточно точно.
Для решения задачи параметрической идентификации рассматриваемой CGE модели на основе
очевидного предположения о несовпадении (в общем случае) точек минимума двух различных функций
предложены два критерия следующего типа:
𝐾𝐴 (𝑝) = √
1
𝑛𝛼 (𝑡2 −𝑡1 +1)
2
𝐴
∑𝑡𝑡=𝑡
∑𝑛𝑖=1
𝛼𝑖 (
1
𝑦 𝑖 (𝑡)−𝑦 𝑖∗ (𝑡)
𝑦 𝑖∗ (𝑡)
2
) , 𝐾𝐵 (𝑝) = √
1
𝑛𝛽(𝑡2 −𝑡1 +1)
2
∑𝑡𝑡=𝑡
∑𝑛𝐵 𝛽 (
1 𝑖=1 𝑖
𝑦 𝑖 (𝑡)−𝑦 𝑖∗ (𝑡)
𝑦 𝑖∗ (𝑡)
2
) . (16)
Здесь {𝑡1 , … , 𝑡2 } – промежуток времени идентификации; 𝑦 𝑖 (𝑡), 𝑦 𝑖∗ (𝑡) – соответственно расчетные и
наблюдаемые значения выходных переменных модели, 𝐾𝐴 (𝑝) – вспомогательный критерий, 𝐾𝐵 (𝑝) –
основной критерий; 𝑛𝐵 > 𝑛𝐴 ; 𝛼𝑖 > 0 и 𝛽𝑖 > 0 – некоторые весовые коэффициенты, значения которых
определяются в процессе решения задачи параметрической идентификации динамической системы;
𝑛𝐵
𝐴
∑𝑛𝑖=1
𝛼𝑖 = 𝑛𝛼 , ∑𝑖=1
𝛽𝑖 = 𝑛𝛽 .
Алгоритм решения задачи параметрической идентификации модели был выбран в виде следующих
этапов.
1. Параллельно, для некоторого вектора начальных значений параметров 𝑝1 ∈ 𝐷, решаются задачи А и
В, в результате находятся точки 𝑝𝐴0 и 𝑝𝐵0 минимума критериев 𝐾𝐴 и 𝐾𝐵 соответственно.
2. Если для некоторого достаточно малого числа 𝜀 верно 𝐾𝐵 (𝑝𝐵0 ) < 𝜀, то задача параметрической
идентификации модели (11), (13), (14) решена.
3. В противном случае, используя в качестве начальной точки 𝑝1 точку 𝑝𝐵0 , решается задача A, и,
используя в качестве начальной точки 𝑝1 точку 𝑝𝐴0 , решается задача B. Переход на этап 2.
Достаточно большое число повторений этапов 1, 2, 3 дает возможность выходить искомым значениям
параметров из окрестностей точек неглобальных минимумов одного критерия с помощью другого критерия
и, тем самым, решить задачу параметрической идентификации.
В результате совместного решения задач A и B согласно указанному алгоритму c использованием с
использованием статистических данных по эволюции экономики Республики Казахстан алгоритма НелдераМидда [6] были получены значения 𝐾𝐴 = 0.015 и 𝐾𝐵 = 0.0063. При этом относительная величина
отклонений расчетных значений переменных используемых в основном критерии от соответствующих
наблюдаемых значений составила менее 0.63%.
В дальнейшем просчет модели за границей периода параметрической идентификации (прогнозный
просчет) с помощью экстраполированных на прогнозный период значений функций 𝑢(𝑡), 𝑎(𝑡) будем
называть базовым просчетом.
Результаты просчета и ретроспективного базового просчета модели на 2008 г., частично
представленные в таблице 1 демонстрируют расчетные (𝑌(𝑡), 𝑌𝑔 (𝑡), 𝑃(𝑡)), наблюдаемые значения и
отклонения расчетных значений основных выходных переменных модели от соответствующих
наблюдаемых значений. Здесь промежуток времени 2000-2007 гг. соответствует периоду параметрической
идентификации модели; 2008г. - период ретропрогноза; Y(𝑡) – валовый выпуск ( × 1012 тенге, в ценах 2000
года); 𝑌𝑔 (𝑡) – ВВП ( × 1012 тенге, в ценах 2000 года); P(𝑡) – индекс потребительских цен в процентах к
предыдущему году; знак «*» соответствует наблюдаемым значениям, знак «Δ» соответствует отклонениям
(в процентах) расчетных значений от соответствующих наблюдаемых значений.
6
Таблица 1. Наблюдаемые, расчетные значения выходных переменных модели и соответствующие
отклонения.
Показатель\ 2000
2001
2002
2003
2004
2005
2006
2007
2008
Год
5.44
6.32
6.47
6.86
7.72
8.52
9.25
9.69
9.84
𝑌 ∗ (𝑡)
𝑌(𝑡)
5.38
6.32
6.47
6.86
7.72
8.52
9.27
9.64
9.82
Δ𝑌(𝑡)
-1.22
-0.02
0.00
0.00
0.05
0.08
0.21
-0.51
-0.26
𝑌𝑔∗ (𝑡)
2.45
2.78
3.05
3.36
3.72
4.09
4.55
5.01
5.18
Yg (𝑡)
2.47
2.78
3.05
3.35
3.72
4.09
4.55
5.01
5.20
Δ𝑌𝑔 (𝑡)
0.88
0.07
-0.04
-0.02
-0.02
-0.02
-0.04
-0.15
0.38
𝑃 (𝑡)
106.4
106.6
106.8
106.7
107.5
108.4
118.8
109.5
𝑃(𝑡)
107.6
106.8
106.9
106.7
107.3
108.2
118.6
109.4
1.13
0.18
0.08
-0.05
-0.23
-0.22
-0.24
-0.05
∗
Δ𝑃(𝑡)
3.3 Оценка некоторых положений макроэкономической теории на базе CGE модели отраслей
экономики
Оценка некоторых положений макроэкономической теории на базе CGE модели
В данном разделе описываются результаты вычислительных экспериментов по оценке отдельных
положений макроэкономической теории о циклических колебаниях макроэкономических показателей при
шоковых изменениях спросов на конечные и инвестиционные товары.
В рамках макроэкономической теории циклические колебания экономических процессов могут
возникнуть при [7]:
- линейной зависимости между текущим доходом и объемом потребительских расходов и
- линейной зависимости между приростом доходов и инвестициями.
С целью проверки этих положений был проведен ряд вычислительных экспериментов по просчету
следующих сценариев a)-f) изменения указанных спросов на конечные и инвестиционные товары с
помощью варьирования соответствующих долей бюджетов экономических агентов:
𝑂1 (𝑡) - доля бюджета домашних хозяйств, идущая на покупку конечных товаров в год t.
𝑂2𝑖 (𝑡) - доля бюджета i-ой отрасли, идущая на покупку инвестиционных товаров в год t. (𝑖 = 1, … ,15).
a) 𝑂1 (𝑡) = 𝑂1 (𝑡 − 1) + 𝑎(𝑌𝑔 (𝑡) − 𝑌𝑔 (𝑡 − 1));
(17)
b) 𝑂2𝑖 (𝑡) = 𝑂2𝑖 (𝑡 − 1) + 𝑏(𝑌𝑔 (𝑡 − 3) − 2𝑌𝑔 (𝑡 − 2) + 𝑌𝑔 (𝑡 − 1)), (𝑖 = 1, … ,15);
(18)
c) Совместное применение сценариев a) и b);
d) Увеличение долей 𝑂1 (𝑡) в 𝑘 раз по сравнению с базовым вариантом;
e) Увеличение долей 𝑂2𝑖 (𝑡) (𝑖 = 1, … ,15) в 𝑙 раз по сравнению с базовым вариантом;
f) Совместное применение сценариев d) и e).
Здесь: 𝑡 = 2010, … , 2015 - время в годах); 𝑎, 𝑏, 𝑙, 𝑘 – положительные константы (𝑙 > 1, 𝑘 > 1).
Применение сценариев a) и b) означает расчет модели при значениях параметров 𝑂1 и 𝑂2𝑖 ,
определяемых формулами (17) и (18) соответственно начиная с 𝑡 = 2010. При этом значения всех
остальных экзогенных переменных модели для сценариев a), b), c), d), e), f) соответствуют ее базовому
просчету. В результате проведения вычислительных экспериментов на базе модели по реализации сценариев
a), b), c), d), e), f) во всех случаях наблюдались циклические колебания переменной 𝑃(𝑡) (индекс
потребительских цен (см. рис. 1, 2). При этом для значений реальных показателей (в частности 𝑌𝑔 (𝑡))
колебательных явлений обнаружено не было.
Результаты экспериментов подтверждают соответствующие положения макроэкономической теории
[7].
3.4.
Подавление
циклических
параметрического регулирования
колебаний
макроэкономических
показателей
методами
В вычислительных экспериментах рассматривалась следующая задача подавления циклических
колебаний уровня потребительских цен, возникших при использовании сценария с) развития экономической
системы, предусматривающего линейную зависимость между объемом потребительских расходов и
текущим доходом и линейную зависимость между инвестициями и приростом доходов.
7
На базе вычислимой модели секторов экономики с использованием сценария c) (см. предыдущий
𝑗
пункт) где 𝑎 = 3 ∙ 10−13 , 𝑏 = 10−13 ) найти значения долей 𝑂𝑖 (𝑡) бюджетов j-го агента-производителя
идущих на покупку товаров и услуг, производимых i-ым агентом-производителем для 2010-2015 г.г.; 𝑖, 𝑗 =
1, … ,16, которые обеспечивали бы минимум следующего критерия 𝐾𝑟 , характеризующего отклонение
расчетных значений индекса потребительских цен 𝑃(𝑡) от соответствующих желаемых значений 𝑃𝑠 (𝑡).
𝐾𝑟 = ∑2015
𝑡=2010 (
𝑃(𝑡)−𝑃𝑠 (𝑡) 2
𝑃𝑠 (𝑡)
) .
Здесь в качестве желаемых значений 𝑃𝑠 (𝑡) использовались расчетные базовые значения индекса
потребительских цен модели без параметрического регулирования. Значение критерия для базового
варианта 𝐾𝑟 = 0.424.
Ограничения на регулируемые параметры:
𝑗
𝑗
𝑗
0.5 ≤ 𝑂𝑖 (𝑡)/ 𝑂̅𝑖 ≤ 1.5; ∑16
𝑖=1 𝑂𝑖 (𝑡) ≤ 1; 𝑖, 𝑗 = 1, … ,16; 𝑡 = 2010, … ,2015.
𝑗
Здесь 𝑂̅𝑖 - базовые значения указанных долей, полученные в результате решения задачи параметрической
идентификации модели по данным 2000 - 2008 г.г.
Ограничения на рост макроэкономического показателя: 𝑌𝑔 (𝑡) ≥ 0.95𝑌𝑔𝑠 (𝑡). Здесь 𝑌𝑔𝑠 (𝑡) расчетные
сценарные значения суммы ВДС отраслей без параметрического регулирования, 𝑌𝑔 [𝑡] - расчетные значения
суммы ВДС отраслей с параметрическим регулированием.
Численное решение данной задачи было получено с помощью алгоритма Нелдера-Мида.
Оптимальное значение критерия 𝐾𝑟 при применении найденного закона параметрического регулирования
оказалось равным 𝐾𝑟 = 0.000844. Базовые значения макроэкономических показателей 𝑌𝑔𝑠 (𝑡) и 𝑃𝑠 (𝑡);
значения, полученные при применении сценария c) и значения, полученные при применении оптимального
закона параметрического регулирования приведены на рисунках 1 и 2.
Анализ результатов вычислительных экспериментов, представленных на рисунках 1, 2 показывает,
что при использовании найденного оптимального закона параметрического регулирования, индекс
потребительских цен на регулируемом периоде практически полностью совпал с желаемыми значениями,
при этом значения суммы ВДС отраслей экономики оказались несколько ниже (за исключением 2015 года)
соответствующих базовых значений.
1.4E+13
1.2E+13
1E+13
8E+12
Базовый
Сценарный
6E+12
Подавление
4E+12
2E+12
0
2011
2012
2013
2014
2015
Рис. 1. Расчетные значения суммарной ВДС отраслей экономики (в тенге, в ценах 2000 года).
8
180
160
140
120
Базовый
100
80
Сценарный
60
Подавление
40
20
0
2011
2012
2013
2014
2015
Рис. 2. Расчетные значения индекса потребительских цен (в процентах к предыдущему году).
Заключение
Приведены результаты по развитию теории параметрического регулирования в рамках двух ее
компонентов: методов оценки оптимальных значений экономических инструментов и исследование влияний
неуправляемых факторов на результаты решения задачи выбора оптимальных значений экономических
инструментов.
Развитые положения теории параметрического регулирования проиллюстрированы на подклассе
вычислимых моделей общего равновесия (CGE моделях).
Показана эффективность предложенного метода параметрической идентификации большеразмерных
CGE моделей.
Приведены результаты вычислительных экспериментов по оценке отдельных положений
макроэкономической теории о циклических колебаниях макроэкономических показателей при шоковых
изменениях спросов на конечные и инвестиционные товары.
Показана эффективность применения теории параметрического регулирования в решении задачи
подавления циклических колебаний уровня потребительских цен.
Полученные результаты могут быть использованы при разработке и осуществлении эффективной
ациклической государственной экономической политики.
[1]
[2]
[3]
[4]
[5]
[6]
[7]
REFERENCES
D. Acemoglu, Introduction to Modern Economic Growth. Princeton, New Jersey, USA: Princeton University
Press. 2008. Р. 1008.
A.A. Ashimov, K. A.Sagadiyev, Yu.V. Borovskiy, N.A. Iskakov, and As.A. Ashimov, “On the market economy
development parametrical regulation theory,” Kybernetes, The Intl. J. of Cybernetics, Systems and Management
Sciences, vol. 37, no. 5, 2008, pp. 623–636.
A.A. Ashimov, B.T. Sultanov, Zh.M. Adilov, Yu.V. Borovskiy, D.A. Novikov, R.M. Nizhegorodtsev, As.A.
Ashimov, Macroeconomic Analysis and Economic Policy Based on Parametric Control. New York: Springer,
2012.
I. Ekeland, R. Temam. Convex analysis and variational problems. Amsterdam: North-Holland publishing
company. 1976.
V. L. Makarov, A. R. Bakhtizin, and S. S. Sulashkin, The Use of Computable Models in Public Administration.
Moscow: Scientific Expert, 2007 (in Russian).
J. A. Nelder and R. Mead, “A simplex method for function minimization,” The Computer Journal, no. 7, 1965,
pp. 308–313.
E.S. Tumanova, N.A. Shagas, Macroeconomics. Мoscow: Infra-M, 2004. (in Russian).
9
Download