資源簡(jiǎn)介
數(shù)值解與理論解對(duì)比可知,四階龍格-庫(kù)塔法的精度已經(jīng)很高,用它來(lái)解一般常微分方程已經(jīng)足夠了。 有程序運(yùn)行說(shuō)明
代碼片段和文件信息
function?y?=?DEEuler(f?haby0varvec)
%f:一階常微分方程的一般表達(dá)式的右端函數(shù)
%h:積分步長(zhǎng)
%a:自變量取值下限
%b:自變量取值上限
%y0:函數(shù)初值
%varvec:常微分方程的變量組
format?long;
N?=?(b-a)/h;
y?=?zeros(N+11);
x?=?a:h:b;
y(1)?=?y0;
for?i=2:N+1
????y(i)?=?y(i-1)+h*(-y(i-1)+x(i-1)+1);
end
format?short;
?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件????????319??2009-05-19?19:38??5\DEEuler.m
?????文件????????501??2009-05-19?19:54??5\DELGKT4_lungkuta.m
?????文件????????423??2009-05-19?19:44??5\DEModifEuler.m
?????文件??????61440??2010-11-04?23:07??5\Given?the?initial.doc
?????目錄??????????0??2010-11-04?23:07??5
-----------?---------??----------?-----??----
????????????????62683????????????????????5
- 上一篇:MATLAB小車跟蹤
- 下一篇:matlab光學(xué)分析程序
評(píng)論
共有 條評(píng)論