Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- matriz=[5.1,3.5,1.4,0.2,0
- 4.9,3.0,1.4,0.2,0
- 4.7,3.2,1.3,0.2,0
- 4.6,3.1,1.5,0.2,0
- 5.0,3.6,1.4,0.2,0
- 5.4,3.9,1.7,0.4,0
- 4.6,3.4,1.4,0.3,0
- 5.0,3.4,1.5,0.2,0
- 4.4,2.9,1.4,0.2,0
- 4.9,3.1,1.5,0.1,0
- 5.4,3.7,1.5,0.2,0
- 4.8,3.4,1.6,0.2,0
- 4.8,3.0,1.4,0.1,0
- 4.3,3.0,1.1,0.1,0
- 5.8,4.0,1.2,0.2,0
- 5.7,4.4,1.5,0.4,0
- 5.4,3.9,1.3,0.4,0
- 5.1,3.5,1.4,0.3,0
- 5.7,3.8,1.7,0.3,0
- 5.1,3.8,1.5,0.3,0
- 5.4,3.4,1.7,0.2,0
- 5.1,3.7,1.5,0.4,0
- 4.6,3.6,1.0,0.2,0
- 5.1,3.3,1.7,0.5,0
- 4.8,3.4,1.9,0.2,0
- 5.0,3.0,1.6,0.2,0
- 5.0,3.4,1.6,0.4,0
- 5.2,3.5,1.5,0.2,0
- 5.2,3.4,1.4,0.2,0
- 4.7,3.2,1.6,0.2,0
- 4.8,3.1,1.6,0.2,0
- 5.4,3.4,1.5,0.4,0
- 5.2,4.1,1.5,0.1,0
- 5.5,4.2,1.4,0.2,0
- 4.9,3.1,1.5,0.1,0
- 5.0,3.2,1.2,0.2,0
- 5.5,3.5,1.3,0.2,0
- 4.9,3.1,1.5,0.1,0
- 4.4,3.0,1.3,0.2,0
- 5.1,3.4,1.5,0.2,0
- 5.0,3.5,1.3,0.3,0
- 4.5,2.3,1.3,0.3,0
- 4.4,3.2,1.3,0.2,0
- 5.0,3.5,1.6,0.6,0
- 5.1,3.8,1.9,0.4,0
- 4.8,3.0,1.4,0.3,0
- 5.1,3.8,1.6,0.2,0
- 4.6,3.2,1.4,0.2,0
- 5.3,3.7,1.5,0.2,0
- 5.0,3.3,1.4,0.2,0
- 7.0,3.2,4.7,1.4,1
- 6.4,3.2,4.5,1.5,1
- 6.9,3.1,4.9,1.5,1
- 5.5,2.3,4.0,1.3,1
- 6.5,2.8,4.6,1.5,1
- 5.7,2.8,4.5,1.3,1
- 6.3,3.3,4.7,1.6,1
- 4.9,2.4,3.3,1.0,1
- 6.6,2.9,4.6,1.3,1
- 5.2,2.7,3.9,1.4,1
- 5.0,2.0,3.5,1.0,1
- 5.9,3.0,4.2,1.5,1
- 6.0,2.2,4.0,1.0,1
- 6.1,2.9,4.7,1.4,1
- 5.6,2.9,3.6,1.3,1
- 6.7,3.1,4.4,1.4,1
- 5.6,3.0,4.5,1.5,1
- 5.8,2.7,4.1,1.0,1
- 6.2,2.2,4.5,1.5,1
- 5.6,2.5,3.9,1.1,1
- 5.9,3.2,4.8,1.8,1
- 6.1,2.8,4.0,1.3,1
- 6.3,2.5,4.9,1.5,1
- 6.1,2.8,4.7,1.2,1
- 6.4,2.9,4.3,1.3,1
- 6.6,3.0,4.4,1.4,1
- 6.8,2.8,4.8,1.4,1
- 6.7,3.0,5.0,1.7,1
- 6.0,2.9,4.5,1.5,1
- 5.7,2.6,3.5,1.0,1
- 5.5,2.4,3.8,1.1,1
- 5.5,2.4,3.7,1.0,1
- 5.8,2.7,3.9,1.2,1
- 6.0,2.7,5.1,1.6,1
- 5.4,3.0,4.5,1.5,1
- 6.0,3.4,4.5,1.6,1
- 6.7,3.1,4.7,1.5,1
- 6.3,2.3,4.4,1.3,1
- 5.6,3.0,4.1,1.3,1
- 5.5,2.5,4.0,1.3,1
- 5.5,2.6,4.4,1.2,1
- 6.1,3.0,4.6,1.4,1
- 5.8,2.6,4.0,1.2,1
- 5.0,2.3,3.3,1.0,1
- 5.6,2.7,4.2,1.3,1
- 5.7,3.0,4.2,1.2,1
- 5.7,2.9,4.2,1.3,1
- 6.2,2.9,4.3,1.3,1
- 5.1,2.5,3.0,1.1,1
- 5.7,2.8,4.1,1.3,1
- 6.3,3.3,6.0,2.5,2
- 5.8,2.7,5.1,1.9,2
- 7.1,3.0,5.9,2.1,2
- 6.3,2.9,5.6,1.8,2
- 6.5,3.0,5.8,2.2,2
- 7.6,3.0,6.6,2.1,2
- 4.9,2.5,4.5,1.7,2
- 7.3,2.9,6.3,1.8,2
- 6.7,2.5,5.8,1.8,2
- 7.2,3.6,6.1,2.5,2
- 6.5,3.2,5.1,2.0,2
- 6.4,2.7,5.3,1.9,2
- 6.8,3.0,5.5,2.1,2
- 5.7,2.5,5.0,2.0,2
- 5.8,2.8,5.1,2.4,2
- 6.4,3.2,5.3,2.3,2
- 6.5,3.0,5.5,1.8,2
- 7.7,3.8,6.7,2.2,2
- 7.7,2.6,6.9,2.3,2
- 6.0,2.2,5.0,1.5,2
- 6.9,3.2,5.7,2.3,2
- 5.6,2.8,4.9,2.0,2
- 7.7,2.8,6.7,2.0,2
- 6.3,2.7,4.9,1.8,2
- 6.7,3.3,5.7,2.1,2
- 7.2,3.2,6.0,1.8,2
- 6.2,2.8,4.8,1.8,2
- 6.1,3.0,4.9,1.8,2
- 6.4,2.8,5.6,2.1,2
- 7.2,3.0,5.8,1.6,2
- 7.4,2.8,6.1,1.9,2
- 7.9,3.8,6.4,2.0,2
- 6.4,2.8,5.6,2.2,2
- 6.3,2.8,5.1,1.5,2
- 6.1,2.6,5.6,1.4,2
- 7.7,3.0,6.1,2.3,2
- 6.3,3.4,5.6,2.4,2
- 6.4,3.1,5.5,1.8,2
- 6.0,3.0,4.8,1.8,2
- 6.9,3.1,5.4,2.1,2
- 6.7,3.1,5.6,2.4,2
- 6.9,3.1,5.1,2.3,2
- 5.8,2.7,5.1,1.9,2
- 6.8,3.2,5.9,2.3,2
- 6.7,3.3,5.7,2.5,2
- 6.7,3.0,5.2,2.3,2
- 6.3,2.5,5.0,1.9,2
- 6.5,3.0,5.2,2.0,2
- 6.2,3.4,5.4,2.3,2
- 5.9,3.0,5.1,1.8,2];
- c1=[5.1,3.5,1.4,0.2];
- c2=[4.9,3.0,1.4,0.2];
- c3=[4.7,3.2,1.3,0.2];
- converg1=[1,1,1,1];
- converg2=[1,1,1,1];
- converg3=[1,1,1,1];
- x = 0;
- while converg1~=[0,0,0,0] & converg2~=[0,0,0,0] & converg2~=[0,0,0,0] do
- x = x + 1;
- classe1 = zeros(1,4);
- classe2 = zeros(1,4);
- classe3 = zeros(1,4);
- for i=1:150
- d1 = sqrt((matriz(i:i,1:1)-c1(1:1,1:1))^2 + (matriz(i:i,2:2)-c1(1:1,2:2))^2 + (matriz(i:i,3:3)-c1(1:1,3:3))^2 + (matriz(i:i,4:4)-c1(1:1,4:4))^2);
- d2 = sqrt((matriz(i:i,1:1)-c2(1:1,1:1))^2 + (matriz(i:i,2:2)-c2(1:1,2:2))^2 + (matriz(i:i,3:3)-c2(1:1,3:3))^2 + (matriz(i:i,4:4)-c2(1:1,4:4))^2);
- d3 = sqrt((matriz(i:i,1:1)-c3(1:1,1:1))^2 + (matriz(i:i,2:2)-c3(1:1,2:2))^2 + (matriz(i:i,3:3)-c3(1:1,3:3))^2 + (matriz(i:i,4:4)-c3(1:1,4:4))^2);
- if d1<=d2 & d1<=d3
- classe1 = cat(1,classe1,matriz(i:i,1:4));
- elseif d2<=d1 & d2<d3
- classe2 = cat(1,classe2,matriz(i:i,1:4));
- elseif d3<d1 & d3<d2
- classe3 = cat(1,classe3,matriz(i:i,1:4));
- end;
- end;
- if isempty(classe1)
- classe1 = [0,0,0,0];
- elseif isempty(classe2)
- classe2 = [0,0,0,0];
- elseif isempty(classe3)
- classe3 = [0,0,0,0];
- end
- [m1,n1] = size(classe1);
- [m2,n2] = size(classe2);
- [m3,n3] = size(classe3);
- converg1 = c1 - mean(classe1(2:m1,:),'r');
- converg2 = c2 - mean(classe2(2:m2,:),'r');
- converg3 = c3 - mean(classe3(2:m3,:),'r');
- c1 = mean(classe1(2:m1,:),'r');
- c2 = mean(classe2(2:m2,:),'r');
- c3 = mean(classe3(2:m3,:),'r');
- setosa = classe1(2:m1,:);
- versicolor = classe2(2:m2,:);
- virginica = classe3(2:m3,:);
- end
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement