python中矩陣的基本運算學習記錄
阿新 • • 發佈:2018-12-20
矩陣運算:
NumPy系統是Python的一種開源的數值計算擴充套件。這種工具可用來儲存和處理大型矩陣,比Python自身的巢狀列表(nested list structure)結構要高效的多(該結構也可以用來表示矩陣(matrix))
import numpy as np
from numpy.linalg import * // 矩陣運算
行列式:
a=np.mat([[1,23,3],[1,6,4],[13,12,14]])
a_det=det(a)
求逆:
a_inv=a.I
伴隨矩陣:
a_adju=a_det*a_inv
求兩個數的最大公因數:
def gcd(a,b):
while a!=0:
a,b = b%a,a
return b
求a模m的逆:
def exgcd(a,m):
if gcd(a,m)!=1:
return None
u1,u2,u3 = 1,0,a
v1,v2,v3 = 0,1,m
while v3!=0:
q = u3//v3
v1,v2,v3,u1,u2,u3 = (u1-q*v1),(u2-q*v2),(u3-q*v3),v1,v2,v3
return u1%m
乘法:
C=a*b
加法:
C=a+b
轉置:
a_tra=a.T