python3.7中矩阵相乘np.dot和*的区别

作为一名菜鸟级初学者,在进行python numpy模块练习的时候发现

import numpy as np
from numpy import *

A=[[8,1,6],[3,5,7],[4,9,2]]
aa=linalg.inv(A)*A
print (aa)

运行结果:

[[ 1.17777778 -0.14444444  0.38333333]
 [-0.18333333  0.11111111  0.73888889]
 [-0.07777778  1.7        -0.20555556]]

当时觉得世界崩塌了,这种计算的结果难道不是单位阵吗?但是博主坚信电脑是不会错的,只可能是自已的问题。为了发现问题

博主做了print (linalg.inv(A))的操作
发现结果和matlab算的一样,居然是对的。

[[ 0.14722222 -0.14444444  0.06388889]
 [-0.06111111  0.02222222  0.10555556]
 [-0.01944444  0.18888889 -0.10277778]]

这就很奇怪了。当时觉得这个编绎器可能有bug,于是换了另一个,但是结果还是一样的。
再后来,想到了python里还有np.dot的运算,于是将程序改为

import numpy as np
from numpy import *

A=[[8,1,6],[3,5,7],[4,9,2]]
aa=linalg.inv(A)*A
print (aa)

bb=np.dot(linalg.inv(A),A)
print(bb)

运行结果为

[[ 1.17777778 -0.14444444  0.38333333]
 [-0.18333333  0.11111111  0.73888889]
 [-0.07777778  1.7        -0.20555556]]
[[ 1.00000000e+00  0.00000000e+00  5.55111512e-17]
 [-5.55111512e-17  1.00000000e+00 -5.55111512e-17]
 [ 5.55111512e-17  9.71445147e-17  1.00000000e+00]]

可以看出np.dot计算的是传统矩阵相乘,而*号运算计算的不知道是什么。。。
但是python能这样计算肯定有它的道理,经过仔细观察后,博主发现前边的结果是矩阵对应位置处两个数相乘的结果,也就是matlab中的点乘操作,而后者才是真正的矩阵相乘的结果。
另外,经过验证发现如果在矩阵定义时采用np.mat函数,那么也可以用*号实现矩阵的相乘
程序修改为:

import numpy as np
from numpy import *

A=np.mat([[8,1,6],[3,5,7],[4,9,2]])
aa=linalg.inv(A)*A
print (aa)

bb=np.dot(linalg.inv(A),A)
print(bb)

运行结果为:

[[ 1.00000000e+00  0.00000000e+00  5.55111512e-17]
 [-5.55111512e-17  1.00000000e+00 -5.55111512e-17]
 [ 5.55111512e-17  9.71445147e-17  1.00000000e+00]]
[[ 1.00000000e+00  0.00000000e+00  5.55111512e-17]
 [-5.55111512e-17  1.00000000e+00 -5.55111512e-17]
 [ 5.55111512e-17  9.71445147e-17  1.00000000e+00]]

可以看出,np.dot和*号计算结果相同。因此在使用numpy进行矩阵计算时一定要注意二者的区别(一个菜鸟的忠告)。

  • 6
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值