以前用VB,还有C++写了画Mandelbrot分形(Fractal)图的程序, 但代码都比较长。新学Python, 发现可以用非常段的代码就可以画出来, 而且速度很快, 分享如下:
import numpy as np
from functools import reduce
import matplotlib.pyplot as plt
x,y=np.ogrid[-2.2:0.8:200j,-1.3:1.3:200j]
c=x+y*1j
maxIt = 30
z = reduce(lambda x, y: x**2 + c, [1] * maxIt, c)
plt.imshow(np.angle(z).T)
plt.show()
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)