歐拉法解常微分方程-文庫(kù)吧
2025-07-09 00:27 本頁(yè)面
【正文】 性和隱形Eualar方法求解下列一階線性微分方程組的近似數(shù)值解： (2) 求解過(guò)程 Eular方法： 一階線性微分方程初值問(wèn)題 （1）方程離散化：差分和差商 （2） 通過(guò)初始值，依據(jù)遞推公式（2）逐步算出就為顯性的Eular方法。隱形Eular方法： （3） 公式（3）即為隱式Eular公式。（三）程序算法1. 利用顯式Eular法方求解 利用MATLAB進(jìn)行求解，編寫腳本文件如下：文件名： %顯性Eular方法 f0=1。 g0 =0。z0=0 delta=。 time=1。 t=0:delta:time。 f=zeros(size(t))。 g=zeros(size(t))。 z=zeros(size(t))。 f1=zeros(size(t))。 g1=zeros(size(t))。 z1=zeros(size(t))。 f(1)=f0。 g(1)=g0。 z(1)=z0。 for i=2:length(t) f1(i1) = *f(i1) + 10000*f(i1)*g(i1)。 f(i)=f(i1)+f1(i1)*delta。 g1(i1) = *f(i1) 10000*f(i1)*g(i1)3*10^7*g(i1)^2。 g(i)=g(i1)+g1(i1
法律信息相關(guān)推薦
鄂ICP備17016276號(hào)-1