گمز 63565

---

شرح پروژه

#کد_پروژه: 63565

موضوع: پروژه گمز

میخواستم عرض کنم که بنده به برنامه ای رو نوشتم توی گمز،یه مشکل میگیره در برنامه و ران نمیتونم بگیرم امکانش هست بنده برنامه ای که نوشتم رو برایتون ارسال کنم و اینکه راهنمایی کنید کجای برنامه مشکل دارهخیلی ممنون بابت راهنمایی $kill

$clear

sets

t  /1*720/

j  /1*77/

*-------------------Defining Parameters----------------*

parameters

MaxOpt

Ta(t),q(t),Ev(t)

Tw

RRmax,MinD

a0,a1

b0,b1,b2,b3

SSmin,SSmax,Spmax

PPC

S0(j)

aa

;

*-------------------Calling Data----------------------*

ppc=520        ;

Tw=172.6       ;

RRmax=1.176043 ;

SSmax=2749.8   ;

SSmin=740      ;

a0=0.0201      ;

a1=15.275      ;

b0=0.000000003 ;

b1=0.00002     ;

b2=0.0626      ;

b3=269.59      ;

Spmax=9.36     ;

MinD=241.553334;

$call gdxxrw Hourly_Apr.xlsx par=S0 rng=S0! rdim=1 cdim=0

$gdxin Hourly_Apr.gdx

$load S0

$gdxin

$call gdxxrw Hourly_Apr.xlsx par=q rng=q! rdim=1 cdim=0

$gdxin Hourly_Apr.gdx

$load q

$gdxin

$call gdxxrw Hourly_Apr.xlsx par=Ev rng=Ev! rdim=1 cdim=0

$gdxin Hourly_Apr.gdx

$load Ev

$gdxin

$call gdxxrw Hourly_Apr.xlsx par=Ta rng=Ta! rdim=1 cdim=0

$gdxin Hourly_Apr.gdx

$load Ta

$gdxin

MaxOpt=226200;

aa=0.9994;

*-------------------Defining Variables----------------*

variables

Z

z1

z2(t)

z3

;

positive variables

Rm,Pm

Income(t),P(t),S(t),L(t),R(t),Sp(t),H(t),A(t)

H1,H2,Se,SumSp

;

binary variables

bi(t)

;

*--------------------Initial Condition----------------*

S.fx('1')=S0('1')

;

*-------------------Defining Equations----------------*

equations

obj

eq0,eq1,eq2,eq3

eq4,eq5,eq6,eq7

eq9,eq10,eq11,eq12

eq13,eq14,eq15,eq16,eq17,eq18,eq19

eq20,eq21

eq22,eq23,eq24,eq25

eq26,eq27,eq28

;

*-----------------------Equations---------------------*

obj .. Z=e=aa*sqr(1-z1/MaxOpt)+(1-aa)*z3;

*-----------------------------------------------------*

eq0..    z1=e=sum(t,Income(t));

eq1..    z3=e=sum(t,z2(t));

eq2(t).. z2(t)=e=sqr(1-720*(R(t)/(Rm+0.0000001)));;

eq3(t).. Income(t)=e=P(t)*Ta(t);

*-----------------------------------------------------*

eq4.. H1=e=b0*S('1')3-b1*S('1')2+b2*S('1')+b3;

eq5.. H2=e=b0*Se3-b1*Se2+b2*Se+b3;

eq6.. Rm=e=sum(t,q(t))+S('1')-Se-sum(t,L(t))-Sum(t,Sp(t));

eq7.. Pm=e=0.003785*[((H1+H2)/2)-Tw]*Rm*0.9;

*eq8.. SumSp=e=sum(t,Sp(t));

*-----------------------------------------------------*

eq9(t)$(ord(t) NE 720).. P(t)=e=0.003785*30*24*[((H(t)+H(t+1))/2)-Tw]*R(t)*0.9;

eq10(t)$(ord(t) EQ 720).. P(t)=e=0.003785*30*24*[((H(t)+H2)/2)-Tw]*R(t)*0.9;

eq11(t)$(ord(t) NE 720).. S(t+1)=e=S(t)+q(t)-R(t)-SP(t)-L(t);

eq12(t)$(ord(t) EQ 720).. S(t)+q(t)-R(t)-SP(t)-L(t)=e=Se;

*-----------------------------------------------------*

eq13(t).. S(t)=g=SSmin;

eq14(t).. S(t)=l=SSmax;

eq15(t).. Rm=g=MinD;

eq16(t).. R(t)=l=RRmax;

eq17(t).. P(t)=l=PPC;

eq18(t).. P(t)=l=1.6*Pm;

eq19(t).. P(t)=g=0.4*Pm;

*-----------------------------------------------------*

eq20.. Se=g=S0('77');

eq21.. Se=l=SSmax;

*-----------------------------------------------------*

eq22(t).. A(t)=e=a0*S(t)+a1;

eq23(t).. H(t)=e=b0*S(t)3-b1*S(t)2+b2*S(t)+b3;

eq24(t)$(ord(t) NE 720).. L(t)=e=(Ev(t)*(A(t)+A(t+1))/2)/1000;

eq25(t)$(ord(t) EQ 720).. L(t)=e=(Ev(t)*(A(t)+(a0*Se+a1))/2)/1000;

*-----------------------------------------------------*

eq26(t)$(ord(t) NE 720).. S(t+1)/SSmax=g=bi(t);

eq27(t)$(ord(t) EQ 720).. Se/SSmax=g=bi(t);

eq28(t).. SP(t)=l=Spmax*bi(t);

*-----------------------------------------------------*

model ccc /all/;

ccc.nodLim = 1e6;

ccc.resLim = 1e6;

Option MINLP = SBB;

Option IterLim = 1e9;

option optcr=0.00001;

solve ccc using MINLP minimizing Z;

*-----------------------------------------------------*

display Z.l,Z1.l,z3.l,z2.l,Rm.l,Pm.l,P.l,R.l,S.l,H.l,Sp.l,bi.l

;

*-----------------------------------------------------*

Execute_unload 'outdata.gdx';

Execute 'GDXXRW outdata.gdx var=S rng=Sheet1!a29';

Execute 'GDXXRW outdata.gdx var=H rng=Sheet2!a29';

Execute 'GDXXRW outdata.gdx var=R rng=Sheet3!a29';

Execute 'GDXXRW outdata.gdx var=sp rng=Sheet4!a29';

Execute 'GDXXRW outdata.gdx var=p rng=Sheet5!a29';

Execute 'GDXXRW outdata.gdx var=z1 rng=Sheet6!a29';

Execute 'GDXXRW outdata.gdx var=se rng=Sheet7!a29';  بنده مهندسی عمران هستم شاخه مهندسی و مدیریت منابع اب

ایتا کاظم انصاری

مهارت ها و تخصص های مورد نیاز

---

---