计算方法_Newton插值法(差商)

clear all; %函数插值——Newton插值

n=input('请输入插值节点个数（正整数）:');

for i=1:n

F(i,1)=input('输入节点值:');

F(i,2)=input('输入对应节点函数值值:');

end

for j=3:n+1

for i=j-1:n

F(i,j)=(F(i,j-1)-F(i-1,j-1))/(F(i,1)-F(i+2-j,1));%差商

end

end

F %矩阵F即为差商表

运行结果：

请输入插值节点个数（正整数）:4

输入节点值:>> -1.1

输入对应节点函数值值:-1.12

输入节点值:-0.5

输入对应节点函数值值:0

输入节点值:0.9

输入对应节点函数值值:1.8

输入节点值:1.7

输入对应节点函数值值:2.2

F =

-1.1000 -1.1200 0 0 0 -0.5000 0 1.8667 0 0

0.9000 1.8000 1.2857 -0.2905 0

1.7000

2.2000 0.5000 -0.3571 -0.0238