低分辨率Mandelbrot分形不是......分辨率足够高?

时间:2012-12-25 09:50:28

标签: c fractals mandelbrot

好的,这是一个奇怪的问题,我遇到了问题(用gcc btw编译)

以下是用于命令提示符的Mandelbrot分形生成器的源代码。我以前做过这个,我想加快测试自己,看看我能在命令提示符下生成实际生成Mandelbrot分形所需的代码有多快。我经常这样做,以便有趣地测试自己

无论如何我遇到了一个新问题,我无法弄清问题是什么。当分形渲染无论迭代次数或我设置的escapeValue时,它总是显示为椭圆形!它不应该这样做。

对于你们所有的mandelbrot / cpp极客们,你们可以帮助我找出为什么我没有得到更多的“mandelbrot”造型吗?

#include <stdio.h>
#include <math.h>

#define DOSWidth 80
#define DOSHeight 25

int iterations = 1024;
float escapeValue = 3.0f;

struct ivar {
    ivar(float _x, float _i) {
        x = _x;
        i = _i;
    }
    void log() {printf("(%g%c%gi)", x, (i<0)?'-':'+', fabs(i));}
    float magnitude() {return sqrtf(x*x+i*i);}
    ivar square() {return ivar(x, i)*ivar(x, i);}

    ivar operator + (ivar v) {return ivar(x+v.x, i+v.i);};
    ivar operator - (ivar v) {return ivar(x-v.x, i-v.i);};
    ivar operator * (ivar v) {return ivar(x*v.x-(i*v.i), x*v.i+i*v.x);};

    float x, i;
};

struct rect {
    rect(float _x, float _y, float _width, float _height) {
        x = _x;y = _y;width = _width;height = _height;
    }

    void setCenter(float cx, float cy) {
        x = cx-width/2.0f;
        y = cy-width/2.0f;
    }

    void log() {printf("(%f, %f, %f, %f)", x, y, width, height);}

    float x, y;
    float width, height;
};

int main() {
    rect region = rect(0, 0, 2.5f, 2.0f);
    region.setCenter(0, 0);
    float xSize = region.width / (float)DOSWidth;
    float ySize = region.height / (float)DOSHeight;
    for(int y=0;y<DOSHeight;y++) {
        for(int x=0;x<DOSWidth;x++) {
            ivar pos = ivar(x*xSize+region.x, y*ySize+region.y);
            bool escapes = false;
            for(int i=0;i<iterations;i++) {
                if(pos.magnitude() > escapeValue) {
                    escapes = true;
                    break;
                }
                pos = pos.square();
            }
            if(escapes)printf(" ");
            else printf("X");
        }
    }
}

如果你有这么远的话,谢谢,谢谢你的帮助!

2 个答案:

答案 0 :(得分:3)

您只是递归地对pos进行平方,直到其幅度超出限制。这不会产生分形;它会产生一个单位圆。

每次迭代后,您需要将(x,y)坐标添加到平方值。请参阅Wikipedia

编辑:一些小改动和voila

答案 1 :(得分:0)

您的转义价值太低,应该是4.00f。