标签：根据 img return info src ++ matrix 一个 style
关于计算矩阵行列式有两个主要方法：
1.根据居住行列式的定义式用递归计算（就是本文所讲）
2.先做矩阵行变换，转化为上三角矩阵，再求行列式。
（我先是思考了行变换转化为三角矩阵，但中途遇到了些问题，所以先把递归的方法写下来，之后会继续更新另外一种方法。）
线性代数里我们已经了解了递归求矩阵行列式的方法。下图：
然后每一个代数余子式又可以看做相对于“n阶母矩阵”的“n-1阶子矩阵”，再次对这个子矩阵按照行或列展开，这就是递归求矩阵行列式的思想。
二维数组不用多说，它和矩阵是一一对应的，表示完全相同。
三元组是指形如（（r，c），d）的集合，我们规定(r,c)是三元组中的一个数在二维数组中的对应位置，d表示数据的值。
三元组的数据结构如下：
1 typedef struct 2 { 3 int r; 4 int c; 5 int d; 6 }TupNode; 7 typedef struct 8 { 9 int rows; 10 int cols; 11 int nums; 12 TupNode data[Maxsize]; 13 }TSMatrix;
具体的计算行列式代码如下：
1 //计算矩阵行列式 2 int DatMat(TSMatrix t){ 3 if(t.cols!=t.rows){ 4 printf("该矩阵无法求行列式！"); 5 return 0; 6 } 7 else{ 8 int n=t.cols; 9 int a[n][n]; 10 //将三元组转化为二维数组 11 for(int i=0;i<n;i++){ 12 for(int j=0;j<n;j++){ 13 a[i][j]=0; 14 } 15 } 16 for(int k=0;k<t.nums;k++){ 17 int i = t.data[k].r; 18 int j = t.data[k].c; 19 a[i][j] = t.data[k].d; 20 } 21 if (n == 1){ 22 return a[0][0]; 23 } 24 else{ 25 int b[n-1][n-1];//创建n-1阶的代数余子式阵bb 26 int c[(n-1)*(n-1)]; 27 int sum = 0;//sum为行列式的值 28 TSMatrix t1; 29 30 //以第一列为基础，求行列式 31 for(int l=0;l<n;l++){ 32 int m1=0; 33 int m2=0; 34 for(int i =0;i<n;i++){ 35 for(int j=0;j<n;j++){ 36 if(i!=l&&j!=0){ 37 c[m1]=a[i][j]; 38 m1++; 39 } 40 } 41 } 42 for(int i =0;i<n-1;i++){ 43 for(int j=0;j<n-1;j++){ 44 b[i][j]=c[m2]; 45 m2++; 46 } 47 } 48 49 //把二维数组转化为三元组 50 t1.rows=n-1; 51 t1.cols=n-1; 52 t1.nums=0; 53 for(int i=0;i<n-1;i++){ 54 for(int j=0;j<n-1;j++){ 55 if(b[i][j]!=0){ 56 t1.data[t1.nums].r=i; 57 t1.data[t1.nums].c=j; 58 t1.data[t1.nums].d=b[i][j]; 59 t1.nums++; 60 } 61 } 62 } 63 sum+=a[l][0]*DatMat(t1)*pow(-1,l);//通过递归来求行列式的值 64 } 65 return sum; 66 } 67 } 68 }
标签：根据 img return info src ++ matrix 一个 style
原文地址：https://www.cnblogs.com/OrdinaryMan/p/10012283.html
