第1章_插值.pptVIP

结束 * 结束 * Matlab程序代码如下： x0=[x1, … xn+1]；y0=[y1,…, yn+1]；x=?； n=length(x0);Q=zeros(n); %用于储存差商表 Q(:,1)=y0’ ; %0阶差商，即函数值 for i=2:n for j=2:i Q(i,j)=(Q(i,j-1)-Q(i-1,j-1))/(x(i)-x(i-j+1)); end end b(n)=Q(n,n) for k=n:-1:2 %另一种迭代 b(k-1)=Q(k-1,k-1)+b(k)*(a-x(k-1)); end b(1), Q 结束 当Mn+1随着n的增大变化不是太大时,|Rn(x)|会随n的增大而减小,但当Mn+1也随n的增大而增大时,就很难保证|Rn(x)|随n的增大而减小了. 先看一个典型的例子： §1.5 分 段 插 值 在一般情况下,似乎可以靠增加插值节点的数目来改善插值的精度，但是否在任何情况下,插值节点的数目越多,插值多项式与被插函数的逼近程度就越好呢?答案是否定的.事实上,从插值公式的余项可以看出： 若在［-1,1］上取等距节点x