MATLAB牛顿插值法例题与程序.pdfVIP

题 目 一 ： 多 项 式 插 值 某气象观测站在 8：00 （AM）开始每隔 10 分钟对天气作如下观测，用三次多项式插值 函数（Newton）逼近如下曲线，插值节点数据如上表， 并求出 9 点 30 分该地区的温度 （x=10）。 x 1 2 3 4 5 6 7 8 y 22.5 23.3 24.4 21.70 25.2 28.5 24.8 25.4 二、数学原理 假设有 n+1 个不同的节点及函数在节点上的值（ x 0 ，y 0 ），……（x n ，y n ），插值多 项式有如下形式： （ ） （ ） （ ） （ ）（ ） （ ） P x x - x x - x (x x ) x - x x x x x n 0 1 0 2 0 1 n 0 1 n （1） 其中系数 （i=0,1,2 ……n）为特定系数，可由插值样条 P （x ） y （i=0,1,2 ……n） i n i i 确定。 根据均差的定义，把 x 看成 [a,b] 上的一点 , 可得 f(x)= f （ x ）+f[ x ，x ] （x - x ） 0 0 1 0 f[x, x ]= f[ x ，x ]+f[x, x ， x ] （x - x ） 0 0 1 0 1 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 ]+f[ x ，x ] （x - x ）+ f[ x ，x ，x ] （x - x ）（x - x ）+ 0 0 1 0 0 1 2 0 1 …+ f[x, ] （ )+ f[x, , x ] （x ） （x ）+ （x ） x 0 ，…x n x - x 0 ）…(x-x n -1 x 0 ，…x n n 1 = Nn Rn 其中 N （x ）= f[ x ]+f