分形图程序、分形图基本图形以及源程序、分形图形学实验指导.doc

（1）Koch曲线程序koch.m

functionkoch(a1,b1,a2,b2,n)

%koch(0,0,9,0,3)

%a1,b1,a2,b2为初始线段两端点坐标,n为迭代次数

a1=0;b1=0;a2=9;b2=0;n=3;

%第i-1次迭代时由各条线段产生的新四条线段的五点横、纵坐标存储在数组A、B中

[A,B]=sub_koch1(a1,b1,a2,b2);

for?i=1:n

???forj=1:length(A)/5;

???????w=sub_koch2(A(1+5*(j-1):5*j),B(1+5*(j-1):5*j));

???????fork=1:4?