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