function [f1bis, f2bis] = separation1(f1,f2,T) ; 

% function [f1bis, f2bis] = separation1(f1,f2,T); 
% 
% Etant donnés deux tableaux f1 et f2 d'entrées positives (tous les deux de
% même taille) et un seuil strictement positif T, renvoie en sortie le
% signal f2bis = f1 là où f1 domine T*f2 (f2bis=0 étant nul ailleurs) et
% f1bis = f2 aux autres points (f1bis étant nul ailleurs). Cete routine est
% utilisée aux fins de la séparation de deux sources via l'étude comparée
% de leurs spectrogrammes.

[M,N]= size(f1); 
A = sign(f1-T*f2); 
f2bis = max(f1.*A, zeros(M,N)); 
f1bis = max(f2.*(-A),zeros(M,N)); 


end

