程序的（旋轉三角形）異常終止

Question

我有以下的C++程序。我必須以某個角度逆時針方向旋轉一個三角形。我通過以下代碼實現了這一點。我的邏輯是我首先輸入旋轉發生的點，旋轉角度，然後座標爲三角形。之後，我應用矩陣公式

M(final coordinates)=M(translate triangle to original position)*M(Rotate triangle)*M(translate triangle to origin)*M(original coordinates) 

其中M指上述公式中的矩陣。
通過乘以這些矩陣以相反的順序，即< - 某人得到了結果。

但是當我運行此代碼程序導致異常終止。

我的代碼是：

#include<iostream.h> 
#include<conio.h> 
#include<graphics.h> 
#include<math.h> 
#define pi 3.14159265 
void rotation(float tx,float ty,float redi,float t[][3]) 
{ 
clrscr(); 
int g=DETECT,gm; 
initgraph(&g,&gm,"C:\\TC\\BGI"); 
setbkcolor(8); 
int i,j,k; 
float x[3][3],c[3][3],r[3][3],mf[3][3],res[3][3]; 
for(j=0;j<3;j++) 
{ 
    for(k=0;k<3;k++) 
    { 
     if(j==k) 
     { 
      r[j][k]=1; 
      x[j][k]=1; 
     } 
     else 
     { 
      r[j][k]=0; 
      x[j][k]=0; 
     } 
    } 
} 
float co,si; 
co=cos(redi); 
si=sin(redi); 
r[0][0]=co; r[0][1]=-si; r[1][0]=si; r[1][1]=co; 

x[0][2]=-tx; 
x[1][2]=-ty; 
for(int mat=0;mat<3;mat++) 
{ 
    for(i=0;i<3;i++) 
    { 
     for(j=0;j<3;j++) 
     { 
      if(mat==0) c[i][j]=0; 
      if(mat==1) mf[i][j]=0; 
      if(mat==2) res[i][j]=0; 
      for(k=0;k<3;k++) 
      { 
       if(mat==0) 
        c[i][j]+=r[i][k]*x[k][j]; 
       if(mat==1) 
       { 
        x[0][2]=tx; 
        x[1][2]=ty; 
        mf[i][j]+=x[i][k]*c[k][j]; 
       } 
       if(mat==2) 
        res[i][j]+=mf[i][k]*t[k][j]; 

      } 
      cout<<res[i][j]; 
     } 
    } 
} 
line(res[0][0],res[1][0],res[0][1],res[1][1]); 
line(t[0][0],t[1][0],t[0][1],t[1][1]); 
line(res[0][2],res[1][2],res[0][1],res[1][1]); 
line(t[0][2],t[1][2],t[0][1],t[1][1]); 
line(res[0][2],res[1][2],res[0][0],res[1][0]); 
line(t[0][2],t[1][2],t[0][0],t[1][0]); 
for(i=0;i<3;i++) 
{ 
    putpixel(res[0][i],res[1][i],RED); 
    putpixel(t[0][i],t[1][i],RED); 
} 
closegraph(); 
getch(); 

} 
void main() 
{ 
clrscr(); 
float tx,ty,t[3][3]; 
float deg,redi; 
cout<<"Rotation about point : "; 
cin>>tx>>ty; 
cout<<"Enter angle of rotation :"; 
cin>>deg; 
for(int i=0;i<3;i++) 
{ 
    cout<<"Enter co-ordintes "<<i+1<<" of triangle: "; 
    for(int j=0;j<3;j++) 
    { 
     if(j==2) 
      t[j][i]=1; 
     else 
      cin>>t[j][i]; 
    } 
} 
redi=(deg*pi)/180.0; 
rotation(tx,ty,redi,t); 
getch(); 
} 

我用TC++ IDE。 請幫助......

Answer 1

我認爲有一些問題，在您的這些循環檢查出來

for(int mat=0;mat<3;mat++) 
{ 
    for(i=0;i<3;i++) 
    { 
    for(j=0;j<3;j++) 
    { 
     if(mat==0) c[i][j]=0; 
     if(mat==1) mf[i][j]=0; 
     if(mat==2) res[i][j]=0; 
     for(k=0;k<3;k++) 
     { 
      if(mat==0) 
       c[i][j]+=r[i][k]*x[k][j]; 
      if(mat==1) 
      { 
       x[0][2]=tx; 
       x[1][2]=ty; 
       mf[i][j]+=x[i][k]*c[k][j]; 
      } 
      if(mat==2) 
       res[i][j]+=mf[i][k]*t[k][j]; 

     } 
     cout<<res[i][j]; 
    } 
    } 
} 

乘以你的矩陣逐一（使用每乘不同的環路）並運行再次編程。

OR

使用此修改後的代碼：

#include<iostream.h> 
#include<conio.h> 
#include<graphics.h> 
#include<math.h> 
void main() 
{ 
clrscr(); 
int gd=DETECT,gm; 
initgraph(&gd,&gm,"C:\\TC\\BGI"); 
float x[3][3],y[3][3],z[3][3],p[3][3],q[3][3]; 
float angle,ptx,pty; 
int i,j,k; 
cout<<"Enter coordinates of triangle:"; 
for(i=0;i<3;i++) 
{ 
    for(j=0;j<2;j++) 
    { 
     cin>>y[j][i]; 
    } 
} 
y[2][0]=1; 
y[2][1]=1; 
y[2][2]=1; 
cout<<"Enter the point about:"; 
cin>>ptx>>pty; 
for(i=0;i<3;i++) 
{ 
    for(j=0;j<3;j++) 
    { 
     if(i==j) 
      x[i][j]=1; 
     else 
      x[i][j]=0; 
    } 
} 
x[0][2]=-ptx; 
x[1][2]=-pty; 
for(i=0;i<3;i++) 
{ 
    for(j=0;j<3;j++) 
    { 
     z[i][j]=0; 
     for(int k=0;k<3;k++) 
     { 
      z[i][j]+=x[i][k]*y[k][j]; 
     } 
    } 
} 
cout<<"Enter angle of rotation:"; 
cin>>angle; 
angle*=(3.14/180); 
for(i=0;i<3;i++) 
{ for(j=0;j<3;j++) 
    { 
     if(i==j) 
      x[i][j]=1; 
     else 
      x[i][j]=0; 
    } 
} 
x[0][0]=cos(angle); 
x[0][1]=-sin(angle); 
x[1][0]=sin(angle); 
x[1][1]=cos(angle); 
for(i=0;i<3;i++) 
{ 
    for(j=0;j<3;j++) 
    { 
     p[i][j]=0; 
     for(int k=0;k<3;k++) 
     { 
      p[i][j]+=(x[i][k]*z[k][j]); 
     } 
    } 
} 
for(i=0;i<3;i++) 
{ 
    for(j=0;j<3;j++) 
    { 
     if(i==j) 
      x[i][j]=1; 
     else 
      x[i][j]=0; 
    } 
} 
x[0][2]=ptx; 
x[1][2]=pty; 
for(i=0;i<3;i++) 
{ 
    for(j=0;j<3;j++) 
    { 
     q[i][j]=0; 
     for(int k=0;k<3;k++) 
     { 
      q[i][j]+=(x[i][k]*p[k][j]); 
     } 
    } 
} 
clrscr(); 
line(y[0][0],y[1][0],y[0][1],y[1][1]); 
line(y[0][0],y[1][0],y[0][2],y[1][2]); 
line(y[0][1],y[1][1],y[0][2],y[1][2]); 

line(q[0][0],q[1][0],q[0][1],q[1][1]); 
line(q[0][0],q[1][0],q[0][2],q[1][2]); 
line(q[0][1],q[1][1],q[0][2],q[1][2]); 
getch(); 
}