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