B题论文(武汉大学杨双簧)(下)

baima908

< >六、实例求解结果<B><p></p></B></P>
< ><FONT face="Times New Roman">1</FONT>、对问题一的解答<p></p></P>
< ><FONT face="Times New Roman">(1)</FONT>、模型①<FONT face="Times New Roman"> (</FONT>带常数项<FONT face="Times New Roman">) </FONT>的线性系数及其检验结果<p></p></P>
<P >表<FONT face="Times New Roman">1(略)  </FONT></P>
<P ><FONT face="Times New Roman">(2)</FONT>、模型②<FONT face="Times New Roman"> (</FONT>不带常数项<FONT face="Times New Roman">) </FONT>的线性系数及其检验结果<p></p></P>
<P >表<FONT face="Times New Roman">2(略)</FONT><FONT face="Times New Roman">2</FONT>、对问题二的解答<p></p></P>
<P >本文设计了两种阻塞费用计算规则，并对其优劣做了详细讨论。具体细节见模型建立部分（<FONT face="Times New Roman">page7</FONT>）。<p></p></P>
<P ><FONT face="Times New Roman">3</FONT>、对问题三的解答<p></p></P>
<P ><FONT face="Times New Roman">   </FONT>预报负荷需求为<FONT face="Times New Roman">982.4</FONT><v:shape><FONT face="Times New Roman"> <v:imagedata></v:imagedata></FONT></v:shape>时的分配预案如表<FONT face="Times New Roman">3</FONT>。<p></p></P>
<P >表<FONT face="Times New Roman">3                       </FONT>分配预案<p></p></P>
<DIV align=center>
<TABLE  cellSpacing=0 cellPadding=0 border=1>

<TR >
<TD  width=168>
<P  align=center>机组号</P></TD>
<TD  width=36>
<P  align=center><FONT face="Times New Roman">1</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">2</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">3</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">4</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">5</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">6</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">7</FONT></P></TD>
<TD  width=45>
<P  align=center><FONT face="Times New Roman">8</FONT></P></TD></TR>
<TR >
<TD  vAlign=top width=168>
<P  align=center>初始出力（兆瓦）</P></TD>
<TD  width=36>
<P  align=center><FONT face="Times New Roman">120</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">73</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">180</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">80</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">125</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">125</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">81.1</FONT></P></TD>
<TD  width=45>
<P  align=center><FONT face="Times New Roman">90</FONT></P></TD></TR>
<TR >
<TD  vAlign=top width=168>
<P  align=center>出力分配预案（兆瓦）</P></TD>
<TD  width=36>
<P  align=center><FONT face="Times New Roman">150</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">79</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">180</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">99.5</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">125</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">140</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">95</FONT></P></TD>
<TD  width=45>
<P  align=center><FONT face="Times New Roman">113.9</FONT></P></TD></TR>
<TR >
<TD  vAlign=top width=168>
<P  align=center>报价（元）</P></TD>
<TD  width=36>
<P  align=center><FONT face="Times New Roman">252</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">300</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">233</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">302</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">215</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">252</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">260</FONT></P></TD>
<TD  width=45>
<P  align=center><FONT face="Times New Roman">303</FONT></P></TD></TR>
<TR >
<TD  vAlign=top width=168>
<P  align=center>允许爬坡值（兆瓦）</P></TD>
<TD  width=36>
<P  align=center><FONT face="Times New Roman">33</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">15</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">48</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">19.5</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">27</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">30</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">21</FONT></P></TD>
<TD  width=45>
<P  align=center><FONT face="Times New Roman">27</FONT></P></TD></TR>
<TR >
<TD  vAlign=top width=168>
<P >实际爬坡值（兆瓦）</P></TD>
<TD  width=36>
<P  align=center><FONT face="Times New Roman">30</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">6</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">0</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">19.5</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">0</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">15</FONT></P></TD>
<TD  width=48>
<P  align=center><FONT face="Times New Roman">13.9</FONT></P></TD>
<TD  width=45>
<P  align=center><FONT face="Times New Roman">23.9</FONT></P></TD></TR>
<TR >
<TD  vAlign=top width=168>
<P  align=center>清算价（元）</P></TD>
<TD  width=369 colSpan=8>
<P ><FONT face="Times New Roman">303</FONT></P></TD></TR></TABLE></DIV>
<P ><p><FONT face="Times New Roman"> </FONT></p></P>
<P ><p><FONT face="Times New Roman"> </FONT></p></P>
<P ><FONT face="Times New Roman">4</FONT>、对问题四的解答<p></p></P>
<P >预报负荷需求为<FONT face="Times New Roman">982.4</FONT><v:shape><FONT face="Times New Roman"> <v:imagedata></v:imagedata></FONT></v:shape>时的分配预案（表<FONT face="Times New Roman">3</FONT>），其阻塞检查结果及调整后的最终出力分配方案分别如表<FONT face="Times New Roman">4</FONT>，表<FONT face="Times New Roman">5</FONT>所示。<p></p></P>
<P >表<FONT face="Times New Roman">4                       </FONT>阻塞检验结果<p></p></P>
<P><FONT face="Times New Roman">biao5          </FONT>调整后的方案<p></p></P>
<P>
<TABLE  cellSpacing=0 cellPadding=0 border=1>

<TR >
<TD >
<P  align=center>机组出力（兆瓦）<p></p></P></TD>
<TD  vAlign=top>
<P ><FONT face="Times New Roman">[153,86.87, 228, 90.1124, 152, 95.3222, 60.1,117]<p></p></FONT></P></TD></TR>
<TR >
<TD >
<P  align=center>允许爬坡值（兆瓦）<p></p></P></TD>
<TD  vAlign=top>
<P ><FONT face="Times New Roman">[33,15,48,19.5,27,30,21,27]<p></p></FONT></P></TD></TR>
<TR >
<TD >
<P  align=center>实际爬坡值（兆瓦）<p></p></P></TD>
<TD  vAlign=top>
<P ><FONT face="Times New Roman">[33,15,48,10.11,27,-29.69,-21,27]<p></p></FONT></P></TD></TR>
<TR >
<TD >
<P  align=center>潮流限值（兆瓦）<p></p></P></TD>
<TD  vAlign=top>
<P ><FONT face="Times New Roman">[165,150,160,155,132,162]<p></p></FONT></P></TD></TR>
<TR >
<TD >
<P  align=center>实际潮流分布（兆瓦）<p></p></P></TD>
<TD  vAlign=top>
<P ><FONT face="Times New Roman">[165,150,155.26,124.51,131.51,159.53]<p></p></FONT></P></TD></TR>
<TR >
<TD >
<P  align=center>清算价（元）<p></p></P></TD>
<TD  vAlign=top>
<P ><FONT face="Times New Roman">495<p></p></FONT></P></TD></TR>
<TR >
<TD >
<P  align=center>阻塞费用<FONT face="Times New Roman">1</FONT>（元）<p></p></P></TD>
<TD  vAlign=top>
<P ><FONT face="Times New Roman">3183.1<p></p></FONT></P></TD></TR>
<TR >
<TD >
<P  align=center>阻塞费用<FONT face="Times New Roman">2</FONT>（元）<p></p></P></TD>
<TD  vAlign=top>
<P ><FONT face="Times New Roman">47155<p></p></FONT></P></TD></TR></TABLE></P>
<P >其中<FONT size=3>阻塞费用<FONT face="Times New Roman">1</FONT>、<FONT face="Times New Roman">2</FONT>分别是用阻塞费用规则一、二求得的。<p></p></FONT></P>
<P ><p><FONT face="Times New Roman"> </FONT></p></P>
<P ><FONT face="Times New Roman">5</FONT>、问题五的求解结果<p></p></P>
<P ><FONT face="Times New Roman">   </FONT>预报负荷需求为<FONT face="Times New Roman">1052.8</FONT><v:shape><FONT face="Times New Roman"> <v:imagedata></v:imagedata></FONT></v:shape>时分配预案、阻塞检验结果及最终出力分配方案分别如表<FONT face="Times New Roman">6</FONT>，表<FONT face="Times New Roman">7</FONT>，表<FONT face="Times New Roman">8</FONT>所示。<p></p></P>
<P >表<FONT face="Times New Roman">6                        </FONT>分配预案</P>
<P >表<FONT face="Times New Roman">7                         </FONT>阻塞检验</P>
<P >表<FONT face="Times New Roman">8               </FONT>调整后的出力分配方案方案<p></p></P>
<P>
<TABLE  cellSpacing=0 cellPadding=0 border=1>

<TR >
<TD >
<P  align=center>出力方案<FONT face="Times New Roman">(</FONT>兆瓦）<p></p></P></TD>
<TD  vAlign=top>
<P ><FONT face="Times New Roman">[153,88,228,99.5,152,155,60.3,117]<p></p></FONT></P></TD></TR>
<TR >
<TD >
<P  align=center>最高报价<FONT face="Times New Roman">(</FONT>元<FONT face="Times New Roman">)<p></p></FONT></P></TD>
<TD  vAlign=top>
<P ><FONT face="Times New Roman">[489,495,356,302,510,380,120,303]<p></p></FONT></P></TD></TR>
<TR >
<TD >
<P  align=center>允许爬坡值<FONT face="Times New Roman">(</FONT>兆瓦<FONT face="Times New Roman">)<p></p></FONT></P></TD>
<TD  vAlign=top>
<P ><FONT face="Times New Roman">[33,15,48,19.5,27,30,21,27]<p></p></FONT></P></TD></TR>
<TR >
<TD >
<P  align=center>实际爬坡值<FONT face="Times New Roman">(</FONT>兆瓦<FONT face="Times New Roman">)<p></p></FONT></P></TD>
<TD  vAlign=top>
<P ><FONT face="Times New Roman">[33,15,48,19.5,27,30,20.8,27]<p></p></FONT></P></TD></TR>
<TR >
<TD >
<P  align=center>潮流限值<FONT face="Times New Roman">(</FONT>兆瓦<FONT face="Times New Roman">)<p></p></FONT></P></TD>
<TD  vAlign=top>
<P ><FONT face="Times New Roman">[165,150,160,155,132,162]<p></p></FONT></P></TD></TR>
<TR >
<TD >
<P  align=center>实际潮流分布<FONT face="Times New Roman">(</FONT>兆瓦<FONT face="Times New Roman">)<p></p></FONT></P></TD>
<TD  vAlign=top>
<P ><FONT face="Times New Roman">[173.4093,143.5833,155.2113,124.6828,135.2969,160.4221]<p></p></FONT></P></TD></TR>
<TR >
<TD >
<P  align=center>相对安全裕度利用率<p></p></P></TD>
<TD  vAlign=top>
<P ><FONT face="Times New Roman">[39.2%,12.3%,4.25%,-85.16%,30.23%,19.29%]<p></p></FONT></P></TD></TR>
<TR >
<TD >
<P  align=center>清算价<FONT face="Times New Roman">(</FONT>元<FONT face="Times New Roman">)<p></p></FONT></P></TD>
<TD  vAlign=top>
<P ><FONT face="Times New Roman">510<p></p></FONT></P></TD></TR>
<TR >
<TD >
<P  align=center>阻塞费用<FONT face="Times New Roman">1(</FONT>元<FONT face="Times New Roman">)<p></p></FONT></P></TD>
<TD  vAlign=top>
<P ><FONT face="Times New Roman">1962.3<p></p></FONT></P></TD></TR>
<TR >
<TD >
<P  align=center>阻塞费用<FONT face="Times New Roman">2(</FONT>元<FONT face="Times New Roman">)<p></p></FONT></P></TD>
<TD  vAlign=top>
<P ><FONT face="Times New Roman">40533<p></p></FONT></P></TD></TR></TABLE></P>
<P>七、总结评价</P>
<P><FONT face="Times New Roman">优点(略)   </FONT></P>
<P >主要缺点表现在：因专业知识匮乏，没有结合现行的几种典型的电力市场运营模式<FONT face="Times New Roman">[6]</FONT>（如<FONT face="Times New Roman">Pool</FONT>模式、双边交易模式等）的特点给出更合理的阻塞管理办法；<p></p></P>
<P >参考文献：<p></p></P>

dianslm

<b>B题求解-----最终定论，大家可对比！
</b><>B题求解-----最终定论，大家可对比！</P><>1、回归---局部精度高是关键</P><>2、阻塞费用：多种算法，结果较多，但数量级不应差太远。</P><P>3、排表4序，加表3量&gt;w=982.4，....</P><P>4、w=982.4时：  5300元/15分钟（等价于21200元/小时），303，0*裕度向量，即不用裕度</P><P>5、w=1052.8时：1555元/15分钟（等价于6000元/小时），356，0.4*裕度向量（安全最好），即要用裕度，
      若用0.55*裕度向量或以上，则阻塞费用&lt;200元，
           若用0.58*裕度向量或以上，则阻塞费用=0元，
           阻塞费用随裕度向量的使用比例增大而减小，可到0，
      阻塞费用的最大值在1555元/15分钟左右，（356，0.4*裕度向量，此时安全度最好）</P><P>（标准的双目标最优化，安全优先，LINGO解较易，matlab解很难---因初始可行点在15分钟内很难找到，即使找到也常常不具一般性，95%以上的队在此编程是错的，当然不少队可能查觉不到自己的错误，而知此难度的队因此不得不使用遮羞布了----即伪造结果数踞，只要能哄过评委就行，实际上大多数评委也不知第4、5问的结果的数量级，这更进一步说明写好论文成为关键，正确计算成为次要！实际上，有不少模型和结果是错误的队冲上了评全国一等奖的序列。）</P><P>阻塞费用-----预定利润和实际利润的差（不要求精确，大体公平就行）
---------------------------------------------------------------------------------
1------    Y=AX+b
其中b=[110.2029  131.3297 -108.7072   77.5706  132.9927  120.4706]’;</P><P>A=  [ 0.0830   -0.0545   -0.0695   -0.0344    0.0005    0.2378
            0.0487    0.1259    0.0602   -0.1032    0.2444   -0.0592
            0.0531    0.0000   -0.1566    0.2052   -0.0646   -0.0776
            0.1200    0.0333   -0.0099   -0.0208   -0.0412    0.0931
            -0.0252    0.0870    0.1239   -0.0120   -0.0662    0.0472
              19   -0.1124    0.0022    0.0060    0.0703    0.0003
           0.1212   -0.0191   -0.0026    0.1448   -0.0041    0.1653
           -0.0010  0.0988   -0.2013    0.0764   -0.0085    0.0012]’;</P><P>3-------当预报负荷为W=982.4 MW时找到下一时段八个机组的出力为[150.0000   79.0000  180.0000   99.5000  125.0000  140.0000   95.0000  113.9000]；清算价303
当预报负荷W=1052.8 MW时找到下一时段八个机组的出力为[150.0000  81.0000  218.2000  99.5000  135.0000  150.0000  102.100 117.0000],  清算价356
</P>

dianslm

<>min R    （R为裕量的统一升缩因子，在0~1中变化）</P><>min F      （F堵塞费用）</P><>sum（dx(i)）=0         dx（i）为x（i）的增量</P><P>0&lt;=x（i）&lt;=x（i）的上界（即第i台机组最大发电量）</P><P>|y（i）|&lt;y（i）的上限*（1+R*裕度T（i））  （6个不等式）      </P><P>vlb&lt;dx&lt;ulb     ulb（i）----爬坡对dx（i）的上限限制</P><P>------------------------------------------------------------------------------------------------------------
上为求解最精巧的模型，其中：</P><P>sum（x（i））=1052.8      </P><P>x（i）=x0（i）+dx（i） </P><P>  x=（x（1） x（2）...x（8））  </P><P><a href="mailt用Lingo8.0解，须会使用@USE()函数，即定义阶跃函数，此外，还须懂一些特定的小计巧才行" target="_blank" >用Lingo8.0解，须会使用@USE()函数，即定义阶跃函数，此外，还须懂一些特定的小计巧才行</A>。</P>

dianslm

<>“matlab解很难---因初始可行点在15分钟内很难找到，即使找到也常常不具一般性”--------当然若用两阶段的思想，即先用很简单的函数来取代原复杂的目标函数，求初始可行解，则可过此关。</P><>       但接下来的求解还有不少难点，95%以上的队编程是错的，当然不少队可能查觉不到自己的错误，而知此难度的队因此不得不使用遮羞布了----即伪造结果数踞，只要能哄过评委就行，实际上大多数评委也不知第4、5问的结果的数量级（4问最小费用=5300元/15分钟；5问费用=1500），这更进一步说明写好论文成为关键，正确计算成为次要！实际上，有不少模型和结果是错误的队冲上了评全国一等奖的序列。</P>

cyh

<>湖北区成绩如何查到？</P>

future_fly

能否将你的完整论文发给我，以便交流，邮件：<a href="mailtfuture_fly@163.com" target="_blank" >future_fly@163.com</A> 万分感谢！！！

飞云0616

<>我们的论文也是上全国了，结果还没出来</P>
<>呵呵</P>
<>可以先交流一下么？</P>

sunzx

<DIV class=quote><B>以下是引用<I>dianslm</I>在2004-10-28 11:53:10的发言：</B>

<>“matlab解很难---因初始可行点在15分钟内很难找到，即使找到也常常不具一般性”--------当然若用两阶段的思想，即先用很简单的函数来取代原复杂的目标函数，求初始可行解，则可过此关。</P>
<>       但接下来的求解还有不少难点，95%以上的队编程是错的，当然不少队可能查觉不到自己的错误，而知此难度的队因此不得不使用遮羞布了----即伪造结果数踞，只要能哄过评委就行，实际上大多数评委也不知第4、5问的结果的数量级（4问最小费用=5300元/15分钟；5问费用=1500），这更进一步说明写好论文成为关键，正确计算成为次要！实际上，有不少模型和结果是错误的队冲上了评全国一等奖的序列。</P></DIV>
<>hehe,你这个量级是对的
<P>不过数值我算出来第一个是5100左右,第二个是1200左右. 这个也许也是个局部最优.
<P>不知道还有没有更优的解.</P>

whuxiaoxiao

<>楼主好象没有获奖啊</P><>哈哈</P><>在这里瞎卖</P>[em06][em06][em06]

dianslm

<>

1. &lt;SCRIPT language=JavaScript&gt;
   
2. var actioninfo3='单帖屏蔽';document.write (dvbbs_show_topic('65803','3','sunzx','','','sunzx@pku.edu.cn','||||||||||||||||||','','Images/userface/image1.gif','','','','☆nudter☆','<DIV class=quote><B>以下是引用<I>dianslm</I>在2004-10-28 11:53:10的发言：</B><P>“matlab解很难---因初始可行点在15分钟内很难找到，即使找到也常常不具一般性”--------当然若用两阶段的思想，即先用很简单的函数来取代原复杂的目标函数，求初始可行解，则可过此关。</P><P>       但接下来的求解还有不少难点，95%以上的队编程是错的，当然不少队可能查觉不到自己的错误，而知此难度的队因此不得不使用遮羞布了----即伪造结果数踞，只要能哄过评委就行，实际上大多数评委也不知第4、5问的结果的数量级（4问最小费用=5300元/15分钟；5问费用=1500），这更进一步说明写好论文成为关键，正确计算成为次要！实际上，有不少模型和结果是错误的队冲上了评全国一等奖的序列。</P></DIV><P>hehe,你这个量级是对的
   
3. <P>不过数值我算出来第一个是5100左右,第二个是1200左右. 这个也许也是个局部最优.
   
4. <P>不知道还有没有更优的解.</P>','2004-10-30 0:37:36','face1.gif','*.*.*.*','8957',1,0,37087,3,'新手上路','32','32','1','2004-10-30','1002','201','31','','1',0,'0','level0.gif',4,'2004-10-30 2:18:45','0',7,'tablebody2','0','0'));
   
5. &lt;script&gt;

复制代码

sunzx， 不过数值我算出来第一个是5100左右------若你真按评分标准定义的费用函数算，那决不会比5300低，比5300略大，比如5600那也不错，若比5300低，则你的计算还是有些问题！</P><>第二个是1200左右. 这个也许也是个局部最优.-----的确如此，不过从1000到1500都算不错！而低于1000（费用太看重），或高于1700（一定是算错）那都有问题！
</P><!--帖子浏览主体--><TABLE fixed; WORD-BREAK: break-all" cellSpacing=1 cellPadding=5 align=center><TR><TD class=tablebody2 vAlign=top width=175><TABLE cellSpacing=0 cellPadding=4 width="100%"><TR><TD glow(color=#9898BA,strength=2)" vAlign=center width=*></TD></TR></TABLE></TD></TR></TABLE>