抓取网站访客qq,网站制作公司业务发展方案,wordpress 无法升级,企业策划是什么意思本程序代码参考EI期刊论文《考虑电动汽车可调度潜力的充电站两阶段市场投标策略》#xff0c;程序中基于历史数据评估可调度潜力#xff0c;由联合报价模型确定节点边际电价#xff0c;作为报价的参考#xff0c;包含个体竞价模式#xff0c;纳什博弈竞价#xff0c;算例…本程序代码参考EI期刊论文《考虑电动汽车可调度潜力的充电站两阶段市场投标策略》程序中基于历史数据评估可调度潜力由联合报价模型确定节点边际电价作为报价的参考包含个体竞价模式纳什博弈竞价算例丰富注释清晰干货满满下面对程序代码和文章进行简要介绍
运行平台MATLABCplexYalmip或者MATLABGurobiYalmip
有很多公证号介绍了该程序代码差别从运行平台开始懂得都懂
使用两种求解器对模型进行求解均可以得到良好效果 在电力市场环境下,充电站优化投标策略能降低电力成本甚至通过售电获取收益。文中考虑了电动汽车成为柔性储荷资源的潜力提出了日前电力市场和实时电力市场下充电站的投标策略。首先基于闵可夫斯基加法提出了充电站内电动汽车集群模型的压缩方法并建立了日前可调度潜力预测模型和实时可调度潜力评估模型。同时考虑充电站间的非合作博弈建立了电力零售市场下充电站的策略投标模型并基于驻点法将其转化为一个广义Nash均衡问题。然后提出了基于日前报价和实时报量的两阶段市场交易模式并与合作投标模式、价格接受模式和集中调度模式进行对比。最后基于一个38节点配电系统进行了仿真。仿真结果表明所提出的可调度潜力计算方法能够将电动汽车集群封装为广义储能设备从而降低了模型的维度。基于可调度潜力的策略投标模型能够挖掘电动汽车的储荷潜力实现电动汽车与电网的有序互动。
算例共分为7节限于篇幅介绍前3节内容。
第一节算例结果 第二节算例结果 第三节算例结果 介绍到这里您看出和其他相关介绍的区别了吗
部分代码
% %% 第七节日前可调度潜力预测算法对比
% %算法1个体
% clear
% clc
% close all
% load data_EV
% X_CS1data_CS1.X;X_CS2data_CS2.X;X_CS3data_CS3.X;X_CS4data_CS4.X;
% S0_CS1data_CS1.S0;S0_CS2data_CS2.S0;S0_CS3data_CS3.S0;S0_CS4data_CS4.S0;
% m1length(S0_CS1);m2length(S0_CS2);m3length(S0_CS3);m4length(S0_CS4);
% %充电站投标问题
% pch_CS1sdpvar(m1,24);pdis_CS1sdpvar(m1,24);%电动汽车个体的充放电功率
% pch_CS2sdpvar(m2,24);pdis_CS2sdpvar(m2,24);
% pch_CS3sdpvar(m3,24);pdis_CS3sdpvar(m3,24);
% pch_CS4sdpvar(m4,24);pdis_CS4sdpvar(m4,24);
% s_CS1sdpvar(m1,24);s_CS2sdpvar(m2,24);s_CS3sdpvar(m3,24);s_CS4sdpvar(m4,24);%电动汽车个体的SOC
% Pchsdpvar(4,24);%充电站充电功率
% Pdissdpvar(4,24);%充电站放电功率
% Ccs[Pch(1,:)sum(pch_CS1),Pch(2,:)sum(pch_CS2),Pch(3,:)sum(pch_CS3),Pch(4,:)sum(pch_CS4),
% Pdis(1,:)sum(pdis_CS1),Pdis(2,:)sum(pdis_CS2),Pdis(3,:)sum(pdis_CS3),Pdis(4,:)sum(pdis_CS4),
% 0pch_CS16.6*X_CS1,0pch_CS26.6*X_CS2,0pch_CS36.6*X_CS3,0pch_CS46.6*X_CS4,
% 0pdis_CS16.6*X_CS1,0pdis_CS26.6*X_CS2,0pdis_CS36.6*X_CS3,0pdis_CS46.6*X_CS4,
% s_CS1(:,1)S0_CS10.95*pch_CS1(:,1)-pdis_CS1(:,1)/0.95,
% s_CS2(:,1)S0_CS20.95*pch_CS2(:,1)-pdis_CS2(:,1)/0.95,
% s_CS3(:,1)S0_CS30.95*pch_CS3(:,1)-pdis_CS3(:,1)/0.95,
% s_CS4(:,1)S0_CS40.95*pch_CS4(:,1)-pdis_CS4(:,1)/0.95,
% s_CS1(:,2:24)s_CS1(:,1:23)0.95*pch_CS1(:,2:24)-pdis_CS1(:,2:24)/0.95,
% s_CS2(:,2:24)s_CS2(:,1:23)0.95*pch_CS2(:,2:24)-pdis_CS2(:,2:24)/0.95,
% s_CS3(:,2:24)s_CS3(:,1:23)0.95*pch_CS3(:,2:24)-pdis_CS3(:,2:24)/0.95,
% s_CS4(:,2:24)s_CS4(:,1:23)0.95*pch_CS4(:,2:24)-pdis_CS4(:,2:24)/0.95,
% 32*0.15s_CS132*0.9,32*0.15s_CS232*0.9,32*0.15s_CS332*0.9,32*0.15s_CS432*0.9
% sum(0.95*pch_CS1-pdis_CS1/0.95,2)32*0.9-S0_CS1,
% sum(0.95*pch_CS2-pdis_CS2/0.95,2)32*0.9-S0_CS2,
% sum(0.95*pch_CS3-pdis_CS3/0.95,2)32*0.9-S0_CS3,
% sum(0.95*pch_CS4-pdis_CS4/0.95,2)32*0.9-S0_CS4
% ];%约束条件
% % %市场出清问题
% Loadcurve[0.955391944564747,0.978345604157644,1,0.995019488956258,0.972932005197055,0.970333477695972,0.930489389346037,0.890428757037679,0.902771762667822,0.941966219142486,0.911000433087917,0.862061498484192,0.840190558683413,0.831095712429623,0.756604590731919,0.671719359029883,0.611520138588133,0.582936336076224,0.572542226071893,0.574707665656128,0.587267215244695,0.644218276310091,0.755521870939801,0.884798614118666];
% PL_base[5.704;5.705;5.631;6.518;4.890;5.705;5.847]*1000;%负荷分布
% PLPL_base*Loadcurve;%基础负荷
% Pfsdpvar(7,24);%馈线功率
% Pf(1,:)PL(1,:)Pch(1,:)-Pdis(1,:);Pf(2,:)PL(2,:);Pf(3,:)PL(3,:);Pf(4,:)PL(4,:)Pch(2,:)-Pdis(2,:);Pf(5,:)PL(5,:)Pch(3,:)-Pdis(3,:);Pf(6,:)PL(6,:);Pf(7,:)PL(7,:)Pch(4,:)-Pdis(4,:);%馈线功率组成
% Pf_limit1000*[40,40,40,40,40,40,40];%馈线功率限制
% Pgsdpvar(10,24);%发电商分段电量
% Pg_step1000*[20,5,3,2,2,2,2,2,2,100];%报价区间
% Price_DSO[3:12]*0.1;%分段电价
% Lagrant_G_leftsdpvar(10,24);%发电商电量下界
% Lagrant_G_rightsdpvar(10,24);%发电商电量上界
% b_Lagrant_G_leftbinvar(10,24);%发电商电量下界布尔变量
% b_Lagrant_G_rightbinvar(10,24);%发电商电量上界布尔变量
% Lagrant_L_leftsdpvar(7,24);%线路功率下界
% Lagrant_L_rightsdpvar(7,24);%线路功率上界
% b_Lagrant_L_leftbinvar(7,24);%线路功率上界布尔变量
% b_Lagrant_L_rightbinvar(7,24);%线路功率下界布尔变量
% Lagrant_Gsdpvar(1,24);%平衡节点电价
% DLMPsdpvar(7,24);%DLMP
% Ckkt[sum(Pg)sum(Pf),
% 0PgPg_step*ones(1,24),
% Price_DSO*ones(1,24)-Lagrant_G_leftLagrant_G_right-ones(10,1)*Lagrant_G0,
% DLMPones(7,1)*Lagrant_GLagrant_L_right-Lagrant_L_left,
% Pg1E6*b_Lagrant_G_left,0Lagrant_G_left(1-b_Lagrant_G_left),
% Pg_step*ones(1,24)-Pg1E6*b_Lagrant_G_right,0Lagrant_G_right(1-b_Lagrant_G_right),
% -Pf_limit*ones(1,24)PfPf_limit*ones(1,24),
% PfPf_limit*ones(1,24)1E6*b_Lagrant_L_left,0Lagrant_L_left(1-b_Lagrant_L_left),
% -PfPf_limit*ones(1,24)1E6*b_Lagrant_L_right,0Lagrant_L_right(1-b_Lagrant_L_right)];%市场出清问题的KKT条件
% Objsum(sum((Price_DSO*ones(1,24)).*Pg))-sum(sum(DLMP.*PL))sum(sum((Pg_step*ones(1,24)).*Lagrant_G_right))sum(sum((Lagrant_L_leftLagrant_L_right).*(Pf_limit*ones(1,24))));
% %求解问题
% C[Ccs,Ckkt];%约束条件
% opssdpsettings(solver,cplex);
% resultsolvesdp(C,Obj,ops);
% DLMPdouble(DLMP);
% Lagrant_Gdouble(Lagrant_G);
% Pfdouble(Pf);
% Pchdouble(Pch);
% Pdisdouble(Pdis);
% Pgdouble(Pg);
% %实时调度
% load data_RT_potential
% load result_order
% Pch1RT_CS1(96).Pch;Pdis1RT_CS1(96).Pdis;Smin1RT_CS1(96).Smin;Smax1RT_CS1(96).Smax;dS1RT_CS1(96).dS;%取最后一个时刻的可调度潜力
% Pch2RT_CS2(96).Pch;Pdis2RT_CS2(96).Pdis;Smin2RT_CS2(96).Smin;Smax2RT_CS2(96).Smax;dS2RT_CS2(96).dS;
% Pch3RT_CS3(96).Pch;Pdis3RT_CS3(96).Pdis;Smin3RT_CS3(96).Smin;Smax3RT_CS3(96).Smax;dS3RT_CS3(96).dS;
% Pch4RT_CS4(96).Pch;Pdis4RT_CS4(96).Pdis;Smin4RT_CS4(96).Smin;Smax4RT_CS4(96).Smax;dS4RT_CS4(96).dS;
% Pchmax[Pch1;Pch2;Pch3;Pch4];
% Pdismax[Pdis1;Pdis2;Pdis3;Pdis4];
% Smin[Smin1;Smin2;Smin3;Smin4];
% Smax[Smax1;Smax2;Smax3;Smax4];
% dS[dS1;dS2;dS3;dS4];
% Linkzeros(24,96);%时段换算矩阵(日前1h换算为实时15min)
% for i1:24
% Link(i,4*i-3:4*i)1;
% end
% DA_PsPch-Pdis;%日前调度数据
% DA_PsDA_Ps*Link;%时段换算
% DA_PfPf*Link;%日前总功率
% %充电站投标问题
% Pchsdpvar(4,96);%充电站充电功率
% Pdissdpvar(4,96);%充电站放电功率
% Ssdpvar(4,96);%充电站广义储能SOC
% Ccs[0PchPchmax,0PdisPdismax,SminSSmax,S(:,1)0.95*0.25*Pch(:,1)-0.25*Pdis(:,1)/0.95dS(:,1),S(:,2:96)S(:,1:95)0.25*0.95*Pch(:,2:96)-0.25*Pdis(:,2:96)/0.95dS(:,2:96)];%约束条件
% % %市场出清问题
% Loadcurve[0.955391944564747,0.978345604157644,1,0.995019488956258,0.972932005197055,0.970333477695972,0.930489389346037,0.890428757037679,0.902771762667822,0.941966219142486,0.911000433087917,0.862061498484192,0.840190558683413,0.831095712429623,0.756604590731919,0.671719359029883,0.611520138588133,0.582936336076224,0.572542226071893,0.574707665656128,0.587267215244695,0.644218276310091,0.755521870939801,0.884798614118666];
% LoadcurveLoadcurve*Link;%换成96个时段
% PL_base[5.704;5.705;5.631;6.518;4.890;5.705;5.847]*1000;%负荷分布
% PLPL_base*Loadcurve;%基础负荷
% Pfsdpvar(7,96);%馈线功率
% Pf(1,:)PL(1,:)Pch(1,:)-Pdis(1,:);Pf(2,:)PL(2,:);Pf(3,:)PL(3,:);Pf(4,:)PL(4,:)Pch(2,:)-Pdis(2,:);Pf(5,:)PL(5,:)Pch(3,:)-Pdis(3,:);Pf(6,:)PL(6,:);Pf(7,:)PL(7,:)Pch(4,:)-Pdis(4,:);%馈线功率组成
% Pf_limit1000*[40,40,40,40,40,40,40];%馈线功率限制
% Pgsdpvar(10,96);%发电商分段电量
% Pg_step1000*[20,5,3,2,2,2,2,2,2,100];%报价区间
% Price_DSO[3:12]*0.1;%分段电价
% Lagrant_G_leftsdpvar(10,96);%发电商电量下界
% Lagrant_G_rightsdpvar(10,96);%发电商电量上界
% b_Lagrant_G_leftbinvar(10,96);%发电商电量下界布尔变量
% b_Lagrant_G_rightbinvar(10,96);%发电商电量上界布尔变量
% Lagrant_L_leftsdpvar(7,96);%线路功率下界
% Lagrant_L_rightsdpvar(7,96);%线路功率上界
% b_Lagrant_L_leftbinvar(7,96);%线路功率上界布尔变量
% b_Lagrant_L_rightbinvar(7,96);%线路功率下界布尔变量
% Lagrant_Gsdpvar(1,96);%平衡节点电价
% DLMPsdpvar(7,96);%DLMP
% Ckkt[sum(Pg)sum(Pf),
% 0PgPg_step*ones(1,96),
% 0.25*Price_DSO*ones(1,96)-Lagrant_G_leftLagrant_G_right-ones(10,1)*Lagrant_G0,
% DLMPones(7,1)*Lagrant_GLagrant_L_right-Lagrant_L_left,
% Pg1E6*b_Lagrant_G_left,0Lagrant_G_left(1-b_Lagrant_G_left),
% Pg_step*ones(1,96)-Pg1E6*b_Lagrant_G_right,0Lagrant_G_right(1-b_Lagrant_G_right),
% -Pf_limit*ones(1,96)PfPf_limit*ones(1,96),
% PfPf_limit*ones(1,96)1E6*b_Lagrant_L_left,0Lagrant_L_left(1-b_Lagrant_L_left),
% -PfPf_limit*ones(1,96)1E6*b_Lagrant_L_right,0Lagrant_L_right(1-b_Lagrant_L_right)];%市场出清问题的KKT条件
% Obj0.25*sum(sum((Price_DSO*ones(1,96)).*Pg))-sum(sum(DLMP.*PL))sum(sum((Pg_step*ones(1,96)).*Lagrant_G_right))sum(sum((Lagrant_L_leftLagrant_L_right).*(Pf_limit*ones(1,96))))0.1*0.25*sum(sum(abs(Pch-Pdis-DA_Ps)));
% %求解问题
% C[Ccs,Ckkt];%约束条件
% resultsolvesdp(C,Obj);
% DLMPdouble(DLMP);
% Lagrant_Gdouble(Lagrant_G);
% Pfdouble(Pf);
% Pchdouble(Pch);
% Pdisdouble(Pdis);
% Pgdouble(Pg);
% Sdouble(S);
% double(Obj)
% figure(1)%实时日前结果对比
% hold on
% plot(sum(Pf),r.-)
% plot(sum(DA_Pf),b--);
% Costsum(sum((Pf-PL).*DLMP));%总用电费用
% Cost_bsum(sum(((Pf-PL)0).*(Pf-PL).*DLMP))%总购电费用
% Cost_s-sum(sum(((Pf-PL)0).*(Pf-PL).*DLMP))%总售电收入
% Cost_r0.1*0.25*sum(sum(abs(Pch-Pdis-DA_Ps)))%平衡费用
% MAPE100*sum(abs((sum(DA_Pf)-sum(Pf))./sum(DA_Pf))/96);%平均百分误差
% %% 第七节日前可调度潜力预测算法对比
% %算法2群体
% clear
% clc
% close all
% load data_EV
% X_CS1data_CS1.X;X_CS2data_CS2.X;X_CS3data_CS3.X;X_CS4data_CS4.X;
% S0_CS1data_CS1.S0;S0_CS2data_CS2.S0;S0_CS3data_CS3.S0;S0_CS4data_CS4.S0;
% m1length(S0_CS1);m2length(S0_CS2);m3length(S0_CS3);m4length(S0_CS4);
% X[sum(X_CS1);sum(X_CS2);sum(X_CS3);sum(X_CS4)];
% Energy[32*0.9*m1-sum(S0_CS1);32*0.9*m2-sum(S0_CS2);32*0.9*m3-sum(S0_CS3);32*0.9*m4-sum(S0_CS4)];
% %充电站投标问题
% Pchsdpvar(4,24);%充电站充电功率
% Pdissdpvar(4,24);%充电站放电功率
% Ccs[0Pch6.6*X,
% 0Pdis6.6*X,
% sum(Pch-Pdis,2)Energy
% ];%约束条件
% % %市场出清问题
% Loadcurve[0.955391944564747,0.978345604157644,1,0.995019488956258,0.972932005197055,0.970333477695972,0.930489389346037,0.890428757037679,0.902771762667822,0.941966219142486,0.911000433087917,0.862061498484192,0.840190558683413,0.831095712429623,0.756604590731919,0.671719359029883,0.611520138588133,0.582936336076224,0.572542226071893,0.574707665656128,0.587267215244695,0.644218276310091,0.755521870939801,0.884798614118666];
% PL_base[5.704;5.705;5.631;6.518;4.890;5.705;5.847]*1000;%负荷分布
% PLPL_base*Loadcurve;%基础负荷
% Pfsdpvar(7,24);%馈线功率
% Pf(1,:)PL(1,:)Pch(1,:)-Pdis(1,:);Pf(2,:)PL(2,:);Pf(3,:)PL(3,:);Pf(4,:)PL(4,:)Pch(2,:)-Pdis(2,:);Pf(5,:)PL(5,:)Pch(3,:)-Pdis(3,:);Pf(6,:)PL(6,:);Pf(7,:)PL(7,:)Pch(4,:)-Pdis(4,:);%馈线功率组成
% Pf_limit1000*[40,40,40,40,40,40,40];%馈线功率限制
% Pgsdpvar(10,24);%发电商分段电量
% Pg_step1000*[20,5,3,2,2,2,2,2,2,100];%报价区间
% Price_DSO[3:12]*0.1;%分段电价
% Lagrant_G_leftsdpvar(10,24);%发电商电量下界
% Lagrant_G_rightsdpvar(10,24);%发电商电量上界
% b_Lagrant_G_leftbinvar(10,24);%发电商电量下界布尔变量
% b_Lagrant_G_rightbinvar(10,24);%发电商电量上界布尔变量
% Lagrant_L_leftsdpvar(7,24);%线路功率下界
% Lagrant_L_rightsdpvar(7,24);%线路功率上界
% b_Lagrant_L_leftbinvar(7,24);%线路功率上界布尔变量
% b_Lagrant_L_rightbinvar(7,24);%线路功率下界布尔变量
% Lagrant_Gsdpvar(1,24);%平衡节点电价
% DLMPsdpvar(7,24);%DLMP
% Ckkt[sum(Pg)sum(Pf),
% 0PgPg_step*ones(1,24),
% Price_DSO*ones(1,24)-Lagrant_G_leftLagrant_G_right-ones(10,1)*Lagrant_G0,
% DLMPones(7,1)*Lagrant_GLagrant_L_right-Lagrant_L_left,
% Pg1E6*b_Lagrant_G_left,0Lagrant_G_left(1-b_Lagrant_G_left),
% Pg_step*ones(1,24)-Pg1E6*b_Lagrant_G_right,0Lagrant_G_right(1-b_Lagrant_G_right),
% -Pf_limit*ones(1,24)PfPf_limit*ones(1,24),
% PfPf_limit*ones(1,24)1E6*b_Lagrant_L_left,0Lagrant_L_left(1-b_Lagrant_L_left),
% -PfPf_limit*ones(1,24)1E6*b_Lagrant_L_right,0Lagrant_L_right(1-b_Lagrant_L_right)];%市场出清问题的KKT条件
% Objsum(sum((Price_DSO*ones(1,24)).*Pg))-sum(sum(DLMP.*PL))sum(sum((Pg_step*ones(1,24)).*Lagrant_G_right))sum(sum((Lagrant_L_leftLagrant_L_right).*(Pf_limit*ones(1,24))));
% %求解问题
% C[Ccs,Ckkt];%约束条件
% opssdpsettings(solver,cplex);
% resultsolvesdp(C,Obj,ops);
% DLMPdouble(DLMP);
% Lagrant_Gdouble(Lagrant_G);
% Pfdouble(Pf);
% Pchdouble(Pch);
% Pdisdouble(Pdis);
% Pgdouble(Pg);
% %实时调度
% load data_RT_potential
% load result_order
% Pch1RT_CS1(96).Pch;Pdis1RT_CS1(96).Pdis;Smin1RT_CS1(96).Smin;Smax1RT_CS1(96).Smax;dS1RT_CS1(96).dS;%取最后一个时刻的可调度潜力
% Pch2RT_CS2(96).Pch;Pdis2RT_CS2(96).Pdis;Smin2RT_CS2(96).Smin;Smax2RT_CS2(96).Smax;dS2RT_CS2(96).dS;
% Pch3RT_CS3(96).Pch;Pdis3RT_CS3(96).Pdis;Smin3RT_CS3(96).Smin;Smax3RT_CS3(96).Smax;dS3RT_CS3(96).dS;
% Pch4RT_CS4(96).Pch;Pdis4RT_CS4(96).Pdis;Smin4RT_CS4(96).Smin;Smax4RT_CS4(96).Smax;dS4RT_CS4(96).dS;
% Pchmax[Pch1;Pch2;Pch3;Pch4];
% Pdismax[Pdis1;Pdis2;Pdis3;Pdis4];
% Smin[Smin1;Smin2;Smin3;Smin4];
% Smax[Smax1;Smax2;Smax3;Smax4];
% dS[dS1;dS2;dS3;dS4];
% Linkzeros(24,96);%时段换算矩阵(日前1h换算为实时15min)
% for i1:24
% Link(i,4*i-3:4*i)1;
% end
% DA_PsPch-Pdis;%日前调度数据
% DA_PsDA_Ps*Link;%时段换算
% DA_PfPf*Link;%日前总功率
% %充电站投标问题
% Pchsdpvar(4,96);%充电站充电功率
% Pdissdpvar(4,96);%充电站放电功率
% Ssdpvar(4,96);%充电站广义储能SOC
% Ccs[0PchPchmax,0PdisPdismax,SminSSmax,S(:,1)0.95*0.25*Pch(:,1)-0.25*Pdis(:,1)/0.95dS(:,1),S(:,2:96)S(:,1:95)0.25*0.95*Pch(:,2:96)-0.25*Pdis(:,2:96)/0.95dS(:,2:96)];%约束条件
% % %市场出清问题
% Loadcurve[0.955391944564747,0.978345604157644,1,0.995019488956258,0.972932005197055,0.970333477695972,0.930489389346037,0.890428757037679,0.902771762667822,0.941966219142486,0.911000433087917,0.862061498484192,0.840190558683413,0.831095712429623,0.756604590731919,0.671719359029883,0.611520138588133,0.582936336076224,0.572542226071893,0.574707665656128,0.587267215244695,0.644218276310091,0.755521870939801,0.884798614118666];
% LoadcurveLoadcurve*Link;%换成96个时段
% PL_base[5.704;5.705;5.631;6.518;4.890;5.705;5.847]*1000;%负荷分布
% PLPL_base*Loadcurve;%基础负荷
% Pfsdpvar(7,96);%馈线功率
% Pf(1,:)PL(1,:)Pch(1,:)-Pdis(1,:);Pf(2,:)PL(2,:);Pf(3,:)PL(3,:);Pf(4,:)PL(4,:)Pch(2,:)-Pdis(2,:);Pf(5,:)PL(5,:)Pch(3,:)-Pdis(3,:);Pf(6,:)PL(6,:);Pf(7,:)PL(7,:)Pch(4,:)-Pdis(4,:);%馈线功率组成
% Pf_limit1000*[40,40,40,40,40,40,40];%馈线功率限制
% Pgsdpvar(10,96);%发电商分段电量
% Pg_step1000*[20,5,3,2,2,2,2,2,2,100];%报价区间
% Price_DSO[3:12]*0.1;%分段电价
% Lagrant_G_leftsdpvar(10,96);%发电商电量下界
% Lagrant_G_rightsdpvar(10,96);%发电商电量上界
% b_Lagrant_G_leftbinvar(10,96);%发电商电量下界布尔变量
% b_Lagrant_G_rightbinvar(10,96);%发电商电量上界布尔变量
% Lagrant_L_leftsdpvar(7,96);%线路功率下界
% Lagrant_L_rightsdpvar(7,96);%线路功率上界
% b_Lagrant_L_leftbinvar(7,96);%线路功率上界布尔变量
% b_Lagrant_L_rightbinvar(7,96);%线路功率下界布尔变量
% Lagrant_Gsdpvar(1,96);%平衡节点电价
% DLMPsdpvar(7,96);%DLMP
% Ckkt[sum(Pg)sum(Pf),
% 0PgPg_step*ones(1,96),
% 0.25*Price_DSO*ones(1,96)-Lagrant_G_leftLagrant_G_right-ones(10,1)*Lagrant_G0,
% DLMPones(7,1)*Lagrant_GLagrant_L_right-Lagrant_L_left,
% Pg1E6*b_Lagrant_G_left,0Lagrant_G_left(1-b_Lagrant_G_left),
% Pg_step*ones(1,96)-Pg1E6*b_Lagrant_G_right,0Lagrant_G_right(1-b_Lagrant_G_right),
% -Pf_limit*ones(1,96)PfPf_limit*ones(1,96),
% PfPf_limit*ones(1,96)1E6*b_Lagrant_L_left,0Lagrant_L_left(1-b_Lagrant_L_left),
% -PfPf_limit*ones(1,96)1E6*b_Lagrant_L_right,0Lagrant_L_right(1-b_Lagrant_L_right)];%市场出清问题的KKT条件
% Obj0.25*sum(sum((Price_DSO*ones(1,96)).*Pg))-sum(sum(DLMP.*PL))sum(sum((Pg_step*ones(1,96)).*Lagrant_G_right))sum(sum((Lagrant_L_leftLagrant_L_right).*(Pf_limit*ones(1,96))))0.1*0.25*sum(sum(abs(Pch-Pdis-DA_Ps)));
% %求解问题
% C[Ccs,Ckkt];%约束条件
% resultsolvesdp(C,Obj);
% DLMPdouble(DLMP);
% Lagrant_Gdouble(Lagrant_G);
% Pfdouble(Pf);
% Pchdouble(Pch);
% Pdisdouble(Pdis);
% Pgdouble(Pg);
% Sdouble(S);
% double(Obj)
% figure(1)%实时日前结果对比
% hold on
% plot(sum(Pf),r.-)
% plot(sum(DA_Pf),b--);
% Costsum(sum((Pf-PL).*DLMP));%总用电费用
% Cost_bsum(sum(((Pf-PL)0).*(Pf-PL).*DLMP))%总购电费用
% Cost_s-sum(sum(((Pf-PL)0).*(Pf-PL).*DLMP))%总售电收入
% Cost_r0.1*0.25*sum(sum(abs(Pch-Pdis-DA_Ps)))%平衡费用
% MAPE100*sum(abs((sum(DA_Pf)-sum(Pf))./sum(DA_Pf))/96);%平均百分误差以上就是本次介绍的主要内容欢迎关注下方公众号获取完整版代码小编会继续推送更有质量的学习资料、文章和程序代码为您的科研加油助力
