计算机图形学实验--完整版-带结果--vc++实现

计算机图形学实验报告信息学院计算机专业20081060183 周建明
综括：
利用计算机编程语言绘制图形，主要实现以下内容：
（1）、中点算法生成任意斜率直线，并设置线型线宽。

（2）、中点算法生成圆
（3）、中点算法生成椭圆
（4）、扫描算法实现任意多边形填充
（5）、Cohen_Sutherland裁剪
（6）、自由曲线与曲面的绘制
（7）、二维图形变换
（8）、三视图变换
实验一、直线的生成
一、实验内容
根据提供的程序框架，修改部分代码，完成画一条直线的功能（中点画线法或者Bresenham画线法任选一），只要求实现在第一象限内的直线。

二、算法原理介绍
双击直线生成.dsw打开给定的程序，或者先启动VC++，文件（file）→打开工作空间（open workspace）。

打开直线生成view.cpp，按注释改写下列函数：
1.void CMyView::OnDdaline() （此为DDA生成直线）
2.void CMyView::OnBresenhamline()（此为Bresenham画直线）
3.void CMYView::OnMidPointLine()（此为中点画线法）
三、程序源代码
1.DDA生成直线画法程序：
float x,y,dx,dy,k;
dx=(float)(xb-xa);
dy=(float)(yb-ya);
k=dy/dx;
x=xa;
y=ya;
if(abs(k)<1)
{
for (x=xa;x<=xb;x++)
{
pdc->SetPixel(x, int(y+0.5),COLOR);
y=y+k;
}
}
if(abs(k)>=1)
{
for(y=ya;y<=yb;y++)
{
pdc->SetPixel(int(x+0.5),y,COLOR);
x=x+1/k;
}
}
//DDA画直线结束
}
2.Bresenham画直线源程序：
float b,d,xi,yi;
int i;
float k;
k=(yb-ya)/(xb-xa);
b=(ya*xb-yb*xa)/(xb-xa);
if(k>0&&k<=1)
for(i=0;i<abs(xb-xa);i++)
{ d=ya+0.5-k*(xa+1)-b;
if(d>=0)
{ xi=xa+1;
yi=ya;
xa++;
ya=ya+0.5;
}
if(d<0)
{ xi=xa+1;
yi=ya+1;
xa++;
ya=ya+1.5;
}
pdc->SetPixel(xi,yi,COLOR);
}
//BresenHam画直线结束}
3.中点画线法源程序：
float b,d,xi,yi;
int i;
float k;
k=(yb-ya)/(xb-xa);
b=(ya*xb-yb*xa)/(xb-xa);
if(k>0&&k<=1)
for(i=0;i<abs(xb-xa);i++)
{ d=ya+0.5-k*(xa+1)-b;
if(d>=0)
{ xi=xa+1;
yi=ya;
xa++;
ya=ya+0.5;
}
if(d<0)
{ xi=xa+1;
yi=ya+1;
xa++;
ya=ya+1.5;
}
pdc->SetPixel(xi,yi,COLOR); }
//BresenHam画直线结束}
四、实验结果
1、DDA生成直线
2、Bresenham画直线
3、中点画线法
实验二、bresenham画圆
一、实验内容
根据提供的程序框架，修改部分代码，用Bresenham画法画一段圆弧或者画圆。

二、算法原理介绍
双击experiment.dsw打开，或者先启动VC++，文件（file）→打开工作空间（open workspace）。

打开直线生成view.cpp，按注释改写下列函数：
void CMyView::OnBresenhamCircle()
三、程序源代码
bresenham画圆源程序
void Bresenhan_circle()
;int x,y=180,color=1,xcenter,ycenter;x=0;
int d=3-2*180;
while(x<y)
{
if(d<0){d=d+4*x+6;x=x+1;}
else{d=d+4*(x-y)+10;x=x+1;y=y-1;}
pdc->SetPixel(0+x,0+y,COLOR);
pdc->SetPixel(0+x,0-y,COLOR);
pdc->SetPixel(0-x,0+y,COLOR);
pdc->SetPixel(0-x,0-y,COLOR);
pdc->SetPixel(0+y,0+x,COLOR);
pdc->SetPixel(0+y,0-x,COLOR);
pdc->SetPixel(0-y,0+x,COLOR);
pdc->SetPixel(0-y,0-x,COLOR);
}
四、实验结果
实验三、二维图形的几何变换
一、实验内容
以一条直线段为例，完成目标的平移、绕任一点旋转。

（缩放等其它变换选做）
二、实验步骤和方法
1、打开二维变换.dsw
2、改写二维变换view.cpp里的void CMyView::Onrecycle()函数（需要改写的地
方我已经做了说明）。

3、生成直线的函数采用VC里的函数。

4、函数的主要任务是计算出变换后的坐标。

三、程序源代码
1.目标的平移的源程序
xa1=xa+dx;
ya1=ya+dy;
xb1=xb+dx;
yb1=yb+dy;
2.绕任意点旋转的源程序
xa1=cos(angle/57.1)*xa-sin(angle/57.1)*ya+x-x*cos(angle/57.1)+y*sin(angle/57.1);y a1=sin(angle/57.1)*xa+cos(angle/57.1)*ya+y-y*cos(angle/57.1)-x*sin(angle/57.1); xb1=cos(angle/57.1)*xb-sin(angle/57.1)*yb+x-x*cos(angle/57.1)+y*sin(angle/57.1);y b1=sin(angle/57.1)*xb+cos(angle/57.1)*yb+y-y*cos(angle/57.1)-x*sin(angle/57.1);
四、实验结果
1、目标的平移的源程序
2、绕任意点旋转的源程序
实验四、用扫描线算法实现多边形填充：
（1）算法思想：
用水平扫描线从上到下扫描由点线段构成的多段定义的多边形。

每根扫描线与多边形各边产生一系列交点，将这些交点按照x坐标进行排序，将排序后的交点成对取出，作为两个端点，用所需填充的色彩画水平直线。

多边形被扫描完毕，则填充结束。

实现步骤：
（1）建立边的分类表ET；
（2）将扫描线纵坐标y的初值为ET中非空元素的最小序号；
（3）置活化边表AEL为空；
（4）执行下列步骤直至ET和AEL都为空；
A、如果ET中的第y类非空，则将其中的所有边取出并插入AEL中，在插入过程忠进行排序；
B、对AEL中的边两两配对，将每对边中x坐标按规则取整，获得有效的填充区段，再填充；
C、将当前扫描线纵坐标y值递增1，即y=1；
D、将AEL中满足y=ymax边删去；
E、对AEL中剩下的每一条边的x递增deltax，即x=x+deltax；
（2）程序实现
typedef struct
{
int y_top;
float x_int;
int delta_y;
float x_change_per_scan; }EACH_ENTRY;
EACH_ENTRY sides[MAX_POINT];
int x[MAX_POINT],y[MAX_POINT];
int side_count,first_s,last_s,scan,bottomscan,x_int_count;
void fill_area(int count)
{
sort_on_bigger_y(count); first_s=1; last_s=1;
for(scan=sides[1].y_top;scan>=bottomscan;scan--)
{ update_first_and_last(count,scan);
process_x_intersections(first_s,last_s);
draw_lines(scan,x_int_count,first_s);
update_sides_list(); } }
sort_on_bigger_y(int n)
{ int k,x1,y1;int tem;
side_count=0; y1=y[n];
x1=x[n]; bottomscan=y[n];
for(k=1;k<n+1;k++)
{ if(y1!=y[k])
{ side_count++; if(k!=n) tem=y[k+1];
else tem=y[1];
put_in_sides_list(side_count,x1,y1,x[k],y[k],tem);
}
else { setcolor(13); line(x1,y1,x[k],y[k]); }
if(y[k]<bottomscan) bottomscan=y[k];
y1=y[k];x1=x[k]; }
}
put_in_sides_list(int entry,int x1,int y1,int x2,int y2,int next_y) { int maxy; float x2_temp,x_change_temp;
x_change_temp=(float)(x2-x1)/(float)(y2-y1);
x2_temp=x2; if((y2>y1)&&(y2<next_y))
{ y2--; x2_temp-=x_change_temp; }
else { if((y2<y1)&&(y2>next_y))
{ y2++; x2_temp+=x_change_temp; } }
maxy=(y1>y2)? y1:y2;
while((entry>1)&&(maxy>sides[entry-1].y_top))
{ sides[entry]=sides[entry-1];
entry--; }
sides[entry].y_top=maxy;
sides[entry].delta_y=abs(y2-y1)+1;
if(y1>y2) sides[entry].x_int=x1;
else sides[entry].x_int=x2_temp;
sides[entry].x_change_per_scan=x_change_temp;
}
update_first_and_last(int count,int scan)
{
while((sides[last_s+1].y_top>=scan)&&(last_s<count))
last_s++; while(sides[first_s].delta_y==0)
first_s++; }
process_x_intersections(int first_s,int last_s)
{ int k; x_int_count=0;
for(k=first_s;k<last_s+1;k++)
{ if(sides[k].delta_y>0)
{ x_int_count++;
sort_on_x(k,first_s);
}
}
}
sort_on_x(int entry,int first_s)
{ while((entry>first_s)&&(sides[entry].x_int<sides[entry-1].x_int)) { swap(&sides[entry],&sides[entry-1]); entry--;
}
}
swap(EACH_ENTRY *x,EACH_ENTRY *y)
{ int i_temp;
float f_temp;
i_temp=x->y_top;
x->y_top=y->y_top;
y->y_top=i_temp;
f_temp=x->x_int;
x->x_int=y->x_int;y->x_int=f_temp;
i_temp=x->delta_y;x->delta_y=y->delta_y;
y->delta_y=i_temp;
f_temp=x->x_change_per_scan;
x->x_change_per_scan=y->x_change_per_scan;
y->x_change_per_scan=f_temp;
}
draw_lines(int scan,int x_int_count,int index)
{
int k,x,x1,x2; for(k=1;k<(int)(x_int_count/2+1.5);k++) { while(sides[index].delta_y==0)
index++; x1=(int)(sides[index].x_int+0.5);
index++;
while(sides[index].delta_y==0)
index++;
x2=(int)(sides[index].x_int+0.5);
setcolor(13);
line(x1,scan,x2,scan); index++;
}
}
update_sides_list() {
int k;
int temp;
for(k=first_s;k<last_s+1;k++)
{ if(sides[k].delta_y>0)
{ sides[k].delta_y--;
sides[k].x_int-=sides[k].x_change_per_scan;
}
}
}
实验五
用Cohen-Sutherland裁剪算法实现直线段
裁剪：
5.1编码算法：
5.1.1算法思想：
本算法分为三个步骤：判断线段两端是否都在窗口内，如果是，线段完全可见；否则判断线段是否显然不可见，如果是，裁剪结束。

否则求线段与窗口边延长线的交点，此线段将线段分为两段，其中一段不可见，舍弃。

对余下的线段继续进行递归裁剪。

算法的编码是进行二进制位运算，若顶点在窗口内，则二进制编码为0000；否则进行编码。

5.1.2实现步骤：
（1）这种算法利用编码的方法，延长窗口边线，使得它们把包含未经裁剪图形的窗口平面区域分成九个区域：
（2）每个区域用一个4位编码CtCbCrCl来表示，代码中每一位分别是0或1，是按照窗口边线来确定的，下面给出具体的编码规则，其中最右边的位Cl是第一位，依次Cr第二、Cb第三、Ct第四位。

A、当两端点P1(x1,y1)和P2(x2,y2)在区域0000中，即满足点的裁剪不等式：
B、当两个端点在窗口边线外的同侧位置，则他们的四位代码中，有一相同位，同时位“1”，显然两个端点代码C1和C2按位与运算C1&C2≠0 .由此可检查判断直线在窗口外，应全部舍弃。

C、果直线两端点不满足上述两种情况，不能简单地全部保留或全部舍弃直线时，则需要计算出直线与窗口边线的交点，将直线分段后继续进行检查判断。

这样可以逐段地舍弃位于窗口外地线段，保留剩余在窗口内的线段。

5.1.3程序实现
/*********定义裁剪函数************/
void line_clip(int x1, int y1, int x2, int y2, int left,int top,int right, int bottom)
{
void getcode(int x,int y,int d[4]); /*定义获得端点的代码函数*/
int i,x11,y11; /*定义交点坐标*/
int aa=1;
getcode(x1,y1,a);
getcode(x2,y2,b);
/*1：裁剪循环开始*/
while(aa!=0)
{
if((a[0]+a[1]+a[2]+a[3]==0)&&(b[0]+b[1]+b[2]+b[3]==0)) /*第一种情况线段完全可见*/
{
aa=0;
return;
}
else if((a[0]&&b[0])+(a[1]&&b[1])+(a[2]&&b[2])+(a[3]&&b[3])!=0) /*线段完全不可见*/ {
setcolor(0);
setwritemode(0); /*设置画线的输出模式为覆盖方式*/
line(x1,y1,x2,y2); /*进行裁剪也就是进行覆盖*/
aa=0;
return;
}
/*2：线段即不完全可见，也不完全不可见，即与边有交点的线段的裁剪处理*/
else
{
if(a[0]+a[1]+a[2]+a[3]==0) /*寻找不可见点*/
{
for(i=0;i<4;i++)
c[i]=b[i];/*x2,y2为不可见点*/
}
else
{
for(i=0;i<4;i++)
c[i]=a[i];/*x1,y1为不可见点*/
}
/*3：直线与窗口边作求交运算，求出交点赋值给x11、y11*/
if(c[0]==1)
{
x11=left;
y11=(int)((y2-y1)*(left-x1)/(x2-x1)+y1);
}
else if(c[1]==1)
{
x11=right;
y11=(int)((y2-y1)*(right-x1)/(x2-x1)+y1);
}
else if(c[2]==1)
{
x11=(int)((x2-x1)*(bottom-y1)/(y2-y1)+x1);
y11=bottom;
}
else if(c[3]==1)
{
x11=(int)((x2-x1)*(top-y1)/(y2-y1)+x1);
y11=top;
} /*3：求交运算结束*/
if((c[0]==a[0])&&(c[1]==a[1])&&(c[2]==a[2])&&(c[3]==a[3]))
{
setcolor(0);
setlinestyle(0,0,3);
setwritemode(0);
line(x1,y1,x11,y11);
x1=x11; y1=y11;
getcode(x11,y11,a);
}
else if((c[0]==b[0])&&(c[1]==b[1])&&(c[2]==b[2])&&(c[3]==b[3])) {
setcolor(0);
setlinestyle(0,0,3);
setwritemode(0);
line(x2,y2,x11,y11);
x2=x11;y2=y11;
getcode(x11,y11,b);
}
} /*2：与边有交点的线段裁剪处理结束*/
} /*1：裁剪循环结束*/
}
/*获得端点代码*/
void getcode(int x,int y,int d[4])
{
d[0]=0;d[1]=0;d[2]=0;d[3]=0;
if(x<left) d[0]=1;
if(x>right) d[1]=1;
if(y>bottom) d[2]=1;
if(y<top) d[3]=1;
return;
}
实验六自由曲线与曲面的绘制
6.1、生成Bezier曲线
（1）算法思想
利用Bezier曲线造型时，如果其次数太高，固然能表示复杂的形状，但同时造成计算复杂度增加，而二次Bezier曲线表示能力有限，故下面我们仅讨论三次Bezier曲线的矩阵表示和计算机生成。

(1) 矩阵表示
三次Bezier曲线的定义式为：
将其分解为两个矢量的点积为：
取Bezier曲线的几何矩阵，则有：
得到三次Bezier曲线的基矩阵为：
(2) 三次Bezier曲线的生成
Bezier曲线用的参数表示法，因此将区间[0，1]分成N等份，确定了控制点后，每取一个都对应一个三维的点。

这样，在每两个和之间画一条直线，再将该这些直线段投影到二维平面，就可以得到我们需要的Bezier曲线了。

（2）程序实现：
void bezier_3(int color, double p[4][2])
{
double t,t1,t2,xt,yt;
int rate=200,x,y;
setcolor(color);
moveto(p[0][0],p[0][1]);
for (t=0;t<=1;t+=1.0/rate)
{
yt=1-t;
t1=yt*yt;
t2=3*yt*t;
xt=p[0][0]*t1*yt+p[1][0]*t2*yt+p[2][0]*t2*t+p[3][0]*t*t*t;
yt=p[0][1]*yt*t1+p[1][1]*t2*yt+p[2][1]*t2*t+p[3][1]*t*t*t;
x=(int)(xt);
y=(int)(yt);
lineto(x,y);
}
}
6.2绘制Bezier曲面
（1）算法思想：
Bezier曲线是一条与控制多边形顶点位置有严格关系的曲线，Bezier曲线形状趋向于特征多边形的形状，而且阶数由控制多边形顶点的个数决定。

Bezier曲面则是由Bezier曲线拓广而来，它也是以Bernstein函数作为基函数，是由Bernstein基函数构造空间点阵列的位置来控制的。

（2）算法实现步骤：
1.给出控制顶点。

2.根据如下公式，编程绘制Bezier曲面。

计算各点坐标
3. 用画直线的方法首尾相接的画直线。

（3）、程序实现：
main()
{
char str[80];
int driver,mode;
int i,k,m,j;
int x0=300,y0=300;
int p[4][4][3]={{{12,145,50},{24,120,50},{45,117,50},{65,140,50}},
{{24,65,34},{56,45,34},{87,56,34},{99,75,34}},
{{43,120,20},{65,98,20},{87,100,20},{111,115,20}},
{{15,160,0},{34,120,0},{55,135,0},{75,170,0}}};
float u0,u1,u2,u3,v0,v1,v2,v3;
float q10,q11,q12,q13,du,dv;
float q[n+1][n+1][3],q20[3],q21[3],q22[3],q23[3],q30[3],q31[3],q32[3],q33[3];
float x,y;
printf("input j:(0---bezier, 1---B_spline)\n");
scanf("%d",&j);
driver=DETECT;mode=0;
initgraph(&driver,&mode,"D:\\TC");
setbkcolor(1);
setcolor(15);
for(i=0;i<=3;i++)
{for(k=0;k<=3;k++)
{x=(-0.7071*p[i][k][0]+0.7071*p[i][k][1]+200)*2-300;
y=(0.4082*p[i][k][0]+0.4082*p[i][k][1]-0.8165*p[i][k][2]+200)*2-300;
if(k==0)moveto(x,y);lineto(x,y);
}
}
for(i=0;i<=3;i++)
{for(k=0;k<=3;k++)
{x=(-0.7071*p[k][i][0]+0.7071*p[k][i][1]+200)*2-300;
y=(0.4082*p[k][i][0]+0.4082*p[k][i][1]-0.8165*p[k][i][2]+200)*2-300;
if(k==0)moveto(x,y);lineto(x,y);
}
}
du=dv=1.0/n;
for(i=0;i<=n;i++)
{
u0=(i*du)*(i*du)*(i*du);u1=(i*du)*(i*du);u2=i*du;u3=1;
if(j==0)
{q10=u0*(-1)+u1*3+u2*(-3)+u3*1;q11=u0*3+u1*(-6)+u2*3;
q12=u0*(-3)+u1*3;q13=u0*1;}
if(j==1)
{q10=(u0*(-1)+u1*3+u2*(-3)+u3*1)/6; q11=(u0*3+u1*(-6)+u3*4)/6;
q12=(u0*(-3)+u1*3+u2*3+u3*1)/6;q13=(u0*1)/6;}
for(k=0;k<=2;k++)
{
q20[k]=q10*p[0][0][k]+q11*p[1][0][k]+q12*p[2][0][k]+q13*p[3][0][k]; q21[k]=q10*p[0][1][k]+q11*p[1][1][k]+q12*p[2][1][k]+q13*p[3][1][k]; q22[k]=q10*p[0][2][k]+q11*p[1][2][k]+q12*p[2][2][k]+q13*p[3][2][k]; q23[k]=q10*p[0][3][k]+q11*p[1][3][k]+q12*p[2][3][k]+q13*p[3][3][k]; if(j==0)
{q30[k]=q20[k]*(-1)+q21[k]*3+q22[k]*(-3)+q23[k]*1;
q31[k]=q20[k]*3+q21[k]*(-6)+q22[k]*3;
q32[k]=q20[k]*(-3)+q21[k]*3;
q33[k]=q20[k]*1;}
if(j==1)
{q30[k]=(q20[k]*(-1)+q21[k]*3+q22[k]*(-3)+q23[k]*1)/6;
q31[k]=(q20[k]*3+q21[k]*(-6)+q22[k]*3)/6;
q32[k]=(q20[k]*(-3)+q22[k]*3)/6;
q33[k]=(q20[k]*1+q21[k]*4+q22[k]*1)/6;}
}
for(m=0;m<=n;m++)
{
v0=(m*dv)*(m*dv)*(m*dv);v1=(m*dv)*(m*dv);v2=m*dv;v3=1;
for(k=0;k<=2;k++)
{q[i][m][k]=q30[k]*v0+q31[k]*v1+q32[k]*v2+q33[k]*v3;}
}
}
setcolor(4);
for(i=0;i<=n;i++)
{for(m=0;m<=n;m++)
{x=(-0.7071*q[i][m][0]+0.7071*q[i][m][1]+200)*2-300;
y=(0.4082*q[i][m][0]+0.4082*q[i][m][1]-0.8165*q[i][m][2]+200)*2-300;
if(m==0)moveto(x,y);
lineto(x,y);
}
}
for(m=0;m<=n;m++)
{for(i=0;i<=n;i++)
{x=(-0.7071*q[i][m][0]+0.7071*q[i][m][1]+200)*2-300;
y=(0.4082*q[i][m][0]+0.4082*q[i][m][1]-0.8165*q[i][m][2]+200)*2-300;
if(m==0)moveto(x,y);
lineto(x,y);
}
}
if(j==0)
{ printf("3Bezier curved surface");}
if(j==1)
{ printf("3B_Splline curved surface");}
outtextxy(270,270,"this is zhouman's homework");
setcolor(3);
settextstyle(TRIPLEX_FONT,HORIZ_DIR,1);
sprintf(str,"%s","----graphics:zhouman--time:2009.12.30----");
outtextxy(120,430,str);
getch();
closegraph();
}
实验七二维图形的缩放、旋转，平移，组
合变换
（1）思想算法：
在齐次坐标理论下，二维图形几何变换矩阵可用下式表示：
平移变换：[x* y* 1] =[x y 1] * =[t*x s*y 1]
比例变换：[x* y* 1]=[x y 1] * =[m+x n+y 1]
旋转变换：在平面上的二维图形饶原点逆时针旋转Ө角，变换矩阵为
[x* y* 1]=[x y 1] * = [x*cosө-y*sinө]
复合变换：以上各种变换矩阵都是以原点为参照点，当以任意参照点进行变换的时
候，我们就要用到复合变换矩阵。

(2)、程序实现：
用直线命令画出一个齿（或六边形的一半）→利用旋转变换或对称变换矩阵实现对其余部分的绘制→调试运行程序→输出图形→分析结果→结束。

void tf1(float sq,float cq,int xp,int yp)
{ t[0][0]=cq; t[0][1]=sq; t[0][2]=0;
t[1][0]=-sq;t[1][1]=cq; t[1][2]=0;
t[2][0]=xp; t[2][1]=yp; t[2][2]=1;
}
void tf2()
{ float x,y;
x=pp[j][0];y=pp[j][1];
pp[j][0]=x*t[0][0]+y*t[1][0]+t[2][0];
pp[j][1]=x*t[0][1]+y*t[1][1]+t[2][1];
}
main()
{
int gdriver,gmode,xp=300,yp=200;
float p[][2]={{80,1},{80,8},{99,17},{97,26},{76,26}};
float i,t[3][3];
gdriver=DETECT;
registerbgidriver(EGA VGA_driver);
initgraph(&gdriver,&gmode,"d:\\tc");
setbkcolor(3);
cleardevice();
setcolor(9);
circle(xp,yp,70);
circle(xp,yp,20);
circle(xp,yp,40);
circle(xp,yp,5);
circle(xp,yp,2);
circle(xp,yp,45);
setlinestyle(0,0,6);
for(i=0;i<6;i=i+0.3142)
{tf1(sin(i),cos(i),xp,yp);
for(j=0;j<5;j++)
{pp[j][0]=p[j][0];pp[j][1]=p[j][1];tf2();}
for(j=0;j<4;j++)
line(pp[j][0],pp[j][1],pp[j+1][0],pp[j+1][1]);
delay(10000);
}
getch();
closegraph();
}
8．图形三视图变换
编写三维变换算法程序→检查程序的正确性→分段调试程序→输入给出的三维形体各顶点的坐标→执行变换→对算法程序进行必要的调整→更换不同的形体数据继续变换→结束。

程序实现：
#include<graphics.h>
#include<math.h>
int a[14][4]={{30,0,0,1},{30,40,0,1},{0,40,0,1},{0,40,10,1},{0,30,30,1},
{0,0,30,1},{30,0,30,1},{30,10,30,1},{10,10,30,1},{10,30,30,1},{10,40,10,1},
{10,10,10,1},{30,10,10,1},{30,40,10,1}};
float t[4][4],p[14][4];
void a400()
{int i,j;
for(i=0;i<4;i++)
for(j=0;j<4;j++)
t[i][j]=0;
}
void a500()
{ int k,i,j;
for(i=0;i<14;i++)
{for(j=0;j<4;j++)
{p[i][j]=0;
for(k=0;k<4;k++)
p[i][j]=p[i][j]+a[i][k]*t[k][j];}
p[i][0]=p[i][0]+280;
p[i][1]=-p[i][1]+180;
}
setcolor(9);
moveto(p[0][0],p[0][1]);
for(i=0;i<14;i++)
lineto(p[i][0],p[i][1]);
line(p[6][0],p[6][1],p[0][0],p[0][1]);
line(p[7][0],p[7][1],p[12][0],p[12][1]);
line(p[8][0],p[8][1],p[11][0],p[11][1]);
line(p[9][0],p[9][1],p[4][0],p[4][1]);
line(p[10][0],p[10][1],p[3][0],p[3][1]);
line(p[13][0],p[13][1],p[10][0],p[10][1]);
line(p[1][0],p[1][1],p[13][0],p[13][1]);
getch();
}
main()
{int driver,mode,i,j;
driver=DETECT;
initgraph(&driver,&mode,"d:\\tc");
setbkcolor(3);
a400();
t[0][0]=0.7071*3;
t[0][1]=-0.4082*3;
t[1][0]=-0.7071*3;
t[1][1]=-0.4082*3;
t[2][1]=0.8165*3;
t[3][3]=1;
a500();
closegraph();
}
实验总结
在实验过程中，尽管过程中任由许多不会的地方，而且有待于今后的提高和改进，但我加深了对书本上知识的理解与掌握，同时也学到了很多书本上没有东西，并积累了一些宝贵的经验，这对我以后的学习与工作是不无裨益的。

由于本实验采用C语言编程，而C语言在图形界面的制作方面较其他面向对
象语言而言没有优势，但是在实验基本的图形绘制方面，只要选择好算法并加以实施也能够实现。

本实验只是完成了一部分主要的功能，而每种功能都只选用了其中一种算法，但是在完成的过程中，我也很深刻地了解和认识了各种算法，留下了印象。

但是未能实现整个画图板界面还是很遗憾，我将继续努力把它完善。

我需要完成的功能主要有以下几点：
（1）、汉字的编写功能（2）、投影变换（3）、隐藏面消隐（4）、三维图形变换等。