歐拉法解常微分方程-免費閱讀
【正文】 y1+y2+y339。)。) figureplot(t,z,39。xlabel(39。ylabel(39。 z1(i1)=3*10^7*g(i1)^2。z1=zeros(size(t))。z0=0delta=。說明了顯性Eualr方法的收斂性是與步長的選擇是相關。迭代三次時，結果與方程描述內(nèi)容基本相符。 ylabel(39。y339。ty2變化圖39。)。)。 g(i)=g(i1)+g1(i1)*delta。 g1=zeros(size(t))。 g0 =0。歐拉方法是一類重要的數(shù)值解法。這類方法回避解y(x)的函數(shù)表達式，而是尋求它在一系列離散節(jié)點上的近似值，相鄰的兩個節(jié)點的間距稱作步長。z0=0 delta=。 z1=zeros(size(t))。 z1(i1)=3*10^7*g(i1)^2。 ylabel(39。 xlabel(39。) figure plot(t,z,39。)。y1+y2+y339。迭代1000次時，模擬結果已經(jīng)嚴重脫離事實，該迭代方法不收斂。這就對我們們選擇步長造成了困難，由于選擇的步長不合適有可能得出錯誤的結論。 time=。 f(1)=f0。 z(i)=z(i1)+z(i1)*delta。y139。t39。o39。title(39。)。ylabel(39。y339。ty2變化圖39。)。)。 g(i)=g(i1)+g1(i1)*delta。g1=zeros(size(t))。 g0 =0。 時間與個變量直接的變化關系如圖所示：從圖中能夠清晰看出，當h=，模擬結果與方程
法律信息相關推薦
鄂ICP備17016276號-1