MO_conspect

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

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

71

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

Правила выбора разрешающего элемента при поиске опорного плана:

1.При условии отсутствия “0-строк” (ограничений-равенств) и “свободных” переменных (т.е. переменных, на которые не наложено требование неотрицательности).

Если в столбце свободных членов симплексной таблицы нет отрицательных элементов, то опорный план найден.

Есть отрицательные элементы в столбце свободных членов, например bi 0 . В такой строке ищем отрицательный коэффициент

ail , и этим самым определяем разрешающий столбец l . Если не найдем отрицательный ail , то система ограничений несовмест-

на (противоречива).

В качестве разрешающей выбираем строку, которой соответствует минимальное отношение:

где r – номер разрешающей строки. Таким образом, arl – разрешающий элемент.

72

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

2.В случае присутствия ограничений-равенств и “свободных” переменных поступают следующим образом.

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

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

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

Рассматривая симплекс-метод, мы предполагали, что задача линейного программирования является невырожденной, т.е. каждый опорный план содержит ровно m положительных компонент, где m – число ограничений в задаче. В вырожденном опорном плане число положительных компонент оказывается меньше числа ограничений: некоторые базисные переменные, соответствующие данному опорному плану, принимают нулевые значения.

Используя геометрическую интерпретацию для простейшего случая, когда n m 2 (число небазисных переменных равно 2), легко отличить вырожденную задачу от невырожденной. В вырожденной задаче в одной вершине многогранника условий пересекается более двух прямых, описываемых уравнениями вида xi 0 . Это значит, что одна или несколько сто-

рон многоугольника условий стягиваются в точку.

Аналогично при n m 3 в вырожденной задаче в одной вершине пересекается более 3-х плос-

костей xi 0 .

В предположении о невырожденности задачи находилось только одно значение, соответствующее минимальному значению

73

зиса вектора условий (выводимой из числа базисных переменной).

скольких индексах сразу (для нескольких строк). В этом случае в находимом опорном плане несколько базисных переменных будут нулевыми.

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

Один из приемов борьбы с вырожденностью состоит в преобразовании задачи путем “незначительного” изменения вектора правых частей системы ограничений на величины i , таким образом, чтобы задача стала

невырожденной и, в то же время, чтобы это изменение не повлияло реально на оптимальный план задачи.

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

Пусть переменную x j необходимо сделать базисной. Рассмотрим множество индексов E0 , состоящее из тех i , для которых достигается

то из базиса исключается вектор условий Ai (переменная xi делается небазисной).

дится переменная xk . В противном случае составляется множество E2 и

т.д.

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

74

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

8.4.1. Введение

Рассмотрим задачу максимизации линейной формы (1) и, одновременно, задачу минимизации (2):

Задача (2) называется двойственной по отношению к прямой (1) (и наобо-

рот!) [8].

Пример: Предприятие выпускает три вида продукции. Каждая продукция требует обработки на трех различных типах установок. Ресурс времени каждого типа установок ограничен. Известна прибыль от единицы каждого вида продукции p1 , p2 , p3 . Если количество выпускаемой про-

дукции каждого вида x1 , x2 , x3 , тогда необходимо максимизировать прибыль

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

a11x1 a12 x2 a13 x3 b1 , a21x1 a22 x2 a23 x3 b2 , a31x1 a32 x2 a33 x3 b3 ,

x 0 ,

где b1, b2 , b3 – ресурсы времени установок первого, второго и третьего типов. Величины aij определяют количество ресурса времени установки i -го типа, которое необходимо для выпуска одной единицы продукции j -го

вида.

Двойственная к ней задача будет иметь вид

W u b1u1 b2u2 b3u3 min

при ограничениях:

a11u1 a21u2 a31u3 p1 ,a12u1 a22u2 a32u3 p2 ,a13u1 a23u2 a33u3 p3 , .

75

Здесь u1 - это оценка (цена), соответствующая одной единице ограничен-

ного ресурса по первой установке. И она равна величине, на которую могла бы увеличиться суммарная прибыль, если бы количество этого ограниченного ресурса увеличилось на единицу, и если это увеличение было бы использовано оптимально. Иными словами, u1 – это количество прибыли,

недополученной из-за нехватки единицы ограниченного ресурса b1 . Аналогичным образом можно интерпретировать смысл величин u2 и u3 .

8.4.2.Преобразования при решении прямой и двойственной задач

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

Для ограничений прямой задачи симплексная таблица имеет вид:

Пусть arS - разрешающий элемент. Сделаем шаг модифицированного жорданова исключения и получим таблицу:

76

где bij aij arS arj aiS , и всю данную таблицу следует разделить еще на

arS .

Симплексную таблицу для двойственной задачи запишем, развернув ее на 90 . Получаем:

Пусть arS - направляющий элемент. Сделаем шаг обыкновенного

жорданова исключения (отличие от модифицированного состоит в том, что элементы в разрешающей строке меняют знаки, а в столбце знаки сохраняются; в остальном преобразование остается тем же):

где bij aij arS arj aiS , и всю данную таблицу также следует разделить еще

на arS .

Замечание: Не следует забывать при преобразованиях, что в данном случае у нас таблица развернута.

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

77

Можно показать, что, решая основную задачу линейного программирования, решаем и двойственную к ней. И наоборот. Причем, max Q min W .

8.4.3. Теоремы двойственности

8.4.3.1. Основная теорема двойственности линейного программирования

Теорема [8]. Пусть рассматривается пара двойственных задач:

Если одна из этих задач обладает оптимальным решением, то и двойственная к ней задача также имеет оптимальное решение. Причем экстремальные значения соответствующих линейных форм равны: max Q min W .

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

Доказательство: Пусть основная задача (1) имеет конечное решение и получена окончательная симплексная таблица:

78

Так как данная таблица, по предположению, соответствует опти-

y1 ... ys xs 1 ... xn 0.

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

v1 ... vs us 1 ... um 0 ,

найдем u1 q1 0, …, us qs 0, vs 1 qs 1 0, …, vn qn 0 . Следовательно, получено опорное решение:

u1 q1 , …, us qs , us 1 0 , …, um 0.

Из последнего столбца,

W b1,n 1v1 ... bs,n 1vs bs 1,n 1us 1 ... bm,n 1um q0

в точке

v1 ... vs us 1 ... um 0

будет минимальным в силу того, что bi,n 1 0 i , i 1, m . Следовательно, max Q min W .

Пусть теперь линейная форма прямой задачи неограничена, т.е. для некоторой верхней переменной, например, ys соответствующий коэффициент qs 0, а все коэффициенты этого столбца симплексной таблицы неположительны: b1,s 0 , b2,s 0 , …, bm,s 0 . Тогда из таблицы для

двойственной задачи:

us b1,sv1 ... bs,svs bs 1,sus 1 ... bm,sum qs qs 0 ,

то есть система ограничений двойственной задачи противоречива. Так как из неотрицательности v1, ..., vs , us 1, ..., um следует неположительность us (нельзя сделать ее положительной). То есть, система несовместна.

Теорема доказана.

8.4.3.2. Вторая теорема двойственности

Теорема [8]. Если хотя бы одно оптимальное решение одной из двойственных задач обращает i -е ограничение этой задачи в строгое неравенство, то i -я компонента (т.е. xi или ui ) каждого оптимального решения вто-

рой двойственной задачи равна нулю.

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

Т.е. оптимальные решения x* и u * пары двойственных задач удовлетворяют условиям

79

80