Ф. X. Асхакова
ВЕКТОРНАЯ ОПТИМИЗАЦИЯ В БАЛАНСОВОЙ МОДЕЛИ ЛЕОНТЬЕВА - ФОРДА, УЧИТЫВАЮЩЕЙ УТИЛИЗАЦИЮ ВРЕДНЫХ
ОТХОДОВ
Работа представлена кафедрой информатики и вычислительной математики Карачаево-Черкесского государственного университета им. У. Д. Алиева.
Научный руководитель - доктор физико-математических наук, профессор Е. А. Семенчин
В статье в рамках модели Леонтьева - Форда, учитывающей утилизацию вредных отходов производства, поставлена и исследована задача одновременной максимизации общего объема производственной продукции и минимизации общего объема экологически вредных отходов производства, подлежащих утилизации. Приведены конкретные примеры решения данной задачи с помощью офисной программы Microsoft Excel.
The author of the article states and investigates the problem of simultaneous maximization of total amount of industrial production and minimization of total amount of ecologically hazardous waste liable to recycling within the limits of the Leontiev-Ford model, which considers utilization of hazardous production waste. Certain examples of this problem’s solution are given by means of the office programme Microsoft Excel.
1. Постановка задачи. Экономико-математическая модель Леонтьева - Форда, учитывающая утилизацию вредных отходов производства, имеет вид1 :
х = Аи х + А12 у + Ь1 - ^13 у , (1)
у = А21х + А22у - Ь2 + А23У, (2)
х >в, у >в , (3)
где х е Яп - вектор валового выпуска продукта, называемого полезным, Я" - п - мерное вещественное пространство,
у е Ят - вектор вредных отходов, выбрасываемых в окружающую среду, которые возникают в процессе производства и подлежат уничтожению с целью поддержания необходимого уровня экологического состояния окружающей среды;
Ъ1 е Я" - вектор, характеризующий величину (объем) спроса на выпускаемый продукт;
Ъ2 е Ят - вектор, характеризующий остаточный уровень вредных отходов, т. е. отходов, которые не могут быть ликвиди-рованы;
А11 - технологическая матрица размера п х т (матрица прямых затрат);
А12- матрица размера п х т, характеризующая затраты при уничтожении вредных отходов;
А13- матрица размера п х т, характеризующая затраты при утилизации отходов;
А21 - матрица размера п х т, характеризующая объем вредных отходов, получаемых при выпуске полезного продукта;
А22 - матрица размера п х т, характеризующая объемы вновь получаемых вредных веществ при уничтожении старых;
А23 - матрица размера п х т, характеризующая объемы вновь получаемых вредных веществ при утилизации старых.
в - нулевой вектор (требуемой размерности: либо п, либо т).
Модель (1)-(3) обобщает известную модель Леонтьева, которая в наших обозначениях имеет вид2.
х = Ап х + А12 у + Ъ1, х >в
и призвана ответить на вопрос3: «... можно ли в условиях данной технологии удовлетворить конечный спрос?» Модель (1)-(2) призвана ответить на вопрос: можно ли в условиях данной технологии, учитывающей загрязнение окружающей
среды вредными отходами производства и возможность утилизации вредных отходов, удовлетворить конечный спрос b1 при заданном (допустимом) остаточном уровне вредных отходов b2.
Модель (1)-(2) впервые появилась в работе4 как результат обобщения известной модели Леонтьева - Форда5
x = A11 x + A1 2 y + b1 ,
У = A2ix + A22 У - b2 .
В работе6 авторов данной статьи предложена методика построения неотрицательного решения модели (1)-(3) «методом простой итерации и ее программная реализация “The productivity of model”»7.
В данной работе в рамках модели (1)-(3) ставится и изучается задача: найти
П
F (x ) = £ xi ^ max , (4)
i=1
m
F2 (y) = Z yi ^ min (5)
i=1
в предположении, что x e Rn, y e Rm удовлетворяют условиям (1)-(3).
П
Очевидно, в данной задаче ZX -
i=1
общий объем продукции, производимой
m
всеми п отраслями, ^ у - это общий объем
i=1
вредных отходов, возникающих в процессе производства и подлежащих утилизации.
2. Методы решения изучаемой задачи. Поставленная в п. 1 задача (4), (5), (1)-(3) представляет собой векторную (двухкритериальную) оптимизационную задачу с критерием
(F1 (4 -F2 (у)), (6)
которой требуется максимизировать при ограничениях (1)-(3). Для решения подобных задач удобно прибегать к методам свертки векторных критериев.
2.1. Линейная сверстка критерия (6). Вместо векторного критерия (6) можно
рассмотреть скалярный c1F1(x) — c2F2(y), Cj = const > 0, c2 = const > 0, c1 + c2 = 1 в котором c1, c2 задаются (определяются) экспертами. Тогда двухкритериальная задача (4), (5), (1)-(3) сводится к задаче линейного программирования:
n m
c i Zx -c 2 Z у ^ max (7)
i=1 i =1
при ограничениях (1)-(3).
Данная задача может быть решена с помощью надстройки Microsoft Excel/Поиск решения.
2.2. Квадратичная свертка критерия (6). Допускаем, что задачи (4), (1)-(3) и (5), (1)— (3), представляющие собой по отдельности задачи линейного программирования, имеют решения. Тогда для свертки критерия (6) можно предположить следующий метод.
1) находим решение
(x(i), у»>)=(хг,...,xi1, у,01,...,ym)
задачи (4), (1)-(3);
2) находим решение
(x<2->, у<2>)=(*р>, ...,Х<;>, ypi,...,ym') задачи (5), (1)-(3);
3) строим скалярный критерий
2 2 n / m , sf
F (x, y)=£ (x- xi") (y, - yp>). (8)
i = 1 i = 1
Тогда двухкритериальная задача (4), (5), (1)-(3) сводится к задаче квадратичного программирования8
2
n / \ m / ч-
Z (x,-x!") +£(у, -уГ-')2 ^ min (9)
i = 1 i = 1
при ограничениях (1)-(3).
Задача квадратичного программирования (9), (1)—(3) может быть решена с помощью надстройки Microsoft Excel/Поиск решения.
п р имер 1. Пусть
А 1 =
(0,49 0,15 0,4 ^
0,1 6 0,09 0,03
,0,008 0,07 0,04у
А1 2
(0,0003 0,0003^ 0,0003 0,0 0 0 3 0,0004 0,0002
А13
(0,0001 0,0001^1
0 0, 0 0 0 1
, Ъ1 =
( 90 ^ 79
V 40 у
А21
(0,001 0,004 0,001^1 (0,07 0,003^
’ ’ , 42 = ’ ’
0,002 0,003 0,001, 0,001 0,002
V7 7 У V7 7у
А23
(0,01 0,001^1 Г0.34^
0,001 0
Ъ2 =
V0'23/
с1 = 0,3 с2 = 0,7
Находим решение задачи (7), (1)-(3): (257,9^
х =
13 3,9 V 53,6 у
, у =
' 0,6л V0,7 У
, значение целевой функции (7) равно 62,3.
П ример 2. Пусть Аи, А12, А13, А21, А22, А2Ъ, Ъ1, Ъ2 имеют тот же вид, что и в примере 1.
Решая задачу (4), (1)-(3), находим (257,9^
133, 9 V 53,6 у
у(1 > =
( 0>|
V °у
значение целевой функции (4) равно 445,4.
Решая в нашем случае задачу (5), (1)-(3), находим
( 0 ^
84,97
0
, значение целевой функции (5) равно 0,025.
Строим критерий
(х-х«)2-(уг -у(2))2 = (х - 257.9)2 + (х2 -133.9)2 + (х3 -53.6)2 + (у1 -0)2 + (у2 -0.025)2
(т. е. критерий (8)) и находим его минимум при ограничениях (1)-(3):
(257,9^
х =
1 3 3,9
V 53,6 ,
у =
, значение целевой функции (9) равно 43683,26.
Сравнивая решения примеров 1 и 2, убеждаемся (и это вполне естественно), что они совпадают.
ПРИМЕЧАНИЯ
1 Данилов Н. Н. Курс математической экономики: Учеб. Пособие. М.: Высшая школа, 2006. С. 239.
2 Там же. С. 239.
3 Там же. С. 239.
4 Омарова А. Д. Модель, учитывающая возможности утилизации вредных отходов. Вузовская наука - Северо-Кавказскому региону: Материалы IV per. науч.-тех. конф. Ставрополь: СевКав-ГТУ, 2000. С. 37.
5 Леонтьев В., Форд Д. Межотраслевой анализ воздействия структуры экономики на окружающую среду // Экономика и математические методы. Т. VIII, 1972. Вып. 3. С. 370-399.
6 Семенчин Е. А., Асхакова Ф. X. Методика построения численными методами решения балансовой модели Леонтьева-Форда, учитывающей утилизацию отходов. Нелокальные краевые задачи и родственные проблемы математической биологии, информатики и физики: Материалы III Международной конференции. Нальчик. 2006. С. 257-259.
7 Программный продукт «Комплекс программ “The productivity of model”». Свидетельство об официальной регистрации программы на ЭВМ регистрационный номер № 2007611644. Зарегистрировано в Реестре программ для ЭВМ от 19.04.2007.
8 Кузнецов Ю. Н., Кузубов В. И., Валощенко А. Б. Математическое программирование. М.:
Высшая школа, 1976. С. 216-248.