Badania operacyjne 1 Badania Operacyjne Operations Research Management
Badania operacyjne 1
Badania Operacyjne (Operations Research, Management Science) • Badania Operacyjne (BO) należą do matematycznych nauk interdyscyplinarnych zajmujących się efektywnym wykorzystaniem środków przez różnego typu organizacje. • Istotne znaczenie w BO ma interakcja pomiędzy człowiekiem a technologią – nacisk na praktyczne zastosowania metod matematycznych. 2
Badania operacyjne – zakres metod BO korzystają z narzędzi, m. in. : • Rachunku prawdopodobieństwa, • Statystyki, • Ekonometrii, • Metod optymalizacji, • Teorii podejmowania decyzji i teorii gier, • Teorii kolejek (masowej obsługi), • Teorii grafów, • Symulacji. 3
PROGRAMOWANIE LINIOWE 4
Wstęp do Programowania Liniowego (PL) • Model PL ma na celu poszukiwanie maksimum bądź minimum funkcji liniowej przy liniowych ograniczeniach • Elementy modelu PL: – Zbiór zmiennych decyzyjnych – Funkcja kryterium. – Układ ograniczeń. 5
Model Programowania Liniowego w postaci klasycznej 6
Model Programowania Liniowego w postaci klasycznej 7
Wstęp do PL Zastosowania modeli LP w różnych dziedzinach: • Produkcja • Finanse • Rolnictwo • Marketing i reklama, itd. . 8
Wstęp do PL • Istotna rola Programowania Liniowego – Efektywne algorytmy obliczeniowe gwarantujące znalezienie rozwiązania optymalnego – Możliwa analiza wrażliwości rozwiązania optymalnego – co by było, gdyby. . . ? . 9
Wstęp do PL • Założenia modelu PL: – Znane wartości parametrów, – Funkcja kryterium i ograniczenia mają własność stałych przyrostów (constant returns to scale) – ten sam co do wielkości przyrost zmiennej , bez względu na początkowy poziom, powoduje zawsze taki sam przyrost wartości funkcji – Addytywność efektów związanych ze zmiennymi, – Zmienne decyzyjne mają charakter ciągły – mogą przyjąć każdą wartość z określonego przedziału liczbowego (inne modelowanie dla zmiennych całkowitoliczbowych czy też binarnych), – Zakłada się nieujemność zmiennych decyzyjnych. 10
Firma „Puchatek” – problem optymalnego planu produkcji • Firma produkuje dwa rodzaje zabawek plastikowych - samochodzików - dla dzieci powyżej 1 roku: – ciężarówka. – traktor. • Występują ograniczone zasoby dwóch środków produkcji: – 1000 kg specjalnego plastiku. – Czas produkcji w ciągu tygodnia 11
Firma „Puchatek” – problem optymalnego planu produkcji • Wymagania rynkowe – Wielkość produkcji nie może przekroczyć 7000 szt. – Liczba ciężarówek nie może przekroczyć liczby traktorów o więcej niż 3500 szt. • Informacja technologiczna – Ciężarówka wymaga 20 dkg plastiku i 0, 3 minut czasu produkcji, – Traktor wymaga 10 dkg plastiku i 0, 4 minut 12 czasu pracy.
Firma „Puchatek” – problem optymalnego planu produkcji • Obecna strategia planowania produkcji: – Produkować jak najwięcej produktu bardziej zyskownego (Ciężarówka – zysk jedn. 8 zł za dziesięć sztuk), – Pozostałe środki przeznaczyć na produkt mniej zyskowny (Traktor – zysk jedn. 5 zł za dziesięć sztuk), pamiętając o zaleceniach działu marketingu. • Obecny tygodniowy plan produkcji: Ciężarówka Traktor Szacowany zysk 8(450) + 5(100) = 4500 sztuk = 1000 sztuk = 4100 zł tygodniowo 13
Firma szuka rozwiązania, które może przynieść zwiększenie zysku 14
Model PL dla firmy „Puchatek” • Zmienne decyzyjne: – X 1 = tygodniowa wielkość produkcji ciężarówek (w 10 szt. ) – X 2 = tygodniowa wielkość produkcji traktorów (w 10 szt. ) • Funkcja kryterium: – maksymalizacja zysku tygodniowego 15
Model PL dla firmy „Puchatek” Max z(x) = 8 X 1 + 5 X 2 zł) (zysk tygodniowy w przy ograniczeniach: 2 X 1 + 1 X 2 <= 1000 3 X 1 + 4 X 2 <= 2400 (plastik w kg) (czas produkcji w minutach) X 1 + X 2 <= 700 (wielkość produkcji w 10 szt. ) X 1 - X 2 <= 350 (Mix) Xj 0, j = 1, 2 (nieujemność zmiennych 16
Analiza graficzna zadania PL Zbiór punktów, które spełniają wszystkie ograniczenia to ZBIÓR ROZWIĄZAŃ DOPUSZCZALNYCH 17
Analiza graficzna – zbiór rozwiązań dopuszczalnych X 2 Ograniczenia na nieujemność zmiennych X 1 18
Analiza graficzna – zbiór rozwiązań dopuszczalnych X 2 1000 700 500 Plastik 2 X 1+X 2 <= 1000 Produkcja całkowita: X 1+X 2 <=700 (nieistotne) Niedopuszczalne Czas produkcji Dopuszczalne 3 X 1+4 X 2 <=2400 500 700 X 1 19
Analiza graficzna – zbiór rozwiązań dopuszczalnych X 2 1000 700 500 Plastik 2 X 1+X 2 <= 1000 Produkcja całkowita: X 1+X 2 <=700 (nieistotne) Niedopuszczalne Mix X 1 -X 2 <= 350 Czas produkcji Dopuszczalne 3 X 1+4 X 2<= 2400 500 700 X 1 Punkty brzegowe Punkty wierzchołkowe Punkty wewnętrzne. • Trzy rodzaje rozwiązań 20
Poszukiwanie rozwiązania optymalnego 21
Poszukiwanie rozwiązania optymalnego X 2 1000 700 500 Ustalamy dowolną wielkość zysku, np. = 2000 zł, i rysujemy odpowiadającą izokwantę funkcji kryterium. (Izokwanta liniowej funkcji kryterium to prosta mająca tę własność, że dla wszystkich punktów tej prostej wartość funkcji jest =4360 jednakowa) Zysk zł Zwiększamy zysk tak dalece jak to mo. . . i kontynuujemy, dopóki jest to dopus X 1 500 22
Podsumowanie rozwiązania optymalnego Ciężarówki = 3200 szt. Traktory = 3600 szt. Zysk maksymalny = 4360 zł – Rozwiązanie optymalne wykorzystuje cały zasób surowca – plastik oraz czasu produkcji – ograniczenia wiążące. – Produkcja całkowita to 6800 szt. (a nie max 7000 szt. ) – Ograniczenie na Mix produktów spełnione jako nierówność: 320 - 360 = -40 < 350 23
Punkty wierzchołkowe a rozwiązanie optymalne – Jeżeli problem PL posiada rozwiązanie optymalne, to jest nim punkt wierzchołkowy, przynajmniej jeden. 24
Punkty wierzchołkowe a rozwiązanie optymalne Ø Jeżeli dokonany zostanie wybór rozwiązania optymalnego, to proste przecinające się w punkcie wierzchołkowym, będącym rozwiązaniem optymalnym, odpowiadają ograniczeniom wiążącym, tj. spełnionym jakoØrównania. W problemie firmy „Puchatek” ograniczeniami wiążącymi są: zapas plastiku oraz czas produkcji. Oznacza to, że cały zapas surowca jest wykorzystany. Również czas produkcji ØPozostałe ograniczenia są niewiążące – wykorzystany w 100% na wielkość obserwujemy zapas wjest ograniczeniu produkcji oraz mix produktów. Zapas – różnica między wartością prawej i lewej strony ograniczenia 25
Niejednoznaczne rozwiązanie optymalne • W przypadku niejednoznaczności rozwiązania optymalnego, izokwanta funkcji kryterium jest równoległa do jednego z ograniczeń. • W przypadku niejednoznacznosci każda liniowa kombinacja (średnia ważona) optymalnych rozwiązań wierzchołkowych jest również optymalna 26
Analiza wrażliwości rozwiązania optymalnego • Jak wrażliwe jest rozwiązanie optymalne na zmiany parametrów modelu? • Powody przeprowadzania analizy wrażliwości: – Założenie o znanych wartościach parametrów nie jest prawdziwe – znamy tylko wartości ocen statystycznych lub eksperckich parametrów – możliwy błąd szacunku, 27
Wrażliwość rozwiązania na zmiany parametrów funkcji kryterium. • Przedział optymalności – Rozwiązanie optymalne pozostaje niezmienne tak długo jak • Parametr funkcji kryterium należy do przedziału optymalności • Nie obserwujemy zmian innych parametrów modelu. – Wartość funkcji kryterium ulegnie zmianie, jeżeli analizowany parametr dotyczy zmiennej, której wartość jest 28
Wrażliwość rozwiązania na zmiany parametrów funkcji kryterium. X 2 1000 M 8 X 1 + 5 X 500 ax 2 M ax x 3 4 X. 75 1 + X 5 X 2 Max 2 X 1 + 5 X 2 X 1 500 800 29
Wrażliwość rozwiązania na zmiany parametrów funkcji kryterium. X 2 1000 1 8 X ax M + Przedział optymalności: [3. 75, 10] 5 X 3. 7 1 1 +5 X 2 5 X +5 Ma x 0 X x 1 Ma 2 500 X 2 400 600 800 X 1 30
Wrażliwość rozwiązania na zmiany parametrów funkcji kryterium. Interpretacja przedziału optymalności dla parametru c 1: Zakładając, że inne elementy modelu (parametry, ograniczenia) nie ulegną zmianie, to zmiana zysku jednostkowego (w 10 szt. ) dla ciężarówek w przedziale [3, 75 ; 10] zł nie spowoduje utraty optymalności przez uzyskane rozwiązanie. Maksymalny zysk odpowiada produkcji 3200 ciężarówek i 3600 szt. traktorów. Oczywiście, zmiana zysku jednostkowego dla ciężarówek spowoduje zmianę wartości maksymalnego zysku, np. dla c 1=9 zł/10 szt. maksymalny zysk wyniesie 320*9+360*5= 4680 zł. 31
Analiza wrażliwości rozwiązania na zmianę prawych stron ograniczeń • Jak zmieni się optymalna wartość funkcji kryterium (np. maksymalny zysk), jeżeli prawa strona wybranego ograniczenia wzrośnie o jednostkę? • Dla jak dużych przyrostów bądź spadków wartości prawej strony ograniczenia, wyznaczona wartość przyrostu funkcji kryterium pozostanie niezmieniona? 32
Analiza wrażliwości rozwiązania na zmianę prawych stron ograniczeń • Każda zmiana wartości prawej strony ograniczenia wiążącego spowoduje zmianę rozwiązania optymalnego. • Dowolna zmiana prawej strony ograniczenia niewiążącego, mniejsza od wielkości zapasu, nie spowoduje zmiany rozwiązania optymalnego, 33
Dualizm w programowaniu liniowym 34
Własności zadania dualnego: • Jeżeli jedno z pary zadań nie posiada skończonego rozwiązania optymalnego, to drugie z zadań jest sprzeczne, • Jeżeli jedno z pary zadań jest sprzeczne, to drugie może być sprzeczne bądź nie posiadać skończonego rozwiązania optymalnego, • Każda ze zmiennych dualnych odpowiada konkretnemu ograniczeniu zadania prymalnego, 35
Interpretacja wycen dualnych • interpretacja wynika z własności równości optymalnych wartości funkcji kryterium obu zadań: przyrost optymalnej wartości funkcji kryterium zadania prymalnego spowodowany marginalnym przyrostem prawej strony odpowiadającego ograniczenia (pamiętamy, że zmiana wartości prawej strony ograniczenia powoduje, w ogólnym przypadku, zmianę wartości zmiennych zadania PL). 36
Wyceny dualne (Shadow Prices) Zakładając, że nie występują zmiany żadnych innych parametrów wejściowych modelu, zmiana optymalnej (max albo min) wartości funkcji kryterium na jednostkę przyrostu wartości prawej strony ograniczenia nazywana jest wyceną (ceną) dualną (najczęściej, wyceną dualną zasobu) 37
Plastik Wyceny dualne – ilustracja graficzna X 2 1000 Max zysk = 4360 zł x 2 +1 01 10 <= 00 x 2 10 +1 <= 2 X 1 Jeżeli dostępna jest większa ilość plastiku (ograniczenie na zasób plastiku będzie rozluźnione), wzrasta wartość prawej strony ograniczenia 500 Max zysk = 4363. 4 zł Wycena dualna = 4363. 40 – 4360. 00 = 3. 40 Czas produkcji X 1 500 38
Wyceny dualne – interpretacja c. d. • Zmienna dualna posiada miano, wynikające ze sposobu pomiaru wartości funkcji kryterium i wartości ograniczenia, np. y 1=3, 4 $/kg (dla 1. ograniczenia na zasób plastiku) • Jeżeli zapas plastiku zwiększy się o 1 kg to maksymalny zysk (odpowiadający nowemu rozwiązaniu optymalnemu) zwiększy się o 3, 4 $ i wyniesie 4360 + 39
Własności zadania dualnego c. d. : • W przypadku modeli PL o mieszanych warunkach ograniczających, zmienne dualne odpowiadające ograniczeniom o przeciwnych znakach niż dla symetrycznej pary (max „ ” oraz min „ ”) są niedodatnie; w przypadku ograniczeń równościowych nie można przewidzieć znaku zmiennej dualnej. 40
Przedział dopuszczalności • Zakładając brak zmian wartości innych parametrów wejściowych modelu, przedziałem dopuszczalności nazywamy: Przedział wartości prawej strony ograniczenia, w zakresie którego nie ulegają zmianie wyceny dualne. • W obrębie przedziału dopuszczalności, zmianę optymalnej wartości funkcji kryterium możemy wyznaczyć następująco: Zmiana wartości f. kryterium = [wycena dualna]x[zmiana wartości prawej 41 strony ograniczenia]
Przedział dopuszczalności Plastik X 2 2 X 1 x 2 +1 1000 00 10 <= Produkcja całkowita X 1 + X 2≤ 700 Zwiększanie zasobu plastiku przynosi efekt tylko do czasu, aż pojawi się nowe ograniczenie Nowe ograniczenie wiążące 500 To jest rozwiązanie niedopuszczaln Czas produkcji X 1 500 42
Plastik Przedział dopuszczalności X 2 2 X 1 x 2 +1 1000 0 0 10 <= Zauważmy, jak zmienia się zysk, gdy rośnie zasób plastiku. 500 Czas produkcji X 1 500 43
Przedział dopuszczalności X 2 Zasób plastiku zmniejsza się (ograniczenie jest bardziej 1000 Rozwiązanie niedopuszczalne restrykcyjne). Zysk zmniejsza się 500 2 X 1 + 1 X 2 <= 1100 Nowe ograniczenie wiążące X 1 500 44
„Puchatek” – wprowadzanie danych w programie Win. QSB Wprowadzanie danych w programie Win. QSB 45
„Puchatek” – rozwiązanie graficzne w programie Win. QSB Zapas plastik u Ilość wyrobó w mix Czas pracy 46
„Puchatek” – rozwiązanie w programie Win. QSB Przedziały optymalności Wyceny dualne Przedziały dopuszczalno ści Zapas/nadmia r 47
• Możliwe, inne niż jednoznaczne, wyniki optymalizacji Sprzeczność zadania: Zbiór rozwiązań dopuszczalnych jest pusty. Powodem są zbyt restrykcyjne ograniczenia. • Nieograniczoność: Funkcja kryterium może być dowolnie duża. Powodem jest brak istotnego ograniczenia w modelu. • Rozwiązanie niejednoznaczne: Więcej niż jeden punkt odpowiada optymalnej wartości funkcji kryterium 48
Zadanie PL jest sprzeczne . 2 3 1 49
Rozwiązanie nieograniczone fu nk cji ro Zb do zwią iór h pus zań zc za ln y Ma kr yte ks ym riu ali m za c c 50 ja
- Slides: 50