wi=diag(1|G(:,j)||), j ∈ [1, 2, ..., p] is a diagonal matrix for deeper source compensation. G(:, j) is the jth column of G. The algorithm is initialized with the minimum norm solution D^MNE, that is, W0=diag(D^MNE)=diag(D^0(1),D^0(2),...,D^0(3p)), where D^0(n) represents the nth element of vector D^0. If continued long enough, FOCUSS converges to a set of concentrated solutions equal in number to the number of electrodes.
