用matlab编写koch雪花的代码

以下是用 MATLAB 编写 Koch 雪花的代码：

function koch_snowflake(n)
% n: 迭代次数

% 初始三角形的三个顶点
p1 = [0, 0];
p2 = [1, 0];
p3 = [cos(pi/3), sin(pi/3)];

% 画初始三角形
line([p1(1), p2(1)], [p1(2), p2(2)]);
line([p2(1), p3(1)], [p2(2), p3(2)]);
line([p3(1), p1(1)], [p3(2), p1(2)]);

% 递归绘制
for i = 1:n
    % 计算当前三角形的三个顶点
    p4 = (2 * p1 + p2) / 3;
    p5 = (p1 + p2) / 2 + (p3 - p1) * sqrt(3) / 6;
    p6 = (p1 + 2 * p2) / 3;
    
    % 画当前三角形的三条线段
    line([p4(1), p5(1)], [p4(2), p5(2)]);
    line([p5(1), p6(1)], [p5(2), p6(2)]);
    line([p6(1), p4(1)], [p6(2), p4(2)]);
    
    % 更新三角形顶点
    p1 = p4;
    p2 = p6;
    p3 = p5;
end

% 闭合三角形
line([p1(1), p2(1)], [p1(2), p2(2)]);
line([p2(1), p3(1)], [p2(2), p3(2)]);
line([p3(1), p1(1)], [p3(2), p1(2)]);

在 MATLAB 中执行 `koch_snowflake(5)`，可以绘制迭代 5 次的 Koch 雪花。你可以尝试调整迭代次数，观察 Koch 雪花的变化。