Python实现求矩阵的伴随矩阵

Python实现求矩阵的伴随矩阵

本文介绍如何使用python实现矩阵的输入并求出矩阵的伴随矩阵。

伴随矩阵

设有一矩阵A

设Aij是矩阵A中元素aij的代数余子式，那么矩阵A*

称为矩阵A的伴随矩阵

基本思路

先输入初始矩阵的阶数，再进行矩阵的录入
利用循环分别求矩阵每个元素对应的代数余子式，将所求结果按原对应位置组成新的矩阵，再进行转置，所求即初始矩阵的伴随矩阵

代码实现

import numpy as np

#矩阵输入
def arry_input(n,m):
    a = []
    b = []
    for i in range(n):
        for j in range(m):
            b.append(float(input()))
        a.append(b)
        b = []
    return a
    
#矩阵转置
def trans(m):
    a = [[] for i in m[0]]
    for i in m:
        for j in range(len(i)):
            a[j].append(i[j])
    return a
    
#矩阵复制
def copy(a):
    b = []
    c = []
    for i in range(len(a)):
        for j in range(len(a)):
            b.append(a[i][j])
        c.append(b)
        b = []
    return c

n = int(input())#矩阵的阶数
a = arry_input(n,n)#开始进行矩阵的输入
b = []
c = []
d = []

for i in range(n):
    for j in range(n):
        b = copy(a)
        for x in range(n):
            b[i][x] = 0
        for y in range(n):
            b[y][j] = 0
        #把aij元素的值变为1，其同行同列的其他元素变为0
        b[i][j] = 1       
        #求代数余子式的值，保留两位小数
        index = round(np.linalg.det(b),2)
        d.append(index)
    c.append(d)
    d = []
    
print(trans(c))

需要注意的是

python中矩阵的复制不能简单地通过‘=’进行赋值，否则改变一个矩阵中元素的值，另一个矩阵中也会做出对应的改变
本代码还将每个元素的代数余子式的输出结果保留了两位小数，如对输出精度要求较高，可以自行调整