Задача линейного программирования. Проверка оптимальности текущего плана, страница 6

20) Так как ∆j ≥ 0 при j=1,2,3,4,5,значит текущий план Х*1=(2, 0, 6, 0, 1) оптимален, а максимальное значение функции Z(X) равно max Z(X) =Z(X1) = 4.

21) Проверка:  2 + 0 + 0(6) – 6(0) – 1 = 1 ,0(2) + 0 + 0(6) + 0 + 1(1) = 1 ,0(2)+0+ 6–3(0)–4(1) = 2 , х1 = 2 > 0 ,  x2 = 0 ,  x3 = 6 > 0 ,  x4 = 0 ,  x5 = 1 > 0.

Z(X) = - 1(2) + 0 + 6 = 4.

2°) Если столбцов Aj с положительными (отрицательными) оценками ∆j > 0 (∆j < 0) несколько, то находим сначала θj для каждого столбца с положительной (отрицательной) оценкой по формуле θj =      =    ,  , . . . ,    .

Выгода Vj находится по формуле Vj = ∆j∙θj.

Включить в базис тот столбец Aj, для которого выгода Vj максимальна. При этом ведущим элементом окажется тот элемент, аij, при котором получено число θj, то есть тот элемент,  а′ij, для которого отношение  - самое малое при переборе всех строк таблицы с аij>0.

1.2.5. Единственность оптимального плана. В задаче на минимум  или на максимум, оптимальное решение является единственным, если ∆j ≠ 0 при всех j = m+1, m+2, …, n.Здесь предполагается, что Am+1, Am+2, …, An небазисных столбцов.

Пример 6.        Z (Х) = 5х1 - 2х2 – 3х3 + 0х4 + 0х5 + 0x6 → min,

2x1 + x2 – x3 + x4 + 0x5 + 0x6 = 5,

3x1 +2x2 + x3 +0x4 + x5 + 0x6 = 7,

x1 + 0x2 – 3x3 +0x4 +0x5 +x6 = 4,

xi ≥ 0, i = 1, 2, 3, 4, 5, 6.


Решение. 1)Составим симплексную таблицу.

Б

Сб

В

5

-2

-3

0

0

0

A1

A2

A3

A4

A5

A6

1

A4

с4= 0

b1=5

2

1

-1

1

0

0

   +

2

A5

с5= 0

b2=7

3

2

1

0

1

0

(1) (1) (3)

3

A6

с6= 0

b 3=4

1

0

-3

0

0

1

             +

4

 

Z=0

-5

2

3

0

0

0

2) Матрица. А размера 3×6 имеет 3 базисных столбца А4 = ,  А5 = ,   А6 =  .

3) Значения базисных переменных с теми же номерами, что имеют базисные столбцы, берутся равными соответствующим элементам столбца В, то есть  х4= b1=5,   х5= b2 =7,   х6= b3 =4.

Значения небазисных переменных берутся равными нулю, то есть  х1 = х2 = х3 = 0.

4) Первоначальный план будет X0=(х1= 0,х2 = 0, х3 = 0, х4 = 5, х5 =7, х6 = 4). Оптимальность плана X0 можно проверить по значениям  1, ∆2, ∆3.

Эти оценки находятся по формуле: ∆j = CбAj – cj,  j = 1, 2, 3.

Если  ∆j ≤ 0  при всех значениях  j = 1, 2, 3,  то текущий план Х0 оптимален:  Х0 = Х*, и тогда

min Z(Х)=Z(Х0)=Z(Х*)= 5∙(0)–2∙(0) – 3∙(0) + 0∙(5) + 0∙(7) + 0∙(4) = 0.

Это число Z = 0  стоит внизу под столбцом В. Если же хотя бы одна оценка  ∆j > 0,  то текущий план Х0 не оптимален, его можно улучшить, то есть перейти к другому текущему плану Х1, для которого значение Z(Х1) окажется меньше, чем  Z = 0.Оценки базисных столбцов  А4, А5 и А6  равны  ∆5 = ∆6 = ∆7 = 0. Эти нули можно записывать в таблицу сразу, не вычисляя. Находим оценки небазисных столбцов  А1, А2 и А3.

1 = Cб∙A1 – c1 = (0; 0; 0)∙(2;3;1)–5 = (0)∙(2) + (0)∙(3) + (0)∙(1)–5 = 0 + 0 + 0 –5= -5 < 0,

2 = Cб∙A2 – c2 = (0;0;0)∙(1;2;0) – (-2) =(0)∙(1) + (0)∙(2) + (0)∙(0)–(2) = 0 + 0 + 0 + 2 = 2 > 0,

3 = Cб∙A3 – c3 = (0;0;0)∙(-1;1;-3) – (-3) = (0)∙(-1) + (0)∙(1) + (0)∙(-3)–(-3) =  0+0+0 – (-3) = 3 > 0.

5) Так как  ∆2 = +2 > 0,  ∆3 = +3 > 0, то текущий план Х0 неоптимален. Для улучшения текущего плана Х0  находим сначала θj для каждого столбца с положительной оценкой по формуле