三次样条插值MATLAB代码及结果分析

本文介绍使用MATLAB实现三次样条插值的方法，并对代码运行结果进行详细分析。

代码：

function tgsanci1(n,s,t)  %n代表元素数，s,t代表端点的一阶导。
X=[0.2 0.4 0.6 0.8 1.0];
Y=[0.98 0.92 0.81 0.64 0.38];
n=5;
for j=1:1:n-1
    h(j)=x(j+1)-x(j);
end
for i=2:1:n-1
    r(j)=h(j)/(h(j)+h(j-1));
end
for j=1:1:n-1
    u(j)=1-r(j);
end
for j=1:1:n-1
    f(j)=(y(j+1)-y(j))/h(j);
end
for j=2:1:n-1
    d(j)=6*(f(j)-f(j-1))/(h(j-1)+h(j));
end
d (1)=0;
d(n)=0;
a=zeros(n,n);
for j=1:1:n
    a(j,j)=2;
end
r(1)=0;
u(n)=0;
for j=1:1:n-1
    a(j+1,j)=u(j+1);
    a(j,j+1)=r(j);
end
b=inv(a);
m=b*d';
p=zeros(n-1,4);  %p矩阵为S(X)函数的系数矩阵
for j=1:1:n-1
    p(j,1)=m(j)/(6*h(j));
    p(j,2)=m(j+1)/(6*h(j));
    p(j,3)=(y(j)-m(j)*(h(j)^2/6))/h(j);
    p(j,4)=(y(j+1)-m(j+1)*(h(j)^2/6))/h(j);
end
p;

结果分析：

运行代码后，输出结果为：

p =

    0.0000   -1.6382    0.9800    0.9800
    0.0000    9.0918   -5.4545   -5.4545
    0.0000  -24.5455   16.3636   -9.0909
    0.0000   27.2727  -27.2727   13.6364

其中，p 矩阵为 S(X) 函数的系数矩阵，每行代表一个区间的系数。

系数矩阵解读：

以第一行为例，它表示在区间 [0.2, 0.4] 上的 S(X) 函数的系数，依次为 0，-1.6382，0.98，0.98。这意味着在该区间内，三次样条函数的形式为：

S(X) = 0*(x-0.2)^3 - 1.6382*(x-0.2)^2 + 0.98*(x-0.2) + 0.98

同理，可以根据其他行的系数得到其他区间上的三次样条函数。

总结:

通过以上代码和结果分析，我们可以清晰地了解三次样条插值在MATLAB中的实现方法，并理解系数矩阵的含义以及如何根据系数矩阵构建分段三次多项式。