b0cf534fed95b8fadb1f170091a37923d9eb0a58
[valse.git] / reports / bazar Emilie / pretraitement.m
1 %Dessins prétraitements
2 script
3
4 XC=X2(:,1)-mean(X2(:,1));
5
6 [C1,L1] = wavedec(X2(:,1),4,'haar');
7 xProj1=C(4:12);
8
9 [C2,L2]=wavedec(XC',4,'haar');
10 xProj2=C(1:12);
11
12 Xrecon1 = waverec([zeros(1,3),xProj1',zeros(1,36)],L1,'haar');
13 Xrecon2 = waverec([xProj2',zeros(1,36)],L2,'haar');
14
15 plot(X2(:,1),'LineWidth',2)
16 hold on
17 plot(Xrecon1,'r-o','LineWidth',2)
18 plot(Xrecon2,'g-s','LineWidth',2)
19 xlabel('Instant of the day','FontSize', 30)
20 ylabel(['Load consumption and its reconstructions' sprintf('\n') ' after projecting and preprocessing'],'FontSize', 30)
21 set(gca, 'FontSize', 20, 'fontName','Times');
22 legend('Original signal','Reconstruction after projection and preprocessing 1','Reconstruction after projection and preprocessing 2','Location','NorthWest')