gmm_test.m

patten regnization source从1-14章能运行

clear all
% 1窜拌 
% X1 : 乞闭捞 [1,1], 傍盒魂捞 [0.3 0 ;0 0.2]牢 盒器甫 捞风绰 单捞磐 100俺 积己
muvec = [1 1];
covmat = [0.3 0; 0 0.2];
X1 = generate_gauss(muvec,covmat,100);
% X2 : 乞闭捞 [5,2], 傍盒魂捞 [0.4 0 ;0 0.4]牢 盒器甫 捞风绰 单捞磐 100俺 积己
muvec = [5 2];
covmat = [0.4 0; 0 0.4];
X2 = generate_gauss(muvec,covmat,100);
% X3 : 乞闭捞 [1,5], 傍盒魂捞 [0.3 0 ;0 0.2]牢 盒器甫 捞风绰 单捞磐 100俺 积己
muvec = [1 5];
covmat = [0.3 0; 0 0.2];
X3 = generate_gauss(muvec,covmat,100);

scatter(X1(1,:),X1(2,:),'b.');
hold on 
scatter(X2(1,:),X2(2,:),'r.');
scatter(X3(1,:),X3(2,:),'k.');

display('3俺狼 努贰胶俊 措茄 单捞磐 积己 肯丰. 浚磐虐甫 穿福搁 促澜 窜拌肺!');
pause 

% 2窜拌
% X1 单捞磐俊 措茄 GMM切嚼阑 困茄 檬扁蔼阑 积己(去钦荐: 3)
[mus, sigmas, weights] = initParams(2,3);
% X1 单捞磐俊 措茄 切嚼 矫累
[mus, sigmas, weights] = trainGMM(X1, mus, sigmas, weights, 0.0001, 100, 'diagcov');
% X1 单捞磐俊 切嚼捞 场唱搁 历厘
X1_mus = mus;
X1_sigmas = sigmas;
X1_weights = weights;
% X1 单捞磐俊 措茄 GMM切嚼 肯丰

% X2 单捞磐俊 措茄 GMM切嚼阑 困茄 檬扁蔼阑 积己(去钦荐: 3)
[mus, sigmas, weights] = initParams(2,3);
% X1 单捞磐俊 措茄 切嚼 矫累
[mus, sigmas, weights] = trainGMM(X2, mus, sigmas, weights, 0.0001, 100, 'diagcov');
% X1 单捞磐俊 切嚼捞 场唱搁 历厘
X2_mus = mus;
X2_sigmas = sigmas;
X2_weights = weights;
% X2 单捞磐俊 措茄 GMM切嚼 肯丰

% X3 单捞磐俊 措茄 GMM切嚼阑 困茄 檬扁蔼阑 积己(去钦荐: 3)
[mus, sigmas, weights] = initParams(2,3);
% X3 单捞磐俊 措茄 切嚼 矫累
[mus, sigmas, weights] = trainGMM(X3, mus, sigmas, weights, 0.0001, 100, 'diagcov');
% X3 单捞磐俊 切嚼捞 场唱搁 历厘
X3_mus = mus;
X3_sigmas = sigmas;
X3_weights = weights;
% X3 单捞磐俊 措茄 GMM切嚼 肯丰
display('3俺狼 努贰胶俊 措茄 切嚼 肯丰. 浚磐虐甫 穿福搁 促澜 窜拌肺!');
pause

% 3 窜拌
% 眠沥等 啊快矫救 盒器甫 敲吩 
% X1 单捞磐俊 措茄 1锅掳 去钦荐俊 措茄 敲吩 
m = X1_mus(:,1)';
s = zeros(2,2);
s(1,1) = X1_sigmas(1,1);
s(2,2) = X1_sigmas(2,1);
plotgaus(m,s);
% X1 单捞磐俊 措茄 2锅掳 去钦荐俊 措茄 敲吩 
m = X1_mus(:,2)';
s = zeros(2,2);
s(1,1) = X1_sigmas(1,2);
s(2,2) = X1_sigmas(2,2);
plotgaus(m,s);
% X1 单捞磐俊 措茄 3锅掳 去钦荐俊 措茄 敲吩 
m = X1_mus(:,3);
m = m';
s = zeros(2,2);
s(1,1) = X1_sigmas(1,3);
s(2,2) = X1_sigmas(2,3);
plotgaus(m,s);

% X2 单捞磐俊 措茄 1锅掳 去钦荐俊 措茄 敲吩 
m = X2_mus(:,1)';
s = zeros(2,2);
s(1,1) = X2_sigmas(1,1);
s(2,2) = X2_sigmas(2,1);
plotgaus(m,s);
% X2 单捞磐俊 措茄 2锅掳 去钦荐俊 措茄 敲吩 
m = X2_mus(:,2)';
s = zeros(2,2);
s(1,1) = X2_sigmas(1,2);
s(2,2) = X2_sigmas(2,2);
plotgaus(m,s);
% X2 单捞磐俊 措茄 3锅掳 去钦荐俊 措茄 敲吩 
m = X2_mus(:,3);
m = m';
s = zeros(2,2);
s(1,1) = X2_sigmas(1,3);
s(2,2) = X2_sigmas(2,3);
plotgaus(m,s);

% X3 单捞磐俊 措茄 1锅掳 去钦荐俊 措茄 敲吩 
m = X3_mus(:,1)';
s = zeros(2,2);
s(1,1) = X3_sigmas(1,1);
s(2,2) = X3_sigmas(2,1);
plotgaus(m,s);
% X3 单捞磐俊 措茄 2锅掳 去钦荐俊 措茄 敲吩 
m = X3_mus(:,2)';
s = zeros(2,2);
s(1,1) = X3_sigmas(1,2);
s(2,2) = X3_sigmas(2,2);
plotgaus(m,s);
% X3 单捞磐俊 措茄 3锅掳 去钦荐俊 措茄 敲吩 
m = X3_mus(:,3);
m = m';
s = zeros(2,2);
s(1,1) = X3_sigmas(1,3);
s(2,2) = X3_sigmas(2,3);
plotgaus(m,s);
hold on 
display('3俺狼 努贰胶俊 措茄 啊快矫救 去钦 葛胆傅 敲吩 肯丰. 浚磐虐甫 穿福搁 促澜 窜拌肺!');
pause 

% 4 窜拌
% 阿 努贰胶喊 切嚼 搬苞甫 府胶飘俊 伎阑 捞侩 历厘
mulist = {X1_mus X2_mus X3_mus};
sigmaslist = {X1_sigmas X2_sigmas X3_sigmas};
weightslist = {X1_weights X2_weights X3_weights};
display('3俺狼 努贰胶俊 措茄 啊快矫救 去钦 葛胆傅 颇扼固磐 历厘 肯丰. 浚磐虐甫 穿福搁 促澜 窜拌肺!');
pause 

% 5 窜拌
% 矫氰 单捞磐 1
test=[1.2; 1.3];
%捞力 牢侥秦焊磊  
classes = classify_GMM(test,mulist, sigmaslist, weightslist);
scatter(test(1,:),test(2,:),'kp');
display('矫氰 单捞磐1 俊 措茄 牢侥 搬苞, 促澜 单捞磐甫 牢侥 且妨搁 浚磐虐甫 穿福矫夸');
classes
pause 
% 矫氰 单捞磐 2
%test=[1.2 4.1 1.2 2.5 4.2; 1.3 2.1 5.4 2.2 1.8];
test=[4.1; 2.1];
%捞力 牢侥秦焊磊  
classes = classify_GMM(test,mulist, sigmaslist, weightslist);
scatter(test(1,:),test(2,:),'kp');
display('矫氰 单捞磐1 俊 措茄 牢侥 搬苞, 促澜 单捞磐甫 牢侥 且妨搁 浚磐虐甫 穿福矫夸');
classes

% 矫氰 单捞磐 3
test=[2.1; 4.1];
%捞力 牢侥秦焊磊  
classes = classify_GMM(test,mulist, sigmaslist, weightslist);
scatter(test(1,:),test(2,:),'kp');
display('矫氰 单捞磐3 俊 措茄 牢侥 搬苞, 场');
classes