clear all
close all
load TRAINLPC
s=[];
for i=1:1:3
    s=[s size(TRAINLPC{i},1)];
end
mx=min(s);
DATA=[];
for i=1:1:3
    DATA=[DATA TRAINLPC{i}(1:1:mx,:)];
end
DATA1=[];
for i=1:1:size(DATA,2)
    temp=DATA(:,i);
DATA1=[DATA1 reshape(temp,13,14)];
end
[E,D]=eigs(cov(DATA1'),3);
DATA1=E'*DATA1;
[M,COV,INDEX,m,D,J]=fuzzykmeansforspeech(DATA1,4,200)