1. 拉格朗日多项式插值
了解概念
插值多项式
插值节点
范德蒙特(Vandermonde)行列式
截断误差、插值余项
特点
函数实现
function y=lagrange(x0,y0,x)
n=length(x0);m=length(x);
for i=1:m
z=x(i);
s=0.0;
for k=1:n
p=1.0;
for j=1:n
if j~=k
p=p*(z-x0(j))/(x0(k)-x0(j));
end
end
s=p*y0(k)+s;
end
y(i)=s;
end
设n个节点数据以数组x0,y0输入(注意Matlat的数组下标从1开始),m个插值点以数组x 输入,输出数组y为m个插值。
则可用y = lagrange(x0,y0,x)调用。
2. 牛顿(Newton)插值
了解概念
差商
差分
等距节点插值公式(Newton向前插值公式)
特点
每增加一个节点,插值多项式只增加一项,因而便于递推运算。而且 Newton 插值的计算量小于Lagrange 插值。
函数实现
3. 分段线性插值
了解概念
插值多项式的振荡
特点
将每两个相邻的节点用直线连起来,如此形成的一条折线就是分段线性插值函数。它是为了解决高次插