无论是牛顿插值还是拉格朗日插值,都只能保证在节点处的函数值没有误差。hermite插值更加复杂,可以保证一阶导数也连续,目前常用的是三次样条插值

一、三次样条插值概念

  • 不超过3次
  • 节点处无误差
  • 一阶导数和二阶导数节点处无误差


如果函数值和函数在点的倒数值是已知的,也就是我们有一下的已知条件

这不就是hermite插值的已知条件吗,我们可以利用Hermite插值先进行一次插值:

但是此时其实我们是不知道 m i m_i mi​的值的,所以我们需要求出他,而现在我们只剩下了一个条件,端点处二阶导数连续:令 h i h_i hi​表示 i − i + 1 i-i+1 i−i+1区间的长度,将 s ( x ) s(x) s(x)转化成下面的形式

然后我们求一个二阶导数:

因为我们的二阶导数需要在端点处连续,也就是说在 [ x i , x i − 1 ] [x_i,x_{i-1}] [xi​,xi−1​]区间的 s ′ ′ ( x i ) s''(x_i) s′′(xi​)要等于区间 [ x i − 1 , x i ] [x_{i-1},x_i] [xi−1​,xi​]上二阶导数的对 x i x_i xi​的值,所以我们将上式的下标 i i i换成 i − 1 i-1 i−1,

此时我们根据区间端点处二阶导数连续的定义,得到下面的等式:
变化不复杂,就是简单的移相合并而已



画圈的部分是我们必须化为1的部分,只有这样我们才能吧 m i m_i mi​单独拿出来,进行一系列的化简:

这不就是一个关于 m i m_i mi​的方程组

计算方法:三次样条插值原理相关推荐

  1. 三次样条插值原理及openCV实现三种边界条件(CSDN为数不多的正确版本)

    没有自定义目录标题 前情总结 算法介绍及原理解析 论证边界问题 边界条件介绍 公式推导 方程组 算法步骤 代码实现 前情总结 同事在工作中遇到需要样条插值的情况,帮他找实现代码的时候想根据博客推一遍原 ...

  2. 【20220207】【信号处理】三次样条插值原理详解

    方程组的求解本文不做介绍. 一.三次样条插值 1. 定义 三次样条插值(Cublic Spline Interpolation),简称 Spline 插值,是通过一系列样本点的光滑曲线,数学上通过求解 ...

  3. python 插值_三次样条插值在Python中的实现

    什么是三次样条插值 三次样条插值(Cubic Spline Interpolation)简称Spline插值,是通过一系列形值点的一条光滑曲线,数学上通过求解三弯矩方程组得出曲线函数组的过程. 实际计 ...

  4. matlab spline三次样条插值x,Spline(三次样条插值)

    关于三次样条插值,计算方法比较复杂,但是静下心来仔细研究也是可以理解的. 本文借鉴文章来源:http://www.cnki.com.cn/Article/CJFDTotal-BGZD200611035 ...

  5. Spline(三次样条插值)

    关于三次样条插值,计算方法比较复杂,但是静下心来仔细研究也是可以理解的. 本文借鉴文章来源:http://www.cnki.com.cn/Article/CJFDTotal-BGZD200611035 ...

  6. Cubic spline(三次样条插值)(转载)

    转自:http://blog.csdn.net/lsxpu/article/details/38849775 自己以前上过数值分析这门课,用的是[1]这本教材,三次样条插值这一节,当时似乎看明白了,但 ...

  7. python三次样条插值拟合的树行线_数学建模笔记——插值拟合模型(一)

    啊好像距离上次写作又过了七天,啊好像我之前计划的一周两三篇,啊辣鸡小说毁我青春,啊我是一只可怜的鸽子. 不管怎样,我又回来了,并坚定地更新着hhh.再过两三天就是我们学校数学建模选拔,再过八九天就是期 ...

  8. Easy3D 三次样条插值

    文章目录 一.简介 二.实现代码 三.实现效果 四.参考资料 一.简介 在插值计算中,最简单的分段多项式近似应该是分段线性插值,它由连接一组数据点组成,仅仅只需要将这些点一一用直线进行顺序相连即可.不 ...

  9. 三次样条插值(Cubic Spline Interpolation)及代码实现(C语言)

    样条插值是一种工业设计中常用的.得到平滑曲线的一种插值方法,三次样条又是其中用的较为广泛的一种.本篇介绍力求用容易理解的方式,介绍一下三次样条插值的原理,并附C语言的实现代码. 1. 三次样条曲线原理 ...

最新文章

  1. 强者愈强!疫情拉大“数据资产”贫富差距,顶级公司数据建设靠什么
  2. Web模板引擎—Mustache
  3. 信守承诺:JAX-RS API的基于合同的测试
  4. cv2 imshow窗口大小_cv2.imshow()图片无法显示
  5. PageObjects 设计模式
  6. new malloc
  7. 配置鸿蒙Windows烧录环境 用Hiburn烧录第一个程序
  8. Unity3D实践3:BOSS血条
  9. I2C双向电平转换电路
  10. 如果你想成为一名黑客
  11. 怎样找回win7密钥
  12. linux卷影复制功能,vssadmin 卷影复制服务管理命令行工具
  13. linux 打印文件目录结构,linux基本文件目录的结构及管理
  14. 2016年俄罗斯M2M市场达100亿俄罗斯卢布
  15. vast显示不出来服务器,VAST实用小知识?
  16. Unity:简单的蓄力技能的逻辑实现
  17. 飞利浦Georg Jensen联名系列真无线耳机JT60重磅上市
  18. 一大波无门槛优惠券来袭(仅限300张)
  19. 计算机品牌及介绍,【推广】电脑电源参数和品牌的介绍
  20. Unity游戏画面参数解析与应用:垂直同步、动态模糊、抗锯齿

热门文章

  1. 韶音骨传导耳机值得入手吗?南卡和韶音骨传导耳机实测对比
  2. Ubuntu操作系统安装
  3. PS系列之 -- 通道抠图进阶 : 用通道抠取头发
  4. systemd的日志存在哪里?
  5. 攻略:苹果手机投屏电脑 iPhone镜像投屏怎么操作
  6. Spring 实战(第 5 版)
  7. BM13-判断一个链表是否为回文结构
  8. VisualVM的配置和使用
  9. Linux服务器连接校园网
  10. MySQL学习之一条SQL更新语句的执行