Математическое моделирование динамических процессов системы управления валовым внутренним продуктом Республики Беларусь

ГУМАНИТАРНЫЕ И СОЦИАЛЬНО-ЭКОНОМИЧЕСКИЕ НАУКИ

УДК 338

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ДИНАМИЧЕСКИХ ПРОЦЕССОВ СИСТЕМЫ УПРАВЛЕНИЯ ВАЛОВЫМ ВНУТРЕННИМ ПРОДУКТОМ РЕСПУБЛИКИ БЕЛАРУСЬ*

Докт. техн. наук, проф. КУЛАКОВ Г. Т., канд. экон. наук, доц. НАГОРНОВ В. Н., канд. техн. наук, доц. КУЛАКОВ А. Т., канд. техн. наук ЦИЛИБИНА В. М.

Белорусский национальный технический университет, Институт экономики Национальной академии наук Беларуси

Важнейшими показателями устойчивого инновационного социально-экономического развития Республики Беларусь являются валовой внутренний продукт (ВВП) и темпы его роста. За шесть лет номинальный ВВП в стране увеличился в 3 раза: с 12,33 млрд дол. США в 2000 г. до 36,95 млрд дол. США в 2006 г.

Успешное выполнение Государственной программы инновационного развития Республики Беларусь на 2006-2010 гг. при достигнутом темпе роста ВВП требует дальнейшего совершенствования государственной системы управления экономикой на основе анализа прошлых лет с разработкой различных сценариев кратко- и долгосрочных прогнозов.

В основу построения структурной схемы динамической модели эффективной системы управления ВВП Республики Беларусь положены теории ноосферного устойчивого социально-экономического развития [1] с использованием кибернетических основ управления экономикой [2, 3] и общей теорией оптимальных систем управления [4-6].

Для построения математической модели системы управления номинальным внутренним валовым продуктом Республики Беларусь с 1990 по 2006 г. с прогнозом до 2010 г. необходимо было определить динамику объекта, представленную безынерционной частью объекта управления и его инерционной частью с лагом (рис. 1).

* Работа выполнена при финансовой поддержке Белорусского республиканского фонда фундаментальных исследований.

Рис. 1. Структурная схема динамической модели системы управления ВВП Республики Беларусь [4]: >(Г) - основная регулируемая величина (номинальный ВВП Республики Беларусь в млрд дол. США по курсу Национального банка); ум(Г) - выходная переменная модели инерционного участка объекта регулирования; хзд - заданное (плановое) значение регулируемого параметра; Л - внутреннее возмущение; /2 - крайнее внешнее возмущение; W0Пp) - передаточная функция опережающего участка объекта регулирования без лага (запаздывания); №нн (р) = (р)е ТуР -то же инерционного участка объекта с лагом; №Лм.(р) -полная модель инерционного участка объекта, равная

передаточной функции инерционного участка ^ (р)е х,р ; №нм(р) - неполная модель инерционного участка объекта без запаздывания; Т - лаг (запаздывание) по каналу регулирующего воздействия; Шъ(р) - передаточная функция крайнего внешнего возмущения; ^укв1(р), ^укв2(р) - соответственно передаточные функции упредителей внутренних и внешних угроз; хР(Г) - управляющее воздействие

Для определения численных значений коэффициентов передаточных функций воспользуемся анализом событий (для выделения этапов и участков) и методами аппроксимации переходных характеристик [5, 6]. В качестве основной переходной характеристики объекта управления принимаем данные Министерства статистики и анализа Республики Беларусь номинального ВВП, выраженного в дол. США по курсу Национального банка Республики Беларусь, на 1990-2006 гг.

Динамику первого этапа моделирования

61,9 й у

—Г й

+У(* -ту) ),

где у(0 - оригинал выходной переменной (номинальное значение ВВП в дол. США по курсу Нацбанка Республики Беларусь); х(() - скачок входной переменной, эквивалентный снижению ВВП из-за развала СССР на 43 млрд дол. США (к1 = -43); ту - лаг, соответствующий одному году, принятый при моделировании в ускоренном масштабе времени равным 10 с.

Повышенная колебательность ВВП первого этапа обусловлена малым численным значением коэффициента демпфирования Е = 0, 44, так как (1) можно записать в общем виде

Т2у"а - ту) + 2ТЕу'а - ту) + у

где Т2 = 61,9 с2, или Т = 7,87 с.

Дифференциальное уравнение (1) представим в виде передаточной функции Ж1(р) разомкнутой системы управления

р) =

к1е

-%уР

Г (р) = X (р) Т2 р2 + 2ТЕр +1

(-43)е

-10 р

61,9 р2 + 6,9 р +1

где У(р) - изображение выходной переменной;

Х(р) - то же входной переменной; р = — -

—

оператор Лапласа как символ дифференцирования.

Второй этап моделирования (1996-2001 гг.) -этап стабилизации экономики Республики Беларусь - характеризуется тем, что государственная система управления начала функционировать более эффективно. Это позволило повысить устойчивость системы управления процессом стабилизации экономики, так как коэффициент демпфирования передаточной функции (4), описывающий динамику инерционного участка объекта управления этого этапа, стал равным единице Е = 1:

к2е-ур

Щр) = 2 2

Т2 р2 + 2Т&р+1

1е

-10 р

-10 р

46,24р2 +13,6 р +1 (6,8 р +1)2

10 /

где к2 = 1,0; 2Т1Е = 13,6, т. е. Е = — =

2Т

= 1,0.

2^46,24

Видно, что колебательное звено превратилось в инерционное звено второго порядка с равными постоянными времени 6,8 с.

Динамика третьего этапа (2001-2006 гг.) инерционного участка объекта - период активного роста номинального ВВП - характеризуется еще большей стабилизацией, так как описывается передаточной функцией вида

Ье-10 р

Шр) 2

116 р2 + 22,5 р +1 'Т2 р2 + 2Т£р +1'

где 2Т2Е = 22,5, т. е. Е =

= 1,052 > 1,0; к3 = 1,0.

22,5 22,5 22,5

2Т2 2л/Л6 2 -10,8

Видно, что инерционное звено второго порядка с равными постоянными времени превратилось на третьем этапе в инерционное звено

второго порядка с различными постоянными времени Т3 Ф Т4

-10 р

Щ( р) =-2-

3 116 р2 + 22,5 р +1 = 1е-10 р = (14,5р2 +1)(8р +1),

где к3 = 1; Т3 = 14,5 с; Т4 = 8 с, т. е. система государственного управления на этом этапе не только стабилизировала экономику, но и организовала активный рост ВВП.

При моделировании динамики инерционного участка объекта на четвертом прогнозном этапе развития экономики Беларуси (2006-2010 гг.) использовалась передаточная функция инерционного участка объекта с коэффициентом демпфирования ^ = 1,13

Ж4( р) =

к4е~10 р

-10 р

176 р2 + 30 р +1

,

30

30

= 1,13.

где к4 = 1А 5==2Ш,

При этом динамика опережающего участка объекта управления без лага на трех последних этапах аппроксимировалась передаточной функцией инерционного звена второго порядка

Кп (р) =

-0,43

Динамика инерционного участка объекта управления Жи0н (р)е Хур представлена передаточной функцией разомкнутой системы (3) для первого этапа моделирования, соответственно (4) - для второго, (6) - для третьего и (7) - для четвертого этапа моделирования динамических процессов тенденций изменения номинального ВВП при двух прогнозных сценариях социально-экономического развития Беларуси: инновационно-прорывном (с максимальным ростом ВВП за счет максимального потребления ТЭР при максимальном снижении энергоемкости ВВП) и инновационно-экологическом (при умеренном росте ВВП с минимальным ростом потребления ТЭР и умеренном снижении энергоемкости ВВП).

Передаточная функция крайнего внешнего возмущения Жв (р) принята равной

Ж (р ) =

к

-3

Тв р +1 3 р +1

где кв = -3 млрд дол. США - снижение объема ВВП в Республике Беларусь в 1999 г. после дефолта в России 1998 г.

Передаточная функция регулятора Жр

на втором и третьем этапах моделирования (1995-2006 гг.) соответствовала пропорционально-интегральному алгоритму принятия управленческих решений, а на четвертом этапе (2006-2010 гг.) - пропорционально-интегрально-дифференциальному (рис. 2).

При этом на втором этапе моделирования (1995-2001 гг.) применялась двухконтур-ная система управления: основной контур -

Ж (р) ^ (р) ^ ЖИ0Н (р)е ХуР ^ ГОС ^ Жр (р) (ГОС - главная обратная связь); внутренний контур - Жр(р) ^ Жоп (р) ^ Жн (р) ^ Жр (р), где Жнм (р) - неполная модель, реализованная в виде реального дифференциатора (элемента предвидения последствий принятия управленческих решений).

На третьем и четвертом этапах моделирования использовалась полная структурная схема (рис. 2), где в качестве неполной модели третьего этапа (2001-2006 гг.) применялось пропорциональное звено с коэффициентом передачи, равным единице, а в четвертом этапе - часть передаточной функции (7), которая не содержит лага.

Модель упредителя внутренних возмущений Жукв1 (р) принималась согласно принципу

инвариантности (независимости) равной обратной передаточной функции регулятора. Модель крайнего внешнего возмущения соответствовала реальному дифференцирующему звену с передаточной функцией

к Т

ЖуКв2

у Тд р +1

где кд - коэффициент передачи; Тд - время дифференцирования. В качестве управляющего воздействия анализировались инвестиции в основные фонды.

е

Упредитель /

I

1990-2001 2002-2005 Формирование задания

1995-2001 2001-2006 2006-2010 Полная модель

2001-2005

2001-2 0

Упредитель /2

2001-2010

Рис. 2. Полная структурная схема динамической модели системы управления ВВП Республики Беларусь

Результаты численного моделирования динамической модели эффективной системы управления номинальным ВВП Республики Беларусь при отработке крайнего внешнего возмущения (/2(0 - дефолт 1998 г. в России) приведены на рис. 3.

ВВП, млрд дол. США 3,5 3,0 2,5 2,0 1,5 1,0 0,5 0 -0,5 -1,0 -1,5 -2,0 -2,5 -3,0 -3,5

Рис. 3. Результаты численного моделирования динамической модели эффективной системы управления ВВП Республики Беларусь при отработке крайнего внешнего возмущения (дефолт в России в 1998 г.) при различных значениях времени принятия управленческих решений: 1 -при компенсации внешнего возмущения в виде дифференциатора первого порядка; 2 - то же второго порядка

Видно, что система управления с устройством компенсации крайнего внешнего возмущения в виде дифференциатора первого порядка инерционности отработала это возмущение весьма эффективно - за три года (кривая 1). Вместе с тем уменьшение времени принятия управленческих решений при использовании дифференциатора второго порядка инерционности приводит к повышению быстродействия в отработке крайнего внешнего возмущения на 30 %, т. е. на один год (кривая 2).

Применение упредителя внутренних возмущений ^укв1 (р) позволяет практически инвариантно отрабатывать последние (кривая 1, рис. 4), что существенно улучшает качество управления при минимальной динамической ошибке регулирования по сравнению с управлением без упредителя внутренних возмущений (кривая 2, рис. 4). При этом время отработки внутреннего возмущения сокращается примерно в два раза.

Адекватность предложенной математической модели фактическому изменению номинального ВВП Республики Беларусь подтверждается достаточно хорошим совпадением расчетного и фактического изменений ВВП до 2006 г. (рис. 5):--фактическая; - расчет.

х

зд

Рис. 4. Графики переходных процессов при единичном скачкообразном внутреннем возмущении динамической модели эффективной системы управления ВВП Республики Беларусь: 1 - с применением упредителя внутренних возмущений; 2 - без применения

ВВП, млрд дол. США

Год

Рис. 5. График переходных процессов расчетного, прогнозного и фактического изменений номинального ВВП Республики Беларусь: -- факт; — расчет;--------

прогноз; 1 - инновационно-прорывной сценарий; 2 - инновационно-экологический

Ускоренные темпы прогнозного роста ВВП четвертого этапа (2006-2010 гг.) достигаются за счет использования пропорционально-интегрально-дифференциального алгоритма принятия управленческого решения и уменьшения инерционных свойств упредителем внешних возмущений. При инновационно-прорывном сценарии (линия 1, рис. 5) ВВП Беларуси достигнет 78,3 млрд дол. США, при инновационно-экологическом сценарии (линия 2, рис. 5) -59,8 млрд дол. США. Это позволит в первом сценарии к 2011 г. поднять уровень среднемесячной заработной платы в республике до 700 дол. США, во втором - до 550 дол. США.

В Ы В О Д Ы

1. Распад Советского Союза, воспринимаемый как скачок входного воздействия на объект управления - экономику Республики Беларусь (на примере ВВП), позволил экспериментально определить динамические характеристики объекта в виде колебательного звена, которое при усилении системы государственного управления валовым внутренним продуктом, начиная с 1995 г., превратилось в инерционное звено с большей устойчивостью.

2. Моделирование различных вариантов принятия управленческих решений в системе государственного регулирования позволило установить, что для сохранения высоких темпов роста номинального ВВП Беларуси имеются потенциальные возможности.

3. Адекватность разработанной динамической модели фактическому изменению ВВП Республики Беларусь с 1990 по 2006 г. позволила сформировать два сценария прогноза роста ВВП до 2010 г. Реализация конкретного сценария будет зависеть от темпов снижения энергоемкости ВВП в условиях неуклонного роста стоимости природного газа.

Л И Т Е Р А Т У Р А

1. Никитенко, П. Г. Модель устойчивого социально-экономического развития Беларуси: проблемы формирования и эволюции / П. Г. Никитенко. - Минск: Право и экономика, 2000. - 312 с.

2. Национальная безопасность Республики Беларусь. Современное состояние и перспективы / М. В. Мяснико-вич [и др]. - Минск: Право и экономика, 2003. - 562 с.

0,65

0,60

0,55

0,50

США 0,45 0,40

ч о и и 0,35 0,30

0,25

ВВП, 0,20 0,15 0,10 0,05 0,00

-0,05

3. Кулаков, Г. Т. Общие принципы управления экономической безопасностью государства с точки зрения экономической кибернетики / Г. Т. Кулаков, Е. Г. Моисе-енко // Социально-экономическая безопасность государства: оценка, обеспечение: материалы междунар. науч.-практ. конф., Минск, 11-12 дек. 2001 г.: в 2 т. - Минск: ИНБ Республики Беларусь, 2003. - Т. 1. - С. 120-126.

4. Никитенко, П. Г. Динамическая модель системы управления валовым внутренним продуктом Республики Беларусь с прогнозом до 2010 г. / П. Г. Никитенко, Г. Т. Кулаков, В. М. Цилибина // Белорусская модель социально-экономического устойчивого инновационного

развития: формирование и пути реализации: сборник материалов междунар. науч.-практ. конф., Минск 19-22 апр. 2006 г. / Право и экономика; науч.-ред. совет: М. В. Мяс-никович [и др.]. - Минск, 2006. - С. 522-525.

5. Кулаков, Г. Т. Анализ и синтез систем автоматического регулирования: учеб. пособие / Г. Т. Кулаков. -Минск: УП «Технопринт», 2003. - 135 с.

6. Кузмицкий, И. Ф. Теория автоматического управления: учеб. пособие / И. Ф. Кузмицкий, Г. Т. Кулаков. -Минск: БГТУ, 2006. - 486 с.

Поступила 6.11.2007

УДК 658.51

ФОРМИРОВАНИЕ КОНЦЕПЦИИ УСТОЙЧИВОГО РАЗВИТИЯ ПРЕДПРИЯТИЯ

БАРМУЦКИЙ Р. И.

Белорусский национальный технический университет

Устойчивое развитие предприятия (УРП) представляет собой траекторию совершенствования производства и повышения эффективности управления, которая отличается от достигнутой точки равновесия экономического состояния и сложившейся структуры хозяйствования. Смена точки равновесия производственной системы заключается в смене траектории или направления ее совершенствования и устойчивого развития.

На современном этапе динамичного развития экономики Республики Беларусь предприятия становятся перед выбором: продолжить свое развитие в рамках индустриальных технологий или переходить на использование инновационных технологий. Новый информационный экономический уклад характеризуется быстрым устареванием капитала, а также определенными финансовыми потерями за риск, инновации, повышение качества продукции и проведение сервиса. Однако даже при такой экономической ситуации упущенная выгода вследствие неиспользования возможностей перехода на выпуск продукции нового поколения перекрывает минимизацию издержек, достигнутых за счет эффективного использования

устаревших технологий и изношенного оборудования.

Основным экономическим условием, способствующим УРП, является расширение рынка продаж продукции, своевременное обеспечение сырьем и материалами, научная организация труда и заработной платы, формирование финансовой стратегии и улучшения социального положения работников. Повышение эффективности производства и управления обеспечивает положительную обратную связь, которая стимулирует экономический рост, процесс развития рынка продаж и повышения заработной платы работников предприятий. Положительная обратная связь оказывает стимулирующее действие на увеличение результатов работы предприятия. Развитие рынков продаж выпускаемой продукции играет ведущую роль для устойчивого функционирования процесса производства.

На основе емкости рынка синхронизируется движение потоков материальных трудовых и денежных ресурсов, согласовываются интересы потребителя, производителя и поставщика материалов и комплектующих, определяется величина экономического эффекта и прибыли.