Python实现矩阵加法和乘法的方法分析

Python实现矩阵加法和乘法的方法分析

本文实例讲述了Python实现矩阵加法和乘法的方孙毕法。分享给大家供大家参考，具体如下：

本来以为python的矩阵用list表示出来应该很简单可以搞洞腊。。其实发现有大学问。

这里贴出我写的特别不pythonic的矩阵加法，作为反例。

def add(a, b):

rows = len(a[0])

cols = len(a)

c = []

for i in range(rows):

temp = []

for j in range(cols):

temp.append(a[i][j] + b[i][j])

c.append(temp)

return c

然后搜索纳凯滑了一下资料，果断有个很棒的，不过不知道有没有更棒的。

矩阵加法

def madd(M1, M2):

if isinstance(M1, (tuple, list)) and isinstance(M2, (tuple, list)):

return [[m+n for m,n in zip(i,j)] for i, j in zip(M1,M2)]

矩阵乘法

def multi(M1, M2):

if isinstance(M1, (float, int)) and isinstance(M2, (tuple, list)):

return [[M1*i for i in j] for j in M2]

if isinstance(M1, (tuple, list)) and isinstance(M2, (tuple, list)):

return [[sum(map(lambda x: x[0]*x[1], zip(i,j)))

for j in zip(*M2)] for i in M1]

1.numpy的导入和使用

data1=mat(zeros((

)))

#创建一个3*3的零矩阵，矩阵这里zeros函数的参数是一个tuple类型(3,3)

data2=mat(ones((

)))

#创建一个2*4的1矩阵，默认是浮点型的数据，如果需要时int类型，可以使用dtype=int

data3=mat(random.rand(

))

#这里拿派弊的random模块使用的是消族numpy中的random模块，random.rand(2,2)创建的是一个二维数组，需要将其转换成#matrix

data4=mat(random.randint(

10

,size=(

)))

#生成一个3*3的0-10之间的随机整数矩阵，如果需要指定下界则羡粗可以多加一个参数

data5=mat(random.randint(

,size=(

))

#产生一个2-8之间的随机整数矩阵

data6=mat(eye(

,dtype=

int

))

#产生一个2*2的对角矩阵

a1=[

]a2=mat(diag(a1))

#生成一个对角线为1、2、3的对角矩阵

---