Em算法（使用C++编写）

Em算法(使用C++编写)，使用Em算法对多维向量进行聚类。程序本身还存在一些问题，比如协方差矩阵的行列式为0的问题(给对角线元素加一个小的数并且取协方差阵为对角阵解决的)，希望大家能讨论下。

#include “em.h“
#include “memory.h“
#include “stdlib.h“
#include “stdio.h“
#include “math.h“
#include “time.h“

namespace ker
{
    #define PI 3.1415926
    #define CONST_E 0.0000001
    
    cEm::cEm（int nVec int nDim int nPat）
    {
        _nVec = nVec;
        _nDim = nDim;
        _nPat = nPat;

        //_pplfZ[i][j]表示第i个向量属于第j类的概率初始化为0
        _pplfZ = new double *[_nVec];
        for （int i = 0; i < _nVec; i++）
        {
            _pplfZ[i] = new double[_nPat];
            memset（_pplfZ[i] 0 sizeof（double） * _nPat）;
        }
        
        //_pplfU[i]表示第i类的均值向量，其初始化的值有可能导致陷入局部极值
        //因此其初始化的值根据待分类数据来给出
        _pplfU = new double *[_nPat];
        for （int i = 0; i < _nPat; i++）
        {
            _pplfU[i] = new double[_nDim];
            //for （int j = 0; j < _nDim; j++） _pplfU[i][j] = i * 2;
        }
        
        //_ppplfDelta[i]表示第i类的协方差矩阵并把协方差矩阵初始化成单位阵
        _ppplfDelta = new double **[_nPat];
        for （int i = 0; i < _nPat; i++）
        {
            _ppplfDelta[i] = new double *[_nDim];
            for （int j = 0; j < _nDim; j++）
            {
                _ppplfDelta[i][j] = new double[_nDim];
                memset（_ppplfDelta[i][j] 0 sizeof（double） * _nDim）;
                _ppplfDelta[i][j][j] = 1;
            }
        }
        
        //_plfPi[i]表示第i类的先验概率初始化为1 / _nPat
        _plfPi = new double [_nPat];
        for （int i = 0; i < _nPat; i++） _plfPi[i] = 1.0 / _nPat;
    }

    cEm::~cEm（）
    {
        for （int i = 0; i < _nVec; i++） delete []_pplfZ[i];
        delete []_pplfZ;

        for （int i = 0; i < _nPat; i++） delete[]_pplfU[i];
        delete []_pplfU;
        
        for （int i = 0; i < _nPat; i++）
        {
            for （int j = 0; j < _nDim; j++） delete[]_ppplfDelta[i][j];
            delete []_ppplfDelta[i];
        }
        delete []_ppplfDelta;
        delete []_plfPi;
    }
    
    double **cEm::Inverse（double **pplfMatSrc）
    {
        //复制原矩阵
        double **pplfMat = new double *[_nDim];
        for （int i = 0; i < _nDim; i++）
        {
            pplfMat[i] = new double [_nDim];
            memcpy（pplfMat[i] pplfMatSrc[i] sizeof（double） * _nDim）;
        }
        
        //创建一个单位阵
        double **pplfI = new double *[_nDim];
        for （int i = 0; i < _nDim; i++）
        {
            pplfI[i] = new double [_nDim];
            memset（pplfI[i] 0 sizeof（double） * _nDim）;
            pplfI[i][i] = 1;
        }
        
        //Gaussian消元法求逆矩阵--正向消元
        for （int i = 0; i < _nDim; i++）
        {
            double lfTmp = pplfMat[i][i];
            for （int j = 0; j < _nDim; j++）
            {
                pplfMat[i][j] /= lfTmp;
                pplfI[i][j] /= lfTmp;
            }
            
            for （int j = i + 1; j < _nDim; j++）
            {
                double lfTmp = -pplfMat[j][i];
                for （int k = 0; k < _nDim; k++

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

     文件        271  2009-02-06 21:38  EM\EM.cfp

     文件      10968  2009-02-06 21:34  EM\em.cpp

     文件      76833  2009-02-06 21:38  EM\EM.exe

     文件       2040  2009-02-06 17:14  EM\em.h

     文件      24897  2009-02-06 21:34  EM\em.o

     文件        739  2009-02-06 21:38  EM\main.cpp

     文件       6163  2009-02-06 21:38  EM\main.o

     文件     335872  2009-02-06 00:46  em.ppt

     目录          0  2009-02-06 21:38  EM

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

               457783                    9