참조
Alireza Soroudi, 슬롯 나라의 전력 시스템 최적화 모델링, 모델DED-PB(Gcode4.5) 장동적 경제 파견, 2017
메인파일 : DED-PB.gms
$title 가격 기반 동적 경제 부하 파견

$onText
자세한 내용은 다음 책의 4장(Gcode4.5)을 참조하세요.
소루디, 알리레자. 슬롯 나라의 전력 시스템 최적화 모델링. 스프링거, 2017.
--------------------------------------------------------------------------------
모델 유형: QCP
--------------------------------------------------------------------------------
작성자:
알리레자 소루디 박사
IEEE 수석 회원
이메일: alireza.soroudi@gmail.com
우리는 개발된 슬롯 나라 코드를 사용하여 파생된 출판물을 요청합니다.
인용함으로써 그 사실을 명시적으로 인정한다.
소루디, 알리레자. 슬롯 나라의 전력 시스템 최적화 모델링. 스프링거, 2017.
DOI: doi.org/10.1007/978-3-319-62350-4
$offText

세트
   t '시간' / t1*t24 /
   i '열 단위' / p1*p4 /;

스칼라 lim / inf /;

테이블 gendata(i,*) '발전기 비용 특성 및 한계'
       a b c d e f Pmin Pmax RU0 RD0
   p1 0.12 14.80 89 1.2 -5 3 28 200 40 40
   p2 0.17 16.57 83 2.3 -4.24 6.09 20 290 30 30
   p3 0.15 15.55 100 1.1 -2.15 5.69 30 190 30 30
   p4 0.19 16.21 70 1.1 -3.99 6.2 20 260 50 50;

테이블 데이터(t,*)
        람다 부하
   t1 32.71 510
   t2 34.72 530
   t3 32.71 516
   t4 32.74 510
   t5 32.96 515
   t6 34.93 544
   t7 44.9 646
   t8 52 686
   t9 53.03 741
   t10 47.26 734
   t11 44.07 748
   t12 38.63 760
   t13 39.91 754
   t14 39.45 700
   t15 41.14 686
   t16 39.23 720
   t17 52.12 714
   t18 40.85 761
   t19 41.2 727
   t20 41.15 714
   t21 45.76 618
   t22 45.59 584
   t23 45.56 578
   t24 34.72 544;

변수
   OF '목표(수익)'
   costThermal '열 단위 비용'
   p(i,t) '화력발전소에서 생산된 전력'
   EM '배출량 계산';

p.up(i,t) = gendata(i,"Pmax");
p.lo(i,t) = gendata(i,"Pmin");

방정식 Genconst3, Genconst4, costThermalcalc, Balance, EMcalc, EMlim, Benefitcalc;

costThermalcalc.. costThermal =e= sum((t,i), gendata(i,'a')*power(p(i,t),2)
                                           + gendata(i,'b')*p(i,t) + gendata(i,'c'));

Genconst3(i,t).. p(i,t+1) - p(i,t) =l= gendata(i,'RU0');

Genconst4(i,t).. p(i,t-1) - p(i,t) =l= gendata(i,'RD0');

Balance(t).. sum(i,p(i,t)) =l= data(t,'load');

EMcalc.. EM =e= sum((t,i), gendata(i,'d')*power(p(i,t),2)
                                  + gendata(i,'e')*p(i,t) + gendata(i,'f'));

EMlim..EM =l= lim;

이익 계산.. OF =e= sum((i,t), 1*data(t,'lambda')*p(i,t)) - costThermal;

모델 DEDPB / 모두 /;
qcp 최대화를 사용하여 DEDPB를 해결합니다.

임베디드 코드 연결:
- 슬롯 나라리더:
    기호: [  이름: p  ]
- 프로젝션:
    이름: p.l(i,t)
    새로운 이름: p_l(i,t)
- 엑셀작성기:
    파일: DEDPB.xlsx
    기호: [  이름: p_l, 범위: 보고서!A1  ]
endEmbeddedCode