1.4.4 分形_mandelbrot集_c程序

《1.4.4 分形_mandelbrot集_c程序》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、1.4.4Mandelbrot集的C语言程序设计（王积社）作为参考，下面给出一个较为简单的用C语言设计的绘制Mandelbrot集的程序．回忆Mandelbrot集的迭代关系：在中固定初值Ｚ0，取使迭代序列(1)有界的参数c，记则(1)变为只要实现（２）即可　　由于算法较简单，故未多言，下面直接给出程序，且对程序作出相应的说明．#include"graphics.h"/*图形屏幕函数头文件*/#include"stdlib.h"/*其它函数头文件*/#include"math.h"/*数学函数头文件*/intmain()3{floatm,dx,dy,x,y,x_n,y_n,C

2、x,Cy;intn,i,j,L=4;intgdriver=DETECT,gmode;/*DETECT；用于硬件测试.本句是测试显示器*/initgraph(&gdriver,&gmode,"");/*initgraph()：初始化图形系统；头文件：；*//*参数gdriver:图形驱动器；gmode：图形模式的；*//*""：指定图形文件的路径，为空时意为默认路径*/setbkcolor(0);/*setbkcolor(1)：设置背景颜色；头文件：*//*参数1表示背景颜色为：蓝色*/dx=4.0/800;/*行微元，可调整*/

3、dy=3.0/600;/*列微元，可调整*/for(i=0;i<800;i++){Cx=-1.9+i*dx;/*计算常数p*/for(j=0;j<600;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*/3m=x_n*x_n;/*计算强行退出标志“M”，此处M=,当然也可换为其它*/if(m>L)break;x=x_n;y=y_n;}putpixel(i,j,(int)(0.4*m)%16);/*在屏幕上画一个点，

4、原型是：voidfarputpixel(intx,inty,intcolor)*/}}getch();closegraph();/*关闭图形工作方式*/}可以改变程序中的有关数据，观察图形变化．3