2011 学年第一学期数值逼近数值逼近上机实验报告上机实验报告院系：文正学院 专业：信息与计算科学姓名：许志杰 学号：0917417001年级：2009 实验场地及名称：文正学院机房实验时间：9 月 26 日实验成绩:指导教师：陈旻昕上机实验一（上机实验一（Runge 现象）现象）在-5,5区间上，作图显示函数的四阶，六阶，十阶 Lagrange 多项211)(xxf式插值。答：程序如下：1、ex1.m：clearclose allx=-5:0.1:5y=1./(1+x.2);plot(x,y)xi=-5,-3,0,3,5;yi=1./(1+xi.2);xx=-5:0.1:5;yy=langrage(xi,yi,xx);hold onplot(xx,yy,-r)xi=-5,-3,-1,0,1,3,5;yi=1./(1+xi.2);xx=-5:0.1:5;yy=langrage(xi,yi,xx);hold onplot(xx,yy,-g)xi=-5,-4,-3,-2,-1,0,1,2,3,4,5;yi=1./(1+xi.2);xx=-5:0.1:5;yy=langrage(xi,yi,xx);hold onplot(xx,yy,-k)2、f.mfunction y=f(x)y=1./(1+x.2);3、Langrage.mfunction yy=langrage(xi,yi,xx)yy=0for k=1:length(xi)s=1;for i1=1:length(xi)if(i1=k)s=s.*(xx-xi(i1)./(xi(k)-xi(i1)endendyy=yy+yi(k)*send程序运行结果：原函数图像原函数图像原函数图像原函数图像+四阶插值图像四阶插值图像原函数图像原函数图像+四阶插值四阶插值+六阶插值图像六阶插值图像原函数图像原函数图像+四阶插值四阶插值+六阶插值六阶插值+十阶插值图像十阶插值图像