VC分形学-可以无穷放大的MandelbrotSet(曼德布洛特集).docVIP

// 程序名称：分形学 - 可以无穷放大的 Mandelbrot Set 曼德布洛特集 // 编译环境：Visual C++ 6.0，EasyX 2011惊蛰版 // 最后更新：2010-9-9 // #include #include // 定义常量 #define ITERATIONS 1000 // 迭代次数，越高，图像越精细 #define MAXCOLOR 64 // 颜色数 ///////////////////////////////////////////////// // 定义复数及乘、加运算 ///////////////////////////////////////////////// // 定义复数 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; ///////////////////////////////////////////////// // 定义颜色及初始化颜色 ///////////////////////////////////////////////// // 定义颜色 int Color[MAXCOLOR]; // 初始化颜色 void InitColor // 使用 HSL 颜色模式产生角度 h1 到 h2 的渐变色 int h1 240, h2 30; for int i 0; i MAXCOLOR/2; i++ Color[i] HSLtoRGB float h1, 1.0f, i * 2.0f / MAXCOLOR ; Color[MAXCOLOR-1-i] HSLtoRGB float h2, 1.0f, i * 2.0f / MAXCOLOR ; ///////////////////////////////////////////////// // 绘制 Mandelbrot Set 曼德布洛特集 ///////////////////////////////////////////////// void Draw double fromx, double fromy, double tox, double toy COMPLEX z, c; int x, y, k; // 定义循环变量 for x 0; x 640; x++ c.re fromx + tox - fromx * x / 640.0 ; for y 0; y 480; y++ c.im fromy + toy - fromy * y / 480.0 ; z.re z.im 0; for k 0; k ITERATIONS; k++ if z.re * z.re + z.im * z.im 4.0 break; z z * z + c; putpixel x, y, k ITERATIONS ? 0 : Color[k % MAXCOLOR] ; ///////////////////////////////////////////////// // 主函数 ///////////////////////////////////////////////// void main // 初始化绘图窗口及颜色 initgraph 640, 480 ; InitColor ; // 初始化 Mandelbrot Set 曼德布洛特集 坐标系 double fromx, fromy, tox, toy; fromx -2.1; tox 1.1; fromy -1.2; toy 1.2; Draw fromx, fromy, tox, toy ; // 捕获鼠标操作，实现放大鼠标选中区域 MOUSEMSG m; bool isLDown false; int selfx, selfy, seltx, selty; // 定义选区 while !kbhit m GetMouseMsg