why mvnpdf returns complex values?
3 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
Nicola Donelli
am 14 Jan. 2015
Beantwortet: Nicola Donelli
am 14 Jan. 2015
Hi, can anybody tell me why sometimes (i.e. with some Variance Matrices) the mvnpdf returns complex values? As far as I know the multivariate normal density is a function from R^d to R...
Thanks
0 Kommentare
Akzeptierte Antwort
Torsten
am 14 Jan. 2015
% FUNCTION Y = MVNPDF(X, [MU], [SIGMA])
%
% X & MU: vectors of same size
% sigma: square matrix
%
% default MU = zeros(size(X))
% SIGMA = diag(ones(size(X)));
function y = mvnpdf(x,mu,sigma)
if nargin < 2
mu = zeros(size(x));
end
if nargin < 3
sigma = diag(ones(size(x)));
end
x = x(:); % Column vector
mu = mu(:); % Column vector
% Check dimensions
n = length(x);
if size(mu,1) ~= n | size(mu,2) ~= 1 | size(sigma,1) ~= n | ...
size(sigma,2) ~= n
error('Parameter dimensions must agree');
end
a = (2*pi)^-(n/2) * det(sigma)^-.5;
b = exp(-.5 * (x-mu).' * inv(sigma) * (x-mu));
y = a * b;
Just check by your own which Expression (a or b) becomes complex in your application.
Best wishes
Torsten.
1 Kommentar
Torsten
am 14 Jan. 2015
Are you sure the covariance matrices you supply are all positive definite ?
Otherwise, the factor "a" may become complex.
Best wishes
Torsten.
Weitere Antworten (1)
Siehe auch
Kategorien
Mehr zu Arduino Hardware finden Sie in Help Center und File Exchange
Produkte
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!