15.11.2014 Views

MATLAB Mathematics - SERC - Index of

MATLAB Mathematics - SERC - Index of

MATLAB Mathematics - SERC - Index of

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

5 Differential Equations<br />

dde23 produces a solution that is continuous on [ t 0 , t f ]. You can use the<br />

function deval and the output <strong>of</strong> dde23 to evaluate the solution at specific<br />

points on the interval <strong>of</strong> integration.<br />

dde23 tracks discontinuities and integrates the differential equations with the<br />

explicit Runge-Kutta (2,3) pair and interpolant used by ode23. The<br />

Runge-Kutta formulas are implicit for step sizes longer than the delays. When<br />

the solution is smooth enough that steps this big are justified, the implicit<br />

formulas are evaluated by a predictor-corrector iteration.<br />

DDE Solver Basic Syntax<br />

The basic syntax <strong>of</strong> the DDE solver is<br />

sol = dde23(ddefun,lags,history,tspan,options)<br />

The input arguments are<br />

ddefun<br />

Handle to a function that evaluates the right side <strong>of</strong> the<br />

differential equations. The function must have the form<br />

dydt = ddefun(t,y,Z)<br />

where the scalar t is the independent variable, the column<br />

vector y is the dependent variable, and Z(:,j) is yt ( – τ j ) for<br />

τ j = lags(j). See “Function Handles” in the <strong>MATLAB</strong><br />

Programming documentation for more information.<br />

lags A vector <strong>of</strong> constant positive delays τ 1 , …,<br />

τ k .<br />

history Handle to a function <strong>of</strong> t that evaluates the solution y()<br />

t for<br />

t ≤ t 0 . The function must be <strong>of</strong> the form<br />

S = history(t)<br />

where S is a column vector. Alternatively, if yt () is constant,<br />

you can specify history as this constant vector.<br />

If the current call to dde23 continues a previous integration to<br />

t0, use the solution sol from that call as the history.<br />

tspan<br />

The interval <strong>of</strong> integration as a two-element vector [t0,tf]<br />

with t0 < tf.<br />

5-52

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!