如何在Matlab中绘制一个圆？ （最小二乘法）

Question

我试图在 x 和 y 上绘制圆的方程回归，但我不知道如何进行。 有什么建议？ （我想要一个圆通过最小二乘解决方案连接点）

x = [5; 4; -1; 1];
y = [3; 5; 2; 1];

% circle's equation: x^2+y^2 = 2xc1+2yc2+c3
a = [2.*x,2.*y,ones(n,3)]
b = [x.^2 + y.^2];
c = a\b;

在此之后我如何绘制圆圈

Answer 1

在 matlab 中有两种绘制圆的方法：

• 绘制一条线，其中数据点形成一个圆圈
• 使用绘图中的'o'标记和'MarkerSize'名称-值对来设置圆的半径
• 您可以使用vscircle 函数绘制圆形图像在您的情况下，我会选择第一个选项，因为您可以控制圆形大小。
• 使用rectangle(...,'Curvature',[1 1]) 函数[编辑：感谢@Cris Luengo]

所以这是一个绘图函数

function circle(x,y,r)
%x and y are the coordinates of the center of the circle
%r is the radius of the circle
%0.01 is the angle step, bigger values will draw the circle faster but
%you might notice imperfections (not very smooth)
ang=0:0.01:2*pi+.01; 
xp=r*cos(ang);
yp=r*sin(ang);

plot(x+xp,y+yp);
end

所以用你的（更正的）代码，它看起来像这样

x = [5; 4; -1; 1];
y = [3; 5; 2; 1];

% circle's equation:  x^2+y^2 = 2xc1+2yc2+c3
a = [2.*x,2.*y,ones(length(x),1)];
b = [x.^2 + y.^2];
c = a\b;

x_m = c(1)/2;
y_m = c(2)/2;
r = sqrt(x_m^2 + y_m^2 -c(3));

% plot data points
plot(x,y,'o')
hold on
% plot center
plot(x_m,y_m,'+')
% plot circle
circle(x_m,y_m,r)
hold off