Файл: Исследование операций и методы оптимизации Заданиевариант 2 Выполнена обучающимся группы о. Издтв 23. 1Б22.docx

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 10.11.2023

Просмотров: 130

Скачиваний: 9

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Определим максимальное значение целевой функции Z(Y) = y1+y2 при следующих условиях-ограничений.
y1+8000y2≤1
6y1+6000y2≤1
Для построения первого опорного плана систему неравенств приведем к системе уравнений путем введения дополнительных переменных (переход к канонической форме).
y1+8000y2+y3 = 1
6y1+6000y2+y4 = 1
Решим систему уравнений относительно базисных переменных: y3, y4
Полагая, что свободные переменные равны 0, получим первый опорный план:
Y0 = (0,0,1,1)

Базис

B

y1

y2

y3

y4

y3

1

1

8000

1

0

y4

1

6

6000

0

1

Z(Y0)

0

-1

-1

0

0


Переходим к основному алгоритму симплекс-метода.
Итерация №0.
Текущий опорный план неоптимален, так как в индексной строке находятся отрицательные коэффициенты.
В качестве ведущего выберем столбец, соответствующий переменной y2, так как это наибольший коэффициент по модулю.
Вычислим значения Di по строкам как частное от деления: bi / ai2
и из них выберем наименьшее:
min (1 : 8000 , 1 : 6000 ) = 1/8000
Следовательно, 1-ая строка является ведущей.
Разрешающий элемент равен (8000) и находится на пересечении ведущего столбца и ведущей строки.

Базис

B

y1

y2

y3

y4

min

y3

1

1

8000

1

0

1/8000

y4

1

6

6000

0

1

1/6000

Z(Y1)

0

-1

-1

0

0






Формируем следующую часть симплексной таблицы. Вместо переменной y3 в план 1 войдет переменная y2.

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

Базис

B

y1

y2

y3

y4

y2

1/8000

1/8000

1

1/8000

0

y4

1/4

21/4

0

-3/4

1

Z(Y1)

1/8000

-7999/8000

0

1/8000

0



Итерация №1.
Текущий опорный план неоптимален, так как в индексной строке находятся отрицательные коэффициенты.
В качестве ведущего выберем столбец, соответствующий переменной y1, так как это наибольший коэффициент по модулю.
Вычислим значения Di по строкам как частное от деления: bi / ai1
и из них выберем наименьшее:
min (1/8000 : 1/8000 , 1/4 : 51/4 ) = 1/21
Следовательно, 2-ая строка является ведущей.
Разрешающий элемент равен (51/4) и находится на пересечении ведущего столбца и ведущей строки.

Базис

B

y1

y2

y3

y4

min

y2

1/8000

1/8000

1

1/8000

0

1

y4

1/4

21/4

0

-3/4

1

1/21

Z(Y2)

1/8000

-7999/8000

0

1/8000

0





Формируем следующую часть симплексной таблицы. Вместо переменной y4 в план 2 войдет переменная y1.

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

Базис

B

y1

y2

y3

y4

y2

1/8400

0

1

1/7000

-1/42000

y1

1/21

1

0

-1/7

4/21

Z(Y2)

401/8400

0

0

-999/7000

7999/42000



Итерация №2.
Текущий опорный план неоптимален, так как в индексной строке находятся отрицательные коэффициенты.
В качестве ведущего выберем столбец, соответствующий переменной y3, так как это наибольший коэффициент по модулю.
Вычислим значения Di по строкам как частное от деления: bi / ai3
и из них выберем наименьшее:
min (1/8400 : 1/7000 , - ) = 5/6
Следовательно, 1-ая строка является ведущей.
Разрешающий элемент равен (1/7000) и находится на пересечении ведущего столбца и ведущей строки.

Базис

B

y1

y2

y3

y4

min

y2

1/8400

0

1

1/7000

-1/42000

5/6

y1

1/21

1

0

-1/7

4/21

-

Z(Y3)

401/8400

0

0

-999/7000

7999/42000





Формируем следующую часть симплексной таблицы. Вместо переменной y2 в план 3 войдет переменная y3.

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

Базис

B

y1

y2

y3

y4

y3

5/6

0

7000

1

-1/6

y1

1/6

1

1000

0

1/6

Z(Y3)

1/6

0

999

0

1/6



Конец итераций: индексная строка не содержит отрицательных элементов - найден оптимальный план
Среди значений индексной строки нет отрицательных. Поэтому эта таблица определяет оптимальный план задачи.
Окончательный вариант симплекс-таблицы:

Базис

B

y1

y2

y3

y4

y3

5/6

0

7000

1

-1/6

y1

1/6

1

1000

0

1/6

Z(Y4)

1/6

0

999

0

1/6


Оптимальный план можно записать так:
y1 = 1/6, y2 = 0
Z(Y) = 1*1/6 + 1*0 = 1/6

Используя последнюю итерацию прямой задачи найдем, оптимальный план двойственной задачи.
x1=0, x2=1/6
Это же решение можно получить, применив теоремы двойственности.
Из теоремы двойственности следует, что X = C*A-1.
Составим матрицу A из компонентов векторов, входящих в оптимальный базис.

A = (A3, A1) =







1

1

0

6














Определив обратную матрицу D = А-1 через алгебраические дополнения, получим:

D = A-1 =







1

-1/6

0

1/6