Министерство образования Р.Ф.
Тюменский государственный нефтегазовый университет
Институт нефти и газаКафедра менеджмента
В отраслях ТЭК
Контрольная работа по
Дисциплине «Экономическая математика, методы и модели»
Вариант №4Выполнил: студент гр.
МОс2 Ваганова А.Р.
Проверил: Захаров А.В
Тюмень 2002 г.
Метод ветвей и границ. Рассмотрим задачу, состоящую в определении
максимального значения функции
[pic]
при условиях
[pic]
Как и при решении сформулированной задачи методом Гомори,
первоначально находим симплексным методом искусственного базиса оптимальный
план задачи без учета целочисленности переменных. Пусть им является план
X0. Если среди компонент этого плана нет дробных чисел, то тем самым
найдено искомое решение задачи и [pic].
Если же компонент плана Х0 имеются дробные числа, то Х0 не
удовлетворяет условию целочисленности и необходимо осуществить
упорядоченный переход к новым планам, пока не будет найдено решение задачи.
Покажем, как это можно сделать, предварительно отметив, что [pic] для
всякого последующего плана Х.
Предполагая, что найденный оптимальный план Х0 не удовлетворяет
условию целочисленности переменных, тем самым считаем, что среди его
компонент есть дробные числа. Пусть например переменная [pic] приняла в
плане Х0 дробное значение. Тогда в оптимальном целочисленном плане её
значение будет по крайней мере либо больше, либо меньше или равно
ближайшему меньшему целому числу [pic], либо больше или равно ближайшему
большему целому числу [pic]. Определяя эти числа, находим симплексным
методом решение двух задач линейного программирования:
[pic]Найдем рассмотренными выше методами решение задач линейного
программирования (I) и (II). Очевидно, здесь возможен один из следующих 4:
1. Одна из задач неразрешима, а другая имеет целочисленный оптимальный
план. Тогда этот план и значение целевой функции на нём и дают
решение исходной задачи.
2. Одна из задач неразрешима, а другая имеет оптимальный план, среди
компонент которого есть дробные числа. Тогда рассматриваем вторую
задачу и в ее оптимальном плане выбираем одну из компонент,
значение которой равно дробному числу, и строим две задачи,
аналогичные задачам (I) и (II).
3. Обе задачи разрешимы. Одна из задач имеет оптимальный целочисленный
план, а в оптимальном плане второй задачи есть дробные числа. Тогда
вычисляем значение целевой функции на этих планах и сравниваем их
между собой. Если на целочисленном оптимальном плане значение
целевой функции больше, или равно ее значению на плане, среди
компонент которого есть дробные числа, то данный целочисленный план
является оптимальным для исходной задачи и он вместе со значением
целевой функции на нем дает искомое решение.
Если же значение целевой функции больше на плане, среди
компонент которого есть дробные числа, то следует взять одно из
таких чисел и для задачи, план которой рассматривается, необходимо
построить две задачи, аналогичные (I) и (II).
4. Обе задачи разрешимы, и среди оптимальных планов обеих задач есть
дробные числа. Тогда выделяем значение целевой функции на данных
оптимальных планах и рассматриваем ту из задач, для которой
значение целевой функции является наибольшим. В оптимальном плане
этой задачи выбираем одну из компонент, значение которой является
дробным числом, и строим две задачи, аналогичные (I) и (II).
Таким образом, описанный выше интеграционный процесс может быть
представлен в виде некоторого дерева, на котором исходная вершина отвечает
оптимальному плану Х0 задачи (32)-(34), а каждая соединенная с ней ветвью
вершина отвечает оптимальным планам задач (I) и (II). Каждая из этих вершин
имеет свои ветвления. При этом на каждом шаге выбирается та вершина, для
которой значение является наибольшим. Если на некотором шаге будет получен
план, имеющий целочисленные компоненты, и значение функции на нем окажется
больше или равно, чем значение функции в других возможных для ветвления
вершинах, то данный план является оптимальным планом исходной задачи
целочисленного программирования и значение целевой функции на нем является
максимальным.
Итак, процесс нахождения решения задачи целочисленного
программирования (32)-(35) методом ветвей и границ включает следующие
этапы:
10 Находят решение задачи линейного программирования (32)-(34)
20 Составляют дополнительные ограничения для одной из переменных,
значение которой в оптимальном плане задачи (32)-(34) является дробным
числом.
30 Находят решения задач (I) и (II), которые получаются из задачи (32)-
(34) в результате присоединения дополнительных ограничений.
40 В случае необходимости составляют дополнительные ограничения
для переменной, значение которой является дробным, формулируют задачи,
аналогичные задачам (I) и (II), и находят их решение. Интеграционный
процесс продолжают до тех пор, пока не будет найдена вершина,
соответствующая целочисленному плану задачи (32)-(34) и такая, что значение
в этой вершине больше или равно значению функции в других возможных для
ветвления вершинах.
Описанный выше метод ветвей и границ имеет более простую логическую
схему расчетов, чем рассмотренный выше метод Гомори. Поэтому в большинстве
случаев для нахождения решения конкретных задач целочисленного
программирования с использованием ЭВМ применяется именно этот метод. В
частности в рассмотренном выше ППП «Линейное программирование в АСУ» для
отыскания целочисленного решения конкретных задач используется метод
ветвления и границ.
2.51 Методом ветвей и границ найти решение задачи, состоящей в
определении максимального значения функции
[pic]
при условиях[pic]
[pic]
xj – целые (j=[pic])
Р е ш е н и е. Находим решение сформулированной задачи симплексным
методом без учета условия целочисленности переменных. В результате
устанавливаем, что такая задача имеет оптимальный план Х0= (18/5, 3/5, 0,
0, 24/5). При этом плане F= (X0)=39/5.
Так как в плане Х0 значения трех переменных являются дробными числами,
то Х0 не удовлетворяет условию целочисленности, и следовательно, не
является оптимальным планом исходной задачи.
Возьмем какую-нибудь переменную, значение которой является дробным
числом, например х1. Тогда эта переменная в оптимальном плане исходной
задачи будет принимать значение, либо меньшее или равное трём:[pic], либо
больше или равно четырём: [pic].
Рассмотрим две задачи линейного программирования:
(I)[pic] (II)[pic]
Задача (I) имеет оптимальный план [pic] на котором значение целевой
функции [pic]. Задача (II) неразрешима.
Исследуем задачу (I). Так как среди компонент оптимального плана этой
задачи есть дробные числа, то для одной из переменных, например x2, вводим
дополнительные ограничения:
[pic]
Рассмотрим теперь следующие две задачи:
[pic](III) [pic]
(IV) [pic]
Задача (IV) неразрешима, а задача (III) имеет оптимальный план [pic](3,
1, 3, 3, 3), на котором значение целевой функции задачи [pic]
Таким образом исходная задача целочисленного программирования имеет
оптимальный план Х*= (3, 1, 2, 3, 3). При этом плане целевая функция
принимает максимальное значение [pic].
Схему реализованного выше вычислительного процесса можно представить в
виде дерева, ветвями которого являются соответствующие ограничения на
переменные, а вершинами – решения соответствующих задач линейного
программирования (рис 2.5).
Дадим геометрическую интерпретацию решения задачи (50)-(53). На рис.
2.6 показана область допустимых решений задачи (50)-(52).
Из него видно, что данная задача имеет оптимальный план[pic](18/5, 3/5,
0, 0, 24/5). В то же время [pic] не является планом задачи (50)-(53),
поскольку три переменные имеют дробные значения. Возьмем переменную х1 и
рассмотрим задачи (I) и (II).
Как видно из рис. 2.7задача (I) имеет оптимальный план [pic](3, 3/2,
0, 9/2, 3/2), а из рис.2.8 следует, что задача (II) неразрешима.
Поскольку среди компонент плана [pic]есть дробные числа, выберем
переменную х2 и рассмотрим задачи (III) (IV). Задача (III) имеет
оптимальный план[pic](3, 1, 2, 3, 3) (рис. 2.9), а задача (IV) неразрешима
(рис. 2.10).
Итак, Х*= (3, 1, 2, 3, 3) является оптимальным планом задачи (50)-(53).
При этом плане [pic].
Решение задачи, правые части которых содержат параметр.
Алгоритм решения задачи (60)-(62) подобен рассмотренному выше алгоритму
решения задачи (57)-(59).
Полагая значение параметра t равным некоторому числу t0, находим
решение полученной задачи линейного программирования (60)-(62). При данном
значении параметра t0 либо определяем оптимальный план, либо устанавливаем
неразрешимость задачи. В первом случае найденный план является оптимальным
для любого, где
[pic]
и числа qi и pi определены компонентами оптимального плана и зависят от
t0:
[pic]
Если при t = t0 задача (60)-(62) неразрешима, то, либо целевая
функция задачи (60) не ограничена на множестве планов, либо система
уравнений не имеет неотрицательных решений. В первом случае задача
неразрешима для всех [pic], а во втором случае определяем все значения
параметра [pic], для которых система уравнений (61) несовместна, и
исключаем их из рассмотрения.
После определения промежутка, в котором задача (60)-(62) имеет один и
тот же оптимальный план или неразрешима, выбираем новое значение параметра
t, не принадлежащее найденному промежутку, и находим решение полученной
задачи линейного программирования. При этом решение новой задачи ищем с
помощью действенного симплекс-метода. Продолжая итерационный процесс, после
конечного числа шагов получаем решение задачи (60)-(62).
Итак, процесс нахождения задачи (60)-(62) включает следующие основные
этапы:
10. Считая значение параметра t равным некоторому числу [pic],
находят оптимальный план или устанавливают неразрешимость полученной задачи
линейного программирования.
20. Находят значения параметра [pic], для которых задача (60)-(62)
имеет один и тот же оптимальный план или неразрешима. Эти значения
параметра t исключают из рассмотрения.
30. Выбирают значения параметра t из оставшейся части промежутка
[pic] и устанавливают возможность определения нового оптимального плана
находят его двойственным симплекс-методом.
40. Определяют множество значений параметра t, для которых задача
имеет один и тот же новый оптимальный план или неразрешима. Вычисления
проводят до тех пор, пока не будут исследованы все значения параметра
[pic].
2.66. Для каждого значения параметра [pic]найти максимальное значение
функции
[pic]
при условиях
[pic]
Р е ш е н и е . Считая значение параметра t в системе уравнений (81) равным
нулю, находим решение задачи (80)-(82) (табл. 2. 41).
Таблица 2.41
|i |Базис |Сб |Р0 |3 |-2 |5 |0 |-4 |
| | | | |Р1 |Р2 |Р3 |Р4 |Р5 |
|1 |Р3 |5 |12+t |1 |1 |1 |0 |0 |
|2 |Р4 |0 |8+4t |2 |-1 |0 |1 |0 |
|3 |Р5 |-4 |10-6t |-2 |2 |0 |0 |1 |
|4 | | |20+29t |10 |-1 |0 |0 |0 |
|1 |Р3 |5 |7+4t |2 |0 |1 |0 |-Ѕ |
|2 |Р4 |0 |13+t |1 |0 |0 |1 |Ѕ |
|3 |Р2 |-2 |5-3t |-1 |1 |0 |0 |Ѕ |
|4 | | |25+26t |9 |0 |0 |0 |Ѕ |Как видно из табл. 2.41, [pic]при t =0 есть оптимальный план задачи.
Однако [pic] является оптимальным планом и тогда среди его компонентов не
окажется отрицательных чисел, т.е. при 5-3t[pic]0; 7+4t[pic]0;
13+t[pic] или при [pic] Таким образом, если [pic] то[pic]- оптимальный план
задачи (80)-(82), при котором [pic]
Исследуем теперь, имеет ли задача оптимальные планы при [pic][pic].
Если [pic], то 5-3t<0 и следовательно, X=(0,5 – 3t, 7+4t, 13+t, 0) не
является планом задачи. Поэтому при [pic] нужно перейти к новому плану,
который был в то же время оптимальным. Это можно сделать в том случае,
когда в строке вектора Р2 имеются отрицательные числа [pic]. В данном
случае это условие выполняется. Поэтому переходим к новому опорному плану,
для чего введем в базис вектор Р1 и исключаем из него вектор Р2 (табл.
2.42).
Таблица 2.42[pic]
|i |Базис |Сб |Р0 |3 |-2 |5 |0 |-4 |
| | | | |Р1 |Р2 |Р3 |Р4 |Р5 |
|1 |Р3 |5 |17+2t |0 |2 |1 |0 |Ѕ |
|2 |Р4 |0 |18-2t |0 |1 |0 |1 |1 |
|3 |Р1 |3 |-5+3t |1 |-1 |0 |0 |-Ѕ |
|4 | | |70-t |0 |9 |0 |0 |5 |Как видно из табл. 2.42, [pic]-оптимальный план задачи для всех t, при
которых [pic] Следовательно, если [pic]является оптимальным планом исходной
задачи, причем [pic].
Если t>17/2, то [pic]не является планом задачи, так как третья
компонента 17 – 2t есть отрицательное число. Поскольку среди элементов 1-й
строки табл. 2.42 нет отрицательных при t>17/2 исходная задача неразрешима.
Исследуем теперь разрешимость задачи при t< -7/4. В этом случае Х=
(0,5 -3t, 7+4t, 13+t, 0) (см. табл.2.41) не является планом задачи, так как
третья компонента 7+4t есть отрицательное число. Чтобы при данном значении
параметра найти оптимальный план (это можно сделать, так как в строке
вектора Р3 стоит отрицательное число -1/2), нужно исключить из базиса
вектор Р3 и ввести в базис вектор Р5 (табл. 2.43).
Таблица 2.43
i |Базис |Сб |Р0 |3 |-2 |5 |0 |-4 | | | | | |Р1 |Р2 |Р3 |Р4 |Р5 | |1 |Р5 |-
4 |-14-8t |-4 |0 |-2 |0 |1 | |2 |Р4 |0 |20+5t |3 |0 |1 |1 |0 | |3 |Р2 |-2
|12+t |1 |1 |1 |0 |0 | |4 | | |32+30t |11 |11 |1 |0 |0 | |Как видно из
табл. 2.43, [pic] является оптимальным планом задачи для всех значений
параметра t, при которых [pic]
Таким образом, если [pic], то задача (80)-(82) имеет оптимальный план
[pic], при котором [pic]
Из табл. 2.43 так же видно, что при t<4 задача неразрешима, поскольку
в строке вектора Р4 нет отрицательных элементов.
Итак, если [pic], то задача не имеет оптимального плана; если
[pic]оптимальный план, а [pic] если [pic], то [pic]- оптимальный план, а
[pic]если [pic], то [pic]- оптимальный план, а [pic] если [pic], то задача
неразрешима.
