[已解决]小白求助 如何计算矩阵行列式的值？

eobeom

请问计算矩阵行列式的值的程序该怎么写？

最佳答案

月排行榜 / 总排行榜

风过无痕1989

2020-10-14 17:24:48

1. #include <stdio.h>
   
2. #include <math.h>
   
3. 
   
4. #define MAX 100
   
5. 
   
6. double D(double arr[][MAX], int n);
   
7. void Display(double arr[][MAX], int n);
   
8. 
   
9. void main()
   
10. {
    
11.     double arr[MAX][MAX];
    
12.     double bak[MAX][MAX];
    
13.     double value;
    
14.     int n;
    
15.     int i,j,k,c;
    
16. 
    
17.     printf("Please enter matrix size n(1 <= n < 20):");
    
18.     scanf("%d",&n);
    
19. 
    
20.     printf("Please input matrix line by line:\n");
    
21.     for(i = 0;i < n;i++)
    
22.     {
    
23.         for(j = 0;j < n;j++)
    
24.         {
    
25.             scanf("%lf",&arr[i][j]);
    
26.         }
    
27.     }
    
28. 
    
29.     printf("matrix a:\n");
    
30.     Display(arr,n);
    
31. 
    
32.     /* 保存D(0,i)的值到bak中 */
    
33.     for(i = 0;i < n;i++)               /* 处理第1行的每一列元素 */
    
34.     {
    
35.         for(j = 1;j < n;j++)           /* 逐行处理 */
    
36.         {
    
37.             for(c = 0,k = 0;k < n;k++) /* 逐列处理 */
    
38.            {
    
39.                 if(k == i)
    
40.                 {
    
41.                     continue;
    
42.                 }
    
43.                 bak[j - 1][c++] = arr[j][k];
    
44.             }
    
45.         }
    
46.         value = D(bak,n-1);
    
47. 
    
48.         printf("Submatrix:\n");
    
49.         Display(bak,n-1);
    
50. 
    
51.         printf("DValue of the Submatrix is %f\n",D(bak,n - 1));
    
52.     }
    
53. 
    
54.     value = D(arr,n);
    
55.     printf("result = %f\n",value);
    
56. }
    
57. 
    
58. double D(double arr[][MAX],int n)
    
59. {
    
60.     double bak[MAX][MAX];
    
61.     int i,j,k,c;
    
62.     double sum = 0;
    
63. 
    
64.     if(n == 1)
    
65.     {
    
66.         return arr[0][0];
    
67.     }
    
68. 
    
69.     for(i = 0;i < n;i++)                /* 处理第1行的每一列元素 */
    
70.     {
    
71.         for(j = 1;j < n;j++)            /* 逐行处理 */
    
72.         {
    
73.             for(c = 0,k = 0;k < n;k++)  /* 逐列处理 */
    
74.             {
    
75.                 if(k == i)
    
76.                 {
    
77.                     continue;
    
78.                 }
    
79.                 bak[j - 1][c++] = arr[j][k];
    
80.             }
    
81.         }
    
82. 
    
83.         /* 计算arr(row,col) * D(row,col)的值 */
    
84.         sum += (i % 2 == 0 ? 1 : -1) * arr[0][i] * D(bak,n - 1);
    
85.     }
    
86.     return sum;
    
87. }
    
88. 
    
89. void Display(double arr[][MAX],int n)
    
90. {
    
91.     int i,j;
    
92.     printf("");
    
93.     for(i = 0;i < n;i++)
    
94.     {
    
95.         for(j = 0;j < n;j++)
    
96.         {
    
97.             printf("%10.1lf",arr[i][j]);
    
98.         }
    
99.         printf("\n");
    
100.     }
     
101. }

复制代码

跳转到最佳答案楼层

kogawananari

你到底是求行列式的值 还是求增广矩阵对应的线性方程组的解

行列式的值np.linalg.det
增广矩阵求解np.linalg.solve

风过无痕1989

1. #include <stdio.h>
   
2. #include <math.h>
   
3. 
   
4. #define MAX 100
   
5. 
   
6. double D(double arr[][MAX], int n);
   
7. void Display(double arr[][MAX], int n);
   
8. 
   
9. void main()
   
10. {
    
11.     double arr[MAX][MAX];
    
12.     double bak[MAX][MAX];
    
13.     double value;
    
14.     int n;
    
15.     int i,j,k,c;
    
16. 
    
17.     printf("Please enter matrix size n(1 <= n < 20):");
    
18.     scanf("%d",&n);
    
19. 
    
20.     printf("Please input matrix line by line:\n");
    
21.     for(i = 0;i < n;i++)
    
22.     {
    
23.         for(j = 0;j < n;j++)
    
24.         {
    
25.             scanf("%lf",&arr[i][j]);
    
26.         }
    
27.     }
    
28. 
    
29.     printf("matrix a:\n");
    
30.     Display(arr,n);
    
31. 
    
32.     /* 保存D(0,i)的值到bak中 */
    
33.     for(i = 0;i < n;i++)               /* 处理第1行的每一列元素 */
    
34.     {
    
35.         for(j = 1;j < n;j++)           /* 逐行处理 */
    
36.         {
    
37.             for(c = 0,k = 0;k < n;k++) /* 逐列处理 */
    
38.            {
    
39.                 if(k == i)
    
40.                 {
    
41.                     continue;
    
42.                 }
    
43.                 bak[j - 1][c++] = arr[j][k];
    
44.             }
    
45.         }
    
46.         value = D(bak,n-1);
    
47. 
    
48.         printf("Submatrix:\n");
    
49.         Display(bak,n-1);
    
50. 
    
51.         printf("DValue of the Submatrix is %f\n",D(bak,n - 1));
    
52.     }
    
53. 
    
54.     value = D(arr,n);
    
55.     printf("result = %f\n",value);
    
56. }
    
57. 
    
58. double D(double arr[][MAX],int n)
    
59. {
    
60.     double bak[MAX][MAX];
    
61.     int i,j,k,c;
    
62.     double sum = 0;
    
63. 
    
64.     if(n == 1)
    
65.     {
    
66.         return arr[0][0];
    
67.     }
    
68. 
    
69.     for(i = 0;i < n;i++)                /* 处理第1行的每一列元素 */
    
70.     {
    
71.         for(j = 1;j < n;j++)            /* 逐行处理 */
    
72.         {
    
73.             for(c = 0,k = 0;k < n;k++)  /* 逐列处理 */
    
74.             {
    
75.                 if(k == i)
    
76.                 {
    
77.                     continue;
    
78.                 }
    
79.                 bak[j - 1][c++] = arr[j][k];
    
80.             }
    
81.         }
    
82. 
    
83.         /* 计算arr(row,col) * D(row,col)的值 */
    
84.         sum += (i % 2 == 0 ? 1 : -1) * arr[0][i] * D(bak,n - 1);
    
85.     }
    
86.     return sum;
    
87. }
    
88. 
    
89. void Display(double arr[][MAX],int n)
    
90. {
    
91.     int i,j;
    
92.     printf("");
    
93.     for(i = 0;i < n;i++)
    
94.     {
    
95.         for(j = 0;j < n;j++)
    
96.         {
    
97.             printf("%10.1lf",arr[i][j]);
    
98.         }
    
99.         printf("\n");
    
100.     }
     
101. }

复制代码

乐乐学编程

帮你顶

风过无痕1989

如果满意了，回来选个最佳答案哈~~你这个举手之劳的动作不需要你付出任何(鱼币、荣誉值、积分)，由系统给我们进行奖励：鱼币+3，荣誉+2，技术+1

乐乐学编程

还没满意？帮你顶上去

风过无痕1989

还有什么不满意的，你倒是说话呀

eobeom

风过无痕1989 发表于 2020-10-17 14:40
还有什么不满意的，你倒是说话呀

不好意思我才看到，我明白了谢谢

风过无痕1989

eobeom 发表于 2020-10-17 19:48
不好意思我才看到，我明白了谢谢

明白了，就给个最佳答案呀