demohermint.m

可以方便的调用这个函数

%函数demoHermite使用hermint函数对y=xe^(-x)进行分段三阶插值
function demoHermint(n)
if nargin<1,n=4;end
x=linspace(0,8,n);
y=x.*exp(-x);
yp=(1-x).*exp(-x);
xi=linspace(min(x),max(x));
ye=xi.*exp(-xi);
yi=hermint(x,y,yp,xi);
err=norm(yi-ye(:));
fprintf('error=%12.2e with %d knots\n',err,n);
plot(x,y,'bo',xi,ye,'b-',xi,yi,'r--');
legend('Given','x*exp(-x)','Hermite');
text(4,0.2,sprintf('%d knots',n));
axis([0 8 0 0.5]);