65号资源-源程序：论文可在知网下载《面向综合能源园区的三方市场主体非合作交易方法》本人博客有解读

function fobjplot(x)%本篇代码的目的是响应main函数中的fobjplot(zbest); T=24; pz=reshape(x,T,[]);%假如x是一个m*n的矩阵，则pz为T*（（m*n）/T）的矩阵，矩阵的内容不会改变 peopr=pz(:,1)*0.6+0.5;%主导商向用户、EV的电费报价 phopr=pz(:,2)*0.2+0.5;%主导商向用户的热费报价 peusr=pz(:,3)*0.4+0.6;%用户向EV的电费报价 Wbest=ones(1,T)*0.5;%建立一个1行24列初值均为0.5的矩阵（这里的意思是刚开始假设充电站的电能一半由主导商供应，另一半由用户供应，w为供应比例） for ic=1:T Wseq=0:0.025:1;%Wseq是一个循循渐进的1行40列的矩阵，随着迭代次数的增加，矩阵不断变长，这里是逐渐改变EV的供电占比，得到最优解 for ii=1:length(Wseq)%ii的值从1到40不断增加 W=Wbest;W(ic)=Wseq(ii);%w为用户或主导商向EV的供电占比 outz(ii)=fobjsub(peopr,phopr,peusr,W,ic);%在fobjsub函数中对outz(1行41列)进行求解（主导商向用户、EV报电价，主导商向用户报热价，用户向EV报电价，EV策略合集） end [~,Ind]=max(outz);%Ind为一个数值，代表了outz矩阵中最大值所在的列数 Wbest(ic)=Wseq(Ind);%wbest的值随着外循环的进行而不断改变，进而不断改变W的值 end %基础数值的图像输出 T=24;%设置初始值 Pepvusr=[0 0 0 0 0 0 10 40 200 600 700 750 850 900 840 780 650 380 100 50 10 0 0 0];%光伏出力 Peacusr=[150 160 170 180 190 200 220 400 550 600 600 550 580 540 570 550 450 400 300 220 200 190 180 170];%电负荷 Phopr=[140 140 160 170 180 190 400 700 950 900 850 800 850 900 850 800 500 350 280 210 190 180 170 160];%热负荷（运营商向用户供热功率） pgrid=[0.4 0.4 0.4 0.4 0.4 0.4 0.7 0.7 1.25 1.25 1.25 0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.7 1.25 1.25 1.25 0.7 0.4];%分时电价表 pgas=2.92;%天然气价格 ceopr=0.002;%配电变压器维护成本 ceoprchp=0.07;%CHP机组的维护成本（电） cgbhopr=0.002;%燃气锅炉维护成本 choprchp=0.07;%CHP机组的维护成本（热） palow=0.3;%光伏发电补贴 cepvusr=0.006;%光伏发电板维护成本 ceacusr=0.002;%制热空调维护成本 Cchs=250*T*2; Pech=ones(1,T)*250;%Pech为1行24列初值为250的矩阵 Pechopr=W.*Pech;%电网给充电站供能 注意这个W!,W影响到了下面所有的参数！ Pechusr=(1-W).*Pech;%用户给充电站供能（供应车辆充电负荷） Peuopr=max(Peacusr-Pepvusr+Pechusr,0);%电网供应用户电负荷 Pgrid = min(Peuopr+Pechopr,500);%电网供电功率 Peopr = Peuopr+Pechopr;%变压器输出功率（运营商向园区内部供能功率） Pechpopr = max(Peopr-Pgrid,0);%CHP机组供电大小 Phchpopr=Pechpopr*1.3;%热电比为1.3，这里是CHP机组供热大小 Phgbopr=min(max(Phopr-Phchpopr,0),500);%燃气供给功率（燃气热输出） Qgas = Phgbopr/10;%热网供热功率？？？ Phusr=Phopr-Phgbopr;%电空调热能（热负荷减去燃气锅炉） figure;hold on; plot(Pgrid,'-x'); plot(Phgbopr,'-o'); legend('电网供电功率','燃气供给功率') xlabel('时间');ylabel('功率 kw'); figure;hold on; plot(Pechopr,'-x'); plot(Pechusr,'-o'); legend('电网给充电站供能','用户给充电站供能') xlabel('时间');ylabel('功率 kw'); figure;hold on; plot(Pepvusr,'-x'); plot(Peacusr,'-o'); plot(Phopr,'->'); legend('光伏最大出力','电负荷','热负荷（运营商向园区供热功率）'); xlabel('时间');ylabel('功率 kw'); figure;hold on; plot(Pepvusr,'-x'); plot(Peacusr,'-o'); legend('光伏输出','空调供电') xlabel('时间');ylabel('需求 kw'); figure;hold on; bar(Peopr,'FaceColor',[1 1 1]); bar(Pgrid,'FaceColor',[0 0 0]); legend('MT输出（运营商向园区内部供能功率）','电网输出') xlabel('时间');ylabel('功率 kw'); figure;hold on; bar(Phchpopr+Phgbopr,'FaceColor',[1 1 1]); bar(Phgbopr,'FaceColor',[0 0 0]); legend('MT热输出','燃气热输出') xlabel('时间');ylabel('功率 kw'); figure;hold on; bar(Pepvusr+Peuopr,'FaceColor',[1 1 1]); bar(Pepvusr,'FaceColor',[0 0 0]); legend('外部输入','光伏输出') xlabel('时间');ylabel('需求 kw'); figure;hold on; bar(Phopr+Phusr,'FaceColor',[1 1 1]); bar(Phusr,'FaceColor',[0 0 0]); legend('外部输入热能','电空调热能') xlabel('时间');ylabel('需求 kw'); figure;hold on; plot(Pechusr,'-x'); plot(Peacusr,'-o'); plot(Phopr,'->'); legend('用户供应车辆充电负荷','自身电负荷','自身热负荷（运营商向园区供热功率）'); xlabel('时间');ylabel('功率 kw'); figure;hold on; plot(Pechopr,'-x'); plot(Peuopr,'-o'); plot(Phgbopr,'->'); legend('电网供应充电负荷','电网供应用户电负荷','热网供应用户热负荷'); xlabel('时间');ylabel('功率 kw'); Csopr=zeros(1,T);%运营商售能收入 Cbopr=zeros(1,T);%购电和购热费用之和 Ceoopr=zeros(1,T);%售电维护成本 Choopr=zeros(1,T);%售热维护成本 for ii=1:T Csopr(ii)=peopr(ii)*Peopr(ii)+phopr(ii)*Phopr(ii);%运营商售能收入 Cbopr(ii)=pgrid(ii)*Pgrid(ii)+pgas*Qgas(ii);%购电和购热费用之和 Ceoopr(ii)=ceopr*Peopr(ii)+ceoprchp*Pechpopr(ii);%售电维护成本 Choopr(ii)=cgbhopr*Phgbopr(ii)+choprchp*Phchpopr(ii);%售热维护成本 end Coopr=sum(Ceoopr)+sum(Choopr);%总维护成本 Ropr=sum(Csopr)-sum(Cbopr)-Coopr;%能源运营商的总收益 Csusr=zeros(1,T);%用户售能收入 Cbusr=zeros(1,T);%用户从运营商购能收入 Cousr=zeros(1,T);%含分布式光伏的用户维护成本 for ii=1:T Csusr(ii)=peusr(ii)*Pechusr(ii)+palow*Pepvusr(ii);%用户售能收入 Cbusr(ii)=peopr(ii)*Peuopr(ii)+phopr(ii)*Phopr(ii);%用户从运营商购能收入 Cousr(ii)=cepvusr*Pepvusr(ii)+ceacusr*Peacusr(ii);%含分布式光伏的用户维护成本 end Rusr=sum(Csusr)-sum(Cbusr)-sum(Cousr);%用户的总收益 Cbch=zeros(1,T);%充电站用能成本 for ii=1:T Cbch(ii)=peopr(ii)*Pechopr(ii)+peusr(ii)*Pechusr(ii);%充电站用能成本 end Rch=Cchs-sum(Cbch);%EV代理商效用函数 aa=1; fprintf('运营商售能收入:%0.2f\n',sum(Csopr)); fprintf('用户售能收入:%0.2f\n',sum(Csusr)); fprintf('运营商总成本:%0.2f\n',sum(Cbopr)+Coopr); fprintf('用户用能成本:%0.2f\n',sum(Cbusr)+sum(Cousr)); fprintf('充电站用能成本:%0.2f\n',sum(Cbch));