基于二阶锥规划的主动配电网最优潮流求解（CPLEX二阶锥规划考虑Wind+CB+SVG+OLTC+ESS多时段24h）

function mpc = IEEE33BW %% 初始各节点有功负荷容量 mpc.pload_prim = [100;90;120;60;60;200;200;60;60;45;60;60;120;60;60;60;90;90;90;90;90;90;420;420;60;60;60;120;200;150;210;60;0]; mpc.qload_prim = [60;40;80;30;20;100;100;20;20;30;35;35;80;10;20;20;40;40;40;40;40;50;200;200;25;25;20;70;600;70;100;40;0]; %这里的第33号节点的初始值PQ为0，因为它是平衡节点。 %% 风机总出力(24小时) mpc.wind = [0.9 1.0 0.95 0.9 0.8 0.79 0.76 0.75 0.55 0.35 0.45 0.43 0.35 0.32 0.30 0.33 0.35 0.43 0.55 0.72 0.75 0.79 0.82 0.95]; % plot(mpc.wind) %% 有功负荷需求（这里可能指的是风电场本地有功负荷需求，但数值不应该比发电还大） mpc.pload = [1.8 1.35 1.2 1.15 1.1 1.35 1.45 1.5 1.65 1.9 2.0 2.1 2.35 2.55 2.25 2.15 2.3 2.6 3.35 3.85 4.0 3.6 3.25 2.7]; % plot(mpc.pload) %% 节点负荷（标幺值,24小时值）(33*24) mpc.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]; %33*24 mpc.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.