歐拉法解常微分方程(存儲(chǔ)版)
【正文】 所表示的顯示意義相吻合。迭代第二次時(shí)，結(jié)果與方程描述內(nèi)容基本相符。)。 ylabel(39。 title(39。o39。t39。 g1(i1) = *f(i1) 10000*f(i1)*g(i1)3*10^7*g(i1)^2。 f1=zeros(size(t))。（三）程序算法1. 利用顯式Eular法方求解 利用MATLAB進(jìn)行求解，編寫腳本文件如下：文件名： %顯性Eular方法 f0=1。求解從實(shí)際問題當(dāng)中歸結(jié)出來的微分方程主要靠數(shù)值解法。假定步長為定數(shù)。 time=1。 f(1)=f0。 z(i)=z(i1)+z(i1)*delta。y139。t39。o39。 title(39。)。時(shí)間與個(gè)變量直接的變化關(guān)系如圖所示：從上述圖形可以明顯看出，在迭代的不斷進(jìn)行時(shí)，各變量與時(shí)間的變化越來越大，且嚴(yán)重脫離了方程所描述的現(xiàn)實(shí)意義?！緦?shí)驗(yàn)小結(jié)】（收獲體會(huì)） 軟件使用 在寫MATLAB語言的時(shí)候要深刻理解題的意圖，整理好思緒再做題目，在我運(yùn)算的過程中，h取值取得越小、越細(xì)微，曲線逼近的越好。t=0:delta:time。g(1)=g0。 Fun=f+g+zend figureplot(t,f,39。)。)。)。ty3變化圖39。title(39。)。ylabel(39。title(39。o39。t39。 g1(i1) = *f(i1) 10000*f(i1)*g(i1)3*10^7*g(i1)^2。f1=zeros(size(t))。三、指導(dǎo)教師評(píng)語及成績(jī)：評(píng)
法律信息相關(guān)推薦
鄂ICP備17016276號(hào)-1