Matlab绘制分形图

%函数功能： %模拟Koch曲线 %程序设计：JianCai Hu %设计时间：2010-05-17 %参考文献: %[1] 段俊生、葛素侨，用MATLAB模拟分形[J].天津商学院学报，2005.6(11),p60-64 clc clear a = 10; k = 5; angle = pi/3; figure; axis equal; A = zeros(2,2*4^k); B = A; A(:,1:2) = [0 3*a;0 0]; C = [cos(angle) -sin(angle);sin(angle) cos(angle)]; D = inv(C); for n = 1:k B = A/3; A(:,1:2*4^(n-1)) = B(:,1:2*4^(n-1)); A(:,2*4^(n-1)+1:4^n) = C*B(:,1:2*4^(n-1)) + [a;0]*ones(1,2*4^(n-1)); A(:,4^n+1:6*4^(n-1)) = D*B(:,1:2*4^(n-1)) + [3*a/2;sqrt(3)*a/2]*ones(1,2*4^(n-1)); A(:,6*4^(n-1)+1:2*4^n) = B(:,1:2*4^(n-1)) + [2*a;0]*ones(1,2*4^(n-1)); end for i = 1:4^k L = line(A(1,2*i-1:2*i),A(2,2*i-1:2*i)); set(L,'LineWidth',2); end