[分形学]JuliaSet(茱莉亚集)VC源代码.docxVIP

关于 Julia Set (茱莉亚集) 的介绍什么的我就不多说了，网上一大堆。执行效果如图：关于 Julia Set，可以通过设置复数 c 的初值，显示出不同的图案，比如，大家可以试试以下几组：c.re = 0.45, c.im = -0.1428;c.re = 0.285, c.im = 0.01;c.re = 0.285, c.im = 0;c.re = -0.8, c.im = 0.156;c.re = -0.835, c.im = -0.2321;c.re = -0.70176, c.im = -0.3842;随便用哪行替换掉源程序中的“c.re = -0.75, c.im = 0;”就可以看到不同的图案了。循环变量 k 是迭代次数，在某些参数下需要高一些会更精细。为了美观，还需要修改一下颜色部分，目前代码中的颜色是这样的：HSLtoRGB((float)((k5) % 360), 1.0, 0.5)这行代码中的 k 的取值范围是 0~180，将其映射到 HSL 颜色空间中的色相上(360 度)。全部代码如下：// 程序名称：分形学 - Julia Set (茱莉亚集)// 编译环境：Visual C++ 6.0，EasyX 2011惊蛰版// 最后更新：2010-9-9//#includegraphics.h#includeconio.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;}/////////////////////////////////////////////////// 主函数/////////////////////////////////////////////////voidmain(){// 初始化绘图窗口initgraph(640,480);/////////////////////////////////////////////////// 绘制 Julia Set (茱莉亚集)/////////////////////////////////////////////////COMPLEX z, c;c.re =-0.75, c.im =0;// 设置迭代初值int x, y, k;// 定义循环变量for(x =0; x 640; x++){for(y =0; y 480; y++){z.re =-1.6+3.2*(x /640.0);z.im =-1.2+2.4*(y /480.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();}+1