matlab代码：
参考文献：主动配电网多源协同运行优化研究_乔珊
摘要：最优潮流研究在配 电网规划运行 中不可或缺 ， 且在大量分布式能源接入 的主动配 电网环境下尤 为重要 。传统的启发式算法 在全局最优 解和求解 速度上均 无法满足主动配电网运行要求 ， 而基于线性化的最优潮流方法在高阻抗的配 电网中适用性也较 弱。 基于此，文章建立 了基 于二阶锥规划的动态最优潮流模型框架，力图将原非线性规划模型松弛转化为SOPC进行快速求解 。 首先 ， 给 出了基于二阶锥松弛 的配 电网动态最优潮流基本模 型； 然后,对主动配 电网中各重要参与元素进行相应 的线性化建模处理 以便 高效求解，如主动管理设备、配 电网重构、需求响应及综合负荷等；同时，分析 了松弛模型和近似等效 的准确性。
程序出图：
部分代码：
%多时段+SVC+CB+OLTC+DG SOCP_OPF Sbase=1MVA, Ubase=12.66KV
%目标函数如果只有网损，那么OLTC永远是高挡位，电压越高，网损越小，因此需进一步考虑目标函数如主网购电，或者电压平衡
%%
%有载调压变压器的位置在那个节点
%%
clear
clc
tic
warning off
%% 1.设参
mpc = IEEE33BW;
wind = mpc.wind;
pload = mpc.pload;
pload_prim = mpc.pload_prim/1000; %化为标幺值
qload_prim = mpc.qload_prim/1000;
a = 3.715; %单时段所有节点有功容量,MW
b = 2.3; %单时段所有节点无功容量,MW
pload = pload/a;%得到各个时段与单时段容量的比例系数
qload = pload/b;%假设有功负荷曲线与无功负荷变化曲线相同
pload = pload_prim*pload; %得到33*24的负荷值，每一个时间段每个节点的负荷
qload = qload_prim*qload;
branch = mpc.branch;
branch(:,3) = branch(:,3)*1/(12.66^2);%求阻抗标幺值
R = real(branch(:,3));
X = imag(branch(:,3));
T = 24;%时段数为24小时
nb = 33;%节点数
nl = 32;%支路数
nsvc = 3;%SVC数 静止无功补偿器 Static Var compensator
ncb = 2;%CB数 分组投切电容器组 (capacitorbanks，CB)
noltc = 1;%OLTC数 有载调压变压器 ( on—load tap changer，OLTC ） transformer
nwt = 2;%2个风机
ness = 2;%ESS数
upstream = zeros(nb,nl);
dnstream = zeros(nb,nl);
for i = 1:nl
upstream(i,i)=1;
end
for i = [1:16,18:20,22:23,25:31]
dnstream(i,i+1)=1;
end
