55号资源-源程序：论文可在知网下载《计及电动汽车和需求响应的多类电力市场下虚拟电厂竞标模型》本人博客有解读

%% 需求响应包括可转移负荷、可中断负荷 %% 经济成本最低，改善负荷曲线 %% 程序初始化 clc clear close all %% 定义变量 xs=1.05*[100,100,100,100,100,100,380,380,380,800,800,800,800,800,380,380,380,800,800,800,380,380,100,100]; %% 需求响应部分参数 kil=[150,200,400];%中断负荷补偿费用 cil=[0.1,0.08,0.05];%最大中断负荷比例 %% 负荷值-有更新 pload=[4.2 4.35 4.8 4.2 3.6 3.15 4.95 6.3 6.15 18.8 19.6 19.2 18.4 17.6 16.5 15.6 14.6 13.2 12.9 9.3 4.95 5.4 4.2 3.6]; %% 定义变量sdpvar/binvar pil=sdpvar(3,24);%中断负荷 cshift=sdpvar(1,24);%定义转移负荷补偿费用 shiftl=sdpvar(1,24);%转移到别处的负荷 shiftq=sdpvar(1,24);%从别处转移过来的负荷 ushiftl=binvar(1,24);%0-1变量 ushiftq=binvar(1,24);%0-1变量 pmgb=sdpvar(1,24);%用户用电费用 cxuefeng=sdpvar(1,24);%削峰填谷收益 %% 约束条件 C=[];%初始化约束 %% 需求响应资源——可中断负荷约束 for m=1:3 for t=1:24 C=[C, 0<=pil(m,t)<=cil(m)*pload(t), %各级中断负荷约束 ]; end end for m=1:3 for t=2:24 C=[C, pil(m,t)+ pil(m,t-1)<=0.2*pload(t), %连续性中断负荷约束 ] ; end end %% 可转移负荷约束 for t=1:24 C=[C, ushiftl(t)+ushiftq(t)<=1;%设置互斥约束，不能同时转移和转接 shiftl(t)<= ushiftl(t)*pload(t)*0.12;%每个时刻最大转移量 shiftq(t)<=ushiftq(t)*pload(t)*0.12;%每个时刻最大转接量 0<=shiftl(t); 0<=shiftq(t); ]; end C=[C,sum(shiftl)==sum(shiftq)];%转移前后总负荷应当保持不变 %% 功率平衡约束 for t=1:24 C=[C, pload(t)-shiftl(t)+shiftq(t)-sum(pil(:,t))==pmgb(t), ]; end for t=1:24 C=[C, cshift(t)==50*shiftl(t)+30*shiftq(t) %计算转移负荷的补偿费用 ]; end %% 费用计算 F1=sum(xs.*pmgb) F2=sum(pil(1,:)*kil(1))+sum(pil(2,:)*kil(2))+sum(pil(3,:)*kil(3))%可中断负荷补偿费用 F3=sum(cshift)%可转移负荷补偿费用 F=F1+F2+F3 %总目标函数 %% 求解器设置 ops=sdpsettings('solver','cplex','verbose',2,'usex0',0); ops.cplex.mip.tolerances.mipgap=1e-6; result=optimize(C,F,ops); if result.problem == 0 % problem =0 代表求解成功 else error('求解出错'); end %% 结果数值读取 F1=value(F1) F2=value(F2) F3=value(F3) shiftl=value(shiftl); shiftq=value(shiftq); pil=value(pil); % 可转移负荷转移特性 figure plot(pload-shiftl+shiftq,'r-*','LineWidth',1.5) hold on plot(pload,'g--*','LineWidth',1.5) xlabel('时间/h'); ylabel('负荷功率/MW'); title('可转移负荷响应结果'); legend('转移后的负荷曲线','转移前的负荷曲线'); legend('boxoff'); box off % 转移负荷结果 figure bar(shiftl-shiftq,'stacked') xlabel('时间/h'); ylabel('转移的负荷功率/MW'); yyaxis right plot(xs,'r-o','LineWidth',1.5); hold on ylabel('购电电价/（元/MWh）'); title('可转移负荷响应结果'); legend('负荷转移值','分时电价值'); legend('boxoff'); box off % 总负荷变化曲线 figure plot(pload,'g-*','LineWidth',1.5) hold on xlabel('时间/h'); plot(pload-shiftl+shiftq-sum(pil(:,:)),'c-*','LineWidth',1.5) ylabel('负荷功率/MW'); hold on yyaxis right plot(xs,'r-.','LineWidth',1.5); ylabel('分时电价/元') title('可转移负荷响应结果'); legend('原负荷曲线','最终负荷曲线','分时电价'); legend('boxoff'); box off %% 可中断负荷响应结果 % 8、中断负荷调度结果 figure bar(pil',0.6,'stacked'); hold on ylabel('功率/MW'); yyaxis right plot(xs,'r--*','LineWidth',1.5); hold on ylabel('购电电价/（元/MWh）'); xlabel('时间/h'); title('可中断负荷调度结果'); legend('一级负荷','二级负荷','三级负荷','市场电价'); legend('boxoff'); box off