Задача оптимального управления ресурсами промышленного предприятия с учетом взаимодействия со см

На правах рукописи







Старинец Дмитрий Владимирович





Задача оптимального управления ресурсами промышленного предприятия с учетом взаимодействия со смежными предприятиями.



			   			

Специальность 05.13.01. –	Системный анализ, управление и 

обработка информации (промышленность).







АВТОРЕФЕРАТ

диссертации на соискание ученой степени

кандидата физико-математических наук















Москва - 2009

Работа выполнена в Московском физико-техническом институте (Государственный университет)



Научный руководитель:

доктор физико-математических наук, профессор

Дикусар Василий Васильевич

Официальные оппоненты:

доктор физико-математических наук 

Березнев Валентин Александрович



кандидат физико-математических наук	

Дарьин Александр Николаевич

Ведущая организация: Центральный экономико-математический институт РАН



Защита диссертации состоится   02  апреля 2009г. в   14       часов на заседании диссертационного совета Д002.017.03 при Учреждении Российской академии наук Вычислительный Центр им. А.А. Дородницына РАН по адресу: 119991, г. Москва, ул. Вавилова, д. 40 в конференц-зале. 



С диссертацией можно ознакомиться в библиотеке Вычислительного центра им. А.А. Дородницына РАН



Автореферат разослан                                           2009 г.



Ученый секретарь совета по защите докторских

 и кандидатских диссертаций 
кандидат физико-математических наук 			                 Мухин А.В.





ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ



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

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

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



Цель работы. Целью работы является:

построение модели взаимодействия промышленных предприятий в условиях кризиса производства; 

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

решение задач линейного программирования большой размерности методом продолжения решения по параметру; 

на основании проведенных исследований – предоставить возможность выработки обоснованных эффективных управленческих решений для оптимального развития промышленного производства в условиях кризиса.

Методы исследования. Основным инструментом для решения поставленных задач является принцип максимума (схема Дубовицкого-Милютина) и метод продолжения решения по параметру. Поставленные задачи (за счет дискретизации обыкновенных дифференциальных уравнений) сводятся к задачам линейного и нелинейного программирования большой размерности. Применение принципа максимума в дискретном варианте сводит первоначальную задачу к задаче линейного программирования большой размерности. В качестве параметра выступает время. Это позволяет сначала на малом отрезке решать задачу малой размерности и затем полученное приближение используется при его продолжении по параметру.



Научная новизна.

Решена новая важная задача эффективного управления ресурсами с учетом  взаимодействия двух промышленных предприятий в условиях кризиса производства.

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



Обоснованность научных положений. 

Теоретические положения и выводы диссертации сформулированы в виде утверждений и теорем, которые строго доказаны.



Практическая ценность.

Модели, методы и алгоритмы, разработанные в диссертации, применялись для решения практических задач взаимодействия промышленных предприятий, а также в учебном процессе в Московском Физико-Техническом институте и в Вычислительном Центре РАН. Предложенные методы продолжения решения по параметру, а также методы регуляризации вырожденных задач могут быть использованы в теоретических исследованиях при решении прикладных задач оптимального управления. Был  адаптирован  пакет прикладных программ БАЛАНС-2 для решения задачи ЛП и использован  для практических численных расчетов показателей эффективности производства на модельном примере (с применением метода продолжения решения по параметрам). 



Апробация работы.

Основные положения исследования докладывались и обсуждались на международной конференции в Черногории (International Conference «Nonlinear Analysis and Optimization Problems», Montenegrin Academy of Sciences and Arts, Petrovac, Montenegro, October 06th – 10th , 2008), на 14-ой Байкальской школе-семинаре СО РАН «Методы оптимизации и их приложения». (Иркутск-Байкал 2-8-го июля 2008г.) и на научных семинарах в МФТИ и в ВЦ РАН.

Личный вклад. 

1) Проведен качественный и количественный анализ задачи эффективного управления взаимодействием двух промышленных предприятий.

2) Разработан прямой численный метод построения  гипотезы по определению множества активных индексов для задачи управления с ограничениями типа неравенств (геометрия оптимальной траектории).

3) Предложена регуляризация вырожденного случая принципа максимума.

4)   Разработан явный эффективный численный метод решения жестких систем ОДУ.

5) Автором адаптирован  пакет прикладных программ БАЛАНС-2, использование которого позволяет выработать обоснованные управленческие решения.

Публикации. Основные результаты исследования отражены в восьми публикациях. Список работ приведен в конце автореферата. В совместных с соавторами работах автору принадлежит 50% результатов.

Структура и объем работы. 

Диссертация состоит из введения и четырех глав и двух приложений. 







ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

Во «Введение» приведены цели исследования, актуальность, обзор литературы по данному вопросу и основные результаты, выносимые на защиту. Дана характеристика научной новизны, практической значимости и указаны апробации полученных результатов.

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

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

Динамическая модель взаимодействия двух предприятий имеет вид:

Задача А:

                                                            

а ограничения на управления и фазовые переменные соответственно

                                                     

                     ,                  

                       найти    при условиях

                                             



Здесь:

  - общий объем ресурсов, передаваемый от предприятия –   

 донора предприятию – акцептору;

     - интенсивность поставки ресурсов (управление);

     - суммарный объем производства промышленной продукции;

      - интенсивность промышленного производства (управление);

       - прибыль от реализации произведенной продукции;

       Остальные переменные являются параметрами  модели.  

Доказаны теоремы:

1. Если параметры и управления задачи А измеримы и для правых частей  дифференциальных уравнений выполнены условия Филиппова А.Ф. и   существует хотя бы одна допустимая пара, удовлетворяющая всем условиям задачи А, то оптимальное решение существует и единственно. 

2.  Принцип максимума для задачи А выполняется тривиально.     

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





График 1

На графиках 1 и 2 приведены характерные решения для фиксированных параметров  и . Приведенные графики  построены с помощью графического пакета и адаптированного пакета прикладных программ  БАЛАНС — 2. 

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

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





График 2





Для рассматриваемой модели в качестве примера параметр-функции задавались в виде: 

	

Где  , ,  - константы.

График 2 иллюстрирует динамику управляющих функций (для модельного примера). Значения на оси ординат слева характеризуют интенсивности, а справа — параметр-функции. А по оси абсцисс отложено характерное время.

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

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

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

,



при условиях

,





;



,





необходимые условия оптимальности имеют вид:

, ,



,



,



где вектор  является решением системы уравнений





с условиями

.



Далее рассматриваются достаточные условия оптимального управления, основанные на методике сведения линейной задачи ОУ со смешанными ограничениями к задаче ЛП в банаховых пространствах, предложенной А.М. Тер-Крикоровым. Далее рассматриваются две задачи:

Задача  1.

Найти управления , дающие максимум линейному функционалу



(1)

при следующих ограничениях:

,

(2)



, ,

(3)



, .

(4)



Матрицы , ,  и  и векторы ,  имеют ограниченные измеримые компоненты, которые выражают обобщенные    технологические и весовые  показатели. Соответствующие матрицы и векторы имеют следующие размеры: , , , , , , , , . Векторы с символом  являются строками, без – столбцами.

Задача  2.

Найти управления , , дающие минимум линейному функционалу



(5)

при следующих ограничениях:

,

(6)

,

(7)

, .

(8)



Достаточные условия оптимальности задач 1 и 2 даются следующей теоремой:

Теорема 1 (Тер-Крикоров). Пусть для некоторых допустимых управлений  и ,  задач 1 и 2 выполнены условия

, ;

(9)

, ;

(10)

, ,

(11)

причем первые два равенства выполняются почти при всех . Тогда ,  будет оптимальным решением задачи 1, а , ,  будет оптимальным решением задачи 2.

Необходимые условия оптимальности для задачи 1 формулируются в терминах принципа максимума Понтрягина с использованием сопряженных переменных . Связь сопряженных переменных  и переменных задачи 2 дается следующими леммами:

       Лемма 1. Если при допустимом управлении задачи 1 существует вектор сопряженных переменных , константа  и векторы множителей Лагранжа , , удовлетворяющие дифференциальным уравнениям и краевым условиям для , условиям Блисса и условиям дополняющей нежесткости для , , то  и  являются допустимыми управлением и фазовым вектором задачи 2.

      Лемма 2. Если существуют допустимые управления , ,  задач 1 и 2, и они удовлетворяют условиям (2.5.9)-(2.5.11), то вектор траектории  задачи (2.5.5)-(2.5.8), соответствующей управлению , является вектором сопряженных переменных  задачи (2.5.1)-(2.5.4) при .

На основании лемм 1 и 2  теорема 1 переформулируется следующим образом:

Теорема 2. Если при данном допустимом управлении  задачи 1 существуют число , кусочно-гладкая вектор-функция , измеримые вектор-функции ,  и вектор  такие, что выполняются условия (2.5.12)-(2.5.15), то  – оптимальное управление задачи 1.

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

Третья глава посвящена вопросу нахождения первого приближения геометрии оптимальной траектории при смешанных ограничениях, типа неравенств.

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

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

Суть рассматриваемой схемы выделения множества активных ограничений заключается в дискретизации времени и сведении исходной задачи 1-2 к вспомогательной задаче математического программирования с конечным числом переменных. Дифференциальные уравнения при этом заменяются конечно-разностными по схеме Эйлера первого или второго порядка точности. Подобные задачи рассмотрены в трудах Ю.Г. Евтушенко. Решение данной вспомогательной задачи рассматривается как некоторое приближение к решению исходной, и на его основании производится выделение подмножества активных ограничений.

В четвертой главе излагаются различные формы задач линейного программирования (ЛП), куда входят также несобственные задачи. Здесь для решения задачи ЛП предлагается метод введения параметра в целевую функцию. Это дало возможность получить эффективную оценку решения задачи ЛП. Кроме того использовался адаптированный пакет прикладных программ БАЛАНС - 2, обеспечивающий многократное формирование условий нахождения решения и создания необходимых для анализа выходных файлов. Была  использована реализация для ОС Windows 2K-XP базовой версии алгоритма анализа неполных математических моделей (разработанная в 1985 году в IIASA, в рамках проекта Regional Development, на языке "Fortran-IV" для ПЭВМ Altus-2. Авторы: Ким К.В. и Умнов А.Е.), адаптированная для языка С++ на кафедре высшей математики МФТИ в рамках совместных исследований с ЗАО «Оптимизационные системы и технологии». В комплекс программных средств решения задач ЛП были включены модули диагностики и анализа качества (получаемых на основе найденных решений) гипотез об оптимальной геометрии фазовых траекторий. Специальные программные средства были разработаны для решения сопряженных задач, проверки формализма Понтрягина-Дубовицкого-Милютина и прямой проверки оптимальности решения на множестве допустимых вариаций.

В приложении 1 изложен новый эффективный метод интегрирования жестких систем обыкновенных дифференциальных уравнений на базе параметризации явных схем.

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

Выводы:

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

	Предложены  явные численные методы для интегрирования жестких систем обыкновенных дифференциальных уравнений. 

	Предложен метод оценки геометрии оптимальной траектории. 

          Предложен метод регуляризации вырожденного принципа максимума в задаче взаимодействия двух промышленных предприятий.

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



Основные результаты диссертации опубликованы в работах:

Дикусар В.В., Старинец Д. В.Управление риском портфеля ценных бумаг  Труды ИСА РАН. Т.31(1) Динамика неоднородных систем. 2007г.С.14-22.

Старинец Д. В. Методы продолжения при решении краевых задач оптимального управления. Труды ИСА РАН. Т.31(1) Динамика неоднородных систем. 2007г.С.74-80.

Дикусар Э.В., Чекарев Д.А., Старинец Д.В.  Достаточные условия экстремума в линейной задаче оптимального управления. Труды ИСА РАН. Т.32(1) Динамика неоднородных систем. 2008г С. 16-23.

Дикусар Э.В., Чекарев Д.А., Старинец Д.В. Сходимость дискретных аппроксимаций. Труды ИСА РАН. Т.32(1) Динамика неоднородных систем. 2008г. С. 101-110.

Дикусар Э.В., Чекарев Д.А., Старинец Д.В. Численно-аналитический метод решения задач оптимального управления со смешанными ограничениями. Труды ИСА РАН. Т.32(1) Динамика неоднородных систем. 2008г.С. 111-122.

Дикусар В.В., Старинец Д.В. Методы интегрирования жестких систем явными методами. Труды 14-ой Байкальской школы-семинара «Методы оптимизации и их приложения». Иркутск-Байкал 2-8-го июля 2008г. т.3 ИСЭМ СО РАН 2008. С. 77-85.

Dikusar V.V., Starinets D.V. Continuation methods for solving boundary value problems. Abstracts of  International Conference «Nonlinear Analysis and Optimization Problems», Montenegrin Academy of Sciences and Arts, Petrovac, Montenegro, October 06th – 10th , 2008. P.37.

 Dikusar V.V., Starinets D.V. Determined portfolio dynamic problem. Abstracts of  International Conference «Nonlinear Analysis and Optimization Problems», Montenegrin Academy of Sciences and Arts, Petrovac, Montenegro, October 06th – 10th , 2008. P.38.



2