[分形几何]IFS系统的实现.

最新推荐文章于 2023-03-28 17:09:08 发布
最新推荐文章于 2023-03-28 17:09:08 发布
阅读量2.9k 收藏 2
点赞数
文章标签： string function encoding random 下载工具 iterator
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/EmilMatthew/article/details/530264
版权

[分形几何]IFS系统的实现.

                                                                                                                   EmilMatthew

                                                                                                                          05/11/15

       IFS(Iterator Function System)是分形几何中一个十分生动而且吸引人的领域.

       用这个系统，可以生成很多自然景观，如树等.

       先看用这个系统可以生成的一些效果:

       http://www.flash8.net/UploadFile/2005-11/20051115211728611.swf

       很有意思，不是吗?

       那么，这样的图形是如何通过程序进行控制的呢，其实从应用的角度去理解，还是相当好懂的。

       那就是仿射坐标变换.

       何谓仿射坐标变换,便是旋转，扭曲，平移三种效果的迭加。

       数学上对应的变换矩阵为:

       http://www.flash8.net/UploadFile/2005-11/20051115211745488.jpg

       所以，只要能根据我们最后所需要的迭代效果，确定出a,b,c,d,e,f的具体取值（当然，这同时也是最难的），根据我们所定出的需要显示的像素点，便能达到显示的效果.

       下面摘录一些我所收集的IFS系统的相应参数(xml文件数据格式)

       <data>中的参数依次为a,b,c,d,e,f,p(p指的是这种迭代效果出现的概率)

 

LEVY曲线:

<inputDatas iterLimit='10000'  fps='20' offsetX='100' offsetY='250' conditionTimes='2' scale='100'>                                                         

       <data>0.5,-0.5,0.5,0.5,0,0,0.5</data>

       <data>0.5,0.5,-0.5,0.5,0.5,0.5,1</data>

</inputDatas>

       http://www.flash8.net/UploadFile/2005-11/20051115211813246.bmp

 

分形树:

<inputDatas iterLimit='10000'  fps='20' offsetX='100' offsetY='250' conditionTimes='4' scale='300'>                                                         

       <data>0,0,0,0.5,0,0,0.05</data>

       <data>0.42,-0.42,0.42,0.42,0,0.2,0.45</data>

       <data>0.42,0.42,-0.42,0.42,0,0.2,0.85</data>

       <data>0.1,0,0,0.1,0,0.2,1</data>

</inputDatas>

      

羊齿草:

<?xml version="1.0" encoding="gb2312" ?>

 <inputDatas iterLimit="20000" fps="20" offsetX="150" offsetY="500" conditionTimes="4" scale="50">

  <data>0,0,0,0.16,0,0,0.01</data>

  <data>0.85,0.04,-0.04,0.85,0,1.6,0.86</data>

  <data>0.2,-0.26,0.23,0.22,0,1.6,0.93</data>

  <data>-0.15,0.28,0.26,0.24,0,0.44,1</data>

  </inputDatas>

 

FLAMBOYENT皇冠:

<?xml version="1.0" encoding="gb2312"?>

<inputDatas iterLimit='20000'  fps='20' offsetX='150' offsetY='500' conditionTimes='5' scale='200'>                                                         

       <data>0.25,0,0,0.5,0,0,0.154</data>

       <data>0.5,0,0,0.5,-0.25,0.5,0.461</data>

       <data>-0.25,0,0,-0.25,0.25,1,0.539</data>

       <data>0.5,0,0,0.5,0,0.75,0.846</data>

       <data>0.5,0,0,-0.25,0.5,1.25,1</data>

</inputDatas>

       http://www.flash8.net/UploadFile/2005-11/20051115211843705.bmp

 

下面给出的是AS2中的具体实现(只列出核心部分):

function mainF(inTransXml:XML)

{

       var i:Number=0;

       var j:Number=0;

       var tmpStr:String=new String();

       var tmpArr:Array=new Array();

      

       _root.gIteratorLimit=parseInt(inTransXml.firstChild.attributes.iterLimit);

       _root.gFps=parseInt(inTransXml.firstChild.attributes.fps);

 

       _root.gXOffset=parseInt(inTransXml.firstChild.attributes.offsetX);

       _root.gYOffset=parseInt(inTransXml.firstChild.attributes.offsetY);

      

       _root.gConditionTimes=parseInt(inTransXml.firstChild.attributes.conditionTimes);

      

       _root.gScale=parseInt(inTransXml.firstChild.attributes.scale);

      

       for(i=0;i<_root.gConditionTimes;i++)

       {

              tmpStr=new String(inTransXml.firstChild.childNodes[i].childNodes[0]);

              tmpArr=tmpStr.split(",");

      

              _root.a[i]=Number(tmpArr[0]);

              _root.b[i]=Number(tmpArr[1]);

              _root.c[i]=Number(tmpArr[2]);

              _root.d[i]=Number(tmpArr[3]);

              _root.e[i]=Number(tmpArr[4]);

              _root.f[i]=Number(tmpArr[5]);

              _root.gCondition[i]=Number(tmpArr[6]);

       }

             

       _root.Xn=_root.Xn_1=_root.Yn=_root.Yn_1=0;

       _root.gIterTimes=0;

       _root.gIntervalID=setInterval(IFSFun,_root.gFps);

}

 

function  IFSFun():Void

{

       var curIndex:Number=0;

       var stepLen:Number=100;

       var i:Number=0;

      

       i=0;

       //trace(String(_root.a[0])+String(_root.b[0])+String(_root.c[0])+String(_root.d[0])+String(_root.e[0])+String(_root.f[0])+String(_root.gCondition[0]));

       //trace(String(a[1])+String(b[1])+String(c[1])+String(d[1])+String(e[1])+String(f[1])+String(_root.gCondition[1]));

       while(i<stepLen&&_root.gIterTimes<_root.gIteratorLimit)

       {

              curIndex=Math.random();

              //trace(_root.gIterTimes);

              //trace(curIndex);

              switch(_root.gConditionTimes)

              {

              case 2:

                     if(curIndex>=0&&curIndex<=_root.gCondition[0])

                     {

                            _root.Xn_1=_root.a[0]*_root.Xn+_root.b[0]*_root.Yn+_root.e[0];

                            _root.Yn_1=_root.c[0]*_root.Xn+_root.d[0]*_root.Yn+_root.f[0];

                            //trace("===level1==");

                           

                            //trace("Xn_1"+String(_root.Xn_1));

                            //trace("Yn_1"+String(_root.Yn_1));

                     }

                     else

                     {

                            _root.Xn_1=_root.a[1]*_root.Xn+_root.b[1]*_root.Yn+_root.e[1];

                            _root.Yn_1=_root.c[1]*_root.Xn+_root.d[1]*_root.Yn+_root.f[1];

                            //trace("===level2==");

                           

                            //trace("Xn_1"+String(_root.Xn_1));

                            //trace("Yn_1"+String(_root.Yn_1));

                     }

                     break;

      

              case 3:

              if(curIndex>=0&&curIndex<=_root.gCondition[0])

              {

                            _root.Xn_1=_root.a[0]*_root.Xn+_root.b[0]*_root.Yn+_root.e[0];

                            _root.Yn_1=_root.c[0]*_root.Xn+_root.d[0]*_root.Yn+_root.f[0];

              }

              else if(curIndex<=_root.gCondition[1])

              {

                            _root.Xn_1=_root.a[1]*_root.Xn+_root.b[1]*_root.Yn+_root.e[1];

                            _root.Yn_1=_root.c[1]*_root.Xn+_root.d[1]*_root.Yn+_root.f[1];

              }

              else

              {

                            _root.Xn_1=_root.a[2]*_root.Xn+_root.b[2]*_root.Yn+_root.e[2];

                            _root.Yn_1=_root.c[2]*_root.Xn+_root.d[2]*_root.Yn+_root.f[2];

              }

              break;

      

              case 4:

              if(curIndex>=0&&curIndex<=_root.gCondition[0])

              {

                            _root.Xn_1=_root.a[0]*_root.Xn+_root.b[0]*_root.Yn+_root.e[0];

                            _root.Yn_1=_root.c[0]*_root.Xn+_root.d[0]*_root.Yn+_root.f[0];

              }

              else if(curIndex<=_root.gCondition[1])

              {

                            _root.Xn_1=_root.a[1]*_root.Xn+_root.b[1]*_root.Yn+_root.e[1];

                            _root.Yn_1=_root.c[1]*_root.Xn+_root.d[1]*_root.Yn+_root.f[1];

              }

              else if(curIndex<=_root.gCondition[2])

              {

                            _root.Xn_1=_root.a[2]*_root.Xn+_root.b[2]*_root.Yn+_root.e[2];

                            _root.Yn_1=_root.c[2]*_root.Xn+_root.d[2]*_root.Yn+_root.f[2];

              }

              else

              {    

                            _root.Xn_1=_root.a[3]*_root.Xn+_root.b[3]*_root.Yn+_root.e[3];

                            _root.Yn_1=_root.c[3]*_root.Xn+_root.d[3]*_root.Yn+_root.f[3];

              }

              break;

      

              case 5:

              if(curIndex>=0&&curIndex<=_root.gCondition[0])

              {

                            _root.Xn_1=_root.a[0]*_root.Xn+_root.b[0]*_root.Yn+_root.e[0];

                            _root.Yn_1=_root.c[0]*_root.Xn+_root.d[0]*_root.Yn+_root.f[0];

              }

              else if(curIndex<=_root.gCondition[1])

              {

                            _root.Xn_1=_root.a[1]*_root.Xn+_root.b[1]*_root.Yn+_root.e[1];

                            _root.Yn_1=_root.c[1]*_root.Xn+_root.d[1]*_root.Yn+_root.f[1];

              }

              else if(curIndex<=_root.gCondition[2])

              {

                            _root.Xn_1=_root.a[2]*_root.Xn+_root.b[2]*_root.Yn+_root.e[2];

                            _root.Yn_1=_root.c[2]*_root.Xn+_root.d[2]*_root.Yn+_root.f[2];

              }

              else if(curIndex<=_root.gCondition[3])

              {    

                            _root.Xn_1=_root.a[3]*_root.Xn+_root.b[3]*_root.Yn+_root.e[3];

                            _root.Yn_1=_root.c[3]*_root.Xn+_root.d[3]*_root.Yn+_root.f[3];

              }

              else

              {    

                            _root.Xn_1=_root.a[4]*_root.Xn+_root.b[4]*_root.Yn+_root.e[4];

                            _root.Yn_1=_root.c[4]*_root.Xn+_root.d[4]*_root.Yn+_root.f[4];

              }

             

                     break;

 

              default:

                     trace3("error while random num produce/n");

                     trace("error while random num produce/n");

                     break;

              }

              //switch ends

              //draw new node

              drawNode2(_root.gXOffset+Xn_1*_root.gScale,_root.gYOffset-Yn_1*_root.gScale,_root.gColorArr[13]);

              //update data.

              _root.Xn=_root.Xn_1;

              _root.Yn=_root.Yn_1;

              i++;

              _root.gIterTimes++;

       }    

      

       if(_root.gIterTimes==_root.gIteratorLimit)

              clearInterval(_root.gIntervalID);

}

 

function drawNode2(x:Number,y:Number,nodeColor:Number):Void

{

       //trace("invoke");

       with(eval("_root.gBrush"))

       {

              lineStyle(1,nodeColor,100);

              moveTo(x-0.5,y);

              lineTo(x+0.5,y);

       }

}

 

源码下载:

http://free3.e-168.cn/as2forward/downloads/IFSSystem.rar

//如果上面这个链接无法响应下载（有可能是被网站给屏蔽掉了），则可使用下载工具（如迅雷等）下载。欢迎提出批评与指正意见.

EmilMatthew
关注
  • 0
    点赞
  • 踩
  • 2
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
用vc++实现IFS分形算法画一棵树
chaos的专栏
02-23 2770
void IFS() { double a[][8]={ 0,0.195,-0.488 ,0.344 ,0.433 ,0.4431 ,0.2452 ,0.25 ,\ 0,0.462 ,0.414 ,-0.252 ,0.361 ,0.2511 ,0.5692 ,0.25,\ 0,-0.058 ,-0.07 ,0.453 ,-0.111 ,0.5976 ,0.0969 ,0.25 ,\
IFS分形
weixin_45028999的博客
01-09 735
分形
经典的分形算法
栖客
11-13 1万+
转载自：http://www.douban.com/note/230496472/ 被誉为大自然的几何学的分形（Fractal）理论，是现代数学的一个新分支，但其本质却是一种新的世界观和方法论。它与动力系统的混沌理论交叉结合，相辅相成。它承认世界的局部可能在一定条件下，在某一方面（形态，结构，信息，功能，时间，能量等）表现出与整体的相似性，它承认空间维数的变化既可以是离散的也可以是连续的，因而拓
微信运动排行榜.html
六桥风月IT随笔
06-17 4万+
<!DOCTYPE html> <html> <head> <meta charset="UTF-8"> <meta id="viewport" name="viewport" content="width=device-width, initial-scale=1.0, maximum-scale=1.0, user-scalable=0" /> <meta content
91.vido.ws v.php,"token" parameter not in video info for unknown reason;
weixin_35307279的博客
03-29 77万+
youtube-dl -v --dump-pages https://www.youtube.com/watch?v=USg3NR76XpQ output[debug] System config: [][debug] User config: [u'--add-metadata', u'--user-agent', u'Mozilla/5.0 (X11; Ubuntu; Linux x86_64...
http://syy7.com/a/25.php,vue_element/vue_elementui.sql at 3ffdcb477ce3fbbc24cf1a4d7e3481a187023a02 ·...
热门推荐
weixin_34064233的博客
04-13 85万+
INSERT INTO `user_tb` VALUES ('1', '女', 'burtyang', 'C041580E1C6192B6084E03CE76D5C05C', '18', null, 'data:image/jpeg;base64,/9j/4AAQSkZJRgABAQEASABIAAD/2wBDAAwICQoJBwwKCQoNDAwOER0TERAQESMZGxUdKiUsKykl...
MATLAB IFS 分形算法
weixin_33782386的博客
03-12 505
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); fo...
分形造型的常用模型，C构造简单的IFS图形
qq_43606914的博客
06-23 3581
分形造型的常用模型随机插值模型迭代函数系统（IFS）L系统 随机插值模型 该模型是1982年由Alain Fournier, DonFussell和Loren Carpenter提出的，它能有效地模拟海岸线和山等自然景象。 为了克服传统模型技术中依赖观察距离的局限性，随机插值模型不是事先决定各种图素和尺度，而是用一个随机过程的采样路径作为构造模型的手段。例如，构造二维海岸线的模型可以 选择控制大致形状的若干初始点； 在相邻两点构成的线段上取其中点，并沿垂直连线方向随机偏移一个距离； 再将偏移后的点与该线段
matlab的ifs程序,[分形几何]IFS系统的实现.
weixin_29338423的博客
03-24 370
[分形几何]IFS系统的实现.EmilMatthew05/11/15IFS(Iterator Function System)是分形几何中一个十分生动而且吸引人的领域.用这个系统，可以生成很多自然景观，如树等.先看用这个系统可以生成的一些效果:http://www.flash8.net/UploadFile/2005-11/20051115211728611.swf很有意思，不是吗?那么，这样的图...
IFS 分形树 分形所有图形
12-27
〖IFS演示程序 V1.4〗 ----------------------------------------------- 作者：赵翔鹏 ([email protected]) 指导：李水根教授 ([email protected]) ----------------------------------------------- ◎ 简介 本软件基于IFS理论，可以产生各种二维分形图形。 可以从http://www.is.pku.edu.cn/~zxp/下载最新版本。 ◎ 更新 与V1.3相比，增加了“全屏”功能，可以绘制大幅的图形。 同时重写了部分代码，使之更精练。 ◎ 操作说明 注意：如果在调色板中选择白色， 则会使用“彩色”方式绘制图形。（等价于单击“颜色设置”－“彩色”。） 如果在调色板中选择其他颜色， 则会使用“单色”方式绘制。 其他功能都是显而易见的，请自己尝试。
IFS分形艺术
04-18
NULL 博文链接：https://huangfeihome.iteye.com/blog/1847204
ifs.zip_ifs_timeiuv_分形_分形几何_迭代
07-14
用迭代函数ifs画分形几何，非常有用，可以仿真画自然界的几何事物
基于IFS分形树的模拟
06-27
根据IFS模型构建分形图形的方法和原理,通过观察树木等自然景物的特征,抽象出一种自然界树木的形状,利用拼贴的方法计算出该树木的IFS码,并用VC++作为工具实现对树木的绘制。通过树木绘制实例详细介绍了绘制的过程、...
IFS.rar_ifs_分形 地形_地形 分形
09-21
一个分形几何中利用IFS系统编写的基于地形的生成
shuye_wcj.rar_分形几何_C/C++_
08-12
TC编写的利用iFS迭代系统产生 一片蕨类树叶
PocketIFS:创建和可视化 IFS l 系统分形
06-30
轻松创建 IFS / l 系统分形（树木、贝壳、海岸线、螺旋线）。 包括 PocketC 源代码。 需要最新版本的 PocketC 运行时（免费）或编译器（共享软件），可从 “PocketIFS 是一款适用于业余应用数学家的不错的软件。它...
vue pc端打印二维码
jingruoannan的博客
07-12 37万+
import Vue from "vue"; let a = "data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAeAAAAKoCAYAAACvNywJAAAABHNCSVQICAgIfAhkiAAAIABJREFUeF7svQmcZVdV73+6q7tr6Bq6q+fupDN15lETZEgkiIkSME4BgYeAQp78ow/Rv/8HweERnnwk8p5/RP+apwQUkI9IiCI8jRoQ0BDgETAkIWNn6KTn6qGqurrmqv6v7.
【个人】网页同步
波波诺夫的专栏
05-27 8万+
<!-- This is an automatically generated file. It will be read and overwritten. DO NOT EDIT! --> Bookmarks Bookmarks 书签栏 BBS 虎踞龙蟠 百度贴吧
基于Jupyterlab的 csv 3D 2D 模型转换
qq_53779256的博客
03-28 2万+
【代码】基于Jupyterlab的 csv 3D 2D 模型转换。
matlab分形几何书籍
最新发布
08-27
关于Matlab分形几何方面的书籍，我推荐以下几本...这些书籍对于想了解和学习Matlab分形几何的人来说都是不错的选择，通过阅读和实践，可以帮助读者理解分形几何的基本概念和算法，并掌握在Matlab中实现分形图形的方法。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值