设为首页 加入收藏

TOP

三元组矩阵行列式的计算(用递归)
2018-11-25 10:08:23 】 浏览:46
Tags:三元 矩阵 行列 计算

1.具体思想:

关于计算矩阵行列式有两个主要方法:

1.根据居住行列式的定义式用递归计算(就是本文所讲)

2.先做矩阵行变换,转化为上三角矩阵,再求行列式。

 

(我先是思考了行变换转化为三角矩阵,但中途遇到了些问题,所以先把递归的方法写下来,之后会继续更新另外一种方法。)

线性代数里我们已经了解了递归求矩阵行列式的方法。下图:

然后每一个代数余子式又可以看做相对于“n阶母矩阵”的“n-1阶子矩阵”,再次对这个子矩阵按照行或列展开,这就是递归求矩阵行列式的思想。

2.三元组和二维数组

二维数组不用多说,它和矩阵是一一对应的,表示完全相同。

三元组是指形如((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 } 

 


编程开发网
】【打印繁体】【投稿】【收藏】 【推荐】【举报】【评论】 【关闭】 【返回顶部
上一篇C语言 经典编程100 下一篇编程测试计算机存储的大小端模式

评论

帐  号: 密码: (新用户注册)
验 证 码:
表  情:
内  容:

array(4) { ["type"]=> int(8) ["message"]=> string(24) "Undefined variable: jobs" ["file"]=> string(32) "/mnt/wp/cppentry/do/bencandy.php" ["line"]=> int(214) }