Question

关于像空间坐标系与像空间辅助辅助坐标系的转换的编程的错误需要大家改正一下

#include "stdio.h"
#include "math.h"
#include "stdlib.h"
double *MatrixOpp(double A[],int m,int n)
{
int i,j,x,y,k;
double *SP=NULL,*AB=NULL,*B=NULL,X,*C;
SP=(double *)malloc(m*n*sizeof(double));
AB=(double *)malloc(m*n*sizeof(double));
B=(double *)malloc(m*n*sizeof(double));

X=Surplus(A,m,n);
X=1/X;

for(i=0;i <m;i++)
for(j=0;j <n;j++)
{
for(k=0;k <m*n;k++)
B[k]=A[k];
{
for(x=0;x <n;x++)
B[i*n+x]=0;
for(y=0;y <m;y++)
B[m*y+j]=0;
B[i*n+j]=1;
SP[i*n+j]=Surplus(B,m,n);
AB[i*n+j]=X*SP[i*n+j];
}
}
C=MatrixInver(AB,m,n);

return C;
}
main()
{
double a[3][1],b[3][1],c[3][3],m[1][3];
double w,k,t;
int i,j;
int y,u,p;
for(i=0,i<3,i++)
for(j=0,j<1,j++)
scanf("%10lf",&a[i][j]);
for(i=0,i<3,i++)
for(j=0,j<1,j++)
scanf("%10lf",&b[i][j]);
for(i=0,i<3,i++)
{
printf("\n");
for(j=0,j<1,j++)
printf("%10lf",a[i][j]);
}
for(i=0,i<3,i++)
{
printf("\n");
for(j=0,j<1,j++)
printf("%10lf",b[i][j]);
}
m[1][3]= MatrixOpp(b,i,j);
for (y=0;y<3;y++)
{
for (u=0;u<3;u++)
{
c[y][u] = 0;
for (p=0;p<0;p++)
c[y][u]+=a[y][p]*m[p][u];
}
}
for (y=0;y<3;y++)
{
for (u=0;u<3;u++)
printf("c[y][u]的值:%10lf",c[y][u]);
printf("\n");
}
return 0;
}
w=atan(c[0][2]/c[1][2])*180/3.14;
t=acos(c[2][2])*180/3.14;
k=atan(c[0][1]/c[0][2])*180/3.14;
printf("w,k,t的值是:%10lf%10lf%10lf",w,k,t);
}

Answer 1

#include "stdio.h"
#include "math.h"
#include "stdlib.h"
double *MatrixOpp(double A[],int m,int n)
{
int i,j,x,y,k;
double *SP=NULL,*AB=NULL,*B=NULL,X,*C;
SP=(double *)malloc(m*n*sizeof(double));
AB=(double *)malloc(m*n*sizeof(double));
B=(double *)malloc(m*n*sizeof(double));

X=Surplus(A,m,n);
X=1/X;

for(i=0;i <m;i++)
for(j=0;j <n;j++)
{
for(k=0;k <m*n;k++)
B[k]=A[k];
{
for(x=0;x <n;x++)
B[i*n+x]=0;
for(y=0;y <m;y++)
B[m*y+j]=0;
B[i*n+j]=1;
SP[i*n+j]=Surplus(B,m,n);
AB[i*n+j]=X*SP[i*n+j];
}
}
C=MatrixInver(AB,m,n);

return C;
}
main()
{
double a[3][1],b[3][1],c[3][3],m[1][3];
double w,k,t;
int i,j;
int y,u,p;
for(i=0,i<3,i++)
for(j=0,j<1,j++)
scanf("%10lf",&a[i][j]);
for(i=0,i<3,i++)
for(j=0,j<1,j++)
scanf("%10lf",&b[i][j]);
for(i=0,i<3,i++)
{
printf("\n");
for(j=0,j<1,j++)
printf("%10lf",a[i][j]);
}
for(i=0,i<3,i++)
{
printf("\n");
for(j=0,j<1,j++)
printf("%10lf",b[i][j]);
}
m[1][3]= MatrixOpp(b,i,j);
for (y=0;y<3;y++)
{
for (u=0;u<3;u++)
{
c[y][u] = 0;
for (p=0;p<0;p++)
c[y][u]+=a[y][p]*m[p][u];
}
}
for (y=0;y<3;y++)
{
for (u=0;u<3;u++)
printf("c[y][u]的值:%10lf",c[y][u]);
printf("\n");
}
return 0;
}
w=atan(c[0][2]/c[1][2])*180/3.14;
t=acos(c[2][2])*180/3.14;
k=atan(c[0][1]/c[0][2])*180/3.14;
printf("w,k,t的值是:%10lf%10lf%10lf",w,k,t);
}