拉格朗日插值法拟合多项式

最新推荐文章于 2021-12-29 22:14:28 发布

水能zai舟

最新推荐文章于 2021-12-29 22:14:28 发布

阅读量2.4k

点赞数

CC 4.0 BY-SA版权

本文链接：https://blog.youkuaiyun.com/easylovecsdn/article/details/84865660

本文介绍了使用拉格朗日插值法来拟合最高次为2次的多项式，提供了随着n变化的通用代码实现。同时，作者表达了对提升图表绘制质量的求助。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

咱们先来拟合一个最高次为2次的多项式（此时固定n = 3）

Code

#拉格朗日插值公式

x = []
y = []

n = 3

for i in range(n) :
    v = input()
    v = int(v)
    x.append(v)

for i in range(n) :
    v = input()
    v = int(v)
    y.append(v)

def f(X) :
    return ((X - x[1]) * (X - x[2]) / (x[0] - x[1]) * (x[0] - x[2])) * y[0] + ((X - x[0]) * (X - x[2]) / (x[1] - x[0]) * (x[1] - x[2])) * y[1] + ((X - x[0]) * (X - x[1]) / (x[2] - x[0]) * (x[2] - x[1])) * y[2]

xx = input()
xx = float(xx)

print(f(xx))

当n不定时，将代码普式化有（n自定义输入）

Code

#拉格朗日插值公式

x = []
y = []

n = input()
n = int(n)

for i in range(n) :
    v = input()
    v = int(v)
    x.append(v)

for i in range(n) :
    v = input()
    v = int(v)
    y.append(v)

def f(X) :
    sum = 0