解释一下这段代码gamma=[0.37*ones(1,32),1.36*ones(1,16),0.82*ones(1,20),1.36*ones(1,16),0.82*ones(1,12)]; %电网的"峰-平-谷"电价 delta=0.33*ones(17,96); %从共享储能电站购电的价格 delta_s=0.25*ones(1,96); %向共享储能电站售电的价格 P_load(1,:)=xlsread('fuhe.xlsx','sheet1','A1:CR1'); P_load(2,:)=xlsread('fuhe.xlsx','sheet1','A2:CR2'); P_pv(1,:)=xlsread('fuhe.xlsx','sheet1','A3:CR3'); P_pv(2,:)=xlsread('fuhe.xlsx','sheet1','A3:CR3'); P_ess_s(1,:)=xlsread('P-ess-s.xlsx','sheet1','A1:CR1'); P_ess_s(2,:)=xlsread('P-ess-s.xlsx','sheet1','A2:CR2'); P_ess_b(1,:)=xlsread('P-ess-b.xlsx','sheet1','A1:CR1'); P_ess_b(2,:)=xlsread('P-ess-b.xlsx','sheet1','A2:CR2'); P_grid(1,:)=xlsread('P-grid.xlsx','sheet1','A1:CR1'); P_grid(2,:)=xlsread('P-grid.xlsx','sheet1','A2:CR2'); g=0.03; i0=0.1; P_max=2161.4878 ; E_max=39700.8022 ; %%约束条件 %%目标函数 B1=(sum(delta_s.*P_ess_s(1,:))+sum(delta_s.*P_ess_s(2,:)))*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; %年用户群从储能电站购电 Cinv= [1000*sum(P_max)+1100*sum(E_max)] + [0] * (19); Copr=72*sum(P_max);%年运维成本 max_t = -1; max_v = -Inf; for t = 1:20 V_t = sum((1+g)^t/(1+i0)^t.*(B1+C0-C1-C2-Cinv-Copr)); %求解净现值 if V_t > 0 if V_t > max_v max_v = V_t; max_t = t; end end end if max_t == -1 % 净现值始终小于0 disp('净现值始终小于0'); else % 净现值大于0 disp(['净现值大于0的年份为: ', num2str(max_t),' 年']); % 输出相应的变量值 fprintf('销售弃光电价: %f\n', B1); fprintf('用户配置储能前从电网购电费用: %f\n', C0); fprintf('用户配置储能后从电网购电费用: %f\n', C1); fprintf('储能总投资成本: %f\n', Cinv(1)); fprintf('年运维成本: %f\n', Copr); end
时间: 2024-02-05 17:03:20 浏览: 39
这段代码是为了计算光伏储能系统的经济效益,并对其进行评估。具体来说,该代码中的变量含义如下:
- gamma:电网的电价,其中0.37代表“平”电价,1.36代表“峰”电价,0.82代表“谷”电价。
- delta:从共享储能电站购电的价格,0.33代表每度电的价格。
- delta_s:向共享储能电站售电的价格,0.25代表每度电的价格。
- P_load:负荷,即用户的用电需求。
- P_pv:光伏发电量。
- P_ess_s:储能电池的售电功率。
- P_ess_b:储能电池的购电功率。
- P_grid:用户从电网购电的功率。
- g:折旧率。
- i0:贴现率。
- P_max:储能电池的最大功率。
- E_max:储能电池的最大容量。
在这些变量的基础上,该代码计算了光伏储能系统在不同年份内的净现值,并输出相关变量的值,以判断该系统是否具有经济价值。其中,约束条件和目标函数的具体含义可以参考代码中的注释。
相关问题
w1=gamma*(r0)^3*omega*f1;
这个公式中,w1表示某个物理量,gamma、r0、omega、f1都是参量,可以看作是常数或者已知变量。你可以把它们带入公式中计算出w1的值。具体来说,gamma是一个常数,r0表示某个长度,omega表示某个角速度,f1表示某个力,它们的单位分别为无量纲、米、弧度/秒、牛顿。因此w1的单位为牛顿·米。
修改这段代码,将向共享储能电站售电的价格改为一个从0.2到0.37间隔为0.1的变量 gamma=[0.37*ones(1,32),1.36*ones(1,16),0.82*ones(1,20),1.36*ones(1,16),0.82*ones(1,12)]; %电网的"峰-平-谷"电价 delta=0.33*ones(1,96); %从共享储能电站购电的价格 delta_s=0.25*ones(1,96); %向共享储能电站售电的价格 P_load(1,:)=xlsread('fuhe.xlsx','sheet1','A1:CR1'); P_load(2,:)=xlsread('fuhe.xlsx','sheet1','A2:CR2'); P_pv(1,:)=xlsread('fuhe.xlsx','sheet1','A3:CR3'); P_pv(2,:)=xlsread('fuhe.xlsx','sheet1','A3:CR3'); P_ess_s(1,:)=xlsread('P-ess-s.xlsx','sheet1','A1:CR1'); P_ess_s(2,:)=xlsread('P-ess-s.xlsx','sheet1','A2:CR2'); P_ess_b(1,:)=xlsread('P-ess-b.xlsx','sheet1','A1:CR1'); P_ess_b(2,:)=xlsread('P-ess-b.xlsx','sheet1','A2:CR2'); P_grid(1,:)=xlsread('P-grid.xlsx','sheet1','A1:CR1'); P_grid(2,:)=xlsread('P-grid.xlsx','sheet1','A2:CR2'); g=0.03; i0=0.1; P_max=2161.4878 ; E_max=39700.8022 ; %%约束条件 %%目标函数 B1=(sum(delta_s.*P_ess_s(1,:))+sum(delta_s.*P_ess_s(2,:)))*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 if max_t == -1 % 净现值始终小于0 disp('净现值始终小于0'); else % 净现值大于0 disp(['净现值大于0的最初年份为: ', num2str(max_t),' 年']); end
修改后的代码如下:
```
gamma=0.2:0.1:0.37;
gamma=[gamma,1.36*ones(1,16),0.82*ones(1,20),1.36*ones(1,16),0.82*ones(1,12)]; %电网的"峰-平-谷"电价
delta=0.33*ones(1,96); %从共享储能电站购电的价格
delta_s=gamma; %向共享储能电站售电的价格
P_load(1,:)=xlsread('fuhe.xlsx','sheet1','A1:CR1');
P_load(2,:)=xlsread('fuhe.xlsx','sheet1','A2:CR2');
P_pv(1,:)=xlsread('fuhe.xlsx','sheet1','A3:CR3');
P_pv(2,:)=xlsread('fuhe.xlsx','sheet1','A3:CR3');
P_ess_s(1,:)=xlsread('P-ess-s.xlsx','sheet1','A1:CR1');
P_ess_s(2,:)=xlsread('P-ess-s.xlsx','sheet1','A2:CR2');
P_ess_b(1,:)=xlsread('P-ess-b.xlsx','sheet1','A1:CR1');
P_ess_b(2,:)=xlsread('P-ess-b.xlsx','sheet1','A2:CR2');
P_grid(1,:)=xlsread('P-grid.xlsx','sheet1','A1:CR1');
P_grid(2,:)=xlsread('P-grid.xlsx','sheet1','A2:CR2');
g=0.03;
i0=0.1;
P_max=2161.4878 ;
E_max=39700.8022 ;
%%约束条件
%%目标函数
B1=(sum(delta_s.*P_ess_s(1,:))+sum(delta_s.*P_ess_s(2,:)))*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
if max_t == -1 % 净现值始终小于0
disp('净现值始终小于0');
else % 净现值大于0
disp(['净现值大于0的最初年份为: ', num2str(max_t),' 年']);
end
```
修改后的代码将 gamma 设置为从 0.2 到 0.37,间隔为 0.1 的变量,然后将向共享储能电站售电的价格 delta_s 设置为 gamma,即向共享储能电站售电的价格与电网的“峰-平-谷”电价相同。
相关推荐
![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)