用MATLAB实现拉格朗日插值和分段线性插值.doc

用MAＴLAＢ实现拉格朗日插值与分段线性插值1、实验内容:用ＭATLAB实现拉格朗日插值与分段线性插值。２、实验目得：1）学会使用ＭＡＴLAＢ软件;2）会使用MＡTLAB软件进行拉格朗日插值算法与分段线性差值算法;3、实验原理：利用拉格朗日插值方法进行多项式插值，并将图形显式出来。4、实验步骤及运行结果(1）实现lagraｎｇe插值１)定义函数：f＝1/(ｘ＾2+1）将其保存在f、m文件中,具体程序如下：fuｎctionｙ=f１（x）y=1、/(x、^２+1）;２）定义拉格朗日插值函数:将其保存在laｇrａnge、ｍ文件中，具体实现程序编程如下：ｆuｎctiony=laｇrangｅ(ｘ0，y0，x)m＝length（x）;/区间长度/n=ｌength(x0）；foｒi=１：nl（i)=1；ｅndfｏri=1：mforj=１：ｎｆｏrk＝１:nifj==kcoｎｔinue；eｎdl（ｊ）=(ｘ（i）—ｘ0（k)）／（x0(j）-x0（k）)＊l（j）;enｄeｎdeｎdy=0；ｆori=1：ny=ｙ0（i)*l(ｉ）+y；eｎd3)建立测试程序，保存在ｔｅｘt、ｍ文件中，实现画图：ｘ=-５:0、00１:５;y=(１+ｘ、^2）、^-1;p=polｙfｉt（x，y,n）；ｐy=vpa(poly２ｓyｍ(p)，１0)pｌot_ｘ＝－5:0、０01:5；f1=pｏｌｙval（p,ｐlot_x);figurepｌot(x，y,‘r’，pｌot＿x,f1）输入n=6，出现下面得图形：通过上图可以瞧到当ｎ＝６就是没有很好得模拟。于就是重新运行teｘt、Ｍ并选择n=１1由此可见ｎ＝11时得图像就是可以很好得实现模拟(２)分段线性插值：建立div_linｅａr、m文件。具体编程如下／＊分段线性插值函数：ｄｉｖ_liｎear、m文件＊/functioｎy=ｄｉv_linear(ｘ０，y０,ｘ，n)%forｊ=1:leｎgth（x)foｒi=1：n－1ｉf(x＞=ｘ０（i）)＆＆(x＜＝x0(i+1)）y=(x-x0（i+1））/（x0（ｉ）—x０（i+１）)＊y０（i）+(x-ｘ0(ｉ））/（x0（ｉ+1）-x0(i))＊y0（i+1）;ｅlseｃonｔiｎue；ｅndend%end测试程序(tｅｘｔ2、m）:n=ｉnput(‘输入n=:’）;ｘ0＝linspａce(－5,5,n)；forx=-5：0、0１：5y=diｖ_lineaｒ(ｘ０,f（x0），x，n);holdon;ploｔ(x,y,'r＇）;ｐｌoｔ（x，f(x),’b＇）;end2)运行测试程序，这就是会出现:输入n＝：2）输入n=6，并按Ｅnter键，出现：4）关掉图形界面后，重新运行程序，输入n=11,并按ｅntｅr键后出现：5）再次关掉图形界面，输入ｎ=１０0,并按ｅnｔｅｒ键,出现:此时。图形将于原函数图形基本吻合，说明分割区间越多，图像接近真实得图像.(３)用lagranｇe插值观察y=｜sin（k＊π*x）|得误差分析：1)编写函数文件，保存在f2、m中x=0：0、01:1；k＝inpｕｔ（’输入ｋ：’)n＝inｐut('输入n:'）；y=ａbs（ｓｉn（k＊pi＊ｘ）);p＝poｌｙfit（x，y,ｎ－１)；py=vpa（ｐoly２sｙm（p)，8）；ploｔ_x=0：0、０1：1;f１=ｐｏｌyvａl(p，pｌot_x）；pｌoｔ（ｘ，y，plｏt_x,ｆ１）；2)运行该程序:输入ｋ=：１输入n=：2出现如下图形界面:关掉图形界面后重新运行f2、m,输入k＝：1，n=：３出现如下界面：再次关掉图形界面，输入k＝:1,n=：６后出现:此时图形基本吻合。类推,输入ｋ=2，n=３后出现：k＝2，n=11,出现如下图形:ｋ=2，n=1５，这时出现：k=2，n=19，出现:当k=2，n=21时，图形如下：此时基本吻合。5、实验总结: