如何根据缩放程度动态设置Mandelbrot集的迭代次数？

问题背景 #

我用Java实现了Mandelbrot集的缩放动画，当前固定最大迭代次数为1000，虽能保证近距离缩放的清晰度，但运行速度极慢；而远景状态下仅需约100次迭代即可保证画面准确。需要一个近似函数f(x)，根据屏幕宽度（对应缩放程度）动态设置迭代次数，兼顾清晰度与性能。

用户原实现代码：

private double[] target = {-1.256640565451168862869, -0.382386428889165027247};

// Returns an integer RGB value (0xRRGGBB) representing the colour which should be drawn at a certain position on the screen
private int getMandleRGB(int x, int y, int w, int h) {
    Complex c = new Complex();
    c.r = lerp(lerp(-2, target[0], 1-zoom), lerp(1, target[0], 1-zoom), x/(double) w);
    c.i = lerp(lerp(-1.5, target[1], 1-zoom), lerp(1.5, target[1], 1-zoom), y/(double) h);
    Complex z = new Complex();
    z.r = c.r;
    z.i = c.i;

    int i = 0;
    for (i = 0; i < 1000; i++) {
        z = Complex.add(Complex.multiply(z, z), c);

        if (z.i*z.i+z.r*z.r > 4) {
            double t = Math.log(i)/Math.log(1000d);
            return (int) lerp((gradient[0] & 0xff0000) >> 16,
                              (gradient[1] & 0xff0000) >> 16, t)*0x10000
                 + (int) lerp((gradient[0] & 0xff00) >> 8,
                              (gradient[1] & 0xff00) >> 8, t)*0x100
                 + (int) lerp((gradient[0] & 0xff),
                              (gradient[1] & 0xff), t);
        }
    }

    return 0x000000; // black
}

动态迭代次数的近似函数 #

Mandelbrot集的迭代需求与像素对应的复平面范围直接相关：缩放倍数越高（画面越近），像素对应复平面范围越小，需要的迭代次数越多；反之则越少。推荐以下两种实用方案：

方案1：基于像素尺寸的对数函数 #

通过计算单个像素对应的复平面宽度，用对数关系映射到迭代次数：

// 计算当前屏幕单个像素对应的复平面范围
double realMin = lerp(-2, target[0], 1-zoom);
double realMax = lerp(1, target[0], 1-zoom);
double pixelSize = (realMax - realMin) / w;

// 动态计算迭代次数，限制在100-1000区间，系数可根据测试调整
int maxIterations = (int) Math.max(100, Math.min(1000, 600 * Math.log(1 / pixelSize)));

逻辑：像素尺寸越小（放大越近），1/pixelSize越大，对数结果越高，迭代次数越多；远景时像素尺寸大，迭代次数自动降到100。

方案2：基于缩放因子的对数映射 #

如果直接用zoom变量（值越大放大倍数越高），可以用：

// 假设zoom=0对应远景，zoom趋近1对应最大放大
double zoomLevel = 1 / (1 - zoom);
// 把缩放级别映射到100-1000的迭代区间
int maxIterations = (int) Math.max(100, Math.min(1000, 100 + 900 * Math.log10(zoomLevel) / Math.log10(1000)));

这里假设最大缩放倍数为1000倍（对应迭代次数1000），远景为1倍（对应100次），通过对数映射实现平滑过渡。

修改后的完整代码 #

将固定迭代次数替换为动态计算值，同时修正颜色渐变的基准值：

private double[] target = {-1.256640565451168862869, -0.382386428889165027247};

// Returns an integer RGB value (0xRRGGBB) representing the colour which should be drawn at a certain position on the screen
private int getMandleRGB(int x, int y, int w, int h) {
    // 计算像素对应的复平面尺寸
    double realMin = lerp(-2, target[0], 1-zoom);
    double realMax = lerp(1, target[0], 1-zoom);
    double pixelSize = (realMax - realMin) / w;
    
    // 动态计算最大迭代次数
    int maxIterations = (int) Math.max(100, Math.min(1000, 600 * Math.log(1 / pixelSize)));

    Complex c = new Complex();
    c.r = lerp(realMin, realMax, x/(double) w);
    c.i = lerp(lerp(-1.5, target[1], 1-zoom), lerp(1.5, target[1], 1-zoom), y/(double) h);
    Complex z = new Complex();
    z.r = c.r;
    z.i = c.i;

    int i = 0;
    for (i = 0; i < maxIterations; i++) {
        z = Complex.add(Complex.multiply(z, z), c);

        if (z.i*z.i+z.r*z.r > 4) {
            // 用当前maxIterations替代固定的1000，保证颜色渐变连续
            double t = Math.log(i)/Math.log(maxIterations);
            return (int) lerp((gradient[0] & 0xff0000) >> 16,
                              (gradient[1] & 0xff0000) >> 16, t)*0x10000
                 + (int) lerp((gradient[0] & 0xff00) >> 8,
                              (gradient[1] & 0xff00) >> 8, t)*0x100
                 + (int) lerp((gradient[0] & 0xff),
                              (gradient[1] & 0xff), t);
        }
    }

    return 0x000000; // black
}

调优提示 #

• 可以调整对数函数的系数（示例中的600），根据实际画面效果和运行速度找到平衡点。
• 若动画卡顿，可在缩放过程中临时降低迭代次数，缩放完成后再补全高迭代次数的画面（渐进式渲染）。
• 极端近距离缩放时，可适当放宽迭代次数上限，但需注意性能开销。

内容的提问来源于stack exchange，提问作者Matthew Miles