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
%在坐標系內隨機分布50個待定位節點，以紅色星號表示
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;
%以dm為半徑，尋找每個待定位節點周圍的錨節點
for?i=1:1:m
d=norm（（xy（:i）-SS）2）;
%%%%%%%%%%%設定權值%%%%%%%%%%%%

if?d<=dm
%xy是一個2行m列的矩陣，所有橫坐標為第一行，縱坐標為第二行；xx（j，i）
%表示第j個待定位節點的第i個錨節點的橫坐標
p（ji）=1/d?;
xx（ji）=xy（1i）;
yy（ji）=xy（2i）;
k=k+1;
else
p（ji）=0;??
xx（ji）=0;
yy（ji）=0;
end
end
%計算范圍內所有錨節點的平均坐標，即質心，以空心圓表示
if?k~=0
cent（:j）=[sum（xx（j:））;sum（yy（j