經(jīng)典Runge-kutta方法數(shù)值分析實驗報告_第1頁
經(jīng)典Runge-kutta方法數(shù)值分析實驗報告_第2頁
經(jīng)典Runge-kutta方法數(shù)值分析實驗報告_第3頁
經(jīng)典Runge-kutta方法數(shù)值分析實驗報告_第4頁
經(jīng)典Runge-kutta方法數(shù)值分析實驗報告_第5頁
全文預(yù)覽已結(jié)束

下載本文檔

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

文檔簡介

1、數(shù)值分析實驗報告姓名學(xué)號日期實驗項目經(jīng)典Runge-kutta方法指導(dǎo)教師一、 上機實驗的問題和要求（需求分析）：實驗要求掌握求解一階常微分方程初值問題經(jīng)典Runge-kutta方法編程運算實驗內(nèi)容：用經(jīng)典Runge-kutta方法在步長選為1/10時計算一階常微分方程初值問題（該問題的精確解為），將誤差（精確解與數(shù)值解的差）列表， 并畫出精確解與數(shù)值解的函數(shù)圖象.2、 程序設(shè)計的基本思想，原理和算法描述：主要利用相應(yīng)公式建立經(jīng)典Runge-kutta方法的RungeK函數(shù)；3、 主要程序代碼或命令：1建立m文件，構(gòu)建函數(shù)RungeKfunction y=RungeK(f,h,a,b,y0,v

2、ar)format long;N=(b-a)/h;y=zeros(N+1,1);y(1)=y0;g(1)=y0;x=a:h:b;for i=2:N+1 k(1)=subs(f,var,x(i-1) y(i-1); k(2)=subs(f,var,x(i-1)+h/2 y(i-1)+h/2*k(1); k(3)=subs(f,var,x(i-1)+h/2 y(i-1)+h/2*k(2); k(4)=subs(f,var,x(i-1)+h y(i-1)+h*k(3); y(i)=y(i-1)+h/6*(k(1)+2*k(2)+2*k(3)+k(4);endplot(x,y,'-rs')在命令窗口輸入如下語句 syms u v>> z=-v+u+1;>> m=RungeK(z,0.1,0,1,1,u v)2 再畫出方程準(zhǔn)確解的圖像，程序語句如下在命令窗口輸入hold onx=0:0.1:1;y=zeros(11,1);for i=1:11 y(i)=exp(-x(i)+x(i);endplot(x,y,'*b')3 計算誤差在命令窗口輸入m-y四、調(diào)試和運行程序過程中產(chǎn)生的問題及采取的措施：開始畫出的圖像，忘記了用不同的顏色區(qū)分開來；經(jīng)修改后，圖像客觀性增強；5、 運行輸出結(jié)果及分析：結(jié)果

溫馨提示

  • 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. 人人文庫網(wǎng)僅提供信息存儲空間，僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理，對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容，請與我們聯(lián)系，我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論