查看: 2769|回复: 0

拉格朗日插值 python实现

字体大小: 正常放大

1186 主题	4 听众	2922 积分

该用户从未签到

电梯直达

1^#

发表于 2024-3-14 11:02 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

让我们逐步解析和详细解释这段代码。
1. 导入必要的库

import numpy as np ]5 S: t/ Q. z7 ^
' Z% r* |' n7 y$ U5 M/ K& [
from scipy.interpolate import lagrange
/ g% i/ r+ M. n( J) [, `4 Z% w
\" A/ B; D, f, U& e Z- f
import matplotlib.pyplot as plt

复制代码

1.numpy: 一个强大的数学库，用于处理数组和矩阵等高级数学运算。
2.scipy.interpolate.lagrange: 从 SciPy 库中导入的，用于执行拉格朗日多项式插值的函数。
3.matplotlib.pyplot: 用于数据可视化的库，可以绘制图表和图形。

2. 源数据

x = np.array([4, 5, 6])6 z' x( ?! H\" W
7 t8 |! R8 [6 b) @ c( H: d
y = np.array([10, 5.25, 1])

复制代码

这里定义了两个数组 x 和 y，代表了三个数据点的坐标：(4, 10), (5, 5.25), 和 (6, 1)。这些数据点将用于拉格朗日插值。
3. 拉格朗日插值

poly = lagrange(x, y)2 M# e2 p% O, O$ J4 i) L0 c
9 g& o9 Y$ A; y4 ]9 g7 R# d/ @% E
res = poly(18)

复制代码

4.使用 lagrange(x, y) 创建一个拉格朗日插值多项式 poly。这个多项式是基于给定的数据点 x 和 y 计算出来的，它能够逼近这些点构成的曲线。
5.poly(18) 计算了在 x = 18 时的插值结果 res。这是通过将 x = 18 代入到拉格朗日多项式中得到的。

4. 打印多项式和插值结果

print('poly =\n', poly)
) J/ Z$ ^( I/ V, E8 _' k
' `$ u5 p) G) l0 b, g
print('res =', res)

复制代码

这里打印了插值多项式 poly 和在 x = 18 处的插值结果 res。
5. 数据可视化

x1 = [i for i in range(-5, 20)]
4 H8 x) Z/ x4 T\" Y: O
: T9 q* H# l H& ~! _& F1 R
y1 = []
. n' S% L# y3 h9 ], u2 w5 i
) a1 A: k9 D- t- m n* C
for each in x1:8 ~5 e I/ W r: K0 k' h- L
; I! J' G/ Z& A- w: t
y1.append(poly(each))
T+ R8 D1 }- m# {
( r4 }/ v# m% _3 F9 D- I
T0 q7 B4 X1 `8 m4 h
# L+ Q! C. S' X- p: t
plt.plot(x1, y1)
/ F/ N' |3 I5 ~) j' u8 w2 d3 O
1 z5 `. `5 [1 V\" h9 R% T
plt.scatter(x, y)8 {\" E8 J3 m. _ i# [\" G, B
9 X% b4 o d0 X3 m: h( T1 R
plt.scatter(18, res)

复制代码

6.首先，创建了一个 x 值列表 x1，范围是从 -5 到 19。然后，对于 x1 中的每个 x 值，都计算对应的拉格朗日插值多项式的 y 值，并将这些 y 值存储在列表 y1 中。
7.使用 plt.plot(x1, y1) 绘制插值多项式形成的曲线。
8.使用 plt.scatter(x, y) 在图上标记原始数据点。
9.使用 plt.scatter(18, res) 在图上标记 x = 18 处的插值结果点。

这样，你就得到了一幅图形，显示了原始数据点、插值曲线，以及在 x = 18 处的插值点。这个过程展示了如何使用拉格朗日插值方法根据一组给定的数据点来估计函数在特定点上的值，并将结果可视化。