dan@0: function lsp = lse_stat_par(sig, krnls, w)
dan@0:   nr_krnls = size(krnls,1);
dan@0:   sig_mat = repmat(sig,nr_krnls,1);
dan@0:   w_mat   = repmat(w  ,nr_krnls,1);
dan@0:   kk = sum(  krnls.*conj(krnls).*w_mat,2);
dan@0:   ks = sum(sig_mat.*conj(krnls).*w_mat,2);
dan@0:   lsp = log(ks/kk*2);
dan@0: end