naspline.txt

来自「方便初学者！插值与拟合的matlab实现！ 可方便同学做课设！」· 文本 代码 · 共 16 行

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