naspline.txt

方便初学者！插值与拟合的matlab实现！ 可方便同学做课设！

function m=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=nachase(lemda,2*ones(1:n-1),mu,g);

m=[dy0;dy;dyn];

if nargin>=5