Newton插值的matlab实现.docVIP

Newton插值的matlab实现 成员： 黄 全 P092314746 28% 付 吉 P091712737 24% 颜学俭 P091712716 24% 罗国庭 P091712739 24% 指导老师：刘华 2012年5月2日 目录 Newton插值的matlab实现 - 1 - 一 过程整理 - 3 - Newton插值的基本原理 - 3 - 二 流程图 - 4 - 三 算法设计 - 6 - 3.1、Newton插值的matlab实现 - 6 - 3.2、程序 - 6 - 3.3、例题 - 6 - 3．4、命令执行图 - 7 - 四 参考文献 - 7 - 一 过程整理 Newton插值的基本原理 假设有n+1个不同的节点及函数在节点上的值（x，y），……（x，y），插值多项式有如下形式： （1） 其中系数（i=0,1,2……n）为特定系数，可由插值样条（i=0,1,2……n）确定。 根据均差的定义，把x看成[a,b]上的一点,可得 f(x)= f（）+f[]（） f[x, ]= f[]+f[x,] （） …… f[x, ，…x]= f[x, ，…x]+ f[x, ，…x](x-x) 综合以上式子，把后一式代入前一式，可得到： f(x)= f（）+f[]（）+ f[]（）（）+ …+ f[x, ，…x]（）…(x-x)+ f[x, ，…x]= N（x）+其中 N（x）= f（）+f[]（）+ f[]（）（）+ …+ f[x, ，…x]（）…(x-x)（2） = f(x)- N（x）= f[x, ，…x]（3） =（）…(x-x) Newton插值的系数（i=0,1,2……n）可以用差商表示。一般有 [] （k=0,1,2，……，n ）（4） 把（4）代入（1）得到满足插值条件N（i=0,1,2，……n）的n次Newton插值多项式 N（x）=f（）+f[]（）+f[]（）（）+……+f[]（）（）…（）. 其中插值余项为： 介于之间。 二 流程图 三 算法设计 3.1、Newton插值的matlab实现 构造计算经过点（）=（k=0,1,2，……，n ）的次数小于等于n的牛顿插值多项式： N（x）=f（）+f[]（）+f[]（）（）+……+f[]（）（）…（） 3.2、程序 function[c,d]=newploy(x,y) n=length(x); d=zeros(n,n); d(:,1)=y; 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 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 3.3、例题 用程序求已知数据点（1，-3），（2，0），（3，15），（4，48），（5，105），（6，162）的牛顿插值多项式。 在命令输入： newpoly([1 2 3 4 5 6 ],[-3,0,15,48,105,192]) 回车得到结果： ans = 0 0 1 0 -4 0 由此得出牛顿插值多项式为： 3．4、命令执行图 例如建立M-文件： t=0:0.1:6; z=t.^3-4*t; x=[1 2 3 4 5 6]; y=[-3,0,15,48,105,192]; plot(t,z,x,y,ro);grid 命令执行得到如下图形： 牛顿插值与多项式数据点之间的关系 四 参考文献 《MATLAB程序设计》 王建卫 取中水 凌宾 《MATLAB数值分析与应用》 张德丰 《数值分析》 黄明游 刘播 徐涛 - 1 - 开 始 int s,int t t=s+1 N Y f(s,t)=(d[t].y-d[s].y)/(d[t].x-d[s].x) f(s,t)=f(s+1,t)-f(s,t-1))/(d[t].x-d[s].x) 输入插值次数n float t=1.0, y=d[0].y, yt=0.0,int j=1 y=(x-d[j-1].x)*t+f(0,j)*t 返回y j+1=j Y j=n N 结 束