NumPy:
1、NumPy 是一个功能强大的第三方库(需要自己安装),主要用于对多维数组执行计算;
它提供了大量的库函数和操作,可以帮助程序员更轻松地进行数值计算
2、可以和另外两个第三方库 SciPy 和 Matplotlib 一起使用从而在一定程度上替换对 Matlab 的使用
3、主要应用:
①数学运算:NumPy 对于执行各种数学运算非常有用,如数值积分、微分、内插、外推以及矩阵的内积、外积、特征向量等。
②图像处理和计算机图形学:计算机中的图像表示为多维数字数组。NumPy 成为同样情况下最自然的选择。
实际上,NumPy 提供了一些优秀的库函数来快速处理图像。例如,镜像图像、按特定角度旋转图像等
③机器学习模型:在编写机器学习算法时,需要对矩阵进行各种数值计算。例如矩阵乘法、换位、加法等
NumPy 数组:
1、NumPy 提供最重要的数据结构是一个称为 ndarray 的 N 维数组类型。 它描述相同类型的元素集合,NumPy 数组是通常的 Python 数组的扩展
ndarray 配备了大量的函数和运算符,可以帮助我们快速编写各种类型计算的高性能代码,每个元素在内存中使用相同大小的块
2、NumPy 数组的创建方法:
①从其他 python 数据类型(如:列表、元组等)转换过来
②NumPy 原生数组的创建(通过 arange、ones、zeros 等创建)
③使用特殊库函数(例如,random)创建
④从磁盘读取数组,无论是标准格式还是自定义格式
⑤通过使用字符串或缓冲区从原始字节创建数组
import numpy as np lst = [1, 2, 3] arr = np.array(lst) # 输出 Type: <class 'numpy.ndarray'> # arr: [1 2 3] print("Type: ", type(arr), " arr: ", arr) # numpy.arange(star, end, step) 返回 [star, end) 隔 step-1 的各个元素 # 输出 [1 3 5 7 9] arr = np.arange(1, 10, 2) print("arr", arr, sep=" = ") # numpy.random.randn(size) 创建一个长度为 size 的服从标准正态分布的随机数组 arr = np.random.randn(5) print("arr: ", arr) # numpy.random.random(size) 创建一个范围为 [0, 1),长度为 size 的服从均匀分布的随机数组 arr = np.random.random(5) print("arr: ", arr) # numpy.linspace(star, end, size) 创建一个长度为 size,范围为 [star, end] 的数组 # 输出 [0. 2.5 5. 7.5 10.] arr = np.linspace(0, 10, 5) print("arr: ", arr) # numpy.zeros(shape) 创建一个用 0 填充,形状为 shape 的数组 a = np.zeros((2, 2)) print("a", a, sep=" = ") b = np.ones((2, 2)) print("b", b, sep=" = ") c = np.full((2, 2), 7) # 创建一个 2*2 矩阵,并填充 7 print("c", c, sep=" = ") d = np.eye(3) # 创建一个 3*3 的单位矩阵 print("d", d, sep=" = ")
注:np.array(object, dtype = None, copy = True, order = None, subok = False, ndmin = 0) 里,
object 表示任何暴露数组接口方法的对象都会返回一个数组或任何(嵌套)序列;
dtype 为数组的数据类型,可选;copy 可选,默认为true,表示对象是否被复制;order 里 C(按行)、F(按列)或A(任意,默认);
subok 默认情况下,返回的数组被强制为基类数组。 如果为true,则返回子类;ndmin 指定返回数组的最小维数
import numpy as np # 本来是一个一维数组,但通过 ndmin=2 使得数组最小维度为2维 # 输出 arr: [[1 2 3 4]] 注意为 2 维数组 arr = np.array([1, 2, 3, 4], ndmin=2) print("arr", arr, sep=": ") # 输出 arr: [[1 2] [3 4]],因为其维度大于最小维度 arr = np.array([[1, 2], [3, 4]], ndmin=1) print("arr", arr, sep=": ") # 指定数据类型为复数 arr = np.array([1, 2, 3], dtype=complex) # 输出 arr: [1.+0.j 2.+0.j 3.+0.j] print("arr", arr, sep=": ")
3、数组属性:
①shape:返回数组形状,如 (2, 3) 表示 2 行 3 列的 2 维数组
②ndim:返回数组维度
③size:返回数组里元素个数
④itemsize:返回数组中每个元素的字节单位长度
⑤T:转置矩阵,但不会改变原矩阵
⑥flags:返回对象的内存信息
⑦real:返回元素的实部
⑧imag:返回元素的虚部
⑨data:包含实际数组元素的缓冲区,由于一般通过数组的索引获取元素,所以通常不需要使用这个属性
import numpy as np # shape 这一数组属性返回一个包含数组维度的元组,它也可以用于调整数组大小 arr = np.array([[1, 2, 3], [4, 5, 6]]) arr_t = arr[:, 2] # 输出列号为 2 的所有元素,即 [3, 6] print("arr_t[][2]: ", arr_t) # 输出数组形状,这里为 (2, 3) print("The arr's shape: ", arr.shape) arr.shape = (3, 2) # 修改形状为 3 行 2 列,注意形状和元素个数 print("The new arr: ", arr) # 也可以通过 reshape() 函数来调整大小,但不会改变原数组 temp = arr.reshape(2, 3) print(temp) # 2 维数组,且数组里有 3 个元素 # ndim 属性返回数组维度 arr = np.arange(24) # 创建 [0...23] 的数组 print(arr.ndim) # 输出 1 arr.shape = (2, 4, 3) # 修改为 3 维数组 print("arr's ndim: ", arr.ndim) # 输出 3 # itemsize 这一数组属性返回数组中每个元素的字节单位长度 arr = np.array([1, 2, 3, 4, 5]) print("size: ", arr.size) print("itemsize: ", arr.itemsize) # 输出 4(整数默认为 4 字节) # dtype 返回数组元素类型 print("dtype: ", arr.dtype) # 输出 dtype: int32(即 4 字节) arr.dtype = np.int8 # 修改整数类型 print("itemsize: ", arr.itemsize) # 输出 1 arr = np.array([[1, 2, 3], [4, 5, 6]]) # 调用矩阵的 T 属性,可以转置这个矩阵,但不会改变原矩阵 print("translate arr: ", arr.T) print("arr: ", arr) # 输出元素的实虚部 arr = np.array([1+2j, 2+3j, 3+4j]) print("arr's real: ", arr.real) # 输出 [1. 2. 3.] print("arr's imag: ", arr.imag) # 输出 [2. 3. 4.] # 输出对象的内存信息 print(arr.flags)
4、数组的简单运算:
大部分数学运算均只对在相应元素进行,如运算四则运算加减乘除,并且运算支持广播
(即不同形状的数组如果符合某种条件则可以进行运算)!
若要进行矩阵乘法则需要函数 dot() 或运算符 @
import numpy as np arr_one = np.array([[1.0, 2.0], [3.0, 4.0]]) arr_two = np.array([[5.0, 6.0], [7.0, 8.0]]) # 四则运算时对应元素进行运算 sum = arr_one + arr_two dif = arr_one - arr_two mul = arr_one * arr_two quo = arr_one / arr_two matrix_product = arr_one.dot(arr_two) # 矩阵乘法 print(("sum: ", sum), ("dif: ", dif), ("mul: ", mul), ("quo: ", quo), sep=" ") print("matrix_product: ", matrix_product) # 大小比较也是通过对应元素进行的 print("arr_one < arr_two = ", arr_one < arr_two) print("arr_one ** 2", arr_one ** 2, sep=" = ") # @ 可以用于矩阵乘法 print("matrix_product: ", arr_one @ arr_two)