MATLAB作业拉格朗日三阶样条插值函数牛顿插值.pdfVIP

Lagrange插值M函数

symsx

x0=[0,1,2];

y0=[1,2,3];

n=length(x0);

fori=1:n

a=1;

forj=1:n

ifj~=i

a=expand(a*(x-x0(j)));

end

end

b=1;

fork=1:n

ifk~=i

b=b*(x0(i)-x0(k));

end

end

A(i)=expand(a/b);

end

L=0;

forp=1:n

L=L+y0(p)*A(p);

end

L

Language

L=x+1

三阶样条插值M函数

functionm=naspline(x,y,dy0,dyn,xx)

n=length(x)-1;

h=diff(x);lemda=h(2/n)./(h(1:n-1)+h(2:n));mu=1-lemda;

g=3*(lemda.*diff(y(1:n))./h(1:n-1)+mu.*diff(y(2:n+1))./h(2:n));

g(1)=g(1)-lemda(1)*dy0;g(n-1)=g(n-1)-mu(n-1)*dyn;

dy=nac