spe.m

公司分拆问题: 聚类算法与图论算法 matlab 的实现

function [ts,xsr,xsb,sums,tk,xkr,xkb,sumk]= spe(filename)
x0 = dlmread(filename);
x1 = x0(:,[2,3]);

ts=cputime;
dm = dist(x1');
m=mean(mean(dm));
W=exp(-dm.^2/m^2);
D = diag(sum(W));
I = eye(size(D));
L=I-D^(-1)*W;
[v,d]=eig(L);
%f=sign(v(:,2));
xsr = x1((v(:,2)>0),:);
xsb = x1((v(:,2)<=0),:);
figure(1);
plot(xsr(:,2),-xsr(:,1),'.r',xsb(:,2),-xsb(:,1),'.b')
sums=innersum(xsr)+innersum(xsb);
ts= cputime -ts;

tk=cputime;
[IDX,C,sumk] = kmeans(x1,2,'replicates',5);
xkr = x1((IDX==1),:);
xkb = x1((IDX==2),:);
figure(2);
plot(xkr(:,2),-xkr(:,1),'.r',xkb(:,2),-xkb(:,1),'.b')
sumk=sum(sumk);
tk= cputime -tk;