Witam, choć to może mało związane z takim prawdziwym programowaniem (Xpress
Mosel),
to na pewno siedzą tu ścisłe mózgi i może podpowiedzą mi, gdzie ten cholerny
błąd siedzi, bo
już z dwie godziny próbuję coś zmieniać w kodzie i nic nie wychodzi (ciągle
dostaje komunikat, że typy zmiennych są niekompatybilne):


model ModelName
uses "mmxprs";
parameters
!liczba oczyszczalni N
N=3
M=3
! liczba wariantów
v=6
end-parameters

declarations
oczyszczalnie=1..N
przekroje=1..M
warianty=1..v
! czy w danej oczyszczalni używamy danego wariantu?
y:array(oczyszczalnie, warianty) of mpvar
BZT:array(przekroje) of linctr
x:array(oczyszczalnie) of real
xwar:array(warianty) of real
q:array(oczyszczalnie) of real
Q:array(przekroje) of real
p:array(przekroje) of real
B:array(przekroje) of real
u:array(oczyszczalnie) of real
k:array(przekroje,przekroje) of real
kb:array(przekroje,przekroje) of real
c:array(oczyszczalnie) of real
cwar:array(warianty) of real
sumy:array(oczyszczalnie) of integer
t:array(przekroje,przekroje) of real
end-declarations
!
!
Q::[15, 18, 25]
q::[0.8, 0.7, 0.5]
p::[220, 250, 265]
B::[4, 4, 3]
xwar::[0.6, 0.7, 0.75, 0.8, 0.9, 0.95]
cwar::[10, 12, 13, 16, 22, 30]
!u to xmax
u::[0.85, 0.9, 0.9]
t::[0, 1, 3,
0, 0, 2,
0, 0, 0]
kb::[0, 0.25, 0.25,
0, 0, 0.27,
0, 0, 0]
forall (j in oczyszczalnie, z in warianty)
y(j,z)=0 or y(j,z)=1
!tylko w jednym miejscu dla danej oczyszczalnie jedynka, dlatego:
forall (j in oczyszczalnie) do
b(j):=(sum(z in warianty) y(j,z))
b(j)=1
end-do
j:=1
while j<=v do
x(j):=(sum(z in warianty) y(j,z)*xwar(z))
!wartości C
c(j):=(sum(z in warianty) y(j,z)*cwar(z))
end-do
forall (i in przekroje, j in przekroje)
k(i,j):=exp(-kb(i,j)*t(i,j))

!ograniczenia
forall (i in przekroje) do
BZT(i):=(sum(j in 1..i) (1-x(j))*p(j)*q(j)*k(j,i))/(Q(i)+sum(j in 1..i)
q(j))
BZT(i)<=B(i)
BZT(i)>=0
end-do
!ograniczenia dot. zakresu stopnia redukcji zanieczyszczen
forall (j in oczyszczalnie) do
x(j)>=0
x(j)<=u(j)
end-do
!Funkcja celu:
z:=sum(j in oczyszczalnie) c(j)*x(j)
minimize(z)

writeln("Niezbędne środki wynoszą: ", getobjval)
forall (i in przekroje) do
writeln(" Stopień redukcji w oczyszczalni nr (",i,") wynosi
",getsol(x(i)))
writeln(" BZT (",i,") wynosi ",getsol(BZT(i)))
end-do
end-model