% clc
% close all
% clear all
load('Mum_NRJ_Data_matlab.mat', 'USEDWATTS', 'GENERATEDWATTS', 'BATTERY')
Start=9;
End=10;
BatterySize= 7 ;
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;
t=(1:length(UW))*30;


for i=1:length(UW)
B(i)=GW(i)-UW(i);
B(1)=300000;
     if  (BS*0.1 <= sum(B))&&(sum(B) < BS*0.9 )
          Grid(i)=0;
     elseif sum(B) >= BS*0.9
          B(i)=0;
         Grid(i)= GW(i)-UW(i);
     elseif sum(B)< BS*0.1
          B(i)=0;
         Grid(i)= GW(i)-UW(i);
     end
end

data=[UW/1000,GW/1000,-Grid/1000];
figure
[hAx,hLine1,hLine2]=plotyy(t/3600,data,t/3600,cumsum(B/(120*70)));
title('Simple Domestic Battery Simulation')
xlabel('Time(hours)')
ylabel(hAx(1),'Power (kW)')
ylabel(hAx(2),'SOC (%)')
legend('Used Power','Generated Power','Power exchange with grid','Battery SOC','Location','northwest')

print -depsc SimpleDomesticModel