【正文】 y(2)=x1= y1=x2= y2=x3= y3=y(3)=x1= y1=x2= y2=x3= y3=y1 = 全區(qū)間拉格朗日插值 y2=lagrange(x,y,x0)y(1)=y(2)=y(3)=y2 = 6 討論與結(jié)論使用tic,toc函數(shù)計(jì)算下列四種方法計(jì)算上述問題所運(yùn)行的時(shí)間Functionlagrange(x0,y0,x)piece_linear(x0,y0,x)piece_square(x0,y0,x)運(yùn)行時(shí)間(s)從三次實(shí)驗(yàn)結(jié)果可知，三個(gè)程序的運(yùn)行時(shí)間都很短。 end作圖比較上圖為三種方法的插值曲線，由圖可得，三種曲線非常接近，這說明我們用拉格朗日插值計(jì)算所給點(diǎn)函數(shù)值的近似值時(shí)，引起的誤差還是比較小的。 x0=[ ]。,i)。 for s=l1:l+1 if s~=k a=a*(zx0(s))/(x0(k)x0(s))。m=length(x)。 l=0。分別通過分段線性插值、分段二次插值和全區(qū)間上拉格朗日插值計(jì)算f()，f()，f()的近似值。已知函數(shù)表如下：xf(x)3 算法原理與流程圖（1）原理設(shè)函數(shù)y=在插值區(qū)間[a,b]上連續(xù)，且在n+1個(gè)不同的插值節(jié)點(diǎn)a≤x0,x1,…,xn≤b上分別取值y0,y1,…,yn。%給l賦初值,根據(jù)x的值確定l if zx0(1)|zx0(n) fprintf(39。% n，p，m分別為x0，y0，x長度if n~=p fprintf(39。 end end sum=sum+y0(k)*a。 break。 y2=lagrange(x,y,x0)y(1)=Error!x(2) is out of range!y2 = 首先輸入函數(shù)變及待求點(diǎn) x=[ ]。參考文獻(xiàn)[1] 易大義，沈云寶，李有法. 計(jì)算方法(