主要内容
该模型为三微网电能交互共享模型，以运行成本和环境成本为目标，考虑负荷需求响应（可削减负荷和可转移负荷）、储能约束、风电约束和功率平衡约束等，利用ADMM算法进行迭代求解，实现联盟群效益的最大化，程序采用matlab+cplex求解，基本做到句句注释，具有较好的参考价值。
部分代码
%% ADMM算法参数设置 lambda_12=zeros(1,24);lambda_13=zeros(1,24);lambda_21=zeros(1,24); lambda_23=zeros(1,24);lambda_31=zeros(1,24);lambda_32=zeros(1,24);%拉格朗日乘子初始化 max_k=1000; %最大迭代次数 tolerant=1e-3;%收敛精度 rho=1e-4;%惩罚因子 k=1;%迭代次数初始化 tao=[]; %残差 P_12=zeros(max_k+1,24);P_21=zeros(max_k+1,24);P_13=zeros(max_k+1,24); P_31=zeros(max_k+1,24);P_23=zeros(max_k+1,24);P_32=zeros(max_k+1,24);%交互量 %% 迭代 while 1 if k==max_k break; end [P_wt_1,P_buy_1,P_batd_1,P_batc_1,P_sell_1,P_grid_1,P_bat_1,L_e_1,L_e0_1,E_co2_1,P_12(k+1,:),P_13(k+1,:),Obj_MG1(k)]=xin1(P_21(k+1,:),P_31(k+1,:),lambda_12,lambda_13,rho); [P_pv_2,P_buy_2,P_batd_2,P_batc_2,P_sell_2,P_grid_2,P_bat_2,L_e_2,L_e0_2,E_co2_2,P_21(k+1,:),P_23(k+1,:),Obj_MG2(k)]=xin2(P_12(k+1,:),P_32(k+1,:),lambda_21,lambda_23,rho); [P_pv_3,P_buy_3,P_batd_3,P_batc_3,P_sell_3,P_grid_3,P_bat_3,L_e_3,L_e0_3,E_co2_3,P_31(k+1,:),P_32(k+1,:),Obj_MG3(k)]=xin3(P_13(k+1,:),P_23(k+1,:),lambda_31,lambda_32,rho); lambda_12=lambda_12+rho*(P_12(k+1,:)+P_21(k+1,:)); lambda_13=lambda_13+rho*(P_13(k+1,:)+P_31(k+1,:)); lambda_21=lambda_21+rho*(P_21(k+1,:)+P_12(k+1,:)); lambda_23=lambda_23+rho*(P_23(k+1,:)+P_32(k+1,:)); lambda_31=lambda_31+rho*(P_31(k+1,:)+P_13(k+1,:)); lambda_32=lambda_32+rho*(P_32(k+1,:)+P_23(k+1,:)); tao=[tao,norm(P_12(k+1,:)-P_12(k,:))+norm(P_13(k+1,:)-P_13(k,:))+norm(P_23(k+1,:)-P_23(k,:))]; %残差计算 if tao(k)<=tolerant break; %判断收敛条件 end k=k+1; end toc disp(['微网1的二氧化碳排放量 : ', num2str(E_co2_1),' kg']); disp(['微网2的二氧化碳排放量 : ', num2str(E_co2_2),' kg']); disp(['微网3的二氧化碳排放量 : ', num2str(E_co2_3),' kg']); %% 画图 figure(1) plot(Obj_MG1,'m','LineWidth',1.5); hold on plot(Obj_MG2,'r','LineWidth',1.5); hold on plot(Obj_MG3,'b','LineWidth',1.5); hold on legend('微网1','微网2','微网3'); xlabel('迭代次数/k'); ylabel('成本/元'); legend('boxoff');
结果一览
