计算机图形学--分形图的生成.docVIP

实验六分形图的生成

班级08信计2班学号20080502056姓名高金玲分数

一、实验目的和要求：

1、掌握分形图的根本原理，熟练其真正的数学内容，解释图形生成的过程。

2、通过对茱莉亚集、谢宾斯基三角形和贝塞尔集的设计，熟练掌握分形图的建立生成。

3、利用茱莉亚集、谢宾斯基三角形和贝塞尔集生成的图形直观了解分形图的内容与过程。

二、实验内容：

1、编程实现茱莉亚集和谢宾斯基三角形的生成图形算法，生成整体有限的图形。

2、进一步编写贝塞尔集生成图形的算法，将其功能放大。

3、将茱莉亚集、谢宾斯基三角形和贝塞尔集的生成图形在应用程序中实现具体的操作要求，观看其图形的结构。

三、程序执行与运行结果

1.茱莉亚集的生成图形程序代码：

#includegraphics.h

#includeconio.h

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

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

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

//定义复数

structCOMPLEX

{

doublere;

doubleim;

};

//定义复数“乘”运算

COMPLEXoperator*(COMPLEXa,COMPLEXb)

{

COMPLEXc;

c.re=a.re*b.re-a.im*b.im;

c.im=a.im*b.re+a.re*b.im;

returnc;

}

//定义复数“加”运算

COMPLEXoperator+(COMPLEXa,COMPLEXb)

{

COMPLEXc;

c.re=a.re+b.re;

c.im=a.im+b.im;

returnc;

}

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

//主函数

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

voidmain()

{

//初始化绘图窗口

initgraph(640,480);

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

//绘制JuliaSet(茱莉亚集)

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

COMPLEXz,c;

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

intx,y,k; //定义循环变量

for(x=0;x640;x++)

{

for(y=0;y480;y++)

{

z.re=-1.6+3.2*(x/640.0);

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

for(k=0;k180;k++)

{

if(z.re*z.re+z.im*z.im4.0) break;

z=z*z+c;

}

putpixel(x,y,(k=180)?0:HSLtoRGB((float)((k5)%360),1.0,0.5));

}

}

//按任意键退出

getch();

closegraph();

}

运行结果：

茱莉亚集的生成图形

2.贝塞尔集的生成图形代码

#includegraphics.h

#includeconio.h

#defineWIDTH640

#defineHEIGHT480

#defineNUMPTS6

HDChdc;

floatanimpts[NUMPTS*2];

floatdeltas[NUMPTS*2];

voidInit()

{

for(inti=0;iNUMPTS*2;i+=2)

{

animpts[i]=(float)rand()/RAND_MAX*WIDTH;

animpts[i+1]=(float)rand()/RAND_MAX*HEIGHT;

deltas[i]=(float)rand()/RAND_MAX*4+2;

deltas[i+1]=(float)rand()/RAND_MAX*4+2;

i