資源簡介
采用歐拉法、改進歐拉法、龍格庫塔法(經典RK法)求解常微分方程初值問題的自編MATLAB代碼。所有函數均獨立成文件便于移植。代碼的使用結合一個具體題目說明,題目來源為浙江大學數值計算方法作業。

代碼片段和文件信息
f_h=@ODE_fun;
x0=0;
y0=1;
xn=1.00001;
h=0.1;
[~Y1]=MethodClassicalRK(f_hx0y0xnh);
[~Y2]=MethodImprovedEuler(f_hx0y0xnh);
[XY3]=MethodEuler(f_hx0y0xnh);
Y=(1+X.*X).^(1/3);
plot(XY1‘kx‘XY2‘go‘XY3‘rx‘XY‘b‘);
legend(‘經典R-K‘‘改進歐拉方法‘‘歐拉方法‘‘準確解‘);
xlswrite(‘result.xlsx‘[X;Y1;Y2;Y3;Y]1‘B1‘);
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????342??2018-12-15?14:36??常微分方程求解\main.m
?????文件?????????481??2016-06-05?13:20??常微分方程求解\MethodClassicalRK.m
?????文件?????????300??2018-12-15?14:35??常微分方程求解\MethodEuler.m
?????文件?????????371??2018-12-15?14:35??常微分方程求解\MethodImprovedEuler.m
?????文件?????????116??2018-12-15?14:32??常微分方程求解\ODE_fun.m
?????文件???????56901??2018-12-15?14:38??常微分方程求解\示例題目.docx
評論
共有 條評論