数值分析实验报告matlab精选 .doc

数值分析实验报告matlab 一、实验内容 设区间[-1,1]上函数 2in f? 11?25x 2 ，考虑区间[-1,1]的一个等距划分，分点为 n xi??1? ，i=0,1,2,...,n，则拉格朗日插值多项式为Ln ? ? i?0 11?25xi 2 li .其中， li，i=0,1,2,...,n是Lagrange插值基函数. 1) 选择不断增大的分点数目n=2,3,...，画出原函数f及插值多项式函数Ln 在[-1,1]上的图像，比较并分析实验结果. 2) 选择其他的函数，例如定义在区间[-5,5]上的函数h x1?x 4 ，g ?arctanx ， 重复上述的实验看其结果如何. 二、实验程序 function chapter2 promps={请选择试验函数，若选f，请输入f，若选好h，请输入h，若选g，请输入g：}; result=inputdlg; Nb_f=char; if errordlg; return; end result=inputdlg; Nd=str2num); if errordlg; return; end switch Nb_f casef f=inline);a=-1;b=1; caseh f=inline);a=-5;b=5; caseg f=inline);a=-5;b=5; end x0=linspace;y0=feval; x=a::b;y=Lagrange; clf; fplot; hold on; plot; xlabel;ylabel x and y=Ln --); 函数 function y=Lagrange n=length; m=length; for i=1:m z=x; s=0; for k=1:n p=; for j=1:n if p=p.*)/-x0); end end s=s+p*y0; end y=s; end 三、实验结果及分析 1)选择不断增大的分点数目n，原函数f及插值多项式函数Ln在[-1,1]上的 图像。 1 随着提高插值多项式次数，可以提高逼近的精度，但是次数的增加，在区间两端点附近与原函数偏离很远，即出现了Runge现象。 2) 选择不断增大的分点数目n，原函数h及插值多项式函数Ln在[-1,1]上的图像。 2 选择不断增大的分点数目n，原函数g及插值多项式函数Ln在[-1,1]上的图像。 3 同样，随着提高插值多项式次数，可以提高逼近的精度，但是次数的增加，在区间两端点附近与原函数偏离很远，即出现了Runge现象。 4 学 生 实 验 报 告 实验课程名称 数值分析 开课实验室 数学与统计学院实验室 学 院 数学与应用数学专业班 学 生 姓 名 学 号 开 课 时 间 2012 至 2013 学年第 一 学期 1 2 3 4 实验 多项式插值的震荡现象 问题提出：考虑在一个固定的区间上用插值逼近一个函数。显然Lagrange插值中使 用的节点越多，插值多项式的次数越高，我们自然关心插值多项式的次数增加时，Ln是否也更加靠近被逼近的函数。Runge 给出的一个例子是极著名并富有启发性的。设区间[-1,1]上函数 f? 11?25x2in 2 . 实验内容：考虑空间[-1,1]的一个等距划分，分点为 xi??1?则拉格朗日插值多项式为 n ， i?0，1,2 ...,n, Ln? ?1?25x i?0 1 2i . 其中，li,i?0,1,2,...,n是n次Lagrange插值基函数。 实验要求： 选择不断增大的分点数目n?2,3,...,画出原函数f及插值多项式函数Ln在[-1,1]上的图像，比较并分析实验结果。 选择其他的函数，例如定义在区间[-5,5]上的函数 h? x1?x 4 , g?arctan， 重复上述的实验看其结果如何。 首先编写拉格朗日插值函数的Matlab实现： Matlab程序为： function y=lagrange