分形之海岸线

      分形理论是当今世界十分风靡和活跃的新理论、新学科。分形的概念是美籍数学家曼德布罗特(B.B.Mandelbort)首先提出的。1967年他在美国权威的《科学》杂志上发表了题为《英国的海岸线有多长?》的著名论文。海岸线作为曲线,其特征是极不规则、极不光滑的,呈现极其蜿蜒复杂的变化。我们不能从形状和结构上区分这部分海岸与那部分海岸有什么本质的不同,这种几乎同样程度的不规则性和复杂性,说明海岸线在形貌上是自相似的,也就是局部形态和整体形态的相似。在没有建筑物或其他东西作为参照物时,在空中拍摄的100公里长的海岸线与放大了的10公里长海岸线的两张照片,看上去会十分相似。事实上,具有自相似性的形态广泛存在于自然界中,如:连绵的山川、飘浮的云朵、岩石的断裂口、布朗粒子运动的轨迹、树冠、花菜、大脑皮层……曼德布罗特把这些部分与整体以某种方式相似的形体称为分形(fractal)。1975年,他创立了分形几何学(fractalgeometry)。在此基础上,形成了研究分形性质及其应用的科学,称为分形理论(fractaltheory)。
  在这一节中,将展示一个海岸线的分形图形。

核心代码:

static void FractalTriangle(const Vector3& vStart, const Vector3& vEnd, Yreal angle, Vector3* pVertices)
{
    Vector3 vSub = vEnd - vStart;
    Yreal len = D3DXVec3Length(&vSub);

    pVertices[0] = vStart;
    pVertices[5] = vEnd;

    Yreal alfa = atan2f(vSub.y, vSub.x);

    pVertices[1] = vStart + vSub*0.25f;
    pVertices[4] = vStart + vSub*0.75f;

    Yreal a0 = alfa + angle;
    Yreal l = len/8.0f/cosf(angle);

    pVertices[2].x = pVertices[1].x + cosf(a0)*l;
    pVertices[2].y = pVertices[1].y + sinf(a0)*l;
    pVertices[2].z = 0.0f;

    pVertices[3].x = pVertices[4].x - cosf(a0)*l;
    pVertices[3].y = pVertices[4].y - sinf(a0)*l;
    pVertices[3].z = 0.0f;
}

软件截图:

软件下载地址:http://files.cnblogs.com/WhyEngine/Fractal.7z

转载于:https://www.cnblogs.com/WhyEngine/p/4020641.html

用递归的方法画分形图 用递归的方法画分形图   分形几何是数学领域里新兴的课题,如果将图形的每个元素按某种规则进行变形,得到新的图形,以此类推,进行若干次变形后得到的图形就是分形图形。Couch曲线是最典型的分形图形:   将一条线段按照图1进行变换,得到图1,再将图1中的每条线段按图1的折线变换得到图2,以此类推,进行6次变换就得到图6,如果进行无限次变换,就得到的是Couch曲线,Couch曲线的维数不是整数维,更多详情请见分形方面的书籍。   用分形图形能画图许多漂亮的图案而被广泛地应用,下面将几个简单的分形图形的代码与图形大家分享。 //Couch曲线的画法 void Couch(CDC *pDC,int x1,int y1,int x2,int y2,int n) { //pDC是画图的设备上下文的指针 //x1,y1,x2,y2是起始的两点 //其中参数n是递归的层数 int x3,y3,x4,y4,x5,y5; //以下是根据空间几何计算出来的坐标 x3=x1+(x2-x1)/3; y3=y1+(y2-y1)/3; x4=x1+(x2-x1)*2/3; y4=y1+(y2-y1)*2/3; x5=x3+(x4-x3)/2+int(sqrt(3)*(y4-y3)/2); y5=y3-int(sqrt(3)*(x4-x3)/2)+(y4-y3)/2; //递归最后一层,递归的出口 if(n==1) { pDC->MoveTo(x1,y1); pDC->LineTo(x3,y3); pDC->LineTo(x5,y5); pDC->LineTo(x4,y4); pDC->LineTo(x2,y2); } else { //递归画图 Couch(pDC,x1,y1,x3,y3,n-1); Couch(pDC,x3,y3,x5,y5,n-1); Couch(pDC,x5,y5,x4,y4,n-1); Couch(pDC,x4,y4,x2,y2,n-1); } } //斯宾斯基篓垫的画法 void Floor(CDC *pDC,int x1, int y1,int x2,int y2,int x3,int y3,int n) { //pDC是画图的设备上下文的指针 //x1,y1,x2,y2,x3,y3是起始的三角形的三点坐标 //其中参数n是递归的层数 int x11,x22,x33,y11,y22,y33; //以下是根据空间几何计算出来的坐标 x11=(x2+x3)/2; y11=(y2+y3)/2; x22=(x1+x3)/2; y22=(y1+y3)/2; x33=(x1+x2)/2; y33=(y1+y2)/2; pDC->MoveTo(x11,y11); pDC->LineTo(x22,y22); pDC->MoveTo(x11,y11); pDC->LineTo(x33,y33); pDC->MoveTo(x22,y22); pDC->LineTo(x33,y33); //递归最后一层,递归的出口 if(n==1) { pDC->MoveTo(x11,y11); pDC->LineTo(x22,y22); pDC->LineTo(x33,y33); pDC->LineTo(x11,y11); } else { //递归画图 Floor(pDC,x1,y1,x33,y33,x22,y22,n-1); Floor(pDC,x33,y33,x2,y2,x11,y11,n-1); Floor(pDC,x22,y22,x11,y11,x3,y3,n-1); } } //分形矩形的画法 void Rect(CDC *pDC,int x1,int y1,int x2,int y2,int n) { //pDC是画图的设备上下文的指针 //x1,y1,x2,y2是起始矩形坐标 //其中参数n是递归的层数 int x3,y3,x4,y4,x5,y5,x6,y6; //以下是根据空间几何计算出来的坐标 x3=x1+(x2-x1)/3; y3=y1+(y2-y1)/3; x4=x1+(x2-x1)*2/3; y4=y1+(y2-y1)*2/3; x5=x3+(y4-y3); y5=y3-(x4-x3); x6=x4-(y3-y4); y6=y4+(x3-x4); pDC->MoveTo(x1,y1); pDC->LineTo(x3,y3); pDC->MoveTo(x4,y4); pDC->LineTo(x2,y2); //递归最后一层,递归的出口 if(n==1) { pDC->MoveTo(x1,y1); pDC->LineTo(x3,y3); pDC->LineTo(x5,y5); pDC->LineTo(x6,y6); pDC->LineTo(x4,y4); } else { //递归画图 Rect(pDC,x3,y3,x5,y5,n-1); Rect(pDC,x5,y5,x6,y6,n-1); Rect(pDC,x6,y6,x4,y4,n-1); } } //分形树的画法,其中参数n是递归的层数 void Tree(CDC *pDC,int x1,int y1,int x2,int y2,int n) { //pDC是画图的设备上下文的指针 //x1,y1,x2,y2是起始矩形坐标 //其中参数n是递归的层数 int x3,y3,x4,y4,x5,y5; //以下是根据空间几何计算出来的坐标 x3=x1+(x2-x1)/3; y3=y1+(y2-y1)/3; x4=x3+int((x1-x3)*cos(5*pi/6))-int((y1-y3)*sin(5*pi/6)); y4=y3+int((x1-x3)*sin(5*pi/6))+int((y1-y3)*cos(5*pi/6)); x5=x3+int((x1-x3)*cos(5*pi/6))+int((y1-y3)*sin(5*pi/6)); y5=y3-int((x1-x3)*sin(5*pi/6))+int((y1-y3)*cos(5*pi/6)); pDC->MoveTo(x1,y1); pDC->LineTo(x2,y2); //递归最后一层,递归的出口 if(n==1) { pDC->MoveTo(x3,y3); pDC->LineTo(x4,y4); pDC->MoveTo(x3,y3); pDC->LineTo(x5,y5); } else { //递归画图 Tree(pDC,x3,y3,x2,y2,n-1); Tree(pDC,x3,y3,x4,y4,n-1); Tree(pDC,x3,y3,x5,y5,n-1); } }   上述的代码及图形附有Visual C++源代码,并在Windows XP和Visual C++6.0下调试成功。更多的分形图形及图形坐标空间几何的计算方法请与作者联系。 最新评论 [发表评论] [文章投稿] 查看所有评论 推荐给好友 打印 改成sqrt(3.0)就好了,vs的问题,参数多类型重载没自己识别出来 ( maxint 发表于 2008-9-4 21:27:00) skinfeature界面产品特点介绍 SkinFeature完全支持各种常用控件及窗口 。对话框,单文档界面,多文档界面的全部支持。 完全支持PNG、TGA、bmp 图像格式。 支持皮肤(.rss)文件加密。保护美工图像不被别人非法使用,进一步保护自身知识产权。 支持从msstyle wba 等主题文件中自动转化为rss皮肤文件。 SkinFeature完全支持 VC,Delphi,C#,VB.Net,Visual Basic,C++ Bulider, PowerBuilder, Win32 SDK Supported WTL, ATL, 第三方厂商SDK, OUTLOOK 等等界面换肤。 支持各Windows平台 Windows 9X/NT/2000/2003/XP/Vista 完全支持ANSI UNICODE 文字编码格式。 完全支持动态换肤(*.rss),支持不规则窗体,支持美工设计自定义界面生成方案。 完全多线程,自定义窗口类换肤支持。 内含15种专业皮肤(.rss)文件。 内含SkinDesigner 皮肤文件开发工具。 http://www.skinfeature.com ( skinfeature 发表于 2008-6-18 17:30:00) change sqrt(3) to sqrt(3.0) ( Eric_wang 发表于 2007-10-30 9:59:00) x5=x3+(x4-x3)/2+int(sqrt(3)*(y4-y3)/2); y5=y3-int(sqrt(3)*(x4-x3)/2)+(y4-y3)/2; 运行是错误的啊? 错误 1 error C2668: 'sqrt' : ambiguous call to overloaded function d:\wordplay\c++\范例\fractral递归方法画分图形\stdafx.cpp 19 错误 2 fatal error C1903: unable to recover from previous error(s); stopping compilation d:\wordplay\c++\范例\fractral递归方法画分图形\stdafx.cpp 19
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值