可以这么干,前提是d>r1+r2:
d = sqrt((x2-x1)^2+(y2-y2)^2); dm = (d-r1-r2)/2+r1;%中心点到r1圆的距离 xm = x1 + dm*(x2-x1)/d; % 中点x坐标 ym = y1 + dm*(y2-y1)/d; % 中点y坐标