MATLAB龙格现象插值实验报告含代码图像结果可直接套用.docxVIP

实验龙格现象的发生、防止和插值效果的比较

【实验目的】

1．掌握用MATLAB编程语言的编程训练，初步体验算法的软件实现；

2．对不同方案结果的精度比较，有无龙格现象的发生。

【实验内容】

将区间【-5,5】等分，xi=-5+ih(i=0,1,…,n),h=10/n，对下列函数进行插值，计算其在xk=-5+0.25k(k=0,1,…,40)上的值，并绘制插值函数的图形。

（1）y=x/(1+x^4)(2)y=arctgx

(1)做n=10,20拉格朗日插值；

(2)做n=10,20分段线性插值；

(3)做n=10,20I型三次样条插值

【程序如下】:

（1）y=x/(1+x^4)

clearall;

closeall;

clc;

N=input(输入插值次数);

xi=zeros(1,N+1);

h=10/N;

fork=1:N+1

xi(k)=-5+(k-1)*h;

end

yi=xi./(1+xi.^4);

xk=zeros(1,41);

fork=1:41

xk(k)=-5+0.25*(k-1)

end

ykl=zeros(1,41);

fori=1:41

ykl(i)=0;

fork=1:N+1

t=1;

forj=1:N+1

ifj~=k

t=(xk(i)-xi(j))/(xi(k)-xi(j))*t;