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

我的博客

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

 
 
 

日志

 
 

分形艺术画(5)——菜花黄  

2016-03-15 15:41:42|  分类: 分形艺术画 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
春风和煦菜花
分形艺术画(5) - liu-xiaoyao - 我的博客
 
分形艺术画(5) - liu-xiaoyao - 我的博客
 
分形艺术画(5) - liu-xiaoyao - 我的博客
 
分形艺术画(5) - liu-xiaoyao - 我的博客

分形艺术画(5) - liu-xiaoyao - 我的博客
 
分形艺术画(5) - liu-xiaoyao - 我的博客
 
分形艺术画(5) - liu-xiaoyao - 我的博客
 
分形艺术画(5) - liu-xiaoyao - 我的博客
 
分形艺术画(5) - liu-xiaoyao - 我的博客
 
%L系统植物----春风和煦菜花黄
clear;clc;
%%%(一)
XQ=[1 2 3 2 3 1 4 2 3 1 2 4 ];XQ=[XQ XQ*1.2 XQ*1.1];
YQ=[1 1 1 1 2 3 1 2 1 2 2 3 ];YQ=[YQ YQ*0.8 YQ*0.5];
ZQ=-[1 1.1 1.4 1.2 1.3 1.2 1.1 1 1 1 1 1.2];ZQ=[ZQ 1.3*ZQ 1.2*ZQ];
%mesh(XQ,YQ,ZQ);
SI=size(XQ);
for Q=1:SI(2)   
    X1=XQ(Q);
    Y1=YQ(Q);
    Z1=[0 ZQ(Q)];
    n=1.5;
    JJ=20;
    TT=1;  
for J=1:JJ                     %J画1层
    X3=zeros(1,2^J);Y3=X3;Z3=X3;
        if J==1; 
            plot3([X1(1) X1(1)],[Y1(1) Y1(1)],[Z1(2) Z1(1)],....
                    'g-','linewidth',4);hold on;
        end
        lk=1;if J>15;lk=J;end
        X2(1)=X1(1)+cos((J-1)*2.2/3*pi)/(lk*3);
        X2(2)=X1(1);
        %            
        Y2(1)=Y1(1)+sin((J-1)*2.3/3*pi)/(lk*3);
        Y2(2)=Y1(1);
        %
        Z2(1)=Z1(1)+1/J;
        Z2(2)=Z1(1)+1/J;
        %
       
          for L=1:2
            if L==1
                plot3([X1(1) X2(L)],[Y1(1) Y2(L)],[Z1(1) 1.1*Z2(L)],....
                    'g-','linewidth',2);hold on;
                plot3([X2(L)],[Y2(L)],1.1*[Z2(L)],....
                    'y+','linewidth',5);hold on;
            else
                plot3([X1(1) X2(L)],[Y1(1) Y2(L)],[Z1(1) Z2(L)],....
                    'g-','linewidth',4);hold on;
                
            end
             if J==JJ
                plot3([X1(1)],[Y1(1)],[Z1(1)],....
                    'y*','linewidth',10);hold on;
             end
          end       
        %
        X3(1:2)=X2(1:2);
        Y3(1:2)=Y2(1:2);
        Z3(1:2)=Z2(1:2);
        Z1=Z3;
    end%J
end%Q
grid;
for J=1:5;
    X(:,J)=[0 1 2 3 4 5 6];
    Z(:,J)=-[1.0 0.9 0.8 0.7 0.7 0.9 1.0]*(1+J/10);
    Z(:,1)=-1.8*[0.9 0.8 0.7 0.8 0.9 0.9 1];
end
for I=1:7
Y(I,:)=[0;1;2;3;4];
end
surf(X,Y,Z);
  评论这张
 
阅读(42)| 评论(4)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

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

页脚

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