15.11.2014 Views

MATLAB Mathematics - SERC - Index of

MATLAB Mathematics - SERC - Index of

MATLAB Mathematics - SERC - Index of

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Initial Value Problems for DDEs<br />

Using a History to Specify the Solution <strong>of</strong> Interest<br />

In an initial value problem, we seek the solution on an interval [ t 0 , t f ]. with<br />

t 0 < t f . The DDE shows that y′ () t depends on values <strong>of</strong> the solution at times<br />

prior to t. In particular, y′ ( t 0 ) depends on yt ( 0 – τ 1 ),…,<br />

yt ( 0 – τ k ). Because <strong>of</strong><br />

this, a solution on [ t 0 , t f ] depends on its values for t ≤ t 0 , i.e., its history St ().<br />

Propagation <strong>of</strong> Discontinuities<br />

Generally, the solution yt () <strong>of</strong> an IVP for a system <strong>of</strong> DDEs has a jump in its<br />

first derivative at the initial point t 0 because the first derivative <strong>of</strong> the history<br />

function does not satisfy the DDE there.<br />

– +<br />

S′ ( t 0 ) ≠ y′ ( t 0 ) = ft ( 0 , yt ( 0 ), St ( 0 – τ 1 ), …,<br />

St ( 0 – τ k ))<br />

A discontinuity in any derivative propagates into the future at spacings <strong>of</strong><br />

τ 1 , τ 2 , …,<br />

τ k .<br />

For reliable and efficient integration <strong>of</strong> DDEs, a solver must track<br />

discontinuities in low order derivatives and deal with them. For DDEs with<br />

constant lags, the solution gets smoother as the integration progresses, so after<br />

a while the solver can stop tracking a discontinuity. See “Discontinuities” on<br />

page 5-57 for more information.<br />

DDE Solver<br />

This section describes:<br />

• The DDE solver, dde23<br />

• DDE solver basic syntax<br />

The DDE Solver<br />

The function dde23 solves initial value problems for delay differential<br />

equations (DDEs) with constant delays. It integrates a system <strong>of</strong> first-order<br />

differential equations<br />

y′ () t = ftyt ( , ()yt , ( – τ 1 ),…,<br />

yt ( – τ k ))<br />

on the interval [ t 0 , t f ], with t 0 < t f and given history yt () = St () for t ≤ t 0 .<br />

5-51

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

Saved successfully!

Ooh no, something went wrong!