计算机图形学--分形图的生成算法的实现

实验六：分形图的生成算法的实现

班级 08信计2班学号 20080502063 姓名李宁分数

一、实验目的与要求：

1、了解WIN-TC1.72等编程环境中常用控件命令与绘图函数，初步掌握在实验设

计集成环境下进行图形处理程序的设计方法。

2、熟练掌握分形图的算法：分形学- Mandelbrot Set (曼德布洛特集)，分形学-

Julia Set (茱莉亚集)。

二、实验内容：

1、在WIN-TC1.72环境中设计程序，利用消息处理函数，搭建能运行图形算法程

序的平台。

2、根据教材中给定的算法，实现分形图生成算法：分形学- Mandelbrot Set (曼

德布洛特集)，分形学- Julia Set (茱莉亚集)

三、实验结果分析：

1、实验程序

分形学- Mandelbrot Set (曼德布洛特集)程序实现的完整源程序

#include

#include

/////////////////////////////////////////////////

// 定义复数及乘、加运算

/////////////////////////////////////////////////

// 定义复数

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;

}

/////////////////////////////////////////////////

// 主函数

/////////////////////////////////////////////////

void main()

{

// 初始化绘图窗口

initgraph(640, 480);

/////////////////////////////////////////////////

// 绘制Mandelbrot Set (曼德布洛特集)

/////////////////////////////////////////////////

COMPLEX z, c;

int x, y, k; // 定义循环变量

for(x = 0; x < 640; x++)

{

c.re = -2.1 + (1.1 - -2.1) * (x / 640.0);

for(y = 0; y < 480; y++)

{

c.im = -1.2 + (1.2 - -1.2) * (y / 480.0);

z.re = z.im = 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();

}

分形学- Julia Set (茱莉亚集)

#include

#include

/////////////////////////////////////////////////

// 定义复数及乘、加运算

/////////////////////////////////////////////////

// 定义复数

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;

}

/////////////////////////////////////////////////

// 主函数

/////////////////////////////////////////////////

void main()

{

// 初始化绘图窗口

initgraph(640, 480);

/////////////////////////////////////////////////

// 绘制Julia Set (茱莉亚集)

/////////////////////////////////////////////////

COMPLEX z, c;

c.re = -0.75, c.im = 0; // 设置迭代初值