Технические науки /12. Автоматизированные системы управления на производстве

К.э.н. Сизиков А.П.

Самарский государственный экономический университет, Россия

АЛГОРИТМ ОПТИМИЗАЦИИ ДВУХУРОВНЕВОЙ ПРОИЗВОДСТВЕННОЙ СИСТЕМЫ

 

Рассматривается двухуровневая модель некоторой системы. Пер­вый (верх­ний) уро­вень представлен линейной моделью с пере­мен­ными параметрами. Это, если речь идет о производственно-техно­ло­ги­ческих объектах, как правило, мо­дель «затраты - выпуск» и/или пре­д­­с­тав­ленная в матричном виде сетевая пото­ко­вая модель. Соот­ветст­венно, пере­менные име­ют смысл интен­сив­ностей «тех­но­ло­гий» и/или материальных пото­ков. К этим переменным можно с помо­щью соот­ветствующих коэф­фи­цие­нтов (в общем слу­чае перемен­ных) при­вя­зать лю­бые производственные фак­торы и, соот­ветст­вен­но, сфор­му­ли­ровать целе­вую функ­цию всей системы. Для удобства даль­ней­шего изло­жения запишем модель верх­него уровня в обоб­щен­­ной канонической форме:

.                                              (1)

Вычислительная схема строится с использованием метода гене­ри­­ро­вания столб­­цов [1]. Столб­цы гене­ри­­­ру­ются в ре­зуль­тате реше­ния локальных задач (задач второго уровня). Модель центра играет роль коор­ди­натора локальных решений. Эту функцию она осу­щест­вляет путем рас­чета симп­лекс­ных муль­ти­пли­каторов, кото­рые трак­ту­ются в данном слу­чае как «цены» ингредиентов (ресурсов, продуктов) сис­­те­мы. Каж­дая локальная модель с уче­том этой инфор­мации опти­ми­зи­рует собст­вен­ную «тех­но­­логию» и соот­­ветст­венно коэф­фи­циен­­ты пред­­став­­ляющего ее столб­ца «зат­ра­ты - выпуск» и передает их в базо­вую модель.  

Каж­дый столбец модели верхнего уровня является пред­ста­ви­телем (агре­га­том) модели соот­вет­ст­вую­щего элемента, или, иначе гово­ря, вектор-столбцом , , где - мно­жест­во воз­мож­ных «технологий» узла. Решение обобщенной задачи сво­дит­ся к формированию опти­­маль­ного ба­зи­са путем гене­ри­рования столбцов, являю­щих­ся край­­ни­­ми точками допус­ти­мых мно­­жеств тех­нологий.

Алгоритм:

1. Формируем начальный вариант задачи, включая в матрицу усло­вий по одно­му пред­ста­ви­телю множеств , . Находим условно-опти­маль­ный базис  и соот­вет­ст­вую­щую ему обрат­ную мат­­рицу.

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

3. Генерируем столбцы с минимальными оценками замещения, т.е. нахо­­дим реше­ния локальных задач

, .                           (2)

4. Проверяем условие . Если оно выполняется, то теку­щий базис оптимален, решение получено. В противном случае переходим в следую­щий пункт.

5. Вводим в базис столбец . Формируем новую мат­рицу . Стол­бец, выводимый из базиса, уда­ляем из задачи. Пере­хо­­дим в п. 2.

Эту вычислительную схему можно проиллюстрировать на при­мере опти­мизации сис­те­мы, состоящей из  взаи­­мосвя­­­зан­ных произ­водств. Согласно Леон­ть­еву для нее мож­но запи­сать урав­нения ма­те­­­риального баланса:

,                                                        (3)

где  - матрица коэффициентов прямых материальных затрат;  - вектор вало­во­го продукта;  - вектор конечного продукта.

Следуя неоклассической теории, будем полагать, что -й про­из­водст­вен­ный объект опи­са­н степенной производственной функ­цией:

,                                                 (4)

где - масштабный множитель;  - коэффициент затрат труда;  - коэф­фи­ци­­ент затрат -го про­дук­та;  - соот­­ветст­вую­щий коэф­фи­ци­­ент эластичности, , , .

Деля обе части (4) на , получаем:

.        (5)

Таким образом, из (5) следует, что в балансовой модели -й объект может быть пред­­став­лен столб­­цом пара­мет­ры, кото­­рого при­над­­лежат мно­жеству ре­ше­ний по  урав­не­ния

.                                                 (6)

где ,   - коэффициент прямых затрат -го продукта .

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

                                                      (7)

где

.                                            (8)

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

.                       (9)

где , .

Функция Лагранжа:

.                   (10)

Условия стационарности:

      

      

Решение:

.                              (11)

Далее идет цикл: 1) расчет муль­­типликаторов по формуле  и 2) расчет тех­но­логических коэф­фици­ен­тов по фор­­му­ле (11). И так до получения прием­ле­мого по точ­ности ре­зу­ль­та­та. Рас­че­ты пока­зывают, что процесс моно­тон­но сходится за число итераций, соиз­­ме­ри­мое с раз­мер­ностью за­да­чи.

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

 

Литература:

1.     Лэсдон Л. Оптимизация больших систем. Главная редакция физико-математической литературы изд-ва «Наука», М., 1975.