load('FeederSolar.mat')
 load('AllData.mat')


n=1;



AllData=cell2mat(AllData);
AllData=AllData(:,[40:94]);
AllDataMat=AllData;
p=q;
DomesticBatterySize=154000;
SolarHouses=11;

while n<=SolarHouses
    House=AllDataMat(:,p(n));
    Grid=zeros(length(FeederSolar),1);
    B=zeros(length(FeederSolar),1);
    for i=1:length(FeederSolar)
         B(i)= ((FeederSolar(i)*1.14)-House(i))*0.8;
         B(1)=0.5*(DomesticBatterySize*2/SolarHouses);
         if  (0.2*(DomesticBatterySize*2/SolarHouses) <= sum(B))&&(sum(B) < 0.8*(DomesticBatterySize*2/SolarHouses) )
            Grid(i)=0;
         elseif sum(B) >= (DomesticBatterySize*2/SolarHouses)*0.8
            B(i)=0;
            Grid(i)= (FeederSolar(i)*1.14)-House(i);
         elseif sum(B)< (DomesticBatterySize*2/SolarHouses)*0.2
            B(i)=0;
            Grid(i)= (FeederSolar(i)*1.14)-House(i);
     end
     AllDataMat(:,p(n))=-Grid;
end
    n=n+1;
end

figure
plot(Feeder_date_time,sum(AllDataMat,2)/1000, Feeder_date_time,(sum(AllData,2)-FeederSolar*12.54)/1000)
ylabel('Power (kW)')
xlabel('Date and Time')
legend('Feeder Demand with Batteries','Feeder Demand without Batteries')