[分形学] Mandelbrot Set (曼德布洛特集) VC 源代码
作者:yangw80

新版地址:http://www.codebus.cn/yw80/post/mandelbrot-set

关于 Mandelbrot Set (曼德布洛特集) 的介绍什么的我就不多说了,网上一大堆。执行效果如图:

为了美观,可以修改一下颜色部分,目前代码中的颜色是这样的:

HSLtoRGB((float)((k<<5) % 360), 1.0, 0.5)

这行代码中的 k 的取值范围是 0~180,将其映射到 HSL 颜色空间中的色相上(360 度)。

全部代码如下:

// 程序名称:分形学 - Mandelbrot Set (曼德布洛特集)
// 编译环境:Visual C++ 6.0,EasyX 2011惊蛰版
// 最后更新:2010-9-9
//
#include <graphics.h>
#include <conio.h>

/////////////////////////////////////////////////
// 定义复数及乘、加运算
/////////////////////////////////////////////////

// 定义复数
struct COMPLEX
{
	double re;
	double im;
};

// 定义复数“乘”运算
COMPLEX operator * (COMPLEX a, COMPLEX b)
{
	COMPLEX c;
	c.re = a.re * b.re - a.im * b.im;
	c.im = a.im * b.re + a.re * b.im;
	return c;
}

// 定义复数“加”运算
COMPLEX operator + (COMPLEX a, COMPLEX b)
{
	COMPLEX c;
	c.re = a.re + b.re;
	c.im = a.im + b.im;
	return c;
}


/////////////////////////////////////////////////
// 主函数
/////////////////////////////////////////////////
void main()
{
	// 初始化绘图窗口
	initgraph(640, 480);

	/////////////////////////////////////////////////
	// 绘制 Mandelbrot Set (曼德布洛特集)
	/////////////////////////////////////////////////
	COMPLEX z, c;
	int x, y, k;	// 定义循环变量
	for(x = 0; x < 640; x++)
	{
		c.re = -2.1 + (1.1 - -2.1) * (x / 640.0);
		for(y = 0; y < 480; y++)
		{
			c.im = -1.2 + (1.2 - -1.2) * (y / 480.0);
			z.re = z.im = 0;
			for(k = 0; k < 180; k++)
			{
				if ( z.re * z.re + z.im * z.im > 4.0 )	break;
				z = z * z + c;
			}
			putpixel(x, y, (k >= 180) ? 0 : HSLtoRGB((float)((k << 5) % 360), 1.0, 0.5));
		}
	}

	// 按任意键退出
	getch();
	closegraph();
}
更新时间:2010/9/9