好吧,这里有一个奇怪的,我有问题(与gcc btw编译)
下面是用于命令提示符的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");
}
}
}
谢谢你走了这么远,感谢你的帮助!
你只是递归地平方pos
,直到它的大小超过极限。这不会产生分形;它将产生一个单位圆。
您需要在每次迭代后将(x,y)坐标添加到平方值。看到维基百科。
EDIT:一些小改动,瞧!
您的escapedvalue太低了,应该是4.00f。