| 分形与混沌
自然界的很多事物，例如树木、云彩、山脉、闪电、雪花以及海岸线等等都呈现出传统的几何学不能
描述的形状。这些形状都有如下的特性：
有着十分精细的不规则的结构
整体与局部相似，例如一根树杈的形状和一棵树很像
分形几何学就是用来研究这样一类的几何形状的科学，借助计算机的高速计算和图像显示，使得我们
可以更加深入地直观地观察分形几何。
Mandelbrot集合
Mandelbrot(曼德布洛特)集合是在复平面上组成分形的点的集合。
Mandelbrot集合的定义(摘自维基百科) Mandelbrot集合可以用下面的复二次多项式定义：
其中c是一个复参数。对于每一个c，从z=0开始对函数 fc(z)
进行迭代。 序列 (0,fc(0),fc(fc(0)),fc(fc(fc(0))),…) 的值或者延伸到无限大，或者只停留在有限半径的圆盘
内。 Mandelbrot集合就是使以上序列不发散的所有c点的集合。
从数学上来讲，Mandelbrot集合是一个复数的集合。一个给定的复数c或者属于Mandelbrot集合，或
者不是。
用程序绘制Mandelbrot集合时不能进行无限次迭代，最简单的方法是使用逃逸时间(迭代次数)进行绘
制，具体算法如下：
判断每次调用函数 fc(z) 得到的结果是否在半径R之内，即复数的模小于R
记录下模大于R时的迭代次数
迭代最多进行N次
不同的迭代次数的点使用不同的颜色绘制
输出图如下:
代码如下:
#
-*- coding: utf-8 -*-
# by whyx 2016/6
import numpy as np
import pylab as pl
from math import log
import time
from matplotlib import cm
def iter_point(c):
z = c
for i in xrange(1, 100): # 最多迭代100次
if abs(z)>2.0: break # 半径大于2则认为逃逸
z = z*z+c
return i # 返回迭代次数
escape_radius = 2.0
iter_num = 100
def smooth_iter_point(c):
z = c
for i in xrange(1, iter_num):
if abs(z)>escape_radius: break
z = z*z+c
absz = abs(z)
if absz > 2.0:
mu = i - log(log(abs(z),2),2)
else:
mu = i
return mu # 返回正规化的迭代次数
def draw_mandelbrot(cx, cy, d):
# 绘制点(cx, cy)附近正负d的范围的Mandelbrot
(x0, x1, y0, y1) = (cx-d, cx+d, cy-d, cy+d)
y, x = np.ogrid[y0:y1:200j, x0:x1:200j]
c = x + y*1j
start = time.clock()
mandelbrot = np.frompyfunc(smooth_iter_point,1,1)(c).astype(np.float)
print "time=",time.clock() - start
#时间单位 秒
pl.imshow(mandelbrot, cmap=cm.Blues_r, extent=[x0,x1,y0,y1])
pl.gca().set_axis_off()
x,y = 0.27322626, 0.595153338
pl.subplot(231)
draw_mandelbrot(-0.5, 0, 1.5)
for i in range(2,7):
pl.subplot(230+i)
draw_mandelbrot(x, y, 0.2**(i-1))
pl.subplots_adjust(0.02, 0, 0.98, 1, 0.02,
0)
pl.show() |
|
