分形与小波应用matlab画分形树.doc

第 1 页 共 2 页一、实验内容一、实验内容应用 matlab 画分形树二、使用环境二、使用环境Matlab四、调试过程四、调试过程程序为：function tree(n,a,b) n=8;a=pi/8;b=pi/8; x1=0;y1=0; x2=0;y2=1; plot([x1,x2],[y1,y2]) hold on [X,Y]=tree1(x1,y1,x2,y2,a,b); hold on W=tree2(X,Y); w1=W(:,1:4); w2=W(:,5:8); w=[w1;w2]; for k=1:nfor i=1:2^k[X,Y]=tree1(w(i,1),w(i,2),w(i,3),w(i,4),a,b);W(i,:)=tree2(X,Y); end w1=W(:,1:4);w2=W(:,5:8);w=[w1;w2]; end function [X,Y]=tree1(x1,y1,x2,y2,a,b) L=sqrt((x2-x1)^2+(y2-y1)^2); if (x2-x1)==0a=pi/2;课程课程:分型与小波实验日期实验日期:2012/10/20成员人数成员人数:3专业：专业： 姓名：姓名：学号：学号：专业：专业： 姓名：姓名： 学号：学号：专业：专业： 姓名：姓名： 学号：学号：第 2 页 共 2 页else if (x2-x1)<0a=pi+atan((y2-y1)/(x2-x1));elsea=atan((y2-y1)/(x2-x1)); end end x3=x2+L*2/3*cos(a+b); y3=y2+L*2/3*sin(a+b); x4=x2+L*2/3*cos(a-b); y4=y2+L*2/3*sin(a-b); a=[x3,x2,x4]; b=[y3,y2,y4]; plot(a,b) axis equal hold on X=[x2,x3,x4]; Y=[y2,y3,y4]; function w=tree2(X,Y) a1=X(1);b1=Y(1); a2=X(2);b2=Y(2); a3=X(1);b3=Y(1); a4=X(3);b4=Y(3); w=[a1,b1,a2,b2,a3,b3,a4,b4];输出结果为: 。