Mercurial > hg > camir-aes2014
diff toolboxes/FullBNT-1.0.7/netlab3.3/gauss.m @ 0:e9a9cd732c1e tip
first hg version after svn
author | wolffd |
---|---|
date | Tue, 10 Feb 2015 15:05:51 +0000 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/toolboxes/FullBNT-1.0.7/netlab3.3/gauss.m Tue Feb 10 15:05:51 2015 +0000 @@ -0,0 +1,33 @@ +function y = gauss(mu, covar, x) +%GAUSS Evaluate a Gaussian distribution. +% +% Description +% +% Y = GAUSS(MU, COVAR, X) evaluates a multi-variate Gaussian density +% in D-dimensions at a set of points given by the rows of the matrix X. +% The Gaussian density has mean vector MU and covariance matrix COVAR. +% +% See also +% GSAMP, DEMGAUSS +% + +% Copyright (c) Ian T Nabney (1996-2001) + +[n, d] = size(x); + +[j, k] = size(covar); + +% Check that the covariance matrix is the correct dimension +if ((j ~= d) | (k ~=d)) + error('Dimension of the covariance matrix and data should match'); +end + +invcov = inv(covar); +mu = reshape(mu, 1, d); % Ensure that mu is a row vector + +x = x - ones(n, 1)*mu; +fact = sum(((x*invcov).*x), 2); + +y = exp(-0.5*fact); + +y = y./sqrt((2*pi)^d*det(covar));