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