clear all

% load('Feeder_houses.mat')
load('FeederSolar.mat')
load('AllData.mat')
load('Feeder_date_time.mat')
load('q.mat')


SolarHouses=11;
DomesticBatterySize=140*10^3;
n=1;

AllData=cell2mat(AllData);
AllData=AllData(:,[40:94]);
AllDataMat=AllData;
%q=randperm(55,11);

while n<=SolarHouses
    p=q(n);
    House=AllDataMat(:,p);
    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;
         if  (0.2*DomesticBatterySize <= sum(B))&&(sum(B) < DomesticBatterySize*0.8 )
            Grid(i)=0;
         elseif sum(B) >= DomesticBatterySize*0.8
            B(i)=0;
            Grid(i)= (FeederSolar(i)*1.14)-House(i);
         elseif sum(B)< DomesticBatterySize*0.2
            B(i)=0;
            Grid(i)= (FeederSolar(i)*1.14)-House(i);
     end
     AllDataMat(:,p)=-Grid;
end
    n=n+1;
end

% vector=sum(AllDataMat,2);
% for k = 1:1:length(vector)
%    if k<49
%         Min_Mean_Ratio(k)=0;
%    else
%         Min_Mean_Ratio(k)=min(vector(k-48:k))/mean(vector(k-48:k));
%    end 
% end  
   
figure
plot(Feeder_date_time,sum(AllData,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')