歐拉法解常微分方程(參考版)
2024-08-04 00:27本頁面
【正文】 ty1+y2+y3變化圖39。)。ylabel(39。t39。) figureplot(t,Fun)。title(39。y339。)。xlabel(39。o39。ty2變化圖39。)。ylabel(39。t39。)。) figureplot(t,g,39。title(39。y139。)。xlabel(39。o39。 z(i)=z(i1)+z(i1)*delta。 g(i)=g(i1)+g1(i1)*delta。 f(i)=f(i1)+f1(i1)*delta。z(1)=z0。 f(1)=f0。g1=zeros(size(t))。z=zeros(size(t))。 f=zeros(size(t))。 time=。 g0 =0。 實(shí)驗(yàn)感想 在這次上機(jī)實(shí)驗(yàn)中，我掌握了解決常微分方程的基本方法，同時(shí)學(xué)會(huì)使用計(jì)算機(jī)軟件對(duì)兩種不同方法得到的結(jié)果進(jìn)行判斷，對(duì)我們以后對(duì)數(shù)據(jù)進(jìn)行分析很有幫助。歐拉法的缺點(diǎn) 簡單地取切線的端點(diǎn)作為下一步的起點(diǎn)進(jìn)行計(jì)算，當(dāng)步數(shù)增多時(shí)，誤差會(huì)因積累而越來越大。這就對(duì)我們們選擇步長造成了困難，由于選擇的步長不合適有可能得出錯(cuò)誤的結(jié)論。 時(shí)間與個(gè)變量直接的變化關(guān)系如圖所示：從圖中能夠清晰看出，當(dāng)h=，模擬結(jié)果與方程所表示的顯示意義相吻合。跌三次迭代結(jié)果，并未產(chǎn)生誤差。=，模擬結(jié)果如下： 迭代第一次， 與A中結(jié)果相同。迭代1000次時(shí)，模擬結(jié)果已經(jīng)嚴(yán)重脫離事實(shí)，該迭代方法不收斂。迭代第二次時(shí)，結(jié)果與方程描述內(nèi)容基本相符。)【實(shí)驗(yàn)結(jié)論】 A步長h=。 title(39。y1+y2+y339。)。 xlabel(39。t
法律信息相關(guān)推薦
鄂ICP備17016276號(hào)-1