目錄

1?主要內容

目標函數

數據說明

節點系統圖

2?部分代碼

3?程序結果

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 下載鏈接?