
第十三章 工具衔接与迁移
NumPy 作为 Python 科学计算的核心库,不仅功能强大,还与众多工具无缝衔接。本文从两个关键方向系统讲解:
MATLAB 用户迁移到 NumPy 时,最大的挑战是思维转换(从矩阵中心到数组中心)和语法差异。以下是核心对照表。
ndarray,不是 matrix → 所有运算默认是逐元素的@ 或 np.dot(),不是 *MATLAB | NumPy | 说明 |
|---|---|---|
|
| 行向量 |
|
| 列向量 |
|
| 注意:元组参数 |
|
| |
|
| 单位矩阵 |
|
| 相同 |
|
| 均匀分布 [0,1) |
|
| 展平为一维 |
💡 MATLAB 的
[1,2,3]是行向量,NumPy 的np.array([1,2,3])是一维数组(无行列之分),但在广播中可视为行或列。
MATLAB | NumPy | 说明 |
|---|---|---|
|
| 索引从 0 开始 |
|
| 第 2 行 |
|
| 第 3 列 |
|
| 从第 2 行到末尾 |
|
| 花式索引(注意:返回副本) |
|
| 布尔索引(相同!) |
⚠️ MATLAB 的
end 在 NumPy 中用-1 表示:A(:, end) →A[:, -1]A(2:end-1, :) →A[1:-1, :]
操作 | MATLAB | NumPy |
|---|---|---|
矩阵乘法 |
|
|
逐元素乘法 |
|
|
转置 |
|
|
求逆 |
|
|
解方程 |
|
|
特征值 |
|
|
FFT |
|
|
最大值 |
|
|
拼接 |
|
|
.py 文件)def func(...):)clear, clc → 不需要(Python 自动管理内存)% MATLAB
function y = myfunc(x)
y = x.^2 + 2*x + 1;
end# Python
def myfunc(x):
return x**2 + 2*x + 1问题 | MATLAB 行为 | NumPy 行为 | 解决方案 |
|---|---|---|---|
| 矩阵乘法 | 逐元素乘法 | 用 |
向量维度 | 明确区分行/列 | 一维数组无行列 | 需要列向量时用 |
整数除法 |
|
| 无问题 |
索引越界 | 报错 | 报错 | 相同 |
默认浮点类型 |
|
| 相同 |
Matplotlib 是 Python 最主流的绘图库,与 NumPy 天然集成。
import numpy as np
import matplotlib.pyplot as plt
# 1. 用 NumPy 生成数据
x = np.linspace(0, 2*np.pi, 100)
y = np.sin(x)
# 2. 用 Matplotlib 绘图
plt.plot(x, y, label='sin(x)')
plt.xlabel('x')
plt.ylabel('y')
plt.title('Sine Wave')
plt.legend()
plt.grid(True)
plt.show()✅ 关键点:Matplotlib 直接接受 NumPy 数组作为输入,无需转换。
x = np.random.randn(100)
y = np.random.randn(100)
colors = np.random.rand(100)
sizes = 1000 * np.random.rand(100)
plt.scatter(x, y, c=colors, s=sizes, alpha=0.5)
plt.colorbar()
plt.show()# 生成 2D 数组(如热力图、图像)
data = np.random.rand(10, 10)
plt.imshow(data, cmap='hot', interpolation='nearest')
plt.colorbar()
plt.show()data = np.random.normal(0, 1, 1000)
plt.hist(data, bins=30, alpha=0.7)
plt.show()from mpl_toolkits.mplot3d import Axes3D
x = np.linspace(-5, 5, 50)
y = np.linspace(-5, 5, 50)
X, Y = np.meshgrid(x, y) # ← 广播生成网格!
Z = np.sin(np.sqrt(X**2 + Y**2))
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
ax.plot_surface(X, Y, Z, cmap='viridis')
plt.show()🔥
np.meshgrid()是连接 NumPy 与 2D/3D 绘图的关键函数!
# 结构化数据
dt = np.dtype([('city', 'U10'), ('temp', 'f4'), ('humidity', 'f4')])
weather = np.array([
('Beijing', 30.0, 40.0),
('Shanghai', 35.0, 70.0),
('Guangzhou', 33.0, 80.0)
], dtype=dt)
# 绘制温度柱状图
plt.bar(weather['city'], weather['temp'])
plt.ylabel('Temperature (°C)')
plt.show() # ❌ 慢
for i in range(len(x)):
plt.plot(x[i], y[i], 'o')
# ✅ 快
plt.plot(x, y, 'o')plt.plot(..., rasterized=True) 或降采样工具 | 与 NumPy 的关系 |
|---|---|
Pandas | DataFrame 底层是 NumPy 数组; |
SciPy | 基于 NumPy,提供高级科学计算(优化、积分、信号处理等) |
scikit-learn | 输入要求为 |
TensorFlow / PyTorch | 张量可与 NumPy 数组互转( |
OpenCV | 图像即 NumPy 数组(H×W×C),可直接操作像素 |
* 改为 @(矩阵乘法)A(1,:) → A[0,:])np.zeros((m,n)) 而非 zeros(m,n)np.linalg.solve(A,b) 代替 A\bnp.arange() 或 np.linspace() 代替 :np.newaxis 或 reshape 处理维度np.save() / np.load() 代替 save / load📌 建议:新项目直接使用 NumPy + Matplotlib + Jupyter Notebook,体验远超 MATLAB!
本文主要讲述了numpy数据持久化。python过渡项目部分代码已经上传至gitee,后续会逐步更新,主要受时间原因限制,当然自己也可以克隆到本地学习拓展。
公众号:咚咚王
《Python编程:从入门到实践》
《利用Python进行数据分析》
《算法导论中文第三版》
《概率论与数理统计(第四版) (盛骤) 》
《程序员的数学》
《线性代数应该这样学第3版》
《微积分和数学分析引论》
《(西瓜书)周志华-机器学习》
《TensorFlow机器学习实战指南》
《Sklearn与TensorFlow机器学习实用指南》
《模式识别(第四版)》
《深度学习 deep learning》伊恩·古德费洛著 花书
《Python深度学习第二版(中文版)【纯文本】 (登封大数据 (Francois Choliet)) (Z-Library)》
《深入浅出神经网络与深度学习+(迈克尔·尼尔森(Michael+Nielsen)》
《自然语言处理综论 第2版》
《Natural-Language-Processing-with-PyTorch》
《计算机视觉-算法与应用(中文版)》
《Learning OpenCV 4》
《AIGC:智能创作时代》杜雨+&+张孜铭
《AIGC原理与实践:零基础学大语言模型、扩散模型和多模态模型》
《从零构建大语言模型(中文版)》
《实战AI大模型》
《AI 3.0》
原创声明:本文系作者授权腾讯云开发者社区发表,未经许可,不得转载。
如有侵权,请联系 cloudcommunity@tencent.com 删除。
原创声明:本文系作者授权腾讯云开发者社区发表,未经许可,不得转载。
如有侵权,请联系 cloudcommunity@tencent.com 删除。