資源簡介
數值解與理論解對比可知,四階龍格-庫塔法的精度已經很高,用它來解一般常微分方程已經足夠了。 有程序運行說明
代碼片段和文件信息
function?y?=?DEEuler(f?haby0varvec)
%f:一階常微分方程的一般表達式的右端函數
%h:積分步長
%a:自變量取值下限
%b:自變量取值上限
%y0:函數初值
%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;
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????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光學分析程序
評論
共有 條評論