Файл: Исследование операций и методы оптимизации Заданиевариант 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) = |
| |
Определив обратную матрицу D = А-1 через алгебраические дополнения, получим:
D = A-1 = |
| |