矩阵转逆及弧度角度变换程序.docVIP

矩阵的逆程序： #includestdio.h void main() { int N; printf(输入不超过10的矩阵的阶数N:\n); scanf(%d,N); float a[10][10],b[10][20],c[10][10],t; int i,j,m; printf(请输入行列式不为0的矩阵A(%d阶):\n,N); for(i=0;iN;i++) for(j=0;jN;j++) scanf(%f,a[i][j]); for(i=0;iN;i++) for(j=0;jN;j++) b[i][j]=a[i][j]; for(i=0;iN;i++) for(j=N;j2*N;j++) b[i][j]=0; for(i=0;iN;i++) b[i][N+i]=1; for(m=0;mN;m++) { t=b[m][m]; i=m; while(b[m][m]==0) { b[m][m]=b[i+1][m]; i++; } if(im) { b[i][m]=t; for(j=0;jm;j++) { t=b[m][j]; b[m][j]=b[i][j]; b[i][j]=t; } for(j=m+1;j2*N;j++) { t=b[m][j]; b[m][j]=b[i][j]; b[i][j]=t; } } for(i=m+1;iN;i++) for(j=2*N-1;j=m;j--) b[i][j]-=b[i][m]*b[m][j]/b[m][m]; for(j=2*N-1;j=m;j--) b[m][j]/=b[m][m]; } printf(第一步变换后得到的增广矩阵为：\n); for(i=0;iN;i++) { for(j=0;j2*N;j++) printf(%3.5f ,b[i][j]); printf(\n); } m=N-1; while(m0) { for(i=0;im;i++) for(j=2*N-1;j=m;j--) b[i][j]-=b[i][m]*b[m][j]; m--; } printf(最后得到的增广矩阵为：\n); for(i=0;iN;i++) { for(j=0;j2*N;j++) printf(%3.5f ,b[i][j]); printf(\n); } for(i=0;iN;i++) for(j=0;jN;j++) c[i][j]=b[i][N+j]; printf(故逆矩阵为：\n); for(i=0;iN;i++) { for(j=0;jN;j++) printf(%3.5f ,c[i][j]); printf(\n); } } 弧度角度转换程序： #include math.h #include iostream #define PI 3.1415926 //i是起点，j是终点 double Azimuth(double xi,double yi,double xj,double yj) { double Dx,Dy,S,T; Dx=xj-xi; Dy=yj-yi; S=sqrt(Dx*Dx+Dy*Dy); T=asin(Dy/S); if(Dx0Dy0) T=T; else if(Dx0Dy0) T+=PI/2; else if(Dx0Dy0) T=2*PI/3-T; else T+=2*PI/3; return T; } using namespace std; int main() { double xi,yi,xj,yj; coutplease input xi:endl; cinxi; coutplease input yi:endl; cinyi; coutplease input xj:endl; cinxj; coutplease input yj:endl; cinyj; double azimuth=Azimuth(xi,yi,x