插值算法(数学建模)

定义

在这里插入图片描述
在这里插入图片描述
注: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’用来指定标注显示的位置。

  • 5
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 1
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

每天一道题

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值