Существование решения озлп в способы его нахождения

Рассмотрим основную задачу линейного программирования (ОЗЛП): найти неотрицательные значения переменных х₁, x₂,..., х_n, удовлетворяющие т условиям-равенствам:

(9.1)

 

и обращающие в максимум линейную функцию этих переменных:

 

L = с₁х₁ + с₂х₂ + … + с_nx_n => max. (9.2)

 

Для простоты предположим, что все условия (9.1) линейно независимы (r = m), и будем вести рассуждения в этом предположении ¹).

Назовем допустимым решением ОЗЛП всякую совокупность неотрицательных значений х₁, х₂,..., x_n, удовлетворяющую условиям (9.1). Оптимальным назовем то из допустимых решений, которое обращает в максимум функцию (9.2). Требуется найти оптимальное решение.

Всегда ли эта задача имеет решение? Нет, не всегда. Во-первых, может оказаться, что уравнения (9.1) вообще несовместны (противоречат друг другу).

 

 

¹) Если среди условий-равенств (9.1) есть некоторые лишние, вытекающие из других, то это обнаружится автоматически в процессе решения ОЗЛП (см. [4, 6]).

 

 

Может оказаться и так, что они совместны, но не в области неотрицательных решений, т. е. не существует ни одной совокупности чисел х₁ 0, х₂ 0,..., х_n 0, удовлетворяющей условиям (9.1). Наконец, может быть и так, что допустимые решения ОЗЛП существуют, но среди них нет оптимального: функция L в области допустимых решений не ограничена сверху. Все эти опасности подстерегают нас, главным образом, в «придуманных», искусственно поставленных задачах, хотя иногда легкомысленное планирование (неполный учет имеющихся ресурсов) приводит к неразрешимым задачам линейного программирования.

Чтобы представить себе принципиальную сторону ОЗЛП, обратимся к геометрической интерпретации. Пусть число уравнении т па два меньше числа переменных п (п - т = k = 2). Такой частный случай дает возможность геометрической интерпретации ОЗЛП па плоскости.

Мы знаем, что т линейно независимых уравнении (9.1) всегда можно разрешить относительно каких-то т базисных переменных, выразив их через остальные, свободные, число которых равно п - т = k (в нашем случае k = 2). Предположим, что свободные переменные — это х₁ и х₂ (если это не так, то всегда можно заново перенумеровать переменные), а остальные: х₃, х₄,..., х_n— базисные. Тогда вместо т уравнений (9.1) мы получим тоже т уравнений, но записанных в другой форме, разрешенных относительно х₃, х₄,...:

(9.3)

 

Будем изображать пару значений свободных переменных точкой с координатами х₁, х₂ (рис. 9.1). Так как переменные х₁, х₂ должны быть неотрицательными, то допустимые значения свободных переменных лежат только выше оси Ox₁ (на которой х₂ = 0) и правее оси Ох₂ (на которой х₁ = 0). Это мы отметим штриховкой, обозначающей «допустимую» сторону каждой оси.

Теперь построим па плоскости х₁0х₂ область допустимых решений или же убедимся, что ее не существует. Базисные переменные х₃,x₄,..., x_n тоже должны быть неотрицательными и удовлетворять уравнениям (9.3). Каждое такое уравнение ограничивает область допустимых решений. Действительно, положим

 

 

в первом уравнении (9.3) х₃ = 0; получим уравнение прямой линии:

 

 

На этой прямой х₃ =0; по одну сторону от нее х₃ >0, по другую — х₃ < 0. Отметим штриховкой ту сторону (полуплоскость), где х₃ > 0 (рис. 9.2). Пусть эта сторона оказалась правее и выше прямой х₃ = 0. Значит, вся область допустимых решений (ОДР) лежит в первом координатном угле, правее и выше прямой х₃ = 0. Аналогично поступим и со всеми остальными условиями (9.3). Каждое из них изобразится прямой со штриховкой, указывающей «допустил1ую» полуплоскость, где только и может лежать решение (рис. 9.3).

Таким образом, мы построили п прямых: две оси координат {Ox₁ и Ox₂) и п — 2 прямых х₃ = 0, х₄ = 0,..., х_n = 0. Каждая из них определяет «допустимую» полуплоскость, где может лежать решение. Часть первого координатного угла, принадлежащая одновременно всем этим полуплоскостям, и есть ОДР. На рис. 9.3 показан случай, когда ОДР существует, т. е. система уравнений (9.3) имеет неотрицательные решения, Заметим, что этих решении— бесконечное множество, так как любая пара значений свободных переменных, взятая из ОДР, «годится», а из х₁ и х₂ могут быть определены и базисные переменные.

Может оказаться, что область допустимых решений не существует, и значит, уравнения (9.3) несовместны в области неотрицательных значений. Такой случай показан на рис. 9.4, где нет области, лежащей одновременно по «нужную» сторону от всех прямых. Значит, ОЗЛП не имеет решения.

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

Для этого дадим геометрическую интерпретацию условию (9.2)

L => max. Подставив выражения (9.3) в формулу (9.2), выразим L через свободные переменные х₁, x₂. После приведения подобных членов получим:

 

L= , (9.4)

 

где , — какие-то коэффициенты, — свободный член, которого в первоначальном виде у функции L не было; теперь, при переходе к переменным х₁, х₂, он мог и появиться. Однако мы его тут же и отбросим: ведь максимум линейной функции L достигается при тех же значениях х₁, х₂, что и максимум линейной однородной функции (без свободного члена):

 

L^’ = . (9.5)

 

Посмотрим, как изобразить геометрически условие L^’ => шах. Положим сначала L^’ = 0, т.е. и построим на - плоскости х₁Ох₂ прямую с таким уравнением; очевидно, он проходит через начало координат (рис. 9.5). Назовем ее «опорной прямой». Если мы будем придавать L^’ какие-то значения С₁, С₂, С₃,..., прямая будет перемещаться параллельно самой себе; при перемещении в одну сторону L^’ будет возрастать, в другую — убывать. Отметим на рис. 9.5 стрелками, поставленными у опорной прямой, то направление, в котором L' возрастает. На рис. 9.5 это оказалось направление «направо — вверх», но могло быть и наоборот: все зависит от коэффициентов . Теперь изобразим опорную прямую и ОДР на одном чертеже (рис. 9.6). Давайте будем мысленно двигать опорную прямую параллельно самой себе в направлении стрелок (возрастания L^’) Когда L^’ достигнет максимума? Очевидно, в точке А (крайней точке ОДР в направлении стрелок). В этой точке свободные переменные принимают оптимальные значения a из них можно по формулам (9.3) найти и оптимальные значения всех остальных (базисных) переменных ¹). Заметим, что максимум L^’ достигается в одной из вершин ОДР, где, по крайней мере, две из базисных переменных (в нашем случае это х₃ и х₅) обращаются в нуль. Могло бы обращаться в нуль и больше базисных переменных, если бы через точку А проходило более двух прямых х_i = 0.

 

 

А может ли оказаться, что оптимального решения не существует? Да, может, если в ОДР функция L^’ (а значит и L) не ограничена сверху. Пример такого ненормального случая показан на рис. 9.7 (в разумно поставленных задачах обычно такого недоразумения не возникает).

На рис. 9.6 оптимальное решение существовало и было единственным. А сейчас рассмотрим случай, когда оптимальное решение существует, но не единственно (их бесконечное множество).

 

¹) Если бы стрелки были направлены иначе («налево — вниз»), то крайней точкой ОДР в направлении стрелок была бы точка В.

 

 

Это случай, когда максимум L^’ достигается не в одной точке А, а на целом отрезке АВ, параллельном опорной прямой (рис. 9.8). Такой случай встречается на практике, но нас он

 

 

не должен волновать. Все равно и в этом случае максимум U достигается в какой-то из вершин ОДР (4 или В — безразлично), и в поисках оптимального решения можно ограничиться вершинами ОДР.

Итак, мы рассмотрели в геометрической интерпретации случай п - т = k=2 и убедились в следующем: оптимальное решение (если оно существует) всегда достигается в одной из вершин ОДР, в точке, где по крайней мере две из переменных х₁, х₂,..., х_n равны нулю.

Оказывается, аналогичное правило справедливо и в случае n - m = k > 2 (только геометрическая интерпретация теряет в этом случае свою наглядность). Обойдемся без доказательства, просто сформулируем это правило.

Оптимальное решение ОЗЛП (если оно существует) достигается при такой совокупности значений переменных х₁, х₂,..., х_n, где по крайней мере k из них обращаются в нуль, а остальные неотрицательны.

При k = 2 такая совокупность значений изображается точкой на плоскости, лежащей в одной из вершин многоугольника допустимых решений (ОДР). При k = 3 ОДР представляет собой уже не многоугольник, а многогранник, и оптимальное решение достигается в одной из его вершин. При k > 3 геометрическая интерпретация теряет наглядность, но все же геометрическая терминология остается удобной. Мы будем продолжать говорить о «многограннике допустимых решений» в k -мерном пространстве, а оптимальное решение (если оно существует) будет достигаться в одной из вершин этого многогранника, где, по крайней мере, k переменных равны нулю, а остальные — неотрицательны. Будем для краткости называть такую вершину «опорной точкой», а вытекающее из нее решение — «опорным решением».

Отсюда вытекает идея, лежащая в основе большинства рабочих методов решения ОЗЛП,— идея «последовательных проб». Действительно, попробуем разрешить уравнения (9.1) относительно каких-нибудь m базисных переменных и выразим их через остальные k свободных. Попробуем положить эти свободные переменные равными нулю — авось повезет, наткнемся на опорную точку. Вычислим базисные переменные при нулевых значениях свободных. Если все они оказались неотрицательными, значит, нам повезло, мы сразу же получим допустимое (опорное) решение, и его остается только оптимизировать. А если нет? Значит, данный выбор свободных и базисных переменных допустимого решения не дает; точка лежит не на границе, а вне ОДР. Что делать? Надо «переразрешить» уравнения относительно каких-то других базисных переменных, но не как попало, а так, чтобы это приближало нас к области допустимых решений (для этого в линейном программировании существуют специальные приемы, на которых мы останавливаться не будем). Пусть, наконец, несколько раз повторив такую процедуру, мы нашли опорное решение ОЗЛП. Браво! Но это еще не все. Тут надо поставить вопрос: а является ли это решение оптимальным? Выразим функцию L через последние получившиеся свободные переменные и попробуем увеличивать их сверх нуля. Если от этого значение L только уменьшается, значит, нам повезло, и мы нашли оптимальное решение, ОЗЛП решена. А если нет? Снова «пере разрешаем» систему уравнений относительно других базисных переменных, и снова не как попало, а так, чтобы, не выходя за пределы допустимых решений, приблизиться к оптимальному. И опять-таки для этого в линейном программировании существуют специальные приемы, гарантирующие, что при каждом новом «пере разрешении» мы будем приближаться к оптимальному решению, а не удаляться от него. На этих приемах мы тоже здесь не будем останавливаться. После конечного числа таких шагов цель будет достигнута — оптимальное решение найдено. А если его не существует? Алгоритм решения ОЗЛП сам покажет вам, что решения нет.

Читатель, может быть, спросит: только-то и всего? Зачем было придумывать хитрые методы решения ОЗЛП, может быть, надо просто перебрать, одну за другой, все возможные комбинации k свободных переменных, полагая их равными нулю, пока, наконец, не будет найдено оптимальное решение?

Действительно, для простых задач, где число переменных невелико, такой «слепой перебор» может привести к решению, и довольно быстро. Но на практике часто встречаются задачи, в которых число переменных (и наложенных условий) очень велико, порядка сотен и даже тысяч. Для таких задач простой перебор становится практически невозможным: слишком велико число комбинаций свободных и базисных переменных. Пример: только при п = 30 и m = 10 число возможных комбинаций свободных переменных с базисными равно = 30045015, т. е. свыше 30 миллионов! А эта задача — далеко не из сложных.

Разработанные в теории линейного программирования вычислительные методы («симплекс-метод», «двойственный симплекс-метод» и другие, см. [4,6,8]) позволяют находить оптимальное решение не «слепым» перебором, а «целенаправленным», с постоянным приближением к решению. Добавим, что современные ЭВМ, как правило, снабжены подпрограммами для решения задач линейного программирования, так что лицу, желающему их решать, нет даже особой надобности обучаться решению таких задач «вручную» — труд крайне неприятный и изнурительный.