数学实验 实验十二 迭代 (2) --分形
迭代与分形
实验名称：迭代与分形专业：信息工程班级：09级四班姓名：序号：29，38提交日期：2011年4月29日一、实验目的与要求1．认识Fibonacci数列，体验发现其通项公式的过程；2．了解matlab软件中进行数据显示与数据拟合的方式；3．掌握matlab软件中plot, polyfit等函数的基本用法；4．提高对数据进行分析与处理的能力。
二、问题描述几何学研究的对象是客观世界中物体的形状。
传统欧氏几何学的研究对象，都是规则并且光滑的，比如：直线、曲线、曲面等。
但客观世界中物体的形状，并不完全具有规则光滑等性质，因此只能近似当作欧氏几何的对象，比如：将凹凸不平的地球表面近似为椭球面。
虽然多数情况下通过这样的近似处理后，能够得到符合实际情况的结果，但是对于极不规则的形态，比如：云朵、烟雾、树木等，传统的几何学就无能为力了。
如何描述这些复杂的自然形态？如何分析其内在的机理？这些就是分形几何学所面对和解决的问题。
三、问题解决（1）对一个等边三角形，每条边按照Koch曲线的方式进行迭代，产生的分形图称为Koch雪花。
编制程序绘制出它的图形，并计算Koch雪花的面积，以及它的分形维数。
（2）自己构造生成元（要有创意），按照图形迭代的方式产生分形图，用计算机编制程序绘制出它的图形，并计算其分形维数。
1、程序如下：function plottrkoch(a,k)％函数，a为迭代0次的三角形的边长，k为迭代次数p=[0 0;a 0;a/2 a/2*sqrt(3);0 0];n=3;A=[cos(pi/3) -sin(pi/3);sin(pi/3) cos(pi/3)];for s=1:kj=0;for i=1:nq1=p(i,:);q2=p(i+1,:);d=(q2-q1)/3;j=j+1;r(j,:)=q1;j=j+1;r(j,:)=q1+d;j=j+1;r(j,:)=q1+d+d*A';j=j+1;r(j,:)=q1+2*d;endn=4*n;clear pp=[r;q2];endmianji=sqrt(3)*(1+3*(1-(4/9)^k)/5)/4*a^2％计算迭代k次后的面积大小weishuD=log(4)/log(3)％计算维数plot(p(:,1),p(:,2))axis equal当k=1时当k=3时当k=5时因为是Koch曲线，所以维数d=ln4÷ln3=1.2619计算面积：假设最初的正三角形边长为1，则其周长为，面积为。
迭代·混沌·分形
迭代·混沌·分形柴文斌(四川省遂宁中学校629000)一、课例背景在20世纪下半叶，计算机的“魔杖”不断制造出新的数学分支，它最拿手的迭代计算引出了“混沌学”，接着又导致了分形几何的产生. 分形的思想和方法在模式识别，自然图象的模拟，信息讯号的处理，以及金融模型，艺术的制作等领域都取得了极大的成功.二、教学目标①本课例按《新课标》的要求，通过分形为载体，引起学生深厚的兴趣，在探究过程中，浅介数学新思想、新发展，同时让学生发现数学美，激发他们勇敢地追求美，主动地创造美，从而陶冶他们的情操，培养他们创新的精神.②总结平常练习过的从迭代、分形为背景数学试题，让他们用联系、发展的眼光，体会“背景深刻，方法独到”高考压轴题设计意图，明白“基础扎实，能力到位”明确要求.三、教学重点①应用计算机让学生感受分形图之美妙及形成数学原理.②分析分形为背景数学试题，形成高观点下审视数学问题.四、教学难点①迭代、混沌、分形定义度的把握.②Julia集、Mandelbrot集及其特征.五、教学过程(一)美丽的分形图形运用多媒体展放《孔雀开屏》等11幅分形艺术作品.师：这些美丽图形自然而优美，纷繁而有序，放射出诱人的色彩，在绚丽的色彩变化背后有几分神秘，似乎没有人会对这些图形无动于衷，你们相信，这些美妙的图形是运用数学知识，通过计算机构造出来的吗？是如何构造的呢？我们还得从函数迭代说起！(二)函数的迭代问题1：计算：①x n n sin lim ∞→ ②=∞→x n n cos lim 问题2：211n n x x +=-11=x轨道：1，0，－1，0，－1，……5.02=x轨道：0.5，―0.75，―0.4375，―0.80859，…―1，0，―1，0，－1问题3：①有没有这样一个初态把它代入211-+=n n x x ，结果不变吗？· ·A B251- 251+ ②618.11=x 写出系统轨道③619.11=x 写出系统轨道问题4：二次函数2)(z z f =进行迭代 ①i z 211=，写出系统轨道 ②i z +=11，写出系统轨道问题5：2)(z z f =且1||0<z求证：1|)(|0<z f证明：i y x z 000+=且1||0<z|2||)(|0020200i y x y x z F +-=20022020)2()(y x y x +-=22020)(y x +=20||z = 因此，在区域1||00<<z 中，1|||)(|00<<z z F ，这就意味着2)(z z F =的每一次迭代，即21n n z z =+都会使z 向靠近0的方向移动，我们说z 向0收敛，或是z 的吸引号，若1||0>z 近似于上面的结论，会发现，经过迭代z 会趋向于∞.若1||0=z ，很明显，z 是平面上单位圆上的点. 于是我们发现复平面上可分为两个区域，一个区域便落在其中的点向0吸引与逼近，而另一个区域便落在其中的点∞逃逸，它们分界线便是1)(0=z F 的单位圆，就是Julia 集.(三)混沌①C z z f +=2)(，0≠C 时，其吸引子不再是0，而是一个区域被吸进去的点会遍整个区域，我们称这个区域为混沌区. 同时，分界线不再是1|)(|0=z F 的单位圆，它是一个不规则不光滑的分界线，就像一个孤岛的海岸线一样.②《三五历经》中说：“天地混沌如鸡子，盘古生其中，万人千岁，天地开辟，阳清为天，阴浊为地，盘古在其中，一日九变；神于天，圣于地. 天月高一丈，地日厚一丈，盘古月长一丈，如此万人千岁，天数极高，地数极深，盘古极长.”③宇宙起源的问题.(四)分形不使系统发散的那么初态的集合组成“内集”，其他的“初态”组成“外集”，内集与外集的边界叫做Julia 集.问题6：运用多媒体展示：i z z f 12.0765.0)(2+-=(一个完全不连通)i z z f +=2)((连通) 特点：处处不光滑，自相似性、精细结构②Mandelbrot 集我们看到，当C 在复平面变化时它的Julia 集也在复平面内变化，而且这些集合可以分成连通与不连通两类. 如果参数C 所对应的Julia 集是连通的，我们就将这个C 染成黑色，否则染上白色，这样得到的黑色集就叫做以参数C 为元素的Mandelbrot 集.问题7：运用多媒体展示Mandelbrot 集，可以看出它有非常复杂的结构，这一结构的明显特征是一个主要心形图与一系列圆盘形的“芽苞”连接在一起，并且，每一个芽苞又被一细节更细小的“芽苞”所环绕，以至无穷. 同时，这些精细的芽苞分支都带有与整体曼德布罗特集相似的微型拷贝.(五)试题研究问题8：将一个单位正三角形一分为四，且挖去中间一个小正三角形，然后再上面三个小三角形中重复上面的步骤. 设初始三角形的面积为1. C n 、S n 分别表示第n 次操作各图形的周长和面积.①求C n 、S n 的表达式.②n 趋于无穷时，C n 、S n 趋于什么？问题9：记P 0表示面积为1的等边三角形，P k+1是对P k 进行如下操作得到：将P k 的每条边三等份，以每边中间部分的线段为边，向外作等边三角形，再将它中间部分的线段去掉，记S n 为曲线P n 所围成图形的面积.①求数列S n 的通项公式.②n n S ∞→lim 问题10：一种树形图为：第一层是一条与水平线垂直的线段长度为1，第二层与第一层线段的前端作两条与线段都成135°角的线段，长度为其一半，第三层按第二层的方法滚动，在第二层线段前端生成两条线段，重复前面的作法，作图到第n 层，称水平线到第n 层最高点的距离为到树形图的第n 层高度，试求：①树形图的第三层及第四层总高度②树形图的第n 层总高度h n③n n h ∞→lim问题11：⎪⎩⎪⎨⎧-+)1(221)(x x x f 121210≤<≤≤x x 定义 *∈=N n x f f f x f n ))(()(①求)152(2007f ②]}1,0[,)(|{15∈==x x x f x B求证：B 中至少含有9个元素问题12：如右图是某计算机的程序框图.(I)求打印出来的x 的值；(II)求打印出来的z 的值；(III)若将程序框图中的语句(9)“n=2007？”改为“94≥z ?”，则张三同学说这是死循环(即一直无限算下去而没有结果)，而李四说不会是死循环，你认为哪个同学说的正确？并说出你的理由.问题13：用牛顿迭代法求根.17世纪，牛顿创立了一种依靠简单计算求解方程根的方法.假设你知道某一方程0)(=x f 的近似解为0x ，此0x 接近于你还不知道的真正解x ，从而可以计算出相应的0)(=x f 及其导数0)(0=x f 的值. 由于0x 接近于x ，所以导数)(0x f '可近似写成00)()(x x x f x f --. 又因为0)(=x f 所以此导数为：000)()(x x x f x f --='] 于是有 )()(000x f x f x x '-=-则修正一次后的近似解为)()(001x f x f x x '-= 重复这个过程得到序列数n x ，它会从极快的速度收敛于此方程的真正解.请你用上述方法 81)(3-=x x f 6.00=x 时①求2x ，3x ，4x . ②100001|21|<-n x 时，n 的范围. 问题14：用多媒体展示基于牛顿迭代法的01=-n z 迭代图形.问题15：(角谷猜想) 任给一个自然数，若它是偶数则将它除以2；若它是奇数，则将它乘3再加1，反复这样运算，经有限步之后其结果必为1. 问题16：分形几何上物理学是怎样？六、课例设计反思：1．数学≠数学题. 数学教育，我想不仅要让学生认识到数学是一门科学，数学是工具，数学是技术，而且应当让他们认识到：作为人类精神、智慧与理性的最高代表之一，数学不仅是文化的重要组成部分，还在文化发展中占据着举足轻重的地位，数学是美的，数学是有意思的.2．Shirley(1986)提出，数学分为形式和非形式，应用和纯粹的，我们平常看到多数中小学讲授的数学知识是形式纯数学，这对学生形成完善的数学的文化观有缺陷，新课标模块设计也充分考虑到这一缺陷,本课例对非形式化教学，研究性学习作些探讨.。
数学实验综合实验报告
一、实验目的：1、初步认识迭代，体会迭代思想的重要性。
2、通过在mathematica 环境下编写程序，利用迭代的方法求解方程的根、线性方程组的解、非线性方程组的解。
3、了解分形的的基本特性及利用mathematica 编程生成分形图形的基本方法， 在欣赏由mathematica 生成的美丽的分形图案的同时对分形几何这门学科有一个直观的了解。
从哲理的高度理解这门学科诞生的必然性，激发读者探寻科学真理的兴趣。
4、从一个简单的二次函数的迭代出发，利用mathematica 认识混沌现象及其所 蕴涵的规律。
5、.进一步熟悉Mathematic 软件的使用，复习总结Mathematic 在数学作图中的应用，为便于研究数学图像问题提供方便，使我们从一个新的视角去理解数学问题以及问题的实际意义。
6、在学习和运用迭代法求解过程中，体会各种迭代方法在解决问题的收敛速度上的异同点。
二、实验的环境：学校机房，mathematica4环境三、实验的基本理论和方法：1、迭代（一）—方程求解函数的迭代法思想：给定实数域上光滑的实值函数)(x f 以及初值0x 定义数列1()n n x f x +=， ,3,2,1,0=n ， (1)n x ， ,3,2,1,0=n ，称为)(x f 的一个迭代序列。
（1）方程求根给定迭代函数)(x f 以及初值0x 利用（1）迭代得到数列n x ， ,3,2,1,0=n .如果数列收敛到某个*x ，则有)(**x f x =. (2)即*x 是方程)(x f x =的解。
由此启发我们用如下的方法求方程0)(=x g 的近似解。
将方程0)(=x g 改写为等价的方程)(x f x =， (3) 然后选取一初值利用（1）做迭代。
迭代数列n x 收敛的极限就是方程0)(=x g 的解。
为了使得迭代序列收敛并尽快收敛到方程0)(=x g 的某一解的条件是迭代函数)(x f 在解的附近的导数将的绝对值尽量小，因此迭代方程修订成x x f x h x )1()()(λλ-+== (4) 选取λ使得|)(|x h '在解的附近尽量小. 为此, 我们可以令,01)()(=-+'='λλx f x h得)(11x f '-=λ. 于是 1)()()(-'--=x f x x f x x h . 特别地,如果取x x g x f +=)()(, 则可得到迭代公式 .,1,0,)()(1 ='-=+n x g x g x x n n n n (5) （2）线性方程组的数值解的迭代求解理论与矩阵理论给定一个n 元线性方程组⎪⎩⎪⎨⎧=++=++,,1111111n n nn n n n b x a x a b x a x a (6)或写成矩阵的形式,b Ax = (7) 其中)(ij a A =是n 阶方阵,T n x x x x ),,(21 =及T n b b b b ),,,(21 =均为n 维列向量.熟知，当矩阵A 的行列式非零时，以上的方程组有唯一解.如何有效，快速地寻求大型的线性方程组的数值解释科学工程计算中非常重要的任务.而迭代法常常是求解这些问题的有效方法之一。
分形、混沌与迭代
试验十二 分形、混沌——迭代一、试验目的：1、Koch 曲线、Sierpinski 三角形、Cantor 集的计算机实现2、掌握用迭代、递归生成分形3、用Matlab 观察分岔与混沌现象二、分形相关程序：1、从一条直线段开始，将线段中间的三分之一部分用一个等边三角形的两边代替，形成山丘形图形如下在新的图形中，又将图中每一直线段中间的三分之一部分都用一个等边三角形的两条边代替，再次形成新的图形如此迭代，形成Koch 分形曲线。
算法分析：考虑由直线段（2个点）产生第一个图形（5个点）的过程。
图1中，设1P 和5P 分别为原始直线段的两个端点，现需要在直线段的中间依次插入三个点2P ，3P ，4P 。
显然2P 位于线段三分之一处，4P 位于线段三分之二处，3P 点的位置可看成是由4P 点以2P 点为轴心，逆时针旋转600而得。
旋转由正交矩阵 ⎪⎪⎪⎪⎭⎫ ⎝⎛-=)3cos()3sin()3sin()3cos(ππππA 实现。
算法根据初始数据（1P 和5P 点的坐标），产生图1中5个结点的坐标。
结点的坐标数组形成一个25⨯矩阵，矩阵的第一行为1P 的坐标，第二行为2P 的坐标……，第五行为5P 的坐标。
矩阵的第一列元素分别为5个结点的x 坐标，第二列元素分别为5个结点的y 坐标。
进一步考虑Koch 曲线形成过程中结点数目的变化规律。
设第k 次迭代产生的结点数为k n ，第1+k 次迭代产生的结点数为1+k n ，则k n 和1+k n 中间的递推关系为341-=+k k n n 。
实验程序及注释：p=[0 0;10 0]; %P 为初始两个点的坐标,第一列为x 坐标,第二列为y 坐标 n=2; %n 为结点数A=[cos(pi/3) -sin(pi/3);sin(pi/3) cos(pi/3)]; %旋转矩阵for k=1:4d=diff(p)/3; %diff 计算相邻两个点的坐标之差,得到相邻两点确定的向量%则d 就计算出每个向量长度的三分之一,与题中将线段三等分对应m=4*n-3; %迭代公式q=p(1:n-1,:); %以原点为起点,前n-1个点的坐标为终点形成向量p(5:4:m,:)=p(2:n,:); %迭代后处于4k+1位置上的点的坐标为迭代前的相应坐标p(2:4:m,:)=q+d; %用向量方法计算迭代后处于4k+2位置上的点的坐标p(3:4:m,:)=q+d+d*A'; %用向量方法计算迭代后处于4k+3位置上的点的坐标p(4:4:m,:)=q+2*d; %用向量方法计算迭代后处于4k 位置上的点的坐标 n=m; %迭代后新的结点数目endplot(p(:,1),p(:,2)) %绘出每相邻两个点的连线axis([0 10 0 10])实验数据记录：由上面的程序，可得到如下的Koch 分形曲线：2、由四边形的四个初始点出发，对于四边形的每条边，生成元如下：可得到火焰般的图形。
3-迭代与分形
试验二迭代与分形一、实验目的与要求1．了解分形几何的基本情况；2．了解通过迭代方式，产生分形图的方法；3．了解matlab软件中简单的程序结构；4．掌握matlab软件中plot, fill等函数的基本用法；二、问题描述几何学研究的对象是客观世界中物体的形状。
传统欧氏几何学的研究对象，都是规则并且光滑的，比如：直线、曲线、曲面等。
但客观世界中物体的形状，并不完全具有规则光滑等性质，因此只能近似当作欧氏几何的对象，比如：将凹凸不平的地球表面近似为椭球面。
虽然多数情况下通过这样的近似处理后，能够得到符合实际情况的结果，但是对于极不规则的形态，比如：云朵、烟雾、树木等，传统的几何学就无能为力了。
如何描述这些复杂的自然形态？如何分析其内在的机理？这些就是分形几何学所面对和解决的问题。
三、问题分析在我们的世界上，存在着许多极不规则的复杂现象，比如：弯弯曲曲的海岸线、变化的云朵、宇宙中星系的分布、金融市场上价格的起伏图等，为了获得解释这些极端复杂现象的数学模型，我们需要认识其中蕴涵的特性，构造出相应的数学规则。
曼德尔布罗特（Mandelbrot）在研究英国的海岸线形状等问题时，总结出自然界中很多现象从标度变换角度表现出对称性，他将这类集合称作自相似集，他发现维数是尺度变换下的不变量，主张用维数来刻划这类集合。
Mandelbrot 将这类几何形体称为分形(fractal)，意思就是不规则的、分数的、支离破碎的，并对它们进行了系统的研究，创立了分形几何这一新的数学分支。
Mandelbrot 认为海岸、山峦、云彩和其他很多自然现象都具有分形的特性，因此可以说：分形是大自然的几何学。
分形几何体一般来说都具有无限精细的自相似的层次结构，即局部与整体的相似性，图形的每一个局部都可以被看作是整体图形的一个缩小的复本。
早在19世纪就已经出现了一些具有自相似特性的分形图形，比如：瑞典数学家科赫（von Koch）设计的类似雪花和岛屿边缘的一类曲线，即Koch曲线；英国植物学家布朗通过观察悬浮在水中的花粉的运动轨迹，提出来的布朗运动轨迹。
分形
k
即是所谓的Koch曲线 Koch曲线的修改规则R是将每一条直线段F0用一 条折线F1替代，称F1为该分形的生成元。分形图形的 基本特性完全由生成元决定。
分形树 分形树的生成方法是：将一条线段三等分，在 等分点上各画一条长度为原线段长度三分之一的线 段，该线段与元线段成固定的夹角。
尽管分形几何的提出只有三十年左右的时间， 但它已经在自然科学的各个领域如数学、物理、化 学、地理、天文、材料、生命乃至经济、社会、艺 术等极其广泛的领域有着广泛的应用。 该实验的目的是以迭代的观点介绍分形的基本 特性以及生成分形图形的基本方法。
生成元产生的分形图形 由IFS(迭代函数系)所生成的分形图形
Sierpinski垫 w a 1 0.5 2 0.5 3 0.5 b 0 0 0 c 0 0 0 d e f P 0.5 0 0 0.333 0.5 0.25 0.433 0.333 0.5 0.5 0 0.333
用IFS迭代绘制分形图形的方法如下： 首先设定迭代的可视区域为 V {( x , y ) | x min x x max , y min y y max } 再按分辨率大小的要求将V分成a*b的网格 网格点为 ( x i , y i )
1 生成元
由生成元产生的分形是一种规则分形，是数学 家按一定规则构造出来的，相当于物理中的模型。 这类图形的构造方式都有一个共同的特点：
最终图形F都是按照一定的规则R通过 对初始图形F0不断修改得到的。
最具代表性的图形有Koch雪花 曲线、分形树
Koch曲线的构造方式是： 给定一条直线段F0，将该直线三等分，并将中间 的一段用以该线段为边的等边三角形的另外两边代替， 得到图形F1，然后，再对图形F1中的每一小段都按上 述方式修改，以至无穷。则最后得到的极限曲线
数学实验迭代：分形
迭代：分形姓名：学号：班级：数学与应用数学4班实验报告实验目的：以迭代的观点介绍分形的基本特性以及生成分形图形的基本方法，使读者在欣赏美丽的分形图案的同时对分形几何这门学科有一个直观的了解，并从哲理的高度理解这门学科诞生的必然，激发读者探寻科学真理的兴趣。
实验环境：Mathematica软件实验基本理论和方法：在19世纪末及20世纪初，一些数学家就构造出一些边界形状极不光滑的图形，而这类图形的构造方式都有一个共同的特点，即最终图形F都是按照一定的规则R通过对初始图形不断修改得到的。
其中最有代表性的图形是Koch曲线，Koch曲线的构造方式是：给定一条直线段，将该直线段三等分，并将中间的一段用以该线段为边的等边三角形的另外两条边代替，得到图形，然后再对图形中的每一小段都按上述方式修改，以至无穷。
则最后得到的极限曲线即是所谓的Koch曲线。
生成元：Koch曲线的修改规则R是将每一条直线段用一条折线代替，我们称为该分形的生成元。
分形的基本特性完全由生成元确定，因此，给定一个生成元，我们就可以生成各种各样的分形图形。
Julia集绘制方法：（1）设定初值p，q，一个最大的迭代次数N，图形的分辨率的大小a，b，和使用的颜色数(如K=16)（或者给定灰度级L）；（2）设定一个上界值；（3）将矩形区域分成的网格，分别以每个网格点，，，，作为初值利用riter做迭代（实际上，只需对满足的初值点做迭代）。
如果对所有，，则将图形的像素点用黑色显示，否则，如果从迭代的某一步开始有，则用modK种颜色显示相应像素（或者用相应的灰度级显示）。
Mandelbrot集绘制方法：设定一个最大的迭代次数N，图形的分辨率的大小a，b，和使用的颜色数(如K=16)（或者给定灰度级L）；（2）设定一个上界值；（3）将矩形区域分成的网格，分别以每个网格点，，，，作为参数值利用riter做迭代（实际上，只需对的初值点做迭代），每次迭代的初值均取为。
ppt02迭代分形图形
%程序1 u=2.6:0.001:4; x=0.6; for j=1:150, x=u.*(x-x.^2); end for i=1:100 x=u.*(x-x.^2); plot(u,x,’r.’) hold on; end 程序运行时间 0.28秒。不保留旧 的X 值，而是直接用它画图。
%程序2 u=2.6:0.001:4; X=ones(250,1401); X(1,:)=0.6*X(1,:); for j=1:250 X(j+1,:)=u.*(X(j,:)-X(j,:).^2); end plot(u,X(150:end,:),’r.’) 运行时间0.15秒，比程序1快。 保 留 所 有 X 值 ， 每 次 计 算 的X 值 生成矩阵的一行元素，最后用矩 阵X 的后150行作图，程序可读性 强。
u=uuu; subplot(3,3,n+1); end
plot(u)， axis([-0.5,0.5,0,1])
3.1.3 树程序3 以矩阵元素赋值的方法计算新点坐标，减少一重循环。 new=[0,i]; subplot(3,3,1); plot(new); axis([-0.5,0.5,0,1 ]); for k=1:8 old=new; n=length(old)/2-1; diff=(old(2:2:end)-old(1:2:end-1))/3; p1=old(1:2:end-1)+diff; p2=p1+diff; lp=p1+diff*(sqrt(3)/2+1/2*i); rp=p2+diff*(sqrt(3)/2-1/2*i); new(2:10:10*n+2)=p1; %第一次循环计算第2,12,22,...点 new(3:10:10*n+3)=p1; %第一次循环计算第3,13...点 new(4:10:10*n+4)=lp; new(5:10:10*n+5)=p1; new(6:10:10*n+6)=p2; new(7:10:10*n+7)=p2; new(8:10:10*n+8)=rp; new(9:10:10*n+9)=p2; new(10:10:10*(n+1))=old(2:2:end); %10,20,... new(1:10:10*n+1)=old(1:2:end-1); %1,11,21.... subplot(3,3,k+1)， plot(new) axis([-0.5,0.5,0,1 ]); end 3.1.4 树程序4(运行时间0.21秒) 每个图形缩小为原来的1/3,再移动位置, 添加到原图形。参看右边的电影。 u=[0,i]; subplot(3,3,1); plot(u) for k=2:8 m=u/3;
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
@8 D @ D @< 8 < @ @8 @ DD @D@L D @D H @ D @ 8 @ D @D @ @ @D@D D D< < @< D< D D @@@ D 88 @ @D D
@@
;
D
4、龙曲线
dragon ptlist_List := Block
ln01 =
ListPlot Nest dragon, ln01, 12 , PlotJoined ® True, AspectRatio ® 1 1.6, Axes ® None
@ D @ 8 @D < @HHL @D @ D L 8 @D @D < @ @D H@ L D D D D8< < 8 D @@ D
, i, pnum = Length ptlist 3 , 3 i+ 1 , ; , 3 i+ 2 3 i+ 3 Show Graphics tmp , AspectRatio ® 1 GoldenRatio
pol =
-1, 0 , 1, 0 , 0, Sqrt 3
showsierpinski Nest redosierpinski, pol, 6
@ D@ D @8 8< < @ @8 @ D @@ D DD H@D @D @D @ D D @L H@D @D @D @D L H@ D L @D @D @D @D H@D @D @D @D L H@D @D @D @D L H@ D @ D @D L @D @< DD D
• 分形的维数
1、相似维数：设分形 F 是自相似的，即 F 由 m 个子集构成，每个子集放大 c 倍后同 F一样 ，则定义 F 的维数为
d ( F ) log( m) / log( c)
例如，对于Cantor集， d ( F ) log 2 / log 3 对于Van Koch 雪花曲线，d ( F ) log 4 / log 3
实验十二 迭代 (2) --分形
实验内容
• 什么是分形？ • 图形迭代 • 函数迭代 • IFS迭代 • 分形的应用
1、什么是分形
• 分形发展简史 欧氏几何、解析几何、微分几何—正则 微积分，复变函数---光滑 反例 1，Cantor集合
F0
F1 F2
Cantor 集合 F中点数不可数（比有理数 还多！），但其区间长度为零！
ptlist
@8 D @ D 8 @< 8 < @H @D @L DD 8 @D @D D @< @ @ @@ D 8 < D @ @ 8 @ D@ D D @ D @ @ D D @ D @ D @ @@ D DD @ D @ D @ D
, i, j, pnum = Length ptlist , tmp1, tmp2, d = ,
反例 2，Weierstrass函数
W ( x)
n 0 ( s 2) n
sin( x)
n
其中 1<s<2 且 1 ，W(x) 是处处连续、
处处不可微的函数。对应 s=1.4, 2
的图象是
lambda = 2; nmax = 30; s = 1.5; Plot Sum lambda^ k, 1, nmax , x, - 10, 10
- ptlist i 5; - ptlist i + 1, 2 ,
i + 1, 1
ptlist
d = Join d,
-1, - 1, 1, - 1, 1, 1, -1, -1, 1, 1, - 1, -1 ; tmp = Join tmp, ptlist i , ptlist i - tmp2, ptlist i - 2 tmp2, ptlist i + tmp1 - 2 tmp2, ptlist i + tmp1 - tmp2, ptlist i + 2 tmp1 - tmp2, ptlist i + 2 tmp1 - 2 tmp2, ptlist i + 2 tmp1 - 3 tmp2, ptlist i + tmp1 - 3 tmp2, ptlist i + tmp1 - 4 tmp2, ptlist i + tmp1 - 5 tmp2, ptlist i + 2 tmp1 - 5 tmp2, ptlist i + 2 tmp1 - 4 tmp2, ptlist i + 3 tmp1 - 4 tmp2, ptlist i + 4 tmp1 - 4 tmp2, ptlist i + 5 tmp1 - 4 tmp2
2;
- ptlist i, 1
+ ptlist i + 1
2 + tmp1,
; tmp
;
0, 0 , 1, 2
;
5、Hilbert曲线
Hilbert ptlist_List := Block tmp =
For i = 1, i < pnum, i ++, tmp1 = tmp2 = ptlist i+ 1 ptlist i, 2
, i, pnum = Length ptlist 3 i +1 , 3 ,
+ ptlist + ptlist + ptlist + ptlist + ptlist + ptlist
3 i +2 3 i +3 3 i +2 3 i +3 3 i +3 3 i +3
2, 2, 2, 2, 2, 2,
ptlist
i+ 1 1, 0
; tmp ;
0, 0 ,
;
3、Sierpinski地毯
redosierpinskiቤተ መጻሕፍቲ ባይዱptlist_List := Block tmp = For i = 0, i < pnum, i ++, tmp = Join tmp, ptlist ptlist ptlist ptlist ptlist ptlist ptlist 3 i+ 1 3 i+ 1 3 i+ 1 3 i+ 2 3 i+ 1 3 i+ 2
tmp = Join tmp,
ln01 =
ListPlot Nest redominkowski, ln01, 3 , PlotJoined ® True, AspectRatio ® 1 GoldenRatio, Axes ® None
@ D @D @8 8< < @8 @D @ D @D @DD D @D @< @8 @D @ D@ @D@D DD @D@D @ @D D @ @D @D @@D D @D @D @@D D @D @D @@D D @D @D @@D D @D @D @@D D @D @D D< @ @ DD < D< < 8@ 88 @ DD
- ptlist i, 1 5; If direct i Š 1, d = Join d, 1, 1, - 1, -1, 1, 1, - 1, - 1, 1, - 1, 1, 1, -1, 1, -1, 1, 1, - 1, 1, - 1, -1, 1, 1, - 1, - 1 ; tmp = Join tmp, ptlist i , ptlist i + tmp1, ptlist i + 2 tmp1, ptlist i + 2 tmp1 + tmp2, ptlist i + tmp1 + tmp2, ptlist i + tmp1 + 2 tmp2, ptlist i + 2 tmp1 + 2 tmp2, ptlist i + 2 tmp1 + 3 tmp2, ptlist i + tmp1 + 3 tmp2, ptlist i + 3 tmp2, ptlist i + 4 tmp2, ptlist i + tmp1 + 4 tmp2, ptlist i + 2 tmp1 + 4 tmp2, ptlist i + 3 tmp1 + 4 tmp2,
tmp = ptlist, tmp1, i, pnum = Length ptlist ptlist i, 2 , For i = 1, i < pnum, i ++, tmp1 = -1 ^ i - 1 ptlist tmp = Insert tmp, ptlist i 2i i + 1, 1
- ptlist i + 1, 2 ,
• 对于一条直线段，将它等分，每段长度为 原来的1/N，共分为N段。 • 将一个正方形每边等分成N段，共有N2个小 正方形。 • 将一个立方体每边等分成N段，共有N3个小 立方体。 • 一般地，设一图形可分解为m个与之相似的 子图形，每个子图形是原来的1/c. 则图形的 维数D满足：cD=m.
2、盒子维数：设 F R是有界集合， 其中 R 是正方形。将 R 分成边长为 的 子正方形。记 N ( ) 为子正方形中包含 F 中点的子正方形的个数。定义 F 的盒子 维数为 ln N ( ) d ( F ) lim 0 ln( 1 / ) 例如，对于 Weierstrass处处连续、处处 不可微的函数，其分形维数为 s.
Koch 曲线的关系，提出了一门描述大自
然的几何形态的学科---分形(Fractal)