拉格朗日插值實(shí)驗(yàn)報(bào)告-文庫(kù)吧在線文庫(kù)
【正文】 值多項(xiàng)式（對(duì)于全區(qū)間上的插值，n取函數(shù)表的長(zhǎng)度）（2） 流程圖 分段線性插值 分段二次插值 全區(qū)間拉格朗日插值4 程序代碼及注釋分段線性插值%分段線性插值function y=piece_linear(x0,y0,x)% x0，y0為已知點(diǎn)，x為待求點(diǎn)n=length(x0)。 l=0。 for s=l1:l if s~=k a=a*(zx0(s))/(x0(k)x0(s))。m=length(x)。Error!x(%d) is out of range!\n39。 for s=l1:l+1 if s~=k a=a*(zx0(s))/(x0(k)x0(s))。m=length(x)。,i)。,i,y(i))。 x0=[ ]。m=length(x)。 end作圖比較上圖為三種方法的插值曲線，由圖可得，三種曲線非常接近，這說(shuō)明我們用拉格朗日插值計(jì)算所給點(diǎn)函數(shù)值的近似值時(shí)，引起的誤差還是比較小的。if zx0(1)|zx0(n) fprintf(39。注：分段線性插值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ì)算上述問(wèn)題所運(yùn)行的時(shí)間Functio
職業(yè)教育相關(guān)推薦
鄂ICP備17016276號(hào)-1