Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Bolshakova_Lineynoe_programmirovanie.pdf
Скачиваний:
103
Добавлен:
27.03.2015
Размер:
6 Мб
Скачать
☆

Ее начальный опорный план будет иметь вид

X0 = (0; 0; K; 0; b1; b2 ; K; bm ) ,

1442443

n

w(X0 ) = −b1 −b2 −K−bm .

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

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

Рассмотрим задачу линейного программирования симметричного вида относительно двух переменных:

z = c1x1 +c2 x2 → max,

a11x1 + a12 x2 ≤ b1,

a

21x1 + a22 x2 ≤ b2 ,

 

 

 

 

................................

 

 

 

 

a

x + a

x

≤ b ,

 

m1 1

m2 2

m

xj ≥ 0 , j = 1,2.

1.4.1. Геометрическая интерпретация области допустимых значений

(1.13)

(1.14)

(1.15)

1.Любое из неравенств (1.14) на плоскости Ох1х2 определяет некоторую полуплоскость.

2.Система неравенств (1.14)−(1.15) определяет выпуклое множество (выпуклый многоугольник, неограниченную выпуклую многоугольную область, пустую область или точку), которое совпадает с многогранником решений D.

16

1.4.2. Геометрическая интерпретация целевой функции

1.

Уравнение z = c1x1 +c2 x2

при

фиксированном значении

z = z0

определяет на плоскости

Ox1x2

прямую z0 = c1x1 +c2x2 .

При изменении z получают семейство параллельных прямых, назы-

ваемых линиями уровня.

2. Вектор коэффициентов целевой функции c = (c1; c2 ) называется градиентом функции. Он перпендикулярен линиям уровня.

3.

Градиент функции c = (c1;

c2 ) показывает направление наи-

большего возрастания целевой функции.

4.

Антиградиент −c = (−c1;

−c2 ) показывает направление

наибольшего убывания целевой функции.

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

Суть графического метода решения задач ЛП основывается на следующих утверждениях:

1) совокупность опорных планов задачи ЛП совпадает с системой вершин многогранника решений;

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

Для практического решения задачи (1.13) – (1.15) необходимо: 1) построить с учетом системы ограничений область допустимых

решений D (многогранник планов);r 2) построить вектор градиента c ;

3)построить перпендикулярно к нему в области допустимых решений одну из прямых семейства z = const;

4)искомая точка экстремума Xopt найдется параллельным пе-

ремещением вспомогательной прямой z = const в направлении вектора cr (если ищется zmax ) и в направлении вектора − c (если

ищется zmin );

5) координаты точки Xopt можно определить, решив совместно уравнения прямых, пересекающихся в этой точке, или по чертежу.

17

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

комбинацией этих вершин. Выпуклой линейной комбинацией точек X1 , X2 , …, Xk

k

 

 

 

 

называется сумма λ1X1 + λ2 X2 +K+ λk Xk , где ∑λi =1

,

λi ≥ 0 , i = 1, k .

i=1

 

 

 

 

1.4.4. Задача со многими переменными

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

1)выделить некоторый базис переменных в системе ограничительных уравнений;

2)опустить базисные переменные и перейти к эквивалентной системе неравенств;

3)выразить целевую функцию через свободные переменные;

4)полученную двумерную задачу рушить обычным графическим способом;

5)найдя две координаты оптимального решения, подставить их

вограничительные уравнения и определить остальные координаты оптимального плана.

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

z = 3x1 + 2x2 + x3 → max ,

x

+

2x

≤12,

1

+

2

= 4,

x1

+ x3

 

 

 

≤14,

2x1 + 2x2

xj ≥ 0, j =1,3.

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

выразим базисную переменную x3 :

x3 = 4 − x1 ,

18

подставим ее значение в целевую функцию:

z = 3x1 + 2x2 + x3 = 3x1 + 2x2 + 4 − x1 = 2x1 + 2x2 + 4 → max .

Так как x3 ≥ 0 , то 4 − x1 ≥ 0 , что равносильно неравенству

x1 ≤ 4 . Следовательно, система ограничений примет следующий симметричный вид:

x + 2x

≤12,

 

1

 

2

≤ 4,

x1

+ 2x

2x

≤14,

 

1

 

2

 

 

x1 ≥ 0 ,

x2 ≥ 0 .

Полученную двумерную задачу решим обычным графическим способом.

1. Так как x1 ≥ 0 , x2 ≥ 0 , то область допустимых решений будет находиться в первой координатной четверти. На плоскости Ox1x2 (рис. 1.1) построим прямые, порождаемые уравнениями:

x

+ 2x

=12

 

x1

+

 

x2

 

=1

−

прямая l

,

 

 

 

 

 

 

 

1

2

12

 

 

 

6

 

 

 

 

 

1

 

x

 

 

 

 

 

 

 

 

 

 

 

 

= 4

x1

=1

 

 

 

 

 

−

прямая l

,

 

 

 

 

 

 

 

1

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

x1

 

 

 

 

 

 

 

 

2x1 + 2x2

=14

 

+

 

 

 

 

 

=1

−

прямая l3.

7

 

7

 

 

(Уравнение прямой в отрезках

 

x1

 

+

x2

=1 означает, что прямая

 

a

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

линия отсекает на оси Ox1 отрезок длиной a, а на оси Ox2

− отре-

зок длиной b.)

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

19

x

+ 2x

≤12,

1

2

≤ 4,

x1

 

 

 

≤14.

2x1 + 2x2

Чтобы определить полуплоскость, соответствующую 1-му нера-

венству x1 +

2x2 ≤12 , возьмем

точку, не лежащую на прямой

x1 + 2x2 =12

(например, (0; 0)),

и подставим ее в неравенство

x1 +2x2 ≤12 ( 0 + 2 0 ≤ 12 0 ≤ 12 ). Если неравенство выполняется, то точка принадлежит данной полуплоскости, в противном

случае − не

принадлежит. В нашем

примере неравенство

x1 +2x2 ≤12

определяет полуплоскость,

лежащую ниже прямой

l1 : x1 +2x2 =12 .

 

3. Построим вектор градиента c = (2;

2) , т.е. соединим точки

начала (0; 0) и конца (2; 2) стрелкой. Перпендикулярно к нему построим одну из прямых семейства z = 2x1 +2x2 +4 = const . На-

пример, при x1 = x2 =1 получим z = 8. Заметим, что эта прямая параллельна l3 : 2x1 +2x2 =14 .

4. Вспомогательную прямую z = 8 будем параллельно перемещать в направлении вектора c (так как ищется zmax ) до последней точки пересечения с областью допустимых решений D. В нашем случае прямая zmax совпадает с прямой l3 : 2x1 +2x2 =14 . Поэто-

му искомых точек экстремума Xopt будет множество:

все точки

отрезка [ Xopt1 ;

Xopt2 ]. Координаты точек

Xopt1

и Xopt2

можно оп-

ределить по рис. 1.1: Xopt1 = (2; 5) и Xopt2

= (4; 3). Но так как пер-

воначальная задача имеет три неизвестные, то

x3 найдем, подста-

вив значение переменной x1 во 2-е ограничительное уравнение:

 

1) для X1

= (2; 5) получим: x = 4 − x = 4 −2 = 2 X1

=

opt

3

1

 

opt

 

= (2; 5; 2);

20

2) для Xopt2 = (4; 3) получим: x3 = 4 − x1 = 4 −4 = 0 Xopt2 = = (4; 3; 0).

Рис. 1.1

5. Для аналитической записи любого оптимального решения Xopt необходимо составить выпуклую линейную комбинацию

опорных решений

Xopt1 и

Xopt2 : Xopt = λ1Xopt1 + λ2 Xopt2 , где

λ1 +λ2 =1 и λ1 ≥ 0 , λ2 ≥ 0 . Следовательно,

 

Xopt = λ1(2; 5; 4) + λ2 (4; 3; 0) = (2λ1 + 4λ2 ; 5λ1 +3λ2 ; 2λ1) ,

где λ1 + λ2 =1 и λ1 ≥ 0,

λ2 ≥ 0 .

 

Тогда

 

 

 

z = z( X1 ) = 3 2 + 2 5 + 1 2 = 18 = z( X2

).

max

opt

opt

 

21

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]