发明名称 居民小区电动汽车充放电控制方法
摘要 本发明公开了一种居民小区电动汽车充放电控制方法,它采用双层离散粒子群算法,通过底层粒子群算法求解各电动汽车满足所有约束条件的充放电计划,再利用顶层粒子群算法优化得出居民小区所有电动汽车充放电功率的控制方法。该方法不但能够有效改善居民小区的负荷特性指标,提高居民小区的电网设备利用率,同时,还能显著降低电动汽车的充电费用,易于实施和推广。
申请公布号 CN105068419A 申请公布日期 2015.11.18
申请号 CN201510471445.2 申请日期 2015.08.04
申请人 西南交通大学 发明人 杨健维;黄宇;王湘;何正友
分类号 G05B13/02(2006.01)I 主分类号 G05B13/02(2006.01)I
代理机构 成都博通专利事务所 51208 代理人 陈树明
主权项 一种居民小区电动汽车充放电控制方法,其步骤为:A、记录电动汽车i的动力电池容量E<sub>i</sub>、最大充电功率P<sub>max,i</sub>、动力电池可用容量比例k<sub>i</sub>;并将当日电动汽车i的最初出行时刻、最后返回时刻、日行驶里程,分别设定为次日电动汽车i的最初出行时刻t<sub>i,s</sub>、最后返回时刻t<sub>i,e</sub>、日耗电量S<sub>i</sub>;其中i代表电动汽车的编号,i=1,2,3…I;I为居民小区内电动汽车总数;B、根据小区的历史基础电力负荷数据,预测出次日第j小时(j=1,2,3…24)内的基础负荷L<sub>j</sub>,进而得到次日基础负荷的日负荷率F<sub>2min</sub>、次日基础负荷的峰谷差F<sub>1max</sub>、次日基础负荷的波动均方差F<sub>3max</sub>;C、计算次日车主的最大总费用F<sub>4max</sub>将电动汽车i在次日的最后返回时刻t<sub>i,e</sub>,即以最大充电功率P<sub>max,i</sub>充电直至充满,所产生的总费用定义为次日车主的最大总费用F<sub>4max</sub>:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>F</mi><mrow><mn>4</mn><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>I</mi></munderover><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><msub><mi>t</mi><mrow><mi>i</mi><mo>,</mo><mi>e</mi></mrow></msub><mo>+</mo><mn>1</mn></mrow><mrow><msub><mi>t</mi><mrow><mi>i</mi><mo>,</mo><mi>e</mi></mrow></msub><mo>+</mo><msub><mi>t</mi><mrow><mi>i</mi><mo>,</mo><mi>c</mi></mrow></msub></mrow></munderover><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>&CenterDot;</mo><msub><mi>cp</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000774451460000011.GIF" wi="1139" he="156" /></maths>其中:cp<sub>j</sub>为第j小时的充电电价;t<sub>i,c</sub>为电动汽车i在最大充电功率P<sub>max,i</sub>下,将动力电池一次性充满所需的充电小时数,可由下式算出:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>t</mi><mrow><mi>i</mi><mo>,</mo><mi>c</mi></mrow></msub><mo>=</mo><mo>&lsqb;</mo><mfrac><msub><mi>S</mi><mi>i</mi></msub><msub><mi>P</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi><mo>,</mo><mi>i</mi></mrow></msub></mfrac><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000774451460000012.GIF" wi="998" he="147" /></maths>其中:[]为四舍五入取整;D、充放电功率约束条件和目标函数的确定设电动汽车i在次日第j小时内的分时充放电功率为p<sub>i,j</sub>,则其满足约束条件式(3)‑(4):<img file="FDA0000774451460000013.GIF" wi="1362" he="178" />其中,H为高电价时段集合,Z为整数集;<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mo>{</mo><mrow><mtable><mtr><mtd><mrow><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mn>1</mn></munderover><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>+</mo><msub><mi>S</mi><mi>i</mi></msub><mo>&lt;</mo><msub><mi>k</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msub><mi>E</mi><mi>i</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mn>2</mn></munderover><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>+</mo><msub><mi>S</mi><mi>i</mi></msub><mo>&lt;</mo><msub><mi>k</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msub><mi>E</mi><mi>i</mi></msub></mrow></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mrow><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mn>24</mn></munderover><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>+</mo><msub><mi>S</mi><mi>i</mi></msub><mo>&lt;</mo><msub><mi>k</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msub><mi>E</mi><mi>i</mi></msub></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></mrow>]]></math><img file="FDA0000774451460000021.GIF" wi="1130" he="550" /></maths>式(3)的含义为:只允许电动汽车i在高电价时段放电;式(4)的含义为:任意时刻电动汽车i的净放电量与日耗电量之和不能超过电动汽车i的可用容量;目标函数为:电动汽车的净充电量与日耗电量的偏差fitness:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>f</mi><mi>i</mi><mi>t</mi><mi>n</mi><mi>e</mi><mi>s</mi><mi>s</mi><mo>=</mo><mrow><mo>|</mo><mrow><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mn>24</mn></munderover><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>-</mo><mrow><mo>&lsqb;</mo><msub><mi>S</mi><mi>i</mi></msub><mo>&rsqb;</mo></mrow></mrow><mo>|</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000774451460000022.GIF" wi="1107" he="167" /></maths>E、以式(5)值最小为目标函数,以式(3)‑(4)为约束条件,通过粒子群算法迭代计算出次日电动汽车i在第j小时的充放电功率p<sub>i,j</sub>的迭代值P<sub>i,j</sub>,则电动汽车i的所有充放电功率迭代值P<sub>i,j</sub>的集合构成电动汽车i的充放电粒子P<sub>i</sub>,P<sub>i</sub>=[P<sub>i,1</sub>,P<sub>i,2</sub>…P<sub>i,j</sub>…P<sub>i,24</sub>];居民小区内所有电动汽车充放电粒子P<sub>i</sub>的集合构成小区充放电粒子P,P=[P<sub>1</sub>,P<sub>2</sub>…P<sub>i</sub>…P<sub>I</sub>];F、次日充放电参数的计算F1、计算归一化的次日负荷峰谷差<img file="FDA0000774451460000026.GIF" wi="84" he="66" />小区次日第j小时内的基础负荷L<sub>j</sub>加上所有电动汽车i在第j小时的充放电功率p<sub>i,j</sub>总和即为L′<sub>j</sub>小区次日第j小时的负荷L′<sub>j</sub>,<img file="FDA0000774451460000023.GIF" wi="387" he="142" />比较小区次日各小时负荷L′<sub>j</sub>,得出次日的最大负荷L′<sub>max</sub>和最小负荷L′<sub>min</sub>,<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mi>L</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mo>&prime;</mo></msubsup><mo>=</mo><mi>m</mi><mi>a</mi><mi>x</mi><mo>&lsqb;</mo><mrow><mo>(</mo><msub><mi>L</mi><mn>1</mn></msub><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>I</mi></munderover><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>,</mo><mrow><mo>(</mo><msub><mi>L</mi><mn>2</mn></msub><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>I</mi></munderover><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mn>2</mn></mrow></msub><mo>)</mo></mrow><mo>...</mo><mo>,</mo><mo>(</mo><mrow><msub><mi>L</mi><mn>24</mn></msub><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>I</mi></munderover><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mn>24</mn></mrow></msub></mrow><mo>)</mo><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000774451460000024.GIF" wi="1446" he="150" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msubsup><mi>L</mi><mi>min</mi><mo>&prime;</mo></msubsup><mo>=</mo><mi>min</mi><mo>&lsqb;</mo><mrow><mo>(</mo><msub><mi>L</mi><mn>1</mn></msub><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>I</mi></munderover><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>,</mo><mrow><mo>(</mo><msub><mi>L</mi><mn>2</mn></msub><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>I</mi></munderover><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mn>2</mn></mrow></msub><mo>)</mo></mrow><mo>...</mo><mo>,</mo><mo>(</mo><mrow><msub><mi>L</mi><mn>24</mn></msub><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>I</mi></munderover><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mn>24</mn></mrow></msub></mrow><mo>)</mo><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000774451460000025.GIF" wi="1442" he="142" /></maths>式中:max表示取最大值,min表示取最小值;进而得到次日负荷的峰谷差F<sub>1</sub>,F<sub>1</sub>=L′<sub>max</sub>‑L′<sub>min</sub>    (8)对次日负荷峰谷差F<sub>1</sub>用次日基础负荷峰谷差F<sub>1max</sub>进行归一化处理,得到归一化的次日负荷峰谷差<img file="FDA00007744514600000317.GIF" wi="85" he="66" /><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msubsup><mi>F</mi><mn>1</mn><mo>*</mo></msubsup><mo>=</mo><msub><mi>F</mi><mn>1</mn></msub><mo>/</mo><msub><mi>F</mi><mrow><mn>1</mn><mi>max</mi></mrow></msub><mo>;</mo></mrow>]]></math><img file="FDA00007744514600000318.GIF" wi="478" he="66" /></maths>F2、计算归一化的次日日负荷率<img file="FDA0000774451460000033.GIF" wi="99" he="83" />计算次日的平均负荷L′<sub>av</sub>,<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msubsup><mi>L</mi><mrow><mi>a</mi><mi>v</mi></mrow><mo>&prime;</mo></msubsup><mo>=</mo><mrow><mo>&lsqb;</mo><mrow><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mn>24</mn></munderover><mrow><mo>(</mo><mrow><msub><mi>L</mi><mi>j</mi></msub><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>I</mi></munderover><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub></mrow><mo>)</mo></mrow></mrow><mo>&rsqb;</mo></mrow><mo>/</mo><mn>24</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000774451460000034.GIF" wi="823" he="147" /></maths>算出次日的平均负荷L′<sub>av</sub>和次日的最大负荷L′<sub>max</sub>的比率,即得到次日的日负荷率F<sub>2</sub>,<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msub><mi>F</mi><mn>2</mn></msub><mo>=</mo><mfrac><msubsup><mi>L</mi><mrow><mi>a</mi><mi>v</mi></mrow><mo>&prime;</mo></msubsup><msubsup><mi>L</mi><mi>max</mi><mo>&prime;</mo></msubsup></mfrac><mo>&times;</mo><mn>100</mn><mi>%</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000774451460000035.GIF" wi="1054" he="146" /></maths>对次日日负荷率F<sub>2</sub>用次日基础负荷日负荷率F<sub>2min</sub>进行归一化处理,得到归一化的次日日负荷率<img file="FDA0000774451460000036.GIF" wi="81" he="79" /><maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><msubsup><mi>F</mi><mn>2</mn><mo>*</mo></msubsup><msub><mrow><mo>=</mo><mi>F</mi></mrow><mrow><mn>2</mn><mi>m</mi><mi>i</mi><mi>n</mi></mrow></msub><mo>/</mo><msub><mi>F</mi><mn>2</mn></msub><mo>;</mo></mrow>]]></math><img file="FDA0000774451460000037.GIF" wi="307" he="84" /></maths>F3、计算归一化的次日负荷波动均方差<img file="FDA0000774451460000038.GIF" wi="94" he="82" /><maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><msub><mi>F</mi><mn>3</mn></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mn>24</mn></munderover><msup><mrow><mo>(</mo><mrow><msub><mi>L</mi><mi>j</mi></msub><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>I</mi></munderover><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>-</mo><msubsup><mi>L</mi><mrow><mi>a</mi><mi>v</mi></mrow><mo>&prime;</mo></msubsup></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>/</mo><mn>24</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000774451460000039.GIF" wi="1197" he="152" /></maths>对次日负荷波动均方差F<sub>3</sub>用次日基础负荷日负荷率F<sub>3max</sub>进行归一化处理,得到归一化的次日负荷波动均方差<img file="FDA00007744514600000310.GIF" wi="80" he="83" /><maths num="0012" id="cmaths0012"><math><![CDATA[<mrow><msubsup><mi>F</mi><mn>3</mn><mo>*</mo></msubsup><mo>=</mo><msub><mi>F</mi><mn>3</mn></msub><mo>/</mo><msub><mi>F</mi><mrow><mn>3</mn><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>;</mo></mrow>]]></math><img file="FDA00007744514600000311.GIF" wi="310" he="84" /></maths>F4、计算归一化的次日车主费用<img file="FDA00007744514600000312.GIF" wi="96" he="79" /><maths num="0013" id="cmaths0013"><math><![CDATA[<mrow><msub><mi>F</mi><mn>4</mn></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>I</mi></munderover><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mn>24</mn></munderover><mrow><mo>(</mo><msub><mi>u</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><msub><mi>cp</mi><mi>j</mi></msub><mo>+</mo><msub><mi>v</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><msub><mi>dp</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>&Delta;</mi><mi>t</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00007744514600000313.GIF" wi="1261" he="150" /></maths>式中:u<sub>i,j</sub>、v<sub>i,j</sub>分别为第i辆电动汽车在第j小时的充电、放电标志符:充电时u<sub>i,j</sub>=1,v<sub>i,j</sub>=0,放电时u<sub>i,j</sub>=0,v<sub>i,j</sub>=1;cp<sub>j</sub>、dp<sub>j</sub>分别为第j小时的充电电价和放电电价;对次日的车主费用F<sub>4</sub>用次日车主的最大总费用F<sub>4max</sub>进行归一化处理,得到归一化的次日的车主费用<img file="FDA00007744514600000314.GIF" wi="82" he="78" /><maths num="0014" id="cmaths0014"><math><![CDATA[<mrow><msubsup><mi>F</mi><mn>4</mn><mo>*</mo></msubsup><mo>=</mo><msub><mi>F</mi><mn>4</mn></msub><mo>/</mo><msub><mi>F</mi><mrow><mn>4</mn><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>;</mo></mrow>]]></math><img file="FDA00007744514600000315.GIF" wi="311" he="78" /></maths>G、优化目标值的确定:将步骤F得到的四个归一化参数采用线性加权求和,得到优化目标值F:<maths num="0015" id="cmaths0015"><math><![CDATA[<mrow><mi>F</mi><mo>=</mo><msub><mi>w</mi><mn>1</mn></msub><msubsup><mi>F</mi><mn>1</mn><mo>*</mo></msubsup><mo>+</mo><msub><mi>w</mi><mn>2</mn></msub><msubsup><mi>F</mi><mn>2</mn><mo>*</mo></msubsup><mo>+</mo><msub><mi>w</mi><mn>3</mn></msub><msubsup><mi>F</mi><mn>3</mn><mo>*</mo></msubsup><mo>+</mo><msub><mi>w</mi><mn>4</mn></msub><msubsup><mi>F</mi><mn>4</mn><mo>*</mo></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00007744514600000316.GIF" wi="1204" he="90" /></maths>式中:w<sub>1</sub>、w<sub>2</sub>、w<sub>3</sub>分别是负荷峰谷差F<sub>1</sub>、日负荷率F<sub>2</sub>、波动均方差F<sub>3</sub>的权系数,其取值均为0.2,w<sub>4</sub>为车主总费用F<sub>4</sub>的权系数,其取值为0.4;H、第一次迭代重复步骤E‑至步骤F的操作L次,得到第一次迭代的L个优化目标值F和L个小区充放电粒子P,令第l个优化目标值为F<sup>1,l</sup>,第l个小区充放电粒子为P<sup>1,l</sup>,每个小区充放电粒子P<sup>1,l</sup>为第l次操作对应的所有电动汽车的次日分时充电功率p<sub>i,j</sub>的集合,即<img file="FDA0000774451460000041.GIF" wi="591" he="86" />最小的F<sup>1,l</sup>值所对应的小区充放电粒子P<sup>1,l</sup>即为全局最优小区充放电粒子P<sup>gbest</sup>,完成第一次迭代;I、第二次及以后各次迭代I1、比较第l个小区充放电粒子历次迭代中的F<sup>k,l</sup>值,若第k′次迭代中F<sup>k′,l</sup>值最小,则第k′次迭代中的P<sup>k′,l</sup>为第l个个体最优小区充放电粒子P<sup>l,pbest</sup>;I2、以步骤D中的式(5)为目标函数,以式(3)‑(4)为约束条件,将上次即k次迭代的全局最优小区充放电粒子P<sup>gbest</sup>以及第l个个体最优小区充放电粒子P<sup>l,pbest</sup>代入粒子群优化算法,迭代求解得到第k+1次,即本次迭代的小区充放电粒子P<sup>k+1,l</sup>,进而得到电动汽车i的充放电粒子<img file="FDA0000774451460000042.GIF" wi="130" he="66" />由<img file="FDA0000774451460000043.GIF" wi="582" he="76" />得到本次迭代第j小时电动汽车i的充放电功率的迭代值P<sub>i,j</sub>;I3、将迭代值P<sub>i,j</sub>重新代入步骤F和步骤G,得到第l个小区充放电粒子P<sup>k+1,l</sup>的优化目标值F<sup>k+1,l</sup>;I4、重复步骤I1‑I3操作L次,最小的F<sup>k+1,l</sup>值所对应的小区充放电粒子P<sup>k+1,l</sup>即为全局最优小区充放电粒子P<sup>gbest</sup>,完成本次迭代;I5、重复步骤I1‑I4操作K次,得到K次迭代时的全局最优小区充放电粒子P<sup>gbest</sup>;次日即控制电动汽车i按电动汽车i在次日第j小时的充放电功率等于P<sup>gbest</sup>的元素p<sub>i,j</sub>的值进行充放电。
地址 610031 四川省成都市二环路北一段111号