擅长:python、mysql、java
<p>使用<a href="https://en.wikipedia.org/wiki/Kronecker_product" rel="nofollow noreferrer">Kronecker product</a>(<a href="https://docs.scipy.org/doc/numpy-1.13.0/reference/generated/numpy.kron.html" rel="nofollow noreferrer">^{<cd1>}</a>):</p>
<pre><code>>>> import numpy as np
>>> A = np.array([[ 1, 2, 3],
... [100, 200, 300]])
>>> B = np.array([[ 10, 20, 30],
... [1000, 2000, 3000],
... [ -10, -20, -2]])
>>> (m,c) = A.shape
>>> (n,c) = B.shape
>>> np.kron(A,np.ones((n,1))) - np.kron(np.ones((m,1)),B)
array([[ -9., -18., -27.],
[ -999., -1998., -2997.],
[ 11., 22., 5.],
[ 90., 180., 270.],
[ -900., -1800., -2700.],
[ 110., 220., 302.]])
</code></pre>