X-Git-Url: https://git.auder.net/?a=blobdiff_plain;ds=sidebyside;f=reports%2Fbazar%20Emilie%2Frecomp.m;fp=reports%2Fbazar%20Emilie%2Frecomp.m;h=b60ed53cb71d1e96fc6ed4fc85dc61dd5c0cedfb;hb=3f62d540d32b70f42b9090087f72426c18cb219e;hp=0000000000000000000000000000000000000000;hpb=51485a7d0aafe7c31c9651fcc2e33ebd2f8a5e82;p=valse.git

diff --git a/reports/bazar Emilie/recomp.m b/reports/bazar Emilie/recomp.m
new file mode 100644
index 0000000..b60ed53
--- /dev/null
+++ b/reports/bazar Emilie/recomp.m	
@@ -0,0 +1,29 @@
+ychap=zeros(n,m,k,94);
+%Ychap=zeros(n,m,l,94);
+
+for LL=1:10
+    for i=1:n
+        for r=1:2
+            ychap(i,:,r,LL)=x(i,:)*(phitrue(:,:,r,LL)*inv(rhotrue(:,:,r,LL)))';
+            Ychap(i,:,r,LL)=waverec(ychap(i,:,r,LL)',L,'sym4');
+        end
+    end
+end
+for i=1:n
+    for r=1:2
+        for LL=1:10
+            RMSE(i,LL,r)=sqrt(sum((donneesC(78:96,i)-Ychap(i,:,r,LL)').^2));
+        end
+    end
+end
+LL=88;
+plot(Ychap(89,:,2,LL),'r')
+hold on
+plot(donnees(49:96,89)/100,'b')
+hold off
+
+
+for LL=1:10
+    Z(LL,:)=principeMAP(Y,X,phiLassoMLE(:,:,:,LL),rhoLassoMLE(:,:,:,LL),piLassoMLE(:,LL),3.14);
+    sum((Z(1:100)==1)+(Z(101:200)==2))
+end
\ No newline at end of file