修改这段代码，是每一个delta_s的值输出一个对应的t，在把输出的t带入P_t得到对应的动态回收期 delta_s=[0.1:0.02:0.37]; %向共享储能电站售电的价格[0.2 0.21 0.22 0.23 0.24 0.25 0.26 0.27 0.28 0.29 0.3 0.31 0.32 0.33 0.34 0.35 0.36 0.37] t_list = zeros(size(delta_s)); % 对应的t值 for i=1:length(delta_s);%1:length(delta_s) delta_1 = delta_s(i); B1=(sum(delta_1.*P_ess_s(1,:))+sum(delta_1.*P_ess_s(2,:)))*365;%年卖弃光收益 B2(i)=B1; C0=(sum(gamma.*P_load(1,:))+sum(gamma.*P_load(2,:)))*365; %年用户群不使用储能的情况下从电网购电费用（不是很懂） C1=(sum(gamma.*P_grid(1,:))+sum(gamma.*P_grid(2,:)))*365; %年用户群从电网购电费用 C2=(sum(delta.*P_ess_b(1,:))+sum(delta.*P_ess_b(2,:)))*365; %年用户群从储能电站购电 Copr=72*sum(P_max);%年运维成本 Cinv = @(t) (t==1)*(1000*sum(P_max)+1100*sum(E_max)); % 在t=1时为f(t)，否则为0 max_t = -1; max_v = -Inf; for t = 1:20 V_t = sum(((1+g).^(1:t)./(1+i0).^(1:t)).*(B1+C0-Cinv(1:t)-C1-C2-Copr));%净现值 V_t = sum(((1+g).^(1:t)./(1+i0).^(1:t))*(B1+C0-Cinv(1:t)-C1-C2-Copr)); if V_t >0 if max_t == -1 % 第一次发现 V_t 大于0 max_t = t; % 记录最初的年份 max_v = V_t; end end end t_list(i) = max_t; for t=1:max_t P_t =t-1+((sum(1+g).^(1:t)./(1+i0).^(1:t)).*(B1+C0-Cinv(1:t)-C1-C2-Copr))/(((1+g).^max_t./(1+i0).^max_t).*(B1+C0-Cinv(1:t)-C1-C2-Copr));%动态回收期 end end disp(['delta_s = ', num2str(delta_s)]); disp(['净现值大于0的最初年份为: ', num2str(t_list),' 年']); disp(['动态回收期为: ', num2str(P_t),' 年']);

修改后的代码如下：

delta_s=[0.1:0.02:0.37]; %向共享储能电站售电的价格[0.2 0.21 0.22 0.23 0.24 0.25 0.26 0.27 0.28 0.29 0.3 0.31 0.32 0.33 0.34 0.35 0.36 0.37]
t_list = zeros(size(delta_s)); % 对应的t值
P_t_list = zeros(size(delta_s)); % 对应的动态回收期
for i=1:length(delta_s)
delta_1 = delta_s(i);
B1=(sum(delta_1.*P_ess_s(1,:))+sum(delta_1.*P_ess_s(2,:)))*365;%年卖弃光收益
B2(i)=B1;
C0=(sum(gamma.*P_load(1,:))+sum(gamma.*P_load(2,:)))*365; %年用户群不使用储能的情况下从电网购电费用
C1=(sum(gamma.*P_grid(1,:))+sum(gamma.*P_grid(2,:)))*365; %年用户群从电网购电费用
C2=(sum(delta.*P_ess_b(1,:))+sum(delta.*P_ess_b(2,:)))*365; %年用户群从储能电站购电
Copr=72*sum(P_max);%年运维成本
Cinv = @(t) (t==1)*(1000*sum(P_max)+1100*sum(E_max)); % 在t=1时为f(t)，否则为0
max_t = -1;
max_v = -Inf;
for t = 1:20
V_t = sum(((1+g).^(1:t)./(1+i0).^(1:t)).*(B1+C0-Cinv(1:t)-C1-C2-Copr));%净现值
if V_t >0
if max_t == -1 % 第一次发现 V_t 大于0
max_t = t; % 记录最初的年份
max_v = V_t;
end
end
end
t_list(i) = max_t;
P_t = 0;
for t=1:max_t
P_t =t-1+((sum(1+g).^(1:t)./(1+i0).^(1:t)).*(B1+C0-Cinv(1:t)-C1-C2-Copr))/(((1+g).^max_t./(1+i0).^max_t).*(B1+C0-Cinv(1:t)-C1-C2-Copr));%动态回收期
end
P_t_list(i) = P_t;
disp(['delta_s = ', num2str(delta_s(i))]);
disp(['净现值大于0的最初年份为: ', num2str(t_list(i)),' 年']);
disp(['动态回收期为: ', num2str(P_t),' 年']);
end