【技术积累】Python中的NumPy库【二】

特务飞鹰

NumPy库的主要类有哪些？

NumPy库的主要类包括：

• ndarray：N维数组对象，是NumPy最重要的类之一。它是Python中数组的基本数据结构，可以进行高效的数学计算和数据处理操作。
  
• ufunc：通用函数对象，是NumPy库中的另一个重要类。它是一种高效的元素级运算工具，提供了基本的数学运算、逻辑运算和位运算等基本的数学和逻辑操作。
  
• linspace：生成一定数量的等差数列，返回一个一维数组。
  
• meshgrid：用于生成二维的坐标矩阵，常用于三维绘图和计算机图形学等领域。
  
• random：随机数生成器，提供了众多生成随机数和随机样本的函数，常用于模型训练和数据分析。
  
• matlib：矩阵库，提供了各种矩阵操作和矩阵运算的函数，包括矩阵乘法、矩阵分解和矩阵求逆等。
  
• fft：快速傅里叶变换（FFT）类，提供了高效的傅里叶变换算法，常用于信号处理和图像处理等领域。
  
• poly：多项式类，提供了多项式求解、多项式拟合和多项式积分等数学运算。常用于数据拟合和模型建立。
  
• linalg：线性代数类，提供了各种线性代数运算的函数，包括矩阵求逆、特征值求解和奇异值分解等。常用于数据分析和机器学习等领域。
  

NumPy库和Python自带的列表有什么区别？

NumPy库是Python中用于科学计算和数学计算的常用库，它提供了多维数组对象和各种数学函数和操作。相比于Python自带的列表，NumPy数组具有以下几点优势：

• 内存占用少：NumPy数组中的所有元素都是相同的数据类型，因此内存占用更小，运行速度更快。
  
• 数组操作方便：NumPy数组可以进行各种数学运算和逻辑运算，例如矩阵乘法、数组切片、数组重塑等操作，这些操作可以在不使用循环的情况下完成。
  
• 广播功能：NumPy数组可以进行广播操作，即对形状不同的数组进行相应的操作，这大大减少了数组形状不同需要循环处理的情况。
  

而Python自带的列表则是一种基本的数据结构，可以存储任意的对象。相较于NumPy数组，其灵活性更高，可以进行各种操作，同时也更容易理解和使用。但列表的操作性能在大型数据集时会比NumPy数组慢。
如何将Python数列转换为NumPy数组？

将Python数列（即列表）转换为NumPy数组可以使用NumPy库中的array()函数。该函数接受一个列表参数，返回一个NumPy数组。
下面是一个示例代码：

1. import numpy as np
3. # 定义一个Python数列
4. lst = [1, 2, 3, 4, 5]
6. # 转换为NumPy数组
7. arr = np.array(lst)
9. # 输出结果
10. print(arr)

复制代码

输出结果为：

1. [1 2 3 4 5]

复制代码

在上面的代码中，我们首先导入了NumPy库。然后，我们定义了一个Python数列lst。接着，我们使用array()函数将lst转换为NumPy数组arr。最后，我们输出了arr的值。
注意，当我们将Python数列转换为NumPy数组时，NumPy会自动根据数列中的值的类型来推断生成的数组的数据类型。例如，上面的代码中lst是一个包含整数的数列，因此生成的数组也是整数类型的。如果数列中包含浮点数，则生成的数组将是浮点数类型的。
NumPy库中的数组是什么？有哪些特征？如何创建

NumPy库中的数组是一种多维数组对象，又称为ndarray。这些数组是用来存储相同数据类型的元素（例如int、float等），它是一个可变的对象，可以在其中进行快速且有效的数据操作。
以下是NumPy数组的一些特征：

• 数组是一种类似于列表的对象，但在NumPy中，它们可以包含更多的维度。
  
• 数组中的元素必须是相同的类型。
  
• NumPy数组的大小是固定的，也就是说，一旦创建，就无法再加入新元素。
  
• 数组可以通过下标访问，也可以使用NumPy提供的高效的函数对数组进行操作。
  

以下是创建NumPy数组的几种方式：
1. 通过列表创建

1. import numpy as np
3. a = np.array([1, 2, 3])  # 一维数组
4. b = np.array([[1, 2, 3], [4, 5, 6]])  # 二维数组

复制代码

2. 通过arange函数创建

1. import numpy as np
3. a = np.arange(10)  # 生成一维数组[0 1 2 3 4 5 6 7 8 9]
4. b = np.arange(0, 10, 2)  # 生成一维数组[0 2 4 6 8]，步长为2
5. c = np.arange(12).reshape(3, 4)  # 生成二维数组

复制代码

3. 使用zeros和ones函数创建

1. import numpy as np
3. a = np.zeros(5)  # 生成一维数组[0. 0. 0. 0. 0.]
4. b = np.zeros((2, 3))  # 生成二维数组
5. c = np.ones((2, 4, 3))  # 生成三维数组，全部元素为1

复制代码

4. 使用random函数创建

1. import numpy as np
3. a = np.random.rand(3)  # 生成一维数组，元素是[0, 1)之间的随机数
4. b = np.random.randn(3, 2)  # 生成二维数组，元素是标准正态分布的随机数
5. c = np.random.randint(1, 10, size=(2, 3))  # 生成二维数组，元素是1-10之间的随机整数

复制代码

 如何使用NumPy数组进行算术运算

1. import numpy as np
3. # 创建两个NumPy数组
4. a = np.array([1, 2, 3])
5. b = np.array([4, 5, 6])
7. # 加法
8. c = a + b
9. print("加法：", c)  # [5 7 9]
11. # 减法
12. c = a - b
13. print("减法：", c)  # [-3 -3 -3]
15. # 乘法
16. c = a * b
17. print("乘法：", c)  # [ 4 10 18]
19. # 除法
20. c = a / b
21. print("除法：", c)  # [0.25 0.4  0.5 ]
23. # 指数
24. c = np.power(a, 2)
25. print("指数：", c)  # [1 4 9]
27. # 平方根
28. c = np.sqrt(a)
29. print("平方根：", c)  # [1.         1.41421356 1.73205081]
31. # 取反
32. c = -a
33. print("取反：", c)  # [-1 -2 -3]

复制代码

如何使用NumPy数组进行三角函数运算？

1. import numpy as np
3. # 创建一个长度为10的一维NumPy数组
4. arr = np.array([0, np.pi/6, np.pi/4, np.pi/3, np.pi/2, np.pi, 3*np.pi/2, 2*np.pi/3, 3*np.pi/4, 5*np.pi/6])
6. # 计算数组中每个元素的正弦值
7. sin_arr = np.sin(arr)
8. print("sin函数结果：", sin_arr)
10. # 计算数组中每个元素的余弦值
11. cos_arr = np.cos(arr)
12. print("cos函数结果：", cos_arr)
14. # 计算数组中每个元素的正切值
15. tan_arr = np.tan(arr)
16. print("tan函数结果：", tan_arr)

复制代码

输出结果：

1. sin函数结果： [ 0.          0.5         0.70710678  0.8660254   1.          0.
2. -1.         -0.8660254  -0.70710678 -0.5       ]
4. cos函数结果： [ 1.00000000e+00  8.66025404e-01  7.07106781e-01  5.00000000e-01
5.   6.12323400e-17 -1.00000000e+00 -1.83697020e-16 -5.00000000e-01
6. -7.07106781e-01 -8.66025404e-01]
8. tan函数结果： [ 0.00000000e+00  5.77350269e-01  1.00000000e+00  1.73205081e+00
9.   1.63312394e+16 -1.22464680e-16  5.44374645e+15  1.73205081e+00
10.   1.00000000e+00  5.77350269e-01]

复制代码

如何使用NumPy数组进行指数和对数运算？

NumPy数组可以使用exponential（指数）和logarithmic（对数）函数进行指数和对数运算。
exponential函数计算每个元素的指数值，logarithmic函数计算每个元素的自然对数，也可以计算底数为其他常数的对数。
下面是一个简单的例子，演示如何使用NumPy数组进行指数和对数操作。

1. import numpy as np
3. arr = np.array([1, 2, 3, 4, 5])
5. # 指数运算
6. print("exponential:", np.exp(arr))
8. # 自然对数
9. print("natural logarithm:", np.log(arr))
11. # 底数为2的对数
12. print("log base 2:", np.log2(arr))
14. # 底数为10的对数
15. print("log base 10:", np.log10(arr))

复制代码

在上面的例子中，我们创建了一个一维NumPy数组，并使用numpy.exp函数计算每个元素的指数值，使用numpy.log函数计算每个元素的自然对数，使用numpy.log2函数计算底数为2的对数，使用numpy.log10函数计算底数为10的对数。
输出结果如下：

1. exponential: [   2.71828183    7.3890561    20.08553692   54.59815003  148.4131591 ]
2. natural logarithm: [ 0.          0.69314718  1.09861229  1.38629436  1.60943791]
3. log base 2: [ 0.          1.          1.5849625   2.          2.32192809]
4. log base 10:[ 0.          0.30103     0.47712125  0.60205999  0.69897   ]

复制代码

因此，使用NumPy进行指数和对数计算非常简单，只需调用相应的函数即可。
如何使用NumPy数组进行统计运算

NumPy是Python中用于数值计算和科学计算的常用库。它提供了一个称为ndarray的多维数组对象，可以使用它进行统计运算。
以下是一个使用NumPy数组进行统计运算的例子：

1. import numpy as np
3. # 创建一个随机数组
4. arr = np.random.randn(100)
6. # 计算数组的均值、中位数、标准差和方差
7. mean = np.mean(arr)
8. median = np.median(arr)
9. std_dev = np.std(arr)
10. variance = np.var(arr)
12. print('均值：', mean)
13. print('中位数：', median)
14. print('标准差：', std_dev)
15. print('方差：', variance)

复制代码

上面的代码创建了一个包含100个随机数的NumPy数组，然后对数组进行了一些统计运算。np.mean()函数计算数组的均值，np.median()函数计算中位数，np.std()函数计算标准差，np.var()函数计算方差。
输出结果可能会因为随机数的不同而有所不同，但通常会接近以下值：

1. 均值： 0.03821741789938476
2. 中位数： 0.04451500104648571
3. 标准差： 0.9616266701218458
4. 方差： 0.9241266652394576

复制代码

如何使用NumPy数组进行线性代数运算？

NumPy是基于Python的库，主要用于科学计算。NumPy中包含大量的数学函数，以支持矩阵和数组的运算。在NumPy中，数组对象被称为ndarray，可以进行各种数学运算。下面是使用NumPy数组进行线性代数运算的案例：

1. import numpy as np
3. # 定义 NumPy 数组
4. A = np.array([[1, 2, 3], [4, 5, 6], [7, 8, 9]])
5. B = np.array([1, 2, 3])
7. # 计算 A 的逆矩阵
8. A_inv = np.linalg.inv(A)
10. # 求解线性方程组 Ax = B
11. x = np.linalg.solve(A, B)
13. # 计算 A 和 B 的乘积
14. C = A.dot(B)
16. # 计算 A 和 B 的点积
17. D = np.dot(A, B)
19. # 计算矩阵 A 的行列式
20. det_A = np.linalg.det(A)
22. # 计算矩阵 A 的特征值和特征向量
23. eigenvalues, eigenvectors = np.linalg.eig(A)
25. # 打印结果
26. print('A:')
27. print(A)
28. print('A_inv:')
29. print(A_inv)
30. print('B:')
31. print(B)
32. print('x:')
33. print(x)
34. print('C:')
35. print(C)
36. print('D:')
37. print(D)
38. print('det_A:')
39. print(det_A)
40. print('eigenvalues:')
41. print(eigenvalues)
42. print('eigenvectors:')
43. print(eigenvectors)

复制代码

这是一个简单的NumPy数组线性代数运算的示例，包括计算逆矩阵、求解线性方程组、计算矩阵乘积、点积、矩阵行列式、特征值和特征向量等等。在NumPy中，还有很多其他的线性代数运算函数可供使用。
如何使用NumPy数组进行傅里叶变换

NumPy中的傅里叶变换函数位于numpy.fft模块中，其中包括傅里叶变换、逆变换和频率域滤波等功能。以下是使用NumPy数组进行傅里叶变换的示例：

1. import numpy as np
2. import matplotlib.pyplot as plt
4. # 生成信号
5. fs = 1000  # 采样率
6. t = np.arange(0, 1, 1/fs)  # 时间轴
7. f1, f2 = 10, 100  # 两个频率成分
8. s = np.sin(2*np.pi*f1*t) + 0.5*np.sin(2*np.pi*f2*t)
10. # 进行傅里叶变换
11. S = np.fft.fft(s)  # 频域信号
12. freqs = np.fft.fftfreq(len(s), 1/fs)  # 频率轴
14. # 绘制原始信号
15. plt.subplot(2, 1, 1)
16. plt.plot(t, s)
17. plt.xlabel('Time (s)')
18. plt.ylabel('Amplitude')
19. plt.title('Original Signal')
21. # 绘制频域信号
22. plt.subplot(2, 1, 2)
23. plt.plot(freqs, np.abs(S))
24. plt.xlabel('Frequency (Hz)')
25. plt.ylabel('Magnitude')
26. plt.title('Frequency Domain Signal')
28. plt.show()

复制代码

在此示例中，我们生成了一个包含两个频率成分的信号，并利用`np.fft.fft()`函数进行了傅里叶变换，得到了信号的频域表示。然后，利用`np.fft.fftfreq()`函数生成了频率轴，最后将时域信号和频域信号绘制在同一张图中进行比较。
需要注意的是，在进行傅里叶变换之前，需要保证信号的采样率是足够高的，以避免出现混淆和重叠的频率成分。此外，傅里叶变换生成的频域信号一般是对称的，因此通常只需要使用其一半进行后续处理和分析。
如何使用NumPy数组进行形状变换？

NumPy提供了多种方法用于数组形状变换。可以使用reshape函数重新构造数组的形状，也可以使用transpose函数交换数组的维度。
示例代码：

1. import numpy as np
3. # 创建一个包含1~9的一维数组
4. arr = np.arange(1, 10)
6. # 将一维数组重塑为3 X 3的二维数组
7. arr_reshaped = arr.reshape((3, 3))
9. # 输出二维数组
10. print("Reshaped array:\n", arr_reshaped)
12. # 交换二维数组的维度
13. arr_transposed = arr_reshaped.transpose()
15. # 输出交换维度后的数组
16. print("Transposed array:\n", arr_transposed)

复制代码

输出结果：

1. Reshaped array:
2. [[1 2 3]
3.   [4 5 6]
4.   [7 8 9]]
5.   
6. Transposed array:
7. [[1 4 7]
8.   [2 5 8]
9.   [3 6 9]]

复制代码

在上面的代码中，我们从1到9创建了一个一维数组，并使用reshape函数将其重新整形为3 X 3的二维数组。然后，我们使用transpose函数交换了数组的维度，并输出了结果。

来源:https://www.cnblogs.com/yyyyfly1/p/17479857.html
免责声明：由于采集信息均来自互联网，如果侵犯了您的权益，请联系我们【E-Mail:cb@itdo.tech】 我们会及时删除侵权内容，谢谢合作！