BatterySizeRange=(0:10000:150*10^3);
HouseNumberRange=11;
p=q;   
resultsPeak = zeros(1,length(BatterySizeRange));

n=1;
Vector=1;
export=0;

for i = 1:length(BatterySizeRange)
    export=0;
    [vector] = FeederDomestic(BatterySizeRange(i),11,q);
for k = 1:1:length(vector)
   if k<49
        Peak_Mean_Ratio(k)=0;
        Min_mean_ratio(k)=0;
   else
        Peak_Mean_Ratio(k)=max(vector(k-48:k))/mean(vector(k-48:k));
        Min_Mean_Ratio(k)=min(vector(k-48:k))/mean(vector(k-48:k));
        MaxVec(k)=max(vector(k-48:k));
        MinVec(k)=min(vector(k-48:k));
   end
   
   if vector(k) < 0
        export = export + vector(k);
   end
   
end
resultsPeakMean(i)=mean(Peak_Mean_Ratio);
resultsMinMean(i)=mean(Min_Mean_Ratio);
resultsMax(i)=mean(MaxVec);
resultsMin(i)=mean(MinVec);
exports(i)=export;
n=n+1;
end

figure
[a,b,c]=plotyy(BatterySizeRange/1000,resultsMin/1000,BatterySizeRange/1000,resultsMax/1000);
ylabel(a(1),'Moving Minimum Mean (kW)')
ylabel(a(2),'Moving Maximum Mean (kW)')
xlabel('Battery Size (kWh)')

figure
plot(BatterySizeRange,exports)
ylabel('Total Exported Energy (Wh)')
xlabel('Domestic Battery Size (kWh)')