33e24c12b6e6b6e79915f1b4365f90c72f2fb307
[valse.git] / reports / bazar Emilie / script.m
1 %clear all
2 load donneesSelec.mat
3 %donnees=[ind100,res100];
4 donnees = BB;
5 Res=mean(donnees,2);
6 for i=1:349
7 X2(:,i)=Res(48*(i-1)+1:48*i);
8 end
9 for i=1:348
10 signal(:,i)=[X2(:,i);X2(:,i+1)];
11 end
12
13 [p1,n1]=size(X2);
14 for i=1:n1
15 for j=1:p1
16 XC(i,j)=X2(j,i)-mean(X2(:,i));
17 end
18
19 [C,L]=wavedec(XC(i,:)',4,'haar');
20 xProj(i,:)=C(1:12);
21 end
22
23 X=xProj(1:end-2,:);
24 Y=xProj(2:end-1,:);
25 for i=1:n1
26 Xrecon(i,:)=waverec([xProj(i,:),zeros(1,36)],L,'haar');
27 end
28 for i=1:348
29 signal3(:,i)=[Xrecon(i,:),Xrecon(i+1,:)];
30 end