УДК 519.216, 519.866
МЕТОДИКА ИДЕНТИФИКАЦИИ НОРМАТИВНОЙ ВЫЧИСЛИМОЙ МОДЕЛИ РЕГИОНАЛЬНОЙ ЭКОНОМИКИ
© В. А. Горбачев
Ключевые слова: ССЕ-модели; вычислимое моделирование; автономные системы дифференциальных уравнений.
В работе проведено исследование модельной системы уравнений в имитационной модели экономического сектора, идентификация вида системы и анализ поведения траекторий решения в малой окрестности точек, соответствующих состоянию равновесия исследуемой экономической системы.
Нормативная балансовая модель общего равновесия региональной экономики (см. [1]) с учетом запасов благ, денег, налогообложения и теневого сектора является частным случаем СОЕ модели. Эта частная модель учитывает потоки благ, по которым имеются статистические временные ряды. Таким образом, модель отражает изменения факторов, действий, предметов потребления, которые можно сравнивать с соответствующей статистикой. Структурно модель представляет собой набор дифференциальных и алгебраических уравнений, часть из которых нелинейны.
При разработке модели предполагалась некая простота математического описания поведения выделенных экономических агентов: в частности, считалось, что материальные и финансовые потоки между агентами определяются некими постоянными долями от запасов материальных благ и запасов денежных средств агентов. С другой стороны, сложность и полнота рассмотрения взаимодействий между агентами привела к тому, что сами уравнения на изменение этих запасов и на начальные данные, как бы просто они ни создавались, довольно громоздки: всего в модели N параметров и М переменных, при этом все они имеют некий понятный экономический смысл. Именно это дает надежду на то, что модель удастся согласовать со статистическими данными при активном использовании современной высокоскоростной вычислительной техники.
В модели рассматриваются следующие экономические агенты: производитель X (перерабатывающий сектор народного хозяйства), производитель У (добывающий сектор); производитель Z (инновационный сектор), соответствующие им теневые части секторов - X', У', Z'; домашние хозяйства (население) Ь ; посредник, или «продавец» Т ; правительство региона (государственные структуры) С; внешние потребители и поставщики (внешний для региона рынок) О; региональная банковская система В (включая долю центрального банка в регионе).
При построении модели сделано важное предположение, что в каждый момент времени Ь весь наличный объем блага г разделен без остатка между агентами, так что в распоряжении каждого агента находится запас блага QV (Ь) . Запас блага QV (Ь) изменяется вследствие его производства (или получения из внешней природной среды) и потребления агентом V, а также вследствие передачи блага от одного агента другому. Последнее предположение
точно записывается в виде уравнения материального баланса:
5У(і) = XV - су - уу - зу - ^(КУ - К%
где XV - объем выпуска (выпуск) или получения из внешней природной среды блага і у агента V ; СУ - конечное потребление блага і агентом V . Потребленное благо не оказывает дальнейшего воздействия на экономические процессы; УУ - текущие затраты, ,1У - материальные затраты. Это тоже форма потребления, но, в отличие от конечного потребления, данные потоки влияют на способность агента производить благо і. К текущим затратам относят блага, которые превращаются в производимый продукт (сырьё, комплектующие), а к капитальным - блага, накопление которых способствует этому превращению (установка оборудования, строительство зданий и сооружений); И,уУ - поток блага і от агента V к агенту ц. Антисимметричная форма члена ^(Ь,уУ - ЬУ^) , описывающего передачи благ, выражает представление о том, что блага не возникают и не исчезают при передаче.
Структурно рассматриваемая модель представляет собой нелинейную автономную систему. В общем виде запишем модельную систему уравнений в следующем виде:
£ = ^ (Т,и),
где — = (Я,\¥,р, в) - набор переменных модели: 5 - вектор запасов материальных благ, IV - вектор денежных запасов, р - вектор цен, в - вектор заработных плат, а и - набор параметров модели. Введем следующие обозначения: 3 = (^1,^2) , где для части переменных статистика известна 3і(іо) = б'0 , а для ^ нет, .^(іо) = Ш2 ; аналогичным образом делится и набор параметров модели Z = (а, и2) : часть параметров Z нам известна из статистики, а часть и2 нет. Без ограничения общности будем считать, что а входит в функцию / . В новых обозначениях мы можем произвести постановку задачи оптимизации параметров и начальных значений модели:
' - = /і(—,иі);
< ^2 = /2(3,и2);
—1 (і0) = в 0, Т2(і0) = и2';
„ ф(—і (іі) - в і) ^ тіп.
Здесь в качестве критерия ф может быть взят индекс близости временных рядов Тейла, или иной критерий близости (среднее квадратичное отклонение, корреляция).
Поставленная задача оптимального управления очень сложна. Даже для небольшого набора уравнений непросто получить решение аналитически, а в нашем случае практически невозможно. Таким образом, единственным способом решения поставленной задачи остается применение численных методов, аналогичных, например, методам из [2]. А именно, предлагается решить данную задачу общим перебором всевозможных вариантов на некоторой достаточно подробной дискретной сети возможных значений искомых параметров с применением параллельных вычислений на суперкомпьютере.
Применение параллельных вычислений на суперкомпьютере обусловлено тем, что такое количество вычислительных операций, которое необходимо для осуществления полного перебора всех значений параметров, не сможет осилить обычный компьютер [3].
Для эффективного применения данного метода необходимо пройти подготовительный этап: при работе с экономикой, не находящейся в переходном периоде, мы можем допустить, что экономика развивается с постоянным темпом роста 7 . С учетом инфляционной составляющей данное предположение можно записать в виде режима квази-сбалансированного роста для исследуемой модели. Для этого введем следующие обозначения — = (5,13) , где
Q - интенсивные переменные (запасы благ и денег у агентов), а P - экстенсивные переменные (цены, заработные платы). В новых обозначениях мы можем записать режим квази-сбалансированного роста:
[ Q = QoeY(t-to);
| P = P0es(t-t0).
Применение режима квази-сбалансированного роста к модели позволяет найти зависимости между параметрами и, выразив зависимые, сократить общее число параметров модели, снизив таким образом нагрузку на вычислительную технику в ходе решения задачи идентификации параметров. В нашем случае применение режима квази-сбалансированного роста позволило исключить 44 параметра и сократить их количество до 12.
Пересчет будет производиться по набору параметров и начальных значений переменных Р2 , которые не могут быть определены из статистики. В качестве критерия ф используем инлекс Тейла. Таким образом, необходимо произвести анализ каждого указанного параметра, чтобы на основании самой сущности параметра, а также основываясь на данных, полученных опытным путем в ходе вычислительного эксперимента для получения работоспособного варианта модели, произвести оценку верхней и нижней допустимых границы ш~, р+ каждого параметра рi, при этом задается дискретная сеть с необходимым шагом разбиения для пересчета параметров Р = (р\, ...,Рп) .
Пересчет будем производить по всевозможным комбинациям значений параметров в рамках установленной сети значений, на каждом шаге проводим проверку критерия оптимальности и сравнение результатов.
Полученный результат можно использовать при изучении последствий изменения политики на рынке, находящемся в режиме сбалансированного роста, предварительно рассчитав начальные условия на исходном и «возмущенном» рынках.
ЛИТЕРАТУРА
1 . Кощеев А.В., Оленев Н.Н. Моделирование взаимодействующих региональных экономических
систем с использованием параллельных вычислений // Труды МФТИ. 2010. Т. 2. № 1 (5).
С. 92-97.
2 . Оленев Н.Н., Фетинина А.И. Моделирование экономики Кировской области с применением
технологий параллельного программирования // Научно-технический вестник СПбГУ ИТМО.
2010. № 1 (65). С. 108-113.
3 . Оленев Н.Н., Фетинина А.И. Параллельные вычисления в идентификации динамической мо-
дели Вятского региона // Вестник Нижегородского университета им. Н.И.Лобачевского. 2009.
№ 6 (1). С. 184-191.
Поступила в редакцию 26 апреля 2011 г.
Gorbachev V.A. The identification of computable model of the regional economy. In this paper we propose the CGE-model of economic sector. This model is represented by a system of differential equations. The analysis of trajectories in the points corresponding to equilibrium of investigated economic system is introduced and studied.
Keywords: CGE-models; computable modeling; independent systems of the differential equations.
Горбачев Владимир Александрович, Российский университет дружбы народов, г. Москва, Российская Федерация, аспирант кафедры нелинейного анализа и оптимизации, e-mail: [email protected]