数值计算笔记之插值(四)三次样条插值

0、定义已知函数在区间上个互异节点,处的函数值为,若构造函数,满足: 在每个小区间上是一个不超过三次的多项式 在上连续则称为的三次样条插值函数。根据定义知道规律为:已知:n+1个数据点[xi, yi],i = 0, 1, …, n 每一分段都是三次多项式函数曲线 节点达到二阶连续 左右两端点处特性(自然边界,固定边界,非节点边界)根据定点,求出每段样条曲线方...
摘要由CSDN通过智能技术生成

0、定义

已知函数f(x)在区间[a,b]n+1个互异节点,a=x_{0}<x_{1}<\cdots < x_{n}=b处的函数值为y_{i}=f(x_{i}),若构造函数s(x),满足:

  1. s(x_{i})=y_{i}(i=0,1,\cdots ,n)
  2. 在每个小区间[x_{i},x_{i+1}]上是一个不超过三次的多项式
  3. s(x),s^{'}(x),s^{''}(x)[a,b]上连续

则称s(x)f(x)的三次样条插值函数。

根据定义知道规律为:

已知:

  • n+1个数据点[xi, yi], i = 0, 1, …, n
  •  每一分段都是三次多项式函数曲线
  •  节点达到二阶连续
  • 左右两端点处特性(自然边界,固定边界,非节点边界)

根据定点,求出每段样条曲线方程中的系数,即可得到每段曲线的具体表达式。

一、推导

插值和连续性:

S_{i}(x_{i}) = y_{i}

S_{i}(x_{i+1}) = y_{i+1}

  其中i = 0, 1,\cdots , n-1.

微分连续性:

S^{'}_{i}(x_{i+1})=S^{'}_{i+1}(x_{i+1})

S^{''}_{i}(x_{i+1}) = S^{''}_{i+1}(x_{i+1})

其中i = 0, 1,\cdots , n-2.

样条曲线的微分式:

S_{i}(x) = a_{i}+b_{i}(x-x_{i})+c_{i}(x-x_{i})^{2}+d_{i}(x-x_{i})^{3}

S^{'}_{i}(x) = b_{i}+2c_{i}(x-x_{i})+3d_{i}(x-x_{i})^{2}

S^{''}_{i} = 2c_{i}+6d_{i}(x-x_{i})

 

将步长 h_{i} = x_{i+1}-x_{i} 带入样条曲线的条件:

  • 14
    点赞
  • 126
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值