matlab实现中点画圆法1/8循环

中点画圆法是一种经典的画圆算法，它可以在计算机屏幕上快速地画出圆形。以下是 Matlab 实现中点画圆法的代码，其中使用了 1/8 循环，可以提高计算效率：

function draw_circle(xc, yc, r)
% xc, yc 为圆心坐标，r 为半径
x = 0;
y = r;
p = 5/4 - r;
plot_points(xc, yc, x, y); % 绘制第一个点和对称点
while x < y
    x = x + 1;
    if p < 0
        p = p + 2 * x + 1;
    else
        y = y - 1;
        p = p + 2 * (x - y) + 1;
    end
    plot_points(xc, yc, x, y); % 绘制其他点和对称点
end

function plot_points(xc, yc, x, y)
% 绘制圆上对称点
plot(xc + x, yc + y, '.');
plot(xc + y, yc + x, '.');
plot(xc - y, yc + x, '.');
plot(xc - x, yc + y, '.');
plot(xc - x, yc - y, '.');
plot(xc - y, yc - x, '.');
plot(xc + y, yc - x, '.');
plot(xc + x, yc - y, '.');

在这段代码中，我们首先定义了一个 `draw_circle` 函数，它接受圆心坐标和半径作为参数。在函数中，我们定义了圆上的第一个点 `(0, r)`，然后使用一个 while 循环不断计算出其他点，并调用 `plot_points` 函数绘制这些点及其对称点。

而 `plot_points` 函数则是负责绘制圆上的点和对称点。我们通过将圆心坐标和当前点的坐标相加或相减，就可以得到这些点的对称点。最后，我们使用 Matlab 的 `plot` 函数来绘制这些点。

需要注意的是，由于 Matlab 中的坐标系原点在左下角，因此我们需要将绘制出的圆进行翻转，才能得到正确的图形。