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

eobeom

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

最佳答案

月排行榜 / 总排行榜

风过无痕1989

2020-10-14 17:24:48

#include <stdio.h>
#include <math.h>
 
#define MAX 100
 
double D(double arr[][MAX], int n);
void Display(double arr[][MAX], int n);
 
void main()
{
    double arr[MAX][MAX];
    double bak[MAX][MAX];
    double value;
    int n;
    int i,j,k,c;
 
    printf("Please enter matrix size n(1 <= n < 20):");
    scanf("%d",&n);
 
    printf("Please input matrix line by line:\n");
    for(i = 0;i < n;i++)
    {
        for(j = 0;j < n;j++)
        {
            scanf("%lf",&arr[i][j]);
        }
    }
 
    printf("matrix a:\n");
    Display(arr,n);
 
    /* 保存D(0,i)的值到bak中 */
    for(i = 0;i < n;i++)               /* 处理第1行的每一列元素 */
    {
        for(j = 1;j < n;j++)           /* 逐行处理 */
        {
            for(c = 0,k = 0;k < n;k++) /* 逐列处理 */
           {
                if(k == i)
                {
                    continue;
                }
                bak[j - 1][c++] = arr[j][k];
            }
        }
        value = D(bak,n-1);
 
        printf("Submatrix:\n");
        Display(bak,n-1);
 
        printf("DValue of the Submatrix is %f\n",D(bak,n - 1));
    }
 
    value = D(arr,n);
    printf("result = %f\n",value);
}
 
double D(double arr[][MAX],int n)
{
    double bak[MAX][MAX];
    int i,j,k,c;
    double sum = 0;
 
    if(n == 1)
    {
        return arr[0][0];
    }
 
    for(i = 0;i < n;i++)                /* 处理第1行的每一列元素 */
    {
        for(j = 1;j < n;j++)            /* 逐行处理 */
        {
            for(c = 0,k = 0;k < n;k++)  /* 逐列处理 */
            {
                if(k == i)
                {
                    continue;
                }
                bak[j - 1][c++] = arr[j][k];
            }
        }
 
        /* 计算arr(row,col) * D(row,col)的值 */
        sum += (i % 2 == 0 ? 1 : -1) * arr[0][i] * D(bak,n - 1);
    }
    return sum;
}
 
void Display(double arr[][MAX],int n)
{
    int i,j;
    printf("");
    for(i = 0;i < n;i++)
    {
        for(j = 0;j < n;j++)
        {
            printf("%10.1lf",arr[i][j]);
        }
        printf("\n");
    }
}

跳转到最佳答案楼层

kogawananari

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

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

风过无痕1989

#include <stdio.h>
#include <math.h>
 
#define MAX 100
 
double D(double arr[][MAX], int n);
void Display(double arr[][MAX], int n);
 
void main()
{
    double arr[MAX][MAX];
    double bak[MAX][MAX];
    double value;
    int n;
    int i,j,k,c;
 
    printf("Please enter matrix size n(1 <= n < 20):");
    scanf("%d",&n);
 
    printf("Please input matrix line by line:\n");
    for(i = 0;i < n;i++)
    {
        for(j = 0;j < n;j++)
        {
            scanf("%lf",&arr[i][j]);
        }
    }
 
    printf("matrix a:\n");
    Display(arr,n);
 
    /* 保存D(0,i)的值到bak中 */
    for(i = 0;i < n;i++)               /* 处理第1行的每一列元素 */
    {
        for(j = 1;j < n;j++)           /* 逐行处理 */
        {
            for(c = 0,k = 0;k < n;k++) /* 逐列处理 */
           {
                if(k == i)
                {
                    continue;
                }
                bak[j - 1][c++] = arr[j][k];
            }
        }
        value = D(bak,n-1);
 
        printf("Submatrix:\n");
        Display(bak,n-1);
 
        printf("DValue of the Submatrix is %f\n",D(bak,n - 1));
    }
 
    value = D(arr,n);
    printf("result = %f\n",value);
}
 
double D(double arr[][MAX],int n)
{
    double bak[MAX][MAX];
    int i,j,k,c;
    double sum = 0;
 
    if(n == 1)
    {
        return arr[0][0];
    }
 
    for(i = 0;i < n;i++)                /* 处理第1行的每一列元素 */
    {
        for(j = 1;j < n;j++)            /* 逐行处理 */
        {
            for(c = 0,k = 0;k < n;k++)  /* 逐列处理 */
            {
                if(k == i)
                {
                    continue;
                }
                bak[j - 1][c++] = arr[j][k];
            }
        }
 
        /* 计算arr(row,col) * D(row,col)的值 */
        sum += (i % 2 == 0 ? 1 : -1) * arr[0][i] * D(bak,n - 1);
    }
    return sum;
}
 
void Display(double arr[][MAX],int n)
{
    int i,j;
    printf("");
    for(i = 0;i < n;i++)
    {
        for(j = 0;j < n;j++)
        {
            printf("%10.1lf",arr[i][j]);
        }
        printf("\n");
    }
}

乐乐学编程

帮你顶

风过无痕1989

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

乐乐学编程

还没满意？帮你顶上去

风过无痕1989

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

eobeom

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

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

风过无痕1989

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

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