迭代次数较少时Mandelbrot Set绘制轮廓异常的问题排查

你的代码存在两个关键问题，其中逃逸判断逻辑错误是导致低迭代次数下出现云状异常的核心原因：

1. 错误的逃逸检测条件 #

Mandelbrot集的逃逸规则是：当复数z的模长平方（即z_real² + z_imag²）大于4时，该点会逃逸到无穷大，可提前终止迭代。但你的代码中单独判断实部或虚部是否超过4：

if(z[0] > 4 || -4 > z[0] || -4 > z[1] || 4 < z[1])
    break;

这个条件会漏掉大量模长平方已经超过4，但单个分量未达到4的情况（比如z=(3,1)，模长平方为10>4，但实部3<4，虚部1<4）。低迭代次数时，这类未被正确检测的逃逸点会继续迭代，导致n的计数错误，最终绘制出的颜色与预期不符，形成云状轮廓。

2. 迭代顺序的小瑕疵（不影响核心，但不符合常规逻辑） #

你的代码在每次循环中先判断当前z是否逃逸，再计算新的z。常规的Mandelbrot迭代流程是先计算新的z，再判断是否逃逸（初始z为0，必然在圆内，第一次判断无意义）。不过这个问题在高迭代次数下会被掩盖，但低迭代时可能加剧误差。

修正后的代码 #

把逃逸条件改为判断模长平方，同时调整迭代顺序让逻辑更清晰：

double getFracPoint(double x,double y){
    // 坐标缩放逻辑保留
    x = x * ((plotEnd.x-plotStart.x) / SCREENWIDTH) + plotStart.x;
    y = y * ((plotEnd.y-plotStart.y) / SCREENHEIGHT) + plotStart.y;
    x /= zoom;
    y /= zoom;

    double z[2] = {0.0, 0.0};
    double c[2] = {x, y};
    int n = 0;

    for(int i = 0; i < ITERS; i++){
        double z_sq = z[0]*z[0] + z[1]*z[1];
        // 正确的逃逸判断：模长平方大于4
        if(z_sq > 4.0){
            break;
        }
        // 计算z² + c
        double t = z[1];
        z[1] = 2 * z[0] * z[1] + c[1];
        z[0] = z[0]*z[0] - t*t + c[0];
        n++;
    }

    return (double)n/(double)ITERS;
}

说明 #

• 模长平方的计算避免了开根号操作，既符合数学逻辑又保证了性能；
• 调整后的迭代顺序先计算新z再判断（初始z为0时跳过无意义的判断），让逻辑更贴合Mandelbrot集的迭代定义；
• 修正后，低迭代次数下的逃逸点会被正确检测，绘制出的图形会符合预期的轮廓样式。

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