gausspdf.m

来自「patten regnization source从1-14章能运行」· M 代码 · 共 24 行

function [z] = gausspdf(x,mu,sigma);

% GAUSSPDF 啊快矫救 犬伏 剐档 窃荐蔼 备窍扁 
%
%    GAUSSPDF(X,MU,SIGMA)绰 乞闭 MUdhk 傍盒魂 SIGMA俊 蝶弗 
%    啊快矫救 橇肺技胶俊 蝶弗 单捞磐 X 趣篮 单捞磐狼 快档甫 府畔茄促. 
%    MU绰 1*D 氦磐(咯扁辑 D绰 橇肺技胶狼 瞒盔)
%    SIGMA绰 D*D 青纺
%    X绰 N*D 青纺肺 阿 青捞 窍唱狼 D-瞒盔 单捞磐 痢捞 等促.  
%
%    See also MEAN, COV, GLOGLIKE

[N,D] = size(x);
if (min(N,D) == 1), x=x(:)'; end;

invSig = inv(sigma);
mu = mu(:)';

x = x-repmat(mu,N,1);

z = sum( ((x*invSig).*x), 2 );

z = exp(-0.5*z) / sqrt((2*pi)^D * det(sigma));