why mvnpdf returns complex values?

Question

Nicola Donelli 2015-1-14

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/169900-why-mvnpdf-returns-complex-values

回答： Nicola Donelli 2015-1-14

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 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Torsten 2015-1-14

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/169900-why-mvnpdf-returns-complex-values#answer_164839

在 MATLAB Online 中打开

% 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;