插值法例题和作业

一、埃尔米特插值pchip

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-') %画出原来的和信生成的

 

 

二、三次样条插值spline

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') %标注显示在东南方

 

三、导入数据