УДК 519.857.3
М.Ю. Андреев1, Н. П. Пильник*, И. Г. Поспелов1
1 Московский физико-технический институт (государственный университет)
2НИУ ВШЭ
Модель межвременного равновесия экономики Республики Казахстан
В статье представлена модель общего экономического равновесия Республики Казахстан. Модель содержит описание восьми макроагентов. Для четырех агентов (домохозяйство, производитель, банк, собственник) сформулированы оптимизационные задачи, поведение остальных агентов задается сценарием. Модель калибруется на квартальных данных с 2005 по 2010 годы и воспроизводит траектории переменных, характеризующих состояние реального и финансового секторов экономики.
Ключевые слова: общее экономическое равновесие, математическое моделирование, рациональные ожидания, банковская система, экономический кризис.
1. Введение
1.1. Структура модели экономики РК
Описанная ниже модель была разработана в отделе математического моделирования экономических систем Вычислительного центра им. А.А. Дородницына РАН в рамках работ по созданию библиотеки моделей национальной экономики РК, проводившихся в 20092011 гг. в Национальной Академии наук РК под руководством академика А. Ашимова [1]. В то же время модель является продолжением работ по системному анализу развивающейся экономики, начатой в ВЦ РАН А.А. Петровым и И.Г. Поспеловым около 40 лет назад [2],
И-
Рассматриваемая ниже модель — это динамическая однопродуктовая модель межвременного равновесия. Однако от классических моделей такого типа, рассматриваемых, например, в [4], наш подход отличается рядом важных методических особенностей, которые уже описывались в журнале (см. [5]). Основная цель настоящей работы - представить эту методику и реализующую ее технологию моделирования более подробно и предметно. По сравнению с моделью экономики РФ, описанной в [5], и предыдущей версии модели РК [6] здесь были использованы новые теоретические результаты относительно описания взаимодействия собственника и его фирмы. Это позволило в рамках единой модели описать динамику реального сектора, кредитно-денежной сферы и фондового рынка. Насколько нам известно, в существующих динамических моделях общего равновесия такого совместного описания получить до сих пор не удавалось.
Модель описывает развитие во времени полного цикла общественного воспроизводства (производство - распределение - потребление) в экономике в разрезе одного агрегированного продукта - реального ВВП - и одного агрегированного ресурса - однородного труда. Движение продукта и ресурса в модели сопровождается и регулируется оборотом 5 финансовых агрегатов: денежной базы (наличность + банковские резервы), безналичных денег (остатки расчетных счетов), кредитных вложений в экономику, сбережений населения и иностранной валюты.
Потоки продукта, ресурса и финансовых агрегатов удовлетворяют системе балансов, над которой надстраивается модель. Уравнения балансов соответствуют строкам таблиц системы национальных счетов (СНС) [13] и сводного баланса банков [14]. Они прослеживают движение продукта и ресурсов от поступления до использования, а движение финансовых агрегатов от эмиссии до погашения.
Величины потоков определяются решениями субъектов экономики - физических и юридических лиц, а также механизмами согласования этих решений, например, рыночным торгом. Модельные агенты представляют агрегированное поведение множества субъектов, исполняющих сходные функции в экономике. По нашим представлениям, это агрегированное поведение оказывается более регулярным («более рациональным»), чем поведение отдельных лиц и организаций [8], [9]. Модель описывает экономическую динамику как результат взаимодействия следующих 8 экономических (макро) агентов.
1) Производитель представляет все нефинансовые коммерческие предприятия в сфере производства и внутренней торговли.
2) Банк представляет коммерческие финансовые учреждения, включая страховые компании и пенсионные фонды.
3) Население представляет все домашние хозяйства.
4) Собственник представляет индивидуальных владельцев предприятий, топ-менеджеров крупных корпораций, а также компании-холдинги и участников фондового рынка.
5) Центральный банк описывает деятельность ЦБ в области обеспечения расчетов, а также монетарного и валютного регулирования.
6) Государство представляет экономическую деятельность органов административного и государственного управления.
7) Экспортер и
8) Импортер представляют посреднические учреждения в сфере внешней торговли.
Первые 4 агента - агрегаты множества субъектов и их поведение описываются принципами оптимальности. Государство и ЦБ рассматриваются как отдельные влиятельные субъекты и их поведение описывается сценариями. Экспортер и Импортер персонифицируют внешнее воздействие на экономику FK и их описание сводится к сценариям изменения внешних факторов.
Модель основывается на принципе рациональных ожиданий [3, 4-6, 11]. Именно, предполагается, что каждый из оптимизирующих свое поведение агентов (1-4) располагает точным прогнозом всех экзогенных для него переменных, в том числе цен, курсов, процентов на весь период расчета,. В конченом счете траектории этих переменных определятся из условий выполнения балансов (условий равновесия), и тем самым для каждого агента получится самосогласованный прогноз конъюнктуры. Можно сказать, что в моделях межвременного равновесия модельные агенты используют для своих прогнозов ту же самую модель, которую строит исследователь.
Исследование модели дает совершенно неожиданный ответ на вопрос о том, как это предположение согласуется с наблюдаемой неопределенностью и неточностью прогнозов экономической конъюнктуры (см. разд. 5.3, а также [11]).
1.2. Технология разработки модели
Мы использовали разработанную в ВЦ FAH под руководством И.Г. Поспелова систему поддержки моделирования экономики ЭКОМОД [2], специально приспособленную для построения моделей межвременного равновесия. Эта система реализована в среде компьютерной алгебры Maple и поддерживает не только численное, но и аналитические исследование модели. В частности, ЭКОМОД
1) автоматически выводит условия оптимальности для задач агентов;
2) проверяет выполнение материальных и финансовых балансов и инвариантность модели относительно выбора единиц измерения [7];
3) позволяет в большинстве случаев автоматически повторить аналитические преобразования модели при модификации исходных соотношений;
4) проводит расчеты по записи модели в стандартной математической нотации.
Практически все формулы, приводимые ниже, взяты непосредственно из записи модели в системе ЭКОМОД.
Нам представляется, что данная технология позволяет реализовать разумный компромисс между произволом предположений, характеризующих имитационные модели (например, [12]), и комбинированием ограниченного числа стандартных хорошо исследованных блоков, что характерно для большинства макроэкономических моделей (например, [1]).
Для удобства аналитических преобразований модель формулируется в непрерывном времени, а перед расчетом дискретизируется (см. разд. 5.1). Формально предполагается, что экономика функционирует на некотором достаточно большом конечном интервале времени £ £ [Ьо,Т]. На этом периоде заданы траектории всех экзогенных переменных модели и на этот же период планируют свои действия модельные агенты. Значения всех фазовых переменных считаются заданными при Ь = ¿о- Для исключения краевых эффектов на кон-
Т
[10], а также конец разд. 2.1).
2. Описание рационального поведения агрегированных агентов 2.1. Блок «Производитель» (J)
Результаты анализа статистических данных [6] показывают, что производство реального ВВП Уз(¿) в РК можно описать линейной производственной функцией
(I) = АМ](I) + Ве(ь(г-г°(I) + ¥зЪ(1), -М3(I) = ^(I) - кМ3(I), (2.1)
где У3ь(Ь) — экзогенно заданная сезонная составляющая выпуска продукции, RJ(¿) — численность занятых, М,](¿) — эффективные производственные фонды, а JJ(¿) - покупки фондообразующего продукта, который в рамках однопродуктовой модели составляет часть Уз (¿). Мы не сопоставляем величину М,] (Ь) с наблюдаемой величиной основных производственных фондов, поскольку эконометрические оценки дают очень большое значение «темпа выбытия» к1.
Мы предполагаем, что все покупатели и продавцы продукта Уз(¿) ориентируются на одну и ту же рыночную цену ру (Ь), а труд оплачивается по ставке вг (¿), поэтому денежное сальдо покупок и продаж продуктов и ресурсов Производителем составляет
МУз(I) = Ру(I) (Уз(I) - JJ(I)) - 8Г(I) RJ(I). (2.2)
Производитель может брать срочные ссуды у Банка на средний срок иод складывающийся на рынке сложный процент г I (¿). Поэтому изменение остатка ссудной задолженности Ьз (Ь) описывается соотношениями
-
-Ьз(I) = Кз(I) - ^Ьз($, Кз(I) > 0, (2.3)
М
где Кз(¿) — поток вновь полученных кредитов.
1 Статистика России также не дает ясной связи отчетной величины основных производственных фондов с объемами выпуска. С другой стороны, в мировой практике нередки модели, в которых вместо фондов (основного капитала) используется величина потока инвестиций, что согласуется с (2.1) при большом значении к.
MAj(t) = (t) - Zj(t). (2.5)
Ссуды в модели (как и в банковском балансе) рассматриваются как особый инструмент, отличный от денег. Денежное сальдо операций по ссудам выражается в виде2
MLj(t) = Кj(t) - (ft + ri(t)) Lj(t). (2.4)
Производитель реализует на рынке свои акции Aje (t), выпущенные по указанию собственника (агента «С», разд. 2.4), по рыночному курсу, выплачивает по ним дивиденды Zj(t) и в результате получает чистый поток платежей:
dAjc, dt
Считается, что ссуды должны быть обеспечены текущей стоимостью акций:
klLj(t) < ва](t) Aje(t).
Рассмотрение акционерного капитала существенно отличает данную модель от наших прежних моделей межвременного равновесия [2, 3, 5, 6, 8]. Как показал недавно Н.П. Пильник, описанный здесь механизм управления фирмой посредством выпуска акций обеспечивает в условиях совершенной конкуренции эффективность равновесия так же, как и использовавшийся нами ранее менее реалистичный механизм управления собственным капиталом.
В модели все налоги условно сведены к «налоговой нагрузке» nyx(t) на добавленную стоимость, поэтому налоговые платежи производителя описываются как
Tax j (t) = nyx(t) Py (t) Yj (t). (2.6)
Предполагается, что все денежные платежи и поступления производителя проводятся через его расчетный счет в Банке, поэтому остаток этого счета Nj(t) изменяется в соответствии с балансом3:
^ (t) = MYj (t) + MLj (t) + MAj (t) - Tax j (t). (2.7)
В модели платежи описываются непрерывными потоками, но в реальности субъекты экономики осуществляют платежи дискретными порциями, и чтобы выплатить такую порцию, ее надо иметь в запасе на расчетном счете. В результате появляется потребность иметь определенную величину среднего запаса денег для обеспечения заданной средней интенсивности платежа. На этом соображении основана «количественная теория денег», а сама связь между потоками и запасами называется ограничением ликвидности. Эконометрическое исследование показало, что для экономики РК стоит учитывать лишь ограничение ликвидности, связанное с выплатой зарплаты sr (t) Rj (t):
Ts sr(t) Rj(t) < Nj(t). (2.8)
Постоянная времени задает характерное время оборота денег. Мы относим ее к настроечным параметрам модели. Наличие эмпирических ограничений типа (2.8) существенно отличает наши модели межвременного равновесия от классических [4].
Соотношения (2.1) - (2.8)4 описывают технологические и институциональные ограничения на деятельность производителя. В рамках этих ограничений Производитель выбирает траектории своих планируемых переменных — тех, которые отмечены индексом J. Переменные, имеющие другие индексы, - это информационные переменные. Их
2Пара модельных соотношений (2.3), (2.4), грубо говоря, соответствует принятой в бухгалтерии двойной проводке операции кредитования производителя. В исходной записи модели в системе ЭКОМОД это соответствие выдерживается еще точнее.
3При формулировке модели в системе ЭКОМОД используется более сложное описание с двойной проводкой, разделяющее обороты инструментов «расчетный счет» и «денежная база», которым банк осуществляет платежи по поручению своего клиента - Производителя.
4 Здесь и ниже такое обозначение имеет в виду все вынесенные в отдельную строку формулы между (2.1) и (2.4), независимо от того, занумерованы они или нет.
Производитель должен прогнозировать и при выборе траекторий планируемых переменных рассматривать как экзогенные. Некоторые из них, например в (2.1), так и останутся экзогенными переменными для всей модели, но основные: цена ру (Ь) ставка зарплаты ,вг (Ь), процент г I (Ь), курс акций 9аь (Ь) и объем их выпуска Ааь (Ь) - определятся из условий выполнения балансов.
Считается, что Производитель действует в интересах своих акционеров и поэтому стремится максимизировать традиционный функционал ожидаемой полезности реальных дивидендов:
где ху(Ь) — описанный ниже дефлятор дивидендов, А — предпочтение времени, а г] — отвращение к риску, которые рассматриваются как настроечные параметры модели.
Легко понять, что задача (2.9) при заданных начальных значениях фазовых переменных М,](Ь0), Ьз(Ь0), N3(Ь0) ^ 0 при ограничениях (2.1) - (2.8) и без каких-либо ограничений на конце горизонта планирования £ = Т неразрешима из-за формальной возможности создания «пирамиды» задолженности. Поэтому, как предложено в [10], мы дополняем условия (2.1) - (2.8) одним краевым условием роста некоторой линейной формы фазовых переменных:
Значения коэффициентов этой формы при Ь = Т однозначно определяются требованием разрешимости задачи Производителя. Аналогичные терминальные условия ставятся ниже и в задачах остальных агрегированных агентов, и коэффициенты в них тоже опре-
одинаковым для всех агентов настроечным параметром.
Описание поведения агента «Производитель» в модели получается добавлением к соотношениям (2.1) - (2.8), (2.10) условий оптимальности5 траекторий планируемых переменных. Эти условия записываются (автоматически) в терминах дополнительных двойственных переменных, которые потом по возможности исключаются. Двойственные переменные относятся к планируемым переменным агента. Примеры записи подобных условий оптимальности можно увидеть в [2, 6, 8, 11].
В силу условий трансверсальности линейная форма в (2.10) приобретает смысл собственного капитала агента, в котором активы и пассивы оцениваются двойственными переменными к их балансам [2]. Поэтому (2.10) мы называем условием роста капитала
2.2. Блок «Банк» (В)
Блок описывает в агрегированном виде деятельность коммерческих финансовых организаций. Сбережения населения в РК в значительной степени состоят из пенсионных накоплений и страховок [6], поэтому модельный Банк, строго говоря, описывает банки в совокупности с пенсионными фондами и страховыми компаниями. Специфику деятельности последних мы, однако, не учитываем. Материальные затраты на финансовые операции считаются несущественными.
Выдавая кредит Производителю, банк исходит из тех же представлений о строке и проценте, что и Производитель, поэтому банк планирует выдачу ссуд в рамках соотношений
(2.9)
аМз(Т)Мз(Т) + аЬз(Т) (Т) + а^(Т)^(Т) > ^ (аМз(ЩМJЩ+аЬ.](Ш) (Щ + а^(Щ^(Щ) е7(т-ш).
(2.10)
[Ю].
(ср. (2.3), (2.4)):
-
(Ъ) = Кв(I) - А (I), Кв(I) > 0, МЬВ(I) = -КВ+ (& + ф)) $). (2.11)
5Достаточных условий в форме Лагранжа [2].
Величины Lb(t) в (2.11) и Lj(t) в (2.3) — это не одно и то же. Первая — то, что банк планирует дать (предложение ссуд), а вторая — то, что Производитель планирует взять (спрос на ссуды). Их естественное равенство оказывается нетривиальным условием равновесия, которое в конечном счете определит величину процента ri(t) (см. (4.1)). Кроме того, надо помнить, что для Производителя Lj(t) — пассив, поэтому, когда он растет, Производитель получает деньги, а для Банка Lb(t) — актив, поэтому, когда он растет, Банк деньги отдает.
Для финансирования ссуд банк привлекает сбережения Sb (t) на средний срок ß-1 иод складывающийся на рынке сложный процент rs (t). Эта операция описывается аналогично операции привлечения ссуд производителем (2.3), (2.4):
d
-Sb(t) = Vb(t) - ßs Sb(t), Vb(t) > 0, MSb(t) = Vb(t) - (ßs + rs(t)) Sb(t). (2.12)
dt
Здесь Vb(t) — поток новых вкладов, a MSb(t) — денежный итог операций с вкладами.
Кроме сбережений населения, Банк привлекает средства в виде остатков расчетных счетов производителей. Поскольку для этих средств нет регулирующей переменной типа процента, Банк, планируя свои операции, должен полагаться на прогноз остатка на расчетных счетах Nj(t) (см. (2.7)). Банк может только отказаться их вести;
d
-Nb(t) = NNb(t), Nb(t) < Nj(t). dt
Банк по указанию собственника выпускает свои акции Abc (t), реализует их на рынке по курсу öab(t), выплачивает то ним дивиденды Z(t) и в результате получает чистый поток платежей:
МАв(t) = eab(t) - Zb(t). (2.13)
Добавленная стоимость Банка состоит из разницы полученных и уплаченных процентов, поэтому налоговые платежи Банка составляют (ср. (2.6))
Tax в (t) = nyx(t) (п (t) Lb (t) - n (t) Sb (t)).
Поскольку банков много и они систематически переводят деньги друг другу через счета своих резервов в ЦБ, платежным средством Банка следует считать именно эти резервы Rb(t). Изменение совокупных остатков расчетных счетов происходит, когда их обналичивают или когда производят инкассацию. Эти операции также проводятся через резервы по требованиям клиентов Банка. Поэтому6
^^ (t) = MLb (t) + MSb (t) + МАв (t) - Тахв (t) + NNb (t).
Резервы Rb (t) представляют собой платежные средства, даже более ликвидные, чем наличность. В модели мы объединяем их с наличностью (см. разд. 2.3) в один инструмент «денежная база».
Банк обязан выполнять резервные требования: остатки Rb (t) должны превосходить установленную нормативно долю привлеченных средств:
с (Sb(t) + Nj(t)) < Rb(t). (2.14)
Ограничения (2.11) - (2.14) на выбор планируемых переменных Банка (переменные с индексом «В») дополняются линейным терминальным ограничением роста капитала, аналогичным (2.10):
ÜLb(T)LB(Т) + ÜSB(Т) SB(Т) + ÜNB(T)NB(Т) + ÜRB(Т) RB(Т) > > (аТв(t0) Lb(t0) + üSb(t0) Sb(to) + üNb(to)NB(to) + üRb(h) Rb(to)) e7 (T-t°>. .
6При формулировке модели в системе ЭКОМОД используется более сложное описание с двойной проводкой, разделяющее обороты инструментов «расчетный счет» и «резервы».
Мы считаем Банк, как и Производителя, частной коммерческой организацией, действующей в интересах своих хозяев, поэтому целью Банка будем считать максимизацию функционала ожидаемой полезности потока реальных дивидендов, аналогичного (2.9):
Г (WГ ^* -max. (2.16)
Параметры zy(t), Д ^ функционалов (2.9) и (2.16) считаются одинаковыми. Хотя исследование более простых моделей с акционерным капиталом показало, что их можно делать разными, не теряя существования и эффективности равновесия, при разных параметрах функционалов у модели отсутствуют решения типа сбалансированного роста, что мешает использовать отработанную технику идентификации и расчета модели.
Описание поведения агента «Банк» получается добавлением к соотношениям (2.11) -(2.14) условий оптимальности траекторий планируемых переменных Банка (переменных с индексом «В»),
2.3. Блок «Население» (Н)
Население служит источником трудовых ресурсов. Для предложения труда в Казахстане (как и ранее для России) нам не удалось найти удовлетворительного замкнутого описания7. В связи с этим в модели занятость задается экзогенно и в блоке «Население» величина занятости Rr(t) — это информационная, а не планируемая переменная.
Продавая ресурсы труда по ставке зарплаты sr (t), получая из бюджета трансферты SBg (t) и оплачивая потребительский продукт Сн (t) по цене py (t), население получает чистые доходы (возможно отрицательные):
МСн(t) = SBG(t) + Sr(t) Rr(t) - Py(t) Сн(t). (2.17)
Предложение сбережений населением Sh(t) описывается симметрично (2.12): d
-Sh(t) = Vh(t) - ps Sh(t), Vh(t) > 0, MSh(t) = -Vh(t) + (f3s + rs(t)) Sh(t). (2.18)
Величину Sh(t) мы сопоставляем с чистыми процентными сбережениями, т. е. суммой банковских депозитов, страховок и пенсионных накоплений за вычетом потребительских кредитов.
Платежным средством населения считаем наличность Wh (t) (которая в модели объединяется с банковскими резервами), поэтому
d
-Wh (t) = МСн (t) + MSh (t).
Чтобы описать потребность населения в наличности, логично было бы поставить ограничение ликвидности, подобное (2.8), Wh(t) ^ тсpy(t) Сн(t). Оно соответствует и статистическим данным, согласно которым Wh(t)/py(t) Сн(t) ~ const. Однако проведенные эксперименты показали, что такая добавка сильно усложняет расчет модели, не добавляя ей существенно точности. Поэтому мы ставим простейшее ограничение:
Wh(t) ^ 0, (2.19)
которое потом особым образом «смягчим» (см. разд. 5.1).
7Рынок определенно не развит - ставки зарплаты сильно различаются по отраслям и регионам. Кроме того, нет надежных данных о фактически выплаченной зарплате.
В задачу Населения добавляется краевое условие роста капитала (ср. (2.10)):
аШн (Т) Шн (Т) + авн (Т) вн (Т) ^ ^ (аШн(1о) Wн(¿с) + аБн(¿с) вн(¿с)) е^(т-Ч (2.Щ
Как и в большинстве моделей межвременного равновесия, мы предполагаем, что
Население максимизирует ожидаемую полезность потребления:
/Т (^)1-" ^ ШаХ.
Описание поведения агента «Население» получается добавлением к соотношениям (2.17) - (2.20) условий оптимальности траекторий планируемых переменных Населения (переменных с индексом «Н»),
2.4. Блок «Собственник» (С)
Обычно в моделях межвременного равновесия собственность ассоциируется с собственностью на факторы производства - труд и основные производственные фонды. Мы считаем, что основные фонды находятся в собственности фирм-производителей, а сбережения домохозяйств отличаются от капитала собственников. Блок модели «Собственник» в агрегированном виде описывает деятельность физических и юридических лиц, управляющих движением акционерного капитала между секторами и за пределы страны. В известных нам макромоделях других авторов такого экономического агента не было.
Мы описываем Собственника как «собрание акционеров» Производителя и Банка. Это собрание принимает решение о выпуске акций, одновременно обязуясь скупить их по текущему рыночному курсу. Решение о выпуске принимается на основе прогноза роста курса и доходности акций. В соответствии с принципом рациональных ожиданий этот прогноз предполагается точным.
Собственник также получает доход в форме прибыли от внешнеторговых операций (фактически ренту за природные ресурсы страны) и вкладывает свой капитал не только в акции Банка и Производителя, но и в иностранные активы (вывоз капитала).
Собственник определяет номинальный акционерный капитал Производителя А] с(1) (число выпущенных и размещенных акций). Акции имеют доходность (Ь), поэтому с учетом их покупки по курсу 0а](Ь) операции с акциями Банка приносят Собственнику чистый доход8:
MAjc(t) = daj(t) ^(t) + daj(t) Ajc(t), Ajc(t) > 0. (2.21)
Совершенно аналогично собственник определяет номинальный акционерный капитал Банка Abc(t), акции которого имеют доходность dab(t) и продаются по курсу, принося Собственнику чистый доход:
MAbc(t) = -dab(t) (t) + dab(t) Abc(t), Abc(t) > 0. (2.22)
В принципе акции являются более доходными, но и более рискованными вложениями, нежели банковские депозиты. В детерминированной модели риск описать нельзя, но результаты его оценки - можно. Будем считать, что надежные депозиты Собственника Sc(t), страхующие вложения в акции, составляют не менее доли ks от текущей стоимости портфеля акций:
8Модельный собственник представляет совокупность многих реальных собственников, которые владеют только частями пакета акций и поэтому должны ориентироваться на их доходность, а не на общую сумму дивидендов Z] (£). Как часто пишут в подобных случаях авторы моделей рациональных ожиданий: «Собственник в модели один, но он об этом не знает».
Se(t) > ks (t9aj(t) Aje(t) + Qab(t) Abe(t)).
Динамика сбережений описывается так же, как у населения (см. (2.18)):
d
-Se(t) = Ve(t) - ps Se(t), Ve(t) > 0, MSe(t) = -Ve(t) + (f3s + rs(t)) Se(t). (2.23)
Статистическим аналогом величины Se (t) мы считаем иностранные депозиты в банках РК. Предполагается, что собственники сначала выводят свои доходы в более безопасные иностранные активы, а потом часть из них возвращают в страну в виде депозитов, страхующих рискованные вложения в акции.
Кроме указанных операций собственник получает известную наперед прибыль Rentax(t) от экспортно-импортных операций (см. разд. 3.3, 3.4).
Считается, что интересы собственника связаны с его иностранными активами, которые формируются за счет вывоза валюты Ce(t). Эту валюту Собственник должен купить на бирже, заплатив за нее по курсу ww (t).
В конечном счете изменение остатка денег We (t) v Собственника описывается балансом:
dWe
d-^(t) = MAje(t) + MAbe(t) + MSe(t) + Rentax(t) - Ww(t) Ce(t), We(t) > 0.
В выписанных балансах нет потребительских расходов собственника. По нашему мнению, эти расходы слабо влияют на поведение собственника и на возможности вложения капитала. Механизмы движения капитала не замыкаются в границах страны, а извлеченную в разных странах прибыль используют не на потребление. Потребительские расходы собственников, правда, могут повлиять на внутренний рынок. Однако это - потребительские расходы в основном богатых, направлены на импортные товары, поэтому могут быть отнесены тоже к покупке валюты. Но остаются еще расходы на услуги внутри страны. Их мы действительно не учитываем.
В задачу Собственника добавляется краевое условие роста капитала:
aWe (Т)We (Т) + aSe (Т)Se (Т) + aAje (Т) Aje (Т) + aAbe (Т) Abe (Т) ^ ^ (aWe (ta)We (Т) + aSe (to)Se (Т) + aAje (ta) Aje (h) + aAbe (ta) Abe (ta)) e7 (T-to).
Предполагается, что Собственник максимизирует ожидаемую полезность потока вывоза валюты Ce (t):
[ (^ )1-,? e-A4t ^ max.
Описание поведения агента «Собственник» получается добавлением к соотношениям (2.21) - (2.24) условий оптимальности траекторий планируемых переменных Населения
C
3. Описания механизмов государственного управления и внешнеэкономических связей
3.1. Блок «Государство» (G)
Механизмы государственного управления экономикой описываются в модели двумя способами: прямо и косвенно. Косвенное (с точки зрения структуры модели) воздействие состоит в установлении различных нормативов: ставок налогов и пошлин, нормы резервирования и пр. Прямое - в реализации планов государственных закупок, выплат населению из бюджета, накоплению валютных резервов, финансированию дефицита бюджета. Оно описывается в блоках «Государство» и «Центральный банк». Экономическая деятельность
государства, как это и принято в большинстве моделей, описывается сценарием, а не решением оптимизационной задачи. Мы связываем эту общепринятую практику с тем, что государство - реальный субъект, наделенный волей и рефлексией, в отличие от рассмотренных выше агрегированных макроагентов, поведение которых можно описывать принципами оптимальности.
В блоке «Государство» описывается исполнение расширенного бюджета, объединяющего национальный, региональные и местные бюджеты, а также внебюджетные фонды.
Доход широкого бюджета формируется за счет налогов и пошлин. Государство назначает ставки, а поступления оно прогнозирует, поэтому величина налоговых поступлений является информационной переменной Тахх(Ь).
Расходы широкого бюджета состоят из выплат населению БВс(Ь), расходов ру(Ь) Сс(Ь) на оплату государственного потребления9 Со(¿) по рыночной цене ру(Ь) и расходов -Мы (¿) ЬйШс (^ на внешние выплаты в валюте Ьс!Шс(1). Величина Ьс!Шс(1) может быть отрицательной, что отвечает поступлениям в бюджет потока денег за счет внешних заимствований.
Положительная разница между доходами и расходами бюджета так или иначе покрывается внутренними заимствованиями, накопление которых образует внутренний государственный долг Вс(Ь). Мы предполагаем, что кредитором государства выступает только Центральный банк, поэтому кредиты можно считать беспроцентными10.
Перечисленные выше статьи доходов, расходов и заимствований изменяют остаток счетов государства Ас(£):
^(1)= ТаХх(Ь) - Ру(I) Ос(1) - ^(I) ЬйШо(1) + ^(I), Аа(1) > 0, Вс(1) > 0.
Эти счета ведет ЦБ, поэтому их можно отнести к модельной денежной базе.
3.2. Блок «Центральный банк» ( СВ)
ЦБ представляет собой часть системы государственного управления, но по закону он обязан обеспечивать собственную безубыточность и обладает большой самостоятельностью маневрирования валютными резервами и кредитами, почему мы и выделяем в модели его описание в отдельный блок. В модели поведение ЦБ описывается сценарием изменения обменного курса и спросом государства на кредиты.
ЦБ покупает или продает валюту на бирже в количестве11 Ус в что изменяет величину золотовалютных резервов Ксв (1), а также дает государству беспроцентный кредит по его запросу. Расходы ЦБ на покупку валюты и кредитование государства покрываются эмиссией Есв (1), которая изменят денежную базу МВсв (1):
СВ (I) = Есв(I), Есв(I) = ^(I) Усв(г) + ^(I). (3.1)
3.3. Блок «Экспортер» (ЕХ)
Международная торговля в модели описывается двумя сравнительно простыми блоками: «Экспортер» и «Импортер».
Экспорт Республики Казахстан в основном представляет собой сырье [6]. Для сырьевого экспорта всегда оказываются существенными спросовые и транспортные ограничения.
9Согласно системе национальных счетов государство покупает по себестоимости и потребляет продукт отрасли «Управление» в составе ВВП, поэтому государственные служащие, военные, полицейские и т. д. входят в число занятых Кг (£), а их заработная платя покрывается расходами ру (?) Со (¿0 •
10Формально проценты есть и выплачиваются, но возвращаются в бюджет как прибыль ЦБ, поскольку в РК, как и во всех других странах, кроме США и Великобритании, ЦБ является государственным учреждением.
11 Эта величина становится отрицательной при валютной интервенции.
По этим причинам мы сочли возможным считать в модели и объем, и цену экспорта (в иностранной валюте) величинами, заданными сценарием12.
Покупка продукта для экспорта Expex (t) на внутреннем рынке требует затрат py (t) Expex (^.Продукт Exp ex (t) затем продается на внешнем рынке по сложившейся мировой цене pe(t), что приносит Экспортеру валютну ю выручку pe(t ) Exp ex (t)- Она продается на бирже по курсу ww (t), принося доход в национальной валюте. Экспортер платит пошлину Tax ex (t) в национальной валюте по переменной ставке пе x(t), а остаток денег отдает Собственнику как ренту Zex (t)'
Zex (t) = -p y (t) Exp ex (t) + ww (t) pe(t) Exp ex (t) - Tax ex (t), Tax ex (t) =пе x(t) Ww (t)p e(t) Exp ex (t).
3.4. Блок «Импортер» (IM)
В рамках односекторной модели потребность в импорте рациональным поведением агентов объяснить невозможно, поэтому объем импорта ImpiM (t) мы задавали найденным эмпирически соотношением
ImPlM(t) = ImpSe(t) + гс Сн(t) + ij Jj(t) + ieExpEx(t) + Ww(t) ebw (t"io>, (3.2)
где Impse(t) - сезонная составляющая спроса на импорт, а гс, ij, iw, bw - постоянные, определенные эконометрическим оцениванием.
Покупка импортного продукта по внешней цене Pi(t) требует затрат в валюте pi(t)lmpim(t) Эта валюта покупается на бирже по курсу ww(t). Продажа продукта ImpiM (t) по внутренней цене приносит импортеру выручку py (t)ImpiM (t) в национальной валюте. Из этой суммы оплачивается импортная пошлина в объеме TaxiM(t) по переменной ставке. Остаток денег передается Собственнику как еще одна составляющая ренты Zim (t)-
Zim (t) = Py (t) ImpiM (t) - ww (t)Pi(t) ImpiM (t) - TaxiM (t),
TaxiM (t) = ni x(t) pi(t) Imp im (t).
4. Описание взаимодействия агентов
Связь блоков агентов описывается условиями равновесия отдельных рынков. На рынке продукта балансируется предложение продукта Производителем и Импортером и спрос со стороны тех же производителей, государства и населения:
0 = Yj (t) - Jj (t) - Gg (t) - Expex (T) + Impm (t) - Сн (t).
На рынке труда балансируется неэластичное предложение труда со стороны Населения со спросом на труд со стороны Производителя:
0 = - Rj (t) + Rr (t).
Равновесие на рынке банковских ссуд выражается в равенстве спроса и предложения
LB (t) = Lj (t). (4.1)
Аналогично описывается и равновесие на рынке сбережений, с учетом того, что сбережения в модели делают два агента: Население (2.18) и Собственник (2.23):
Sb (t) = Sh (t) + Sc (t). (4.2)
12 Для РФ нам удалось описать экспорт и импорт как эндогенные переменные, определяемые рациональным поведением агентов, но только в трехпродуктовой модели [5].
На валютном рынке балансируется спрос на иностранную валюту со стороны Собственника, Импортера, Государства и ЦБ с предложением иностранной валюты экспортерами:
ww (t)pe(t) Exp ex (t) = pi(t) ImpiM (t) + Cc (t) + Vc в (t) + ww (t) LdWG(t).
Условия равновесия на рынке акций состоят в том, что Банк и Производитель передают собственнику сальдо выплаты дивидендов и выручки от продажи акций:
0 = -МАв(t) - MAbc(t), 0 = -MAj(t) - MAjc(t).
Легко видеть, что отсюда в силу (2.5), (4.3), (2.13), (2.21), (2.22) следуют естественные соотношения Zj(t) = daj(t)Aj(t) и Zb(t) = dab(t)Ab(t), определяющие доходности акций Производителя daj (t) и Банк a dab (t).
Еще одно взаимодействие описывает передачу в бюджет собранных налогов и пошлин, а также получение собственником ренты от внешней торговли:
0 = -Тахх (t)+TaxJ (t)+TaxB (t)+TaxEX (t)+TaxIM (t), 0 = -Rentax(t)+ZEX (t)+ZIM (t).
5. Результаты исследования модели
5.1. Окончательный вид соотношений модели
Полная система соотношений модели состоит из подсистем, выписанных в разд. 2-4, дополненных условиями оптимальности поведения агентов. Следует отметить, что это - уже автоматически упрощенная, особенно в части описания взаимодействий, система. Исходная запись в системе ЭКОМОД соответствует бухгалтерским двойным проводкам финансовых потоков и, как и они, содержит много вспомогательных переменных и зависимых соотношений. Но зато такая запись позволяет автоматически проверять корректность системы балансов, которую, как показывает опыт, вручную без ошибок для более чем трех инструментов написать практически невозможно. Соотношения равенства активов и пассивов, например (4.1), (4.2), фактически получаются как первые интегралы исходных уравнений в потоках, которые ЭКОМОД находит автоматически в процессе проверки балансов. Еще одним примером такого интеграла служит вытекающее из (3.1) в силу остальных соотношений равенство
МВсв (t) = Rb (t) + WH (t) + Wc (t) + AG(t),
которое показывает корректность определения эмиссии в модели: денежная база МВсв(t), выпущенная ЦБ, равна остаткам ликвидных средств остальных агентов.
Ограничения-неравенства, такие как (2.19), дают в системе условий оптимальности условия дополнительности вида
W(t) ^ 0, ф) ^ 0, ф) • W(t) = 0, (5.1)
где <¿>(t) - двойственная переменная к ограничению (2.19). Эти условия мы регуляри-зируем: либо считаем, что на всем интервале расчета один из сомножителей равен О, либо заменяем «угол» (5.1) на плоскости W,<р графиком гладкой убывающей функции W = f (ф), содержащей свои настроечные параметры. Такого рода регуляризации следует рассматривать как самостоятельный содержательный этап построения модели [4].
С экономической точки зрения соотношения (5.1) задают вырожденную функцию спроса на деньги. Аргумент этой функции составляет главный предмет дискуссий в экономической теории. Эти споры соотношение (5.1) однозначно решает в пользу величины ^(t), которая в конечном счете окажется некоторой комбинацией показателей экономической конъюнктуры (цен, процентов, курсов и т. п.). Решив вопрос об аргументе, мы, как и в традиционных макроэкономических моделях, подбираем форму функции спроса (регуляризацию (5.1)) эмпирически.
В результате аналитических упрощений, переобозначений13 и регуляризаций модель становится системой из 16 независимых уравнений на 16 эндогенных переменных. Кроме этого остаются соотношения, которые по значениям эндогенных переменных позволяют вычислить траектории изменения других показателей, например, дефицита бюджета или объема международных резервов. Заметим, что среди эндогенных переменных остаются некоторые ненаблюдаемые двойственные переменные из условий оптимальности.
Параметры модели в основном определяются независимыми эконометрическими оценками (например, коэффициенты в (2.1) и в (3.2)). Оставшиеся 15 настроечных параметров определяются в процессе калибровки.
Система уравнений модели обладает группой симметрий масштабной инвариантности, которой при экспоненциальных траекториях экзогенных переменных отвечают автомодельные решения - на них каждая эндогенная переменная растет с постоянным темпом, определяемым ее размерностью в базовых единицах измерения [2, 7]. Сопоставление автомодельных решений со статистикой позволяет определить большинство настроечных параметров и зафиксировать типичные «темпы роста размерностей»: продукта д, труда Л, внутренних цен г и курса валюты ш, а также темп роста капитала 7 = г + д в краевых условиях (2.10), (2.15), (2.20), (2.24). Отметим, что сравнительная статика автомодельных решений содержит известную 1Б-ЬМ модель [6].
Система уравнений нормируется на выбранное автомодельное решение, и получается расчетная система, мотивировать которую без развернутого описания в разд. 2-4 было бы очень затруднительно.
Система дискретизируется с шагом 1 квартал по неявной схеме и решается довольно сложным алгоритмом, подробно описанным в [2, 6]. Траектория рассчитывается не шаг за шагом, а вся сразу, поскольку речь фактически идет о поиске ограниченного решения динамической системы, неустойчивой в обоих направлениях времени.
Расчет ведется от некоторого момента в прошлом до некоторого момента в будущем. На периоде прошлого экзогенные величины задаются в соответствии со статистикой. Результат расчета на этом участке используется для верификации модели. На периоде будущего экзогенные величины задаются независимым прогнозом, результат расчета на этом участке дает прогноз развития экономики.
5.2. Результаты калибровки модели общего равновесия Республики Казахстан
Калибровка модели осуществлялась по восьми представленным ниже переменным на квартальных данных с 1 квартала 2005 года по 4 квартал 2010 года. Черной линией изображены статистические величины, серой — модельные.
На представленных ниже графиках для численной оценки точности модельных оценок приведено значение обратного коэффициента Тейла:
где уг, уг — соответственно статистический и модельный ряды данных. Хотя Кот может принимать отрицательные значения, такие случаи можно игнорировать как результаты калибровки модели с крайне низкой точностью. Этот показатель в некотором смысле схож с популярным коэффициентом детерминации В2, но В2 малопригоден для калибровки моделей САРЭ, в которых наряду с постоянно растущими (а следовательно, нестационарными) показателями, такими как ВВП, присутствуют близкие к стационарным показатели темпов изменения (дефлятор ВВП).
13Система ЭКОМОД следит за корректностью переобозначений и помнит их историю. Кроме того, она помнит, из какого исходного соотношения получено то или иное конечное.
7 15
ras ¡¡от гю7 гооа дни joio joos ¡emt 2007 гам моя 2010
Рис. 1. Рояльный ВВП (КоТ 0.993). реальное потребление домохозяйств (КоТ 0.941). млрд тенге
ioos iow i»7 2 cos гмэ 2010 im 2 сад гоо? гад» im 2010
Рис. 2. Реальное валовое накопление (КоТ 0.817). реальный импорт (КоТ 0.921). млрд тенге
Рис. 3. Темп роста дефлятора ВВП (темп инфляции) (КоТ 0.957). чистые депозиты домашних хозяйств (КоТ 0.85). млрд тенге
Рис. 4. Кредиты фирмам (КоТ 0.903). расчетные счета фирм (КоТ 0.957). млрд тенге
Наихудший показатель точности модели получается для величины чистых сбережений населения S(t). Видимо, это связано с тем, что они описаны в модели как банковские депозиты, а фактически в PK депозиты населения меньше выданных населению потребительских кредитов (см. разд. 2.2), да и сама динамика чистых сбережений весьма нестабильна. В этом пункте модель явно требует совершенствования.
Уже после калибровки параметров модели по указанным 8 переменным была построена модельная оценка индекса фондовохх) рынка KASE, который в рамках модели описывается как курс акций. Полученные результаты представлены на следующем храфикс (черная ли-
ния статистика, серая модельные расчеты). Специально отметим, что статистические данные изображены с дневным шагом.
В течение рассматриваемого периода на казахстанском фондовом рынке надувался один из самых больших пузырей в мире. Модель четко показывает два периода роста индекса: начало 2006 года и начало 2007 года. Максимальный модельный и статистический уровни совпадают достаточно точно. В 2008 году происходит обвал, который тоже описывается моделью, хотя модель дает спад в два раза меньше, чем в реальности. Тем не менее после отскока в конце 2009 года фондовый индекс KASE достигает ровно тех уровней, которые предсказывает модель. Способность модели совместно описывать динамику реального сектора, кредитно-денежной сферы и фондового рынка, насколько нам известно, в мировой практике уникальна.
5.3. Сильный магистральный эффект
Результаты прогнозных расчетов мы не приводим. В модели достаточно много экзогенных переменных, и разработка прогноза их совместного изменения представляет собой отдельную еще не решенную задачу. Главным для нас было показать работоспособность модели, и с этой точки зрения важнейшим результатом работы стало новое подтверждение наличия сильного магистрального эффекта [11].
D-J 0Í5 0.1 0.15 0.1 0.05
о
-o.os -0.1 -0.1Б
Рис. 5. Фондовый индекс KASE
S
0.3 5
-0.15 .1
Рис. 6. Модельные значения чистых депозитов домохозяйства (млн тенге), модельные значения темпов инфляции
Хотя модель представляет собой краевую задачу и в случае обпито положения изменение экзогенных переменных в будущем влияют на всю предысторию, при идентифицированных значениях параметров магистральное свойство оказывается выраженным так сильно, что влияние будущих) на настоящее практически затухает за один шаг, и модель фактически ведет себя как причинная динамическая система. Это хорошо видно на рис. 6, где показана реакция расчетных траекторий на изменение притока иностранных депозитов. Серая линия то же, что на рис. 3, а пунктирная результат расчета в предположении увеличения иностранных депозитов в первом квартале 2010 года.
При других значениях параметров модели отклонения пунктирной лини от расчетной начинаются ранее момента (предвидение агентов).
С математической точки зрения это выражается в том, что все положительные действительные части ляпуновских показателей в результате идентификации становятся достаточно большими, так что в силу асимптотики Лапласа вклад в ограниченное решение будущих возмущений фактически сводится к вкладу их текущих значений.
С содержательной точки зрения это означает, что набор сложившихся институтов и реально используемых инструментов так «зажимает» выбор агентов, что их оптимальное решение фактически оказывается не зависящим от будущих возмущений. Представляется, что так и должна быть устроена работоспособная экономическая система. В такой системе принцип рациональных ожиданий не противоречит наблюдаемой неопределенности прогнозов. Деятельность, сильно зависящая от точности прогнозов, так или иначе элиминируется или локализуется как слишком рискованная.
Исследование выполнено при финансовой поддержке РФФИ в рамках научных проектов Ж№ 11-01-00644, 11-01-12136-офи-м-2011, 12-01-00916-а (2011-2012), 12-01-31333, 1201-31189; при финансовой поддержке РГНФ в рамках научно-исследовательского проекта РГНФ ("Построение многопродуктовой макроэкономической модели на основе дезагрегирования финансовых балансов"), проект №11-02-00241а; ПФИ ОМН РАН №3, проект 3.14; ПФИ Президиум РАН №14, проект 109. Исследование осуществлено в рамках программы фундаментальных исследований НИУ ВШЭ в 2012 году. Расчеты выполнены на суперкомпьютере МВС-ЮОК МСЦ РАН.
Литература
1. Pekarski S.E., Merzlyakov S. Макроэкономическая политика экспортоориентированной экономики (на англ. яз.). W14 2 «Научные доклады лаборатории макроэкономического анализа». - М.: НИУ ВШЭ, 2011. - 32 с.
2. Turnovsky S.J. Methods of Macroeconomic Dynamics. — MIT Press Cambridge, Mass. — 2000. - 687 p.
3. Андреев М.Ю., Врежещ В.П., Хохлов M.А., Пет,ров A.A., Поспелов И.Г. Жукова A.A., Здановская В. С. Опыт моделирования экономической динамики Республики Казахстан в период мирового финансового кризиса. — М.: ВЦ РАН, 2010. — 161 с.
4. Андреев М.Ю., Пильник Н.П., Поспелов П.Г. Моделирование деятельности современной российской банковской системы // Экономический журнал ВШЭ. — 2009. — Т. 13, № 2. - С. 143-171.
5. Андреев М.Ю., Пильник Н.П., Поспелов П.Г. Сильный магистральный эффект в модели рациональных ожиданий современной банковской системы России // Журнал новой экономической ассоциации. — 2009. — Т. 1, № 2. — С. 70-84.
6. Андреев М.Ю., Поспелов И.Г., Поспелова, И.П., Хохлов М.А. Технология моделирования экономики и модель современной экономики России. — М.: МИФИ, 2007.
7. Ашимов A.A., Султанов Б.Т., Адилов Ж.М. [и др.]. О применении теории параметрического регулирования для вычислимых моделей общего равновесия // Экономика и мат. методы. — 2010. — Вып. 3.
8. Кондратов П.А., Поспелова Л.Я., Шананин A.A. Обобщенный непараметрический метод. Применение к анализу товарных рынков // Труды МФТИ. — 2010. — Т. 2, № 3. - С. 32-45.
9. Петров A.A., Поспелов П.Г. Математические модели экономики России // Вестник РАИ. - Т. 79, № 6. - 2009. - С. 492-506.
10. Пильник Н.П., Поспелов П.Г. О естественных терминальных условиях в моделях межвременного равновесия // Экономический журнал ВШЭ. — 2007. Т. 11, № 1. — С. 1-33.
11. Поспелов И.Г. Равновесные модели экономики в период мирового финансового кризиса // Труды МФТИ. - 2009. - Т. 1, № 4. - С. 66-83.
12. Поспелов П.Г., Хохлов М.А. Метод проверки размерности для исследования моделей экономической динамики // Математическое моделирование. — 2006. — Т. 18, № 10. — С.113-122.
13. Сайт Агентства Республики Казахстан по статистике. http://www.stat.kz/vs_scheta/int_scheta/Pages/default.aspx
14. Сайт Национального Банка Республики Казахстан. http://www.nationalbank.kz/?switch=rus
15. Форрестер Дж. Мировая динамика. — М.: Наука, 1979. — 167 с.
Поступим в редакцию 31.12.2012