数值分析用MATLAB实验报告,数值分析MATLAB实验报告

实验 2.1 多项式插值的震荡现象

问题提出：考虑在一个固定的区间上用插值逼近一个函数。显然Lagrange插值中使

用的节点越多，插值多项式的次数越高，我们自然关心插值多项式的次数增加时，Ln(x)是否也更加靠近被逼近的函数。Runge 给出的一个例子是极著名并富有启发性的。设区间[-1,1]上函数

f(x)

11 25x2in

2

.

实验内容：考虑空间[-1,1]的一个等距划分，分点为 xi 1 则拉格朗日插值多项式为

n

， i 0，1,2 ...,n,

Ln(x)

1 25x

i 0

1

2i

(x).

其中，li(x),i 0,1,2,...,n是n次Lagrange插值基函数。

实验要求：

(1)选择不断增大的分点数目n 2,3,...,画出原函数f(x)及插值多项式函数Ln(x)在[-1,1]上的图像，比较并分析实验结果。

(2)选择其他的函数，例如定义在区间[-5,5]上的函数 h(x)

x1 x

4

, g(x) arctan(x)，

重复上述的实验看其结果如何。

首先编写拉格朗日插值函数的Matlab实现： Matlab程序为：

function y=lagrange(x0,y0,x) %Lagrange插值 n=length(x0); m=length(x); for i=1:m z=x(i); s=0.0; for k=1:n p=1.0; for j=1:n if(j~=k)

p=p*(z-x0(j))/(x0(k)-x0(j)); end end

s=s+p*y0(k); end