计算方法MATLAB实验报告Runge现象与方程求根合集.docxVIP

实验Runge现象的演示

【实验目的】

通过上机绘制Runge函数图像，理解高次插值的病态性质。

【实验内容】

在区间[-1,1]上取n=11等距节点对龙格(Runge)函数f(x)=1/(1+25*x^2)做多项式插值，并画出插值函数。

【实验程序】

functiontry_Runge

f=inline(1./(1+25*x.^2));

n=11;

X=linspace(-1,1,n);

Y=f(X);

x=-1:0.01:1;

y=mylagpoly(X,Y,x);

plot(x,f(x),r,X,Y,o,x,y,b)

title(Runge现象)

legend(y=1/(1+25*x^2),插值点,等分的10次插值多项式,0)

function[C,L]=interp_lag(X,Y)

n=length(X);

L=zeros(n,n);

forj=1:n

C1=[X(1:j-1),X(j+1:n)];

C2=poly(C1);