<p>当我计算Vandermonde时,我得到了正确的答案
这个矩阵的系数。然而,输出矩阵是相反的。
它应该是<code>[6,-39,55,27]</code>,而不是<code>[27,55,-39,6]</code></p>
<p>我的Vandermonde矩阵的输出被翻转,最终的解决方案
c、 他被打翻了</p>
<pre><code>import numpy as np
from numpy import linalg as LA
x = np.array([[4],[2],[0],[-1]])
f = np.array([[7],[29],[27],[-73]])
def main():
A_matrix = VandermondeMatrix(x)
print(A_matrix)
c = LA.solve(A_matrix,f) #coefficients of Vandermonde Polynomial
print(c)
def VandermondeMatrix(x):
n = len(x)
A = np.zeros((n, n))
exponent = np.array(range(0,n))
for j in range(n):
A[j, :] = x[j]**exponent
return A
if __name__ == "__main__":
main()
</code></pre>
<p>你可以
<code>print(c[::-1])</code>
这将颠倒c的顺序。
从<a href="https://stackoverflow.com/questions/29806239/how-can-i-flip-the-order-of-a-1d-numpy-array">How can I flip the order of a 1d numpy array?</a></p>