混合高斯模型matlab

MATLAB的相关混合高斯模型，用起来很好，有些瑕疵，适合学习了解混合高斯模型算法初学者

function plotall
%
% Copyright （c） 2008 Sylvain Calinon 
% LASA Lab EPFL CH-1015 Lausanne Switzerland 
% http://www.calinon.ch http://lasa.epfl.ch


nbSamples=3;

[priorsMuSigma] = readGmmFile（‘../outdata/gmm.txt‘）;

Data=[];
for n=1:nbSamples
  Datatmp = load（[‘../outdata/data‘ num2str（n‘%.2d‘） ‘_rescaled.txt‘]）‘;
  Data = [Data  Datatmp];
end
nbVar = size（Datatmp1）;
nbData = size（Datatmp2）;

regr = load（[‘../outdata/gmr_Mu.txt‘]）‘;
rgSigma = load（[‘../outdata/gmr_Sigma.txt‘]）‘;
size（regr）
rgSigma = reshape（rgSigmanbVar-1nbVar-1nbData）;

%% Plot GMM
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Axis parameters
minX = min（Data（1:））;
maxX = max（Data（1:））;
dYmax=1E-10;
for i=1:nbVar-1
  dY（i） = max（Data（i+1:）） - min（Data（i+1:））;
  if dY（i）>dYmax
    dYmax=dY（i）;
  end
end
for i=1:nbVar-1
  minY（i） = min（Data（i+1:））+（dY（i）/2）-1.3*（dYmax/2）; 
  maxY（i） = min（Data（i+1:））+（dY（i）/2）+1.3*（dYmax/2）;
end

figure（‘position‘[20-801200700]‘GMM/GMR‘）;
%Plot GMM
for i=1:nbVar-1
  subplot（2nbVar-1i）; hold on; box on;
  plotGMM（Mu（[1 i+1]:）Sigma（[1 i+1][1 i+1]:）[0.3 0.1 0]1）;
  for n=1:nbSamples
    plot（Data（1（n-1）*nbData+1:n*nbData）Data（i+1（n-1）*nbData+1:n*nbData）‘k-‘）;
  end
  axis（[minX maxX minY（i） maxY（i）]）;
  xlabel（‘t‘）; ylabel（[‘x_‘ num2str（i）]）;
end
%Plot GMR
for i=1:nbVar-1
  subplot（2nbVar-1i+nbVar-1）; hold on; box on;
  plotGMM（regr（[1 i+1]:）rgSigma（ii:）[0 0.3 0.1]3）;
  axis（[minX maxX minY（i） maxY（i）]）;
  xlabel（‘t‘）; ylabel（[‘x_‘ num2str（i）]）;
end

pause;
close all;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2008-04-29 11:25  gmm-gmr\indata\
     文件        7390  2008-04-21 07:22  gmm-gmr\indata\data01.txt
     文件        8530  2008-04-21 07:22  gmm-gmr\indata\data02.txt
     文件        8830  2008-04-21 07:22  gmm-gmr\indata\data03.txt
     文件         592  2008-04-21 07:22  gmm-gmr\Makefile
     目录           0  2008-04-29 11:25  gmm-gmr\matlab_src\
     文件        1568  2008-04-21 07:22  gmm-gmr\matlab_src\plotall.m
     文件        2053  2008-04-21 07:22  gmm-gmr\matlab_src\plotGMM.m
     文件         276  2008-04-21 07:22  gmm-gmr\matlab_src\readGmmFile.m
     目录           0  2008-04-29 11:26  gmm-gmr\MSVCproject\
     文件         886  2008-04-21 07:22  gmm-gmr\MSVCproject\MSVCproject.sln
     文件       14336  2008-04-21 07:22  gmm-gmr\MSVCproject\MSVCproject.suo
     文件        4562  2008-04-21 07:22  gmm-gmr\MSVCproject\MSVCproject.vcproj
     目录           0  2008-04-29 11:25  gmm-gmr\obj\
     目录           0  2008-04-29 11:25  gmm-gmr\outdata\
     文件        1167  2008-04-21 07:22  gmm-gmr\README
     目录           0  2008-04-29 11:26  gmm-gmr\src\
     文件       14162  2008-04-21 07:22  gmm-gmr\src\gmr.cpp
     文件        3036  2008-04-21 07:22  gmm-gmr\src\gmr.h
     文件         141  2008-04-21 07:22  gmm-gmr\src\Macros.cpp
     文件        1231  2008-04-21 07:22  gmm-gmr\src\Macros.h
     文件        3250  2008-04-21 07:22  gmm-gmr\src\main.cpp
     文件         494  2008-04-21 07:22  gmm-gmr\src\MathLib.h
     文件         168  2008-04-21 07:22  gmm-gmr\src\Matrix.cpp
     文件       34823  2008-04-21 07:22  gmm-gmr\src\Matrix.h
     文件         176  2008-04-21 07:22  gmm-gmr\src\Vector.cpp
     文件       11810  2008-04-21 07:22  gmm-gmr\src\Vector.h
     目录           0  2008-04-29 11:25  gmm-gmr\