质心算法matlab仿真程序

这是一个能运行的质心算法matlab仿真程序

clear
clc
for i=1:1:10
    for j=1:1:10
        x（j+（i-1）*10）=（i-1）*10;
        y（j+（i-1）*10）=（j-1）*10;
    end
end
figure
plot（xy‘.‘）
hold on 
axis（[0 100 0 100]）
xy=[x;y]
hold on
xm=90;
ym=90; 
n=50;   
for i=1:1:n
   Sx（i）=rand（11）*xm;
   Sy（i）=rand（11）*ym;
    plot（Sx（i）Sy（i）‘r*‘）
    xlabel（‘x轴‘）
    ylabel（‘y轴‘）
   hold on
end
dm=30  
m=100;  
for j=1:1:n
    SS=[Sx（j）;Sy（j）]; 
    k=0;
    for i=1:1:m
        d=norm（（xy（:i）-SS）2）; 
        if d<=dm  
            xx（ji）=xy（1i）;
            yy（ji）=xy（2i）;
            k=k+1;
        else
            xx（ji）=0;
            yy（ji）=0;
        end
    end
    if k~=0
        cent（:j）=[sum（xx（j:））;sum（yy（j:））]/k;
    else
        cent（:j）=0;
    end
 plot（cent（1j）cent（2j）‘o‘）  
hold on
 plot（[cent（1j） Sx（j）][cent（2j） Sy（j）]‘R‘）  
title（‘Centroid‘）
hold on
MM=[cent（1j）;cent（2j）]
e（j）=norm（（MM-SS）2）/dm
end
figure
axis（[0 n 0 1]）
   j=1:1:n
    plot（je（j） ‘-r.‘）
   hold on
   title（‘Centroid‘）
   E=sum（e）/n

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1059  2010-07-28 17:27  Centroid.m

----------- ---------  ---------- -----  ----

                 1059                    1