научный журнал «Актуальные исследования» #18 (45), май '21

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

В статье сравниваются дискретный и непрерывный варианты управления для односекторной модели экономики. Задача дискретного оптимального управления сводится к поиску условного экстремума функции многих переменных. Для решения системы нелинейных уравнений на каждом временном шаге применяется итерационный метод. Оптимальная траектория достигается вариацией начального значения множителя Лагранжа. Для получения оптимального значения целевой функции используется метод градиентного спуска.

Аннотация статьи
дискретное управление
условный экстремум
разностные уравнения
градиентные методы
целевая функция
Ключевые слова

Введение. Вопросам оптимального управления для моделей экономики посвящены работы многих исследователей [1,3-5]. В основном рассматривается непрерывное управление. Представляет интерес сравнение алгоритмов непрерывного и дискретного методов управлений.

Уравнение Солоу записывается в виде

,      .  (1)

В этом уравнении  – объем основных фондов в момент времени t,  – затраты на трудовые ресурсы. Параметрами уравнения являются: μ>0 – темп выбытия основных фондов, А – коэффициент, учитывающий уровень технического прогресса, 0≤σ≤1 – доля инвестиций на обновление основных фондов, α,β – параметры производственной функции, которая принята в виде функции Кобба-Дугласа:

.

На рис. 1(а) изображена схема финансовых потоков. Баланс финансовых потоков определяется равенством

,

где , .

Рис. 1. (а) Схема финансовых потоков; (б) Дискретное время

В [4] рассмотрена следующая задача: найти такое управление , при котором за время  капитал  вырастет от  до . При этом интегральный функционал

    (2)

примет минимальное значение.

Математическая формулировка непрерывной задачи имеет вид:

,     (3)

,   ,   (4)

.     (5)

Система уравнений Эйлера для расширенного функционала задачи (3)-(5) приводилась к следующему виду:

,   ,

.

Для поиска решения представленной системы уравнений применялась явная схема численного метода Эйлера.

1. Дискретная задача. На практике изменение управляющей переменной L(t) возможно только в определенные моменты времени и принятая величина не меняется некоторый конечный период времени.

При переходе к дискретной задаче выполняются преобразования:

1. Непрерывная область изменения времени заменяется дискретной, то есть [0,T] отрезок заменяется дискретным множеством точек (рис. 1(б)):

, .

2. Производная заменяется разностным отношением и, таким образом, получаем дискретную задачу

,    (6)

,   ,      (7)

 .        (8)

Необходимо найти такое управление , , при котором за время  капитал  вырастет от значения  до значения . При этом функционал (8) примет минимальное значение.

При условии , где  при , интеграл (8) можно преобразовать, тогда вместо (8) получим условие

 .      (9)

Заметим, что   – функция N переменных.  Переменные K1, K2,...,KN  содержатся в уравнении связи (6). Таким образом, сформулирована задача на условный экстремум функции многих переменных. Для поиска решения воспользуемся методом множителей Лагранжа [2, 6].

Составляем функцию Лагранжа:
.

Запишем необходимые условия существования экстремума

,   .

Приравнивая к нулю производные, получаем уравнения

 ,    (10)

 ,       (11)

 .        (12)

Кроме того, должны выполняться условия:

,    (13)

Рассмотрим алгоритм решения системы уравнений (10)-(12).

Из уравнения (10) выражаем Li-1≠0:

   (14)

При заданном значении λ0 находим L0.

Из (11) следует (при i=1,2,...,N-1)

 (15)

Из уравнения (12) находим

При известных значениях Ki, Li-1 и λi-1 строится итерационный процесс:

, ,   (16)

,   .  (17)

Решение задачи (10)-(13) зависит от выбора начального значения . Методом вариации данного начального значения находится такое значение , при котором . Для реализации метода вариации начального значения λ0 составляем целевой функционал

,  (18)

где ω  – постоянная.

Минимум функционала (18) определяется итерационным процессом:

,   .  (19)

2. Алгоритм решения дискретной задачи. Функциональная схема программы в системе Mathcad приведена на рис.2.

Рис. 2. Функциональная схема программы

В программе функция L(Kl) определяет оптимальное значение трудовых ресурсов (14); Dif(NfFλ0) предназначена для решения системы уравнений (10)-(12); f(KL ) – правая часть уравнения (12); F(KLl ) – правая часть уравнения (14); Iter(KLλ) – реализует итерационный процесс на каждом временном шаге для поиска решения системы (10)-(12); U(Nλ0 ,Yν) – вычисляет отклонение второго краевого условия от заданного значения. Управление методом градиентного спуска реализовано в блоках GRADU(… ) и GRAD(…). Результаты вычислений представлена на рис. 3, 4.

Рис. 3. (а) Оптимальный рост К: 1 – непрерывное управление; 2 – дискретное управление; (б) целевая функция в итерационном процессе

Оптимальное управление представлено на рис. 4. Следует отметить, что с увеличением числа шагов разность между непрерывным и дискретным управлением уменьшается.

Рис. 4. Оптимальное управление: 1 – непрерывное, 2 – дискретное; (аN=2 (бN=4

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

Таблица

Отношение оптимальных затрат в зависимости от количества промежутков

Число временных промежутков

при дискретном управлении

2

4

10

20

Отношение оптимальных затрат

0,834

0,926

0,979

0,947

Основные выводы

  1. Составлен алгоритм решения задачи дискретного оптимального управления для односекторной модели экономики. Построен итерационный процесс решения системы нелинейных алгебраических уравнений на каждом временном шаге.
  2. Разработана программа в системе Mathcad, реализующая алгоритмы дискретного и непрерывного управлений. Программа проверена в широком диапазоне параметров.
  3. Сравнительный анализ решений дискретного и непрерывного оптимальных управлений продемонстрировал существенное расхождение оптимальных решений при малом числе шагов для дискретного варианта. Оптимальные затраты трудовых ресурсов оказались меньше на 10…15% при дискретном управлении.
Текст статьи
  1. Горбунов В.К., Львов А.Г. Эффективные производственные фонды и производственные функции малого предпринимательства регионов // Экономика региона. 2018. Т. 14, № 2. С. 502–515.
  2. Демидович Б.П., Марон И.А. Основы вычислительной математики: Учебник для вузов. М.: Наука, 1970. 664с.
  3. Колемаев В.А. Математическая экономика: Учебник для вузов. – 2-е изд., перераб. и доп. М.: ЮНИТИ-ДАНА, 2002. 399с.
  4. Охорзин В.А. Прикладная математика в системе Mathcad: учебное пособие. – 2-е изд., испр. и доп. СПб. «Лань», 2008. 352с.
  5. Слиденко А.М., Агапова Е.А. Численные алгоритмы оптимизации моделей управления экономикой предприятия // Вестник ВГАУ, Т. 12, 3(62), 2019. С. 192-205.
  6. Эльсгольц Л.Э. Дифференциальные уравнения и вариационное исчисление. М.: Наука, 1969. 424 с.
Список литературы
Ведется прием статей
Прием материалов
c 16 октября по 22 октября
Осталось 3 дня до окончания
Публикация электронной версии статьи происходит сразу после оплаты
Справка о публикации
сразу после оплаты
Размещение электронной версии журнала
26 октября
Загрузка в eLibrary
26 октября
Рассылка печатных экземпляров
03 ноября