Курсовая работа: Выбор параметров контроля с использованием метода динамического программирования и метода ветвей и границ
Ls ’ = ∑сj + h1 l ∆b1 ,
j=k+1
l-1
∆b1 = b1 - ∑ a1 j ∙xj - ∑a1 j .
xj ЄSj=k+1
Из условий следует, что Ls ’ не меньше максимального значения величины
∑cj ∙xj
xj ЄGS
при ограничениях
∑ a1j ∙xj b1 - ∑ a1j ∙xj = b1 ’ ,
xj ЄGS xj ЄS
xj Є {0;1}, xj ЄGS ,
Выбор очередной переменной для включения во множество S производится с помощью условия
h1r (xr ) = max h1j (xj )
xj ЄGS
Для выбранной переменной xr определяются величины Hs (xr ) и Hs (xr ), т.е. в S включаются xr = 1 или xr = 0.
Если в процессе решения окажется, что во множестве GS нет элементов, которые могут быть введены во множество S без нарушения ограничения (5), то полученное решение Ls =∑cj ∙xj принимается в качестве первого приближенного xj ЄSрешения L0 .
Все вершины дерева возможных вариантов, для которых выполняются условия
Hs ≤ L0 , из дальнейшего рассмотрения исключаются.
Из оставшихся ветвей выбирается ветвь с максимальным значением Hs , и процесс поиска оптимального варианта продолжается. Если в процессе решения будет найдено Ls = ∑cj ∙xj > L0 , то полученное решение принимается
xj ЄS
в качестве нового приближенного результата. Вычислительная процедура заканчивается, если для оставшихся ветвей выполняется условие Hs ≤ L0 .
2.2 Практическая часть
Ручной счёт
Данные для расчета:
С≤16
Таблица 4
N | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 |
Qi | 0.17 | 0.03 | 0.15 | 0.09 | 0.13 | 0.08 | 0.07 | 0.02 | 0.06 | 0.04 |
с(xi) | 5 | 1 | 4 | 2 | 6 | 3 | 2 | 3 | 1 | 1 |
hj | 0.034 | 0.03 | 0.0375 | 0.045 | 0.0217 | 0.0267 | 0.035 | 0.0067 | 0.06 | 0.04 |
Для удобства расчетов проранжируем таблицу 4 в порядке убывания hj и присвоим новые номера элементам, следующим образом:
Таблица 5
N | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 |
n | 9 | 4 | 10 | 3 | 7 | 1 | 2 | 6 | 5 | 8 |
Qi | 0.06 | 0.09 | 0.04 | 0.15 | 0.07 | 0.17 | 0.03 | 0.08 | 0.13 | 0.02 |
с(xi) | 1 | 2 | 1 | 4 | 2 | 5 | 1 | 3 | 6 | 3 |
hj | 0.06 | 0.045 | 0.04 | 0.0375 | 0.035 | 0.034 | 0.03 | 0.0267 | 0.02167 | 0.0067 |
Для формирования таблицы 6 произведем расчеты:
1)
8
∑сj =19>b1 - ∑ cj ∙xj =16-0=16;