Ritter's求最小包围圆

Ritter's算法如下：

1.从点集中随机选出两个点作为直径对圆进行初始化。

2.判断下一个点p是否在圆中，如果在则继续本步骤，如果不在则进行步骤3。

3.使用p作为新圆的一个边界点，另一个边界点为距离p最远的圆上的点（旧圆心和新的点p连线方向的向量，直径d+r），使用这两个点作为直径构造新圆。

4.继续步骤2，直到遍历完所有点。

n=100;p=rand(n,2); %100个点的坐标，二维坐标p1=p(1,:);p2=p(2,:);r=sqrt((p1(1)-p2(1))^2+(p1(2)-p2(2))^2)/2;  %两点距离的一半是半径cenp=(p1+p2)/2;  %圆心的坐标for i=3:n    newp=p(i,:);        d=sqrt((cenp(1)-newp(1))^2+(cenp(2)-newp(2))^2);      if d>r        r=(r+d)/2;        cenp=cenp+(d-r)/d*(newp-cenp);    %(newp-cenp)/d是这个方向上的单位向量，这个新的圆不是最优的，是直径d+原先r的，新直径就是r=(r+d)/2    end    endhold on;plot(p(:,1),p(:,2),'o');x0=cenp(1);  %迭代出来的最后的圆心坐标y0=cenp(2);  %迭代出来的最后的圆心坐标theta=0:0.01:2*pi;  %参数theta，用来极坐标画圆圈用x=x0+r*cos(theta);y=y0+r*sin(theta);plot(x,y,'-',x0,y0,'.');  %把圆周上的点用-直线段依次连接起来，正多边形无限逼近圆，步长0.01很小，最终的圆心（x0,y0）用一点句号点标出来axis equal   %横纵坐标刻度是一样比例的，防止圆变形成椭圆。axis square命令意思是画出来正方形的坐标图