Python中Quiver绘制磁场矢量图异常问题排查

核心问题排查方向 #

1. 维度切片错误 #

你要绘制y=0平面，但如果代码中对数组的切片维度搞反（混淆y和z维度），就会导致显示的不是目标平面的数据：

• 正确切片：若磁场数组形状为(x, y, z)，y=0平面的x方向分量应取bx[:, 0, :]；如果误写为bx[:, :, 0]，实际取的是z=0平面的数据，这会直接导致你看到的现象。
• 验证方法：执行print(bx[:, 0, 7])查看y=0、z=7位置的x分量是否全0，再执行print(bx[:, 0, 0])确认z=0位置的x分量状态——如果结果和你描述的现象相反，说明切片维度完全搞反了。

2. Quiver参数与坐标网格不匹配 #

Matplotlib的quiver(X, Y, U, V)有严格的参数对应规则：

• X、Y：箭头位置的坐标网格，形状必须和U、V（箭头分量）完全一致。
• U：箭头在X轴方向的分量，V：箭头在Y轴方向的分量。

常见错误：

• 网格索引模式错误：用np.meshgrid(x, z)默认是indexing='xy'，生成的网格形状为(nz, nx)，但你的bx[:,0,:]形状是(nx, nz)，维度不匹配会导致数据错位显示。解决方法是生成网格时指定indexing='ij'，保证网格形状和分量数组一致：

  X, Z = np.meshgrid(x_axis, z_axis, indexing='ij')
  

• 分量-轴对应错误：若目标平面是x-z平面（y=0），U应对应x方向分量bx[:,0,:]，V应对应z方向分量bz[:,0,:]；如果误把V设为y方向分量by[:,0,:]，会显示错误的矢量方向。

3. 数据与网格维度不匹配 #

确保坐标网格（如x、z轴）的长度，和磁场数组对应维度的长度完全一致。例如x轴长度为nx、z轴长度为nz，那么bx[:,0,:]的形状必须是(nx, nz)，网格X、Z的形状也必须是(nx, nz)。

修正示例代码 #

import numpy as np
import matplotlib.pyplot as plt

# 定义轴参数
nx, ny, nz = 10, 15, 10
x = np.linspace(0, 1, nx)
y = np.linspace(-1, 1, ny)
z = np.linspace(0, 1, nz)

# 模拟磁场数据：y=0时z=7处bx全0
bx = np.random.rand(nx, ny, nz)
bz = np.random.rand(nx, ny, nz)
bx[:, 0, 7] = 0  # 目标位置设为全0

# 生成y=0平面的坐标网格（x-z平面）
X, Z = np.meshgrid(x, z, indexing='ij')

# 提取y=0平面的分量
u = bx[:, 0, :]  # x方向分量
v = bz[:, 0, :]  # z方向分量

# 绘制quiver图
plt.quiver(X, Z, u, v)
plt.xlabel('x')
plt.ylabel('z')
plt.title('Magnetic Field at y=0 Plane')
plt.show()

内容的提问来源于stack exchange，提问作者Guillermo Fuentes Morales