拉格朗日插值實(shí)驗(yàn)報(bào)告-文庫吧
2025-07-05 19:03 本頁面
【正文】 于是我們可以得出如下的拉格朗日n次插值多項(xiàng)式（對(duì)于全區(qū)間上的插值，n取函數(shù)表的長度）（2） 流程圖 分段線性插值 分段二次插值 全區(qū)間拉格朗日插值4 程序代碼及注釋分段線性插值%分段線性插值function y=piece_linear(x0,y0,x)% x0，y0為已知點(diǎn)，x為待求點(diǎn)n=length(x0)。p=length(y0)。m=length(x)。% n，p，m分別為x0，y0，x長度if n~=p fprintf(39。Error! Please input again!\n39。)。% x0和y0長度不等時(shí)，報(bào)錯(cuò)elsefor i=1:m z=x(i)。 sum=。 l=0。%給l賦初值,根據(jù)x的值確定l if zx0(1)|zx0(n) fprintf(39。Error!x(%d) is out of range!\n39。,i)。 break。 end%當(dāng)插值點(diǎn)超出范圍時(shí)，報(bào)錯(cuò) for j=2:n if zx0(j) l=j。 end if l~=0 break。 end end%一旦l有非零值，則終止循環(huán),選出合適的l for k=l1:l a=。 for s=l1:l if s~=k a=a*(zx0(s))/(x0(k)x0(s))。 end end sum=sum+y0(k)*a。 end y(i)=sum。 fprintf(39。y(%d)=%f\nx1=%.3f y1=%.5f,x2=%.3f y2=%.5f\n\n39。,i,y(i),x0(l1),y0(l1),x0(l),y0(l))。%輸出插值結(jié)果和所需節(jié)點(diǎn)endendend分段二次插
職業(yè)教育相關(guān)推薦
鄂ICP備17016276號(hào)-1