1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36
|
function [Y] =Residus_SenseiA(Date_prix)
B=Date_prix(:,1);
D=Date_prix(:,2);
C=weekday(B);
durant_week_end=C==1|C==7;
[~,~,~,heures,~,~]=datevec(Date_prix(:,1));
hors_8h_20h = 8>heures | heures>=20;
entre_8h_20h = 8<=heures & heures<20;
D_semaine_entre_8h_20h =Date_prix(~durant_week_end & entre_8h_20h,2);
D_Offpeak =Date_prix(durant_week_end | hors_8h_20h,2);
Moyenne_Offpeak=mean(D_Offpeak);
Moyenne_peak=mean(D_semaine_entre_8h_20h);
%%%%%%%%%%%%%%%%%%%%%%%%%%%
Donnees_Prix_hors_8h_20h=bsxfun(@times,hors_8h_20h,D);
Moyenne_Prix_hors_8h_20h=bsxfun(@times,Moyenne_Offpeak,hors_8h_20h);
Prix_hors_8h_20h=bsxfun(@minus,Donnees_Prix_hors_8h_20h,Moyenne_Prix_hors_8h_20h);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Donnees_Prix_entre_8h_20h=bsxfun(@times,entre_8h_20h,D);
Moyenne_Prix_entre_8h_20h=bsxfun(@times,Moyenne_peak,entre_8h_20h);
Prix_entre_8h_20h=bsxfun(@minus,Donnees_Prix_entre_8h_20h,Moyenne_Prix_entre_8h_20h);
Y=bsxfun(@plus,Prix_hors_8h_20h,Prix_entre_8h_20h);
end |
Partager