注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

我的博客

清谈不觉晓:完全原创、自娱自乐

 
 
 

日志

 
 

分形艺术画(8)——多层树枝干  

2016-04-25 09:47:18|  分类: 分形艺术画 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
6层树
分形艺术画(8)——多层树枝干 - liu-xiaoyao - 我的博客
 5层树分形艺术画(8)——多层树枝干 - liu-xiaoyao - 我的博客
4层变形树
分形艺术画(8)——多层树枝干 - liu-xiaoyao - 我的博客
 
%Fractals1_Snow2
%L系统植物
clear;clc;
CC=0;
%(一)字符解释
%1----按起点角直行一次;
%20\40----(右\左)转角一次;
%3.0X----退回X步;
%2\4----(右\左)转角一次+直行一次;%3----退回一步;
%5.0X----y+X步;
%+-6.0X----x+-X步;
%(二)左右解释
%(1)R=1―-顺时针转角;%(2)R=-1―逆时针转角;
R=1;
%(三)基因的分支解释
Lm=3;am=R*pi/8;%T=6.5
M9=[1 1 1 20 1 1 40 1 3.01 20 1 3.01 3.09 40 1 1 20 1....
    3.01 40 1 3.01 3.09];%BASE%%%%%%
%   1 1 1 [  1 1 {  1  }   [  1  ]    ]   {  1 1 [  1
%   ]    {  1  }    } ;
%(2-2)P:基因分支_梅花
Lw=0;aw=R*pi/8;
W1=[1 20 5.25 -6.002 1 40 40 5.15 6.01 1 40 1 5.26 20 20 1];
sw=size(W1);
%(四)层数:
NP=6;
%Snow
for P=1:NP;P
    if P==1;XYA1=zeros(1,3);end
    XYA1=zeros(1,3);
    if P>1;XYA1=XYAO;end    
    NP=4^(P);
    NNP=4^(P-1);
  for NU=1:NNP
    X(1)=XYA1(P,(NU-1)*3+1);
    Y(1)=XYA1(P,(NU-1)*3+2);
    A(1)=XYA1(P,(NU-1)*3+3);    
    KW=W1(1);
    if KW==20
        A(1)=A(1)+aw;
        XYA1(1,3)=A(1);
    end
    if KW==40
        A(1)=A(1)-aw;
        XYA1(1,3)=A(1);
    end
    if KW>5 && KW<6
        KK=((KW-5)*100);KK=fix(KK+0.001);
        XYA1(1,2)=XYA1(1,2)+KK;  
    end
    if abs(KW)>6 && abs(KW)<7
        Fu=abs(KW)/KW;
        %KK=((abs(KW)-6)*100);KK=Fu*fix(KK+0.001);
        KK=((abs(KW)-6)*100);KK=Fu*KK;
        XYA1(1,1)=XYA1(1,1)+KK;  
    end    
    if KW==1        
        J0=3;%%%%%%%%%%%%%%%%%%{4/22}
        MKJ=[M9];sm=size(MKJ);                
        [XYA]=Fractals2_Y2(MKJ,Lm,am,sm(2),X(1),Y(1),A(1),P);
        sxy=size(XYA);
        for JJ=1:sxy(1)
            H(JJ)=((XYA(JJ,1)-XYA(1,1))^2+(XYA(JJ,2)-XYA(1,2))^2)^0.5;
        end
        [H,N]=(sort(H));
        sHN=size(H);
        XYN(1:4)=N(sHN(2)-0:-1:sHN(2)-3);
        XYA2(1:4,1:3)=XYA(XYN,1:3);
        for IIP=1:4
            for JJP=1:3
                XYAO(P+1,(NU-1)*12+(IIP-1)*3+JJP)=XYA2(IIP,JJP);
            end%JJP
        end%IIP
     end%KW
   end%NU
end%PN
  评论这张
 
阅读(50)| 评论(3)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2017