| 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 |