52牛顿插值.ppt

  1. 1、本文档共16页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
* 24 * 24 * 24 * 24 第二节、牛顿插值 Newton插值多项式的构造 差商的定义及性质 Newton插值方法 算法设计 Lagrange 插值虽然易算,但若要增加一个节点时,全部基函数 li(x) 都需要重新计算。 能否重新在Pn中寻找新的基函数 ? 希望每加一个节点时,只附加一项上去即可。 能否构成Pn的一组基函数 利用插值条件,代入上式,得关于Ak的线性代数方程组 引例 设一个点    唯一确实0次多项式为 代入可得 设二个点       多项式为 第一个点代入,后一项为零,即为一点式 第二个点代入,得 1、2个点的牛顿多项式 三个点: 前两个点代入得,第三项恒为零,即为二点式 第3个点代入得 3点牛顿多项式 General Form 差商(亦称均差) /* divided difference */ 0阶差商 1阶差商 2阶差商 K 阶差商 General Form where Rewriting 利用插值条件和差商,可求出Nn(x)的系数 bi 差商表 xk f ( xk ) 一阶差商 二阶差商 三阶差商 n 阶差商 差商表 例 x 1 2 3 4 5 6 y -3 0 15 48 105 192 根据下表数据,求牛顿插值多项式 xk f [xk] 一阶 二阶 三阶 四阶 五阶 1 2 0 3 15 4 48 5 105 6 192 3 15 33 57 87 6 9 12 15 1 1 1 0 0 0 解:构造差商表 -3 例1 :给定f(x)=lnx的数据表 xi 2.20 2.40 2.60 2.80 3.00 f(xi) 0.78846 0.87547 0.95551 1.02962 1.09861 解:差商表 N2(x)=0.78846 +0.43505(x-2.20) - 0.087375(x-2.20) (x-2.40) N4(x)= 0.78846 +0.43505(x-2.20) - 0.087375(x-2.20)(x-2.40) +0.0225(x-2.20)(x-2.40)(x-2.60) -0.00755(x-2.20)(x-2.40)(x-2.60)(x-2.80) Largran Code X=[1 2 3 4 5 6]; %X is a vector that contains a list of abscissas Y=X.^3-4*X; % Y is a vector that contains a list of ordinates n=length(X); D=zeros(n,n); D(:,1)=Y; %Use fozmula (20) to form the divided-difference table for j=2:n for k=j:n D(k,j)=(D(k,j-1)-D(k-1,j-1))/(X(k)-X(k-j+1)); end end %Determine the coefficients of the Newton interpolating polynomial C=D(n,n); for k=(n-1):-1:1 C=conv(C,poly(X(k))); m=length(C); C(m)=C(m)+D(k,k); end * 24 * 24 * 24 * 24

文档评论(0)

canjiadyy + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档