分形图形学实验指导.docxVIP

分形图形学实验指导 ———————————————————————————————— 作者 : ———————————————————————————————— 日期： ? 分形图形学实验指导 实验一 二维空间上的分形图形生成 实验目的 Ｍaｎdelb ｒ ot 集与 Jｕlia 集的计算机实现 掌握用 L 系统语言生成分形 实验内容及步骤 1. 编写程序生成 Maｎdelbrot 集 在复迭代中影响最大的当属迭代ｚ→ｚ＾ 2＋ c，实际上它只是形式更一 般的复解 析迭代 z＿ ( ｎ＋ 1)=F(z_n)+c 的一种， F 是一个非线性函数。显然 z→z＾ 2+c 也是最 简单的一种 , 它 在复迭代中的地位相当于逻辑斯蒂映射ｘ _(n+ １ )= ａ x_n(1- ｘ_n）在 实迭代中的地位 ( 见第八章）。 考虑一般形式的 F，令ｚ =x+iy,c=c_(  Ｘ）＋ ic_(Y)  ，其中  i  表示虚数  , ｉ =SQRT(-1) 。 分离实部与虚部，具体化迭代关系便有  : x→f( ｘ ,y)+c_ （Ｘ ),y →g（ x, ｙ） +c_( Ｙ ). 通常所说的 M集是迭代二次函数 z→ z^2＋ c 产生的 , 此函数具 体化就是 x→ x^２ - y^2+c＿ ( X） , ｙ→２ ｘ y＋ c_( Y） . 其中 z＝ x+i y， c=ｃ_( X)+i c_( Y ), 以横轴 ｘ记录实数的实部，以纵轴 y 记录实数 的虚部。 M集合实际上是常数 c=（ c_( Ｘ), c_( Ｙ） ) 构成的图象。让 c 从屏幕左上角 开始变化，逐行增加 , 一直变到屏幕右下角。如果取的区域是２０ 0×2 00，则一共要 计算 40，000 个点 , 把计算的结果用不同的颜色标记下来，就得到一幅图象，这 就是 集。对于不同的 c 值 , 如何得到表征迭代性质的不同的结果呢 ?  M 容易知道 , 无穷远处肯定是迭代的一个吸引子， 即对于复平面上相当多的初始条件 , 迭代最终都跑到无穷远处。 但研究发现， 在原点附近还存在一个奇特的区域， 在迭代过程中此区域永远不会跑掉。在非严格的意义上，这个不变的集合就是Ｍ集 , 我们的主要任务就是画出这个集合的边界——实际上边界是分形曲线 , 极其复杂，Ｍ集图象的全部魅力就在这里。 在复平面上 , 我们以原点 (0, ０) 作为参考点，观察迭代过程是否远离原点，以及逃 离原点的速度如何。为此规定一个距离函数 D=x^2+ｙ^２ , 其实 D有许多不同的取法 , 以上取法是最普通的。 可以看出 , 如果 D较大 , 表 明迭代 点离原点较远 , 如果 Ｄ 较小，表明迭代点离原点较近。假设对于任何一个 c, 迭代都从 x_0， y＿ 0) ＝(0, ０ ) 开始 , 我们观察迭代点列 x_1，ｙ_1) ，（ｘ_2，y_2),( ｘ_３ , ｙ_3),，（ x_（ 100），y＿ (10 ０)) ， 的变化状况。每一次都计算一下 D值，即该点与原点 (0 ， 0) 的距离 ( 为方便计，这 里 实际上计算的是距离的平方 ) 。取一个参考距离 Ｒ，不妨取得大一些 , 比如 R=40( 实 际上取２ 0 就足够了 ) 。现在想知道迭代多少次后实际的距离 D大于 R 。在迭代过程中 如果 D小于 R，则继续让计算机迭代 , 要规定一个上限， 比如 30０次。如果迭代了 300 次后结果仍然不跑掉 ( 即 D 仍然小于 Ｒ ), 则可以近 似认为此点属于 集合。 M 对于迭代次数小于 300 次的情况 , 如果迭代 10 次 D就大于 R, 则标记 c 点为白色 ; 如果迭代 35 次 D 开始大于 Ｒ ，则标记 c 点为红色；如果迭代 １10 次 开始大于 Ｒ ， D 则标记 c 点为蓝色，等等。 2. 编写程序生成 Juli ａ集 一个Ｍ集可以对应无数种 J 集, 实际上 M集就是所有 J 集的浓缩 。M集不同部位 的形状 , 反映了对应于该处的 J 集的形状， 于是用 M集 可以对 J 集进行分类， 至少在计 算机图形学的层次上可以这样说。 计算 J 集时，初始迭代点就不能总取 （ 0, ０ ) 了，而是要根据实际位置取实际的 (x  ， 坐标值。仍以迭代 z→z^２ +c 为例说明。先取 定一个 c 值, 比如 c=(1.02 ２ 1,0.2433), 迭代关系化为 x→x^2 -y^2+1.0 ２ 2１,y →2xy+0.243 ３ . 从屏幕左上角算起 , 逐行计算 , 一直算到屏幕右下角。 当然，也可以不取整个屏幕那么大，只选一个２ 0０× 200 的小区域做。 标色的原理与上面讲  M集合时完全类似  , 此略。