定义
注:1.只要n+1个节点互异,满足上述插值条件的多项式是唯一存在的。
2.如果不限制多项式的次数,插值多项式并不唯一。
拉格朗日插值法
值分析中,拉格朗日插值法是以法国十八世纪数学家约瑟夫∙路易斯∙拉
格朗日命名的一种多项式插值方法。在若干个不同的地方得到相应的观测值,拉格朗日插值法可以找到一个多项式,其恰好在各个观测的点取到观测到的值。
龙格现象
高次插值会产生龙格现象,即在两端处波动极大,产生明显的振荡。在不熟悉运动趋势的前提下,不要轻易使用高次插值。
分段线性插值
牛顿插值
分段三次埃尔米特插值
直接使用Hermite插值得到的多项式次数较高,也存在着龙格现象,
因此在实际应用中,往往使用分段三次 Hermite 插值多项式 (PCHIP)。
Matlab有内置的函数(实现过程已经帮我们封装好了,会调用就行了):
p = pchip(x,y, new_x)
x是已知的样本点的横坐标;y是已知的样本点的纵坐标;new_x是要插入处对应的横坐标
例:
x = ‐pi:pi; y = sin(x);
new_x = ‐pi:0.1:pi;
p = pchip(x,y,new_x);
plot(x, y, 'o', new_x, p, 'r‐'
plot函数用法:
plot(x1,y1,x2,y2)
线方式: ‐ 实线 :点线 ‐. 虚点线 ‐ ‐ 波折线
点方式: . 圆点 +加号 * 星号 x x形 o 小圆
颜色: y黄; r红; g绿; b蓝; w白; k黑;
三次样条插值
Matlab有内置的函数:
p = spline (x,y, new_x)
x是已知的样本点的横坐标;y是已知的样本点的纵坐标;new_x是要插入处对应的横坐标
例:
x = -pi:pi;
y = sin(x);
new_x = -pi:0.1:pi;
p1 = pchip(x,y,new_x); %分段三次埃尔米特插值
p2 = spline(x,y,new_x); %三次样条插值
plot(x,y,'o',new_x,p1,'r-',new_x,p2,'b-')
legend('样本点','三次埃尔米特插值','三次样条插值','Location','SouthEast') %标注显示在东南
说明:
legend(string1,string2,string3, …)
分别将字符串1、字符串2、字符串3……标注到图中,每个字符串对应的图标为画图时的图标。‘Location’用来指定标注显示的位置。