MAAB牛顿插值法例题与程序.docx

题 目 一 ： 多 项 式 插 值

某气象观测站在8：00（AM）开始每隔10分钟对天气作如下观测，用三次多项式插值

x

1

2

3

4

5

6

7

8

y

二、数学原理

假设有n+1个不同的节点及函数在节点上的值（x，y），……（x，y ），插值多

0 0 n n

项式有如下形式：

P（x）?? ??（x-x）??（x-x）(x?x)?????（x-x）（x?x）??（x?x）（1）

n 0 1 0 2 0 1 n 0 1 n

其中系数?

i

确定。

（i=0,1,2……n）为特定系数，可由插值样条P（x）?y

n i i

（i=0,1,2……n）

根据均差的定义，把x看成[a,b]上的一点,可得

f(x)=f（x

0

）+f[x，x

0 1

]（x-x ）

0

f[x,x

0

]=f[x，x

0 1

]+f[x,x，x

0 1

……

]（x-x）

1

f[x,x

0

，…x

n-1

]=f[x,x

0

，…x

n

]+f[x,x

0

，…x

n

](x-x)

n

综合以上式子，把后一式代入前一式，可得到：

f(x)=f[x

0

]+f[x，x

0 1

]（x-x

0

）+f[x

，x，x

0 1 2

]（x-x

0

）（x-x）+

1

…+f[x,x ，…x

]（x-x ）…(x-x

)+f[x,x ，…x,x]?

（x）=N

（x）+R（x）

0

n

0

n-1

0

n n

?1 n

n

其中

N（x）=f[x

n 0

]+f[x，x

0 1

]（x-x

0

）+f[x

，x，x

0 1 2

]（x-x

0

）（x-x）+

1

…+f[x,x

0

，…x

n

]（x-x

0

）…(x-x

n-1

)（2）

R（x）=f(x)-N

n n

（x）=f[x,x

0

，…x

n

,x]?

（x）（3）

n?1

?（x）=（x-x ）…(x-x)

n?1 0 n

Newton插值的系数?

i

（i=0,1,2……n）可以用差商表示。一般有

? ?f[x，x

k 0 1

??x

k

]（k=0,1,2，……，n）（4）

把（4）代入（1）得到满足插值条件N（x）?f（x）（i=0,1,2，……n）的n次Newton

插值多项式

n i i

N（x）=f（x

n 0

）+f[x，x

0 1

]（x-x

1

）+f[x

，x，x

0 1 2

]（x -x

1

）（x -x

2

）+……

+f[x，x

0 1

??x

n

]（x -x

1

）（x -x

2

）…（x-x

）.

n-1

其中插值余项为：

?介于x，x

0 1

??x

k

之间。

三、程序设计

function[y,A,C,L]=newdscg(X,Y,x,M)

%y为对应x的值，A为差商表，C为多项式系数，L为多项式

%X为给定节点，Y为节点值，x为待求节点n=length(X);m=length(x);%n为X的长度fort=1:mz=x(t);A=zeros(n,n);A(:,1)=Y;s=;p=;q1=;c1=;

forj=2:nfori=j:n

A(i,j)=(A(i,j-1)-A(i-1,j-1))/(X(i)-X(i-j+1));

end

q1=abs(q1*(z-X(j-1)));c1=c1*j;end

C=A(n,n);q1=abs(q1*(z-X(n)));fork=(n-1):-1:1

C=conv(C,poly(X(k)));

d=length(C);C(d)=C(d)+A(k,k);end

y(k)=polyval(C,z);%输出y值end

L(k,:)=poly2sym(C);%输出多项式

symsM,X=[1,3,5,7];Y=[,,,];x=10;

[y,A,C,L]=newdscg(X,Y,x,M)y=

A=

四、结果分析和讨论

对于不超过三次的插值多项式，x如果选取1,3,5,7这三个点能够得到较好的三次插值多项式L=^^2++。当x=10时，也即9点30分时的温度为度，结果分析知此值应是偏小的。对于选取不同的插值节点，能够得到不同的插值多项式，误差也不尽相同。

五、完成题目的体会与收获

对于牛顿插值法有了更深的了解，合理选择插值节点很重要。加深了对其原理的认识