1. 程式人生 > >python中矩陣的基本運算學習記錄

python中矩陣的基本運算學習記錄

 矩陣運算:

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