博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
matlab练习程序(Ritter‘s最小包围圆)
阅读量:6266 次
发布时间:2019-06-22

本文共 720 字,大约阅读时间需要 2 分钟。

原始算法是sphere,我这里简化为circle了。

Ritter's求最小包围圆为线性算法,因为非常简单,所以应用非常广泛。

该算法求出的圆比最优圆大概会大个5%到20%左右,求最优圆应该可以用Bouncing Bubble算法,以后有机会可以尝试一下。

Ritter's算法如下:

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

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

3.使用p作为新圆的一个边界点,另一个边界点为距离p最远的圆上的点,使用这两个点作为直径构造新圆。

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

结果如下:

matlab代码如下:

clear all;close all;clc;n=100;p=rand(n,2);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);    end    endhold on;plot(p(:,1),p(:,2),'o');x0=cenp(1);y0=cenp(2);theta=0:0.01:2*pi;x=x0+r*cos(theta);y=y0+r*sin(theta);plot(x,y,'-',x0,y0,'.');axis equal

参考:

转载地址:http://pidpa.baihongyu.com/

你可能感兴趣的文章
浏览器 UA 判断
查看>>
理解OAuth 2.0
查看>>
高并发处理思路与手段(三):消息队列
查看>>
Docker+Nginx部署Angular
查看>>
Docker & ASP.NET Core (4):容器间的连接
查看>>
beam 的异常处理 Error Handling Elements in Apache Beam Pipelines
查看>>
将png图片转换为字体图标
查看>>
/var/log/wtmp
查看>>
C# 获取机器码
查看>>
什么是医嘱?医嘱的书写内容?
查看>>
如何通过CSP编程卸载Windows Mobile应用程序
查看>>
用delphi实现完美屏幕截图
查看>>
matlab练习程序(差异演化DE)
查看>>
这就是搜索引擎:核心技术详解
查看>>
加解密技术处理时间对比
查看>>
g++命令行详解 (转)
查看>>
Ubuntu菜鸟入门(九)—— 支付宝支付控件安装
查看>>
什么是 SRS 呢?在我们大部分的音频播放器里都内欠有这种音效。
查看>>
对/etc/rc.d/init.d目录的一点理解(转)
查看>>
c#使用params重载方法
查看>>