拉格朗日插值實(shí)驗(yàn)報(bào)告-預(yù)覽頁
【正文】 的插值函數(shù)類Φ中，求一簡單函數(shù)P(x)，滿足插值條件P(xi)=yi(i=0,1,…,n)，而在其他點(diǎn)x≠xi上，作為f(x)近似值。類似地，我們可以選取距離x最近的三個(gè)節(jié)點(diǎn)xi1，xi與xi+1，然后進(jìn)行二次插值，即得這種分段低次插值叫分段二次插值，又稱分段拋物線插值。m=length(x)。% x0和y0長度不等時(shí)，報(bào)錯(cuò)elsefor i=1:m z=x(i)。Error!x(%d) is out of range!\n39。 end if l~=0 break。 end y(i)=sum。%輸出插值結(jié)果和所需節(jié)點(diǎn)endendend分段二次插值%分段二次插值function y=piece_square(x0,y0,x)% x0，y0為已知點(diǎn)，x為待求點(diǎn)n=length(x0)。Error! Please input again!\n39。 l=0。 break。 end%一旦l有非零值，則終止循環(huán),選出合適的l end if l==0 l=n1。 end y(i)=sum。%輸出插值結(jié)果與所需節(jié)點(diǎn)endendend拉格朗日全區(qū)間插值%拉格朗日全區(qū)間插值function y=lagrange(x0,y0,x)% x0，y0為已知點(diǎn)，x為待求點(diǎn)n=length(x0)。Error! Please input again!\n39。 if zx0(1)|zx0(n) fprintf(39。 end%當(dāng)插值點(diǎn)超出范圍時(shí)，報(bào)錯(cuò) for k=1:n p=。fprintf(39。 y=[2 3 4]。 y2=lagrange(x,y,x0)Error!x(1) is out of range! x=[1 2 3 4]。 y=[ ]。所以加入以下判斷以提高插值的準(zhǔn)確性n=length(x0)。Error! Please input again!\n39。,i)
職業(yè)教育相關(guān)推薦
鄂ICP備17016276號-1