将dj0输出的13个数据导成矩阵即可 后面还需要分别使用13个数进行计算
A=[100 90 5 1200;110 80 5 800;120 30 2 600;130 10 2 450;85 40 5 900;80 80 2 500;95 100 2 400;83 40 2 560;125 8 5 870;75 101 2 490;85 10 2 560;50 70 2 550;30 40 5 850]
n=0
x=0
y=0
for j=1:13
n1=A(j,3)*A(j,4)
n=n+n1
x0=A(j,1)*A(j,3)*A(j,4)
x=x+x0
y0=A(j,2)*A(j,3)*A(j,4)
y=y+y0
end
xd0=x/n
yd0=y/n
t0=0
for j=1:13
DJ0=zeros(13,1);
for i=1:13
dj0=sqrt((xd0-A(j,1))^2+(yd0-A(j,2))^2)
DJ0(i)=dj0;
end
end