2022年 11月 5日

Python计算行列式

 
  1. from itertools import combinations,permutations
  2. #检测
  3. #生成一个排序函数
  4. #计算逆序数
  5. #计算行列式
  6.  
  7. Matrix = [[1,5,3],[7,4,3],[2,4,232]]
  8. #检测
  9. def test(Matrix):
  10.     #检测是否是一个方阵
  11.     for i in Matrix:
  12.         if len(i)==len(Matrix):
  13.             pass
  14.             for ii in i:
  15.                 if type(ii)==type(101):
  16.                     pass
  17.  
  18.                 else:
  19.                     print("请检查您输入的矩阵")
  20.                     print("不可靠的结果")
  21.         else:
  22.             print("请检查您输入的矩阵")
  23.             print("不可靠的结果")
  24.             break
  25.             
  26.     #输出行列式阶数和n!
  27.     sum = 1
  28.     for i in range(len(Matrix)+1):
  29.         sum *= i
  30.         
  31.     return len(Matrix),sum
  32.  
  33.  
  34. def InversionNumber(l):
  35.     sum = 0
  36.     l0 = []
  37.     for i in range(len(l)):
  38.         if l[i] < l[i-1] and i-1 >= 0:
  39.            l0.append(l[i-1])
  40.            sum += 1
  41.  
  42.         for k in range(i-1):
  43.             if l[i] < l[k]:
  44.                 sum += 1
  45.  
  46.     return (-1)**sum
  47.  
  48.  
  49. def Permute(n):
  50.     #返回一个列表,包含整个序列,最后一位表明(-1)^ε
  51.     l = []
  52.     for k in range(1,n+1):
  53.             l.append(k)
  54.     L = list(permutations(l, len(l)))
  55.     for i in L:
  56.         L[L.index(i)]=list(i)
  57.  
  58.     for i in L:
  59.         i.append(InversionNumber(i))
  60.  
  61.     return L
  62.  
  63. def main(Matrix):
  64.     l,S = test(Matrix)
  65.     L = Permute(l)
  66.     sum = 0
  67.     for i in L:
  68.         sum1 = 1
  69.         for ii in range(len(Matrix)):
  70.             sum1 *= Matrix[ii][i[ii]-1]
  71.         sum1 *= i[-1]
  72.         sum += sum1
  73.     return sum
  74.  
  75. if __name__ == "__main__":
  76.     det = main(Matrix)

这个提示一下,咱们得先修改最上头那一段Matrix

写入您想计算的那个哦!