Matlab进行Lagrange多项式插值
拉格朗日插值法对函数y=1./(1+25*x.^2)在区间[-1,1]进行5次、10次、15次插值观察龙格现象
主程序
1.拉格朗日
function [c,l]=lagran(x,y)
%c为多项式函数输出的系数
%l为矩阵的系数多项式
%x为横坐标上的坐标向量
%y为纵坐标上的坐标向量
w=length(x);
n=w-1;
l=zeros(w,w);
for k=1:n+1
v=1;
for j=1:n+1
if k~=j
v=conv(v,poly(x(j)))/(x(k)-x(j)) %对多项式做卷积运算
end
end
l(k,:)=v;
end
c=y*l;
2.在matlab窗口中输入:
x=linspace(-1,1,6);y=1./(1+25*x.^2);
lagran(x,y)
回车可得结果:
ans =
-0.0000 1.2019 -0.0000 -1.7308 -0.0000 0.5673
在matlab窗口中输入:
x=linspace(-1,1,11);y=1./(1+25*x.^2);
lagran(x,y)
回车可得结果:
ans =
-220.9417 0.0000 494.9095 -0.0000 -381.4338 -0.0000 123.3597 0.0000 -16.8552 0.0000 1.0000
在matlab窗口中输入:
x=linspace(-1,1,16);y=1./(1+25*x.^2);
lagran(x,y)
回车可得结果:
ans =