Daniel@0: function W = initializeDiagMKL(X)
Daniel@0:     
Daniel@0:     [d,n,m] = size(X);
Daniel@0: 
Daniel@0:     W = zeros(d,m);
Daniel@0:     for i = 1:m
Daniel@0:         W(:,i) = 1./std(X(:,:,i),1,2);
Daniel@0:     end
Daniel@0:     W(isinf(W)) = 1;
Daniel@0: end