C++求矩阵的逆程序代码(最新整理)

1、求解逆矩阵源代码：#include #include #include #definetiny 1.0e-20void inverse(double*,int);void ludcmp(double*, int, int*, double*); void lubksb(double*, int, int*, double*); double *matrix(int,int,int,int);double *vector(int,int);voidfree_matrix(double*,int,int,int,int); voidfree_vector(double*,int,int);void

2、 inverse(double *mat, int dim)int i,j,*indx; double *y,d,*col;y = matrix(0,dim-1,0,dim-1);indx = (int *)malloc(unsigned)(dim*sizeof(int); col = vector(0,dim-1);ludcmp(mat,dim,indx,&d); for (j=0;jdim;j+)for (i=0;idim;i+) coli = 0.0; colj = 1.0; lubksb(mat,dim,indx,col);for (i=0;idim;i+) yij = coli;fo

3、r (i=0;idim;i+) for (j=0;jdim;j+)matij = yij; free_matrix(y,0,dim-1,0,dim-1); free_vector(col,0,dim-1); free(indx);void ludcmp(double *a, int n, int *indx, double *d)int i,imax,j,k;doublebig,dum,sum,temp; double*vv;vv = (double*)malloc(unsigned)(n*sizeof(double); if (!vv)fprintf(stderr,error allocat

4、ing vector memoryn); exit(1);*d = 1.0;for (i=0;in;i+)big = 0.0;for (j=0;j big) big = temp;if (big = 0.0)fprintf(stderr,singular matrix in routine ludcmpn); for (j=0;jn;j+) printf( %f ,aij); printf(/n);exit(1);vvi = 1.0/big;for (j=0;jn;j+)for (i=0;ij;i+)sum = aij;for (k=0;ki;k+) sum -= aik * akj; aij

5、 = sum;big = 0.0;for (i=j;in;i+)sum = aij;for (k=0;k= big)big = dum; imax = i;if (j != imax)for (k=0;kn;k+)dum = aimaxk; aimaxk = ajk; ajk = dum;*d = -(*d);vvimax = vvj;indxj = imax;if (ajj = 0.0) ajj = tiny;if (j != n-1)dum = 1.0 / ajj;for (i=j+1;in;i+) aij *= dum;free(vv);void lubksb(double *a, in

6、t n, int *indx, double *b)int i,ip,j,ii=-1; doublesum;for (i=0;i=0)for (j=ii;j=0;i-)sum = bi;for (j=i+1;jn;j+) sum -= aij * bj; bi = sum / aii;“”“”at the end, xiao bian gives you a passage. minand once said, people who learn to learn are very happy people. in every wonderful life, learning is an ete

7、rnal theme. as a professional clerical and teaching position, i understand the importance of continuous learning, life is diligent, nothing can be gained, only continuous learning can achieve better self. only by constantly learning and mastering the latest relevant knowledge, can employees from all walks of life keep up with the pace of enterprise development and innovate to meet the needs of the market. this document is also ed