计算机图形学 分形图的生成.doc

实验六 分形图的生成班级 08信计2班 学号 20080502084 姓名 吴涛 分数 一、实验目的和要求：1、掌握希尔宾斯基三角形和Julia Set (茱莉亚集)的基本原理。2、熟悉两个图形的生成算法。3、掌握希尔宾斯基三角形和和Julia Set (茱莉亚集)的绘制.。4、提高分形图形生成的理解应用能力。二、实验内容：1、对于第一个图形在平面内随机的设置种子，并由此而设定三角形的三个顶点。形成初始化模式后，绘制三万个点，使规则传递下去。2、对于第二个图形则运用逃逸时间法后设定一个常数c的值。 3、对两图形分别进行分析对比其局部与整体的自相似性。 三、实验结果分析：1、该程序实现了递归算法和逃逸时间法的图形绘制 2、比较每一小部分与整体的关系： 图形的层次是无限的、分形往往可以从局部“看出”整体、虽然看上去十分复杂，但其背后的规则却是相当简单。四、程序代码：1、Julia Set (茱莉亚集)的程序代码为：#include #include / 定义复数及乘、加运算/ 定义复数struct COMPLEXdouble 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);/ 绘制 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 4.0 )break;z = z * z + c;putpixel(x, y, (k = 180) ? 0 : HSLtoRGB(float)(k 5) % 360), 1.0, 0.5);/ 按任意键退出getch();closegraph();输出的结果为： 2、希尔宾斯基三角形为：#include #include #include void main()srand(unsigned)time(NULL);/ 设置随机种子POINT P3 = 320, 50, 120, 400, 520, 400;/ 设定三角形的三个顶点POINT p = rand() % 640, rand() % 480;/ 随机产生当前点/ 初始化图形模式initgraph(640, 480);/ 绘制三万个点int n;for(int i = 0; i = 30000; i+)n = rand() % 3;p.x = (p.x + Pn.x) / 2;p.y = (p.y + Pn.y) / 2;putpixel(p.x, p.y, GREEN);/ 按任意键退出getch();closegraph();运行结果为：3、自相似图形程序：#include #include #include #include void star1(int x,int y,int r);void star2(int x, int y , int r);void main() int graphdriver,graphmode,x,y,r; graphdriver=DETECT; initgraph(&graphdriver,&graphmode, ); setcolor(LIGHTMAGENTA); star1(160,160,80); getch(); star2(480,320,80);getch(); void star1(int x,int y,int r) if(r0)star1(x-r,y+r,r/2);star1(x+r,y+r,r/2);star1(x-r,y-r,r/2);star1(x+r,y-r,r/2);bar(x-r,y+r,x+r,y-r);void star2(int x, int y , int r) if(r0) star2(x-r