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

需要注意的是

python中矩阵的复制不能简单地通过‘=’进行赋值，否则改变一个矩阵中元素的值，另一个矩阵中也会做出对应的改变

本代码还将每个元素的代数余子式的输出结果保留了两位小数，如对输出精度要求较高，可以自行调整