NumPy之矩阵向量线性代数等操作示例

 更新时间:2023年04月28日 14:26:01   作者:CodeDevMaster  
这篇文章主要为大家介绍了NumPy之矩阵向量线性代数等操作示例详解,有需要的朋友可以借鉴参考下,希望能够有所帮助,祝大家多多进步,早日升职加薪

NumPy矩阵和向量

矩阵

在NumPy中,矩阵可以看作是一个二维数组,其中每个元素都可以通过行列坐标来定位。它表示为一个m×n的矩形网格,其中m表示矩阵的行数,n表示矩阵的列数。在计算机科学中,矩阵通常用数字或符号表示,并且可以进行加、减、乘等运算。

一个M X N的矩阵是一个由M行(row)N列(column)元素排列成的矩形阵列。矩阵里的元素可以是数字、符号或数学式。

以下是一个由 6 个数字元素构成的 2 行 3 列的矩阵:

[1 2 3]
[4 5 6]

注意:2×3矩阵即2行3列,请把左、右两边的多个[与多个]看成一个整体[]

向量

在 NumPy 中,向量是一维数组对象,其所有元素都必须具有相同的数据类型。向量可以通过创建一个一维数组来实现。

向量是一种特殊的矩阵,其中只包含一行或一列元素。向量通常用数字或符号表示,其大小表示向量包含的元素数量。

创建一个包含五个元素的向量:

import numpy as np
v = np.array([1, 2, 3, 4, 5])

创建向量

import numpy as np
# 创建一维数组
v = np.array([1, 2, 3])
print(v)  # 输出:[1 2 3]

创建矩阵

import numpy as np
# 创建二维数组
m = np.array([[1, 2, 3], [4, 5, 6]])
print(m)  # 输出:[[1 2 3]
          #      [4 5 6]]

访问元素

import numpy as np
# 创建一维数组
v = np.array([1, 2, 3])
# 访问元素
print(v[0])  # 输出:1
print(v[1])  # 输出:2
print(v[2])  # 输出:3

转置矩阵

NumPy中除了可以使用numpy.transpose 函数来对换数组的维度,还可以使用 T 属性。

例如有个 m 行 n 列的矩阵,使用 t() 函数就能转换为 n 行 m 列的矩阵。

import numpy as np
# 创建矩阵
A = np.array([[1, 2], [3, 4]])
# 转置矩阵
AT = np.transpose(A)
print(AT)  # 输出:[[1 3]
           #      [2 4]]
print (A.T)    

矩阵加减乘除

# 矩阵加法
np.add(A, B)
# 矩阵减法
np.subtract(A, B)
# 矩阵乘法
# 在进行矩阵乘法时,前一个矩阵的列数必须等于后一个矩阵的行数,才能进行乘法运算
# (M行, N列)*(N行, L列) = (M行, L列)
p.dot(A, B)
# 矩阵除法
np.divide(A, B)
import numpy as np
# 创建矩阵A和B
A = np.array([[1, 2], [3, 4]])
B = np.array([[5, 6], [7, 8]])
print("A + B:")
print(np.add(A, B))
[1, 2]  	[5, 6]		[6, 8]
		+  			=  
[3, 4]		[7, 8]		[10, 12]
print("A - B:")
print(np.subtract(A, B))
[1, 2]  	[5, 6]		[-4, -4]
		-  			=  
[3, 4]		[7, 8]		[-4, -4]
print("A * B:")
print(np.dot(A, B))
[1, 2]  	[5, 6]		[19, 22]
		*  			=  
[3, 4]		[7, 8]		 [43, 50]
print("A / B:")
print(np.divide(A, B))
[1, 2]  	[5, 6]		[0.2,  0.33333333]
		/  			=  
[3, 4]		[7, 8]		[0.42857143, 0.5]

矩阵和矩阵(向量)相乘: (M行, N列)*(N行, L列) = (M行, L列)

矩阵向量乘法

m×n 的矩阵乘以 n×1 的向量,得到的是 m×1 的向量

[1, 2]  	[1]			[19]
		*  			=  
[3, 4]		[1]			 [43]

矩阵求逆

使用numpy.linalg.inv()函数进行矩阵求逆操作

import numpy as np
# 创建矩阵
matrix = np.array([[1, 2], [3, 4]])
# 求逆矩阵
result = np.linalg.inv(matrix)
print(result)

矩阵的迹

使用numpy.trace()函数可以计算矩阵的迹

import numpy as np
# 创建矩阵
matrix = np.array([[1, 2], [3, 4]])
# 计算矩阵的迹
result = np.trace(matrix)
print(result)

向量点积

使用numpy.dot()函数进行向量点积操作

import numpy as np
# 创建两个向量
vector1 = np.array([1, 2])
vector2 = np.array([3, 4])
# 向量点积
result = np.dot(vector1, vector2)
print(result)

向量范数

使用numpy.linalg.norm()函数可以计算向量的范数

import numpy as np
# 创建向量
vector = np.array([1, 2, 3])
# 计算向量的L2范数
result = np.linalg.norm(vector)
print(result)

NumPy线性代数

NumPy的线性代数模块(numpy.linalg)提供了许多矩阵运算函数,如矩阵乘法、求逆、行列式、特征值等,该库包含了线性代数所需的所有功能。

常用的 NumPy 线性代数函数:

函数说明
np.dot(a, b):两个数组的点积,即元素对应相乘
np.matmul(a, b)两个数组的矩阵积
np.linalg.inv(a)计算矩阵的逆
np.linalg.det(a)计算矩阵的行列式
np.linalg.eig(a)计算矩阵的特征值和特征向量
np.linalg.solve(a, b)解线性方程组 ax=b

计算矩阵乘积

import numpy as np
a = np.array([[1, 2], [3, 4]])
b = np.array([[5, 6], [7, 8]])
print(np.dot(a, b))
# [[19 22]
#  [43 50]]

计算矩阵的逆

import numpy as np
a = np.array([[1, 2], [3, 4]])
print(np.linalg.inv(a))
# [[-2.   1. ]
#  [ 1.5 -0.5]]

解线性方程组

import numpy as np
a = np.array([[2, 3], [4, 5]])
b = np.array([5, 6])
x = np.linalg.solve(a, b)
print(x) # [-4.  5.]

解得的线性方程组如下,其解为x1=−4,x2=5

以上就是NumPy之矩阵向量线性代数等操作示例的详细内容,更多关于NumPy 矩阵向量线性代数的资料请关注脚本之家其它相关文章!

相关文章

  • Python 继承,重写,super()调用父类方法操作示例

    Python 继承,重写,super()调用父类方法操作示例

    这篇文章主要介绍了Python 继承,重写,super()调用父类方法,结合完整实例形式详细分析了Python面向对象程序设计中子类继承与重写父类方法的相关操作技巧,需要的朋友可以参考下
    2019-09-09
  • python爬虫 猫眼电影和电影天堂数据csv和mysql存储过程解析

    python爬虫 猫眼电影和电影天堂数据csv和mysql存储过程解析

    这篇文章主要介绍了python爬虫 猫眼电影和电影天堂数据csv和mysql存储过程解析,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下
    2019-09-09
  • Python字节串"b"前缀案例解析

    Python字节串"b"前缀案例解析

    文章主要介绍了Python中的字节串(bytes)类型,包括其定义、特性、应用场景、语法规则、类型转换方法以及实际应用案例,本文结合实例代码给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友参考下吧
    2026-02-02
  • tensorflow实现读取模型中保存的值 tf.train.NewCheckpointReader

    tensorflow实现读取模型中保存的值 tf.train.NewCheckpointReader

    今天小编就为大家分享一篇tensorflow实现读取模型中保存的值 tf.train.NewCheckpointReader,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧
    2020-02-02
  • Pandas 中的join函数应用实现删除多余的空行

    Pandas 中的join函数应用实现删除多余的空行

    这篇文章主要介绍了Pandas 中的join函数应用实现删除多余的空行,str.join也就是sequence要连接的元素序列,下面我们来看看他的作用实现删除多余的空行,需要的小伙伴可以参考一下
    2022-02-02
  • Python之——生成动态路由轨迹图的实例

    Python之——生成动态路由轨迹图的实例

    今天小编就为大家分享一篇Python之——生成动态路由轨迹图的实例,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧
    2019-11-11
  • Python pygame绘制文字制作滚动文字过程解析

    Python pygame绘制文字制作滚动文字过程解析

    这篇文章主要介绍了Python pygame绘制文字制作滚动文字过程解析,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下
    2019-12-12
  • Python中作用域的规则与应用

    Python中作用域的规则与应用

    在 Python 编程中,作用域(Scope) 是指一个变量可以被访问和引用的范围,这篇文章主要为大家介绍了Python中作用域的规则与应用的相关知识,需要的可以了解下
    2025-01-01
  • python二维图制作的实例代码

    python二维图制作的实例代码

    这篇文章主要给大家介绍了关于python二维图制作的相关资料,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧
    2020-12-12
  • python中opencv Canny边缘检测

    python中opencv Canny边缘检测

    这篇文章主要介绍了python中opencv Canny边缘检测,Canny边缘检测是一种使用多级边缘检测算法检测边缘的方法。OpenCV提供了函数cv2.Canny()实现Canny边缘检测。更多相关内容需要的小伙伴可以参考下面文章内容
    2022-06-06

最新评论