求矩阵的行列式

最新推荐文章于 2025-06-05 19:24:07 发布

原创最新推荐文章于 2025-06-05 19:24:07 发布 · 1.5k 阅读

1 ·

CC 4.0 BY-SA版权

测试、分析、逆向专栏收录该内容

15 篇文章

订阅专栏

本文介绍了一个用Python实现的计算矩阵行列式的函数，并通过一个具体示例进行验证。该函数使用递归方式，适用于任意大小的方阵。同时，为了对比准确性，文章还展示了使用NumPy库内置函数计算同一矩阵行列式的结果。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

import numpy as np

mat_1=[
    [2,4,1,54],
    [23,42,4,4],
    [3,1,4,56],
    [6,3,2,6],
]


mat_2=[
    [100,2,3],
    [4,5,6],
    [7,8,9]
]


def det(mat,n):
    A=0
    if n==1:
        A=mat[0][0]
        return  A
    if n==2:
        A=mat[0][0]*mat[1][1]-mat[0][1]*mat[1][0]
        return A
    
    b=[0 for x  in xrange(n) ]
    t1=[ [0 for x  in xrange(n)]  for y in xrange(n)]

    for i in xrange(n):
        b[i]=mat[i][0]*(-1)**(i+2)

    for i in xrange(n):
        for j in xrange(n-1):
            t1[i][j]=mat[i][j+1]

    for i in xrange(n):
        t2 = [[0 for x in xrange(n-1)] for y in xrange(n-1)]
        for j in xrange(n-1):
            for k in xrange(n-1):
                if  j<i:
                    t2[j][k]=t1[j][k]
                else:
                    t2[j][k]=t1[j+1][k]
        A += b[i]*det(t2,n-1)
    return A



print det(mat_1,len(mat_1))
print np.linalg.det(mat_1)