Après mon merge héroïque
[valse.git] / reports / bazar Emilie / scriptEssai.m
1 %clear all
2 load donneesSelec.mat
3 donnees=BB;
4
5 Res=mean(donnees,2);
6 %Ind=mean(donnees(:,1:100),2);
7 for i=1:340
8 X2(:,i)=Res(48*(i-1)+1:48*i);
9 end
10 for i=1:339
11 signal(:,i)=[X2(:,i);X2(:,i+1)];
12 end
13
14 [p1,n1]=size(X2);
15 for i=1:n1
16 [C,L]=wavedec(X2(:,i)',4,'haar');
17 xProj(i,:)=C(4:12);
18 end
19
20 X=xProj(1:end-2,:);
21 Y=xProj(2:end-1,:);
22 for i=1:n1
23 Xrecon(i,:)=waverec([zeros(1,3),xProj(i,:),zeros(1,36)],L,'haar');
24 end
25
26 for i=1:339
27 signal2(:,i)=[Xrecon(i,:),Xrecon(i+1,:)];
28 end