1.4.4 分形_Mandelbrot集_C程序.doc

1.4.4 Mandelbrot集的C语言程序设计 （王积社） 作为参考，下面给出一个较为简单的用C语言设计的绘制Mandelbrot集的程序． 回忆Mandelbrot集的迭代关系： 在 中固定初值Ｚ0，取使迭代序列(1)有界的参数 c， 记 则(1)变为 只要实现（２）即可 　　由于算法较简单，故未多言，下面直接给出程序，且对程序作出相应的说明． #include graphics.h /*图形屏幕函数头文件*/ #include stdlib.h /*其它函数头文件*/ #include math.h /*数学函数头文件*/ int main() { float m,dx,dy,x,y,x_n,y_n,Cx,Cy; int n,i,j,L=4; int gdriver=DETECT,gmode; /* DETECT；用于硬件测试. 本句是测试显示器*/ initgraph(gdriver,gmode, ); /* initgraph()：初始化图形系统；头文件：graphics.h；*/ /*参数gdriver: 图形驱动器；gmode： 图形模式的；*/ /* ：指定图形文件的路径，为空时意为默认路径*/ setbkcolor(0); /* setbkcolor(1)： 设置背景颜色；头文件：graphics.h*/ /*参数1表示背景颜色为：蓝色*/ dx=4.0/800; /*行微元，可调整*/ dy=3.0/600; /*列微元，可调整*/ for(i=0;i800;i++) { Cx= -1.9+i*dx; /*计算常数p*/ for(j=0;j600;j++) { Cy= -1.2+j*dy; /*计算常数q*/ x=y=0; /*取Z0=0*/ for(n=0;n=1000;n++) { x_n=x*x-y*y+Cx; /*计算Xn*/ y_n=2*x*y+Cy; /*计算Yn*/ m=x_n*x_n; /*计算强行退出标志“M”，此处M=,当然也可换为其它*/ if(mL) break; x=x_n; y=y_n; } putpixel(i,j,(int)(0.4*m)%16); /*在屏幕上画一个点，原型是：void far putpixel(int x, int y, int color) */ } } getch(); closegraph(); /*关闭图形工作方式*/ } 可以改变程序中的有关数据，观察图形变化． 1