三次（自然）样条插值及其Python代码实现

最新推荐文章于 2024-12-08 22:47:30 发布

原创

最新推荐文章于 2024-12-08 22:47:30 发布 · 2.1k 阅读

12 ·

CC 4.0 BY-SA版权

文章标签：

#python

本文介绍了如何使用三次样条插值方法对给定点进行插值，并通过Python的numpy和matplotlib库实现。通过计算系数矩阵和常数矩阵，求解插值函数的参数，最后绘制了插值曲线。

前言：

利用一个未知理论推出的结果来证明该未知理论是一种荒谬的行为……

理论：

参考《数值分析》这本书（为什么不参考课上PPT，是不想吗QAQ）

代码：

（强烈建议先看一看理论部分的讲解，尝试自己写一写，其实并不难的）

# 三次样条插值
import numpy as np
from matplotlib import pyplot as plt

''' 
定义插值函数为: S_i(x_i) = y_i + b_i(x - x_i) + c_i(x - x_i)^2 + d_i(x - x_i)^3
'''

#  输入需要插值的点
ls = [(-1.00, -14.58), (0.00, 0.00), (1.27, 0.00), (2.55, 0.00), (3.82, 0.00), (4.92, 0.88), (5.02, 11.17)