function [R,W]=soluciona_escalado(U,W,n)
% U contiene los espectrogramas de las señales separadas, sin problemas de
% permutacion, W las matrices de separacion para cada frecuencia y n el
% numero de señales. Devuelve los espectrogramas y las matrices W con el
% escalado correcto

for f=1:length(U(:,1,1))
    D=diag(diag(pinv(W(:,:,f))));
    W(:,:,f)=D*W(:,:,f);
    for k=1:n
        R(f,:,k)=D(k,k)*U(f,:,k);
    end
end