frac_tree.m

本人收集的一些有关matlab的代码程序设计 也不知道改选什么分类

function [x,y]=frac_tree(x0,y0,v,N)
%FRAC_TREE 函数用来演示函数的编写方法。它以分形树为背景加以
%   介绍，详细内容见书中第 3 章。
%
%   [x,y]=frac_tree(x0,y0,v,N)
%   给定初始点 x0, y0, 随机向量 v 和随机点个数 N, 则可以
%   计算出分形树上的 N 个点坐标向量 x 和 y。

%Designed by Prof D Xue (c) 2000
x=[x0; zeros(N-1,1)]; y=[y0; zeros(N-1,1)];
for i=2:N
   vv=v(i);
   if vv<0.05, y(i)=0.5*y(i-1);   
   elseif vv<0.45, 
      x(i)=0.42*(x(i-1)-y(i-1)); y(i)=0.2+0.42*(x(i-1)+y(i-1));
   elseif vv<0.85, 
      x(i)=0.42*(x(i-1)+y(i-1)); y(i)=0.2-0.42*(x(i-1)-y(i-1));   
   else, 
      x(i)=0.1*x(i-1); y(i)=0.1*y(i-1)+0.2;
   end
end