不同维度矩阵相乘[通俗易懂]

不同维度矩阵相乘[通俗易懂]在深度学习中经常会遇到不同维度的矩阵相乘的情况,本文会通过一些例子来展示不同维度矩阵乘法的过程。总体原则:在高维矩阵中取与低维矩阵相同维度的分片来与低维矩阵相乘,结果再按分片时的顺序还原为高维矩阵。相乘结果的维度与原来的高维矩阵一致。二维乘一维三维乘一维三维乘二维…

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全家桶1年46,售后保障稳定

前言

在深度学习中经常会遇到不同维度的矩阵相乘的情况,本文会通过一些例子来展示不同维度矩阵乘法的过程。

总体原则:在高维矩阵中取与低维矩阵相同维度的子矩阵来与低维矩阵相乘,结果再按子矩阵的排列顺序还原为高维矩阵。相乘结果的维度与原来的高维矩阵一致。
具体来说,当一方为一维矩阵时,另一方取其最后一维子矩阵来做乘法;当两方都是大于等于2维的矩阵时,取各自的最后两维构成的子矩阵来做乘法,其他维度体现结果的拼接信息,不参与运算(为batch训练提供了便利,batch中各样本的顺序在矩阵运算前后保持一致)。

实例:下面我们从低维到高维,依次演示不同维度矩阵相乘的结果。

二维乘一维

二维矩阵依次取出一维的行向量与一维矩阵做内积

#二维乘一维
import numpy as np

a = np.linspace(1,4,4).reshape(2,2)
b = np.array([1,1])
c = np.matmul(a,b)
print('a:\n',a)
print('b:\n',b)
print('ab:\n',c)

Jetbrains全家桶1年46,售后保障稳定

在这里插入图片描述

三维乘一维

三维矩阵包含两个二维矩阵,分别将这两个二维矩阵与一维矩阵相乘(乘积为一维),结果按原来的顺序拼接起来,构成一个二维矩阵

#三维乘一维
import numpy as np

a = np.linspace(1,8,8).reshape(2,2,2)
b = np.array([1,1])
c = np.matmul(a,b)
print('a:\n',a)
print('b:\n',b)
print('ab:\n',c)

在这里插入图片描述

二维乘二维

最常见的矩阵相乘形式

#二维乘二维
import numpy as np

a = np.linspace(1,4,4).reshape(2,2)
b = np.ones((2,2))
c = np.matmul(a,b)
print('a:\n',a)
print('b:\n',b)
print('ab:\n',c)

在这里插入图片描述

三维乘二维

将三维矩阵中的后两维组成的二维子矩阵分别与二维矩阵相乘(二维),结果再按原顺序拼接起来(三维)

#相当于三维矩阵里的二维分量分别与二维矩阵相乘,再拼接起来
import numpy as np

a=np.linspace(1,8,8).reshape(2,2,2)
# print(a)
b = np.array([[1,0],[0,1]]) #单位矩阵
c = np.matmul(a,b)
print('a:\n',a)
print('b:\n',b)
print('ab:\n',c)

在这里插入图片描述

三维乘三维

两个三维矩阵中对应位置的二维子矩阵分别相乘,结果按第0维分量更多的那个矩阵的结构拼接。
注意:,并不是任意两个三维矩阵都能相乘,其必须满足两个条件:

1:两个矩阵的后两个维度构成的二维矩阵之间必须满足二维矩阵相乘的条件,即第一个矩阵的列数等于第二个矩阵的行数
2:两个矩阵的第0维分量数必须相等(每个分量对应相乘) 或 有一方为1(broadcast-广播机制) —-反例见下方第3种情况

1-(2,2,2)*(2,2,2)

#三维乘三维 (2,2,2)*(2,2,2)
#计算时都是二维乘二维,第三维度反映二维矩阵的拼接信息;对应位置二维矩阵相乘
import numpy as np

a=np.linspace(1,8,8).reshape(2,2,2)
# print(a)
e = np.array([[[1,0],[0,1]]])
f = np.array([[[0,1],[1,0]]])
b = np.vstack((e,f))
c = np.matmul(a,b)
print('a:\n',a)
print('b:\n',b)
print('ab:\n',c)

在这里插入图片描述

2-(2,2,2)*(1,2,2)

#三维乘三维 (2,2,2)*(1,2,2)
#广播机制(broadcast)
import numpy as np

a=np.linspace(1,8,8).reshape(2,2,2)
# print(a)
b = np.array([[[0,1],[1,0]]])
c = np.matmul(a,b)
print('a:\n',a)
print('b:\n',b)
print('ab:\n',c)

在这里插入图片描述

3-(3,2,2)*(2,2,2)–失败

第0维分量数不满足条件2,不能相乘

# 三维乘三维--不同形状:(4,2,2)*(2,2,2)
import numpy as np

a=np.linspace(1,16,16).reshape(4,2,2)
# print(a)
e = np.array([[[1,0],[0,1]]])
f = np.array([[[0,1],[1,0]]])
b = np.vstack((e,f))
c = np.matmul(a,b)
print('a:\n',a)
print('b:\n',b)
print('c:\n',c)

在这里插入图片描述

多维乘多维

与三维乘三维类似,可乘条件2改为:除最后两维外,每一维的分量数必须对应相等(每个分量对应相乘) 或 有一方为1(broadcast-广播机制)

#各维度的分量相互对应,最终仍是计算二维乘二维
import numpy as np

a=np.linspace(1,16,16).reshape(2,2,2,2)
b = np.ones((2,2,2,2)) #全1矩阵
c = np.matmul(a,b)
print('a:\n',a)
print('b:\n',b)
print('ab:\n',c)

在这里插入图片描述

参考资料

知乎:多维矩阵相乘的可视化

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/230764.html原文链接:https://javaforall.net

(0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • iOS中什么是superView?(新手概念简述版)[通俗易懂]

    iOS中什么是superView?(新手概念简述版)[通俗易懂]文/司马陶冶(简书作者)/文章有删改原文链接:http://www.jianshu.com/p/8537dac37852著作权归作者所有,转载请联系作者获得授权,并标注“简书作者”。在iPhone

    2022年8月1日
    4
  • js垃圾处理机制_java中垃圾回收有什么目的

    js垃圾处理机制_java中垃圾回收有什么目的文章目录前置知识堆栈栈堆执行上下文与作用域链执行上下文作用域链一、JavaScript中怎么被定义为垃圾使用局部变量使用对象概括二、两种回收策略标记清理引用计数概括三、什么时候执行垃圾回收关于ChromeV8引擎的GC分代回收指针与活跃对象的区分回收的执行周期四、内存问题五、Es6WeakMap参考文章前置知识堆栈栈什么是栈栈其实是一种数据结构,有着先进后出,后进先出的特性,用生活中的事物来理解最形象的就是汉诺塔了。我们在栈中存储的数据就像汉诺塔的盘子一样,最先放进去在最下面,最后放入的盘.

    2022年10月9日
    0
  • Pytest(1)安装与入门「建议收藏」

    Pytest(1)安装与入门「建议收藏」pytest介绍pytest是python的一种单元测试框架,与python自带的unittest测试框架类似,但是比unittest框架使用起来更简洁,效率更高。根据pytest的官方网站介绍,它

    2022年7月31日
    4
  • soapUI工具使用方法、简介、接口测试[通俗易懂]

    soapUI工具使用方法、简介、接口测试[通俗易懂]一、soapUI简介SOAP:   WebService通过Http协议发送请求和接收结果时,发送的请求内容和结果内容都采用XML格式封装,并增加了一些特定的HTTP消息头,以说明HTTP消息头的内容格式,这些特定的HTTP消息头和XML内容格式就是SOAP协议。SOAP提供了标准的RPC方法来调用WebService。 SOAP协议 =HTTP协议 +XML数

    2022年6月17日
    45
  • java主要用来开发什么_浅谈Java开发主要是用来做什么的?

    java主要用来开发什么_浅谈Java开发主要是用来做什么的?java是一种编程语言,所以当然是用来“创造”的啦。它能做的事情非常多,涵盖了编程领域的方方面面。1、桌面级应用,简单来说就是主要功能都在计算机上运行的程序,比如word、excel等办公软件就归为此类。2、企业级应用(应用最广泛),简单来说就是,使用人数较多,数据量较大的应用。它对系统的稳定性、安全性、扩展性等要求较高。包括HR管理软件、营销管理软件、等几乎所有你能想到的应用。3、嵌入式设备及…

    2022年7月7日
    19
  • 最详细eclipse汉化插件安装教程

    最详细eclipse汉化插件安装教程教程作者:阿良,欢迎转载,转载请说明出处!本人QQ:583393588,欢迎广大安卓朋友交流!本文与《最详细安卓ADT装教程》是兄弟篇,欢迎查阅:http://blog.csdn.net/dai_zhenliang/article/details/8490837我鄙视那些中文不好而鄙视我用中文的人离线安装中文包方式一:使用下载单独的语言包比如汉化eclipseindig

    2022年6月1日
    32

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号