拉格朗日插值實(shí)驗(yàn)報(bào)告-wenkub
2022-08-17 19:03:10 本頁(yè)面
【正文】 ror!x(%d) is out of range!\n39。% x0和y0長(zhǎng)度不等時(shí)，報(bào)錯(cuò)elsefor i=1:m z=x(i)。m=length(x)。y(%d)=%f\nx1=%.3f y1=%.5f,x2=%.3f y2=%.5f\n\n39。 for s=l1:l if s~=k a=a*(zx0(s))/(x0(k)x0(s))。 break。 l=0。Error! Please input again!\n39。對(duì)應(yīng)于每一節(jié)點(diǎn)xk(k=0,1,…,n)，都能寫出一個(gè)滿足此條件的多項(xiàng)式，這樣寫出了n+1個(gè)多項(xiàng)式l0(x),l1(x),…,ln(x)，其中；由條件可得于是我們可以得出如下的拉格朗日n次插值多項(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)。在本實(shí)驗(yàn)中，采用拉格朗日插值法。分別通過(guò)分段線性插值、分段二次插值和全區(qū)間上拉格朗日插值計(jì)算f()，f()，f()的近似值。所以我們總是希望根據(jù)已有的數(shù)據(jù)點(diǎn)（或函數(shù)表）來(lái)構(gòu)造某個(gè)簡(jiǎn)單函數(shù)P(x)作為f(x)的近似值。實(shí)驗(yàn)名稱： 實(shí)驗(yàn)一 拉格朗日插值1 引言我們?cè)谏a(chǎn)生活中常常會(huì)遇到這樣的問(wèn)題：某個(gè)實(shí)際問(wèn)題中，函數(shù)f(x)在區(qū)間[a,b]上存在且連續(xù)，但卻很難找到其表達(dá)式，只能通過(guò)實(shí)驗(yàn)和觀測(cè)得到有限點(diǎn)上的函數(shù)表。插值法是解決此類問(wèn)題的一種比較古老的、但卻很常用的方法。已知函數(shù)表如下：xf(x)3 算法原理與流程圖（1）原理設(shè)函數(shù)y=在插值區(qū)間[a,b]上連續(xù)，且在n+1個(gè)不同的插值節(jié)點(diǎn)a≤x0,x1,…,xn≤b上分別取值y0,y1,…,yn。①分段低次插值當(dāng)給定了n+1個(gè)點(diǎn)x0x1…xn上的函數(shù)值y0,y1,…,yn后，若要計(jì)算x≠xi處函數(shù)值f(x)的近似值，可先選取兩個(gè)節(jié)點(diǎn)xi1與xi使x
職業(yè)教育相關(guān)推薦
鄂ICP備17016276號(hào)-1