如何用numpy concatenate实现类似广播的外连接操作？

针对你需要实现的「外连接操作」，可以直接利用NumPy的广播机制结合np.concatenate或np.stack来实现，比meshgrid的方式更简洁高效：

方法1：广播 + np.concatenate #

通过扩展数组维度触发广播，再沿最后一维拼接：

import numpy as np

a = np.array([0.0, 10.0, 20.0, 30.0])
b = np.array([1.0, 2.0, 3.0])

# 扩展a的维度为(4,1,1)，b的维度为(1,3,1)，广播后均变为(4,3,1)
result = np.concatenate([a[:, None, None], b[None, :, None]], axis=-1)

方法2：广播 + np.stack（更简洁） #

np.stack可以直接在指定维度合并数组，省去手动扩展最后一维的步骤：

import numpy as np

a = np.array([0.0, 10.0, 20.0, 30.0])
b = np.array([1.0, 2.0, 3.0])

# 先扩展维度让a和b满足广播条件：a变为(4,1)，b变为(1,3)
a_expanded = a[:, None]
b_expanded = b[None, :]

# 沿最后一维堆叠，自动广播为(4,3)后合并为(4,3,2)
result = np.stack([a_expanded, b_expanded], axis=-1)

这两种方法的优势：

• 直接利用NumPy原生的广播机制，逻辑更直观，避免了meshgrid后的转置操作
• 生成的中间数组更少，计算效率更高
• 代码行数更少，符合Pythonic的简洁风格

内容的提问来源于stack exchange，提问作者Guillaume Mougeot