Графический (геометрический) метод решения задач ЛП.

Графические методы связаны прежде всего с геометрическим изображением функциональной зависимости при помощи линий на плоскости. Графики используются для быстрого нахождения значения функций по соответствующему значению аргумента, для наглядного изображения функциональных зависимостей.
В экономическом анализе применяются почти все виды графиков: диаграммы сравнения, диаграммы временных рядов, кривые распределения, графики корреляционного поля, статистические картограммы. Особенно широко распространены в анализе диаграммы сравнения - для сравнения отчетных показателей с плановыми, предшествующих периодов и передовых предприятий отечественных или зарубежных. Для наглядного изображения динамики экономических явлений (а в анализе с динамическими рядами приходится иметь дело очень часто) используются диаграммы временных рядов.
С помощью координатной сетки строятся графики зависимости, например, уровня издержек от объема произведенной и реализованной продукции, а также. графики, на которых можно изображать и корреляционные связи между показателями. В системе осей координат изображение показывает влияние различных факторов на тот или иной показатель.
Широко применяется графический метод для исследования производственных процессов, организационных структур, процессов программирования и т. д. Например, для анализа эффективности использования производственного оборудования строятся расчетные графики, в том числе графики множественных факторов.

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

В математически формализованной системе анализа, планирования и управления особое место занимают сетевые графики. Они дают большой экономический эффект при строительстве и монтаже промышленных и других предприятий.
Сетевой график (рис. 6.1) позволяет выделить из всего комплекса работ наиболее важные, лежащие на критическом пути, и сосредоточить на них основные ресурсы строительномонтажных организаций, устанавливать взаимосвязь между различными специализированными организациями и координировать их работу. Работы, лежащие на критическом пути, требуют наиболее продолжительного ожидания поступления очередного события. На стадии оперативного анализа и управления сетевой график дает возможность осуществлять действенный контроль за ходом строительства, своевременно принимать меры по устранению возможных задержек в работе.
Применение сетевых графиков анализа, планирования и управления обеспечивает, как показывают многие примеры, сокращение сроков строительства на 20-30%, повышение производительности труда на 15-20%.
При анализе, осуществляемом непосредственно на стройках, использование материалов сетевого планирования и управления способствует правильному определению причин, влияющих на ход строительства, и выявлению предприятий, не обеспечивающих выполнение порученных им работ или поставку оборудования в сроки, установленные графиком.
Разработка сетевого графика в строительстве осуществляется при наличии: норм продолжительности строительства и срока ввода в действие объекта или комплекса объектов, проектно-сметной документации, проекта организации строительства и производства работ, типовых технологических карт, действующих норм затрат труда, материалов и работы машин. Кроме того, при составлении графика используются опыт выполнения отдельных работ, а также данные о производственной базе строительных и монтажных организаций.
На основе всех этих данных составляется таблица работ и ресурсов, где в технологической последовательности производства работ указываются их характеристика, объем, трудоемкость в человеко-днях, исполнитель (организация и бригада), численность рабочих, сменность, потребность в механизмах и материалах, источники их поступления, общая продолжительность выполнения работы в днях, а также предшествующее задание, после окончания которого можно начинать данную работу. Исходя из показателей такой таблицы, подготавливают сетевой график, который может иметь различную степень детализации в зависимости от принятой схемы произ
водства работ и уровня руководства; кроме общего графика исполнители разрабатывают график выполняемых ими работ.
Основные элементы сетевого графика: событие, работа, ожидание, зависимость.
При анализе хода строительства объекта следует устанавливать, правильно ли составлен сетевой график, не допущено ли при этом завышение критического пути, учтены ли при оптимизации графика все возможности его сокращения, нельзя ли какие-либо работы выполнять параллельно или сократить время, затрачиваемое на них, путем увеличения средств механизации и др. Это особенно важно в тех случаях, когда продолжительность работ по графику не обеспечивает окончание строительства в срок.
Основным материалом сетевого планирования, используемого при анализе, является информация о ходе работ по графику, который обычно составляется не реже одного раза в декаду. В качестве примера приводится карта задания и информации о ходе работы по объекту строительства, осуществляемому по сетевому графику (табл. 6.1). По данным карты, критические работы выполнялись в начале месяца с опережением графика, однако затем было допущено отставание монтажа подкрановых балок по ряду Б, а последующая работа - монтаж подкрановых балок по ряду А - закончена с отставанием на один день.
Оптимизация сетевых графиков осуществляется на стадии планирования посредством сокращения критического пути, т. е. минимизации сроков выполнения строительных работ при заданных уровнях ресурсов, минимизации уровня потребления материальных, трудовых и финансовых ресурсов при фиксированных сроках выполнения строительных работ. Возможен и смешанный подход: для одной части работ (более дорогостоящих) - минимизировать уровень потребления ресурсов при фиксированных сроках выполнения работ, для другой - минимизировать сроки при фиксированном уровне ресурсов.
Решение оптимизационных задач существенно облегчается наличием пакетов прикладных программ (ППП), приспособленных к составлению оптимальных сетевых графиков на ЭВМ.
В зарубежной практике системного анализа распространен графо-математический метод, получивший название «дерево решений». Суть этого метода заключается в следующем.
Путем предварительной оценки потребностей, предварительного анализа возможных организационных, технических или технологических условий намечаются все предполагаемые варианты решения данной задачи. Вначале разрабатываются



Задание


Информация

Резерв времени по работам

Чис
тый

Наименование
работ

шифр

дата
начала

дата
оконча

плановая
продол

Ре
зерв
вре

%
тех-

требуемое время для

при
чина

фактическая дата

находя
щимся

не находящимся

резерв времени с


работ

работ
(план)

ния
работ
(план)

житель
ность,
дней

мени

кой
готов
ности

оконча
ния
работ,
дней

задер
жки

оконча
ния
работ

на критическом пути

аа критическом пути

начала месяца, дней

1

2

3

4

5

6

7

8

9

10

11

12

13

Разработка грунта

1-2

1/IV

6/IV

5

0

100

-

-

6/IV

¦-

-

-

Бетонирование фундаментов под котлы

2-3

7/IV

17/1V

9

0

100

14/IV

2

2

Бетонирование фундаментов по ряду А

2-4

7/IV

14/1V

7

2

100

14/IV




То же по ряду Б

2-5

7/IV

14/IV

7

2

100

-

-

14/IV




Устройство трубной разводки

6-18

18/IV

21/IV

4

19

100

-

-

29/IV

-7

Устройство обратной засыпки

6-7

18/IV

19/IV

2

0

100

17/IV

2

2

Монтаж сборных железобетонных ко













лонн:
по ряду Б

7-8

20/IV

22/IV

3

1

100

-

-

22/IV

_

-

-

по ряду А

7-9

20/IV

22/IV

3

1

100

-

-

22/IV

-

-

-

Устройство подкрановых путей и монтаж башенного крана 7-10
Установка опорных рам на фундамент под оборудование 7-16 Монтаж подкрановых балок:
по ряду Б 8-11
20/IV 24/IV 4
20/IV 24/IV 4
24/IV 25/IV 2

по ряду А 10-12 25/IV 26/IV
Монтаж первой части балок и плит покрытия 12-13 27/IV 4/V
Монтаж подкрановых путей мостового lt;3 крана 12-14 27/IV 3/V


6

7

8

9

10

11

12

13

0

100

-

-

22/IV

1

-

1

14

100.

-

-

29/IV

-

-5

-

1

100

за-

27/IV

-2

27/IV -1
держ- ка с поставкой ж/б конструкций
  1. 100 -

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


f = –х 1 + 5х 2 ¾> min ;

4х 1+ 3х 2 £ 24,

х 1– 10х 2 £ 0,

8х 1– 3х 2 ³ 0,

5х 1+ 3х 2 ³ 15,

х 1³0, х 2³ 0. (1)

Совокупность переменных хj , удовлетворяющих условию (1), называется областью допустимых решений. Допустимое решение, обращающее целевую функцию в min или max , называется оптимальным. Для его определения необходимо построить область допустимых решений (область определения). Так как в условии задачи заданы две переменные, то область допустимых решений находится на плоскости х 10х 2. Каждое неравенство (1) определяет полуплоскость, а равенство – прямую. Для построения полуплоскости необходимо найти ее границу и установить, с какой стороны от нее лежит искомая полуплоскость. Перепишем условия (1) в виде равенств (2) и пронумеруем их.

4х 1+ 3х 2 = 24 (I ),
х 1– 10х 2 = 0 (II ),
8х 1– 3х 2 = 0 (III ),
5х 1+ 3х 2 = 15 (IV ). (2)

Введем систему координат х 10х 2 и построим последовательно эти прямые – границы полуплоскостей. Для построения прямой на плоскости необходимо определить любые две точки, лежащие на этой прямой. Если прямая пересекает оси 0х 1и 0х 2, то можно найти координаты точек ее пересечения с осями координат. Определим координаты пересечения прямой (I ) с осью 0х 1: х 1=0; Þ 3х 2= 24; Þ х 2= 8. Соответственно определим координаты второй точки пересечения первой прямой с осью 0х 2: х 2=0; Þ 4х 1= 24; Þ х 1= 6. Следовательно, точки пересечения прямой (I ) с осями координат равны (0,8) и (6,0). Построим эту прямую (рис. 1).

Определим полуплоскость. Для этого подставим в первое неравенство (1) координаты любой точки, не лежащей на данной прямой, например (0,0). Тогда из первого условия следует: 4×0+3×0 £24, значит, неравенство справедливо, откуда следует, что полуплоскость лежит с той стороны прямой, где находится точка с координатами (0,0).


Аналогичным образом строятся и другие полуплоскости. Необходимо учесть, что прямые (II) и (III) проходят через начало координат, т.е. точку (0,0). Координаты второй точки желательно брать пропорционально коэффициентам в уравнении искомой прямой. Например, для второй прямой – точки (0,0) и (10,1), а для третьей – (0,0) и (3,8). После построения всех полуплоскостей область допустимых решений примет следующий вид (рис. 3):



Целевая функция f определяет на плоскости прямую, которая должна проходить через точку или сторону многоугольника и иметь наименьшее значение. Построим направляющий вектор для этой прямой. Данный вектор перпендикулярен искомой прямой, и его направление всегда определяет максимум целевой функции. Противоположное направление вектора определяет минимум. Обозначим этот вектор через . Он проходит через точку (0,0) и (–1,5). Координаты второй точки берут из коэффициентов целевой функции и с их помощью определяют направление вектора. Перпендикулярно ему построим прямую –х 1+ 5х 2=0. Как было сказано выше, вектор всегда показывает направление возрастания значения целевой функции (max ) , противоположный ему вектор –– направление убывания значения целевой функции (min ). Перемещаем прямую –х 1+5х 2=0 по области определения параллельно самой себе в направлении min . Целевая функция f достигнет своего минимального значения в точке С (рис. 4).


Оптимальному решению задачи (1) соответствует точка С , которая лежит на пересечении прямых (I ) и (II ):

4х 1+ 3х 2= 24;

х 1– 10х 2= 0.

Для решения данной системы уравнений умножить второе уравнение на 4 и сложить соответственно по элементам с 1-м уравнением:

4х 1+ 3х 2 = 24;

4х 1– 40х 2 = 0.

Вычтем из первого уравнения второе, получим: 43х2= 24 Þ х 2= 0,56.

Подставив найденное значение х 2во второе уравнение, получим:

х 1= 10х х 1=5,6. Подставив координаты точки С в целевую функцию, получим следующий результат:

f min = – 5,6 + 5×0,56 = – 2,8.

Окончательный результат задачи запишем в следующем виде:

х 1= 5,6, х 2= 0,56;f min = – 2,8.

Решение данного примера на ПЭВМ осуществляется программным комплексом «Блок-3». С его помощью производятся ввод, решение и вывод результативной информации на внешний носитель. Простота и доступность комплекса позволит без труда освоить его и применять на практике.

Задача № 1.1.2.

f = 2х 1+ 3х 2 ¾> max;

2х 1+ 3х 2 £ 12,

2х 1– 5х 2 £ 0,

7х 1– 2х2³ 0,

х 1, х 2³ 0. (3)

Определения и построение области допустимых решений аналогичны заданию 1.1.1. Окончательный вид области допустимых решений представлен на рис. 5 многоугольником АВС (точка А совпадает с точкой 0).

Очевидно, что прямая, определяющая целевую функцию, совпадает с прямой, образующей сторону многоугольника ВС . Отсюда следует, что решением данной ЭММ являются точки, лежащие на стороне ВС много-

угольника АВС . Для записи решения ЭММ необходимо найти координату x 1B – точки В и x 1C – точки С . Определив их, мы сможем найти отрезок, лежащий на оси 0x 1(рис. 6).


Координаты точки В – x1B определяются в результате пересечения прямых 2х 1+ 3х 2 = 12 и 7х 1– 2х 2 = 0. Для этого необходимо решить систему уравнений:

2х 1+ 3х 2= 12 ´ 2 Þ 4х 1+ 6х 2= 24;

7х 1– 2х 2= 0 ´ 3 Þ 21х 1– 6х2= 0.

Сложив два последних уравнения, получим: 25х 1=24, х 1=0,96. Из этого следует, что x 1B =0,96. Координата точки С x 1C определяется в результате пересечения прямых 2х 1+ 3х 2=12 и 2х 1–5х 2=0. Решим систему уравнений:

2х 1+ 3х 2= 12 ´ 5 Þ 10х 1+ 15х 2= 60;

2х 1– 5х 2= 0 ´ 3 Þ 6х 1 – 15х 2= 0.

Сложив два последних уравнения, получим: 16х 1= 60, х 1= 3,75, откуда следует, что x 1C = 3,75.

Значение целевой функции для данной ЭММ равно 12 (так как уравнение прямой, на которой определен отрезок ВС – 2х 1+3х 2= 12).

Таким образом, ответ данной задачи:

x 1Î[x 1B ; x 1C ] Þ x 1Î;

2х 1+ 3х 2=12 Þ 3х 2= 12 – 2х х 2= (12 – 2х 1)/3.

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

x 1Î; x 2= (12 – 2х 1)/3; f max = 12.

Задача № 1.1.3.

f = 2х 1+ 3х 2 ¾> max;

2х 1+ 3х 2 ³ 12,

2х 1– 5х 2 £ 0,

7х 1– 2х 2³ 0,

х 1, х 2 ³0. (4)

Используя схему построения области допустимых решений задач 1.1.1–1.1.2, получим следующий график (рис. 7):


f = 2х 1+ 3х 2 ¾> max ;

х 1+ х2 £ 2,

2х 1+ 3х 2³ 12,

2х 1– 5х 2£ 0,

7х 1– 2х 2³ 0,

х 1, х 2³ 0. (5)

Используя график задачи 1.1.3 и достроив первую полуплоскость х 1+х2£ 2, получим область определения, показанную на рис. 8.


Из графика (рис. 8) видно, что для данной ЭММ области допустимых решений нет. Ответ: нет области допустимых решений.

Задача № 1.1.5.

f = – х 1+ 5х 2 ¾> min;

10х 1+ 3х 2£ 30,

10х 1+ 5х 2³ 50,

2х 1– 6х 2£ 0,

х 1, х 2³ 0. (6)

Область определения ЭММ (6) представлена на рис. 9. Из анализа графика следует, что областью допустимых решений будет являться точка А с координатами (0,10) (10х 1+ 5х 2= 50, х 1= 0, 5х 2= 50, х 2=10). В случае, когда решением ЭММ является единственная точка, целевую функцию можно не строить.

Ответ: x 1= 0; x 2=10; fmin = 0+5×10 = 50.


Таким образом, при решении задач ЭММ ЛП возможны следующие ситуации:

– задача имеет одно оптимальное решение;

– задача имеет бесконечное число оптимальных решений;

– задача не имеет оптимального решения;

– задача не имеет области допустимых решений.

На практике ЭММ ЛП не имеет решений только в том случае, если некорректна постановка задачи.

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

Задача № 1.1.6.

Предприятие может организовать производство своей продукции двумя способами. При первом способе предприятие за месяц выпускает C 1 тыс. изделий, при втором – C 2 тыс. изделий. Расход производственных, людских ресурсов, амортизация оборудования и ограничения ресурсов, приведены ниже в таблице.

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

1) Решить графическим способом;

2) Решить на базе комплекса «Блок-3»;

3) Симплекс-методом.

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

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

Алгоритм графического метода решения злп

Реализацию графического метода решения ЗЛП рассмотрим на примерах.

Пример 2.2.1. Решить ЗЛП графическим методом:

(2.2.1)

max z =x 1 + 4x 2 (2.2.2)

Решение. Для построения области допустимых решений, которая состоит из пересечения полуплоскостей, соответствующих каждому неравенству системы ограничений (2.2.1), запишем уравнения граничных прямых:

l 1: x 1 + 5x 2 = 5; l 2: x 1 + x 2 = 6; l 3: 7x 1 + x 2 = 7.

l 1 к виду (2.2.3.) разделим обе его части на 5:
. Таким образом, прямаяl 1 отсекает на оси Ох 1 5 единиц, на оси Ох 2 1 единицу. Аналогично имеем для l 2:
иl 3:
.

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

Таким образом, первая и вторая искомые полуплоскости расположены в противоположную сторону от начала координат (0 – 5·0– 5; 7·0 + 07), а вторая – в сторону начала координат (0 + 06). Область допустимых решений на рисунке 2.2.1 заштрихована.

Рисунок 2.2.1 – Область допустимых решений

Для нахождения оптимального плана, который будет находиться в вершине многоугольника решений, нужно построить вектор направлений
=(с 1 ,с 2), который указывает направление наибольшего возрастания целевой функцииz =с 1 х 1 +с 2 х 2 .

В данной задаче вектор направлений
= (1, 4): он начинается в точкеО (0,0) и заканчивается в точкеN (1, 4).

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

Таким образом, точкой максимума целевой функции z является точкаА пересечения прямыхl 2 иl 3 .

Для вычисления оптимального значения целевой функции z найдем координаты точки А. Поскольку точка А – это точка пересечения прямых l 2 и l 3 , то ее координаты удовлетворяют системе уравнений, составленной из уравнений соответствующих граничных прямых:



Таким образом, точка А имеет координаты x 1 =1/6, x 2 = 35/6.

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

Подставив координаты точки А в целевую функцию (2.4), получим

max z = 1/6 + 4·(35/6) = 47/2.

Пример 2.2.2. Построить на плоскости область допустимых решений системы линейных неравенств (2.2.4) и найти наибольшее и наименьшее значения целевой функции (2.2.5):

(2.2.4)

z = –2x 1 –x 2 (2.2.5)

Решение. Для построения области допустимых решений, которая состоит из пересечения полуплоскостей, соответствующих каждому неравенству системы ограничений (2.2.4), запишем уравнения граничных прямых:

l 1: 4x 1 – x 2 = 0; l 2: x 1 + 3x 2 = 6; l 3: x 1 – 3x 2 = 6; l 4: x 2 = 1.

Прямая l 1 проходит через точку с координатами (0;0). Для ее построения выразим x 2 через x 1: x 2 = 4x 1 . Найдем еще одну точку, через которую проходит прямая l 1 , например (1;4). Через точку с координатами (0;0) и точку с координатами (1;4) проведем прямую l 1 .

Для приведения уравнения прямой l 2 к виду в отрезках на осях (2.2.3) разделим обе его части на 6:
. Таким образом, прямаяl 2 отсекает на оси Ох 1 6 единиц, на оси Ох 2 - 2 единицы. Аналогично имеем для l 3:
и Прямаяl 4 параллельна оси Ох 1 и проходит через точку с координатами (0;1) .

Для определения полуплоскостей, которые отвечают ограничениям системы (2.2.4) в ограничения нужно подставить координаты какой-либо точки, не лежащей на граничной прямой. В силу ограничений х 1 0, х 2 0, область допустимых решений ЗЛП лежит в первой четверти координатной плоскости.

О
бласть допустимых решений на рисунке 2.2.2 заштрихована.

Рисунок 2.2.2 – Область допустимых решений

Построим вектор направлений
= (–2,–1). Далее строим линию уровня, перпендикулярно к вектору.

Для нахождения наибольшего значения целевой функции передвигаем линию уровня в направлении вектора до последнего пересечения с областью допустимых решений. Таким образом, точкой максимума целевой функцииz является точкаА (пересечение прямыхl 1 иl 2).

Для вычисления оптимального значения целевой функции z найдем координаты точкиА . Поскольку точкаА – это точка пересечения прямыхl 1 иl 2 , то ее координаты удовлетворяют системе уравнений, составленной из уравнений соответствующих граничных прямых:



Таким образом, точка А имеет координаты x 1 =6/13, x 2 = 24/13.

Подставив координаты точки А в целевую функцию (2.2.5), получим оптимальное значение целевой функции

max z = – 2·(6/13) – (24/13) = – 36/13.

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

В результате решения ЗЛП возможны следующие случаи:

    Целевая функция достигает оптимального значения в единственной вершине многоугольника решений;

    Целевая функция достигает оптимальное значение в любой точке ребра многоугольника решений (ЗЛП имеет альтернативные опорные планы с одинаковыми значениями z);

    ЗЛП не имеет оптимальных планов;

    ЗЛП имеет оптимальный план в случае неограниченной области допустимых решений.

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

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

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

Однако для выработки наглядных представлений о решениях задач линейного программирования графический метод представляет определённый интерес. Кроме того, он позволяет геометрически подтвердить справедливость теорем линейного программирования .

Теоретические основы графического метода

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

при которых линейная форма принимает оптимальное значение.

Пример 3.

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

Продолжаем решать задачи графическим методом вместе

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

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

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

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

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

Наиболее простым и наглядным методом решения задачи линейного программирования (ЗЛП) является графический метод. Он основан на геометрической интерпретации задачи линейного программирования и применяется при решении ЗЛП с двумя неизвестными:

Будем рассматривать решение этой задачи на плоскости. Каждое неравенство системы функциональных ограничений геометрически определяет полуплоскость с граничной прямой а п х, + + a j2 х 2 = b n i = 1, т. Условия неотрицательности определяют полуплоскости с граничными прямыми х { = 0, х 2 = 0 соответственно. Если система совместна, то полуплоскости, пересекаясь, образуют общую часть, которая является выпуклым множеством и представляет собой совокупность точек; координаты каждой из этих точек являются решением данной системы. Совокупность этих точек называют многоугольником решений. Он может быть точкой, отрезком, лучом, ограниченным и неограниченным многоугольником.

Геометрически ЗЛП представляет собой отыскание такой угловой точки многоугольника решений, координаты которой доставляют максимальное (минимальное) значение линейной целевой функции, причем допустимыми решениями являются все точки многоугольника решений.

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

Определим, какую часть плоскости описывает неравенство 2х { + Зх 2 12.

Во-первых, построим прямую 2х, + Зх 2 = 12. Она проходит через точки (6; 0) и (0; 4). Во-вторых, определим, какая полуплоскость удовлетворяет неравенству. Для этого выбираем любую точку на графике, не принадлежащую прямой, и подставляем ее координаты в неравенство. Если неравенство будет выполняться, то данная точка является допустимым решением и полуплоскость, содержащая точку, удовлетворяет неравенству. Для подстановки в неравенство удобно использовать начало координат. Подставим х { = х 2 = 0 в неравенство 2х, + Зх 2 12. Получим 2 0 + 3 0

Аналогично графически можно изобразить все ограничения задачи линейного программирования.

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

Необходимо помнить, что область допустимых решений удовлетворяет условиям неотрицательности (Xj > 0, j = 1, п). Координаты любой точки, принадлежащей области определения, являются допустимым решением задачи.

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

Этот вектор показывает направление наискорейшего изменения целевой функции. Прямая c [ x l + с 2 х 2 = f(x 0), перпендикулярная вектору-градиенту, является линией уровня целевой функции (рис. 2.2.2). В любой точке линии уровня целевая функция принимает одно и то же значение. Приравняем целевую функцию постоянной величине а. Меняя значение а, получим семейство параллельных прямых, каждая из которых является линией уровня целевой функции.


Рис. 2.2.2.

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

Графический метод решения ЗЛП состоит из четырех этапов:

  • 1. Строится область допустимых решений (ОДР) ЗЛП.
  • 2. Строится вектор-градиент целевой функции (ЦФ) с началом в точке х 0 (0; 0): V = (с, с 2).
  • 3. Линия уровня CjXj + с 2 х 2 = а (а - постоянная величина) - прямая, перпендикулярная вектору-градиенту V, - передвигается в направлении вектора-градиента в случае максимизации целевой функции f(x v х 2) до тех пор, пока не покинет пределов ОДР. При минимизации /(*, х 2) линия уровня перемещается в направлении, противоположном вектору-градиенту. Крайняя точка (или точки) ОДР при этом движении и является точкой максимума (минимума) f(x p jc 2).

Если прямая, соответствующая линии уровня, при своем движении не покидает ОДР, то минимума (максимума) функции f(x р х 2) не существует.

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

4. Определяются координаты точки максимума (минимума). Для этого достаточно решить систему уравнений прямых, дающих в пересечении точку максимума (минимума). Значение f(x { , х 2), найденное в полученной точке, является максимальным (минимальным) значением целевой функции.

Возможные ситуации графического решения ЗЛП представлены в табл. 2.2.1.

Таблица 2.2.1

Вид ОДР

Вид оптимального решения

Ограниченная

Единственное решение

Бесконечное множество решений

Неограниченная

ЦФ не ограничена снизу

ЦФ не ограничена сверху

Единственное решение

Бесконечное множество решений

Единственное решение

Бесконечное множество решений

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

Намечается выпуск двух видов костюмов - мужских и женских. На женский костюм требуется 1 м шерсти, 2 м лавсана и 1 человекодень трудозатрат; на мужской - 3,5 м шерсти, 0,5 м лавсана и 1 человекодень трудозатрат. Всего имеется 350 м шерсти, 240 м лавсана и 150 человекодней трудозатрат.

Требуется определить, сколько костюмов каждого вида необходимо сшить, чтобы обеспечить максимальную прибыль, если прибыль от реализации женского костюма составляет 10 ден. ед., а от мужского - 20 ден. ед. При этом следует иметь в виду, что необходимо сшить не менее 60 мужских костюмов.

Экономико-математическая модель задачи

Переменные : х, - число женских костюмов; х 2 - число мужских костюмов.

Целевая функция :

Ограничения :

Первое ограничение (по шерсти) имеет вид х { + 3,5х 2 х { + 3,5х 2 = 350 проходит через точки (350; 0) и (0; 100). Второе ограничение (по лавсану) имеет вид 2х { + 0,5х 2 2х х + 0,5х 2 = 240 проходит через точки (120; 0) и (0; 480). Третье ограничение (по труду) имеет вид х у +х 2 150. Прямая х { + х 2 = 150 проходит через точки (150; 0) и (0; 150). Четвертое ограничение (по количеству мужских костюмов) имеет вид х 2 > 60. Решением этого неравенства является полуплоскость, лежащая выше прямой х 2 = 60.

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

На рис. 2.2.3 затенена область допустимых решений (ОДР). Для определения направления движения к оптимуму построим вектор- градиент V, координаты которого являются частными производными целевой функции:

Чтобы построить такой вектор, нужно соединить точку (10; 20) с началом координат. Для удобства можно строить вектор, пропорциональный вектору V. Так, на рис. 2.2.3 изображен вектор (30; 60).

Затем построим линию уровня 10xj + 20х 2 = а. Приравняем целевую функцию постоянной величине а. Меняя значение а , получим семейство параллельных прямых, каждая из которых является линией уровня целевой функции.