一、原因
在实现拉格朗日插值时遇到一个问题,即:
当我们手算时用拉格朗日很容易得到具体表达式,但是当我们用机器实现时,却很难得到,因为机器不能表达出未知数x(虽然可以用TensorFlow的占位符,但杀鸡焉用牛刀),因此我们能很容易的得到具体的近似解,却无法得到表达式。
但是,python的scipy库中的lagrange实现了这一功能,如图:
所以我不禁思考如何实现这一功能,后来我换了一个思路,将系数当成方程组解,用高斯消元法求解。
二、实现代码(python)
import numpy as np
import matplotlib.pyplot as plt
from scipy.interpolate import lagrange
def main_element_gauss(A,n):
for i in range(0,n-