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