数学建模第三讲----插值算法+MAYTLAB绘图

目录

一.适用情况

二.常见但建模不会去用的插值方法

三.比赛会用到的插值方法----分段三次埃尔米特插值&三次样条插值

（1）分段埃尔米特插值

（2） 三次样条插值

(3)n维数据的插值（不常用）

（4）可以使用插值算法来预测数据

四.MATLAB画图

1.plot画图函数用法

2.figure(n)

---

一.适用情况

已知的数据较少，需要根据已知数据再模拟创造出一些新的数据

二.常见但建模不会去用的插值方法

• 拉格朗日插值法
• 牛顿插值法

注意：

1.当插值多项式次数过高时，会出现Runge(龙格)现象；为了解决这一现象，使用分段插值的方法

2.同时这两种方法都无法满足求导的需求，因此建模比赛中不会去用

三.比赛会用到的插值方法----分段三次埃尔米特插值&三次样条插值

（1）分段埃尔米特插值

代码示例：

% 分段三次埃尔米特插值示例
x = -pi:pi; y = sin(x); 
new_x = -pi:0.1:pi;
p = pchip(x,y,new_x);
figure(1); % 在同一个脚本文件里面，要想画多个图，需要给每个图编号，否则只会显示最后一个图哦~
plot(x, y, 'o', new_x, p, 'r-')

（2） 三次样条插值

注：三次样条插值的图像比埃尔米特插值更为平滑，但是使用要求条件更为苛刻。

要使用三次样条插值，必须满足：

代码：

代码示例:

% 三次样条插值和分段三次埃尔米特插值的对比
x = -pi:pi; 
y = sin(x); 
new_x = -pi:0.1:pi;
p1 = pchip(x,y,new_x);   %分段三次埃尔米特插值
p2 = spline(x,y,new_x);  %三次样条插值
figure(2);
plot(x,y,'o',new_x,p1,'r-',new_x,p2,'b-')
legend('样本点','三次埃尔米特插值','三次样条插值','Location','SouthEast')   %标注显示在东南方向
% 说明：
% LEGEND(string1,string2,string3, …)
% 分别将字符串1、字符串2、字符串3……标注到图中，每个字符串对应的图标为画图时的图标。
% ‘Location’用来指定标注显示的位置

(3)n维数据的插值（不常用）

（4）可以使用插值算法来预测数据

注：很少使用，且仅适用于短期预测

四.MATLAB画图

1.plot画图函数用法

plot(x1,y1,'样式1',x2,y2,'样式2')
线方式： - 实线 :点线 -. 虚点线 - - 波折线
点方式： . 圆点  +加号  * 星号  x x形  o 小圆
颜色： y黄； r红； g绿； b蓝； w白； k黑； m紫； c青

想让几组图同时显示在一个表中，可以依次往下写下去

样式表示方法：颜色+线方式+点方式，如果只需点方式，可省略线方式
如蓝色实线，对应值处用小圆圈可表示为

figure(1);
plot(x, y, 'b-o')

运行结果如图所示：

2.figure(n)

在同一个脚本文件里面，要想画多个图，需要给每个图编号,如figure(1)，figure(2)，figure(3)

代码

legend('string1','string2','string3','Location','方位（如SouthEast）')

例如：

plot(x,y,'o',new_x,p1,'r-',new_x,p2,'b-')
legend('样本点','三次埃尔米特插值','三次样条插值','Location','SouthEast')   %标注显示在东南方向

运行结果为：

注：标注方位也可以不写，可以用鼠标直接在生成的图片中拖动

3.横纵坐标标注

xlabel('X的值')%在横轴下方标注横轴的名称
ylabel('y的值')%在纵轴左侧标注纵轴的名称

运行结果为：