%函数功能:
%模拟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
matlab画分形图
5星 · 超过95%的资源 需积分: 44 168 浏览量
2012-04-05
16:37:45
上传
评论 2
收藏 2KB RAR 举报
ronghao1230
- 粉丝: 0
- 资源: 2