数值分析考试复习题 附录：matlab函数 .docVIP

附录：函数 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function p=newpoly(x,y,n,x0) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %p=newpoly(x,y,n) %(x,y)为节点坐标 %p为所求的的差值多项式的系数向量 %n为牛顿插值次数,n=1时为线性插值,n=2时为二次插值... N=length(x); d=[x;y]; if nargin==4 for i=1:N %根据插值误差估计式(根据距x0的远近)重新排列插值节点 for j=i+1:N if abs(d(1,i)-x0)abs(d(1,j)-x0) t=d(:,i); d(:,i)=d(:,j); d(:,j)=t; end end end end dq=zeros(n+1,n+2); dq(:,1)=d(1,1:n+1); dq(:,2)=d(2,1:n+1); for i=2:n+1 for j=3:i+1 dq(i,j)=(dq(i,j-1)-dq(i-1,j-