计算机图形学实验三：圆的生成算法的实现

实验三：圆的生成算法的实现班级08信计学号20080502069姓名分数一、实验目的和要求1、理解圆的基本原理。2、掌握几种常见的圆生成算法。3、利用TurboC实现圆生成的中点画圆算法。4、理解圆生成的基本原理，掌握几种常见的圆生成算法。5、利用VisualC++实现圆生成的中点画圆的算法。6、利用VisualC++实现圆的Bresenham算法。7、简单了解其他算法。二、实验内容：1.利用中点画图算法，在屏幕上生成任意一段圆弧。2.利用图的对称性，将(1)题生成的圆弧扩展为一个整圆。3.利用bresebham算法设计出一段圆弧。三、实验步骤：1.预习教材关于圆的生成原理。2.仿照教材关于圆生成的中点画圆算法和bresenham算法，使用C++实现该算法。3.调试、编译、运行程序。四、实验结果分析：1.这里用圆生成的中点画圆算法为例来绘制一个圆。程序代码：#include<graphics.h>#include<conio.h>voidcirclepoints(intx,inty,intcolor){ intm,n; floatxasp,yasp; floataspectratio; m=200; n=200; getaspectratio(&xasp,&yasp); aspectratio=xasp/yasp; putpixel(x+m,y*aspectratio+n,color); putpixel(y+m,x*aspectratio+n,color); putpixel(-y+m,x*aspectratio+n,color); putpixel(-x+m,y*aspectratio+n,color); putpixel(y+m,-x*aspectratio+n,color);putpixel(x+m,-y*aspectratio+n,color);putpixel(-x+m,-y*aspectratio+n,color); putpixel(-y+m,-x*aspectratio+n,color);}voidmidpointcircle(intr,intc){ intx,y; floatd; x=0;y=r;d=5.0/4-r; circlepoints(x,y,c); while(y>x) { if(d<=0) d+=2.0*x+3; else { d+=2.0*(x-y)+5; y--; } x++; circlepoints(x,y,c); }}voidmain(){ inta,b; intgraphdriver=DETECT; intgraphmode=0; initgraph(&graphdriver,&graphmode,""); cleardevice(); a=200;b=200; midpointcircle(a,b); getch(); closegraph();}运行结果:2.利用bresenham算法生成圆的代码：#include<graphics.h>#include<stdio.h>#include<conio.h>voidBresenhemCircle(intcenterx,intcentery,intradius,intcolor,inttype);voidmain(){intdrive=DETECT,mode;inti,j;initgraph(&drive,&mode,"");BresenhemCircle(300,200,100,150,0);getch();closegraph();}voidBresenhemCircle(intcenterx,intcentery,intradius,intcolor,inttype){intx=type=0;/*初始横坐标为原点*/inty=radius;/*初始纵坐标远离原点*/intdelta=2*(1-radius);intdirection;while(y>=0){getch();if(!type)/*执行*/{/*在上半圆画两点*/putpixel(centerx+x,centery+y,color);putpixel(centerx-x,centery+y,color);/*在下半圆画两点*/putpixel(centerx-x,centery-y,color);putpixel(centerx+x,centery-y,color);getch();}else{line(centerx+x,centery+y,centerx+x,centery-y);line(centerx-x,centery+y,centerx-x,centery-y);getch();}if(delta<0){if((2*(delta+y)-1)<0)direction=1;elsedirection=2;}elseif(delta>0){if((2*(delta-x)-1)>0)direction=3;elsedirection=2;}elsedirection=2;switch(direction){case1:x++;delta+=(2*x+1);break;case2:x++;y--;delt