二维分形作图.ppt

由生成元生成严格自相似分形图形 Koch曲线的生成过程 生成Koch曲线的MATLAB程序 clear;clf; new=[0,1];%定义初始线段的两端 kmax=20; %迭代次数 for k=1:kmax old=new; %保存原有各点的坐标 n=length(old)-1; %计算需要改变图形的线段数目 for m=0:n-1 %计算各新线段两端点的坐标(5点) diff=(old(m+2)-old(m+1))/3; %取新线段长度向量 new(4*m+1)=old(m+1); %新线段第一点坐标 new(4*m+2)=old(m+1)+diff; %新线段第二点坐标 new(4*m+3)= new(4*m+2)+diff*((1+sqrt(3)*i)/2); %第三点 new(4*m+4)=old(m+1)+2*diff; %第四点 end new(4*n+1)=old(n+1); %最后一点 plot(new) ; axis equal; end 运行结果 Cantor集C的生成过程 Levy曲线(C曲线)的生成过程 分形树的MATLAB程序 clear;clf; th=pi/6;%设置角度 u=[0,0;0,1]; %初始线段坐标 rov1=[cos(th),-sin(th);sin(th),cos(th)] ; rov2=rov1; %坐标变换矩阵 kmax=15; %循环次数 for k=1:kmax uuu=[]; %建立一空数组，保存点坐标 n=length(u)/2; %计算需要改变图形的线段数目 for m=0:n-1 %计算各新端点的坐标 p1=(u(2*m+1,:)*2+u(2*m+2,:))/3 ; %a点 p2=(u(2*m+1,:) +u(2*m+2,:)*2)/3; %b点 pp=[(u(2*m+2,1)-u(2*m+1,1)); (u(2*m+2,2)-u(2*m+1,2))]/3; %计算新线段向量 lp=rov1*pp; %第一条树枝偏转向量 lp=p1+lp; %第一条树枝端点C坐标 rp=rov2*pp; %第二条树枝偏转向量 rp=p2+rp; %第二条树枝端点D坐标 uu=[u(2*m+1,:);p1;p1;lp;p1;p2;p2;rp;p2;u(2*m+2,:)]; %按顺序排列新线段两端点的坐标 uuu=[uuu;uu]; end u=[uuu]; plot(u(:,1),u(:,2)) ; axis([-0.5,0.5,0,1] ); end 分形树运行结果 函数迭代产生的分形 z 表示复数，在复平面上定义函数 f(z)。 任意给定初始复数值 ，定义复数序列 对于什么样的初始值 ，复数序列 收敛或有界？ L 系统 生物学家Lindenmayer引进一种字符串迭代重写机制, 称为L系统.一个L系统可表示为一个有序的三元素集合： 复平面上的牛顿迭代法 从复平面上某一点 出发, 利用牛顿迭代法求方程 的根. 其迭代公式为 离 的根越近, 则迭代过程就越少, 离的越远迭代 时间就越长甚至不收敛. 这使得逃逸时间法可行 文法构图算法 字母表：L，R 生成规则：L→R，R→L R 初始字母：R 则有：R→L R→R L R→L R R L R→R L R L R R L R→L R R L R R L R L R R L R→…… 其中： V 是字母表(运动过程集合)， ω是初始字符串(初始元)， P 是生成规则(生成元)。 F 表示向前距离d,并画线 + 表示左转弯角度δ - 表示右转弯角度δ 其中 例如 F 表示向前距离d,并画线 + 表示左转弯角度δ - 表示右转弯角度δ FFF-FF-F-F+F+FF-FF-FFF 字符串 的解释 Koch曲线的L系统文法生成 δ=60o ω：F ( 初始元) P ： F → F ＋ F － － F ＋ F 第0次迭代：F 第1次迭代：F ＋ F － － F ＋ F 第2次迭代：(F ＋ F － － F ＋ F )＋ (F ＋ F － － F ＋ F) － － (F ＋ F － － F ＋ F )＋ (F ＋ F － － F ＋ F) 第3次迭代：F ＋ F － － F ＋ F ＋ F ＋ F － － F ＋ F － － F ＋ F － － F ＋ F ＋ F ＋ F － － F ＋ F ＋ F ＋