該樓層疑似違規(guī)已被系統(tǒng)折疊?隱藏此樓查看此樓

求常微分方程的數(shù)值解 ode45方法的源代碼怎么看懂呢？四百多行 如何理解這些代碼的核心思想 以方便未來自己使用呢？求大神指點(diǎn)迷津 感激不盡

function varargout = ode45(ode,tspan,y0,options,varargin)

%ODE45 Solve non-stiff differential equations, medium order method.

% [TOUT,YOUT] = ODE45(ODEFUN,TSPAN,Y0) with TSPAN = [T0 TFINAL] integrates

% the system of differential equations y' = f(t,y) from time T0 to TFINAL

% with initial conditions Y0. ODEFUN is a function handle. For a scalar T

% and a vector Y, ODEFUN(T,Y) must return a column vector corresponding

% to f(t,y). Each row in the solution array YOUT corresponds to a time

% returned in the column vector TOUT. To obtain solutions at specific

% times T0,T1,...,TFINAL (all increasing or all decreasing), use TSPAN =

% [T0 T1 ... TFINAL].

%

% [TOUT,YOUT] = ODE45(ODEFUN,TSPAN,Y0,OPTIONS) solves as above with default

% integration properties replaced by values in OPTIONS, an argument created

% with the ODESET function. See ODESET for details. Commonly used options

% are scalar relative error tolerance 'RelTol' (1e-3 by default) and vector

% of absolute error tolerances 'AbsTol' (all components 1e-6 by default).

% If certain components of the solution must be non-negative, use

% ODESET to set the 'NonNegative' property to the indices of these

% components.

%

% ODE45 can solve problems M(t,y)*y' = f(t,y) with mass matrix M that is

% nonsingular. Use ODESET to set the 'Mass' property to a function handle

% MASS if MASS(T,Y) returns the value of the mass matrix. If the mass matrix

% is constant, the matrix can be used as the value of the 'Mass' option. If

% the mass matrix does not depend on the state variable Y and the function

% MASS is to be called with one input argument T, set 'MStateDependence' to

% 'none'. ODE15S and ODE23T can solve problems with singular mass matrices

總結(jié)

以上是生活随笔為你收集整理的matlab编程ode45,matlab 中ode45的源代码如何看懂的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。