SelectModel/suppressionmodelesegaux.m

   1 function[B1,B2,glambda,ind,rho,pi]=suppressionmodelesegaux(B1,B2,glambda,rho,pi)
   2
   3     ind=[];
   4     for l=1:length(glambda)
   5         for ll=1:l-1
   6             if B1(:,:,l)==B1(:,:,ll)
   7                 ind=[ind l];
   8             end
   9         end
  10     end
  11     ind=unique(ind);
  12     B1(:,:,ind)=[];
  13     glambda(ind)=[];
  14     B2(:,:,ind)=[];
  15     rho(:,:,:,ind)=[];
  16     pi(:,ind)=[];
  17
  18 end