Границы области

Вид материалаЗадача

Содержание


Границы области
F(x) = const
Подобный материал:
1   2   3

Границы области


Обозначим границы области многоугольника решений.


Целевая функция F(x) → max


Рассмотрим целевую функцию задачи F = 5X1+6X2 → max.
Построим прямую, отвечающую значению функции F = 0: F = 5X1+6X2 = 0. Будем двигать эту прямую параллельным образом. Поскольку нас интересует максимальное решение, поэтому двигаем прямую до последнего касания обозначенной области. На графике эта прямая обозначена пунктирной линией.



Равный масштаб



Конец формы


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

Прямая F(x) = const пересекает область в точке B. Так как точка B получена в результате пересечения прямых (2) и (5), то ее координаты удовлетворяют уравнениям этих прямых:
10x1+6x2≤60
x1=0
Решив систему уравнений, получим: x1 = 0, x2 = 10
Откуда найдем максимальное значение целевой функции:
F(X) = 5*0 + 6*10 = 60

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

Определим максимальное значение целевой функции F(X) = 5x1 + 6x2 при следующих условиях-ограничений.

3x1 + 2x2≥6

10x1 + 6x2≤60

11x1≤11

Для построения первого опорного плана систему неравенств приведем к системе уравнений путем введения дополнительных переменных (переход к канонической форме).

В 1-м неравенстве смысла (≥) вводим базисную переменную x3 со знаком минус. В 2-м неравенстве смысла (≤) вводим базисную переменную x4. В 3-м неравенстве смысла (≤) вводим базисную переменную x5.

3x1 + 2x2-1x3 + 0x4 + 0x5 = 6

10x1 + 6x2 + 0x3 + 1x4 + 0x5 = 60

11x1 + 0x2 + 0x3 + 0x4 + 1x5 = 11

Введем искусственные переменные x: в 1-м равенстве вводим переменную x6;

3x1 + 2x2-1x3 + 0x4 + 0x5 + 1x6 = 6

10x1 + 6x2 + 0x3 + 1x4 + 0x5 + 0x6 = 60

11x1 + 0x2 + 0x3 + 0x4 + 1x5 + 0x6 = 11

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

F(X) = 5x1+6x2 - Mx6 => max

Из уравнений выражаем искусственные переменные:

x6 = 6-3x1-2x2+x3

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

F(X) = (5+3M)x1+(6+2M)x2+(-1M)x3+(-6M) => max

Матрица коэффициентов A = a(ij) этой системы уравнений имеет вид:


Решим систему уравнений относительно базисных переменных:

x6, x4, x5,

Полагая, что свободные переменные равны 0, получим первый опорный план:

X1 = (0,0,0,60,11,6)

Базис

B

x1

x2

x3

x4

x5

x6

x6

6

3

2

-1

0

0

1

x4

60

10

6

0

1

0

0

x5

11

11

0

0

0

1

0

F(X0)

-6M

-5-3M

-6-2M

1M

0

0

0


Переходим к основному алгоритму симплекс-метода.

Итерация №0.

Текущий опорный план неоптимален, так как в индексной строке находятся отрицательные коэффициенты.

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

Вычислим значения Di по строкам как частное от деления: bi / ai1

и из них выберем наименьшее:

min (6 : 3 , 60 : 10 , 11 : 11 ) = 1

Следовательно, 3-ая строка является ведущей.

Разрешающий элемент равен (11) и находится на пересечении ведущего столбца и ведущей строки.

Базис

B

x1

x2

x3

x4

x5

x6

min

x6

6

3

2

-1

0

0

1

2

x4

60

10

6

0

1

0

0

6

x5

11

11

0

0

0

1

0

1

F(X1)

-6M

-5-3M

-6-2M

1M

0

0

0

0



Получаем новую симплекс-таблицу:

Базис

B

x1

x2

x3

x4

x5

x6

x6

3

0

2

-1

0

-3/11

1

x4

50

0

6

0

1

-10/11

0

x1

1

1

0

0

0

1/11

0

F(X1)

5-3M

0

-6-2M

1M

0

5/11+3/11M

0


Итерация №1.

Текущий опорный план неоптимален, так как в индексной строке находятся отрицательные коэффициенты.

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

Вычислим значения Di по строкам как частное от деления: bi / ai2

и из них выберем наименьшее:

min (3 : 2 , 50 : 6 , - ) = 11/2

Следовательно, 1-ая строка является ведущей.

Разрешающий элемент равен (2) и находится на пересечении ведущего столбца и ведущей строки.

Базис

B

x1

x2

x3

x4

x5

x6

min

x6

3

0

2

-1

0

-3/11

1

11/2

x4

50

0

6

0

1

-10/11

0

81/3

x1

1

1

0

0

0

1/11

0

-

F(X2)

5-3M

0

-6-2M

1M

0

5/11+3/11M

0

0



Получаем новую симплекс-таблицу:

Базис

B

x1

x2

x3

x4

x5

x6

x2

11/2

0

1

-1/2

0

-3/22

1/2

x4

41

0

0

3

1

-1/11

-3

x1

1

1

0

0

0

1/11

0

F(X2)

14

0

0

-3

0

-4/11

3+1M


Итерация №2.

Текущий опорный план неоптимален, так как в индексной строке находятся отрицательные коэффициенты.

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

Вычислим значения Di по строкам как частное от деления: bi / ai3

и из них выберем наименьшее:

min (- , 41 : 3 , - ) = 132/3

Следовательно, 2-ая строка является ведущей.

Разрешающий элемент равен (3) и находится на пересечении ведущего столбца и ведущей строки.

Базис

B

x1

x2

x3

x4

x5

x6

min

x2

11/2

0

1

-1/2

0

-3/22

1/2

-

x4

41

0

0

3

1

-1/11

-3

132/3

x1

1

1

0

0

0

1/11

0

-

F(X3)

14

0

0

-3

0

-4/11

3+1M

0



Получаем новую симплекс-таблицу:

Базис

B

x1

x2

x3

x4

x5

x6

x2

81/3

0

1

0

1/6

-5/33

0

x3

132/3

0

0

1

1/3

-1/33

-1

x1

1

1

0

0

0

1/11

0

F(X3)

55

0

0

0

1

-5/11

1M