14.基于双层优化的电动汽车优化调度研究(文章复现)

说明书

MATLAB代码：基于双层优化的电动汽车优化调度研究

关键词：双层优化  选址定容 输配协同 时空优化    

参考文档：《考虑大规模电动汽车接入电网的双层优化调度策略_胡文平》中文版

《A bi-layer optimization based temporal and spatial scheduling for large-scale electric vehicles》完全复现

仿真平台：MATLAB+CPLEX 平台

优势：代码具有一定的深度和创新性，注释清晰，非烂大街的代码，非常精品！

主要内容：代码主要做的是一个双层的电动汽车充放电行为优化问题，具体来讲，输电网上层优化将电动汽车与发电机、基本负荷协调，同时考虑风力发电，从而在时域内优化电动汽车的负荷周期。然后，配电网的下层优化在空间上调度电动汽车负荷的位置。同时代码考虑了风电的出力场景，研究了不同风电出力下电动汽车的适应性，该代码具有一定的创新性，适合新手学习以及在此基础上进行拓展，代码质量非常高，保姆级的注释以及人性化的模块子程序，所有数据均有可靠来源

实现效果：具体如下

%%%%%%配电网潮流优化%%%%%%%%
%%%%%%标幺值SB=100MVA,UB=12.66kV,二阶锥松弛%%%%%%
%%%%%%MISOCP模型,分时段优化,并行计算%%%%%%%%%%%%
clear
clc
tic;
%%%%%IEEE33配电网数据%%%%%%%%%%%;
Pload=[0.0004666666667,0.0005,0.0005666666667,0.0006333333333,0.0006666666667,0.0007333333333,0.0007666666667,0.0008,0.0008666666667,0.0009333333333,0.0009666666667,0.001,0.0009333333333,0.0008666666667,0.0008,0.0007,0.0006666666667,0.0007333333333,0.0008,0.0009333333333,0.0008666666667,0.0007333333333,0.0006,0.0005333333333;0.00042,0.00045,0.00051,0.00057,0.0006,0.00066,0.00069,0.00072,0.00078,0.00084,0.00087,0.0009,0.00084,0.00078,0.00072,0.00063,0.0006,0.00066,0.00072,0.00084,0.00078,0.00066,0.00054,0.00048;0.00056,0.0006,0.00068,0.00076,0.0008,0.00088,0.00092,0.00096,0.00104,0.00112,0.00116,0.0012,0.00112,0.00104,0.00096,0.00084,0.0008,0.00088,0.00096,0.00112,0.00104,0.00088,0.00072,0.00064;0.00028,0.0003,0.00034,0.00038,0.0004,0.00044,0.00046,0.00048,0.00052,0.00056,0.00058,0.0006,0.00056,0.00052,0.00048,0.00042,0.0004,0.00044,0.00048,0.00056,0.00052,0.00044,0.00036,0.00032;0.00028,0.0003,0.00034,0.00038,0.0004,0.00044,0.00046,0.00048,0.00052,0.00056,0.00058,0.0006,0.00056,0.00052,0.00048,0.00042,0.0004,0.00044,0.00048,0.00056,0.00052,0.00044,0.00036,0.00032;0.0009333333333,0.001,0.001133333333,0.001266666667,0.001333333333,0.001466666667,0.001533333333,0.0016,0.001733333333,0.001866666667,0.001933333333,0.002,0.001866666667,0.001733333333,0.0016,0.0014,0.001333333333,0.001466666667,0.0016,0.001866666667,0.001733333333,0.001466666667,0.0012,0.001066666667;0.0009333333333,0.001,0.001133333333,0.001266666667,0.001333333333,0.001466666667,0.001533333333,0.0016,0.001733333333,0.001866666667,0.001933333333,0.002,0.001866666667,0.001733333333,0.0016,0.0014,0.001333333333,0.001466666667,0.0016,0.001866666667,0.001733333333,0.001466666667,0.0012,0.001066666667;0.00028,0.0003,0.00034,0.00038,0.0004,0.00044,0.00046,0.00048,0.00052,0.00056,0.00058,0.0006,0.00056,0.00052,0.00048,0.00042,0.0004,0.00044,0.00048,0.00056,0.00052,0.00044,0.00036,0.00032;0.00028,0.0003,0.00034,0.00038,0.0004,0.00044,0.00046,0.00048,0.00052,0.00056,0.00058,0.0006,0.00056,0.00052,0.00048,0.00042,0.0004,0.00044,0.00048,0.00056,0.00052,0.00044,0.00036,0.00032;0.00021,0.000225,0.000255,0.000285,0.0003,0.00033,0.000345,0.00036,0.00039,0.00042,0.000435,0.00045,0.00042,0.00039,0.00036,0.000315,0.0003,0.00033,0.00036,0.00042,0.00039,0.00033,0.00027,0.00024;0.00028,0.0003,0.00034,0.00038,0.0004,0.00044,0.00046,0.00048,0.00052,0.00056,0.00058,0.0006,0.00056,0.00052,0.00048,0.00042,0.0004,0.00044,0.00048,0.00056,0.00052,0.00044,0.00036,0.00032;0.00028,0.0003,0.00034,0.00038,0.0004,0.00044,0.00046,0.00048,0.00052,0.00056,0.00058,0.0006,0.00056,0.00052,0.00048,0.00042,0.0004,0.00044,0.00048,0.00056,0.00052,0.00044,0.00036,0.00032;0.00056,0.0006,0.00068,0.00076,0.0008,0.00088,0.00092,0.00096,0.00104,0.00112,0.00116,0.0012,0.00112,0.00104,0.00096,0.00084,0.0008,0.00088,0.00096,0.00112,0.00104,0.00088,0.00072,0.00064;0.00028,0.0003,0.00034,0.00038,0.0004,0.00044,0.00046,0.00048,0.00052,0.00056,0.00058,0.0006,0.00056,0.00052,0.00048,0.00042,0.0004,0.00044,0.00048,0.00056,0.00052,0.00044,0.00036,0.00032;0.00028,0.0003,0.00034,0.00038,0.0004,0.00044,0.00046,0.00048,0.00052,0.00056,0.00058,0.0006,0.00056,0.00052,0.00048,0.00042,0.0004,0.00044,0.00048,0.00056,0.00052,0.00044,0.00036,0.00032;0.00028,0.0003,0.00034,0.00038,0.0004,0.00044,0.00046,0.00048,0.00052,0.00056,0.00058,0.0006,0.00056,0.00052,0.00048,0.00042,0.0004,0.00044,0.00048,0.00056,0.00052,0.00044,0.00036,0.00032;0.00042,0.00045,0.00051,0.00057,0.0006,0.00066,0.00069,0.00072,0.00078,0.00084,0.00087,0.0009,0.00084,0.00078,0.00072,0.00063,0.0006,0.00066,0.00072,0.00084,0.00078,0.00066,0.00054,0.00048;0.00042,0.00045,0.00051,0.00057,0.0006,0.00066,0.00069,0.00072,0.00078,0.00084,0.00087,0.0009,0.00084,0.00078,0.00072,0.00063,0.0006,0.00066,0.00072,0.00084,0.00078,0.00066,0.00054,0.00048;0.00042,0.00045,0.00051,0.00057,0.0006,0.00066,0.00069,0.00072,0.00078,0.00084,0.00087,0.0009,0.00084,0.00078,0.00072,0.00063,0.0006,0.00066,0.00072,0.00084,0.00078,0.00066,0.00054,0.00048;0.00042,0.00045,0.00051,0.00057,0.0006,0.00066,0.00069,0.00072,0.00078,0.00084,0.00087,0.0009,0.00084,0.00078,0.00072,0.00063,0.0006,0.00066,0.00072,0.00084,0.00078,0.00066,0.00054,0.00048;0.00042,0.00045,0.00051,0.00057,0.0006,0.00066,0.00069,0.00072,0.00078,0.00084,0.00087,0.0009,0.00084,0.00078,0.00072,0.00063,0.0006,0.00066,0.00072,0.00084,0.00078,0.00066,0.00054,0.00048;0.00042,0.00045,0.00051,0.00057,0.0006,0.00066,0.00069,0.00072,0.00078,0.00084,0.00087,0.0009,0.00084,0.00078,0.00072,0.00063,0.0006,0.00066,0.00072,0.00084,0.00078,0.00066,0.00054,0.00048;0.00196,0.0021,0.00238,0.00266,0.0028,0.00308,0.00322,0.00336,0.00364,0.00392,0.00406,0.0042,0.00392,0.00364,0.00336,0.00294,0.0028,0.00308,0.00336,0.00392,0.00364,0.00308,0.00252,0.00224;0.00196,0.0021,0.00238,0.00266,0.0028,0.00308,0.00322,0.00336,0.00364,0.00392,0.00406,0.0042,0.00392,0.00364,0.00336,0.00294,0.0028,0.00308,0.00336,0.00392,0.00364,0.00308,0.00252,0.00224;0.00028,0.0003,0.00034,0.00038,0.0004,0.00044,0.00046,0.00048,0.00052,0.00056,0.00058,0.0006,0.00056,0.00052,0.00048,0.00042,0.0004,0.00044,0.00048,0.00056,0.00052,0.00044,0.00036,0.00032;0.00028,0.0003,0.00034,0.00038,0.0004,0.00044,0.00046,0.00048,0.00052,0.00056,0.00058,0.0006,0.00056,0.00052,0.00048,0.00042,0.0004,0.00044,0.00048,0.00056,0.00052,0.00044,0.00036,0.00032;0.00028,0.0003,0.00034,0.00038,0.0004,0.00044,0.00046,0.00048,0.00052,0.00056,0.00058,0.0006,0.00056,0.00052,0.00048,0.00042,0.0004,0.00044,0.00048,0.00056,0.00052,0.00044,0.00036,0.00032;0.00056,0.0006,0.00068,0.00076,0.0008,0.00088,0.00092,0.00096,0.00104,0.00112,0.00116,0.0012,0.00112,0.00104,0.00096,0.00084,0.0008,0.00088,0.00096,0.00112,0.00104,0.00088,0.00072,0.00064;0.0009333333333,0.001,0.001133333333,0.001266666667,0.001333333333,0.001466666667,0.001533333333,0.0016,0.001733333333,0.001866666667,0.001933333333,0.002,0.001866666667,0.001733333333,0.0016,0.0014,0.001333333333,0.001466666667,0.0016,0.001866666667,0.001733333333,0.001466666667,0.0012,0.001066666667;0.0007,0.00075,0.00085,0.00095,0.001,0.0011,0.00115,0.0012,0.0013,0.0014,0.00145,0.0015,0.0014,0.0013,0.0012,0.00105,0.001,0.0011,0.0012,0.0014,0.0013,0.0011,0.0009,0.0008;0.00098,0.00105,0.00119,0.00133,0.0014,0.00154,0.00161,0.00168,0.00182,0.00196,0.00203,0.0021,0.00196,0.00182,0.00168,0.00147,0.0014,0.00154,0.00168,0.00196,0.00182,0.00154,0.00126,0.00112;0.00028,0.0003,0.00034,0.00038,0.0004,0.00044,0.00046,0.00048,0.00052,0.00056,0.00058,0.0006,0.00056,0.00052,0.00048,0.00042,0.0004,0.00044,0.00048,0.00056,0.00052,0.00044,0.00036,0.00032;0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0];%节点有功(何立夫数据）;
Qload=[0.00028,0.0003,0.00034,0.00038,0.0004,0.00044,0.00046,0.00048,0.00052,0.00056,0.00058,0.0006,0.00056,0.00052,0.00048,0.00042,0.0004,0.00044,0.00048,0.00056,0.00052,0.00044,0.00036,0.00032;0.0001866666667,0.0002,0.0002266666667,0.0002533333333,0.0002666666667,0.0002933333333,0.0003066666667,0.00032,0.0003466666667,0.0003733333333,0.0003866666667,0.0004,0.0003733333333,0.0003466666667,0.00032,0.00028,0.0002666666667,0.0002933333333,0.00032,0.0003733333333,0.0003466666667,0.0002933333333,0.00024,0.0002133333333;0.0003733333333,0.0004,0.0004533333333,0.0005066666667,0.0005333333333,0.0005866666667,0.0006133333333,0.00064,0.0006933333333,0.0007466666667,0.0007733333333,0.0008,0.0007466666667,0.0006933333333,0.00064,0.00056,0.0005333333333,0.0005866666667,0.00064,0.0007466666667,0.0006933333333,0.0005866666667,0.00048,0.0004266666667;0.00014,0.00015,0.00017,0.00019,0.0002,0.00022,0.00023,0.00024,0.00026,0.00028,0.00029,0.0003,0.00028,0.00026,0.00024,0.00021,0.0002,0.00022,0.00024,0.00028,0.00026,0.00022,0.00018,0.00016;9.333333333e-05,0.0001,0.0001133333333,0.0001266666667,0.0001333333333,0.0001466666667,0.0001533333333,0.00016,0.0001733333333,0.0001866666667,0.0001933333333,0.0002,0.0001866666667,0.0001733333333,0.00016,0.00014,0.0001333333333,0.0001466666667,0.00016,0.0001866666667,0.0001733333333,0.0001466666667,0.00012,0.0001066666667;0.0004666666667,0.0005,0.0005666666667,0.0006333333333,0.0006666666667,0.0007333333333,0.0007666666667,0.0008,0.0008666666667,0.0009333333333,0.0009666666667,0.001,0.0009333333333,0.0008666666667,0.0008,0.0007,0.0006666666667,0.0007333333333,0.0008,0.0009333333333,0.0008666666667,0.0007333333333,0.0006,0.0005333333333;0.0004666666667,0.0005,0.0005666666667,0.0006333333333,0.0006666666667,0.0007333333333,0.0007666666667,0.0008,0.0008666666667,0.0009333333333,0.0009666666667,0.001,0.0009333333333,0.0008666666667,0.0008,0.0007,0.0006666666667,0.0007333333333,0.0008,0.0009333333333,0.0008666666667,0.0007333333333,0.0006,0.0005333333333;9.333333333e-05,0.0001,0.0001133333333,0.0001266666667,0.0001333333333,0.0001466666667,0.0001533333333,0.00016,0.0001733333333,0.0001866666667,0.0001933333333,0.0002,0.0001866666667,0.0001733333333,0.00016,0.00014,0.0001333333333,0.0001466666667,0.00016,0.0001866666667,0.0001733333333,0.0001466666667,0.00012,0.0001066666667;9.333333333e-05,0.0001,0.0001133333333,0.0001266666667,0.0001333333333,0.0001466666667,0.0001533333333,0.00016,0.0001733333333,0.0001866666667,0.0001933333333,0.0002,0.0001866666667,0.0001733333333,0.00016,0.00014,0.0001333333333,0.0001466666667,0.00016,0.0001866666667,0.0001733333333,0.0001466666667,0.00012,0.0001066666667;0.00014,0.00015,0.00017,0.00019,0.0002,0.00022,0.00023,0.00024,0.00026,0.00028,0.00029,0.0003,0.00028,0.00026,0.00024,0.00021,0.0002,0.00022,0.00024,0.00028,0.00026,0.00022,0.00018,0.00016;0.0001633333333,0.000175,0.0001983333333,0.0002216666667,0.0002333333333,0.0002566666667,0.0002683333333,0.00028,0.0003033333333,0.0003266666667,0.0003383333333,0.00035,0.0003266666667,0.0003033333333,0.00028,0.000245,0.0002333333333,0.0002566666667,0.00028,0.0003266666667,0.0003033333333,0.0002566666667,0.00021,0.0001866666667;0.0001633333333,0.000175,0.0001983333333,0.0002216666667,0.0002333333333,0.0002566666667,0.0002683333333,0.00028,0.0003033333333,0.0003266666667,0.0003383333333,0.00035,0.0003266666667,0.0003033333333,0.00028,0.000245,0.0002333333333,0.0002566666667,0.00028,0.0003266666667,0.0003033333333,0.0002566666667,0.00021,0.0001866666667;0.0003733333333,0.0004,0.0004533333333,0.0005066666667,0.0005333333333,0.0005866666667,0.0006133333333,0.00064,0.0006933333333,0.0007466666667,0.0007733333333,0.0008,0.0007466666667,0.0006933333333,0.00064,0.00056,0.0005333333333,0.0005866666667,0.00064,0.0007466666667,0.0006933333333,0.0005866666667,0.00048,0.0004266666667;4.666666667e-05,5.e-05,5.666666667e-05,6.333333333e-05,6.666666667e-05,7.333333333e-05,7.666666667e-05,8.e-05,8.666666667e-05,9.333333333e-05,9.666666667e-05,0.0001,9.333333333e-05,8.666666667e-05,8.e-05,7.e-05,6.666666667e-05,7.333333333e-05,8.e-05,9.333333333e-05,8.666666667e-05,7.333333333e-05,6.e-05,5.333333333e-05;9.333333333e-05,0.0001,0.0001133333333,0.0001266666667,0.0001333333333,0.0001466666667,0.0001533333333,0.00016,0.0001733333333,0.0001866666667,0.0001933333333,0.0002,0.0001866666667,0.0001733333333,0.00016,0.00014,0.0001333333333,0.0001466666667,0.00016,0.0001866666667,0.0001733333333,0.0001466666667,0.00012,0.0001066666667;9.333333333e-05,0.0001,0.0001133333333,0.0001266666667,0.0001333333333,0.0001466666667,0.0001533333333,0.00016,0.0001733333333,0.0001866666667,0.0001933333333,0.0002,0.0001866666667,0.0001733333333,0.00016,0.00014,0.0001333333333,0.0001466666667,0.00016,0.0001866666667,0.0001733333333,0.0001466666667,0.00012,0.0001066666667;0.0001866666667,0.0002,0.0002266666667,0.0002533333333,0.0002666666667,0.0002933333333,0.0003066666667,0.00032,0.0003466666667,0.0003733333333,0.0003866666667,0.0004,0.0003733333333,0.0003466666667,0.00032,0.00028,0.0002666666667,0.0002933333333,0.00032,0.0003733333333,0.0003466666667,0.0002933333333,0.00024,0.0002133333333;0.0001866666667,0.0002,0.0002266666667,0.0002533333333,0.0002666666667,0.0002933333333,0.0003066666667,0.00032,0.0003466666667,0.0003733333333,0.0003866666667,0.0004,0.0003733333333,0.0003466666667,0.00032,0.00028,0.0002666666667,0.0002933333333,0.00032,0.0003733333333,0.0003466666667,0.0002933333333,0.00024,0.0002133333333;0.0001866666667,0.0002,0.0002266666667,0.0002533333333,0.0002666666667,0.0002933333333,0.0003066666667,0.00032,0.0003466666667,0.0003733333333,0.0003866666667,0.0004,0.0003733333333,0.0003466666667,0.00032,0.00028,0.0002666666667,0.0002933333333,0.00032,0.0003733333333,0.0003466666667,0.0002933333333,0.00024,0.0002133333333;0.0001866666667,0.0002,0.0002266666667,0.0002533333333,0.0002666666667,0.0002933333333,0.0003066666667,0.00032,0.0003466666667,0.0003733333333,0.0003866666667,0.0004,0.0003733333333,0.0003466666667,0.00032,0.00028,0.0002666666667,0.0002933333333,0.00032,0.0003733333333,0.0003466666667,0.0002933333333,0.00024,0.0002133333333;0.0001866666667,0.0002,0.0002266666667,0.0002533333333,0.0002666666667,0.0002933333333,0.0003066666667,0.00032,0.0003466666667,0.0003733333333,0.0003866666667,0.0004,0.0003733333333,0.0003466666667,0.00032,0.00028,0.0002666666667,0.0002933333333,0.00032,0.0003733333333,0.0003466666667,0.0002933333333,0.00024,0.0002133333333;0.0002333333333,0.00025,0.0002833333333,0.0003166666667,0.0003333333333,0.0003666666667,0.0003833333333,0.0004,0.0004333333333,0.0004666666667,0.0004833333333,0.0005,0.0004666666667,0.0004333333333,0.0004,0.00035,0.0003333333333,0.0003666666667,0.0004,0.0004666666667,0.0004333333333,0.0003666666667,0.0003,0.0002666666667;0.0009333333333,0.001,0.001133333333,0.001266666667,0.001333333333,0.001466666667,0.001533333333,0.0016,0.001733333333,0.001866666667,0.001933333333,0.002,0.001866666667,0.001733333333,0.0016,0.0014,0.001333333333,0.001466666667,0.0016,0.001866666667,0.001733333333,0.001466666667,0.0012,0.001066666667;0.0009333333333,0.001,0.001133333333,0.001266666667,0.001333333333,0.001466666667,0.001533333333,0.0016,0.001733333333,0.001866666667,0.001933333333,0.002,0.001866666667,0.001733333333,0.0016,0.0014,0.001333333333,0.001466666667,0.0016,0.001866666667,0.001733333333,0.001466666667,0.0012,0.001066666667;0.0001166666667,0.000125,0.0001416666667,0.0001583333333,0.0001666666667,0.0001833333333,0.0001916666667,0.0002,0.0002166666667,0.0002333333333,0.0002416666667,0.00025,0.0002333333333,0.0002166666667,0.0002,0.000175,0.0001666666667,0.0001833333333,0.0002,0.0002333333333,0.0002166666667,0.0001833333333,0.00015,0.0001333333333;0.0001166666667,0.000125,0.0001416666667,0.0001583333333,0.0001666666667,0.0001833333333,0.0001916666667,0.0002,0.0002166666667,0.0002333333333,0.0002416666667,0.00025,0.0002333333333,0.0002166666667,0.0002,0.000175,0.0001666666667,0.0001833333333,0.0002,0.0002333333333,0.0002166666667,0.0001833333333,0.00015,0.0001333333333;9.333333333e-05,0.0001,0.0001133333333,0.0001266666667,0.0001333333333,0.0001466666667,0.0001533333333,0.00016,0.0001733333333,0.0001866666667,0.0001933333333,0.0002,0.0001866666667,0.0001733333333,0.00016,0.00014,0.0001333333333,0.0001466666667,0.00016,0.0001866666667,0.0001733333333,0.0001466666667,0.00012,0.0001066666667;0.0003266666667,0.00035,0.0003966666667,0.0004433333333,0.0004666666667,0.0005133333333,0.0005366666667,0.00056,0.0006066666667,0.0006533333333,0.0006766666667,0.0007,0.0006533333333,0.0006066666667,0.00056,0.00049,0.0004666666667,0.0005133333333,0.00056,0.0006533333333,0.0006066666667,0.0005133333333,0.00042,0.0003733333333;0.0028,0.003,0.0034,0.0038,0.004,0.0044,0.0046,0.0048,0.0052,0.0056,0.0058,0.006,0.0056,0.0052,0.0048,0.0042,0.004,0.0044,0.0048,0.0056,0.0052,0.0044,0.0036,0.0032;0.0003266666667,0.00035,0.0003966666667,0.0004433333333,0.0004666666667,0.0005133333333,0.0005366666667,0.00056,0.0006066666667,0.0006533333333,0.0006766666667,0.0007,0.0006533333333,0.0006066666667,0.00056,0.00049,0.0004666666667,0.0005133333333,0.00056,0.0006533333333,0.0006066666667,0.0005133333333,0.00042,0.0003733333333;0.0004666666667,0.0005,0.0005666666667,0.0006333333333,0.0006666666667,0.0007333333333,0.0007666666667,0.0008,0.0008666666667,0.0009333333333,0.0009666666667,0.001,0.0009333333333,0.0008666666667,0.0008,0.0007,0.0006666666667,0.0007333333333,0.0008,0.0009333333333,0.0008666666667,0.0007333333333,0.0006,0.0005333333333;0.0001866666667,0.0002,0.0002266666667,0.0002533333333,0.0002666666667,0.0002933333333,0.0003066666667,0.00032,0.0003466666667,0.0003733333333,0.0003866666667,0.0004,0.0003733333333,0.0003466666667,0.00032,0.00028,0.0002666666667,0.0002933333333,0.00032,0.0003733333333,0.0003466666667,0.0002933333333,0.00024,0.0002133333333;0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0];%节点无功(何立夫数据）;
%Pload=[0.00058,0.00065,0.00071,0.00083,0.00088,0.00087,0.00084,0.000798,0.000804,0.000798,0.000815,0.00085,0.0009,0.000898,0.000975,0.000993,0.000895,0.00078,0.00069,0.00061,0.00058,0.000582,0.00057,0.000563;0.000522,0.000585,0.000639,0.000747,0.000792,0.000783,0.000756,0.0007182,0.0007236,0.0007182,0.0007335,0.000765,0.00081,0.0008082,0.0008775,0.0008937,0.0008055,0.000702,0.000621,0.000549,0.000522,0.0005238,0.000513,0.0005067;0.000696,0.00078,0.000852,0.000996,0.001056,0.001044,0.001008,0.0009576,0.0009648,0.0009576,0.000978,0.00102,0.00108,0.0010776,0.00117,0.0011916,0.001074,0.000936,0.000828,0.000732,0.000696,0.0006984,0.000684,0.0006756;0.000348,0.00039,0.000426,0.000498,0.000528,0.000522,0.000504,0.0004788,0.0004824,0.0004788,0.000489,0.00051,0.00054,0.0005388,0.000585,0.0005958,0.000537,0.000468,0.000414,0.000366,0.000348,0.0003492,0.000342,0.0003378;0.000348,0.00039,0.000426,0.000498,0.000528,0.000522,0.000504,0.0004788,0.0004824,0.0004788,0.000489,0.00051,0.00054,0.0005388,0.000585,0.0005958,0.000537,0.000468,0.000414,0.000366,0.000348,0.0003492,0.000342,0.0003378;0.00116,0.0013,0.00142,0.00166,0.00176,0.00174,0.00168,0.001596,0.001608,0.001596,0.00163,0.0017,0.0018,0.001796,0.00195,0.001986,0.00179,0.00156,0.00138,0.00122,0.00116,0.001164,0.00114,0.001126;0.00116,0.0013,0.00142,0.00166,0.00176,0.00174,0.00168,0.001596,0.001608,0.001596,0.00163,0.0017,0.0018,0.001796,0.00195,0.001986,0.00179,0.00156,0.00138,0.00122,0.00116,0.001164,0.00114,0.001126;0.000348,0.00039,0.000426,0.000498,0.000528,0.000522,0.000504,0.0004788,0.0004824,0.0004788,0.000489,0.00051,0.00054,0.0005388,0.000585,0.0005958,0.000537,0.000468,0.000414,0.000366,0.000348,0.0003492,0.000342,0.0003378;0.000348,0.00039,0.000426,0.000498,0.000528,0.000522,0.000504,0.0004788,0.0004824,0.0004788,0.000489,0.00051,0.00054,0.0005388,0.000585,0.0005958,0.000537,0.000468,0.000414,0.000366,0.000348,0.0003492,0.000342,0.0003378;0.000261,0.0002925,0.0003195,0.0003735,0.000396,0.0003915,0.000378,0.0003591,0.0003618,0.0003591,0.00036675,0.0003825,0.000405,0.0004041,0.00043875,0.00044685,0.00040275,0.000351,0.0003105,0.0002745,0.000261,0.0002619,0.0002565,0.00025335;0.000348,0.00039,0.000426,0.000498,0.000528,0.000522,0.000504,0.0004788,0.0004824,0.0004788,0.000489,0.00051,0.00054,0.0005388,0.000585,0.0005958,0.000537,0.000468,0.000414,0.000366,0.000348,0.0003492,0.000342,0.0003378;0.000348,0.00039,0.000426,0.000498,0.000528,0.000522,0.000504,0.0004788,0.0004824,0.0004788,0.000489,0.00051,0.00054,0.0005388,0.000585,0.0005958,0.000537,0.000468,0.000414,0.000366,0.000348,0.0003492,0.000342,0.0003378;0.000696,0.00078,0.000852,0.000996,0.001056,0.001044,0.001008,0.0009576,0.0009648,0.0009576,0.000978,0.00102,0.00108,0.0010776,0.00117,0.0011916,0.001074,0.000936,0.000828,0.000732,0.000696,0.0006984,0.000684,0.0006756;0.000348,0.00039,0.000426,0.000498,0.000528,0.000522,0.000504,0.0004788,0.0004824,0.0004788,0.000489,0.00051,0.00054,0.0005388,0.000585,0.0005958,0.000537,0.000468,0.000414,0.000366,0.000348,0.0003492,0.000342,0.0003378;0.000348,0.00039,0.000426,0.000498,0.000528,0.000522,0.000504,0.0004788,0.0004824,0.0004788,0.000489,0.00051,0.00054,0.0005388,0.000585,0.0005958,0.000537,0.000468,0.000414,0.000366,0.000348,0.0003492,0.000342,0.0003378;0.000348,0.00039,0.000426,0.000498,0.000528,0.000522,0.000504,0.0004788,0.0004824,0.0004788,0.000489,0.00051,0.00054,0.0005388,0.000585,0.0005958,0.000537,0.000468,0.000414,0.000366,0.000348,0.0003492,0.000342,0.0003378;0.000522,0.000585,0.000639,0.000747,0.000792,0.000783,0.000756,0.0007182,0.0007236,0.0007182,0.0007335,0.000765,0.00081,0.0008082,0.0008775,0.0008937,0.0008055,0.000702,0.000621,0.000549,0.000522,0.0005238,0.000513,0.0005067;0.000522,0.000585,0.000639,0.000747,0.000792,0.000783,0.000756,0.0007182,0.0007236,0.0007182,0.0007335,0.000765,0.00081,0.0008082,0.0008775,0.0008937,0.0008055,0.000702,0.000621,0.000549,0.000522,0.0005238,0.000513,0.0005067;0.000522,0.000585,0.000639,0.000747,0.000792,0.000783,0.000756,0.0007182,0.0007236,0.0007182,0.0007335,0.000765,0.00081,0.0008082,0.0008775,0.0008937,0.0008055,0.000702,0.000621,0.000549,0.000522,0.0005238,0.000513,0.0005067;0.000522,0.000585,0.000639,0.000747,0.000792,0.000783,0.000756,0.0007182,0.0007236,0.0007182,0.0007335,0.000765,0.00081,0.0008082,0.0008775,0.0008937,0.0008055,0.000702,0.000621,0.000549,0.000522,0.0005238,0.000513,0.0005067;0.000522,0.000585,0.000639,0.000747,0.000792,0.000783,0.000756,0.0007182,0.0007236,0.0007182,0.0007335,0.000765,0.00081,0.0008082,0.0008775,0.0008937,0.0008055,0.000702,0.000621,0.000549,0.000522,0.0005238,0.000513,0.0005067;0.000522,0.000585,0.000639,0.000747,0.000792,0.000783,0.000756,0.0007182,0.0007236,0.0007182,0.0007335,0.000765,0.00081,0.0008082,0.0008775,0.0008937,0.0008055,0.000702,0.000621,0.000549,0.000522,0.0005238,0.000513,0.0005067;0.002436,0.00273,0.002982,0.003486,0.003696,0.003654,0.003528,0.0033516,0.0033768,0.0033516,0.003423,0.00357,0.00378,0.0037716,0.004095,0.0041706,0.003759,0.003276,0.002898,0.002562,0.002436,0.0024444,0.002394,0.0023646;0.002436,0.00273,0.002982,0.003486,0.003696,0.003654,0.003528,0.0033516,0.0033768,0.0033516,0.003423,0.00357,0.00378,0.0037716,0.004095,0.0041706,0.003759,0.003276,0.002898,0.002562,0.002436,0.0024444,0.002394,0.0023646;0.000348,0.00039,0.000426,0.000498,0.000528,0.000522,0.000504,0.0004788,0.0004824,0.0004788,0.000489,0.00051,0.00054,0.0005388,0.000585,0.0005958,0.000537,0.000468,0.000414,0.000366,0.000348,0.0003492,0.000342,0.0003378;0.000348,0.00039,0.000426,0.000498,0.000528,0.000522,0.000504,0.0004788,0.0004824,0.0004788,0.000489,0.00051,0.00054,0.0005388,0.000585,0.0005958,0.000537,0.000468,0.000414,0.000366,0.000348,0.0003492,0.000342,0.0003378;0.000348,0.00039,0.000426,0.000498,0.000528,0.000522,0.000504,0.0004788,0.0004824,0.0004788,0.000489,0.00051,0.00054,0.0005388,0.000585,0.0005958,0.000537,0.000468,0.000414,0.000366,0.000348,0.0003492,0.000342,0.0003378;0.000696,0.00078,0.000852,0.000996,0.001056,0.001044,0.001008,0.0009576,0.0009648,0.0009576,0.000978,0.00102,0.00108,0.0010776,0.00117,0.0011916,0.001074,0.000936,0.000828,0.000732,0.000696,0.0006984,0.000684,0.0006756;0.00116,0.0013,0.00142,0.00166,0.00176,0.00174,0.00168,0.001596,0.001608,0.001596,0.00163,0.0017,0.0018,0.001796,0.00195,0.001986,0.00179,0.00156,0.00138,0.00122,0.00116,0.001164,0.00114,0.001126;0.00087,0.000975,0.001065,0.001245,0.00132,0.001305,0.00126,0.001197,0.001206,0.001197,0.0012225,0.001275,0.00135,0.001347,0.0014625,0.0014895,0.0013425,0.00117,0.001035,0.000915,0.00087,0.000873,0.000855,0.0008445;0.001218,0.001365,0.001491,0.001743,0.001848,0.001827,0.001764,0.0016758,0.0016884,0.0016758,0.0017115,0.001785,0.00189,0.0018858,0.0020475,0.0020853,0.0018795,0.001638,0.001449,0.001281,0.001218,0.0012222,0.001197,0.0011823;0.000348,0.00039,0.000426,0.000498,0.000528,0.000522,0.000504,0.0004788,0.0004824,0.0004788,0.000489,0.00051,0.00054,0.0005388,0.000585,0.0005958,0.000537,0.000468,0.000414,0.000366,0.000348,0.0003492,0.000342,0.0003378;0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0];%节点有功;
%Qload=[0.000348,0.00039,0.000426,0.000498,0.000528,0.000522,0.000504,0.0004788,0.0004824,0.0004788,0.000489,0.00051,0.00054,0.0005388,0.000585,0.0005958,0.000537,0.000468,0.000414,0.000366,0.000348,0.0003492,0.000342,0.0003378;0.000232,0.00026,0.000284,0.000332,0.000352,0.000348,0.000336,0.0003192,0.0003216,0.0003192,0.000326,0.00034,0.00036,0.0003592,0.00039,0.0003972,0.000358,0.000312,0.000276,0.000244,0.000232,0.0002328,0.000228,0.0002252;0.000464,0.00052,0.000568,0.000664,0.000704,0.000696,0.000672,0.0006384,0.0006432,0.0006384,0.000652,0.00068,0.00072,0.0007184,0.00078,0.0007944,0.000716,0.000624,0.000552,0.000488,0.000464,0.0004656,0.000456,0.0004504;0.000174,0.000195,0.000213,0.000249,0.000264,0.000261,0.000252,0.0002394,0.0002412,0.0002394,0.0002445,0.000255,0.00027,0.0002694,0.0002925,0.0002979,0.0002685,0.000234,0.000207,0.000183,0.000174,0.0001746,0.000171,0.0001689;0.000116,0.00013,0.000142,0.000166,0.000176,0.000174,0.000168,0.0001596,0.0001608,0.0001596,0.000163,0.00017,0.00018,0.0001796,0.000195,0.0001986,0.000179,0.000156,0.000138,0.000122,0.000116,0.0001164,0.000114,0.0001126;0.00058,0.00065,0.00071,0.00083,0.00088,0.00087,0.00084,0.000798,0.000804,0.000798,0.000815,0.00085,0.0009,0.000898,0.000975,0.000993,0.000895,0.00078,0.00069,0.00061,0.00058,0.000582,0.00057,0.000563;0.00058,0.00065,0.00071,0.00083,0.00088,0.00087,0.00084,0.000798,0.000804,0.000798,0.000815,0.00085,0.0009,0.000898,0.000975,0.000993,0.000895,0.00078,0.00069,0.00061,0.00058,0.000582,0.00057,0.000563;0.000116,0.00013,0.000142,0.000166,0.000176,0.000174,0.000168,0.0001596,0.0001608,0.0001596,0.000163,0.00017,0.00018,0.0001796,0.000195,0.0001986,0.000179,0.000156,0.000138,0.000122,0.000116,0.0001164,0.000114,0.0001126;0.000116,0.00013,0.000142,0.000166,0.000176,0.000174,0.000168,0.0001596,0.0001608,0.0001596,0.000163,0.00017,0.00018,0.0001796,0.000195,0.0001986,0.000179,0.000156,0.000138,0.000122,0.000116,0.0001164,0.000114,0.0001126;0.000174,0.000195,0.000213,0.000249,0.000264,0.000261,0.000252,0.0002394,0.0002412,0.0002394,0.0002445,0.000255,0.00027,0.0002694,0.0002925,0.0002979,0.0002685,0.000234,0.000207,0.000183,0.000174,0.0001746,0.000171,0.0001689;0.000203,0.0002275,0.0002485,0.0002905,0.000308,0.0003045,0.000294,0.0002793,0.0002814,0.0002793,0.00028525,0.0002975,0.000315,0.0003143,0.00034125,0.00034755,0.00031325,0.000273,0.0002415,0.0002135,0.000203,0.0002037,0.0001995,0.00019705;0.000203,0.0002275,0.0002485,0.0002905,0.000308,0.0003045,0.000294,0.0002793,0.0002814,0.0002793,0.00028525,0.0002975,0.000315,0.0003143,0.00034125,0.00034755,0.00031325,0.000273,0.0002415,0.0002135,0.000203,0.0002037,0.0001995,0.00019705;0.000464,0.00052,0.000568,0.000664,0.000704,0.000696,0.000672,0.0006384,0.0006432,0.0006384,0.000652,0.00068,0.00072,0.0007184,0.00078,0.0007944,0.000716,0.000624,0.000552,0.000488,0.000464,0.0004656,0.000456,0.0004504;5.8e-05,6.5e-05,7.1e-05,8.3e-05,8.8e-05,8.7e-05,8.4e-05,7.98e-05,8.04e-05,7.98e-05,8.15e-05,8.5e-05,9.e-05,8.98e-05,9.75e-05,9.93e-05,8.95e-05,7.8e-05,6.9e-05,6.1e-05,5.8e-05,5.82e-05,5.7e-05,5.63e-05;0.000116,0.00013,0.000142,0.000166,0.000176,0.000174,0.000168,0.0001596,0.0001608,0.0001596,0.000163,0.00017,0.00018,0.0001796,0.000195,0.0001986,0.000179,0.000156,0.000138,0.000122,0.000116,0.0001164,0.000114,0.0001126;0.000116,0.00013,0.000142,0.000166,0.000176,0.000174,0.000168,0.0001596,0.0001608,0.0001596,0.000163,0.00017,0.00018,0.0001796,0.000195,0.0001986,0.000179,0.000156,0.000138,0.000122,0.000116,0.0001164,0.000114,0.0001126;0.000232,0.00026,0.000284,0.000332,0.000352,0.000348,0.000336,0.0003192,0.0003216,0.0003192,0.000326,0.00034,0.00036,0.0003592,0.00039,0.0003972,0.000358,0.000312,0.000276,0.000244,0.000232,0.0002328,0.000228,0.0002252;0.000232,0.00026,0.000284,0.000332,0.000352,0.000348,0.000336,0.0003192,0.0003216,0.0003192,0.000326,0.00034,0.00036,0.0003592,0.00039,0.0003972,0.000358,0.000312,0.000276,0.000244,0.000232,0.0002328,0.000228,0.0002252;0.000232,0.00026,0.000284,0.000332,0.000352,0.000348,0.000336,0.0003192,0.0003216,0.0003192,0.000326,0.00034,0.00036,0.0003592,0.00039,0.0003972,0.000358,0.000312,0.000276,0.000244,0.000232,0.0002328,0.000228,0.0002252;0.000232,0.00026,0.000284,0.000332,0.000352,0.000348,0.000336,0.0003192,0.0003216,0.0003192,0.000326,0.00034,0.00036,0.0003592,0.00039,0.0003972,0.000358,0.000312,0.000276,0.000244,0.000232,0.0002328,0.000228,0.0002252;0.000232,0.00026,0.000284,0.000332,0.000352,0.000348,0.000336,0.0003192,0.0003216,0.0003192,0.000326,0.00034,0.00036,0.0003592,0.00039,0.0003972,0.000358,0.000312,0.000276,0.000244,0.000232,0.0002328,0.000228,0.0002252;0.00029,0.000325,0.000355,0.000415,0.00044,0.000435,0.00042,0.000399,0.000402,0.000399,0.0004075,0.000425,0.00045,0.000449,0.0004875,0.0004965,0.0004475,0.00039,0.000345,0.000305,0.00029,0.000291,0.000285,0.0002815;0.00116,0.0013,0.00142,0.00166,0.00176,0.00174,0.00168,0.001596,0.001608,0.001596,0.00163,0.0017,0.0018,0.001796,0.00195,0.001986,0.00179,0.00156,0.00138,0.00122,0.00116,0.001164,0.00114,0.001126;0.00116,0.0013,0.00142,0.00166,0.00176,0.00174,0.00168,0.001596,0.001608,0.001596,0.00163,0.0017,0.0018,0.001796,0.00195,0.001986,0.00179,0.00156,0.00138,0.00122,0.00116,0.001164,0.00114,0.001126;0.000145,0.0001625,0.0001775,0.0002075,0.00022,0.0002175,0.00021,0.0001995,0.000201,0.0001995,0.00020375,0.0002125,0.000225,0.0002245,0.00024375,0.00024825,0.00022375,0.000195,0.0001725,0.0001525,0.000145,0.0001455,0.0001425,0.00014075;0.000145,0.0001625,0.0001775,0.0002075,0.00022,0.0002175,0.00021,0.0001995,0.000201,0.0001995,0.00020375,0.0002125,0.000225,0.0002245,0.00024375,0.00024825,0.00022375,0.000195,0.0001725,0.0001525,0.000145,0.0001455,0.0001425,0.00014075;0.000116,0.00013,0.000142,0.000166,0.000176,0.000174,0.000168,0.0001596,0.0001608,0.0001596,0.000163,0.00017,0.00018,0.0001796,0.000195,0.0001986,0.000179,0.000156,0.000138,0.000122,0.000116,0.0001164,0.000114,0.0001126;0.000406,0.000455,0.000497,0.000581,0.000616,0.000609,0.000588,0.0005586,0.0005628,0.0005586,0.0005705,0.000595,0.00063,0.0006286,0.0006825,0.0006951,0.0006265,0.000546,0.000483,0.000427,0.000406,0.0004074,0.000399,0.0003941;0.00348,0.0039,0.00426,0.00498,0.00528,0.00522,0.00504,0.004788,0.004824,0.004788,0.00489,0.0051,0.0054,0.005388,0.00585,0.005958,0.00537,0.00468,0.00414,0.00366,0.00348,0.003492,0.00342,0.003378;0.000406,0.000455,0.000497,0.000581,0.000616,0.000609,0.000588,0.0005586,0.0005628,0.0005586,0.0005705,0.000595,0.00063,0.0006286,0.0006825,0.0006951,0.0006265,0.000546,0.000483,0.000427,0.000406,0.0004074,0.000399,0.0003941;0.00058,0.00065,0.00071,0.00083,0.00088,0.00087,0.00084,0.000798,0.000804,0.000798,0.000815,0.00085,0.0009,0.000898,0.000975,0.000993,0.000895,0.00078,0.00069,0.00061,0.00058,0.000582,0.00057,0.000563;0.000232,0.00026,0.000284,0.000332,0.000352,0.000348,0.000336,0.0003192,0.0003216,0.0003192,0.000326,0.00034,0.00036,0.0003592,0.00039,0.0003972,0.000358,0.000312,0.000276,0.000244,0.000232,0.0002328,0.000228,0.0002252;0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0];%节点无功;
Line=[1,1,2;2,2,3;3,3,4;4,4,5;5,5,6;6,6,7;7,7,8;8,8,9;9,9,10;10,10,11;11,11,12;12,12,13;13,13,14;14,14,15;15,15,16;16,16,17;17,5,25;18,25,26;19,26,27;20,27,28;21,28,29;22,29,30;23,30,31;24,31,32;25,1,18;26,18,19;27,19,20;28,20,21;29,2,22;30,22,23;31,23,24;32,33,1];%支路数据;
Pd=0.000018;%放电功率;
Pc=0.000018;%充电功率;
%Ndsum=[0;0;0;0;0;0;0;0;200;400;400;400;400;400;0;0;0;0;0;400;400;0;0;0];%上层调度放电电动汽车数量;
%Ncsum=[950;950;950;950;0;0;0;0;0;0;0;0;0;0;0;0;300;0;0;0;0;0;950;950];%上层调度充电电动汽车数量;
Ndsum=[0;0;0;0;0;0;0;0;160;160;160;160;160;160;0;0;0;0;0;160;80;0;0;0];%上层调度放电电动汽车数量(何立夫数据）;
Ncsum=[380;380;380;360;110;0;0;0;0;0;0;0;0;0;0;0;60;0;0;0;0;0;350;380];%上层调度充电电动汽车数量(何立夫数据）;
Nd_resid=zeros(24,1);Nc_resid=zeros(24,1);Nd_comme=zeros(24,1);Nc_comme=zeros(24,1);Nd_indus=zeros(24,1);Nc_indus=zeros(24,1);%区域电动汽车数量;