clear syms x y s(10)=0; % 已知量 xs(2,10)=0; ys(2,10)=0; for n=1:10 [sx,sy]=solve( ((x+s(n)/2)^2+(y-50)^2)^0.5==60-s(n),(y-50)/(x-s(n)/2)==25/(s(n)/2-(((160-s(1))^2-50^2)^0.5+s(1)))/2 ); xs(:,n)=sx;ys(:,n)=sy; end xs,ys