基于非合作博弈的风-光-氢微网容量优化配置

部分代码：

clear
clc
num_WT=40;
num_PV=2603;
num_EC=22;
num_HST=12;
num_FC=30;

load PV;
load WT;
load Load;
%% 参数区
r=0.06; %折旧率（已换算为百分比）
m=15; %使用年限/a
a=0.95;%储氢效率
V_hst=10;%储氢罐的容量
E=1000;%储氢罐标准大气压下的最大储氢量/Nm3
E_hst_all=num_HST*E;%总的可储氢量/Nm3
E_hst_init=E_hst_all*0.5;%假设初始状态的储氢量为一半
P_FC_max=2;%燃料电池额定功率
P_EC_max=13.62;%电解槽的功率，由文中参数计算所得/kW
W_EC=1.941;%单个电解槽的额定产氢量
eff_p2g=0.75;%电解槽效率
W_p2g=0.19;%每度电的制氢量
P_grid=40;%联络线最大功率/kW
P_pv=0.083;%光伏的单板功率/kw
P_wt=10; %风机额定功率/kw
P_fc=2;%燃料电池单机功率/kw
eff_g2p=0.4;%燃料电池的效率
W_g2p=0.8;%生产每度电的耗氢量
e_wt=4912; %风机单位造价/元
e_pv=1400; %光伏板单位造价/元
e_fc=5000; %燃料电池单位造价/元
e_hst=7000; %储氢罐单位造价/元
e_ec=4000; %电解槽单位造价/元
c_buy