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