function [f1bis, f2bis, f3bis, f4bis] = separation2(f1,f2,T1,T2) ; 

% [f1bis, f2bis, f3bis,f4bis] = separation1(f1,f2,T1,T2); 
% 
% Etant donnés deux tableaux f1 et f2 d'entrées positives tous les deux de
% même taille et deux seuils T1<T2 strictement positifs, renvoie le signal
% f4bis = f1 là où f1 domine T2*f1 (f4bis=0 étant nul ailleurs), f1bis = f2
% là où f1 est dominé par T1*f2 (f1bis étant nul ailleurs), f2bis = f2
% (resp. f3bis =f1) en tous les autres points (f2bis et f3bis état nuls
% ailleurs). Cette routine est utilisée aux fins de la séparation de trois
% sources via l'étude comparée des spectrogrammes.

[M,N]= size(f1); 
A = sign(f1-T2*f2); 
B = sign(T1*f2-f1); 
C = A.*B ; 
f4bis = max(f1.*A, zeros(M,N));
f1bis = max(f2.*B, zeros(M,N)); 
f2bis = max(f2.*C, zeros(M,N)); 
f3bis = max(f1.*C, zeros(M,N));  
end

