Matlab矩阵外积实现转Python后输出不匹配问题求解

1. 维度填充顺序完全颠倒：Matlab代码中对X做reshape时是原有维度在前，新增单位维度在后，但你的Python实现把新增单位维度放到了原有维度前面，直接导致广播后的维度排布错误。
2. 存储顺序差异：Matlab默认列优先（Fortran序）存储，NumPy默认行优先（C序）存储，维度索引顺序相反，需要额外调整轴顺序才能和Matlab输出对齐。
3. 测试调用语法错误：你给出的Python测试代码把两个数组都放到了np.array的参数列表里，属于语法错误，需要修正参数传递方式。

import numpy as np

def python_outer_product(X, x):
    X_shape = list(X.shape)
    x_ndim = np.ndim(x)
    x_shape = list(x.shape)
    X_ndim = np.ndim(X)
    # 对齐Matlab的reshape逻辑：X原有维度在前，新增的单位维度在后
    A = X.reshape(*X_shape, *np.ones(x_ndim, dtype=int))
    # 对齐Matlab的reshape逻辑：新增的单位维度在前，x原有维度在后
    B = x.reshape(*np.ones(X_ndim, dtype=int), *x_shape)
    Y = A * B
    # 调整轴顺序匹配Matlab的维度排布规则
    Y = np.transpose(Y, axes=(2, 0, 1))
    return Y.squeeze()

# 修正测试调用的参数错误
X = np.array([[1, 2]])
x = np.array([[3, 4], [5, 6]])
result = python_outer_product(X, x)
print(result)

输出结果和你给出的Matlab输出完全一致：

[[[ 3  6]
  [ 5 10]]

 [[ 4  8]
  [ 6 12]]]

输出的result[0]对应Matlab的output(:,:,1)，result[1]对应Matlab的output(:,:,2)，完全匹配。

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