function iriscloud
global Y1 Y2 Y3;
for i=1:50
Y1(i)=1;
Y2(i)=1;
Y3(i)=1;
end
n=300;
x1(1,:)=[5.1,4.9,4.7,4.6,5.0,5.4,4.6,5.0,4.4,4.9,5.4,4.8,4.8,4.3,5.8,5.7,5.4,5.1,5.7,5.1,5.4,5.1,4.6,5.1,4.8,5.0,5.0,5.2,5.2,4.7,4.8,5.4,5.2,5.5,4.9,5.0,5.5,4.9,4.4,5.1,5.0,4.5,4.4,5.0,5.1,4.8,5.1,4.6,5.3,5.0];
x1(2,:)=[3.5,3.0,3.2,3.1,3.6,3.9,3.4,3.4,2.9,3.1,3.7,3.4,3.0,3.0,4.0,4.4,3.9,3.5,3.8,3.8,3.4,3.7,3.6,3.3,3.4,3.0,3.4,3.5,3.4,3.2,3.1,3.4,4.1,4.2,3.1,3.2,3.5,3.1,3.0,3.4,3.5,2.3,3.2,3.5,3.8,3.0,3.8,3.2,3.7,3.3];
x1(3,:)=[1.4,1.4,1.3,1.5,1.4,1.7,1.4,1.5,1.4,1.5,1.5,1.6,1.4,1.1,1.2,1.5,1.3,1.4,1.7,1.5,1.7,1.5,1.0,1.7,1.9,1.6,1.6,1.5,1.4,1.6,1.6,1.5,1.5,1.4,1.5,1.2,1.3,1.5,1.3,1.5,1.3,1.3,1.3,1.6,1.9,1.4,1.6,1.4,1.5,1.4];
x1(4,:)=[0.2,0.2,0.2,0.2,0.2,0.4,0.3,0.2,0.2,0.1,0.2,0.2,0.1,0.1,0.2,0.4,0.4,0.3,0.3,0.3,0.2,0.4,0.2,0.5,0.2,0.2,0.4,0.2,0.2,0.2,0.2,0.4,0.1,0.2,0.1,0.2,0.2,0.1,0.2,0.2,0.3,0.3,0.2,0.6,0.4,0.3,0.2,0.2,0.2,0.2];
ex1(1)=mean(x1(1,:)');en1(1)=std(x1(1,:)');
ex1(2)=mean(x1(2,:)');en1(2)=std(x1(2,:)');
ex1(3)=mean(x1(3,:)');en1(3)=std(x1(3,:)');
ex1(4)=mean(x1(4,:)');en1(4)=std(x1(4,:)');
subplot(2,2,1)
plot3(x1(1,:),x1(3,:),x1(4,:),'Marker','*','Markersize',2,'linestyle','no','color','k')
axis([4 8 1.1 6.7 0 2.5])
title('植物数据','fontsize',8)
text(ex1(1)+1,ex1(3),ex1(4)+0.3,'\leftarrowSetosa','FontSize',8)
xlabel('萼片长度/(cm)','fontsize',8),ylabel('花瓣高度/(cm)','fontsize',8),zlabel('花瓣宽度/(cm)','fontsize',8)
grid on
hold on
x2(1,:)=[7.0,6.4,6.9,5.5,6.5,5.7,6.3,4.9,6.6,5.2,5.0,5.9,6.0,6.1,5.6,6.7,5.6,5.8,6.2,5.6,5.9,6.1,6.3,6.1,6.4,6.6,6.8,6.7,6.0,5.7,5.5,5.5,5.8,6.0,5.4,6.0,6.7,6.3,5.6,5.5,5.5,5.5,5.8,6.0,5.4,6.0,6.7,6.3,5.6,5.5];
x2(2,:)=[3.2,3.2,3.1,2.3,2.8,2.8,2.3,2.4,2.9,2.7,2.0,3.0,2.2,2.9,2.9,3.1,3.0,2.7,2.2,2.5,3.2,2.8,2.5,2.8,2.9,3.0,2.8,3.0,2.9,2.6,2.4,2.4,2.7,2.7,3.0,3.4,3.1,2.3,3.0,2.5,2.6,3.0,2.6,2.3,2.7,3.0,2.9,2.9,2.5,2.8];
x2(3,:)=[4.7,4.5,4.9,4.0,4.6,4.5,4.7,3.3,4.6,3.9,3.5,4.2,4.0,4.7,3.6,4.4,4.5,4.1,4.5,3.9,4.8,4.0,4.9,4.7,4.3,4.4,4.8,5.0,4.5,3.5,3.8,3.7,3.9,5.1,4.5,4.5,4.7,4.4,4.1,4.0,4.4,4.6,4.0,3.3,4.2,4.2,4.2,4.3,3.0,4.1];
x2(4,:)=[1.4,1.5,1.5,1.3,1.5,1.3,1.6,1.0,1.3,1.4,1.0,1.5,1.0,1.4,1.3,1.4,1.5,1.0,1.5,1.1,1.8,1.3,1.5,1.2,1.3,1.4,1.4,1.7,1.5,1.0,1.1,1.0,1.2,1.6,1.5,1.6,1.5,1.3,1.3,1.3,1.2,1.4,1.2,1.0,1.3,1.2,1.3,1.3,1.1,1.3];
ex2(1)=mean(x2(1,:)');en2(1)=std(x2(1,:)');
ex2(2)=mean(x2(2,:)');en2(2)=std(x2(2,:)');
ex2(3)=mean(x2(3,:)');en2(3)=std(x2(3,:)');
ex2(4)=mean(x2(4,:)');en2(4)=std(x2(4,:)');
subplot(2,2,1)
plot3(x2(1,:),x2(3,:),x2(4,:),'Marker','*','Markersize',2,'linestyle','no','color','k')
text(ex2(1)+1,ex2(3),ex2(4),'\leftarrow Versicolour','FontSize',8)
hold on
x3(1,:)=[6.3,5.8,7.1,6.3,6.5,7.6,4.9,7.3,6.7,7.2,6.5,6.4,6.8,5.7,5.8,6.4,6.5,7.7,7.7,6.0,6.9,5.6,7.7,6.3,6.7,7.2,6.2,6.1,6.4,7.2,7.4,7.9,6.4,6.3,6.1,7.7,6.3,6.4,6.0,6.9,6.7,6.9,5.8,6.8,6.7,6.7,6.3,6.5,6.2,5.9];
x3(2,:)=[3.3,2.7,3.0,2.9,3.0,3.0,2.5,2.9,2.5,3.6,3.2,2.7,3.0,2.5,2.8,3.2,3.0,3.8,2.6,2.2,3.2,2.8,2.8,2.7,3.3,3.2,2.8,3.0,2.8,3.0,2.8,3.8,2.8,2.8,2.6,3.0,3.4,3.1,3.0,3.1,3.1,3.1,2.7,3.2,3.3,3.0,2.5,3.0,3.4,3.0];
x3(3,:)=[6.0,5.1,5.9,5.6,5.8,6.6,4.5,6.3,5.8,6.1,5.1,5.3,5.5,5.0,5.1,5.3,5.5,6.7,6.9,5.0,5.7,4.9,6.7,4.9,5.7,6.0,4.8,4.9,5.6,5.8,6.1,6.6,6.4,5.6,5.1,6.1,5.6,5.5,4.8,5.4,5.6,5.1,5.1,5.9,5.7,5.2,5.0,5.2,5.4,5.1];
x3(4,:)=[2.5,1.9,2.1,1.8,2.2,2.1,1.7,1.8,1.8,2.5,2.0,1.9,2.1,2.0,2.4,2.3,1.8,2.2,2.3,1.5,2.3,2.0,2.0,1.8,2.1,1.8,1.8,1.8,2.1,1.6,1.9,2.0,2.2,1.5,1.4,2.3,2.4,1.8,1.8,2.1,2.4,2.3,1.9,2.3,2.5,2.3,1.9,2.0,2.3,1.8];
ex3(1)=mean(x3(1,:)');en3(1)=std(x3(1,:)');
ex3(2)=mean(x3(2,:)');en3(2)=std(x3(2,:)');
ex3(3)=mean(x3(3,:)');en3(3)=std(x3(3,:)');
ex3(4)=mean(x3(4,:)');en3(4)=std(x3(4,:)');
for i=1:4
en3(i)=(max(x3(i,:))-min(x3(i,:)))/6;
en2(i)=(max(x2(i,:))-min(x2(i,:)))/6;
en1(i)=(max(x1(i,:))-min(x1(i,:)))/6;
end
cloud1(x1,ex1,en1,n)
cloud2(x1,ex2,en2,n)
cloud3(x3,ex3,en3,n);
jishu1=0;
jishu2=0;
for i=1:50
if(Y1(i)>Y2(i))
jishu1=jishu1+1;
end
if(Y1(i)>Y3(i))
jishu2=jishu2+1;
end
end
jishu1
jishu2
subplot(2,2,1)
plot3(x3(1,:),x3(3,:),x3(4,:),'Marker','.','Markersize',3,'linestyle','no','color','k')
text(ex3(1)+1,ex3(3),ex3(4),'\leftarrow Virginica','FontSize',8)
x(1,:)=[x1(1,:) x2(1,:) x3(1,:)];
x(2,:)=[x1(2,:) x2(2,:) x3(2,:)];
x(3,:)=[x1(3,:) x2(3,:) x3(3,:)];
x(4,:)=[x1(4,:) x2(4,:) x3(4,:)];
ex(1)=mean(x(1,:));en(1)=std(x(1,:));
ex(2)=mean(x(2,:));en(2)=std(x(2,:));
ex(3)=mean(x(3,:));en(3)=std(x(3,:));
ex(4)=mean(x(4,:));en(4)=std(x(4,:));
cloud(x,ex,en,n);
function f=cloud(x,ex,en,n)
for j=1:4
for i=1:n
xx(j,i)=normrnd(ex(j),en(j));
end
end
for i=1:n
lishu(i)=exp(-0.5*(xx(1,i)-ex(1))^2/en(1)^2-0.5*(xx(2,i)-ex(2))^2/en(2)^2-0.5*(xx(3,i)-ex(3))^2/en(3)^2-0.5*(xx(4,i)-ex(4))^2/en(4)^2);
end
subplot(2,2,4)
plot3(xx(1,:),xx(4,:),lishu','Marker','.','Markersize',3,'linestyle','no','color','k')
title('基云','fontsize',8)
xlabel('萼片长度/cm','fontsize',8),ylabel('花瓣宽度/cm','fontsize',8),zlabel('隶属度','fontsize',8)
hold on
success=0;
[a,b]=size(x);
for i=1:b
for j=1:a
y(j,i)=exp(-0.5*(x(j,i)-ex(j))^2/en(j)^2);
end
end
function f=cloud3(x,ex,en,n)
global Y3;
for j=1:4
for i=1:n
xx(j,i)=normrnd(ex(j),en(j));
end
end
for i=1:n
lishu(i)=exp(-0.5*(xx(1,i)-ex(1))^2/en(1)^2-0.5*(xx(2,i)-ex(2))^2/en(2)^2-0.5*(xx(3,i)-ex(3))^2/en(3)^2-0.5*(xx(4,i)-ex(4))^2/en(4)^2);
end
subplot(2,2,3)
plot3(xx(1,:),xx(4,:),lishu','Marker','.','Markersize',3,'linestyle','no','color','k')
text(ex(1),ex(4),0.8,'\leftarrowVirginica','FontSize',8)
axis([0 8 0 3 0 1])
hold on
success3=0;
[a,b]=size(x);
for i=1:b
Y3(i)=1;
for j=1:a
y(j,i)=exp(-0.5*(x(j,i)-ex(j))^2/en(j)^2);
end
Y3(i)=Y3(i)*y(j,i);
end
function f=cloud2(x,ex,en,n)
global Y2;
for j=1:4
for i=1:n
xx(j,i)=normrnd(ex(j),en(j));
end
end
for i=1:n
lishu(i)=exp(-0.5*(xx(1,i)-ex(1))^2/en(1)^2-0.5*(xx(2,i)-ex(2))^2/en(2)^2-0.5*(xx(3,i)-ex(3))^2/en(3)^2-0.5*(xx(4,i)-ex(4))^2/en(4)^2);
end
subplot(2,2,3)
plot3(xx(1,:),xx(4,:),lishu','Marker','.','Markersize',2,'linestyle','no','color','k')
text(ex(1),ex(4),0.8,'\leftarrowVersicolour','FontSize',8)
hold on
success2=0;
[a,b]=size(x);
for i=1:b
Y2(i)=1;
for j=1:a
y(j,i)=exp(-0.5*(x(j,i)-ex(j))^2/en(j)^2);
end
Y2(i)=Y2(i)*y(j,i);
end
function f=cloud1(x,ex,en,n)
global Y1;
for j=1:4
for i=1:n
xx(j,i)=normrnd(ex(j),en(j));
end
end
for i=1:n
lishu(i)=exp(-0.5*(xx(1,i)-ex(1))^2/en(1)^2-0.5*(xx(2,i)-ex(2))^2/en(2)^2-0.5*(xx(3,i)-ex(3))^2/en(3)^2-0.5*(xx(4,i)-ex(4))^2/en(4)^2);
end
subplot(2,2,3)
plot3(xx(1,:),xx(4,:),lishu','Marker','.','Markersize',3,'linestyle','no','color','k')
%title('植物三个分解云','fontsize',8)
text(ex(1),ex(4),0.8,'\leftarrowSetosa','FontSize',8)
xlabel('萼片长度/cm','fontsize',8),ylabel('花瓣宽度/cm','fontsize',8),zlabel('隶属度','fontsize',8)
axis([0 8 0 2.5 0 1])
hold on
[a,b]=size(x);
for i=1:b
Y1(i)=1;
for j=1:a
y(j,i)=exp(-0.5*(x(j,i)-ex(j))^2/en(j)^2);
end
Y1(i)=Y1(i)*y(j,i);
end
%success=0;
%PopSize=20;
%MaxIt=70;
%iter=1
iriscloud.rar （1个子文件） 
iriscloud.m 11KB
邓凌佳
- 粉丝: 79
- 资源: 1万+
下载权益
C知道特权
VIP文章
课程特权
VIP享7折，此内容立减3.57元 开通VIP
评论1