- •Исследованиеопераций
- •1 Литература 3
- •2 Основныепонятияисследованияопераций(исо).Классификациязадач 4
- •3 Основылинейногопрограммирования 8
- •4 Основытеорииграфов 23
- •1 Литература
- •2Основныепонятияисследованияопераций(исо).
- •2.1КлассификациязадачИсо
- •2.2Многокритериальныезадачи
- •3Основылинейногопрограммирования
- •3.1ГрафическийметодрешениязадачЛп
- •3.2 ВозможныевариантырешениязадачЛп
- •3.3 ЭквивалентныепостановкизадачЛп
- •3.4Симплекс-метод
- •3.5ТранспортнаязадачаЛп
- •3.6Задачаоназначениях
- •4Основытеорииграфов
- •4.1Основныепонятия
- •4.2 Алгоритмпостроенияэйлеровацикладлянеориентированно-гографа
- •4.3Заданиеграфов
- •4.4Топологическаясортировка
- •4.5Остовноедерево(остов)минимальноговеса
- •4.6Кратчайшиепути
- •Ij (j, еслиw0
- •Ij 0,еслиw0
- •5Сетевоепланированиеиуправлениепроектами
- •6Задачакоммивояжераиметодветвейиграниц
- •7Имитационноемоделированиенапримереметода
- •8Основытеориисложностивычислений
- •9Динамическоепрограммирование
3.1ГрафическийметодрешениязадачЛп
Вначалерассмотримнаиболеепростойметодрешения,применимыйдляслучая,когдаколичествопеременныхнепревосходитдвух.Рассмотримследующийпример.
Компанияпроизводитпогрузчикиитележки.Отодногопогрузчикакомпанияполучаетдоходвразмере$80иотоднойтележкивразмере$40.Имеетсятриобрабатывающихцен-тра,накоторыхвыполняютсяоперацииметаллообработки,сваркиисборки,необходимыедляпроизводствалюбогоизпродуктов.Дляинтервалапланирования,равногомесяцу,за-данапредельнаяпроизводственнаямощностькаждогообрабатывающегоцентравчасах,атакжеколичествочасов,необходимоенаэтомцентредляпроизводстваодногопогрузчика
иоднойтележки.Этаинформациязаданавтаблице.
-
Центр\Изделие
Погрузчик Тележка
(часынаед.) (часынаед.)
Общаямощность
(часы)
Мет.обработка
Сварка
Сборка
6 4
2 3
9 3
2400
1500
2700
Требуетсясоставитьдопустимыйпланработнамесяцсмаксимальнымдоходом.Матема-тическаямодельзадачиможетбытьзаписанаследующимобразом.
80x1+40x2→max,
6x1
+4x2
≤2400
2x1+3x2≤1500
1 2
9x+3x≤2700x1≥0,x2≥0
гдеx1иx2–количествапроизводимыхпогрузчиковитележексоответственно.
Представимограниченияввидеследующегографика,гдепоосям(0x1)и(0x2)отклады-ваетсяколичествопроизведенныхпогрузчиковитележексоответственно,алиниипред-ставляютограничениянапроизводственныемощности.Допустимаяобластьзадачипред-ставленаввидемногоугольника.
Представимцелевуюфункцию80x1+40x2=zспеременнымзначениемzпрерывистойлинией.Любаяточка(x1,x2)налинии80x1+40x2=zсоответствуетдоходувразмереz.Перемещаяеепараллельносебесамой,получаемразныезначениядохода.
Приx1=0значениеz=40x2.Следовательно,значениеzувеличиваетсяприувеличенииx2,т.е.приперемещениицелевойфункциивверх.Приэтомлиния80x1+40x2=zнедолжнапокинутьдопустимуюобласть.
Нетруднозаметить,чтомаксимальныйдоходсредидопустимыхточекдостигаетсяводнойизвершинмногоугольникаограничений.Вданномслучаевточкепересечениялиний,
соответствующихограничениямнаметаллообработкуисборку,x∗=200,x∗=300.
1 2
Дляопределенияточки,соответствующейоптимальномурешению,нужносравнитьуглынаклонапрямых,соответствующихцелевойфункциииограничению.
3.2 ВозможныевариантырешениязадачЛп
1)ЗЛПимеетединственноерешение
x2
✻
900
❇
❇
❇
❇
❇
600 ❇
❏ ❇
500❏❇
❏
❇
◗◗❇
❏◗
❏◗❇◗
❇
❏◗
•❏❇ ◗
❇ ◗
❆
❆❆✟✟✯
❏ ◗
◗
❇❏ ◗❇❏ ◗
❇❏ ◗
✲
◗
❆
❆ ❇ ❏ ◗300400 750 x1
2)НесуществуетдопустимогорешенияЗЛП
3)ЗЛПимеет∞многооптимальныхрешений,которыеназываютсяальтернативными
4)ЦелеваяфункцияЗЛПнеограничена(вдопустимойобласти).
3.3 ЭквивалентныепостановкизадачЛп
Общаяпостановка:
n
j=1
n
cjxj→max(min),
j=1aijxj≤bi,i∈I1
n
n
j=1aijxj=bi,i∈I2 j=1aijxj≥bi,i∈I3
j 1
x≥0,j∈N,
xj∈R,j∈N2,xj≤0,j∈N3,
гдеcj,aij,bi–заданныечисловыепараметры,множестваI1,I2иI3попарнонепересе-каются,I1∪I2∪I3={1,...,m},R–множестводействительныхчисел,множестваN1,N2иN3попарнонепересекаются,иN1∪N2∪N3={1,...,n}.Cуществуетнесколько
эквивалентныхпостановокзадачЛП:
–Задачимаксимизациииминимизации.Дляпереходаотоднойзадачикдругойкоэффи-циентыцелевойфункцииможноумножитьна-1.
–Задачисограничениямитипа”=”,”≤”и”≥”.Дляпереходаотограничения”≤”кограни-
чению”≥”инаоборотдостаточнодомножатьсоответствующееограничениена-1.Огра-
ничениетипа”=”можетбытьпредставленоввидедвухограничений”≤”и”≥”стойже
левойиправойчастью.Дляпереходаотограниченийтипа”=”кограничениямтипа≤
можноприменитьметодГаусса.Например,рассмотримзадачу
−x1−x2+x5→min
x1+x2=1
x2−2x3=−3
−
x3 x4+x5=1
xj≥0
ПрипомощиметодаГауссавыделяемединичнуюподматрицувматрицеограничений.
-
1
1
0
0
0
1
1
0
2
0
0
4
1
0
0
2
-2
2
0
1
-2
0
0
-3
⇒
0
1
-2
0
0
-3
⇒
0
1
0
-2
2
-1
0
0
1
-1
1
1
0
0
1
-1
1
1
0
0
1
-1
1
1
Эквивалентнаязадача
x1=2−(2x4−2x5)≥0
x2=−1−(−2x4+2x5)≥0x3=1−(−x4+x5)≥0
−2+2x4−2x5−1+2x4−2x5+x5=−3+4x4−3x5→min
2x4
−2x5≤2
−2x4+2x5≤−1
−x4+x5≤1
x4,x5≥0
Дляпереходаотограниченийтипа≤кограничениямтипа=достаточноввестидополни-тельныенеотрицательныепеременные,пооднойдлякаждогоограничения≤.Этипере-
менныевойдутвцелевуюфункциюснулевымикоэффициентами.
Переменнаяxj∈Rможетбытьпредставленаввидедвухновыхнеотрицательныхпере-
менных:xj=x+−x−,x+≥0,x−≥0.
j j j j
ДлякаждойзадачиЛПможносформулироватьдвойственнуюей задачу.Рассмотримзадачумаксимизациисограничениями”≤”.Представимэтузадачувматричнойформе:
cx→max,Ax≤b,
гдеc=(c1,...,cn)–вектор-строка,x=(x1,...,xn)T–вектор-столбец,A=||aij||–матрицаразмерностиm×n,b=(b1,...,bm)T–вектор-столбец.Отметим,чтоограниченияx≥0отсутствуют.ОнимогутбытьвойтивограниченияAx≤bпутемвведенияэквивалентныхограничений−x≤0.Введемврассмотрениедвойственнуюейзадачу:
yb→min,
yA=c,y≥0,
гдеy=(y1,...,ym)–вектор-строкадвойственныхпеременныхилипотенциалов.ПустьAiобозначаетстрокуiматрицыA,аAj–столбецjэтойжематрицы.Соответствиемеждукомпонентамипрямойидвойственнойзадачпредставленовследую-щейтаблице.
-
Прямаязадача
Двойственнаязадача
max{cx}
Aix≤bi,i∈I1
Aix=bi,i∈I2
xj≥0,j∈N1
xj∈R,j∈N2
xj≤0,j∈N3
min{yb}
yi≥0,i∈I1
yi∈R,i∈I2
yAj≥cj,j∈N1
yAj=cj,j∈N2
yAj≤cj,j∈N3
Теорема1(Теоремадвойственности)Справедливоодноизследующихутверждений:
а)обезадачи,прямаяидвойственнаяей,имеютдопустимыерешенияи
max{cx}=min{yb},
б)еслиоднаиззадач,прямаяидвойственнаяей,неимеетдопустимогорешенияадругаяимеет,тоцелеваяфункциядругойзадачинеограничена.
в)обезадачинеимеютдопустимыхрешений.
Теорема2(Одополняющейнежесткости)Пустьxиy–допустимыерешенияпрямойидвойственнойзадачсоответственно.Тогдаследующиеусловияа),б)ив)эквивалент-ны:
а)xиyоптимальныерешенияпрямойидвойственнойзадач,б)cx=yb,
в)(условиедополняющейнежесткости)yi(Aix−bi)=0i=1,...,m,и(yAj−cj)xj=0,
j=1,...,n.
Изэтойтеоремы,например,следует,чтоеслирешенияпрямойидвойственнойзадачдопу-стимыивыполняетсяусловиев)дополняющейнежесткости,тоэтирешенияоптимальны.Оптимальноерешениепрямойзадачиможетбытьполученоизоптимальногорешениядвойственнойзадачи.Например,рассмотримзадачу
4x1−3x2−6x3−x4→max,
−2x1−x2+5x3+2x4≥2−3x1+x2−x3−x4≥1
xi≥0,i=1,2,3,4
Построимдвойственнуюейзадачу
2y1+y2→min,
−2y1−3y2≥4
−y1+y2≥−3
5y1−y2≥−6
2y1−y2≥−1yj≥0,j=1,2
Решимдвойственнуюзадачуграфически:y∗=(−9/4,−21/4),z∗=−39/4.Оптимальное
решениеопределяетсяпересечениемлиний,соответствующихограничению2и3двой-
1
ственнойзадачи.Этиограничениявыполняютсякакстрогиеравенства,аостальные-какнестрогие.Тогдапоусловиюдополняющейнежесткостидолжновыполнятьсяx∗=0,4=0,x2иx3могутбытьнайденыизравенств−x2+5x3=2иx2−x3=1(рассмат-
x∗ ∗ ∗
∗ ∗ ∗ ∗
риваемподматрицуматрицыпрямойзадачи,состоящуюизстолбцов2и3).Получаем
x∗=(0,7/4,3/4,0).
