MATLAB IFS 分形算法

最新推荐文章于 2024-11-13 09:04:53 发布

原创最新推荐文章于 2024-11-13 09:04:53 发布 · 2.8k 阅读

1 ·

CC 4.0 BY-SA版权

图像处理专栏收录该内容

79 篇文章

订阅专栏

本文介绍了一种使用MATLAB或Octave实现的IFS（迭代函数系统）绘制算法。该算法通过定义一组仿射变换矩阵和对应的概率分布，随机选择变换来生成分形图像。代码中详细展示了如何初始化参数、迭代计算点坐标并最终绘制出分形图案。

function IFS_draw(M,p)
N=30000; 
for k=1:length(p);
    eval(['a',num2str(k),'=reshape(M(',num2str(k),',:),2,3);']);
end
xy=zeros(2,N); 
pp=meshgrid(p); 
pp=tril(pp); 
pp=sum(pp,2); 
for k=1:N-1;
    a=rand-pp; 
d=find(a<=0);
    xy(:,k+1)=eval(['a',num2str(d(1)),'(:,1:2)'])*xy(:,k)+...
        eval(['a',num2str(d(1)),'(:,3)']); 
end
P=complex(xy(1,:),xy(2,:)); plot(P,'k.', 'MarkerFaceColor','k','markersize',2); 
axis equal; 

clc;clear;close all; M=[0.06,0,0,0.6,0,0;...         0.04,0,0,-0.5,0,1;...         0.46,-0.34,0.32,0.38,0,0.6;...         0.48,0.17,-0.15,0.42,0,1;...         0.43,-0.26,0.27,0.48,0,1;...         0.42,0.35,-0.36,0.31,0,0.8];  p=[0.1,0.1,0.1,0.23,0.23,0.24];  IFS_draw(M,p); 
axis image; 
set(gcf,'Color','w');