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

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

Рассмотрим задачу ЛП в стандартной форме записи:

Положим n=2, т.е. рассмотрим эту задачу на плоскости. Пусть система неравенств совместна (имеет хотя бы одно решение).

Каждое неравенство этой системы геометрически определяет полуплоскость с граничной прямой a_i1 x₁ + a_i2 x₂ = b_i , i=1,2,…m. Условия неотрицательности определяют полуплоскости, соответственно, с граничными прямыми x1=0,x2 =0. Система совместна, поэтому полуплоскости, как выпуклые множества, пересекаясь, образуют общую часть, которая является выпуклым множеством и представляет собой совокупность точек, координаты каждой из которых являются решением данной системы. Совокупность этих точек называют многоугольником решений. Он может быть точкой, отрезком, лучом, многоугольником, неограниченной многоугольной областью.

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

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

неравенству. Удобной для использования при подстановке в неравенство является начало координат. Подставим х₁=х₂=0 в неравенство 2х₁+3х₂£12. Получим 2´0+3´0£12. Данное утверждение является верным, следовательно, неравенству 2х₁+3х₂£12 соответствует нижняя полуплоскость, содержащая точку (0.0). Это отражено на графике, изображенном на рис.1.

Рис. 1. Неравенству 2х₁+3х₂£12 соответствует нижняя полуплоскость.

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

Решением каждого неравенства системы ограничений ЗЛП является полуплоскость, содержащая граничную прямую и расположенная по одну сторону от нее. Пересечение полуплоскостей, каждая из которых определяется соответствующим неравенством системы, называется областью допустимых решений или областью определения. Необходимо помнить, что область допустимых решений удовлетворяет условиям неотрицательности (x_j ³0, j=1,…,n). Координаты любой точки, принадлежащей области определения являются допустимым решением задачи.

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

.

Этот вектор показывает направление наискорейшего изменения це­левой функции. Прямая , перпендикулярная вектору–градиенту, является линией уровня целевой функции. В любой точке линии уровня целевая функция принимает одно и то же значение. Приравняем целевую функцию постоянной величине “a”. Меняя значение “a”, получим семейство параллельных прямых, каждая из которых является линией уровня.

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

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

Графический метод решения ЗЛП состоит из следующих этапов.

1. Строится многоугольная область допустимых решений ЗЛП – ОДР,

2. Строится вектор-градиент ЦФ в какой-нибудь точке Х₀ принадлежащей ОДР –

.

3. Линия уровня C₁x₁+C₂x₂ = а (а–постоянная величина) - прямая, перпендикулярная вектору –градиенту – передвигается в направлении этого вектора в случае максимизации f(x_1,x₂) до тех пор, пока не покинет пределов ОДР. Предельная точка (или точки) области при этом движении и является точкой максимума f(x_1,x₂).

4. Для нахождения ее координат достаточно решить два уравнения прямых, получаемых из соответствующих ограничений и дающих в пересечении точку максимума. Значение f(x_1,x₂), найденное в получаемой точке, является максимальным.

При минимизации f(x_1,x₂) линия уровня перемещается в направлении, противоположном вектору-градиенту. Если прямая при своем движении не покидает ОДР, то соответствующий максимум или минимум f(x_1,x₂) не существует.

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

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

Задача. о планировании выпуска продукции пошивоч­ному предприятию.

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

Модель задачи.

Введем следующие обозначения: х₁ - число женских костюмов; x₂ - число мужских костюмов.

Прибыль от реализации женских костюмов составляет 10х₁, а от реализации мужских 20х₂, т.е. необходимо максимизировать целевую функцию

f(x) = 10´ х₁ + 20´ х₂ -> max.

Ограничения задачи имеют вид:

х₁ + х₂ £ 150

2 х₁ + 0.5 х₂ £ 240

х₁ + 3.5 х₂ £ 350

х₂³ 60

х₁ ³ 0

Первое ограничение по труду х₁ + х₂ £ 150. Прямая х₁ + х₂ = 150 проходит через точки (150, 0) и (0, 150).

Рис. 2. Решением первого неравенства является нижняя полуплоскость.

Второе ограничение по лавсану 2 х₁ + 0.5 х₂ £ 240. Прямая 2 х₁ + 0.5 х₂ = 240 проходит через точки (120, 0) и (0, 480). Третье ограничение по шерсти х₁ + 3.5 х₂ £ 350. Добавим четвертое ограничение по количеству мужских костюмов х2 ³ 60. Решением этого неравенства является полуплоскость, лежащая выше прямой х2 = 60. На рис.3. заштрихована область допустимых решений.

Рис. 3. Заштрихована область допустимых решений.

Для определения направления движения к оп­тимуму построим вектор-градиент Ñ, координаты которого являются частными производными целевой функции, т.е. = (10;20).

Что бы построить этот вектор, нужно соединить точку (10;20) с началом координат. При макси­мизации целевой функции необходимо двигаться в направ­лении вектора-градиента, а при минимизации — в противо­положном направлении. Для удобства можно строить век­тор, пропорциональный вектору Ñ. Так, на рис. 2.1.6. изобра­жен вектор градиент (30;60).

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

х₁ + 3.5 х₂ = 350

х₁ + х₂ = 150 .

Отсюда лег­ко записать решение исходной ЗЛП: max f(x) = 2300 и дости­гается при x₁=70 и x₂=80 (рис. 4.)

Рис.4. Максимум целевой функции достигается в точке (70, 80).

Задача.Для изготовления двух видов продукции А₁ и А₂ используют три вида ресурсов S₁, S₂, S₃, запасы которых составляют 18, 16 и 5 усл.ед. Расход ресурсов на 1 ед. продукции приведен в таблице:

Виды ресурсов	Запасы ресурсов	Расходы ресурсов на 1 изд.
		А1	А2
S1	18	1	3
S2	16	2	1
S3	5	-	1
	Прибыль	2 руб.	3 руб.

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

Составим экономико-математическую модель (ЭММ) задачи.

Пусть надо выпустить изделий A₁ - x₁ шт., а изделий А₂ - x₂ шт. Тогда прибыль F = 2x₁ + 3x₂ Þ max

x1 + 3x2 £	18
2x1 + x2 £	16
x2 £	5
x1 ³ 0,	x2 ³ 0

Решим задачу графически.

1)	x1 + 3x2 £ 18
	x1 + 3x2 = 18 (0; 6) (18; 0)
	к.т. (0; 0), 0 + 3*0 < 18 (в) – входит
2)	2x1 + x2 £ 16
	2x1 + x2 = 16 (0; 16) (8; 0)
	к.т. (0; 0), 2*0 + 0 < 16 (в) – входит
3)	x2 £ 5
	x2 = 5, x2 < 5 - ниже прямой
4)	x1 ³ 0 - правее ОX2
5)	x2 ³ 0 - выше ОX1
Линия уровня	F = 2x1 + 3 x2 F = 0
	2x1 + 3x2 = 0 (0; 0) (3; -2)
	q = {2; 3} - указывает направление возрастания F.

max F достигается в т. С

т.С	x1 + 3 x2 = 18	Þ -	2 x1 + 6 x2 = 36	Þ	5 x2 = 20	Þ
2x1 + x2 = 16	2 x1 + x2 = 16	x1 + 3 x2 = 18
Þ	x2 = 4
x1 = 6

Таким образом, необходимо выпустить x₁ = 6 шт. изделий А₁, x₂ = 4 шт. изделий А₂, чтобы получить max F = 2*6 + 3*4 = 24 ден.ед.