我有(N0, N1, N2, N3)矩阵V和(N1, N1)矩阵M。 N1通常约为30-50，N0xN1xN2xN3约为1,000,000。我想要一个新的矩阵EV，其中i0, i1, i2, i3条目由：
sum(V(i0, :, i2, i3) * M(i1, :)', 2)
我需要做很多次这个操作，主要关注的是速度。
我前一段时间使用numpy问过同一个问题的解决方案（see here）有人建议我使用np.einsum，这非常有效。我现在需要将整个事情移植到Matlab，我想知道在这种情况下最有效的解决方案是什么。
提前谢谢。