画一个圆，很容易。

发布的布雷特Shoelson，2008年3月14日,

123次观看(过去30天)|0喜欢|19日的评论

你可能已经注意到这个博客的格式最近发生了一些变化。以下是每周两个主题的常规内容。

周二Doug将提供MATLAB教程。
周五的客座博主Jiro, Brett和Bob将重点介绍File Exchange提交的内容。

一个文件不需要很长才有用。布雷特这周的选择，伊扎克·布彻圆适合，不包括注释，只有5行。但我真的很喜欢Izhak的词条，在我下载它之后的几年里，我有很多机会使用它。不知何故，在我的工作中，圆与点的拟合要求似乎以令人费解的频率出现。伊扎克的CIRCFIT很容易做到这一点，而且在大多数情况下都相当强大。那么如何使用伊扎克函数呢?考虑以下为说明目的而创建的噪声数据:

R = 10;x_c = 5;y_c = 8;θ= 0:π/ 64:π;xs = x_c + R*cos();ys = y_c + R*sin();现在添加一些随机噪声，使问题更具挑战性:乘= 0.5;x = x +乘* randn(大小(x));y = y +乘* randn(大小(y));

情节点……

图绘制(x, y,“b”。)轴平等的

.．.然后计算并显示最佳拟合圆

[xfit, yfit Rfit] = circfit (x, y);图绘制(x, y,“b”。)举行在矩形(“位置”(xfit-Rfit yfit-Rfit Rfit * 2, Rfit * 2),.．.“弯曲”[1],“线型”，“- - -”，“edgecolor”，“r”）;标题(sprintf ('最佳拟合:R = %0.1f;Ctr = (% 0.1 f, % 0.1度)，.．.Rfit、xfit yfit));情节(xfit yfit,“g”。([xfit-Rfit-2,xfit+Rfit+2]) ylim([yfit-Rfit-2,yfit+Rfit+2])轴平等的