Python实现拉格朗日插值法

Posted

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了Python实现拉格朗日插值法相关的知识,希望对你有一定的参考价值。

已知sinx的一组x,y对应关系,用拉格朗日插值法估计sin(0.3367)的值.

x x0.32 0.34 0.36
y 0.314567 0.333487 0.352274

/
/

class Interpolation:
    def __init__(self, x, y):
        self.x = x
        self.y = y

    def func(self, X):
        s = 0
        for i in range(len(self.x)):
            W = 1
            w = (X - self.x[i])
            for j in range(len(self.x)):
                W *= (X - (self.x[j]))
                if i == j:
                    continue
                else:
                    w *= (self.x[i] - self.x[j])
            L = self.y[i] * W / w
            s += L
        print(round(s, 8))

while (True):
    x=input(‘x:‘)
    y=input(‘y:‘)
    X=input(‘X:‘)
    if x==‘Q‘ or y==‘Q‘ or X==‘Q‘:
        break
    else:
        p1 = Interpolation(eval(x),eval(y))   # 前面是x的取值,后面是y的取值;
        p1.func(eval(X))   # 现在要计算的是  sin(0.3367)  的值

程序还有好多不足的问题需要改进(^o^),

以上是关于Python实现拉格朗日插值法的主要内容,如果未能解决你的问题,请参考以下文章

Python实现拉格朗日插值法

数值分析Python实现系列—— 拉格朗日插值法

python拉格朗日插值

MATLAB实现拉格朗日插值法

拉格朗日插值与牛顿插值

拉格朗日插值法