function [stdgrid,sumgrid]=battery_sim(BatterySize,End) % clc % close all % clear all load('Mum_NRJ_Data_matlab.mat', 'USEDWATTS', 'GENERATEDWATTS', 'BATTERY') Start=1; %End=3; % BatterySize= 10; UW=USEDWATTS(Start*2880:End*2880); GW=GENERATEDWATTS(Start*2880:End*2880); Grid=zeros(2880*(End-Start), 1 ); B=zeros(2880*(End-Start), 1 ); BS=BatterySize*120000; for i=1:length(UW) B(i)=GW(i)-UW(i); if (0 <= sum(B))&&(sum(B) < BS ) Grid(i)=0; elseif sum(B) >= BS B(i)=0; Grid(i)= GW(i)-UW(i); elseif sum(B)< 0 B(i)=0; Grid(i)= GW(i)-UW(i); end end BatteryCost = BatterySize; stdgrid=std(subplus(-Grid)); sumgrid=sum(subplus(-Grid/120000)); % data = [UW,GW, cumsum(B)/120]; % % figure; % plot(data) % figure; % plot(-Grid)