p2.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;
mn=100;
Xnew=X*cos(pi/4)+Y*sin(pi/4);
Ynew=Y*cos(pi/4)-X*sin(pi/4);
m0new=(-sin(pi/4)+cos(pi/4)*m0)/(cos(pi/4)+sin(pi/4)*m0);
mnnew=(-sin(pi/4)+cos(pi/4)*mn)/(cos(pi/4)+sin(pi/4)*mn);
N=length(Xnew);
M=zeros(1,N);
x=[Xnew(1):0.01:Xnew(N)];
y=cubicspline(Xnew,Ynew,m0new,mnnew,x);
plot(x,y,'r')
hold on
plot(Xnew,Ynew,'o')