目录
1主要内容
上层目标函数考虑光伏和储能的投资成本。
程序采用模块化编程，并有每个模块功能介绍，方便学习。
2部分代码
3程序结果
4 程序结果
1主要内容
该程序主要方法复现《含高比例可再生能源配电网灵活资源双层优化配置》运行-规划联合双层配置模型，上层为光伏、储能选址定容模型，即优化配置，下层考虑弃光和储能出力，即优化调度，模型以IEEE33节点为例，采用粒子群算法求解，下层模型为运行成本和电压偏移量的多目标模型，并采用多目标粒子群算法得到pareto前沿解集，从中选择最佳结果带入到上层模型，最终实现上下层模型的各自求解和整个模型迭代优化。
上层目标函数考虑光伏和储能的投资成本。
parameter; %输入所有的数据 % 各个决策变量的含义 st_pv1=st_pvc(x(1));%光伏1选址cap_pv1=x(2);%光伏1容量st_cn=st_cnc(x(3));%储能选址cap_cn=x(4);%储能容量%目标%投资成本u=(r*(1+r)^year2)/((1+r)^year2-1);C1=u*(cpv*cap_pv1*s_pv+cess*cap_cn*s_cn);%光伏储能投资成本之和%带入参数调用运行层程序[gxbest1,fymin1]=lowerlayer(x);
下层目标函数考虑运行成本和电压偏移量双目标。
%运行成本=弃光成本+运维成本+购电成本+网损成本%弃光成本Cq=0;for i=1:K%不同场景Cq=Cq+cqpv*sum(gailv(i)*(cap_pv1.*center(i,:)-pv1_s(i)));end%运维成本Cy=cpvy*sum(gailv'.*(sum(pv1_s,2)+sum(pv1_s,2)))+cessy*sum(abs(cn));%购电成本 网损成本 电压偏移量Cgc=zeros(K,T);Clossc=zeros(K,T);Vpc=zeros(K,T);%matpower潮流计算
程序采用模块化编程，并有每个模块功能介绍，方便学习。
2部分代码
% ***************导入网络参数******************%parameter;st_pv1=st_pvc(upx(1));%光伏1选址cap_pv1=upx(2);%光伏1容量% st_pv2=st_pvc(upx(3));% cap_pv2=upx(4);st_cn=st_cnc(upx(3));%储能选址cap_cn=upx(4);%储能容量global center;gama=0.9;%光伏最低消纳率% *******************导入结束********************%%% ****************决策空间的设置1**************%maxFun=2;%2个目标函数fff=[0,800000;0,90]; %各个目标函数的最小值和最大值，即绝对正理想解和绝对负理想解%fff=[0,0.4;30,100;10000,20000]; n = 5;% 初始种群个数d = 120;% 空间维数，即决策变量（各设备控制量）个数maxIterations = 2; % 最大迭代次数 wmax=0.9;%maximum of inertia factor，最大惯性系数 wmin=0.4;%minimum of inertia factor，最小惯性系数c1=1.4962; %1.4962; %learning factor1，自我学习因子c2=1.4962; %1.4962; %learning factor2，群体学习因子Xmax=[cap_pv1.*center(1,:),cap_pv1.*center(2,:),cap_pv1.*center(3,:),cap_pv1.*center(4,:),s_cn.*cap_cn.*ones(1,24)];Xmin=[gama.*cap_pv1.*center(1,:),gama.*cap_pv1.*center(2,:),gama.*cap_pv1.*center(3,:),gama.*cap_pv1.*center(4,:),-s_cn.*cap_cn.*ones(1,24)];dX=Xmax-Xmin;Vmax=dX; Vmin=ones(1,120);%***********决策空间设置1结束**********%%% ******种群位置与速度初始化*******% X = repmat(Xmin,n,1)+repmat((Xmax-Xmin),n,1).*rand(n,d);%初始种群的位置 V=repmat(Vmin,n,1)+repmat((Vmax-Vmin),n,1).*rand(n,d);% 初始种群的速度%******变量维数（总共120维）：场景4*1个pv消纳值+储能出力
3程序结果
4 程序结果
点击直达！
