我正在尝试为n阶Koch分形建立代码,定义为:
为此,我在Matlab中实现了一个代码.我已经能够开发出能够绘制二阶分形的代码(零和第一阶很简单).为此,我在这篇文章中基于代码.到目前为止,我已经有了这个,它有效,但显然只是为了第二阶段(或者一个零,只需编辑一下代码)
clear all;
clc;
close all;
hold on;
% Counter-clockwise rotation matrix
R = [cosd(60) -sind(60);sind(60) cosd(60)];
angles = [0 60 -60 0 0]; % main angles
n = 2; % order
% Length of each straight segment
seglength = 3^-(n-1);
% Initialization of variables
a=0;
b=0;
c=1/3;
d=0;
for i=1:4^(n-2)
for j=1:4
p1 = [a;b];
p5 = [c;d];
p2 = (2.*p1+p5)/3;
p4 = (2.*p5+p1)/3;
p3 = p2 + R*(p4-p2);
x = [p1(1) p2(1) p3(1) p4(1) p5(1)]; …Run Code Online (Sandbox Code Playgroud)