Python代码绘制Julia分形集 - 详细解析与示例

以下是绘制Julia分形集的Python代码：

import numpy as np
from PIL import Image

def julia(c, max_iter=255, size=(1000, 1000)):
    '绘制Julia分形集'
    # 创建一个空白图像，并将其转换为numpy数组
    img = Image.new('RGB', size, (0, 0, 0))
    pixels = img.load()
    w, h = size
    # 设置Julia集参数
    re_min, re_max = -2, 2
    im_min, im_max = -2, 2
    cre, cim = c
    # 计算每个像素点是否属于Julia集
    for x in range(w):
        for y in range(h):
            zx = re_min + x * (re_max - re_min) / (w - 1)
            zy = im_min + y * (im_max - im_min) / (h - 1)
            i = 0
            while i < max_iter:
                zxn = zx * zx - zy * zy + cre
                zyn = 2 * zx * zy + cim
                zx, zy = zxn, zyn
                if zx * zx + zy * zy > 4:
                    break
                i += 1
            # 设置像素点的颜色
            r, g, b = i % 4 * 64, i % 8 * 32, i % 16 * 16
            pixels[x, y] = (r, g, b)
    # 显示并保存图像
    img.show()
    img.save('julia_{}.png'.format(c))

# 绘制Julia分形集
julia((-0.70176, -0.3842))

该代码将生成一个大小为1000x1000的Julia分形集，并将其保存为PNG图像。您可以使用不同的参数来生成不同的Julia分形集。