view functions/betaDivergencePerElement.m @ 0:b4e26b53072f tip

Initial commit.
author Holger Kirchhoff <holger.kirchhoff@eecs.qmul.ac.uk>
date Tue, 04 Dec 2012 13:57:15 +0000
parents
children
line wrap: on
line source
function betaDivPerElem = betaDivergencePerElement(x, y, beta)

switch beta
    case 0
        betaDivPerElem = x./y - log(x./y) - 1;

    case 1
        betaDivPerElem = x .* log(x./y) + x - y;

    otherwise
        betaDivPerElem = (x.^beta) / (beta*(beta-1)) ...
                         + (y.^beta) / beta ...
                         - (x .* (y.^(beta-1))) / (beta-1);
end