MATLAB科赫曲线

本程序生成分形图案——Cantor曲线

修改迭代次数（程序中是4），可使图案更精细化。

clear

clc

x=0.1:0.2:10; %初始点集的x坐标

y=zeros(1,50); %初始点集的y坐标

length_0=length(x);

L=[x;y]; %L的没列代表一个点的坐标，总的列数等于点数A1=1/3*[1,0;0,1]; %相乘项只需保持2列，不必扩展匹配

%W1的IFS码中无t1项

A2=1/3*[1,0;0,1];

t2=10*[2/3;0]*ones(1,length_0);

for i=1:4 %迭代次数。次数越大，Koch曲线越复杂m1=A1*L;

x1=m1(1,:);

y1=m1(2,:);

m2=A2*L+t2;

x2=m2(1,:);

y2=m2(2,:);

L=[x1 x2;y1 y2]; %形成新的部分点集的组合

length_0=length_0*2;

t2=10*[2/3;0]*ones(1,length_0);

end

plot(L(1,:),L(2,:),'.')

axis([-1,11,-1,11])