拉格朗日插值實驗報告(已改無錯字)
2022-08-17 19:03:10 本頁面
【正文】 區(qū)間插值%拉格朗日全區(qū)間插值function y=lagrange(x0,y0,x)% x0，y0為已知點，x為待求點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長度不等時，報錯elsefor i=1:m z=x(i)。 s=。 if zx0(1)|zx0(n) fprintf(39。Error!x(%d) is out of range!\n39。,i)。 break。 end%當插值點超出范圍時，報錯 for k=1:n p=。 for j=1:n if j~=k p=p*(zx0(j))/(x0(k)x0(j))。 end end s=p*y0(k)+s。 end y(i)=s。fprintf(39。y(%d)=%.5f\n39。,i,y(i))。%輸出插值結(jié)果endendend5 算例分析 測試示例 x=[1 2 3 4]。 y=[2 3 4]。 y2=lagrange(x,y,x0)Error! Please input again! x=[1 2 3 4]。 y=[2 3 4 5]。 x0=[ ]。 y2=lagrange(x,y,x0)Error!x(1) is out of range! x=[1 2 3 4]。 y=[2 3 4 5]。 x0=[ ]。 y2=lagrange(x,y,x0)y(1)=Error!x(2) is out of range!y2 = 首先輸入函數(shù)變及待求點 x=[ ]。 y=[ ]。 x0=[ ]。注：分段線性插值y0=piece_linear(x,y,x0)y(1)=x1= y1=
職業(yè)教育相關(guān)推薦
鄂ICP備17016276號-1