數(shù)值分析第五版計算實習(xí)題_第1頁
已閱讀1頁，還剩7頁未讀， 繼續(xù)免費閱讀

下載本文檔

版權(quán)說明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請進行舉報或認領(lǐng)

文檔簡介

1、數(shù)值分析計算實習(xí)題數(shù)值分析計算實習(xí)題第二章21程序：clearclcx1=[0.20.40.60.81.0]y1=[0.980.920.810.640.38]n=length(y1)c=y1(:)fj=2:n%求差商fi=n:1:jc(i)=(c(i)c(i1))(x1(i)x1(ij1))endendsymsxdfddf(1)=1d(1)=y1(1)fi=2:n%求牛頓差值多項式df(i)=df(i1)(xx1(i1))d(i)=c(

2、i)df(i)enddisp(4次牛頓插值多項式)P4=vpa(collect((sum(d)))5)%P4即為4次牛頓插值多項式并保留小數(shù)點后5位數(shù)pp=csape(x1y1variational)%調(diào)用三次樣條函數(shù)q=pp.coefsdisp(三次樣條函數(shù))fi=1:4S=q(i:)[(xx1(i))^3(xx1(i))^2(xx1(i))1]S=vpa(collect(S)5)endx2=0.2:0.08:1.08dot=[121

3、112]figureezplot(P4[0.21.08])holdony2=fnval(ppx2)x=x2(dot)y3=eval(P4)y4=fnval(ppx2(dot))plot(x2y2rx2(dot)y3bx2(dot)y4co)title(4次牛頓插值及三次樣條)結(jié)果如下：4次牛頓插值多項式P4=0.52083x^40.83333x^31.1042x^20.19167x0.98三次樣條函數(shù)Dw=zeros(1n)fi=1:n

4、Dw(i)=y1(i)polyval(q(i:)x1(i))%系數(shù)endp=DwqsymsxL8fi=1:nL8(i)=p(ni1)x^(i1)enddisp(8次拉格朗日插值)L8=vpa(collect((sum(L8)))5)xi=0:64yi=polyval(pxi)figureplot(xiyix1y1r)holdontitle(8次拉格朗日插值)結(jié)果如下：8次拉格朗日插值L8=3.2806e10x^86.7127e8x^75

溫馨提示

  • 1. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽，若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 眾賞文庫僅提供信息存儲空間，僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理，對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對任何下載內(nèi)容負責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容，請與我們聯(lián)系，我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論