example_iris_02.m

模糊C均值聚类的一个参考材料,做模糊聚类是用的一个程序

% ----------------------------------------------------------------------
% Example for testing IRID data clustering (variable m and eta)
% 
% ----------------------------------------------------------------------
% Written By: Mahdi Amiri.
% E-Mail: [e-mail]
% Homepage: http://yashil.20m.com/
% June 2003
% ----------------------------------------------------------------------

% ----------------------------------------------------------------------
% Clean up
clc;
clear all;
close all;

% ----------------------------------------------------------------------
% Input data => Xin

load iris_1.dat; Xin = iris_1;

% ----------------------------------------------------------------------
cMarker = ['+' 'o' '*' '.' 'x' 's' 'd' '^' 'v' '>' '<' 'p' 'h'];
cColor =  ['r' 'g' 'b' 'm' 'c' 'y' 'k' 'r' 'g' 'b' 'y' 'm' 'c'];

p = size (Xin, 2);

% ----------------------------------------------------------------------
% Number of clusters => nC
nC = 3;

% ----------------------------------------------------------------------
% Options
term_thr = 0.00001;
max_iter = 100;
info_display = 0;

%init_V = Xin(1:nC, :);
%init_V = rand (nC, p);
init_V = [0.077057 0.408295 0.362549 0.997681;
    0.928650 0.432343 0.231903 0.345886
    0.341949 0.927887 0.713257 0.988281];

% ----------------------------------------------------------------------
% Main loop
m = 1.5;
eta = 3.5;
for i = 1:15

    % ----------------------------------------------------------------------
    % FCM
    [V,U,E] = Yf_FCMC1 (Xin, nC, [m; max_iter; term_thr; info_display; 1], init_V);

    V_FCM = V;
    U_FCM = U;
    E_FCM = E;

    % ----------------------------------------------------------------------
    % PCM
    init_V_PCM = V_FCM;
    %w = Yf_PCMC1_FindWeights1 (Xin, U, V, m, 1)
    w = [.35 .15 .35];
    [V,T,E] = Yf_PCMC1 (Xin, nC, w, [m; max_iter; term_thr; info_display; 1], init_V_PCM);

    V_PCM = V;
    T_PCM = T;
    E_PCM = E;
    
    % ----------------------------------------------------------------------
    % FPCM
    [V,U,T,E] = Yf_FPCMC1 (Xin, nC, [m; eta; max_iter; term_thr; info_display; 1], init_V);

    V_FPCM = V;
    U_FPCM = U;
    T_FPCM = T;
    E_FPCM = E;

    Iter_FCM(i) = length (E_FCM);
    Iter_PCM(i) = length (E_PCM);
    Iter_FPCM(i) = length (E_FPCM);

    ResubE_FCM(i) = Tmpf_FindResubErrorIris1 (U_FCM);
    ResubE_PCM(i) = Tmpf_FindResubErrorIris1 (T_PCM);
    ResubE_FPCM_U(i) = Tmpf_FindResubErrorIris1 (U_FPCM);
    ResubE_FPCM_T(i) = Tmpf_FindResubErrorIris1 (T_FPCM);
    
    %eta = eta + 0.5;
    m = m + .5;
end

Iter_FCM
Iter_PCM
Iter_FPCM

ResubE_FCM
ResubE_PCM
ResubE_FPCM_U
ResubE_FPCM_T
% ----------------------------------------------------------------------
% ----------------------------------------------------------------------
% ----------------------------------------------------------------------