目录
[1 主要内容](#1 主要内容)
[2 部分代码](#2 部分代码)
[3 程序结果](#3 程序结果)
[4 下载链接](#4 下载链接)
1 主要内容
该程序复现《具有源荷不平衡特性的配电网智能软开关和储能联合规划》部分模型，未考虑聚类分析和分布鲁棒部分，就智能软开关和储能联合规划部分进行了方法复现，以IEEE33节点模型为例，实现了对原文中NY配电网配置，网络中分别接入3个光伏和风电节点，安装位置参照文献，以年综合费用最小为目标函数，建立协同规划模型，通过大M法和二阶锥松弛技术，将原始的非线性非凸模型转化为混合整数二阶锥规划模型，采用matlab+gurobi求解方式，储能和sop的配置数量可以设置，默认为分别配置两组，程序运行可靠，作为原文献中的核心非常具有参考价值！
程序的难点在于储能的配置约束表达，还是非常需要技巧的，既要避免非线性带来的求解难题，又要保证储能配置满足各项要求，程序非常精心处理该部分内容，非常具有参考意义！
目标函数
数据说明
程序中电价执行北京市峰谷电价，负荷和电价数据和文献一致，目标函数中部分参数根据经验值暂估设定，因有明确注释，方便修改调试。
节点系统图
2 部分代码
.rtcContent { padding: 30px; } .lineNode {font-size: 10pt; font-family: Menlo, Monaco, Consolas, "Courier New", monospace; font-style: normal; font-weight: normal; }u_data=sqrt(double(x_ui_square))/12.66;result=double(f);% disp(['程序总运行时间：',num2str(etime(clock,t1))]);%程序结果x_p_sop=value(x_p_sop);x_ui_square=value(x_ui_square);u_sop=value(u_sop);x_s_sop=value(x_s_sop);fi=value(fi);fo=value(fo);fsop=value(fsop);fl=value(fl);x_u=value(x_u);n_u=find(x_u>=0.9);x_p_cd=value(x_p_cd);% x_p_dis=value(x_p_dis);soc=value(soc);figure;plot(Pload,'b-*','LineWidth',1)ylabel('功率');yyaxis rightplot(Price,'LineWidth',1);ylabel('电价');ylim([2,3.6]);xlabel('时间');legend('负荷','电价')figure;xn=1:33;yn=1:N;mesh(yn,xn,u_data(:,:,1)),xlabel('时间'),ylabel('节点'),zlabel('电压（p.u.）')figure;bar(x_p_cd(n_u(1),:,1));ylabel('功率');yyaxis right;plot(soc(n_u(1),:,1),'LineWidth',1);ylim([0 1]);xlabel('时间');ylabel('soc');figure;bar(x_p_cd(n_u(2),:,1));ylabel('功率');yyaxis right;plot(soc(n_u(2),:,1),'LineWidth',1);ylim([0 1]);xlabel('时间');ylabel('soc');fprintf('sop选址为')c2=find(u_sop==1);sop_st(c2,:)fprintf('储能选址为')n_ufprintf('容量为%dMVA\n',x_s_sop(c2))fprintf('投资成本%d万元\n',fi)fprintf('运维成本%d万元\n',fo)fprintf('损耗成本%d万元\n',fsop)fprintf('弃风弃光%d万元\n',value(Cdis))fprintf('弃负荷%d万元\n',value(fp))fprintf('购电成本%d万元\n',value(Cpz))3 程序结果
4 下载链接
