function qiujlt(x,method)
close all
A=biaozh(x);
B=qiumhjz(A,method);
R=qiubb(B);
lamda=qiusp(R);
m=length(lamda);
for j=1:m
M(:,:,i)=R>=lamda(i);
cls=qiucls(M(:,:,i));
allcls{i}=cls;
end
rr=[];
for j=1:length(allcls)-1
ind=[];
for j=1:i
ind0=[];
for k=1:i+1
ind0=[ind0 ismember(allcls{i+1}{k}(1),allcls{i}{j})];
end
ind=[ind;ind0];
end
tmp=allcls{i+1};
for 1:i
if sum(ind(r,:))==2
tp=ind(r,:);
s=find(tp==1);
allcls{i+1}{r}=tmp{s(1)};
allcls{i+1}{r+1}=tmp{s(2)};
t1=[1:r-1,r+1:i];
t2=[1:r-1,r+2:i+1];
fot t=1:i-1
allcls{i+1}{t2(t)}=allcls{i}{t1(t)};
end
rr=[rr r];
end
end
end
seq0=allcls{end};seq=[];
for i=1:m
seq=[seq allcls{end}{i}];
end
figure,hold on,title('动态聚类图')
wide=30;high=25;
x0=wide*(m-1);y0=high*m;par=(wide+high)/2;
axis([-2,5*par x0+1.5par-par y0+par])
axis('off')
lbx{l}=wide*[0:m-1];
for i=2:m
lbx{i}(rr(m+1-i))=(lbx{i-1}(rr(m+1-i))+lbx{i-1}(rr(m+1-i)+1))/2;
for u=1:rr(m+1-i)-1
lbx{i}(u)=lbx{i-1}(u);
end
for v=rr(m+1-i)+1:m+1-i
lbx{i}(v)=lbx{i-1}(v+1);
end
for i=1:m
text((i-1)*wide-wide/9,-0.5*high,int2str(allcls{end}{i}))
end
text(-3.45*par,(m+0.5)*high,'分类水平')
text(-3.1*par,0.5*high,num2str(lamda(m)))
text(x0+0.3*par,(m+0.5)*high,'分类数')
text(x0+1.35*par,0.5*high,num2str(m))
for 1:m
line([(i-1)*wide (i-1)*wide],[0 high])
end
for 2:m
text(-3.1*par,(i-0.5)*high,num2str(lamda(m+1-i)))
text(x0+1.35*par,(i-0.5)*high,num2str(m+1-i))
line([lbx{i-1}(rr(m+1-i)) lbx{i-1}(rr(m+1-i)+1)],[i-1 i-1]*high)
for j=1:length(lbx{i})
line([lbx{i}(j) lbx{i}(j)],[i-1 i]*high)
end
end
