|
发表于 2005-9-14 06:09:55
|
显示全部楼层
<>2004电力市场输电阻塞管理</P>
<>要准确计费,必须给出第i台机组出力和对应其当前段价的函数关系式,这是个难点,没几个队能正确写出,好在当年没有最佳答案,大家就可各显神通,错了评委也难知道.</P>
<>想用LINGO解,必用其自定义函数@user(),我们可定义其为阶跃函数,从而解决本问题.</P>
<P> </P>
<P>MODEL:<BR>SETS:<BR> i6/1..6/:y,c;<BR> i8/1..8/:x,x0,u,v,T,TT,q,vv,uu,xx;<BR> i12/1..12/:z6;<BR> i10/1..10/:;<BR> cc(i6,i8):A;<BR> dd(i8,i10):e,h,ee;<BR>ENDSETS<BR>p=0.5;<BR>!j0=303; !r=0;<BR>j0=356; r=0.397; <BR><a href="mailtmin=0.25*@sum(i8(i):p*j0*(-x(i))*@user(-x(i))+T(i)*@user(x(i" target="_blank" >min=0.25*@sum(i8(i):p*j0*(-x(i))*@user(-x(i))+T(i)*@user(x(i</A>)));</P>
<P>@for(i8(i):T(i)=x(i)*(TT(i)-j0));</P>
<P>@for(i8(i):TT(i)=h(i,1)*@user(x0(i)+x(i))+(h(i,2)-h(i,1))*@user(x0(i)+x(i)-ee(i,1))+(h(i,3)-h(i,2))*@user(x0(i)+x(i)-ee(i,2))+(h(i,4)-h(i,3))*@user(x0(i)+x(i)-ee(i,3))+(h(i,5)-h(i,4))*@user(x0(i)+x(i)-ee(i,4))+(h(i,6)-h(i,5))*@user(x0(i)+x(i)-ee(i,5))+(h(i,7)-h(i,6))*@user(x0(i)+x(i)-ee(i,6))+(h(i,8)-h(i,7))*@user(x0(i)+x(i)-ee(i,7))+(h(i,9)-h(i,8))*@user(x0(i)+x(i)-ee(i,8))+(h(i,10)-h(i,9))*@user(x0(i)+x(i)-ee(i,9)));</P>
<P>@sum(i8(i):x(i))=0;<BR>@sum(i8(j):A(1,j)*(x(j)+x0(j)))+c(1)<z6(1)*(1+z6(7)*r);<BR>@sum(i8(j):A(2,j)*(x(j)+x0(j)))+c(2)<z6(2)*(1+z6(8)*r);<BR><a href="mailt-@sum(i8(j):A(3,j)*(x(j)+x0(j)))+c(3)<z6(3)*(1+z6(9)*r" target="_blank" >-@sum(i8(j):A(3,j)*(x(j)+x0(j)))+c(3)<z6(3)*(1+z6(9)*r</A>);<BR>@sum(i8(j):A(4,j)*(x(j)+x0(j)))+c(4)<z6(4)*(1+z6(10)*r);<BR>@sum(i8(j):A(5,j)*(x(j)+x0(j)))+c(5)<z6(5)*(1+z6(11)*r);<BR>@sum(i8(j):A(6,j)*(x(j)+x0(j)))+c(6)<z6(6)*(1+z6(12)*r);</P>
<P>@for(i6(i):y(i)=@sum(i8(j):A(i,j)*(x0(j)+x(j)))+c(i));<BR>@for(i8(i):xx(i)=x0(i)+x(i));<BR>@FOR(i8(i)BND(vv,xx,uu));<BR>@FOR(i8(i)BND(-v,x,u));<BR>@for(i8(i)free(x(i)));<BR>@for(i8(i):@free(xx(i)));<BR>@for(i6(i):@free(y(i)));<BR>@for(i8(i):@free(q(i)));<BR>@for(i8(i):@free(TT(i)));<BR>@for(i8(i):@free(T(i)));</P>
<P>DATA:<BR>A=0.0828, 0.0483, 0.0530, 0.1199, -0.0254, 0.1220, 0.1216, -0.0012,<BR> -0.0546, 0.1279, 0, 0.0333, 0.0868, -0.1124, -0.0189, 0.0987,<BR> -0.0695, 0.0616, -0.1566, -0.0099, 0.1245, 0.0021, -0.0025, -0.2014,<BR> -0.0345, -0.1024, 0.2052, -0.0208, -0.0118, 0.0060, 0.1449, 0.0765,<BR> 0.0005, 0.2433, -0.0646, -0.0411, -0.0652, 0.0703, -0.0043, -0.0089,<BR> 0.2378, -0.0602, -0.0779, 0.0930, 0.0469, 0.0001, 0.1659, 0.0007;</P>
<P><BR>c=110.2965 131.2289 -108.8732 77.4817 132.9745 120.6633;</P>
<P><BR>!v=63.0000 21.0000 48.0000 39.0000 27.0000 45.0000 34.9000 50.9000;<BR>!u=3.0000 ,9.0000 ,48.0000 , 0 , 27.0000 , 15.0000 , 7.1000 , 3.1000;<BR>u=3.0000 7.0000 9.8000 0 17.0000 5.0000 0.0000 0;<BR>v= 63.0000 23.0000 86.2000 39.0000 37.0000 55.0000 42.0000 54.00;</P>
<P>z6=165, 150, 160, 155, 132, 162, 0.13, 0.18, 0.09, 0.11, 0.15, 0.14 ;<BR>!x0=150.0000 79.0000 180.0000 99.5000 125.0000 140.0000 95.0000 113.9000;<BR>x0= 150.0000 81.0000 218.2000 99.5000 135.0000 150.0000 102.100 117.0000;</P>
<P>uu=190 89 280 116 155 180 125 160;<BR>vv=0 0 0 0 0 0 0 0;<BR>!x1=153 86.87 228 90.1124 152 95.3222 60.1 117;</P>
<P>e=70 0 50 0 0 30 0 0 0 40<BR> 30 0 20 8 15 6 2 0 0 8<BR> 110 0 40 0 30 0 20 40 0 40<BR> 55 5 10 10 10 10 15 0 0 1<BR> 75 5 15 0 15 15 0 10 10 10<BR> 95 0 10 20 0 15 10 20 0 10<BR> 50 15 5 15 10 10 5 10 3 2<BR> 70 0 20 0 20 0 20 10 15 5;</P>
<P>h=-505 0 124 168 210 252 312 330 363 489<BR> -560 0 182 203 245 300 320 360 410 495<BR> -610 0 152 189 233 258 308 356 415 500<BR> -500 150 170 200 255 302 325 380 435 800<BR> -590 0 116 146 188 215 250 310 396 510<BR> -607 0 159 173 205 252 305 380 405 520<BR> -500 120 180 251 260 306 315 335 348 548<BR> -800 153 183 233 253 283 303 318 400 800;</P>
<P><BR> ee=70 70 120 120 120 150 150 150 150 190<BR> 30 30 50 58 73 79 81 81 81 89<BR> 110 110 150 150 180 180 200 240 240 280<BR> 55 60 70 80 90 100 115 115 115 116<BR> 75 80 95 95 110 125 125 135 145 155<BR> 95 95 105 125 125 140 150 170 170 180<BR> 50 65 70 85 95 105 110 120 123 125<BR> 70 70 90 90 110 110 130 140 155 160;</P>
<P><BR>ENDDATA<BR> END<BR></P> |
|