拉格朗日插值實(shí)驗(yàn)報(bào)告-閱讀頁
2024-08-08 19:03本頁面
【正文】 0(1)|zx0(n) fprintf(39。,i)。 end%當(dāng)插值點(diǎn)超出范圍時(shí)，報(bào)錯(cuò) for k=1:n p=。 end end s=p*y0(k)+s。fprintf(39。,i,y(i))。 y=[2 3 4]。 y=[2 3 4 5]。 y2=lagrange(x,y,x0)Error!x(1) is out of range! x=[1 2 3 4]。 x0=[ ]。 y=[ ]。注：分段線性插值y0=piece_linear(x,y,x0)y(1)=x1= y1=,x2= y2=y(2)=x1= y1=,x2= y2=y(3)=x1= y1=,x2= y2=y0 = 分段二次插值 y1=piece_square(x,y,x0)y(1)=x1= y1=x2= y2=x3= y3=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í)間都很短。所以加入以下判斷以提高插值的準(zhǔn)確性n=length(x0)。m=length(x)。Error! Please input again!\n39。if zx0(1)|zx0(n) fprintf(39。,i)。 end作圖比較上圖為三種方法的插值曲線，由圖可得，三種曲線非常接近，這說明我們用拉格朗日插值計(jì)算所給點(diǎn)函數(shù)值的近似值時(shí)，引起的誤差還是比較小的
職業(yè)教育相關(guān)推薦
鄂ICP備17016276號-1