我试图画曼德布罗特集时，我的数学出了什么问题？

Question

我目前正在尝试用python编写一个程序，结果是绘制mandelbrot集。分形的形状看起来不错，但轮廓一点也不精确。我正试图弄清楚这个问题是来自我的数学还是我的代码

我试图提高最大迭代次数（我称之为i），看看计算是否太低，但变化不大。

import numpy as np
from PIL import Image

taille = 1000
nb_points = 500
centre_x = taille/2
centre_y = taille/2

fractale = np.zeros((taille, taille, 3), dtype = np.uint8)

for x in range(-nb_points,nb_points):
    x = float(x)
    x = x/250

    for y in range(-nb_points,nb_points):
        y = float(y)
        y = y/250

        X = float(0)
        Y = float(0)
        i = 0

        module_carre = 0

        while module_carre < 4 and i <20 :
            X = float(X**2 - Y**2 + x)
            Y = float(2*X*Y + y)
            i += 1 
            module_carre = float(X*X + Y*Y)

        if module_carre < 4:
            i=0

        couleur = liste_couleur[i]

        fractale[centre_x + x*250,centre_y + y*250] = couleur
imgpil = Image.fromarray(fractale, 'RGB')
imgpil.save("resultat.jpg")

我是法国人，这就是为什么你可能不明白一切。我没有粘贴所有关于定义不同深浅的蓝色等的线条...我不明白为什么我的轮廓这么糟糕。我相信它来自我的数学，但我没有看到任何错误。这是我第一次发布堆栈溢出，我还不明白所有的工具。我没有设法添加输出的图像对不起我的英语，我希望你能帮我修复我的代码

Answer 1

主要问题是：

X = float(X**2 - Y**2 + x)
Y = float(2*X*Y + y)

在Y的计算中，您希望使用X的旧值，但它已经更新，不再可用。你可以这样做：

new_X = float(X**2 - Y**2 + x)
Y = float(2*X*Y + y)
X = new_X

或者您可以并行更新它们：

X, Y = float(X**2 - Y**2 + x), float(2*X*Y + y)

此外，您不需要调用float，因此您真正需要的是：

X, Y = X**2 - Y**2 + x, 2*X*Y + y

顺便说一下，Python确实有一个内置的复数类，所以您可以将z和Z初始化为z=复杂（x， y）和Z=复杂（0.0,0.0）。那么你的循环体会有：

Z = Z**2 + z

我还建议增加最大迭代次数。将其限制在20将给出非常低的分辨率。当生成高分辨率图像时，我通常至少使用1000，至少对于边界附近的点。