Решение оптимизационной задачи линейного программирования (стр. 5 из 7)

Таким образом, для исследования влияния изменения запаса ресурса на оптимальное решение нет необходимости решать задачу заново (с новым ограничением). Для нахождения оптимального решения достаточно по окончательной симплекс-таблице исходной задачи составить уравнения и подставить в них величину изменения запаса ресурса (значение d).

Изменение запасов ресурсов (т.е. правых частей ограничений) может привести к недопустимости оптимального базиса, найденного для исходной задачи. Так как на все переменные, используемые в задаче, накладывается требование неотрицательности, допустимый диапазон изменения запаса ресурса (т. е. диапазон допустимых значений d) находят из системы неравенств. Таким образом, допустимый диапазон изменения запаса времени работы токарного станка, при котором состав переменных в базисе оптимального решения не изменяется, находится из условия:

Х3 = 8 + 1*d > 0

Х6 = 0 – 0,5*d > 0

Х4 = 2,67 + 0,17*d > 0

Х5 = 5,33 + 0,33*d > 0

Решив данную систему неравенств, получим, что –8< d < 0. Таким образом, базис оптимального решения будет состоять из переменных (Х3,Х6,Х4,Х5), если запас времени работы токарного станка будет находиться в диапазоне от 0 до 8 часов. Выход значения d за границы этого диапазона приведет к недопустимости найденного нами оптимального решения, так как минимум одна из базисных переменных окажется отрицательной, и для того, чтобы найти оптимальное решение, нам придется решать задачу заново.

Аналогично выполняется анализ на чувствительность к изменению запаса времени работы станка-автомата.

5.4. АНАЛИЗ НА ЧУВСТВИТЕЛЬНОСТЬ К ИЗМЕНЕНИЯМ КОЭФФИЦИЕНТОВ ЦЕЛЕВОЙ ФУНКЦИИ

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

6. ОПРЕДЕЛЕНИЕ ОПТИМАЛЬНОГО ЦЕЛОЧИСЛЕННОГО РЕШЕНИЯ

Данная задача по своему содержанию является частично целочисленной. Переменные X1 , X2 , X3 , X4 , X5 , X6 ,обозначающие время работы определенного станка над деталями определенного типа, должны принимать целые значения. В то же время, переменные Х7 , Х8, обозначающие время простоя соответственно токарного станка и станка-автомата, могут принимать дробные значения. Для поиска оптимального целочисленного решения воспользуемся методом Гомори для частично целочисленных задач.

6.1. МЕТОД ГОМОРИ ДЛЯ ЧАСТИЧНО ЦЕЛОЧИСЛЕННЫХ ЗАДАЧ

Метод Гомори для нахождения целочисленного решения относится к большой группе методов, называемых методами отсечений. Эти методы основаны на введении в задачу дополнительных ограничений, позволяющих учесть требование целочисленности. Основная идея методов отсечений состоит в том, что на полученное оптимальное нецелочисленное решение накладывается дополнительное ограничение, которое делает это решение недопустимым, но и не отсекает ни одного целочисленного решения от области допустимых решений.

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

L1*W1 + L2*W2 + … +Ln*Wn ≥ {Bi} ,где

Aij, если Aij≥0 и Wjможет быть дробной, (1)

({Bi}*Aij)/({Bi}-1), если Aij<0 и Wjможет быть дробной, (2)

Lj = {Aij}, если{Aij}£{Bi} и Wiдолжна быть целой, (3)

{Bi}*(1-{Aij})/(1-{Bi}), если{Aij}>{Bi} и Wiдолжна быть целой, (4)

j=1,…,n

где Wn – небазисная переменная;

Bi - базисная переменная, имеющая максимальную дробную часть ( дробная часть числа – это разность между этим числом и максимальным целым числом, не превосходящим его);

Aij – коэффициент, стоящий на пересечении строки i-ой базисной переменной и столбца j-ой небазисной переменной;

Далее полученное ограничение приводится к стандартному виду:

-L1*W1 - L2*W2 - … -Ln*Wn + Sr = -{Bi}

где r – номер итерации алгоритма.

Здесь Sr – неотрицательная остаточная переменная, не имеющая никакого содержательного смысла; в оптимальном целочисленном решении эта переменная оказывается равной нулю.

В нашем случае переменная, имеющая максимальную дробную часть – это Х4 ({2,67}=0,67), она должна быть целой, переменные Х7и Х8 могут быть дробными, переменные Х1 и Х2 должны быть целыми, поэтому, согласно выше приведенной формуле, составим новое дополнительное ограничение. Так как все коэффициенты на пересечениях базисной переменной Х4 и небазисных переменных Х1 , Х2 , Х7 , Х8 ≥ 0(0,44≥0, 0,11≥0, 0,17≥0), токоэффициенты при переменных Х1 и Х2 рассчитали по формуле (3):L1={0,44}=0,44, L2={0,11}=0,11, а коэффициенты при переменных Х7 и Х8 рассчитали по формуле (1):L3=0,17, L4=0,17. {В4}={Х4} = {2,67} = 0,67. Ограничение будет иметь вид:

0,44Х1 + 0,11Х2 + 0,17Х7 + 0,17Х8≥ 0,67

Можно убедиться, что это ограничение сделало наше оптимальное решение недопустимым ( если подставить Х1=0, Х2=0, Х7=0, Х8=0, - значения переменных, полученных в оптимальном нецелочисленном решении, то получим 0≥0,67 – неверно).

Приведя ограничение к стандартному виду, имеем:

-0,44Х1 - 0,11Х2 - 0,17Х7 - 0,17Х8 + Х9= -0,67

Добавим к нашей финальной симлекс-таблице строку и столбец, соответствующие построенному ограничению и новой базисной переменной Х9:

БП	X1	X2	X3	X4	X5	X6	X7	X8	X9	БР
E	1,67	1,67	0	0	0	0	2,5	2,5	0	40
X3	1	1	1	0	0	0	1	0	0	8
X6	-0,67	-0,67	0	0	0	1	-0,5	0,5	0	0
X4	0,44	0,11	0	1	0	0	0,17	0,17	0	2,67
Х5	0,22	0,55	0	0	1	0	0,33	0,33	0	5,33
X9	-0,44	-0,11	0	0	0	0	-0,17	-0,17	1	-0,67

Таблица 8. Симплекс-таблица №7.

Как видно, полученная симплекс-таблица содержит недопустимое решение (переменная Х9имеет отрицательное значение). Произведем дальнейший пересчет таблицы, причем ведущую строку определяем максимальным по модулю отрицательным элементом столбца решений, а ведущий столбец – минимальным по модулю отношением элемента Е-строки к отрицательным элементам ведущей строки. Пересчет симплекс-таблицы осуществляется на основе стандартных процедур симплекс-метода.

Итак, переменная, исключаемая из базиса – это X9, т.к. ее значение –0,67 - это максимальный по модулю отрицательный элемент столбца решений. В базис включаем переменную X1, т.к. |1,67/(-0,44)|=3,8, |1,67/(-0,11)|=15,2, |2,5/(-0,17)|=14,7, 3,8 – минимальное по модулю отношение элемента Е-строки к отрицательным элементам ведущей строки. Ведущий элемент равен –0,44. Получим новую симплекс-таблицу:

БП	X1	X2	X3	X4	X5	X6	X7	X8	X9	БР
E	0	1,25	0	0	0	0	1,875	1,875	3,75	37,5
X3	0	0,75	1	0	0	0	0,625	-0,375	2,25	6,5
X6	0	-0,5	0	0	0	1	-0,25	0,75	-1,5	1
X4	0	0	0	1	0	0	0	0	1	2
Х5	0	0,5	0	0	1	0	0,25	0,25	0,5	5
X1	1	0,25	0	0	0	0	0,375	0,375	-2,25	1,5

Таблица 9. Симплекс-таблица №8.

Все значения базисных переменных стали неотрицательными, это означает остановку вычислительного процесса на данной итерации и анализ полученных результатов. Как видно из таблицы, в базис вошла новая переменная Х1, переменные Х3, Х4 и Х5уменьшили свое значение, а переменная Х6увеличилась. Значение целевой функции уменьшилось и стало равно 37,5 , что объясняется тем, что оптимальное нецелочисленное решение было отсечено нашим дополнительным ограничением, и для поиска оптимального целочисленного решения мы ушли вглубь области допустимых решений, где значение целевой функции меньше оптимального. Наше решение все еще нецелочисленное, поэтому составим новое ограничение.

Переменная, имеющая максимальную дробную часть – это Х3 ({6,5}=0,5) (Х1имеет такую же дробную часть, поэтому выбрали любую из них, например, Х3), она должна быть целой, переменные Х7, Х8 и Х9 могут быть дробными, переменная Х2 должна быть целой, поэтому, согласно формуле, составим новое дополнительное ограничение. Так как коэффициенты на пересечениях базисной переменной Х3 и небазисных переменных Х2 , Х7 , Х9 ≥ 0(0,75≥0, 0,625≥0, 2,25≥0), токоэффициент при переменной Х2 рассчитаем по формуле (3):L1={0,75}=0,75, коэффициенты при переменных Х7 и Х9 рассчитаем по формуле (1):L3=0,625, L4=2,25. Так как коэффициент на пересечении базисной переменной Х3 и небазисной переменной Х8<0, то коэффициент при переменной Х8 рассчитаем по формуле (2): L2=({6,5}*(-0,375))/({6,5}-1)=0,375. {В3}={Х3} = {6,5} = 0,5. Ограничение будет иметь вид: