l1.m

来自「计算插值与曲线拟合的MATLAB程序」· M 代码 · 共 10 行

clear
xn=[-5:0.05:5]; %插值点（多个）
fn=1./(1+xn.^2); %函数的准确值
x1=linspace(-5,5,11) ;    y1=1./(1+x1.^2);  %11个样本点 x1=[-5  -4  -3  -2  -1  0  1  2  3  4  5]；
L10=lagrange (x1,y1,xn);  % 用11个样本点的Lagrange插值结果

x2= linspace(-5,5,7) ; y2=1./(1+x2.^2); %用7个样本点
L6=lagrange (x2,y2,xn);  % 用7个样本点的Lagrange插值结果

plot(xn,L10,'r',xn,L6, '*g',xn,fn,'-b',x1,y1,'o', x2,y2,'o')