Линейное программирование

Download Report

Transcript Линейное программирование

Линейное
программирование
Двойственность в линейном
программировании
Виды двойственных задач
Симметричная
 Не симметричная
 Смешанные

Алгоритм построения двойственной
симметричной задачи





L(x) = b1y1 + b2y2 +…+bmym → min
Каждому неравенству системы
ограничений исходной задачи поставим в
соответствие переменную yi
Построим целевую функцию, слагаемыми
которых являются произведения
свободных коэффициентов исходной
системы на соответствующие новые
Составим систему ограничений,
коэффициентами которой являются
коэффициентами матрицы,
транспонированной из коэффициентов
исходной задачи
Свободными членами в системе
ограничений являются коэффициенты
целевой функции исходной задачи. Знаки
неравенства меняются на
противоположные, целевая функция
стремиться к минимуму, все yi ≥0
a1 y1  a21 y2  ...  am1 ym  C1
a y  a y  ...  a y  C
 12 1
22 2
m2 m
2

   

a1m y1  a2 m y2  ...  anm ym  Cm
Не симметричные двойственные
задачи
Не симметричные двойственные задачи – это задачи, в которых
система ограничений является равенством
(задачи в каноническом виде)
Основное неравенство
двойственности
C1x1 +C2x2 +…+Cnxn ≤ b1y1 +b2y2 +…+bmym
Для произвольного допустимого плана исходной и двойственной задачи
Доказательство
Умножим каждое неравенство в двойственной задачи на соответствующую
переменную Xj и сложим все неравенства перегруппировав их
Теоремы
Если одна из двойственных задач имеет оптимальное решение,
то и вторая имеет такое же по значению целевой функции L(x)
оптимальное решение
 Для оптимальности допустимых решений x и y,пары
двойственных задач необходимо и достаточно выполнение
следующих ограничений
система

 yоптi  (ai1 x1  ai 2 x2  ...  ain xn  bi )  0

 xоптj  (a1 j y1  a2 j y2  ...  amj ym  C j )  0
Пример симметричной двойственной задачи

L(x) = x1-x2 → max
 2 x1  x2  2

 x1  2 x2  2
x  x  5
 1 2

xj  0
L(y) = 2y1 + 2y2 + 5y3 → min
 2 y1  y 2  y3  1

 y1  2 y 2  y3  1



yj  0
Пусть найдено решение исходной задачи.
L(x) опт=3
опт (4; 1)
Используя две теоремы двойственности, получим:
L(y) min = L(x) max = 3 по Теореме №1
 y1опт  (2 x1  x2  2)  0
y
 2 опт  ( x1  2 x2  2)  0
 y3опт  ( x1  x2  5)  0
 x  (2 y  y  y 1)  0
1
2
3
 1
 x2  ( y1  2 y 2  y3  1)  0
y1опт (9)  0  y1опт  0
y 2опт  0  0
y3опт  0  0
4  ( y 2  y3  1)  0
1  (2 y 2  y3  1)
y1  0
y 2  y3  1
 2 y1  y3  1
y2  2
3
y3  1
3
Пример симметричной двойственной задачи
Для определения решения двойственной задачи можно использовать способ, основанный
на взаимнооднозначном соответствие переменных в конечном решении симплекс таблицы.
Основными переменными исходной задачи соответствуют балансовые переменные
двойственной задачи и наоборот.
x1
x2
x3
x4
x5
bj
x3
2
1
1
0
0
2
x4
1
2
0
1
0
2
x5

1
1
1
1
0
0
0
0
1
0
5
0
x3
x1
x1
0
1
x2
0
0
x3
1
0
x5
0
1
0

0
0
0
x4
1
2
3
1
3
2
3
x5 b j
1 9
2
4
3
1
1
3
1
3
3
Решением исходной системы является x1=4 (второе уравнение в симплекс таблице), ей
соответствует введенная в исходной системе переменная x4 (y2). Элемент на пересечении
соответствующей строки и столбца дает:
аналогично
y3  1
3
y2  2
3
Экономический смысл двойственной задачи

L(y)= b1 y1  b2 y 2  ...  bm y m  min
L( y )
 yi
bi
L
 yi  L  yi bi
bi
Пусть ∆bi=1 => ∆L ≈ yi – для оптимального решения найденной значение yi обозначает
ценность рассматриваемого продукта
Yi -условная цена единицы i- того ресурса, с его помощью можно определить степень
влияния ограничений на значение целевой функции. Предельное значение (верхняя и
нижняя границы ресурсов) для которых yi остается неизменным определяется следующим
образом
н
b
i
 xj 

 min
d 
 ij 
в
b
i
 xj 

 max
d 
 ij 
x - значение в оптимальном решении
d - соотв. базису матрицы в оптимальном решении
j
ij
Если в план производства включаются новые виды продукции, то эффективность их
m
включения в план можно оценить по формуле
  a y  C
j
i 1
ij
i
j
Если ∆j < 0 , то план улучшается и новую продукцию выгодно выпускать
Если ∆j > 0 производство не выгодно