MATLAB画二维树
具体的要求跟出来的结果如两张图所示，之前查过了好像跟元胞数组有关，但是我自己没能理解，希望诸位大神帮帮忙。function tree(n,a,b)
% tree(13,pi/3,pi/5),n为分形树迭代次数 %a,b为分枝与竖直方向夹角
%x1,y1,x2,y2为初始线段两端点坐标,nn为迭代次数
x1=0;y1=0;
x2=0;y2=1;
plot(,)
hold on
=tree1(x1,y1,x2,y2,a,b);
hold on
W=tree2(X,Y);
w1=W(:,1:4);
w2=W(:,5:8);
% w为2^k*4维矩阵,存储第k次迭代产生的分枝两端点的坐标,
% w的第i(i=1,2,…,2^k)行数字对应第i个分枝两端点的坐标
w=;
for k=1:n
for i=1:2^k
=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=;
end
%由每个分枝两端点坐标(x1,y1),(x2,y2)产生两新点的坐标(x3,y3),(x4,y4),画两分枝图形,并把
%(x2,y2)连同新点横、纵坐标分别存储在数组X,Y中
function =tree1(x1,y1,x2,y2,a,b)
L=sqrt((x2-x1)^2+(y2-y1)^2);
if (x2-x1)==0
a=pi/2;
else if(x2-x1)<0
a=pi+atan((y2-y1)/(x2-x1));
else
a=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*1/2*cos(a-b);
y4=y2+L*1/2*sin(a-b);
a=;
b=;
plot(a,b)
axis equal
hold on
X=;
Y=;
%把由函数tree1生成的X,Y顺次划分为两组,分别对应两分枝两个端点的坐标,并存储在一维%数组w中
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=;
你试一下这个程序，长度已经能控制了，不过角度控制似乎存在问题，如果你能调试出来，欢迎与我交流！
MATLAB画二维树 真棒！
真棒啊！！洋洋！！
6666666666666666
页:
[1]