分形M集绘制方法和其实验结果(11年使用,08级,试卷).docVIP

数学试验课试卷 分形M集的绘制方法及其实验结果 本次数学实验试题的题目是“分形M集的绘制方法及其实验结果”（注：M集是Mandelbrot集的简称），其具体内容与绘制方法可参见《数学实验》（李尚志）第122至124页，更进一步内容可参见《分形几何---数学基础及其应用》（作者：肯尼思.法尔科内，曾文曲、刘世耀译，东北工学院出版社）。 对于给定的复数，考虑如下函数迭代公式： （1） 对于不同的复值参数，当时，此迭代所产生的迭代序列可能有界，也可能无界，令是使迭代序列有界的所有参数值所构成的集合，即有 ={ | 所有使公式（1）迭代序列有界时的值} 当时，一般称为复平面上的Mandelbrot集(即M集)。 如何画出M集在复平面上某个矩形区域内的图像？选定复平面上的某个矩形区域，在此区域内任取一个，并取初值，取迭代函数为，即进行迭代，则此迭代对不同的，有如下三种结果：①发散到无穷大；②收敛到0；③其它情况。实际计算表明，当此迭代所得到复数值的模大于2时，此迭代一般发散到无穷大（这只对成立），但是，不同的值发散到无穷大的速度并不相同。而此迭代对某个收敛到0时，它将快速收敛到0。为了画出M集，根据集合的定义，我们只要用不同的颜色画出给定的矩形区域内迭代发散到无穷大点、收敛到0的点、其它情况点即可，则除了发散到无穷大的点，另外两种类型的点的并集就是M集。为了增加图像的趣味性，我们将此区域内所有的发散到无穷大的速度的快慢用不同的颜色来标记，这样就可将发散到无穷大的的点与另外两种类型点区分开来，因此就可画出一张分形图像（由于计算机计算速度的限制，对，若模大于某一个固定的值，就可认为迭代发散到无穷大，其它迭代函数则可根据实际情况确定的值）。 下面是实现此程序的代码 其中选项ColorFunction是颜色函数，若变为如下，则图形会更加好看。 在上面的函数M[x,y]中，只要z的模大于2，就认为是发散到无穷大，M[x,y]的返回值k是绘图所用的颜色数目（上面的例子是用50种颜色画图，你可以根据实际情况修改此数值），PlotPoints-200表示绘制一个200X200的图像，Mesh-False表示不要网格，ColorFunction表示mathematica绘图时，使用什么样的颜色函数，关于Hue，RGBColor函数的详细说明，请参见mathematica的帮助文件。 但是，上面的程序有一个问题，就是运行速度太慢，为此我们使用mathematica的Compile函数（参见mathematica的帮助文件），将此程序编译后运行，速度会成倍提高（首次运行除外），下面是程序代码。请注意其中的一行是“k]”，不能写成“k];” 不同的矩形区域，不同的color与zmax（注：对Abs[z]zmaxkcolor），以及不同的ColorFunction，你会得到不同的图形。例如 此程序运行后，将会得到如下图形： 而下面的程序生成的图形是： 以上的图形，其绘图的值实际上迭代的次数，如果我们我们选取迭代函数为一个任意复函数，返回值也不用迭代次数，而改用最后的迭代复数值的模，则图形就会丰富多彩，例如 经过试验，mathematica的Compile函数有点小问题，如果在申请局部变量z时，写z=0+0*I，则有时候会出现语法问题，而改用z=x+y*I则不出现问题，对此种改动，两者运行结果一样。下面是另一个例子。 其图形是： 下面也是用mathematica生成的图形，不过是只有黑白两种颜色的图形，其生成方法是，对x,y，根据函数的迭代次数k，将M[x,y]变为只取0或者1两个值的函数，例如，若k为偶数，M[x,y]=0，若k为奇数，m[x,y]=1，或者是，若kK（K为一个选定的整数），则M[x,y]=0，否则m[x,y]=1，然后用DensityPlot绘图。如果生成的是含有0与1的数据矩阵，可用ListDensityPlot绘图，此种绘图方法比DensityPlot速度快。下面是用此种方法生成的三个图形。 （图1） （图2） （图3） 关于本次数学实验试题的内容要求 本次数学实验试题的题目是“分形M集的绘制方法及其实验结果”，下面是本次数学实验考试的试题内容要求： （1）对分形理论及其分形图像的生成方法，做简要介绍。 （2）在复平面上，选定一些区域，画出M集的图形，根据这些图形的形状，你能得出哪些猜测。 （3）将M集的迭代函数替换为其它你能想到的函数，使用上面所提及的方法，或者是你所想到的新方法，画出其分形图形，图形的个数不限，其评判标准是图形视觉上的好坏。 （4）进一步地，你能够通过实验，找到一些或者一类迭代函数，它们所生成的图像，能够模拟自然