WUT_Computer_Science/Programming/WDWR/projekt/old_code/WDWR1.mod

139 lines
4.3 KiB
Modula-2
Raw Permalink Blame History

###########################################################
# WDWR 25406 #
# Planowanie produkcj w warunkach ryzyka. #
# MODEL #
# Autor: Krzysztof Rudnicki #
###########################################################
##########
# Zbiory #
##########
# Produkty
set PRODUCTS = {"P1", "P2", "P3", "P4"};
# Narzedzia
set TOOLS;
# Miesiace
set MONTHS ordered;
# Scenariusze
param scenarioCount = 10000;
set SCENARIOS = {0..scenarioCount};
#############
# Parametry #
#############
# Liczba kazdego z narzedzi
param toolCount {TOOLS} >= 1;
# Dochody ze sprzedazy [pln/szt]
param expectedProfitPerUnit {PRODUCTS} >= 0;
# Scenarios
param scenarioProfitPerUnit {SCENARIOS, PRODUCTS};
# Czasy produkcji [godz]
param toolTimePerUnit {TOOLS, PRODUCTS} >= 0;
# Ograniczenia rynkowe liczby sprzedawanych produktow [szt]
param salesMarketLimit {MONTHS, PRODUCTS} >= 0;
# Ograniczeine liczby magazynowanych produktow [szt]
param storageLimit {PRODUCTS} >= 0;
# Koszt magazynowania produktow [pln/szt per msc]
param storageUnitCost >= 0;
# Aktualny stan magazynowy [szt]
param startingStorage {PRODUCTS} >= 0;
# Pozadany stan magazynowy na koniec symulacji [szt]
param desiredEndStorage {PRODUCTS} >= 0;
# Liczba dni roboczych w miesiacu [d]
param daysPerMonth >= 1;
# Liczba zmian w ciagu jednego dnia roboczego
param shiftsPerDay >= 1;
# Dlugosc zmiany [godz]
param hoursPerShift >= 1;
# Liczba roboczogodzin w miesi<EFBFBD>cu [godz]
param workHoursPerMonth = daysPerMonth*shiftsPerDay*hoursPerShift;
###########
# Zmienne #
###########
# Produkcja produktow
var produced {MONTHS, PRODUCTS} >= 0 integer;
# Sprzedaz produktow w danym miesiacu
var sold {MONTHS,PRODUCTS} >= 0 integer;
var totalSold {p in PRODUCTS} = sum {m in MONTHS} sold[m, p];
# Sprzedaz z podziałem na normalną (do 80% limitu rynku) i z obniżonym zyskiem (powyżej 80% limitu)
var salesNormal {m in MONTHS, p in PRODUCTS} >= 0 integer;
var salesDiscounted {m in MONTHS, p in PRODUCTS} >= 0 integer;
# Iloosc produktow przekazanych do magazynu w danym miesiacu
var stored {m in MONTHS, p in PRODUCTS} = produced[m, p] - sold[m, p];
# Stan magazynowy na koniec danego miesiaca
var storage {m in MONTHS, p in PRODUCTS} =
startingStorage[p] + sum {m2 in MONTHS: ord(m2) <= ord(m)} stored[m2, p];
# Czas pracy narzedzi w danym miesi<EFBFBD>cu
param availableToolTime {t in TOOLS} =
toolCount[t]*workHoursPerMonth;
var usedToolTime {m in MONTHS, t in TOOLS} =
sum {p in PRODUCTS} produced[m,p]*toolTimePerUnit[t,p];
# Koszt magazynowania
var monthlyStorageCost {m in MONTHS} =
(sum {p in PRODUCTS} storage[m, p])*storageUnitCost;
var totalStorageCost = sum {m in MONTHS} monthlyStorageCost[m];
# Zysk dla warto<EFBFBD>ci oczekiwanej
var expectedSalesProfit =
sum {m in MONTHS, p in PRODUCTS} (
salesNormal[m, p] * expectedProfitPerUnit[p] +
salesDiscounted[m, p] * 0.8 * expectedProfitPerUnit[p]
);
var expectedNetProfit =
expectedSalesProfit - totalStorageCost;
#######################
# Ograniczenia modelu #
#######################
# Podział sprzedaży na normalną i z obniżonym zyskiem
subject to TotalSales {m in MONTHS, p in PRODUCTS}:
sold[m, p] = salesNormal[m, p] + salesDiscounted[m, p];
# Ograniczenie normalnej sprzedaży do 80% limitu rynkowego
subject to NormalSalesLimit {m in MONTHS, p in PRODUCTS}:
salesNormal[m, p] <= 0.8 * salesMarketLimit[m, p];
# Ograniczenie rynkowe sprzedazy produktow
subject to SalesMarketLimit {m in MONTHS, p in PRODUCTS}:
sold[m, p] <= salesMarketLimit[m, p];
# Ograniczenie magazynowe sprzedazy produktow
subject to SalesLimit1 {p in PRODUCTS}:
sold[first(MONTHS), p] <= produced[first(MONTHS), p];
subject to SalesLimit2 {m in MONTHS, p in PRODUCTS: m != first(MONTHS)}:
sold[m, p] <= produced [m, p] + storage[m, p];
# Ograniczenie pojemno<EFBFBD>ci magazynowej
subject to StorageLimit {m in MONTHS, p in PRODUCTS}:
storage[m, p] <= storageLimit[p];
# Ograniczenie na po<EFBFBD><EFBFBD>dany stan magazynowy na koniec marca
subject to DesiredStorage {p in PRODUCTS}:
storage[last(MONTHS), p] >= desiredEndStorage[p];
#Ograniczenie czasu pracy narzedzi w miesiacu
subject to ToolWorkTime {m in MONTHS, t in TOOLS}:
usedToolTime[m, t] <= availableToolTime[t];
################
# Funkcje celu #
################
maximize Profit: expectedNetProfit;