% close all
load('Feeder.mat', 'Feeder')

BatterySize= 100*10^3;

Feeder_Power=Feeder*2;
Feeder_Bat=zeros(length(Feeder_Power), 1);
Mov_Avg_Day= tsmovavg(Feeder_Power, 's' , 48, 1);
Mov_Avg_Hour= tsmovavg(Feeder_Power, 's' , 8, 1);
std_day=zeros(length(Feeder_Power),1);

for i= 1:length(Feeder_Power)
    if i<49
        std_day(i)=0;
    else
       std_day(i)=std(Feeder_Power(i-48:i));
    end 
end

Dispatch=Mov_Avg_Day+std_day/1.5;
Charge=Mov_Avg_Day-std_day/1.5;

for i = 1:length(Feeder_Power)
    
if Feeder_Power(i) > Dispatch(i) && sum(Feeder_Bat) > 0
    Feeder_Bat(i)= Dispatch(i) - Feeder_Power(i);
elseif Feeder_Power(i) < Charge(i) && sum(Feeder_Bat) < BatterySize*2
    Feeder_Bat(i) = Charge(i) - Feeder_Power(i) ;
end

end

Grid=Feeder_Power + Feeder_Bat;
figure
plot(Feeder_date_time,Feeder_Power,Feeder_date_time,Grid,Feeder_date_time,cumsum(Feeder_Bat))
title('Feeder Power without Solar')