Mercurial > hg > ddm

view lse_stat_par.m @ 4:72c011ed1977 tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
more elaborate example with non-stat. estimate explanation
author smusevic
date Tue, 30 Jul 2013 09:56:27 +0100
parents a4a7e3405062
children
line wrap: on
line source
function lsp = lse_stat_par(sig, krnls, w)
  nr_krnls = size(krnls,1);
  sig_mat = repmat(sig,nr_krnls,1);
  w_mat   = repmat(w  ,nr_krnls,1);
  kk = sum(  krnls.*conj(krnls).*w_mat,2);
  ks = sum(sig_mat.*conj(krnls).*w_mat,2);
  lsp = log(ks/kk*2);
end