GMM（matlab源码）

GMM只是一个数学模型，只是对数据形态的拟和，但是和你所看到的数据分布存在出入也是正常的，因为用EM估计GMM的那些参数时，一般假设我们所得到的数据是不完备的（也就是说假设我们看到的数据分布不是真正的分布，它在运算时把那部分丢失或者叫隐藏的数据“补”上了）

function y = cmvnpdf（X Mu Sigma）
%CMVNPDF     - （Complex range） multivariate normal probability density function （pdf）.
%   Y = CMVNPDF（X） returns the n-by-1 vector Y containing the probability
%   density of the multivariate normal distribution with zero mean and
%   identity covariance matrix evaluated at each row of the n-by-d matrix
%   X. Rows of X correspond to observations and columns correspond to
%   variables or coordinates.
%
%   Y = CMVNPDF（XMU） returns the density of the multivariate normal
%   distribution with mean MU and identity covariance matrix evaluated
%   at each row of X.  MU is a 1-by-d vector or an n-by-d matrix in which
%   case the density is evaluated for each row of X with the corresponding
%   row of MU.  MU can also be a scalar value which CMVNPDF replicates to
%   match the size of X.
%
%   Y = CMVNPDF（XMUSIGMA） returns the density of the multivariate normal
%   distribution with mean MU and covariance SIGMA evaluated at each row
%   of X.  SIGMA is a d-by-d matrix or an d-by-d-by-n array in which case
%   the density is evaluated for each row of X with the corresponding page
%   of SIGMA i.e. CMVNPDF computes Y（I） using X（I:） and SIGMA（::I）.
%   Pass in the empty matrix for MU to use its default value when you want
%   to only specify SIGMA.
%
%   If X is a 1-by-d vector CMVNPDF replicates it to match the leading
%   dimension of MU or the trailing dimension of SIGMA.
%
%   Example:
%
%      mu = [1 -1];
%      Sigma = [.9 .4; .4 .3];
%      X = mvnrnd（mu Sigma 10）;
%      p = cmvnpdf（X mu Sigma）;
%
%   See also MVNRND NORMPDF.

%   Copyright 1993-2002 The MathWorks Inc.
%   Revision: 1.2   Date: 2002/03/28 16:51:27

%   Modified by Pekka Paalanen LUT 2003
%   
%
%   $Name:  $
%   $Revision: 1.1 $  $Date: 2004/08/16 15:06:44 $

if nargin < 1 | isempty（X）
    error（‘Requires the input argument X.‘）;
elseif ndims（X） > 2
    error（‘X must be a matrix.‘）;
end

% Get size of data.  Column vectors provisionally interpreted as multiple scalar data.
[nd] = size（X）;

% Assume zero mean data are already centered
if nargin < 2 | isempty（Mu）
    X0 = X;
 
 % Get scalar mean and use it to center data
elseif prod（size（Mu）） == 1
    X0 = X - Mu;

% Get vector mean and use it to center data
elseif ndims（Mu） == 2
    [n2d2] = size（Mu）;
    if d2 ~= d % has to have same number of coords as X
        error（‘X and MU must have the same number of columns.‘）;
    elseif n2 == n % lengths match
        X0 = X - Mu;
    elseif n2 == 1 % mean is a single row rep it out to match data
        X0 = X - repmat（Mun1）;
    elseif n == 1 % data is a single row rep it out to match mean
        n = n2;
        X0 = repmat（Xn21） - Mu;
    else % sizes don‘t match
        error（‘X or MU must be a row vector or X and MU must have the same number of rows.‘）;
    end
     
 else
   error（‘MU must be a matrix.‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     目录          0  2008-10-17 16:53  （GMM）matlab源码

     文件       6613  2005-04-06 21:25  （GMM）matlab源码\cmvnpdf.m

     文件       3208  2005-04-06 21:18  （GMM）matlab源码\covfixer2.m

     文件       1104  2005-06-02 16:10  （GMM）matlab源码\estS_generate.m

     文件       1099  2005-04-04 11:37  （GMM）matlab源码\getargs.m

     文件       1008  2005-04-07 13:58  （GMM）matlab源码\gmmb_cmeans.m

     文件      10323  2005-04-07 18:25  （GMM）matlab源码\gmmb_em.m

     文件       1132  2005-04-07 12:25  （GMM）matlab源码\gmmb_em_init_cmeans1.m

     文件       1221  2005-04-06 21:13  （GMM）matlab源码\gmmb_em_init_fcm1.m

     文件       1227  2005-04-04 11:43  （GMM）matlab源码\init_fcm1.m

     文件        474  2005-04-06 21:09  （GMM）matlab源码\warning_wrap.m

     文件        248  2006-10-17 11:23  （GMM）matlab源码\注意.txt

----------- ---------  ---------- -----  ----

                27657                    12