质心定位算法

clear allclc;
for n=6:2:30
x=100*rand（1100）;      %布置10m*10m的网格区域           
y=100*rand（1100）;
w=100*rand（1n）;
z=100*rand（1n）;
plot（xy‘b*‘wz‘rO‘）
axis（[0 100 0 100]）
grid on;
xlabel（‘x‘）ylabel（‘y‘）
title（‘原始点分布‘）
C=0;
X=zeros（1100）;
Y=zeros（1100）;
for  i=1:100  
       m=0; a=0; b=0;
       for k=1:n          
               dist=distance（x（i）y（i）w（k）z（k））;        %节点与锚节点间距离
               if dist<=10                                %通信半径值=2             
                   a=a+w（k）;
                   b=b+z（k）;
                   m=m+1;
               end         
       end
       if m>=2
             X（i）=a/m;          
             Y（i）=b/m; 
             %plot（[x（i）X（i）][y（i）Y（i）]）;
             %hold on;         
       else