matlab插值龙格实验,matlab实现Lagrange多项式插值观察龙格现象

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

-16.8552 0.0000 1.0000

在matlab窗口中输入：

x=linspace(-1,1,16);y=1./(1+25*x.^2);

lagran(x,y)

回车可得结果：

ans = 0.0000 123.3597