clear all

% load('Feeder_houses.mat')
load('FeederSolar.mat')
load('AllData.mat')
load('\\filestore.soton.ac.uk\users\gh10g10\mydocuments\Year1\Summer Project\MATLAB\TVV\Feeder_date_time.mat')

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

B=zeros(length(FeederSolar),1);

AllData=cell2mat(AllData);
AllDataMat=AllData;
AllDataSol=sum(AllDataMat,2)-(0.5*FeederSolar*SolarHouses);
    p= [51];%randi(55,SolarHouses,1);
House1=zeros(length(FeederSolar),5);
    
while n<=SolarHouses
    House=AllDataMat(:,p(n));
    House1(:,n)=(House-1.14*FeederSolar)/1000;
    Grid=zeros(length(FeederSolar),1);
    for i=1:length(FeederSolar)
       B(i)= ((FeederSolar(i)*1.14)-House(i))*0.8;
         if  (0 <= sum(B))&&(sum(B) < DomesticBatterySize*2 )
            Grid(i)=0;
         elseif sum(B) >= DomesticBatterySize*2
            B(i)=0;
            Grid(i)= (-FeederSolar(i)*1.14)+House(i);
         elseif sum(B)< 0
            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), Feeder_date_time,sum(AllData,2), Feeder_date_time,sum(AllDataSol,2))
legend('Modified Data with Solar Power and Domestic Batteries', 'Original Feeder Data')

% Feeder_max=zeros(length(FeederSolar),1);
% Grid_max=zeros(length(FeederSolar),1);
% Grid_mean =zeros(length(FeederSolar),1);
% Peak_Mean_Ratio =zeros(length(FeederSolar),1);
% Feeder_PowerS=sum(AllData,2);
% GridS=sum(AllDataMat,2);
% 
% for i= 1:length(Feeder_PowerS)
%     if i<49
%         Feeder_max(i)=0;
%         Grid_max(i)=0;
%     else
%        Feeder_max(i)=max(Feeder_PowerS(i-48:i));
%        Grid_max(i)=max(GridS(i-48:i));
%        Grid_mean(i)=mean(GridS(i-48:i));
%        Peak_Mean_Ratio(i)= Grid_max(i)/Grid_mean(i);
%     end 
% end
% 
% Peak_Reduction=mean(Grid_max-Feeder_max);
% 
% PMR=mean(Peak_Mean_Ratio);
% 
% Cost_Benefit=(1/(DomesticBatterySize*SolarHouses))*Peak_Reduction;