利用matlab绘制简单IFS图形（Sierpinski三角形和BarnsleyFern巴恩斯利蕨）

hyhhyh21

已于 2023-08-23 10:58:14 修改

阅读量8.4k

点赞数 10

文章标签： matlab IFS图形分形谢尔宾斯基三角巴恩斯利蕨

于 2019-02-17 21:04:41 首次发布

本文链接：https://blog.csdn.net/weixin_42943114/article/details/87552879

版权

利用matlab绘制简单IFS图形（Sierpinski谢尔宾斯基三角形和BarnsleyFern巴恩斯利蕨）

一、SierpinskiTriangle谢尔宾斯基三角形

谢尔宾斯基三角形（英语：Sierpinski triangle）是一种分形，由波兰数学家谢尔宾斯基在1915年提出。它的维度为log(3)/log(2)=1.585。（来自百度）

下面用迭代的方式绘制谢尔宾斯基三角形，方法如下：

在平面上任取三个不共线的点A、B、C作为三角形的3个顶点，之后选取另一点作为当前点。
1计算当前点与随机某一顶点的中点O
2绘制该中点O，并将该点作为新的当前点。重复步骤1
在这里插入图片描述
效果如下：

虽然初始当前点在三角形之外，但是随着迭代，很快便具有三角形的外观。随着时间的增加呈现出谢尔宾斯基三角形。

谢尔宾斯基三角形具有空间的自相似性。比如根据上面的算法得到的正三角形如下图所示（共计2^17个点）。
在这里插入图片描述
对上图进行灰度化显示，按分辨率大小划分区间，将点数多的区间标记为白色，将点数少的区间标记为黑色。逐渐提高图像分辨率，可以得到如下动图：

matlab代码如下：

%谢尔宾斯三角形
%给定初始3点
Tripots=rand(3,2);
Tripots=[0,1;1*sin(2*pi/3),1*cos(2*pi/3);1*sin(4*pi/3),1*cos(4*pi/3)];

%定义SierpinskiTriangle三角形点数
N=2^17;%定义1000个点
SierpinskiTriangle=zeros(N+3,2);
SierpinskiTriangle(1:3,:)=Tripots;
%定义参考点
ReferencePoint=rand(1,2);

for j=1:N
    %选取参考点与随机角点的中点，作为新的参考点
    ReferencePoint=(Tripots(randi([1,3],1),:)+ReferencePoint)/2;
    %记录在SierpinskiTriangle中
    SierpinskiTriangle(3+j,:)=ReferencePoint;
end

scatter(SierpinskiTriangle(:,1),SierpinskiTriangle(:,2),0.1,'filled')
axis off