for i = 1:length(delta_s)
时间: 2024-04-28 15:27:05 浏览: 9
这段代码是一个 for 循环,用来遍历 delta_s 数组中的每个元素。具体来说,i 的取值范围是 1 到 delta_s 数组的长度,即 i 从 1 开始逐渐增加,直到 i 等于 delta_s 数组的长度为止。
在循环中,可以通过 delta_s(i) 访问 delta_s 数组中的第 i 个元素,进而进行相应的操作。整个循环的作用是对于 delta_s 中的每个元素,计算其对应的净现值大于 0 的最初年份。
相关问题
修改这段代码,是每一个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
解释下这段代码delta_s=[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) delta = delta_s(i); B1=(sum(delta_s.*P_ess_s(1,:))+sum(delta_s.*P_ess_s(1,:)))*365;%年卖弃光收益 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; end disp(['delta_s = ', num2str(delta_s)]); disp(['净现值大于0的最初年份为: ', num2str(t_list)],' 年');
这段代码主要是用来计算共享储能电站售电价格 delta_s 对应的净现值大于 0 的最初年份。具体实现步骤如下:
1. 定义共享储能电站售电价格 delta_s 和对应的 t 值。
2. 对每个 delta_s,计算相关的年度收益和成本,包括卖弃光收益、用户群不使用储能的情况下从电网购电费用、用户群从电网购电费用、用户群从储能电站购电费用、年运维成本等。
3. 定义投资成本函数 Cinv,计算在不同年份的投资成本。
4. 对于每个 delta_s,从 t=1 开始逐年计算净现值 V_t,直到 V_t 不再大于 0。记录最初的年份和对应的净现值。
5. 输出每个 delta_s 对应的净现值大于 0 的最初年份。
总的来说,这段代码主要是用来优化共享储能电站售电价格 delta_s,使得在最初的年份净现值大于 0。
相关推荐
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)