您的当前位置:首页正文

矩阵函数

来源:九壹网
/****************************** 矩阵相乘

******************************/

void plus(double *A, double *B, double *C, int m, int n, int R) { int i,j,k; for(i=0;i/****************************** 矩阵转置

矩阵a 为m*n的矩阵, 返回矩阵a_t n*m

******************************/

void transpose(double *a, double*a_t,int m,int n) { int i,j; for(i=0;i*(a_t+i+j*m)=*(a+j+i*n); } return; }

/****************************** 矩阵求逆

矩阵a 为m*n的矩阵, 返回矩阵e 为a的逆矩阵

矩阵a可逆的充要条件是a是方阵且a为非奇异矩阵 ******************************/

void inverse(double *a,double*e,int m,int n) { int i,j; double temp; double temp2; /*若该矩阵不是方阵,则不能求逆*/ if(m!=n) { cout<<\"a不是方阵,不能求逆!\"<{ temp2=*(a+i+k*m); for(j=0;j0;i--) { for(j=0;j=0;k--) { temp2=*(a+i+k*m); for(j=0;j

因篇幅问题不能全部显示,请点此查看更多更全内容