p1.m

自行编写的三次样条插值函数源代码

clc;
clear;
X=[8.125    8.4      9.0      9.485    9.6    9.959    10.166     10.2];
Y=[0.0774   0.099    0.280    0.60     0.708   1.200     1.800     2.177];
m0=0.01087;  %Y0'
mn=100; %Yn'
N=length(X);
M=zeros(1,N); %初始化
x=[8.125:0.001:10.2];
y=cubicspline(X,Y,m0,mn,x);
plot(x,y,'-r')
hold on
plot(X,Y,'o')