西北農(nóng)林科技大學(xué)數(shù)值分析數(shù)值法實(shí)驗(yàn)報(bào)告(存儲(chǔ)版)
【正文】 on。原圖像39。 for j=1:n if j~=k p=p*(zx0(j))/(x0(k)x0(j))。而分段低次插值則精確度較高，擬合效果較好，而三次樣條插值具有良好的收斂性與穩(wěn)定性，與分段低次插值相比較光滑度更高，而且提供的信息也相對(duì)少一些。最值得欣賞的風(fēng)景，是自己奮斗的足跡。不要做金錢、權(quán)利的奴隸；應(yīng)學(xué)會(huì)做“金錢、權(quán)利”的主人。而在區(qū)間[0，64]上從圖32中可以看出藍(lán)色虛線條和綠色線條是幾乎重合的，而紅色線條在[30，70]之間有很大的振蕩，所在在區(qū)間[0，64]三次樣條插值更精確寫。 s=。三次樣條插值39。g39。o39。M5=。 p(j,2)=m(j+1)/(6*h(j))。for j=1:1:n a(j,j)=2。y=[0 1 2 3 4 5 6 7 8]。y4=polyval(S(4,:),x4X(4))。dxn= 。 y3=polyval(S(3,:),x3X(3))。Y=1./(25*X.^2+1)。for k=N2:1:1 M(k+1)=(U(k)C(k)*M(k+2))/B(k)。B(1)=B(1)H(1)/2。r39。 y=polyval(C,x)。 z=1./(1+25*xp.^2)。 [C,D]=newpoly(X,Y)。D=zeros(n,n) D(:,1)=Y39。)hold ongtext(39。plot(x,y)hold on for i=1:1:5y(i)=*(x(i))*(x(i))*(x(i)) *(x(i))*(x(i))*(x(i))*(x(i))*(x(i))*(x(i))*(x(i))endk=[0 1 10 11]x0=+*kfor i=1:1:4y0(i)=*(x(i))*(x(i))*(x(i)) *(x(i))*(x(i))*(x(i))*(x(i))*(x(i))*(x(i))*(x(i))endplot( x0,y0,39。 p(j,3)=(y(j)m(j)*(h(j)^2/6))/h(j)。end r(1)=0。y=[ ]。%39。)。已知n次牛頓插值多項(xiàng)式如下：Pn=f(x0)+f[x0,x1](xx0)+ f[x0,x1,x2](xx0) (xx1)+2）實(shí)驗(yàn)題目實(shí)驗(yàn)一：已知函數(shù)在下列各點(diǎn)的值為xi.f（xi）試用4次牛頓插值多項(xiàng)式P4（x）及三次樣條函數(shù)S（x）（自然邊界條件）對(duì)數(shù)據(jù)進(jìn)行插值。實(shí)驗(yàn)三：下列數(shù)據(jù)點(diǎn)的插值x01491625364964y012345678可以得到平方根函數(shù)的近似，在區(qū)間[0,64]上作圖。+ f[x0,x1，fx=[ ]。for i=2:n for j=i:n FF(j,i)=(FF(j,i1)FF(j1,i1))/(x(j)x(ji+1))。\n39。endfor j=1:1:n1 u(j)=1r(j)。 a(j,j+1)=r(j)。M1= 。o39。原圖像39。for k=(n1):1:1 C=conv(C,poly(X(k))) m=length(C)。 plot(x,y,X,Y,39。)得到插值函數(shù)和f（x）圖形：當(dāng)n=20時(shí)，我們?cè)贑ommand Window中輸入以下的命令：clear,clf,hold on。)。 D=diff(Y)./H。U(N1)=U(N
物理相關(guān)推薦
鄂ICP備17016276號(hào)-1