检查以下代码,指出错误的地方:%% 综合能源系统-日前优化调度模型
% 电制氢(P2G)、碳捕集(CCS)
% 阶梯式碳交易机制
clc
clear
%% 输入原始数据, 光照强度:w/m^2, 风速:m/s, 电热冷负荷:kw
% 光照强度
GZ_num = [31.9200000000000,44.7000000000000,43.3800000000000,63.6600000000000,104.520000000000,537.360000000000,694.620000000000,823.380000000000,913.260000000000,1069.38000000000,1134.66000000000,1150.68000000000,1120.56000000000,1148.82000000000,1131.90000000000,1105.38000000000,1063.14000000000,944.460000000000,511.560000000000,95.4600000000000,33.1200000000000,64.8600000000000,27.8400000000000,17.4600000000000];
% 风速大小
FS_num = [10.7100000000000,10.8000000000000,10.8300000000000,10.3350000000000,10.7700000000000,10.6050000000000,9.64500000000000,8.95500000000000,9.49500000000000,9.84000000000000,9.94500000000000,10.0500000000000,9.97500000000000,10.3800000000000,10.2900000000000,10.2750000000000,10.3350000000000,10.5150000000000,10.2300000000000,9.55500000000000,10.5750000000000,10.5900000000000,10.5450000000000,10.6500000000000];
% 电负荷
P_load_e = [1086.76000000000,1052.42000000000,1044.34000000000,1057.47000000000,1128.17000000000,1225.13000000000,1367.54000000000,1422.08000000000,1443.29000000000,1464.50000000000,1355.42000000000,1338.25000000000,1326.13000000000,1207.96000000000,1207.96000000000,1214.02000000000,1328.15000000000,1346.33000000000,1334.21000000000,1308.96000000000,1266.54000000000,1237.25000000000,1170.59000000000,1128.17000000000];
% 热负荷
P_load_h = [1134.23000000000,1236.24000000000,1257.45000000000,1254.42000000000,1297.85000000000,1207.96000000000,1184.73000000000,1295.83000000000,1032.22000000000,923.140000000000,872.640000000000,821.130000000000,882.740000000000,905.970000000000,842.340000000000,839.310000000000,938.290000000000,1057.47000000000,1065.55000000000,1090.80000000000,1065.55000000000,1061.51000000000,1050.40000000000,1061.51000000000];
% 冷负荷
P_load_c = [503.280000000000,447.840000000000,450.720000000000,471.600000000000,479.520000000000,408.240000000000,441.360000000000,432.720000000000,398.160000000000,426.240000000000,481.680000000000,501.120000000000,524.160000000000,501.840000000000,492.480000000000,522,456.480000000000,442.800000000000,491.760000000000,496.800000000000,466.560000000000,447.840000000000,468.720000000000,460.080000000000];
% 上网购电价格,采用分时电价
C_e_buy = [0.47*ones(1,6),1.24*ones(1,7),0.89*ones(1,6),1.24*ones(1,5)];
% 上网购气价格,天然气的价格是 2.46或 1.57元/立方米,折合成单位热值的价格为 0.348 元/(k·Wh),采用分时气价;
C_g_buy = [2.46*ones(1,6),1.57*ones(1,7),2.46*ones(1,6),1.57*ones(1,5)];
%% 风、光机组设备模型
% 光伏机组模型
S_PV = 10000; % 光伏面积
% 1:计算光伏发电功率
P_pv_max = 0.001*0.157*S_PV*GZ_num; % 光伏发电效率为0.157,(光照强度已载入),0.001是指将光照强度换算单位:w/m^2转换成kw/m2
% 2:计算光伏运维成本
PV_om_fix = 0.02*0.15*S_PV; % 固定运维成本,按天结算,如设备清灰这样的人工费,一平米折合约0.02*0.15元
PV_om_work = 0.039*sum(P_pv_max); % 可变运维成本,与发电量有关,运维系数为:0.039元/kw
PV_om = PV_om_fix+PV_om_work; % 总运维成本
% 风电机组模型
for i=1:24 % 风机,单台额定2MW,切入风速3,额定风速12,切出风速20
if 3<=FS_num(i) && FS_num(i)<=12
Pe_swind(i)=222.2*(FS_num(i)-3); % 2000=2220*(12-3),因为上限容量为2MW,所以风速大于12时满额输出
elseif 12<FS_num(i) && FS_num(i)<=20
Pe_swind(i)=2000; % 单台额定容量2MW
else
Pe_swind(i)=0; % 其余情况风力不发电
end
end
% 1:计算风电发电功率
P_wt_max = Pe_swind;
% 2:计算风电运维成本
WT_om=0.07*sum(P_wt_max); % 风电运维总成本,运维系数为:0.07元/kw
%% 各变量及常量定义
P_g_CCHP = sdpvar(1,24); % CCHP消耗天然气功率
P_CCHP_e = sdpvar(1,24); % CCHP的输出电功率
P_CCHP_h = sdpvar(1,24); % CCHP的输出热功率
P_CCHP_c = sdpvar(1,24); % CCHP的输出冷功率
P_g_GB = sdpvar(1,24); % 输入GB设备的天然气功率
P_GB_h = sdpvar(1,24); % GB设备输出的热功率
P_e_EB = sdpvar(1,24); % 输入EB设备的电功率
P_EB_h = sdpvar(1,24); % EB设备输出的热功率
P_e_EC = sdpvar(1,24); % 输入电制冷机EC的耗电功率
P_EC_c = sdpvar(1,24); % 电制冷机EC输出的制冷功率
P_e_EL = sdpvar(1,24); % 输入电解槽EL的耗电功率
P_EL_H = sdpvar(1,24); % 电解槽EL输出的产氢功率
P_H_MR = sdpvar(1,24); % 输入甲烷反应器MR的氢能功率
P_MR_g = sdpvar(1,24); % 甲烷反应器MR输出的天然气功率
P_e_CCS = sdpvar(1,24); % 输入碳捕集CCS的耗电量
C_MR_cc = sdpvar(1,24); % 碳捕集CCS输出的碳捕集量(等于MR所用的二氧化碳量,单位为t)
P_wt_e = sdpvar(1,24); % 风电功率
P_pv_e = sdpvar(1,24); % 光伏功率
P_buy_e = sdpvar(1,24); % 购电功率
P_buy_g = sdpvar(1,24); % 购气功率
% 多元储能部分(电ES1、热ES2、冷ES3)
P_ES1_cha = sdpvar(1,24);
P_ES2_cha = sdpvar(1,24);
P_ES3_cha = sdpvar(1,24);
% 充放功率
P_ES1_dis = sdpvar(1,24);
P_ES2_dis = sdpvar(1,24);
P_ES3_dis = sdpvar(1,24);
% 各储能的实时容量状态
SOC_1 = sdpvar(1,24);
SOC_2 = sdpvar(1,24);
SOC_3 = sdpvar(1,24);
% 充放标志二进制变量
B_ES1_cha = binvar(1,24);
B_ES2_cha = binvar(1,24);
B_ES3_cha = binvar(1,24);
% 充、放标志
B_ES1_dis = binvar(1,24);
B_ES2_dis = binvar(1,24);
B_ES3_dis = binvar(1,24);
%% 约束条件
Constraints=[];
Constraints=[Constraints,
P_CCHP_e==0.73*P_g_CCHP, % CCHP的气-电能量转换约束
P_CCHP_h==0.35*P_g_CCHP, % CCHP的气-热能量转换约束
P_CCHP_c==0.28*P_g_CCHP, % CCHP的气-冷能量转换约束
0<=P_g_CCHP<=1200, % CCHP消耗的气功率上下限约束
-0.2*1200<=P_g_CCHP(2:24)-P_g_CCHP(1:23)<=0.2*1200, % CCHP的爬坡约束(1-24时段)
];
Constraints=[Constraints,
P_GB_h==0.91*P_g_GB, % GB的气-热能量转换约束
0<=P_g_GB<=800, % GB的出力上下限约束
-0.2*800<=P_g_GB(2:24)-P_g_GB(1:23)<=0.2*800, % GB的爬坡约束(1-24时段)
];
Constraints=[Constraints,
P_EB_h==1.86*P_e_EB, % EB的电-热能量转换约束
0<=P_e_EB<=600, % EB的出力上下限约束
-0.2*600<=P_e_EB(2:24)-P_e_EB(1:23)<=0.2*600, % EB的爬坡约束(1-24时段)
];
Constraints=[Constraints,
P_EC_c==0.73*P_e_EC, % EC的电-冷的约束关系
0<=P_e_EC<=500, % EC上下限约束
-0.2*500<=P_e_EC(2:24)-P_e_EC(1:23)<=0.2*500, % EC的爬坡约束(1-24时段)
];
Constraints=[Constraints,
P_EL_H==0.87*P_e_EL, % EL(电解槽)的电-氢能量转换约束
0<=P_e_EL<=500, % EL的消耗电功率的上下限约束
-0.2*500<=P_e_EL(2:24)-P_e_EL(1:23)<=0.2*500, % EL的爬坡约束(1-24时段)
];
Constraints=[Constraints,
P_MR_g==0.63*P_H_MR, % MR(甲烷反应器)的氢-气能量转换约束
0<=P_H_MR<=350, % MR消耗的氢功率的上下限约束
-0.2*350<=P_H_MR(2:24)-P_H_MR(1:23)<=0.2*350, % MR的爬坡约束(1-24时段)
];
% 带P2g-CCS协同运行部分
Constraints=[Constraints,
C_MR_cc==1.06*P_MR_g, % CCS捕捉到的CO2与其输出的天然气功率的关系
P_e_CCS==0.55*C_MR_cc, % CCS碳捕捉所耗的电能与捕捉到CO2的质量之间的关系
0<=C_MR_cc<=500, % CCS碳捕集捕获能力上下限约束
-0.2*500<=C_MR_cc(2:24)-C_MR_cc(1:23)<=0.2*500, % CCS的出力爬坡约束(1-24时段)
0<=P_e_CCS<=500, % CCS设备的耗电功率上下限约束
];
Constraints=[Constraints,
0<=P_ES1_cha<=B_ES1_cha*0.5*450, % 储电设备的最大充电功率约束
0<=P_ES2_cha<=B_ES2_cha*0.5*500, % 储热设备的最大充热功率约束
0<=P_ES3_cha<=B_ES3_cha*0.5*350, % 储冷设备的最大充气功率约束
0<=P_ES1_dis<=B_ES1_dis*0.5*450, % 储电设备的最大放电功率约束
0<=P_ES2_dis<=B_ES2_dis*0.5*500, % 储热设备的最大放热功率约束
0<=P_ES3_dis<=B_ES3_dis*0.5*350, % 储冷设备的最大放气功率约束
SOC_1(1)==0.3*450, % 储电设备的初始容量
SOC_2(1)==0.3*500, % 储热设备的初始容量
SOC_3(1)==0.3*350, % 储冷设备的初始容量
% 始末状态守恒约束
SOC_1(24)==SOC_1(1),
SOC_2(24)==SOC_2(1),
SOC_3(24)==SOC_3(1),
% 充放状态唯一
B_ES1_cha+B_ES1_dis<=1,
B_ES2_cha+B_ES2_dis<=1,
B_ES3_cha+B_ES3_dis<=1,
% 储能容量上下限约束
0.2*450<=SOC_1<=0.8*450,
0.2*500<=SOC_2<=0.8*500,
0.2*350<=SOC_3<=0.8*350,
% 储能容量变化约束
SOC_1(2:24)==SOC_1(1:23)+0.95*P_ES1_cha(2:24)-P_ES1_dis(2:24)/0.95,
SOC_2(2:24)==SOC_2(1:23)+0.95*P_ES2_cha(2:24)-P_ES2_dis(2:24)/0.95,
SOC_3(2:24)==SOC_3(1:23)+0.95*P_ES3_cha(2:24)-P_ES3_dis(2:24)/0.95,
];
Constraints=[Constraints,
0.53*P_wt_max<=P_wt_e<=P_wt_max, % 风电功率
0.67*P_pv_max<=P_pv_e<=P_pv_max, % 光伏功率
0<=P_buy_e<=2000, % 购电功率
0<=P_buy_g<=2000, % 购气功率
];
Constraints=[Constraints,
P_buy_e+P_wt_e+P_pv_e+P_CCHP_e+P_ES1_dis==P_load_e+P_e_EB+P_e_EL+P_e_CCS+P_e_EC+P_ES1_cha, % 电功率平衡约束
P_CCHP_h+P_GB_h+P_EB_h+P_ES2_dis==P_load_h+P_ES2_cha, % 热功率平衡约束
P_EC_c+P_CCHP_c+P_ES3_dis==P_load_c+P_ES3_cha,% 冷功率平衡约束
P_buy_g+P_MR_g==P_g_GB+P_g_CCHP, % 气功率平衡约束
P_EL_H==P_H_MR, % 氢功率平衡约束
];
%% 阶梯式碳交易机制
% 碳排放权配额模型
E_e_buy = 0.728*sum(P_buy_e); % 购电配额
E_CCHP = 0.102*3.6*sum(P_CCHP_h+6/3.6*P_CCHP_e+6/3.6*P_CCHP_c); % CCHP配额(产生CO2的机组)
E_GB = 0.102*3.6*sum(P_GB_h); % GB配额(产生CO2的机组)
E_IES = E_e_buy+E_CCHP+E_GB; % IES总碳排放配额
E1 = 0.728*(P_buy_e)+0.102*3.6*(P_CCHP_h+6/3.6*P_CCHP_e+6/3.6*P_CCHP_c)+0.102*3.6*(P_GB_h); % 细化每小时碳排放配额
% 实际碳排放模型
E_e_buy_a = 1.08*sum(P_buy_e); % 实际购电产生的等价CO2
E_CCHP_a = 0.789*sum(P_CCHP_h+6/3.6*P_CCHP_e+6/3.6*P_CCHP_c); % 实际CCHP运行过程产生的CO2,0.789kg/kWh
E_GB_a = 0.25*sum(P_GB_h); % 实际GB运行过程产生的CO2,0.25kg/kWh
E_MR_a = 0.91*sum(P_MR_g); % 实际甲烷生成器MR减少的碳排放(合成天然气需要的CO2)
E = E_e_buy_a+E_CCHP_a+E_GB_a-E_MR_a-E_IES; % 实际IES总碳排放
E2 = 1.58*(P_buy_e)+0.789*(P_CCHP_h+6/3.6*P_CCHP_e+6/3.6*P_CCHP_c)+0.25*(P_GB_h)-0.9*(P_MR_g)-E1; % 细化每小时实际碳排放
% 阶梯碳交易成本(分段线性化)
E_v = sdpvar(1,5); % 每段区间内的长度,分为5段,每段长度是2000
lamda = 0.268; % 碳交易基价:0.268元/kg
Constraints=[Constraints,
E==sum(E_v), % 总长度等于E
0<=E_v(1:4)<=2000, % 除了最后一段,每段区间长度小于等于2000
0<=E_v(5)];
Constraints=[Constraints,C_MR_cc<=E2]; % 每个时刻的二氧化碳捕获量必须小于或等于实际碳排放量(超过捕获限制的二氧化碳将释放到大自然)
% 碳交易成本
C_CO2=0;
for v=1:5
C_CO2=C_CO2+(lamda+(v-1)*0.25*lamda)*E_v(v);
end
% 机组运维成本
C_om=0;
for t=1:24
C_om=C_om+P_g_CCHP(t)*0.04+P_g_GB(t)*0.025+P_e_EB(t)*0.021+P_e_EC(t)*0.028+P_e_EL(t)*0.023+P_H_MR(t)*0.031+C_MR_cc(t)*0.012+...
0.022*(P_ES1_cha(1,t)+P_ES1_dis(1,t))+0.022*(P_ES2_cha(1,t)+P_ES2_dis(1,t))+0.022*(P_ES3_cha(1,t)+P_ES3_dis(1,t));
end
C_om=C_om+PV_om+WT_om; % 加上风、光机组运维成本
% 购能成本(电、气)
Cost_e=0;
Cost_g=0;
for t=1:24
Cost_e=Cost_e+C_e_buy(t)*P_buy_e(t)'; % 购电成本
end
for t=1:24
Cost_g=Cost_g+C_g_buy(t)*P_buy_g(t)'; % 购气成本
end
Cost_buy=Cost_e+Cost_g; % (购电+购气)
% 弃风弃光成本
Cost_WT=0;
Cost_PV=0;
for t=1:24
Cost_WT=Cost_WT+0.927.*(P_wt_max(t)-P_wt_e(t))'; % 弃风成本
end
for t=1:24
Cost_PV=Cost_PV+0.927.*(P_pv_max(t)-P_pv_e(t))'; % 弃光成本
end
Cost_FG=Cost_WT+Cost_PV;
%% 总目标函数
obj = Cost_buy+C_om+C_CO2+Cost_FG; % 购能成本+运维成本+环境成本+弃风弃光
%% 模型求解
ops=sdpsettings('solver','cplex','verbose',0,'usex0',0);
ops.cplex.mip.tolerances.mipgap=1e-6;
result=optimize(Constraints,obj,ops);
Cost_F=double(obj);
display(['通过Yalmip求得的日前最优规划值为(总成本/元 ) : ', num2str(Cost_F)]);
Cost_E_buy=value(Cost_e);
display(['购电成本/元 : ', num2str(Cost_E_buy)]);
Cost_Gas_buy=value(Cost_g);
display(['购气成本/元 : ', num2str(Cost_Gas_buy)]);
Cost_C_om=value(C_om);
display(['运维成本/元 : ', num2str(Cost_C_om)]);
Cost_CO2_trade=value(C_CO2);
display(['碳交易成本/元 : ', num2str(Cost_CO2_trade)]);
Cost_Q_WT=value(Cost_WT);
display(['总弃风成本/元 : ', num2str(Cost_Q_WT)]);
Cost_Q_PV=value(Cost_PV);
display(['总弃光成本/元 : ', num2str(Cost_Q_PV)]);
disp(['---------------新能源机组运行结果--------------']);
disp(['风能利用率 : ', num2str(sum(value(P_wt_e))/sum(P_wt_max)*100),' %']);
Cost_Cost_WT1=value(0.01*(sum(P_wt_max)-sum(value(P_wt_e))));
display(['风电消纳量/MWh : ', num2str(Cost_Cost_WT1)]);
disp(['光伏利用率 : ', num2str(sum(value(P_pv_e))/sum(P_pv_max)*100),' %']);
Cost_Cost_PV1=value(0.01*((sum(P_pv_max)-sum(value(P_pv_e)))));
display(['光电消纳量/MWh : ', num2str(Cost_Cost_PV1)]);
disp(['---------------外部联络线交互结果--------------']);
Cost_Cost_e=value(P_buy_e);
display(['电网联络线交互总功率/kWh : ', num2str(sum(Cost_Cost_e))]);
Cost_Cost_g=value(P_buy_g);
display(['天然气网联络线交互总功率/kWh : ', num2str(sum(Cost_Cost_g))]);
disp(['--------------计算机组运行碳排放量-------------']);
Cost_Cost_E2=value(E2);
display(['总系统碳排放量/t : ', num2str(sum(Cost_Cost_E2))]);
%% 日前调度结果存放
P_buy_e_save = [double(P_buy_e')]; % 购电功率(电出力)
P_buy_g_save = [double(P_buy_g')]; % 购气功率(气出力)
P_g_CCHP_save = [double(P_g_CCHP')]; % CCHP运行功率
P_CCHP_e_save = [double(P_CCHP_e')]; % CCHP的电出力
P_CCHP_h_save = [double(P_CCHP_h')]; % CCHP的热出力
P_CCHP_c_save = [double(P_CCHP_c')]; % CCHP的冷出力
P_GB_h_save = [double(P_GB_h')]; % 燃气锅炉的热出力
P_EB_h_save = [double(P_EB_h')]; % 电锅炉的热出力
P_EC_c_save = [double(P_EC_c')]; % 电制冷机的冷出力
P_EL_H_save = [double(P_EL_H')]; % 电解槽的氢出力
P_MR_g_save = [double(P_MR_g')]; % 甲烷反应器的气出力
%% 画图
figure(1); % 电功率平衡
Figure_E_in=[P_buy_e', P_wt_e', P_pv_e', P_CCHP_e',P_ES1_dis'];
bar(Figure_E_in,'stacked');
hold on
Figure_E_out=[-P_ES1_cha',-P_e_EB',-P_e_EL',-P_e_CCS',-P_e_EC'];
a=bar(Figure_E_out,'stacked');
a(3).FaceColor = [.1 .7 .3];
a(4).FaceColor = [.7 .5 .3];
a(5).FaceColor = [.8 .8 .3];
plot(P_load_e','-g*');
xlabel('时间/h');
ylabel('功率/kw');
title('电功率平衡');
legend('上网电量','风电机组','光伏机组','CCHP机组','储能放电','储能充电','电锅炉','电解槽','碳捕集','电制冷机','电负荷');
grid on;
axis([0 25,-2000,3500])
set(gcf,'color','white');
set(gca,'color','none');
set(gca,'XTick',[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25])
set(gca,'XTickLabel',{'0:00','1:00','2:00','3:00','4:00','5:00','6:00','7:00','8:00','9:00','10:00','11:00','12:00','13:00','14:00','15:00','16:00','17:00','18:00'...
'19:00','20:00','21:00','22:00','23:00','24:00'})
figure(2) % 热功率平衡
Figure_H_in=[P_EB_h',P_GB_h',P_CCHP_h',P_ES2_dis'];
bar(Figure_H_in,'stacked');
hold on
Figure_H_out=[-P_ES2_cha'];
b=bar(Figure_H_out,'stacked');
plot(P_load_h','-r*');
xlabel('时间/h');
ylabel('功率/kw');
title('热功率平衡');
legend('电锅炉','燃气锅炉','CCHP机组','储能放热','储能充热','热负荷');
grid on;
axis([0 25,-500,1500])
set(gcf,'color','white');
set(gca,'color','none');
set(gca,'XTick',[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25])
set(gca,'XTickLabel',{'0:00','1:00','2:00','3:00','4:00','5:00','6:00','7:00','8:00','9:00','10:00','11:00','12:00','13:00','14:00','15:00','16:00','17:00','18:00'...
'19:00','20:00','21:00','22:00','23:00','24:00'})
figure(3) % 冷功率平衡
Figure_C_in=[P_EC_c',P_CCHP_c',P_ES3_dis'];
bar(Figure_C_in,'stacked');
c1=bar(Figure_C_in,'stacked');
hold on
Figure_C_out=[-P_ES3_cha'];
c=bar(Figure_C_out,'stacked');
plot(P_load_c','-b*');
xlabel('时间/h');
ylabel('功率/kw');
title('冷功率平衡');
legend('电制冷机','CCHP机组','储能放冷','储能充冷','冷负荷');
grid on;
axis([0 25,-200,700])
set(gcf,'color','white');
set(gca,'color','none');
set(gca,'XTick',[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25])
set(gca,'XTickLabel',{'0:00','1:00','2:00','3:00','4:00','5:00','6:00','7:00','8:00','9:00','10:00','11:00','12:00','13:00','14:00','15:00','16:00','17:00','18:00'...
'19:00','20:00','21:00','22:00','23:00','24:00'})
figure(4)
x=1:24;
plot(x,P_load_e,'-g*');
hold on
plot(x,P_load_h,'-r*');
hold on
plot(x,P_load_c,'-b*');
ylabel('负荷/kW')
xlabel('时间/h')
title('电、热、冷负荷曲线');
legend('电功率','热功率','冷功率');
axis([0 25,0,1800])
grid on
set(gcf,'color','white');
set(gca,'color','none');
set(gca,'XTick',[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25])
set(gca,'XTickLabel',{'0:00','1:00','2:00','3:00','4:00','5:00','6:00','7:00','8:00','9:00','10:00','11:00','12:00','13:00','14:00','15:00','16:00','17:00','18:00'...
'19:00','20:00','21:00','22:00','23:00','24:00'})
figure(5)
x1=1:24;
plot(x1,P_wt_max,'--bo');
hold on
plot(x1,P_pv_max,'--ro');
hold on
plot(x1,value(P_wt_e),'-b*');
hold on
plot(x1,value(P_pv_e),'-r*');
hold on
xlabel('时间/h')
ylabel('功率/kW')
title('可再生能源消纳水平(日前)');
legend('风电预测出力','光伏预测出力','光伏出力功率','风电出力功率');
axis([0 25,0,2500])
grid on
set(gcf,'color','white');
set(gca,'color','none');
set(gca,'XTick',[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25])
set(gca,'XTickLabel',{'0:00','1:00','2:00','3:00','4:00','5:00','6:00','7:00','8:00','9:00','10:00','11:00','12:00','13:00','14:00','15:00','16:00','17:00','18:00'...
'19:00','20:00','21:00','22:00','23:00','24:00'})
figure(6)
Q_N=double(E2);
Q_CC=double(C_MR_cc);
Q_CN=double(E2-C_MR_cc);
bar([Q_N',Q_CC'])
hold on
plot(Q_CN,'-g*');
xlabel('时段');
ylabel('二氧化碳量/t');
legend('碳排放量','碳捕集量','碳封存量');
title('碳轨迹结果图(日前)');
axis([0 25,-100,800])
grid on
set(gcf,'color','white');
set(gca,'color','none');
set(gca,'XTick',[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25])
set(gca,'XTickLabel',{'0:00','1:00','2:00','3:00','4:00','5:00','6:00','7:00','8:00','9:00','10:00','11:00','12:00','13:00','14:00','15:00','16:00','17:00','18:00'...
'19:00','20:00','21:00','22:00','23:00','24:00'})
最新发布