%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%% Effet de moyennage des profils temporels 
%%%% H. BEN AHMED, ENS Rennes, SATIE Lab.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


clc;
close all;
clear all;

load conso;    % profil de consommation en W sur 1 semaine pas 1 min
load lambda;     % profil d'irradiance en W/2 sur 1 semaine pas 1 min

T_echan=120;    % Nouveau pas de moyennage des profils PV et de conso

Npoints=length(conso);


Nmoy=floor(Npoints/T_echan);
temps_moy=zeros(1,Nmoy);
Pmoy=zeros(1,Nmoy);
conso_moy=zeros(1,Nmoy);
lambda_moy=zeros(1,Nmoy);

for k=1:Nmoy
    
    M=0;
    S=0;
    for i=1:T_echan
    S(i)=conso(i+(k-1)*T_echan);
    M(i)=lambda(i+(k-1)*T_echan);
    end
    conso_moy(k)=mean(S);
    lambda_moy(k)=mean(M);
    temps_moy(k)=k*T_echan;
end;

for k=1:Nmoy
Pmoy(k)=mean(conso_moy);
end;
conso_moyenne_en_W=Pmoy(1);
conso_max_en_W=max(conso_moy);

save conso_moy;
save lambda_moy;
save temps_moy;
save Pmoy;
save Nmoy;

figure (1)

subplot(211);
plot(conso)
hold on
plot(temps_moy, conso_moy,'r');
plot(temps_moy, Pmoy,'g');
axis([0 max(temps_moy) 0 max(conso)+0.1*max(conso)]);
hold off

title('Profil de consommation en W sur 1 semaine')
xlabel('temps (min)');
ylabel('Profils (W)');
legend('conso','conso moyennée','Valeur moyenne')

subplot(212);
plot(lambda)
hold on
plot(temps_moy, lambda_moy,'r');
axis([0 max(temps_moy) 0 max(lambda)+0.1*max(lambda)]);
hold off

title('Profil d irradiance en W/m2 sur 1 semaine')
xlabel('temps (min)');
ylabel('Profils (W)');
periode=int2str(T_echan);
titre='lambda moyennée Techan=';
legend('lambda','lambda moyennée')

