Hermite插值可以看作牛顿插值的极限状况.为什么可以这么说呢?我们来看一个实例:

构造一个三次多项式 $p_3$ 使得 $p_3(0)=0$,$p_3(1)=1,p_3'(0)=1,p_3'(1)=0$.

证明:我们进行牛顿插值.不妨构造这么几个插值点:

\begin{equation}
x_0,x_1,x_2,x_3
\end{equation}
其中$x_0=0,x_2=1$.$x_1$和$x_0$足够接近,且$x_1>x_0$.使得
\begin{equation}
f[x_0,x_1]
\end{equation}
约等于1.$x_3>x_2$且$x_3$和$x_2$足够接近,使得
\begin{equation}
f[x_2,x_3]
\end{equation}
约等于0.

这样构造出来的三次牛顿插值多项式和本题所需要的Hermite插值多项式的要求应该不会相差太远.根据牛顿插值公式,可得这样的多项式为

\begin{align*}
f(x)=f(0)+(x-0)f[0,x_1]+(x-0)(x-x_1)f[0,x_1,1]+(x-0)(x-x_1)(x-1)f[0,x_1,1,x_3]
\end{align*}

由于$x_1$和$x_0$足够接近,$x_3$和$x_2$足够接近,因此我们不妨把这个牛顿插值得来的多项式取极限.得到Hermite插值所需要的多项式为

\begin{align*}
f(0)+(x-0)f'(0)+(x-0)^2f[0,0,1]+(x-0)^2(x-1)f[0,0,1,1]
\end{align*}

现在,我们来计算
\begin{align*}
f[0,0,1]=\frac{f[0,0]-f[0,1]}{0-1}=\frac{f'(0)-f[0,1]}{0-1}=0
\end{align*}

\begin{align*}
f[0,0,1,1]=\frac{f[0,0,1]-f[0,1,1]}{0-1}
\end{align*}

为此,我们先计算$f[0,1,1]$.
\begin{align*}
f[0,1,1]=\frac{f[0,1]-f[1,1]}{0-1}=-1
\end{align*}因此,\begin{align*}
f[0,0,1,1]=-1
\end{align*}.

因此可以求得满足条件的hermite多项式为

\begin{align*}
x-x^2(x-1)=-x^3+x^2+x
\end{align*}

一般地,已知这些数据信息:

\begin{align*}
(x_0,f(x_0)),(x_2,f(x_2)),\cdots,(x_{2n},f(x_{2n})),(x_0,f'(x_0)),(x_2,f'(x_2)),\cdots,(x_{2n},f'(x_{2n}))
\end{align*}

求构造一个多项式$p_{2n+1}(x)$满足这些条件.

解答:为此,我们先构造牛顿插值多项式.

\begin{equation}
x_0<x_1<x_2<\cdots <x_{2n}<x_{2n+1}
\end{equation}

其中$x_1$与$x_0$足够接近,$x_3$与$x_2$足够接近,...,$x_{2n+1}$与$x_{2n}$足够接近,使得

\begin{equation}
f[x_{2k+1},x_{2k}]
\end{equation}

约等于$f'(x_{2k})$.$(0\leq k\leq n)$.

根据牛顿插值公式,我们知道牛顿插值多项式为

\begin{align*}
f(x_0)+(x-x_0)f[x_0,x_1]+(x-x_0)(x-x_1)f[x_0,x_1,x_2]+\cdots
+(x-x_0)(x-x_1)(x-x_2)\cdots (x-x_{2n})f[x_0,x_1,x_2,\cdots,x_{2n+1}]
\end{align*}

取极限,可得相应的Hermite插值多项式为

\begin{align*}
f(x_0)+(x-x_0)f'(x_0)+(x-x_0)^2f[x_0,x_0,x_2]+\cdots+(x-x_0)^2(x-x_2)^2\cdots (x-x_{2n-2})^2(x-x_{2n})f[x_0,x_0,x_2,x_2,\cdots,x_{2n-2},x_{2n-2},x_{2n},x_{2n}]
\end{align*}

转载于:https://www.cnblogs.com/yeluqing/archive/2012/12/17/3827566.html

Hermite插值是牛顿插值的极限情形相关推荐

  1. 拉格朗日插值和牛顿插值的龙格现象

    文章目录 一.实验目的 二.实验设备信息 三.实验内容 (一)拉格朗日插值多项式 (二)牛顿插值多项式 四.实验步骤 (一)拉格朗日插值函数实现 (二)牛顿插值函数实现 (三)观察拉格朗日插值和牛顿插 ...

  2. matlab全域基函数,多项式函数插值:全域多项式插值(一)单项式基插值、拉格朗日插值、牛顿插值 [MATLAB]...

    全域多项式插值指的是在整个插值区域内形成一个多项式函数作为插值函数.关于多项式插值的基本知识,见"计算基本理论". 在单项式基插值和牛顿插值形成的表达式中,求该表达式在某一点处的值 ...

  3. Matlab实现线性插值、抛物插值、牛顿插值、拉格朗日插值、分段抛物插值、分段线性插值

    目录 线性插值 原理 流程图 代码 抛物插值 原理 流程图 代码 拉格朗日插值 代码 牛顿插值 原理 代码 分段线性插值 代码 线性插值 原理 流程图 单个点的线性插值代码 X=[0.2 0.4]; ...

  4. 【插值】牛顿插值、拉格朗日插值、三次样条插值的Python代码实现

    插值简介 插值即根据有限的离散点绘制出穿过所有样本点的曲线,从直观上想象似乎画一条穿过n个特定点的曲线有无数种画法,但从数学意义上来说我们希望画出的曲线能够尽量平滑,震荡幅度尽量小能够在非样本点上符合 ...

  5. 插值问题(拉格朗日插值、牛顿插值)

    agui_lagrange.m: function f=agui_lagrange(x0,y0,x) % x0为节点向量,y0为节点上的函数值,x为插值点,f为返回插值 n=length(x0);m= ...

  6. 数值分析之 拉格朗日插值、牛顿插值、分段线性插值实现

    1.拉格朗日插值法 考虑全局信息的比较经典的插值方法,编程简单,计算量大. #coding=utf-8 from matplotlib import pyplot as pltdef Lg(data, ...

  7. 1月16日:拉格朗日中值定理,罗尔定理,柯西中值,拉格朗日插值,牛顿插值,重心插值,拉格朗日乘子法的证明

    拉格朗日中值定理 https://www.bilibili.com/video/BV117411E7kx?from=search&seid=17921778669593975548 拉格朗日中 ...

  8. 插值法:拉格朗日插值、牛顿插值

    拉格朗日插值法 (*以下定义选自维基百科) 算法流程图 算法代码 [cpp] view plaincopy #include<iostream> #include<string> ...

  9. 【数值分析】插值法:拉格朗日插值、牛顿插值

    本科课程参见:<软件学院那些课> 拉格朗日插值法 (*以下定义选自维基百科) 算法流程图 算法代码 #include<iostream> #include<string& ...

最新文章

  1. 开源助力!武汉新型冠状病毒防疫开源信息收集平台
  2. 网站开启https后很慢_HTTPS会影响网站打开速度吗
  3. JavaScript实现网页元素的拖拽效果
  4. 原生Js 两种方法实现页面关键字高亮显示
  5. densenet网络结构_FC-DENSENET用于图像分割学习笔记
  6. UNDO 表空间管理
  7. 解决ssh或ftp下root用户认证失败问题
  8. 最短路问题_Dijkstra算法
  9. python怎样计算增长率_增长率如何计算
  10. TM中拒收自定义表情的设置方法(转)
  11. Cannot find name ‘console‘. Do you need to change your target library?ging the ‘lib‘ compiler option
  12. elasticjob已下线_elasticJob 源码解析之自诊断恢复
  13. 【Java攻城狮宝典】04-for循环
  14. 课时23 YUM软件管理
  15. 处女作真人语音计算器上线了
  16. CCS安装失败的原因
  17. J-LINK 烧录程序
  18. docker容器内pip install 显示 warning
  19. html、css基础---学习笔记
  20. 输出图案-平行四边形

热门文章

  1. 【廖雪峰python进阶笔记】模块
  2. 还需要“attention”吗?一堆“前馈层”在ImageNet上表现得出奇得好
  3. 总结 | 机器学习的通俗讲解!
  4. 无人机航拍记录生活真爽,包邮送一个!
  5. SOTA论文也未必能被接收,谷歌科学家Eric Jang谈顶会审稿标准
  6. 美团架构师开源5万字的《面试手册》PDF免费下载!
  7. k-means聚类算法从入门到精通
  8. 10 年前被删的初恋,凌晨 1 点突然加我…
  9. 视学算法第六轮送书活动,16本高质量书免费拿走!
  10. Nginx源码分析:epoll事件处理模块概述