C语言-三维图形变换

{
float bf,af,ca,cb,cg,sg,xp,yp,zp;
int i;
bf=bta*rad;
af=afa*rad;
cb=cos(bf);
ca=cos(af);
for(i=0;i<mp;i++)
{
if(x[i]==0)
xp=ox-sf*z[i]*cb;
else
xp=ox+sf*(x[i]-z[i]*cb);
void rotationy(float xc,float yc,float zc,float beta,float tranm[][4]);
void rotationz(float xc,float yc,float zc,float gama,float tranm[][4]);
void projectx(float i,float m,float n,float tranm[][4]);
axono2d(1.,320.,240.,-45.,35.2667,u,v,w,ipen,kind,NUMBER);
mirrxoy(0.,0.,-60.,tm);
draw(1.,320.,240.,tm);
mirrxoz(0.,-60.,0.,tm);
draw(1.,320.,240.,tm);
mirryoz(-60.,0.,0.,tm);
{
for(mj=0;mj<4;mj++)
{
mc[mi][mj]=0;
for(mk=0;mk<4;mk++)
mc[mi][mj]=mc[mi][mj]+ma[mi][mk]*mb[mk][mj];
}
}
}
void transform(float x0,float y0,float z0,float *x,float *y,float *z,float tranm[][4])
setwritemode(1);
initm(tm);
for(j=1;j<=15;j++)
{
setcolor(j);
translation(0.,0.,-10.*j,tm);
draw(1.,300.,350.,tm);
}
for(j=1;j<=15;j++)
{
setcolor(j);
translation(-10.*j,0.,0.,tm);
void mirrxoz(float l,float m,float n,float tranm[][4]);
void mirryoz(float l,float m,float n,float tranm[][4]);
void shearing(float a,float b,float c,float d,float e,float f,float tranm[][4]);
tranm[3][0]=l;
tranm[3][1]=m;
tranm[3][2]=n;
}
void scaling(float sx,float sy,float sz,fltm(tranm);
tranm[0][0]=sx;
tranm[1][1]=sy;
tranm[2][2]=sz;
}
void rotationx(float xc,float yc,float zc,float alfa,float tranm[][4])
{
initm(tranm);
tranm[1][1]=cos(rad*alfa);
tranm[1][2]=sin(rad*alfa);
tranm[2][1]=-tranm[1][2];
v[NUMBER]={20.,60.,60.,0.,0.,20.,20.,20.,20.,20.,0.,0.,20.,20.,60.,60.,60.,60.,60.,60.,0.,0.,0.,0.},
w[NUMBER]={30.,30.,0.,0.,90.,90.,30.,30.,90.,90.,90.,90.,90.,30.,30.,0.,0.,30.,30.,0.,0.,0.,0.,90.};
else
{
if(lk==0)
setlinestyle(USERBIT_LINE,0xFFFF,3);
else
setlinestyle(USERBIT_LINE,0xFFFF,3);
lineto((int)x,(int)y);
}
}
void initm(float unitm[][4])
{
int n;
for(n=0;n<4;n++)
{
unitm[n][0]=0;
unitm[n][1]=0;
unitm[n][2]=0;
unitm[n][3]=0;
unitm[n][n]=1;
}
}
void mult4x4(float ma[][4],float mb[][4],float mc[][4])
{
int mi,mj,mk;
for(mi=0;mi<4;mi++)
void projecty(float i,float m,float n,float tranm[][4]);
void projectz(float i,float m,float n,float tranm[][4]);
void mirrxoy(float l,float m,float n,float tranm[][4]);
int ipen[NUMBER]={3,2,2,2,2,2,2,2,2,2,3,2,2,3,2,2,2,3,2,3,2,2,3,2};
int kind[NUMBER]={0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,0,1};
void draw(float sf,float xc,float yc,float tranm[][4]);
getch();
axono2d(sf,xc,yc,-45.,35.2667,u1,v1,w1,ipen,kind,m);
}
void axono2d(float sf,float ox,float oy,float bta,float afa,float x[],float y[],float z[],int kp[],int kl[],int mp)
h=tranm[0][3]*x0+tranm[1][3]*y0+tranm[2][3]*z0+tranm[3][3];
(*x)=xu/h;
(*y)=yv/h;
(*z)=zw/h;
}
void translation(float l,float m,float n,float tranm[][4])
{
initm(tranm);
draw(1.,320.,240.,tm);
getch();
closegraph();
}
void draw(float sf,float xc,float yc,float tranm[][4])
{
int i,m=NUMBER;
float *px,*py,*pz,xu,yv,zw;
float u1[NUMBER],v1[NUMBER],w1[NUMBER];
{
float xu,yv,zw,h;
xu=tranm[0][0]*x0+tranm[1][0]*y0+tranm[2][0]*z0+tranm[3][0];
yv=tranm[0][1]*x0+tranm[1][1]*y0+tranm[2][1]*z0+tranm[3][1];
zw=tranm[0][2]*x0+tranm[1][2]*y0+tranm[2][2]*z0+tranm[3][2];
main()
{
int j,m=NUMBER;
int gd=DETECT,gm;
float tm[4][4]={0.};
initgraph(&gd,&gm,"d:\\tc");
settextstyle(3,0,2);
outtextxy(170,4,"3D GRAPFICS TRANSFORMATION");
draw(1.,320.,240.,tm);
}
for(j=1;j<=15;j++)
{
setcolor(j);
rotationy(0.,0.,10.,24.*j,tm);
draw(1.,320.,240.,tm);
}
for(j=1;j<=15;j++)
{
setcolor(j);
rotationz(0.,0.,10.,24.*j,tm);
void initm(float unitm[][4]);
void mult4x4(float ma[][4],float mb[][4],float mc[][4]);
void transform(float x0,float y0,float z0,float *x,float *y,float *z,float tranm[][4]);
void axono2d(float sf,float ox,float oy,float bta,float afa,float x[],float y[],float z[],int kp[],int kl[],int mp);
void plot(float x,float y,int ip,int lk);