擅长:python、mysql、java
<p>只需从一开始将指数<code>range</code>换成另一种方式,之后就不必翻转:</p>
<pre><code>def VandermondeMatrix(x):
n = len(x)
A = np.zeros((n, n))
exponent = np.array(range(n-1,-1,-1))
for j in range(n):
A[j, :] = x[j]**exponent
return A
</code></pre>
<p>输出:</p>
<pre><code>#A_matrix:
[[64. 16. 4. 1.]
[ 8. 4. 2. 1.]
[ 0. 0. 0. 1.]
[-1. 1. -1. 1.]]
#c:
[[ 6.]
[-39.]
[ 55.]
[ 27.]]
</code></pre>