西北農(nóng)林科技大學(xué)數(shù)值分析數(shù)值法實(shí)驗(yàn)報(bào)告(參考版)
2025-08-07 05:34本頁(yè)面
【正文】 壓力不是有人比你努力，而是那些比你牛幾倍的人依然比你努力。最值得欣賞的風(fēng)景，是自己奮斗的足跡。不要做金錢、權(quán)利的奴隸；應(yīng)學(xué)會(huì)做“金錢、權(quán)利”的主人。贈(zèng)語(yǔ)； 如果我們做與不做都會(huì)有人笑，如果做不好與做得好還會(huì)有人笑，那么我們索性就做得更好，來(lái)給人笑吧！ 現(xiàn)在你不玩命的學(xué)，以后命玩你。6）實(shí)驗(yàn)結(jié)果分析與小結(jié)通過(guò)此次實(shí)驗(yàn)，我對(duì)牛頓多項(xiàng)式插值，三次樣條插值，拉格朗日插值有了更進(jìn)一步的了解，知道了三次樣條的擬合程度在高次的情況下更高，在理論上和應(yīng)用上都有重要意義，在利用計(jì)算機(jī)編程軟件進(jìn)行高次插值的時(shí)候，我們可以多考慮利用三次樣條進(jìn)行插值。而分段低次插值則精確度較高，擬合效果較好，而三次樣條插值具有良好的收斂性與穩(wěn)定性，與分段低次插值相比較光滑度更高，而且提供的信息也相對(duì)少一些。而在區(qū)間[0，64]上從圖32中可以看出藍(lán)色虛線條和綠色線條是幾乎重合的，而紅色線條在[30，70]之間有很大的振蕩，所在在區(qū)間[0，64]三次樣條插值更精確寫。圖31為[0 1]的曲線，圖32為[0 64]區(qū)間上的曲線。 end y(i)=s。 for j=1:n if j~=k p=p*(zx0(j))/(x0(k)x0(j))。 s=。 m=length(x)。拉格朗日插值39。原圖像39。三次樣條插值39。hold on%axis([0 2 0 1])。:b39。hold on。g39。variational39。hold on。r39。o39。y=lagr1(x0,y0,x)。y0=[0 1 2 3 4 5 6 7 8]。M7=。M5=。M3=。M1=。 p(j,4)=(y(j+1)m(j+1)*(h(j)^2/6))/h(j)。 p(j,2)=m(j+1)/(6*h(j))。t=ap=zeros(n1,4)。 a(j,j+1)=r(j)。 u(n)=0。for j=1:1:n a(j,j)=2。endfor j=2:1:n1 d(j)=6*(f(j)f(j1))/(h(j1)+h(j))。endfor j=1:1:n1 u(j)=1r(j)。 n=9for j=1:1:n1 h(j)=x(j+1)x(j)。y=[0 1 2 3 4 5 6 7 8]。 三次樣條可以利用自然邊界條件。（1）用這9各點(diǎn)作8次多項(xiàng)式插值L8(x).（2）用三次樣條（自然邊界條件）程序求S（x）。.39。y4=polyval(S(4,:),x4X(4))。 y3=polyval(S(3,:),x3X(3))。y2=polyval(S(2,:),x2X(2))。y1=polyval(S(1,:),x1X(1))。dxn= 。Y=1./(25*X.^2+1)。.39。y4=polyval(S(4,:),x4X(4))。 y3=polyval(S(3,:),x3X(3))。y2=polyval(S(2,:),x2X(2))。y1=polyval(S(1,:),x1X(1))。dxn= 。Y=1./(25*X.^2+1)。 S(k+1,4)=Y(k+1)。 S(k+1,2)=M(k+1)/2。M(N+1)=3*(dxnD(N))/H(N)M(N)/2。for k=N2:1:1 M(k+1)=(U(k)C(k)*M(k+2))/B(k)。 U(k)=U(k)temp*U(k1)。 for k=2:N1 temp=A(k1)/B(k1)。B(N1)=B(N1)H(N)/2。B(
物理相關(guān)推薦
鄂ICP備17016276號(hào)-1