黑龙江大学数值分析实验一拉格朗日插值算法

本实验代码仅供参考

一、实验目的
为了掌握数值分析的基本思想、原理和方法，要注意数值分析的处理技巧与计算机实现的结合，需将各种数值方法设计成算法，并编制好程序，拿到计算机上实现，最后得到可行性的验证。
1.学习和掌握拉格朗日插值多项式。
2.运用拉格朗日多项式进行运算。
3.验证随着插值节点的增多插值曲线的变化情况。
二、实验内容
利用Lagrange插值多项式求被插值函数f(x)在点x=65处的近似值。建议：画出Lagrange插值多项式的曲线。

如下是实验的代码片

def f1(x):# 返回拉格朗日插值
    s = 0
    lx = [0,10,20,30,40,50,60,70,80,90,100,110,120]
    ly = [5,1,7.5,3,4.5,8.8,15.5,6.5,-5,-10,-2,4.5,7]
    for i,j in zip(lx,ly):#双行遍历
        l = lx[:]
        l.remove(i)#移除Xk
        sum1,sum2 = 1,1
        for each in l:
            sum1 *= x - each
            sum2 *= i - each
        s += j * sum1 / sum2#插值结果
    return s
import matplotlib.pyplot as plt#画图函数
print(f1(65))
lx = [each for each in range(121)]#从0-120的列表
ly = [f1(i) for i in lx]#与x相对应的f(x)值
plt.plot(lx,ly)#图像绘制
plt.title('Lagrange')#标题
plt.xlabel('x')#x轴
plt.ylabel('y')#y轴
plt.show()

代码结果12.859933590888977
图像代码运行可得