--PAGE_BREAK--
Суммарные затраты:
f(x) = 6´18+7´2+10´13+8´32+8´5+8´40+7-23+14-7=1127
Положим U1=0
V2 = U1+C12=7=U4+10 Þ U4 = -3
V3 = U4+8=5; V4=U1+8=8=U2+8 Þ U2=0
V5 = U1+7= 7 = U3+14 Þ U3= -7
V1 = U3+6= -1
dij = (Ui+Cij)-Vj
9 0 9 0 0
(dij) = 11 1 5 0 8
0 -3 -2 -3 0
14 0 0 1 6
Наличие отрицательных оценок свидетельствует о том, что план не является оптимальным. Построим контур перераспределения для клетки (3,2). Наименьшая поставка в вершине контура со знаком “-” равна 2. Произведем перераспределение поставок. Результаты представим в таблице 3.
Мощности
постав-
щиков
140
Мощности потребителей
U i
18
15
32
45
30
30
10
7
14
8/5
7/25
0
40
12
8
10
8/40
15
0
25
6/18
10/2
10
12
14/5
-7
45
16
10/13
8/32
12
16
-7
Vj
-1
7
5
8
7
Суммарные затраты:
f(x) = 6´18+10´2+10´13+8´32+8´5+8´40+7´25+14´7=1119
Положим, U1=0 Þ V4=8, V5=7; V4=U2+8 Þ U2=0
V5 = U3+14 Þ U3= 7-14= -7; V1= -7+6= -1; V2= -7+10= +3
V2=U4+10 Þ U4=3-10= -7; v3= -7+8=1
9 4 13 0 0
(dij) = 13 5 9 0 8
2 0 2 -3 0
10 0 0 -3 2
Наличие отрицательных оценок свидетельствует о том, что план не является оптимальным. Построим контур перераспределения для клетки (3,4).
Наименьшая поставка в клетке со знаком “-” равна 5. Произведем перераспределение поставок результаты представим в таблице 4.
Мощности
постав-
щиков
140
Мощности потребителей
U i
18
15
32
45
30
30
10
7
14
8
7/30
0
40
12
8
10
8/40
15
0
25
6/18
10/2
10
12/5
14
-4
45
16
10/13
8/32
12
16
-4
Vj
2
+6
4
8
7
Суммарные затраты:
f(x) = 7´30+8´40+6´18+10´2+12´5+10´13+8´32=1104
U1=0 Þ V5= 7; U2=0 Þ V4=8=U3+12 Þ U3=-4 Þ
V1= 6-4=2, V2=10-4=+6=U4+10; V3= -4+8= +4
8 1 10 0 0
(dij) = 10 2 6 0 8
0 0 2 0 3
10 0 0 0 5
Матрица оценок (dij) не содержат отрицательных величин Þ данный план является оптимальным, т.к. С34 = 0, а клетка (3,4) не является запятой, то данный план не является единственным. Стоимость перевозок по этому плану, как было рассчитано ранее, равна f(x) = 1104.
3.6. Симплекс-метод решения задач линейного программирования.
Симплекс-метод позволяет отказаться от метода перебора при решении задач линейной оптимизации, является основным численным методом решения задач линейного программирования и позволяет за меньшее число шагов, чем в методе перебора, получить решение.
Реализация алгоритма симплекс-метода.
1. Записать задачу в канонической форме: заменить все ограничения-неравенства с положительной правой;
2. Разделить переменные на базисные и свободные: перенести свободные переменные в правую часть ограничений-неравенств.
3. Выразить базисные переменные через свободные: решить систему линейных уравнений (ограничений-неравенств) – относительно базисных переменных;
4. Проверить неотрицательность базисных переменных: убедиться в неотрицательности свободных членов в выражениях для базисных переменных. Если это не так, вернуться к пункту 2, выбирая другой вариант разделения переменных на базисные и свободные.
5. Выразить функцию цели через свободные переменные: базисные переменные, входящие в функцию, выразить через свободные переменные;
6. Вычислить полученное базисное решение и функцию цели на нем: приравнять к 0 свободные переменные;
7. проанализировать формулу функции цели: если все коэффициенты свободных переменных положительны (отрицательны), то найденное базисное решение будет минимально (максимально) и задача считается решенной;
8. Определить включаемую в базис и исключаемую из базиса переменные: если не все коэффициенты при свободных переменных в функции цели положительны (отрицательны), то следует выбрать свободную переменную, входящую в функцию цели с максимальным по модулю отрицательным (положительным) коэффициентом, и увеличивать ее до тех пор, пока какая-нибудь из базисных переменных не станет равной 0. Свободную переменную рассматриваем как новую базисную переменную (включаемую в базис), а базисную переменную рассматриваем как новую базисную переменную (исключаемую из базиса);
9. Используя новое разделение переменных на базисное и свободное, вернуться к пункту 3 и повторять все этапы до тех пор, пока не будет найдено оптимальное решение.
В заключение отметим, что определение оптимального решения распадается на два этапа:
· Нахождение какого-либо допустимого решения с положительным свободным членом;
· Определение оптимального решения, дающего экстрему целевой функции.
IV. Методы нелинейного программирования.
4.1. Основные понятия, постановка и методы решения задачи нелинейного программирования.
Нелинейное программирование (планирование) – математические методы отыскания максимума или минимума функции при наличии ограничений виде неравенств или уравнений. Максимизируя (минимизируя) функция представляет собой принятый критерий эффективности решения задачи, соответствующий поставленной цели. Он носит название целевой функции. Ограничение характеризует имеющиеся возможности решения задачи.
Целевая функция или хотя бы одно из ограничений нелинейное (т.е. на графиках изображается не прямыми-кривыми-линиями) существо решения задач нелинейного программирования заключается в том, чтобы найти условия, обращающие целевую функцию в минимум или максимум. Решение, удовлетворяющее условию задачи и соответствующее намеченной цели, называется оптимальным планом. Нелинейное программирование служит для выбора наилучшего плана распределения ограниченных ресурсов в целях решения поставленной задачи. В общем виде постановка задачи нелинейного программирования сводится к следующему. Условия задачи представляются с помощью системы нелинейных уравнений или неравенств, выражающих ограничение, налагаемое на использование имеющихся ресурсов.
Z1(X1, X2,...,Xn) ³ 0;
Z2(X1, X2,...,Xn) ³ 0;
…
Zm(X1, X2,...,Xn) ³ 0;
при Xi ³ 0,
где Z1, Z2,…,Zm – соответствующие функции, характеризующие условие решения поставленной задачи (ограничения); Хi – искомые величины, содержащие решение задачи.
Целевая функция задается в виде:
y = f (X1, X2,…, Xn).
Причем по крайней мере одна из функций y, Z1, Z2,…, Zm – нелинейная.
Методами нелинейного программирования решаются задачи распределения неоднородных ресурсов.
Пусть имеется m разнородных ресурсов, которые предполагается реализовать для бизнеса в n регионах страны.
Известны оценочные возможности (вероятности) начать бизнес в j-м регионе (Pj), а также эффективности использования i-го ресурса в n-м регионе (wij).
Распределение ресурсов по регионам характеризуется так называемым параметром управления (hij):
hij = 0, если i-й ресурс не направляется в j-й регион,
1, если i-й ресурс направляется в j-й регион.
Необходимо распределить ресурсы по регионам таким образом (выбирать такие значения hij), чтобы величина полной вероятности достижения цели Рц была максимальной:
Рц = å Pj [1 — P (1-hijwij] = max.
Должно выполняться также ограничение
S hij = 1, i = 1, 2,…m
Ограничение означает, что каждый из m ресурсов обязательно должен назначаться в какой-либо из регионов.
Динамическое программирование (планирование)
Динамическое программирование (планирование) служит для выбора наилучшего плана выполнения многоэтапных действий. Для многоэтапных действий характерно протекание во времени. Кроме действий, естественно носящих многоэтапный характер (например, перспективное планирование), в ряде задач прибегают к искусственному расчленению на этапы, с тем, чтобы сделать возможным применение метода динамического программирования.
В общем виде постановка задачи динамического программирования сводится к следующему:
Имеется некоторая управляемая операция (целенаправленное действие), распадающаяся (естественно или искусственно) на m шагов – этапов. На каждом шаге осуществляется распределение и перераспределение участвующих в операции с целью улучшения ее результата в целом. Эти распределения в динамическом программировании называются управлениями операцией и обозначаются буквой U. Эффективность операции в целом оценивается тем же показателем, что и эффективность ее управления W (U).
При этом эффективность управления W(U) зависит от всей совокупности управлений на каждом шаге операции:
W = W(U) = W(U1, U2, ..., Um).
Управление, при котором показатель W достигает максимума, называется оптимальным управлением. Оптимальное управление обозначается буквой U.
Оптимальное управление многошаговым процессом состоит из совокупности оптимальных шаговых управлений:
U = (U1, U2, ..., Um).
Задача динамического программирования – определить оптимальное управление на каждом шаге Ui (i = 1, 2, …, m) и, тем самым, оптимальное управление всей операцией в целом.
В большинстве практических задач принимается, что показатель эффективности операции W в целом представляет собой сумму эффективности действий на всех этапах (шагах) операции:
W = å wi,
где wi – эффективность операции на i-м шаге.
При этом в случае оптимального управления
W = max åwi
Существо решения динамического программирования заключается в следующем:
- Оптимизация производится методом последовательных приближений (итераций) в два круга; в начале от последнего шага операции к первому, а затем наоборот от первого к последнему;
- На первом круге, идя от последующих шагов к предыдущим, находится так называемое условное оптимальное управление;
- Условное оптимальное управление выбирается таким, чтобы все предыдущие шаги обеспечивали максимальную эффективность последующего шага, иными словами, на каждом шаге имеется такое управление, которое обеспечивает оптимальное продолжение операции; этот принцип выбора управления называется принципом оптимальности;
- Так продолжается до первого шага, но поскольку первый шаг не имеет предыдущего, то полученное для него условное оптимальное управление теряет свой условный характер и становится просто оптимальным управлением, которое мы ищем;
- Второй круг оптимизации начинается с первого шага, для которого оптимальное управление известно.
Имея для всех шагов после него условное оптимальное управление, мы знаем, что необходимо делать на каждом последующем шаге. Это дает нам возможность последовательно переходить от условных к оптимальным управлениям дл всех последующих шагов, что обеспечивает оптимальность операции в целом.
Пусть имеется m типов различных грузов, которыми необходимо загрузить транспортное средство таким образом, чтобы общая ценность груза W была максимальной. Ценность груза является функцией отгрузоподъемности транспортного средства:
W = f (G)
Известны массы грузов i-го типа Рi и их стоимости Ci.
Необходимо загрузить транспортное средство таким образом, чтобы общая ценность груза была максимальной:
W = fm(G) = max å xiCi,
где xi – число предметов груза i-го типа, загружаемых в транспортное средство; xi выступает здесь в качестве управления (Ui=xi)
Ограничивающими условиями являются:
å xi Pi £ G
xi = 0, 1, 2…
Первое условие требует, чтобы общая масса груза не превышала грузоподъемности транспортного средства, а второе – чтобы предметы, составляющие груз различных типов, были неделимы.
Понятие критерия оптимальности
Формулировка критериев экономических систем является необходимой предпосылкой оптимизации плановых решений. В общем случае под критерием оптимальности понимается признак, на основании которого производится оценка, сравнение альтернатив, классификация объектов и явлений. Критерий оптимальности функционирования экономической системы – это один из возможных критериев (признаков) ее качества, а именно тот признак, по которому функционирование системы признается наилучшим из возможных вариантов ее функционирования. В сфере принятия экономических решений критерий оптимальности – это показатель, выражающий предельную меру экономического эффекта принимаемого хозяйственного решения для сравнительной оценки возможных решений выбора наилучшего из них. Наиболее часто используется максимум прибыли или минимум затрат.
Критерий оптимальности обычно носит количественный характер и показывает, насколько один из вариантов лучше ли хуже другого. Порядковый критерий определяет лишь то, что один вариант лучше или хуже другого. Математической формой критерия оптимальности в экономико-математических моделях является целевая функция, экстремальное значение которой характеризует предельно допустимую эффективность деятельности моделируемого объекта.
Если за классифицирующий признак принять уровень общности, то для экономической системы существуют глобальный критерий оптимального развития в масштабе Земли, социально-экономический критерий, а также «глобальный» (обобщенный) и локальный критерий оптимальности в частных системах моделей.
Если за классифицирующий признак взять математическую формулировку, то критерии подразделяются на скалярные и векторные, аддитивные и мультипликативные, интегральные критерии во временном аспекте и интегральные в пространственном аспекте и др.
Возможна классификация моделей по временному аспекту, по способам формирования критериев, по типу применяемых измерителей, по способам использования критериев.
Сущность глобального и локального критериев оптимальности.
Чаще всего термин «глобальный» применяется либо по отношению к критерию одноуровневой модели, либо по отношению к критерию «верхней» модели многоуровневой системы моделей. В последнем случае, наряду с глобальным, фигурируют локальные критерии моделей нижних уровней, отражающие интересы отдельных хозяйственных звеньев, социальных групп.
Разделение критериев на глобальный и локальный может быть отнесено к любой иерархически построенной системе моделей, например модели отрасли или предприятия.
Глобальному критерию может быть дана словесная формулировка, а для решения практических задач планирования и управления такая формулировка детализируется и представляется в виде совокупности более конкретных локальных критериев. Математически глобальный критерий принято формулировать в виде скалярной целевой функции, которая обобщенно выражает все многообразие целей или в виде векторной функции, представляющей собой набор несводимых друг к другу частных целевых функций.
Большинство многоуровневых систем имеют два уровня: верхний и нижний. Система моделей производственной программы предприятия включает в себя модели расчета общезаводских показателей и показателей отдельных цехов. При формировании обобщенных критериев должны учитываться и местные (частные интересы), а локальные критерии – подчинены обобщенному.
Сложность системы целей объясняется многообразием задач общественного развития и развития систем, а также тем, насколько обширны и интенсивны внешние связи данной системы.
Предприятие является элементом более общих систем: отрасли промышленности, эк5ономического региона. Поэтому деятельность предприятия оценивается в рамках любой из этих общих систем по соответствующим показателям. С этой точки зрения предприятие должно наилучшим образом соответствовать целям внешней системы. С другой стороны, само предприятие – сложная система, элементами которой являются коллективы его работников (бригады, отделы, службы, участки и т.д.) и отдельные индивидуумы. Следовательно, деятельность предприятия должна быть направлена на наилучшее обеспечение интересов коллектива и его работников. Система критериев оптимальности деятельности предприятия включают объемы выпуска основных типов продукции высшей категории качества, производительность труда, себестоимость продукции, фонд заработной платы.
Система критериев отраслевой системы включает удовлетворение общественных потребностей производимой продукции, экономию ресурсов, внедрение достижений научно-технического прогресса, обеспечение надежности выполнения плановых заданий. Внешние связи отраслевых систем, а значит, и комплексы их целей, усложняются фактором времени, пространственной организацией, сочетанием различных подходов и аспектов планирования.
Множественность целей развития систем существенно осложняет планирование, особенно, если цели разнонаправленные, и приближение к одним целям удаляет систему от достижения других. Таким образом возникает задача их согласования. Отыскание наилучших решений по нескольким критериям называется многокритериальной или векторной оптимизацией.
Векторная оптимизация
Математическая формулировка задачи векторной оптимизации:
Пусть X = {x1,…, x N} (j = 1,N) — вектор переменных, обычно предполагается неотрицательность вектора переменных X³0, функциональная взаимосвязь переменных устанавливается определенными соотношениями, на которые накладываются ограничения:
gi (X)£bi (i = 1,M).
Функционирование системы оценивается определенными критериями, записываемыми в виде целевых функций fr(X) (r = 1,K). Множество критериев можно представить в виде векторной целевой функции
F(X) = {f1(X),…>fr(X)}.
Чтобы минимизировать частный критерий fr(X), достаточно максимизировать -fr(X), так как min fr(X)=-max (-fr(X)). Поэтому в дальнейшем предполагается, что каждая компонента векторного критерия максимизируется. Задача многоцелевой оптимизации записывается как векторная задача математического программирования (ВЗМП)
F(X) = {f1(X),…>fr(X)} (max),
gi (X)£bi (i = 1,M),
X³0.
Будем рассматривать ВЗМП для случая, когда точки оптимума X*r(r=1,K), полученные при решении задачи по каждому критерию fr(r=1,K) не совпадают (случай их совпадения встречается крайне редко и такая задача не представляет интереса). Поэтому с математической точки зрения задача является некорректной, так как если один из критериев достигает своего оптимума, то улучшение по другим компонентам векторного критерия невозможно. Отсюда вытекает, что решением ВЗМП может быть только какое-то компромиссное решение.
Особенностью задач векторной оптимизации является наличие в области допустимых значений области компромиссов, в которой невозможно одновременное улучшение всех критериев. Принадлежащие области компромиссов планы называют эффективными, или оптимальными по Парето (по имени итальянского экономиста, впервые сформулировавшего проблему векторной оптимизации и принцип оптимальности решения).
Понятие предпочтительности плана. План X° не хуже плана X`, если
fr(X°)³ fr(X`) (r = 1,K). Если среди этих неравенств хотя бы одно строгое, то план X° предпочтительнее (лучше) X`, т.е. при переходе от X° к X`значение ни одного критерия не ухудшилось и хотя бы одного критерия улучшилось. План X° оптимален по Парето (эффективен), если он допустим и не существует другого плана X`, для которого fr(X°)³ fr(X`) (r = 1,K), и хотя бы для одного критерия выполняется строгое неравенство.
К общей формулировке многокритериальной задачи могут сводиться задачи различного содержания, которые можно подразделить на четыре типа.
1. Задачи оптимизации на множестве целей, каждая из которых должна быть учтена при выборе оптимального решения. Примером может служить задача составления плана работы предприятия, в которой критериями служит ряд экономических показателей.
2. Задачи оптимизации на множестве объектов, качество функционирования каждого из которых оценивается самостоятельным критерием. Если качество функционирования каждого объекта оценивается несколькими критериями (векторным критерием), то такая задача называется многовекторной. Примером может служить задача распределения дефицитного ресурса между несколькими предприятиями. Для каждого предприятия критерием оптимальности является степень удовлетворения его потребностей в ресурсе или другой показатель, например, величина прибыли. Для планирующего органа критерием выступает вектор локальных критериев предприятий.
3. Задачи оптимизации на множестве условий функционирования. Задан спектр условий, в которых предстоит работать объекту, и применительно к каждому условию качество функционирования оценивается некоторым частным критерием.
4. Задачи оптимизации на множестве этапов функционирования. Рассматривается функционирование объектов на некотором интервале времени, разбитом на несколько этапов. Качество управления на каждом этапе оценивается частным критерием, а на множестве этапов – общим векторным критерием. Примером может служить распределение квартального плана цеха по декадам. В каждой декаде необходимо обеспечить максимальную загрузку. В результате получится критерий максимизации загрузки в каждой декаде квартала.
Многокритериальные задачи можно также классифицировать по другим признакам: по вариантам оптимизации, по числу критериев, по типам критериев, по соотношениям между критериями, по уровню структуризации, наличию фактора неопределенности.
При разработке методов решения векторных задач приходится решать ряд специфических проблем.
Проблема нормализации возникает в связи с тем, что локальные критерии имеют, как правило, различные единицы и масштабы измерения, и это делает невозможным их непосредственное сравнение. Операция приведения критериев к единому масштабу и безразмерному виду носит название нормирования. Наиболее распространенными способами нормирования является замена абсолютных значений критериев их безразмерными относительными величинами
fr(X) = fr(X),
f*r
или относительными значениями отклонений от оптимальных значений критериев f*r
fr(X) = f*r— fr(X),
f*r
Проблема выбора принципа оптимальности связана с определением свойств оптимального решения и решением вопроса — в каком смысле оптимальное решение превосходит все остальные.
Проблема учета приоритета критериев встает, если локальные критерии имеют различную значимость. Необходимо найти математическое определение приоритета и степень его влияния на решение задачи.
Проблема вычисления оптимума возникает, если традиционные вычислительные схемы и алгоритмы непригодны для решения задач векторной оптимизации.
Решение перечисленных проблем идет в нескольких направлениях. Основные направления:
Методы, основанные на свертывании критериев в единый;
Методы, использующие ограничения на критерии;
Методы целевого программирования;
Методы, основанные на отыскании компромиссного решения;
Методы, в основе которых лежат человеко-машинные процедуры принятия решений (интерактивное программирование).
В методах, основанных на свертывании критериев, из локальных критериев формируется один. Наиболее распространенным является метода линейной комбинации частных критериев. Пусть задан вектор весовых коэффициентов критериев a = {a1,…,ar}, характеризующих важность соответствующего критерия, åar = 1, ar ³ 0 (r = 1,K). Линейная скаляризованная функция представляет собой сумму частных критериев, умноженных на весовые коэффициенты. Задача математического программирования становится однокритериальной и имеет вид
F° = åarfr(X) (max),
qi(X) £ bi (I = 1,M),
X ³ 0.
Критерии в свертке могут быть нормированы. Решение, полученное в результате оптимизации скаляризованного критерия эффективно.
К недостаткам метода можно отнести то, что малым приращениям коэффициентов соответствуют большие приращения функции, т.е. решение задачи неустойчиво, а также необходимость определения весовых коэффициентов.
Направление методов, использующих ограничения на критерии включает два подхода:
1) метод ведущего критерия;
2) методы последовательного применения критериев (метод последовательных уступок, метод ограничений).
В методе ведущего критерия все целевые функции кроме одной переводятся в разряд ограничений. Пусть g = (g2, g3,…, gк-1) – вектор, компоненты которого представляют собой нижние границы соответствующих критериев. Задача будет иметь вид
F = f1 (max)
fr ³ gr (r = 2,K),
qi (X) £ bi (I = 1,M),
X ³ 0.
Полученное этим методом решение может не быть эффективным, поэтому необходимо проверить его принадлежность области компромиссов.
Метод ведущего критерия применяется в таких задачах, как минимизация полных затрат при условии выполнения плана по производству различных видов продукции, максимизация выпуска комплектных наборов при ограничении на потребляемые ресурсы.
Алгоритм метода последовательных уступок:
1. Критерии нумеруются в порядке убывания важности.
2. Определяется значение f*1. Лицом, принимающим решение, устанавливается величина уступки D1 по этому критерию.
3. Решается задача по критерию f2 с дополнительным ограничением f1(X) ³ f*1 — D1.
Далее пункты 2 и 3 повторяются для критерия f2,…, fk.
Полученное решение не всегда принадлежит области компромиссов.
При решении задач методами целевого программирования предполагается приближение значения каждого критерия к определенной величине fr, т.е. достижение определенной цели. В самом общем виде задача целевого программирования формулируется как задача минимизации сумм отклонений целевых функций от целевых значений с нормированными весами.
d(F(X), F) = ( å wR êfR (X) — fR êp) (min),
где F = {f1,...., fR} — вектор целевых значений,
W = {w1,..., wR} — вектор весов, обычно å wR= 1, wR³ 0
(r = 1, K), значения p находятся в пределах 1 £ p £ ¥,
d(.) – расстояние (мера отклонения) между F(X) и F.
Во многих случаях применения целевого программирования полагают p = 1. Например, в линейном целевом программировании функции fR (X) (r=1, K) и