关于用追赶法求解系数矩阵的三次样条插值实验报告。.docVIP

数值分析实验报告 实验名称 关于用追赶法求解系数矩阵的三次样条插值实验报告。 实验目的 会使用matlab 语言编程使用追赶法和三次样条插值 结合求解线性方程组的解。 实验内容 算法如下： function?x=followup(a,b,c,d) n=length(d); a(1)=0; %“追”的过程 L(1)=b(1); y(1)=d(1)/L(1); u(1)=c(1)/L(1); for?i=2:(n-1) ????L(i)=b(i)-a(i)*u(i-1); ????y(i)=(d(i)-y(i-1)*a(i))/L(i); ????u(i)=c(i)/L(i); end L(n)=b(n)-a(n)*u(n-1); y(n)=(d(n)-y(n-1)*a(n))/L(n); %“赶”的过程 x(n)=y(n); for?i=(n-1):-1:1 ????x(i)=y(i)-u(i)*x(i+1); end用matlab编程如下: function[s,y0]=spline3?(x,y,x0) %x,y为数表x0为插值点s表示插值函