擅长:python、mysql、java
<p>让我们试试这个-</p>
<pre class="lang-py prettyprint-override"><code>
# self.t x self.N > self.N x self.t x 1
# arranging one column in one layer
self.Bsolution = np.swapaxes(np.expand_dims(self.Bsolution, axis=-1), 0, 1)
# self.N x self.t x self.N == N columns of Bsolution multiplied with N columns of Asolution
# multiplying 1 column of B with N columns of A in a layer done for N layers
product = np.multiply(self.Bsolution, self.Asolution)
# self.N x self.N > self.N x 1 x self.N
self.Matrix = np.expand_dims(self.Matrix, axis=1)
# integration need to take place along axis=1 i believe
# multiplying each of N x N columns(4 rows in length) with N x N matrix elements
result = np.sum(simps(np.multiply(self.Matrix, product), self.time))
</code></pre>
<p>让我知道它是否有效</p>