蕨型叶分形

目录

要点

基本语句

EraseMode

习题

[1 设置颜色](#1 设置颜色)

[2 旋转蕨型叶图](#2 旋转蕨型叶图)

[3 枝干](#3 枝干)

[4 塞平斯基三角形](#4 塞平斯基三角形)


要点


蕨型叶是通过一个点的反复变换产生的,假设x是一个含有两个分量的向量,可以用来表示平面内的一个点,则可以用Ax+b的形式对其进行变换。


基本语句


语句darkgreen=[0 2/3 0]设置深绿色的颜色变量,在matlab中用红绿蓝分量来表示颜色,这里将绿色分量设置为2/3,其余为0。
采用animatedline来绘制蕨型叶,h为蕨型叶图形的句柄。

复制代码
function fern_a


    shg
    clf reset
    set(gcf,'color','white','menubar','none', ...
        'numbertitle','off','name','Fractal Fern')
    x = [.5; .5];
    h = animatedline('Marker', '.', 'LineStyle', 'none', 'Color', [0 2/3 0]);
    axis([-3 3 0 10])
    axis off  %不显示坐标轴,只显示叶子本身
    stop = uicontrol('style','toggle','string','stop', ...
        'background','white');
    drawnow %绘制

    p  = [ .85  .92  .99  1.00];
    A1 = [ .85  .04; -.04  .85];  b1 = [0; 1.6];
    A2 = [ .20 -.26;  .23  .22];  b2 = [0; 1.6];
    A3 = [-.15  .28;  .26  .24];  b3 = [0; .44];
    A4 = [  0    0 ;   0   .16];

    cnt = 1; %将计数器的初始值设置为1
    tic
    while ~get(stop,'value')
        r = rand;
        if r < p(1)
            x = A1*x + b1;
        elseif r < p(2)
            x = A2*x + b2;
        elseif r < p(3)
            x = A3*x + b3;
        else
            x = A4*x;
        end
        addpoints(h, x(1), x(2));
        drawnow %重新绘制图形
        cnt = cnt + 1; %生成新的点后计数器的数值自增1
    end
    t = toc; %读取秒表数据
    s = sprintf('%8.0f points in %6.3f seconds',cnt,t);
    text(-1.5,-0.5,s,'fontweight','bold');
    set(stop,'style','pushbutton','string','close','callback','close(gcf)')
end

EraseMode


EraseMode 属性不再受支持,而且在以后的版本中会出错。

初始语句:

复制代码
h = plot(x(1),x(2),'.');
darkgreen = [0 2/3 0];
set(h,'markersize',1,'color',darkgreen,'erasemode','none');

修改后:

复制代码
h = animatedline('Marker', '.', 'LineStyle', 'none', 'Color', [0 2/3 0]);

习题


1 设置颜色

在MATLAB中,颜色可以使用多种方式来表示,其中一种常用的方式是使用RGB(Red-Green-Blue)表示法。在RGB表示法中,每个颜色的红、绿、蓝三个分量的取值范围为0到1,表示颜色的深浅程度。

复制代码
function fern_a


    shg
    clf reset
    set(gcf,'color','black','menubar','none', ...
        'numbertitle','off','name','Fractal Fern')
    x = [.5; .5];
    h = animatedline('Marker', '.', 'LineStyle', 'none', 'Color', [1 0.75 0.8]);
    axis([-3 3 0 10])
    axis off
    stop = uicontrol('style','toggle','string','stop', ...
        'background','black','ForegroundColor','white');
    drawnow

    p  = [ .85  .92  .99  1.00];
    A1 = [ .85  .04; -.04  .85];  b1 = [0; 1.6];
    A2 = [ .20 -.26;  .23  .22];  b2 = [0; 1.6];
    A3 = [-.15  .28;  .26  .24];  b3 = [0; .44];
    A4 = [  0    0 ;   0   .16];

    cnt = 1;
    tic
    while ~get(stop,'value')
        r = rand;
        if r < p(1)
            x = A1*x + b1;
        elseif r < p(2)
            x = A2*x + b2;
        elseif r < p(3)
            x = A3*x + b3;
        else
            x = A4*x;
        end
        addpoints(h, x(1), x(2));
        drawnow
        cnt = cnt + 1;
    end
    t = toc;
    s = sprintf('%8.0f points in %6.3f seconds',cnt,t);
    text(-1.5,-0.5,s,'fontweight','bold','Color','red');
    set(stop,'style','pushbutton','string','close','callback','close(gcf)')
end

2 旋转蕨型叶图

复制代码
%修改坐标
addpoints(h, x(2), x(1));

%修改标注位置
s = sprintf('%8.0f points in %6.3f seconds',cnt,t);
text(2.5,-2.9,s,'fontweight','bold','Color','red');

3 枝干

复制代码
A4 = [0, 0; 0, 0.16];
%修改为
A4 = [0, 0; 0, 0.08];
%这将使得蕨型叶的枝干更加瘦长

蕨型叶的起始点均为(0,0)。


4 塞平斯基三角形

复制代码
function sierpinski_triangle_a
    shg
    clf reset
    set(gcf, 'color', 'white', 'menubar', 'none', ...
        'numbertitle', 'off', 'name', 'Fractal Sierpinski Triangle')
    x = [.5; sqrt(3)/6];
    h = animatedline('Marker', '.', 'LineStyle', 'none', 'Color', [0 2/3 0]);
    axis([0 1 0 sqrt(3)/2])
    axis off
    stop = uicontrol('style', 'toggle', 'string', 'stop', ...
        'background', 'white');
    hold on
    plot([0, 1, 0.5, 0], [0, 0, sqrt(3)/2, 0], 'k-*'); % 绘制三角形的顶点
    drawnow

    A1 = [0.5, 0; 0, 0.5];  b1 = [0; 0];
    A2 = [0.5, 0; 0, 0.5];  b2 = [0.5; 0];
    A3 = [0.5, 0; 0, 0.5];  b3 = [0.25; sqrt(3)/4];

    cnt = 1;
    tic
    while ~get(stop, 'value')
        r = rand;
        if r <= 1/3
            x = A1 * x + b1;
        elseif r <= 2/3 & r>1/3
            x = A2 * x + b2;
        else
            x = A3 * x + b3;
        end
        addpoints(h, x(1), x(2));
        drawnow
        cnt = cnt + 1;
    end
    t = toc;
    s = sprintf('%8.0f points in %6.3f seconds', cnt, t);
    text(0.25, -0.05, s, 'fontweight', 'bold');
    set(stop, 'style', 'pushbutton', 'string', 'close', 'callback', 'close(gcf)')
end
相关推荐
松果集33 分钟前
MATLAB基础知识【8】绘图,peaks,sphere,积分
matlab
freexyn2 天前
Matlab自学笔记六十一:快速上手解方程
数据结构·笔记·matlab
ytttr8732 天前
matlab通过Q学习算法解决房间路径规划问题
学习·算法·matlab
梦子要转行2 天前
matlab/Simulink-全套50个汽车性能建模与仿真源码模型9
开发语言·matlab·汽车
Zevalin爱灰灰2 天前
MATLAB GUI界面设计 第六章——常用库中的其它组件
开发语言·ui·matlab
曹勖之12 天前
simuilink和ROS2数据联通,Run后一直卡在Initializting
windows·matlab·simulink·ros2
Zevalin爱灰灰12 天前
MATLAB GUI界面设计 第三章——仪器组件
开发语言·ui·matlab
算法如诗12 天前
基于SOA(海鸥优化算法)的路径规划Matlab实现方案
开发语言·算法·matlab
项目申报小狂人12 天前
2025年中科院三区全新算法,恒星振荡优化器:受自然启发的元启发式优化,完整MATLAB代码免费获取
开发语言·算法·matlab
Zevalin爱灰灰13 天前
MATLAB GUI界面设计 第二章——APP Designer操作正式入门
开发语言·ui·matlab