用matlab计算牛顿插值的例题,MATLAB牛顿插值法例题与程序.doc

MATLAB牛顿插值法例题与程序

题目一：多项式插值

某气象观测站在8：00(AM)开始每隔10分钟对天气作如下观测，用三次多项式插值函数(Newton)逼近如下曲线，插值节点数据如上表，并求出9点30分该地区的温度(x=10)。

x12345678y22.523.324.421.7025.228.524.825.4二、数学原理

假设有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)

=(