CC BY
42№3(15)2009
К. В. Кетова, И. Г. Русяк
Идентификация и прогнозирование обобщающих показателей развития региональной экономической системы
В работе предложена экономико-математическая модель региональной экономической системы, которая наряду с вещественным капиталом учитывает человеческий капитал. Построен алгоритм идентификации неизвестных параметров модели на основе генетического алгоритма с вещественным кодированием и метода Хука— Дживса. Расчеты по идентификации и прогнозированию проведены на примере Удмуртской Республики.
1. Введение. Математическая модель экономической системы региона
При формулировке математической модели экономической системы региона учтем, что региональная экономика взаимодействует с внешней экономической средой посредством кредитов, инвестиций, налогообложения, дотаций, трансфертов и субвенций (рис. 1). При этом далее механизм воздействия кредитов и инвестиций на экономику региона будем отождествлять между собой.
56
Рис. 1. Схема взаимодействия экономики региона с внешней экономической средой
^-
- №3(15)2009
г> *
В качестве показателей макроэкономическом системы региона примем: 2
1) объем произведенной продукции У;
2) стоимость основных производственных фондов (ОПФ) или производственного капи- з: тала К; §"
о
3) стоимость человеческого капитала Н; |
4) объемы инвестиций в производственный капитал I; ^
5) объемы инвестиций в человеческий капитал ]; ^
6) объем потребления С;
7) доходы регионального бюджета С.
Таким образом, соответствующий паспорт неизвестных задачи имеет вид:
(У, К, Н, I, ], С, С) г,
где I — время.
Коэффициенты V и п будут определены ниже. Производственную функцию примем в виде
У = Г(К, Н) = АК а Н3,
где параметры А, а, 3 определяются по статистическим данным экономической системы [Русяк, Кетова (2008)].
Запишем основное балансовое уравнение макроэкономической модели региона, используя схему воспроизводства экономики, приведенную на рис.2:
У + Б( + Т —Ы! — Яг = 1г+1 + Л+1 + С(+1, (1)
где Б( — внешние инвестиции в экономику региона; Я( — внешний долг;
— налоговые отчисления в федеральный бюджет; Т — дотации, трансферты, субвенции. Нижний индекс в данном случае характеризует зависимость переменных от времени.
Рис. 2. Схема цикла воспроизводства региональной экономики
В относительных переменных балансовое уравнение имеет вид:
Во + Я + Вь — = 1,
№3(15)2009
где з0 = ^ — уровень потребления;
зк = -I — норма инвестиций в основные производственные фонды (капиталовложения); Е
зь = — — норма инвестиций в человеческий капитал;
В г
зь = — — уровень внешних инвестиций;^ = У + Т — Я.
>3 §
» 3
г
0
1
я» '3
о X
Л §
X о
3
«I
а §
а а
'3
а
I!
о с
19 г
а ?
о о
«I
3 X а
<8
О &
3 55
0 &
3 §
г
3
■е-
3
1
«I
»о 5
Балансовое уравнение (1) содержит эндогенные и экзогенные параметры. К последним относятся параметры (В, Я, , Т) . При этом внешний долг Я{ определяется динамикой поступления кредитов В{, политикой возврата кредита, или темпом его погашения ст, и процентной ставкой по кредитам р.
Для описания формирования дохода регионального бюджета используем схему, приведенную на рис. 3.
Пусть
Рис. 3. Схема бюджетного взаимодействия региона с внешней средой
N = + МЯ,
где ,МЯ — налоги, собираемые на территории региона, поступающие в федеральный и региональный бюджеты соответственно.
Обозначим
р = Я = N1
Р = N ' Р = N ' где рри рЯ — соответствующие доли отчислений,
Р I „Я
тогда рр +рЯ = 1.
58
^-
- №3(15)2009
Объем налогов определяется через долю п объема реализованной продукции У (рис. 1): N = пУ. Федеральное бюджетное регулирование выразим через пропорцию V (рис.1) возврата средств в виде дотаций, трансфертов, субвенций как долю от уровня региональных налогов [Булгаков В. К., Булгаков О. В. (2006)]:
Т = vNR,
к *
с
X
<8
0
1
оа
где =рЯ пУ.
Консолидированный доход региона можно записать в виде:
С = + Т = (^)рЯ пУ.
Для описания динамики ОПФ, или К, человеческого капитала Н и внешнего долга 2 используем эволюционные уравнения, рассмотренные в работах [Кетова (2007)];[Русяк, Кето-ва (2007)].
Таким образом, общая постановка задачи макроэкономической динамики включает следующие соотношения:
Е = У + Т — — Я = I + ] + С —В, (2)
50 + + 5Ь — = 1, (3)
У = Г (К, Н) = АК а Н3, (4)
С = 50Е, I = 5кЕ, ] = 5ЛЕ, В = 5ьЕ, Т = vрRпУ, = рГпУ, (5)
(РЯ )
Е = иГ(К, Н), ш = 1 + пр Fv^- —1, С = (^)рЯ пУ, (6)
К(Г) = ВкЕ (Г) — цК (Г), К 0 = К (Г 0), Кт = К (Гт ), (7)
Н(Г) = ЁЕ(Г) — ХН(Г), Н0 = Н(Г0), Нт = Н(и), (8)
¿(1) = ВьЕ(Т) — о2(Г), 20 = 2(Г0), 2т = 2(Гт); Я(Г) = (р + о)2(Г), (9)
где ц и х —коэффициенты выбытия производственного и человеческого капитала соответственно;
о — как было указано выше, темп погашения кредита;
ё — усредненная по всем возрастам доля населения, участвующего в общественном производстве.
2. Алгоритм идентификации модели
Систему уравнений (2)—(9) можно представить в виде:
X = ^ (х, г, а, Г), (10)
Г(х, г, а) = 0, (11)
х, г, а)< 0, (12)
№3(15)2009 ^
где х = (К, Н, 7) — вектор фазовых переменных; Г1(Г2 —вектор-функции ограничений;
г = (У, Р, Ь, 1°, 1, 7, С, С, В, Т, Я, , ) — вектор дополнительных переменных (здесь Р — общая численность населения (рори!а^ог);[. — численность трудовых ресурсов (!аЬог),Ч° — численность населения трудоспособного возраста); а = (а/) = (Л, а, р, г|, в, х, ^, у, Рр , ря, ст, р, 5о, 5к, 5ь, 5Ь) вектор параметров системы, причем
0 < а™ < а! < атах,
где |П, отах —минимальные и максимальные значения соответствующих параметров.
Задача идентификации состоит в следующем. Известно (например, из опыта) поведение системы х(г). Необходимо подобрать коэффициенты а/ таким образом, чтобы отклонение поведения системы, которое определяется решением уравнений (10) при наличии ограничений (11) и (12), от заданного поведения было бы минимальным. | Точное решение задачи идентификации во многих случаях не представляется возможным. Алгоритмы ее приближенного решения формально мало отличаются от алгоритмов решения задач математического программирования [Федоренко (1994)]. При практической реализации алгоритмов существуют проблемы, связанные с дифференцируемостью функционалов и наличием многих экстремумов целевой функции.
Реализация вычислительной схемы требует конечно-разностной аппроксимации задачи (10). На отрезке [tо, Гт] вводится сетка tо <11 <...<Г^ — Гт; при этом дифференциальные I уравнения (10) заменяются дискретным аналогом
з <8
8 а
'3
У"+1 =л( yn+1, yn, zn+1, zn, a, tn), n = 0, N -1, (13)
§ где Л — оператор, определяющий метод численного интегрирования; у — сеточная функция, аппроксимирующая х.
Условие минимизации отклонения поведения системы (13) от заданного поведения хс I в заданных точках оси времени t е [Т0,Т1 ] принимается в виде
I (yt )2
| 7(а)= £ ^-^ т'г; а/ е[а Г, а Г ]. (14)
<» Г €[Т0,Т1] (хС) а/
3 X а
<8
0 &
1
8 рования: &
3 §
г
3
Если J(a) ^ 0, то yt ^ x'd, a ^ a0, где a0 решение задачи идентификации (10)—(14). Таким образом, задача идентификации сводится к задаче математического программи-
J(a)^ ext; G^a) = 0, G2 (a) < 0. (15)
Выбор метода оптимизации определяется особенностями решаемой задачи. Хорошо ■§■ разработанный аппарат математического программирования предлагает множество различных алгоритмов, но при их практической реализации существуют такие проблемы, как многоэкстремальность и выбор начального приближения. Для решения оптимизационных задач, имеющих сложности такого рода, предлагается применять гибридные алгоритмы
60
_ №3(15)2009
(ГБА): генетические алгоритмы (ГА) в комбинации с градиентными методами оптимизации § [Курейчик (2000)];[Тененев, Паклин (2003)]. ¡^
При использовании ГА для решения оптимизационных задач параметры оптимизации з: представляются в виде кодированных значений (генов) [Рутковская (2004)]. Совокупность ге- I нов образует хромосому. Из хромосом составляется популяция. Каждой хромосоме ставится § в соответствие функция приспособленности, которая выступает мерой качества решения, ^ описываемого данной особью. В качестве функции приспособленности выступает целевая ^ функция задачи (14).
На каждой итерации ГА путем применения генетических операторов происходит эволюция исходной популяции, т. е. происходит изменение информации, содержащейся в хромосомах. Блок-схема ГА представлена на рис.4.
Рис. 4. Блок-схема генетического алгоритма
Для представления параметра оптимизации в виде гена хромосомы применяются различные способы кодирования [Рутковская (2004)], [Тененев (2006)]. В случае вещественного
кодирования параметр приводится к безразмерному виду с помощью преобразования:
а. — ат'п _
= ^т_а1_ , = 1 N. (16)
а т ах —а т|П
Соответственно обратное преобразование будет иметь вид:
а1 = ат|п + Ь, (атах — ат|п), I = Ш (17)
Пример хромосомы в случае вещественного кодирования представлен на рис.5. После того как каждой особи поставлено в соответствие значение функции приспособленности (14), осуществляется отбор особей с целью выбора тех представителей, которые будут участвовать в создании новой популяции.
№3(15)2009
0,23496 0,78154 0,23955 0,00092
>3 §
«I
3
г
0
1
я» '3
о X
л §
X
0
3
«I
а §
£ з <8
8 а
'3
г
1 I!
о с
19 г
а ?
о о
«I
3 X а
<8
О &
3 55
0 &
3 §
г
3 ■&■
3
1
«I
»о 5
Ь, Ь2 Ь, Ь„
Рис. 5. Пример хромосомы в случае вещественного кодирования
Существует множество типов операторов отбора: например, метод рулетки, турнирный отбор, ранговый отбор. Будем применять турнирный отбор, суть которого заключается в том, что из популяции выбираются т5 особей и среди них определяется особь, имеющая наилучшее значение функции приспособленности. Эта операция продолжается до тех пор, пока не будет получено нужное число родительских особей для формирования следующего поколения. Значение т5 называется размером турнира.
После создания родительской популяции к ней применяются так называемые генетические операторы. Среди генетических операторов, как правило, используются операторы скрещивания, мутации и инверсии.
Результатом реализации оператора скрещивания является формирование потомков на основе родительских особей. Суть действия оператора скрещивания состоит в выборе точки разрыва внутри хромосомы с последующим обменом информации, располагающейся за точкой разрыва, между двумя родительскими особями. Пример механизма действия одноточечного оператора скрещивания для вещественного кодирования приведен на рис.6.
Отметим, что операторов скрещивания, или механизмов формирования особей потом-
Родительские особи
Потомки
Ць]. •лу/н- ~Ь1
То' разр ки ыва
..Ь} к,-
.ы
¿Я2- ■Ь1
Рис. 6. Механизм действия одноточечного оператора скрещивания
ков, существует достаточно много. Воспользуемся оператором скрещивания Б1п2, который приведен в работе [Тененев (2006)]:
с, = аЬ1 + (Ь,2, / = 1, N.
(18)
где с, — значение /-го элемента потомка;
Ь, — значение /-го гена родительской особи;
З = 2—1(2—и) ; а = 1 —(; и € [0; 1] — случайная величина, имеющая равномерный закон распределения;
Я — параметр, отражающий количество разрядов для кодирования вещественного числа.
Следующим генетическим оператором является оператор мутации. Он предназначен для внесения случайных изменений в информацию, содержащуюся в хромосоме, и, как оператор скрещивания, имеет достаточно разнообразные варианты реализации. Критерием ка-
^-
- т3(15)2009
чества оператора мутации служит его способность производить небольшие изменения с высокой долей вероятности и с низкой вероятностью — большие.
Будем использовать оператор мутации, определяемый по формуле [МиЫепЬе1п, БсЬНегкатр-УооБеп (1993)]:
ст = с, + в, г, р,, (19)
где р, = 2~к"; к е [4; 20];
г, е [4; 5] — диапазон шага рекомбинации; в, е {—1,+1} — направление рекомбинации.
В генетических алгоритмах также применяется оператор инверсии. Суть его работы заключается в разбиении хромосомы случайным образом на две части с последующей их перестановкой (рис. 7).
До инверсии После инверсии
с1с2...с, сш...сн
СО ¡с
Точка разбиения
Рис. 7. Действие оператора инверсии
Реализация операторов скрещивания, мутации и инверсии носит вероятностный характер. При практической реализации генетических операторов задают некоторое фиксированное число рг, затем производится генерация случайного числа р из отрезка [0; 1]. Если р < рг, то оператор выполняется. Для оператора скрещивания рг е [0,5; 1], в то время как для операторов мутации и инверсии рг е[0; 0,1] [Рутковская (2004)].
Еще одним важным аспектом ГА является способ формирования следующего поколения особей, который может происходить на основе различных стратегий. В нашем случае была использована элитарная стратегия. При выборе данной стратегии одна или несколько лучших особей переходят в следующую популяцию в случае, если в результате применения генетических операторов были получены особи, среди которых лучшее значение функции приспособленности хуже, чем у лучшей особи (одной или нескольких) предыдущего поколения. Это сделано для того, чтобы предотвратить потерю лучшего решения при переходе от старого поколения к новому.
Для решения задачи идентификации параметров модели (15) был использован генетический алгоритм с вещественным кодированием с применением следующих генетических операторов и стратегий:
1) стратегия отбора — турнирный отбор (размер турнира — 4);
2) оператор скрещивания, который основан на вероятностном применении двух типов операторов — применяется двухточечный оператор, и в данном случае между двумя точками разрыва хромосомы происходит получение новых генов с использованием оператора В1п2 вида (18) (рис. 8);1
3) оператор мутации. Каждый ген подвергается мутации с вероятностью, равной —,
тд
где тд — количество генов. Новое значение гена вычисляется по формуле (19);
№3(15)2009
Родительские особи
,1,1 ,1 ,.1 ,.1 А ,1 С, с2...с( с(+1 • • • СУ . с)+Г"'«
1 очки разрыв а
г2 с2 г2 Ч 2 ' / с2 с2 г2 г2 у+1" N
Потомки
Г1 Г1 Г1 12'/ Вш2 с1 г1
>
Г2 Г2 Г2 С1 С2•■•'/ Вш2 г2 г2
Рис. 8. Действие двухточечного оператора скрещивания совместно с оператором В1п2
з 4) оператор инверсии, который состоит в перестановке частей хромосом, разделенных
| случайным образом;
3 5) стратегия формирования следующего поколения — элитарная стратегия (одна элитная особь).
* В используемой гибридной схеме генетический алгоритм с вещественным кодировани-
| ем является основным методом оптимизации. Дополнительным методом является метод Ху-
| ка—Дживса.
■ИЗ Метод Хука—Дживса имеет ряд преимуществ [Лесин, Лисовец (1995)], а именно:
5 1) в отличие от методов оптимизации, использующих производные, он не требует диф-§ ференцируемости целевой функции;
3 2) в процессе организации вычислительного поиска не требуется решения задач одномерной минимизации (как, например, при методе сопряженных направлений);
I 3) в отличие от генетического алгоритма обладает более высокой скоростью сходимо-
| сти вблизи точки экстремума.
£ Алгоритм метода Хука—Дживса состоит из процедуры, определяющей направление
| убывания минимизируемой функции (исследующий покоординатный поиск), и процедуры
| перемещения в направлении убывания.
| Применяемый гибридный алгоритм основан на последовательной работе генетического ° алгоритма и метода Хука—Дживса. Суть работы гибридного алгоритма заключается в том, | что на каждой (очередной) итерации основного метода предпринимается попытка улучше-| ния решения с помощью дополнительного метода оптимизации. Благодаря этому гибридно ный алгоритм повышает адаптационные свойства каждого из методов, применяемых в данной схеме, и дает в большинстве случаев лучшие результаты, чем отдельные методы. § Тестирование алгоритмов производилось на наборе стандартных функций, таких как ^ функции Розенброка, Растригина, Экли [Евтушенко (1982)] и др. Во всех случаях достигнута § приемлемая точность определения глобального экстремума, при этом подтверждена высо-
8 кая работоспособность данного гибридного алгоритма.
6
з
| 3. Результаты программной реализации модели
з
Рассмотрим представленную модель на примере статистических данных по Удмуртской
| Республике (УР).
§ Использовать статистические данные переходного периода экономики в их первоначальном виде практически невозможно, поскольку они испытывают значительные и трудно-
64
_ №3(15)2009
объяснимые колебания, вызванные перестройкой экономической системы. Поэтому для по- § вышения точности моделирования вначале производилось сглаживание колебаний временных рядов методом скользящей средней. з:
Полученные таким образом исходные временные ряды основных экономических пара- £ метров УР за 1996-2006 гг., построенные по статистическим данным [Удмуртия в цифрах § (1996-2006)];[Закон о бюджете УР (1996-2006)];[0сновные фонды УР (1996-2006)];[Доходы, ^ расходы и потребление домашних хозяйств (2006)];[0тчетность Федеральной налоговой ^ службы (2007)];[0сновные итоги деятельности Государственного совета УР (2007)], представлены в табл. 1.
При решении задачи идентификации, следуя «Отчетности Федеральной налоговой службы» (2007), полагали, что доли отчислений в федеральный и региональный бюджеты до 2002 г. включительно одинаковы: рг = рЯ = 0,5. Начиная же с 2006 г. эти доли составляют: рг = 0,7 и рЯ = 0,3, при этом в данном временном промежутке они изменяются линейным образом. Напомним, что рг +рЯ = 1.
В качестве минимизируемого функционала выбран функционал (14), где вектор моделируемых макроэкономических характеристик, определяемый из системы уравнений (2)—(8), задан в виде у = (У, К, Н, I,}, NF, NЯ, Т); t е[ 1996,2006]. Известное поведение системы хд = (У, К, Н, I, ], NF, NЯ, Т) , в заданных точках оси времени t €[1996,2006] определяется из табл. 1. Вектор варьируемых параметров имел вид а = (а,) = (г|, ё, х, ^, у, 50, , 5^).
Решение задачи идентификации представлено в табл.2 и 3.
В табл.4 приведены средние относительные отклонения расчетных экономических показателей от их статистических значений за период ретропрогноза, определяемые по формуле:
1 I у г' _xt I
Су =- £ у = {У, к, Н, I, л NF, ^, Т, О}, (20)
МХд t €[Т0,Т1] Хд
где МХд — количество заданных точек показателя х на оси времени в периоде ретропрог-ноза.
Из анализа величины Су следует, что для годового производства продукции среднее за 1996-2006 гг. отклонение составляет 4,09%, величина основного регионального капитала — 8,05%, человеческого капитала — 8,11%. Расхождение между расчетными и статистическими данными по уровню суммарного потребления и региональным бюджетом за исследуемый период составляет 11,20 и 7,07% соответственно.
Полученные таким образом значения параметров системы (см. табл. 3) были приняты базовыми в прогнозных расчетах на 2007-2016 гг. Результаты прогнозирования абсолютных и относительных (на демографическую единицу или единицу трудовых ресурсов) параметров экономики УР на основе предложенной модели приведены в табл. 5 и 6.
На рис. 9-11 представлены изменения расчетных макроэкономических параметров с 1996 по 2016 г. Здесь же отражены результаты сравнения расчетных и статистических данных за период ретропрогноза. На рис.12 и 13 приведены кривые изменения соответствующих удельных показателей за этот же период.
Идентификация и прогнозирование обобщающих показателей развития региональной экономической системы
Таблица 1
О»
о»
Исходные данные макроэкономических показателей экономики УР в ценах 2006 г.
Обозначение 1996 г. 1997 г. 1998 г. 1999 г. 2000 г. 2001 г. 2002 г. 2003 г. 2004 г. 2005 г. 2006 г.
У, млн руб. 113971,6 103 819,9 119805,9 130 629,8 140412,9 1 28 794,0 125 684,2 125 957,8 128 024,2 131 817,0 136 655,1
К, млн руб. 775 903,5 721 594,0 675 827,1 599 839,1 500 876,7 409 730,0 381 166,2 371 453,7 373 078,2 386 568,4 399 996,8
Н, млн руб. 67 539,8 70 974,9 77 373,4 86 067,0 93 306,4 94 036,6 88 870,5 83 584,7 82 175,0 86617,0 92 848,9
В, млн руб. — — 5161,8 5144,3 5016,2 2913,5 2944,3 2898,6 1165,9 1 266,7 1707,4
/, млн руб. 24 222,7 20423,8 22 193,3 22516,9 25 522,2 21 949,4 1 9 584,6 17 465,3 21 255,0 26 117,9 32 656,5
млн руб. 7265,2 7977,9 8700,5 9579,9 10 973,1 11 663,2 9781,8 7867,2 8011,8 9638,5 11 668,3
С, млн руб. — — — — — — — — - - -
, млн руб. — — — — — — 22 233,8 23 573,5 29459,1 36 924,2 46 288,3
Ын, млн руб. — — — — — — 21 275,1 20 388,6 1 9 197,8 18 891,5 18 541,3
Т, млн руб. — — — — — — 4485,8 3007,8 2648,9 2268,4 2405,3
О, млн руб. — — — — — — 25 760,9 23 396,5 21 846,7 21 159,9 20 946,7
1°, чел. 933 400 948 100 959 300 973 000 986 800 968437 972 298 979 074 988 658 993 052 989 600
/_, чел. 639 600 669 700 700 800 744 700 767 800 792 900 782 000 781 200 772 900 792 900 766 000
Р, чел. 1 612 600 1 607 700 1 604 000 1 601 400 1 595 600 1 588100 1 578 200 1 568 200 1 560 200 1 552 800 1 544400
1
§
§
1 I
Таблица 2
Расчетные данные макроэкономических показателей экономики УР в ценах 2006 г., млн руб.
Обозначение 1996 г. 1997 г. 1998 г. 1999 г. 2000 г. 2001 г. 2002 г. 2003 г. 2004 г. 2005 г. 2006 г.
У 113 216,9 115 700,7 118071,3 1 20 343,0 122 530,7 124 650,5 126 718,4 128 750,6 130 514,4 132 016,7 1 33 263,8
К 765 214,6 689 227,6 623 095,7 565 575,9 515581,8 472 164,8 434497 401 856,1 373 055,8 347 560,5 324 905,2
н 64 821,9 68106,1 71 374,7 74 622,7 77 845,9 81 041,2 84 206,6 87 340,9 90 263,8 92 965,9 95 439,7
в 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0
1 19 994,7 20433,3 20 852,0 21 253,2 21 639,6 22 013,9 22 379,1 22 172,9 21 903,7 21 576,4 21 195,2
] 8348,1 8531,3 8706,1 8873,6 9034,9 9191,2 9343,7 9257,6 9145,2 9008,5 8849,4
с 66 991,5 68461,2 69 863,9 71 208,1 72 502,6 73 756,8 74 980,4 74 289,4 73 387,6 72 290,8 71 013,8
ые 20 786,6 21 242,6 21 677,9 22 094,9 22 496,5 22 885,8 23 268,1 26 005,4 28 758,1 31 513,3 34 257,9
N° 20 786,6 21 242,6 21 677,9 22 094,9 22 496,5 22 885,8 23 268,1 21 277,1 19 172,1 16 968,7 14681,9
Т 2910,1 2974,0 3034,9 3093,3 3149,5 3204,0 3253,1 2974,7 2680,4 2372,4 2052,7
О 23 696,7 24216,6 24712,8 25 188,2 25 646,0 26 089,8 26 521,2 24 251,9 21 852,6 19341,1 16 734,7
I
Таблица 3
Параметры модели
А а (3 ц £ X и V «Л ш
1,120 0,180 0,820 0,131 0,750 0,046 0,367 0,140 0,704 0,208 0,088 0,760
О» 41
Таблица 4
Средняя погрешность определения макроэкономических показателей экономики УР за период ретропрогноза
Параметр У К Н / } м" Т О
Погрешность, (,у, % 4,09 8,05 8,11 12,86 12,28 11,60 8,97 9,80 7,07
К. В. Кетова, И. Г. Русяк
rn3(15)2009
У-МУ\млнруб. 160
--Расчетные данные
• —Статистические данные
1996 1998 2000 2002 2004 2006 2008 2010 2012 2014 Рис. 9. Динамика изменения валового регионального продукта
2016
ЛИ 0", млн руб. 900
— Расчетные данные
— Статистические данные
1996 1998 2000 2002 2004 2006 2008 2010 2012 Рис. 10. Динамика изменения производственного капитала
2014
2016
НЛОмлн руб. 900
— Расчетные данные
— Статистические данные
t,¿od
1996 1998 2000 2002 2004 2006 2008 2010 2012 Рис. 11. Динамика изменения человеческого капитала
2014
2016
68
№3(15)2009
у', с', тыс. руб. чел.
100 90 80 70 60 50 40
1- 1 1 1 -1- -Г 1 1
1 1 1 1 1 1
1 1 1 1
1 1 1 1 1 1 1 1
1 1 1 1 "Г - - ■ Г - "I - - - г 1 I 1
1 I 1
1 1 1 1 1 I 1 1
Г, год
1996
1998
2000
2002
2004
2006
2008
2010
2012
2014
2016
Рис. 12. Динамика изменения удельных значений валового регионального продукта (1) и общественного потребления (2)
к1, тыс. руб. /чел. 1400
1200 1000 800 600 400
- 1 1 1 -г 1 1 I 1
1 X 1 1 1 I 1
1 1 I
1 1 1 1 1 I 1 I 1
© 1 1 1 I 1 1 1
1 1 1 1 " I 1 1 '
А1, тыс. руб./чел. 180
160 140 120 100 80
{угод
1996 1998 2000 2002 2004 2006 2008 2010 2012 2014 2016 Рис. 13. Динамика изменения фондовооруженности (1) и капиталовооруженности (2) труда
* с
X
<8
0
1 оа ¡с
Анализируя полученные данные, можно сделать следующие выводы:
• экономика Удмуртской Республики находится на этапе стабильного роста всех макроэкономических показателей;
• темп роста экономики невелик;по удельному ВРП составляет около 2%, по уровню удельного общего потребления — около 1% в год.
В табл. 6 приведены удельные показатели (см. формулу (20));верхний индекс р указывает на то, что показатель рассчитан на единицу общей численности населения региона, верхний индекс I — что показатель приведен в расчете на единицу трудовых ресурсов.
41
о
Идентификация и прогнозирование обобщающих показателей развития региональной экономической системы
Таблица 5 3
Прогнозные значения абсолютных макроэкономических показателей экономики УР в ценах 2006 г., млн руб.
Обозначение 2007 г. 2008 г. 2009 г. 2010 г. 2011 г. 2012 г. 2013 г. 2014 г. 2015 г. 2016 г.
У 133 515,8 133 768,9 134 022,4 1 34 276,3 134 530,7 134 785,6 135 041,2 135 297,4 135 554,4 135 319,8
К 321 287,7 317 766,1 314 337,7 311 000,3 307 751,8 304 590,1 301 512,9 298 518,4 295 604,5 287 203,4
н 95 895,2 96 349,8 96 802,7 97 254,0 97 703,7 98151,8 98 598,3 99 043,3 99486,7 99 907,2
в 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0
1 21 267,4 21 307,7 21 348,1 21 388,5 21 429,1 21 469,7 21 510,4 21 551,2 21 592,1 21 554,8
] 8912,0 8928,9 8945,8 8962,8 8979,8 8996,8 9013,8 9030,9 9048,1 9032,4
с 71 094,0 71 228,8 71 363,8 71 499,0 71 634,4 71 770,2 71 906,3 72 042,7 72 1 79,5 72 267,8
ые 34 300,2 34 365,2 34430,3 34495,5 34 560,9 34 626,4 34 692,0 34 757,9 34 823,9 34 763,6
N° 14 700,0 14 727,9 14 755,8 14 783,8 14811,8 14 839,9 14 868,0 14 896,2 14 924,5 14 898,7
Т 2058,0 2061,9 2065,8 2069,7 2073,6 2077,5 2081,5 2085,4 2089,4 2085,8
0 16 758,1 16 789,8 16 821,6 16 853,5 16 885,4 16 917,4 16 949,5 16 981,7 17 013,9 16 984,5
141
Таблица 6
Прогнозные значения относительных макроэкономических показателей экономики УР в ценах 2006 г., тыс. руб/чел.
Обозначение 2007 г. 2008 г. 2009 г. 2010 г. 2011 г. 2012 г. 2013 г. 2014 г. 2015 г. 2016 г.
УР 86,97 87,49 88,03 88,60 89,21 89,87 90,57 91,33 92,14 92,67
к' 405,20 393,46 388,90 385,78 384,20 383,78 384,23 385,70 387,76 383,10
Ь1 120,94 119,30 119,77 120,64 121,97 1 23,67 1 25,65 127,97 130,50 133,26
¡1 26,82 26,38 26,41 26,53 26,75 27,05 27,41 27,84 28,32 28,75
/ 11,24 11,06 11,07 11,12 11,21 11,34 11,49 11,67 11,87 12,05
ср 46,31 46,59 46,87 47,18 47,50 47,85 48,23 48,63 49,06 49,49
10,92 10,98 11,05 11,12 11,20 11,28 11,37 11,46 11,57 11,63
/
1
§
§
1 I
№3(15)2009
Список литературы
Булгакове. К., Булгаков О. В. Моделирование динамики обобщающих показателей развития регио- Н нальных экономических систем России II Экономика и математические методы. 2006. Т.42. № 1.
Доходы, расходы и потребление домашних хозяйств в 2002-2006 годах (по итогам выборочного ^
00
обследования бюджетов домашних хозяйств): Стат. сб. Государственного комитета РФ по статистике, ^
Евтушенко Ю. Г. Методы решения экстремальных задач и их применение в системах оптимизации. М.: Наука, 1982.
Закон о бюджете Удмуртской Республики на 1996,..., 2006 гг.
КетоваК.В. Об одной задаче макроэкономической динамики региона с учетом факторов экономического развития II Периодический научно-теоретический журнал «Вестник ИжГТУ». Ижевск: Изд-во ИжГТУ. 2007. №3(35). С. 33.
КурейчикВ. М. Генетические алгоритмы II Перспективные информационные технологии и интеллектуальные системы. 2000. № 1. С. 18.
Лесин В. В., Лисовец Ю. П. Основы методов оптимизации. М.: Изд-во МАИ, 1995.
Основные итоги деятельности Государственного совета УР третьего созыва по правовому обеспечению социально-экономического развития УР. Управление аналитического обеспечения и информационных ресурсов при аппарате Государственного совета УР. Ижевск. 2007.
Основные фонды Удмуртской Республики: Стат. сб. Госкомстата УР. 1996-2006. Ижевск, 2007
Отчетность Федеральной налоговой службы «Отчет о начислении и поступлениям налогов, сборов и иных обязательных платежей в бюджетную систему Российской Федерации», по состоянию на 1 января 2007 г.
РусякИ.Г., КетоваК.В. Построение производственной функции экономической системы региона с учетом человеческого капитала II Вестник МГУ. 2008. №3. (Экономика).
РусякИ. Г., КетоваК.В. Экономико-математическая модель анализа и прогноза фактора человеческого капитала II Научно-практический журнал «Экономика, статистика, информатика. Вестник УМО», [раздел «Статистика и математические методы в экономике»]. М.: Изд-во ГОУ ВПО МЭСИ. №2. 2007. C. 56-60.
РутковскаяД. Нейронные сети, генетические алгоритмы и нечеткие системы. М.: Горячая линия Телеком, 2004.
ТененевВ.А. Применение генетических алгоритмов с вещественным кроссовером для минимизации функций большой размерности II Интеллектуальные системы в производстве. Ижевск: Изд-во ИжГТУ. 2006. № 1. C. 93-107.
ТененевВ.А., Паклин Н. Б. Гибридный генетический алгоритм с дополнительным обучением лидера II Интеллектуальные системы в производстве. 2003. № 2. С. 181.
Удмуртия в цифрах: стат. сб. Госкомстата УР. 1996-2006. Ижевск, 2007.
Федоренко Р. П. Введение в вычислительную физику. М.: Изд-во МФТИ, 1994.
Muhlenbein H., Schlierkamp-Voosen D. Predictive Models for the Breeder Genetic Algorithm: I. Continuous Parameter Optimization II Evolutionary Computation. 1993. Vol. 1. № 1. P. 25-49.
С.32.
<s о
2006.
