next up previous
Next: 7 Time-domain analysis of Up: guide Previous: 5 Representation of systems

Subsections


6 Linear systems with delays

6.1 Linear time-varying systems

To simulate the linear system with delays

$\displaystyle \dot{x}$(t) = A0(tx(t) + $\displaystyle \sum_{i=1}^{m}$Ai(tx(t - $\displaystyle \tau_{i}^{}$(t)) + $\displaystyle \int\limits_{-\tau (t)}^{0}$G(t, sx(t + s) ds + B(t) , (6.1)

(here A0(t), Ai(t) (i = 1,..., m) are n×n matrices with piece-wise continuous elements, G(t, s) is n×n matrix with piece-wise continuous elements, B(t) is an n-dimensional vector with piece-wise continuous elements) one can use the function dde45lin.

The conditional representation of (6.1) is

$\displaystyle \dot{x}$ = A0(tx + $\displaystyle \sum_{i=1}^{m}$Ai(ty(- $\displaystyle \tau_{i}^{}$(t)) + $\displaystyle \int\limits_{-\tau (t)}^{0}$G(t, sy(s) ds + B(t) , (6.2)

so, to simulate the system it is necessary to define the following finite number of parameters:

- matrices A0(t), Ai(t) (i = 1,..., m), B(t),

- a matrix-function G(t, s),

- an initial time-moment t0,

- an initial point x0,

- an initial prehistory y0( . ) = {y0(s), - $ \tau$ $ \leq$ s < 0}.

6.2 Linear time-invariant systems

To simulate the linear time-invariant system with delay

$\displaystyle \dot{x}$(t) = A0 x(t) + A$\scriptstyle \tau$ x(t - $\displaystyle \tau$) + $\displaystyle \int\limits_{-\tau}^{0}$G(sx(t + s) ds (6.3)

with respect to the basic initial functions one can use the function test.

Here A0, A$\scriptstyle \tau$ are constant n×n matrices, G(s) is n×n matrix with piece-wise continuous on [- $ \tau$, 0] elements.

The conditional representation of (6.3) is

$\displaystyle \dot{x}$ = A0 x + A$\scriptstyle \tau$ y(- $\displaystyle \tau$) + $\displaystyle \int\limits_{-\tau}^{0}$G(sy(s) ds (6.4)

and to simulate the system it is necessary to define the following finite number of parameters:

- matrices A0, A$\scriptstyle \tau$,

- a matrix-function G(s),

- an initial point x0,

- an initial prehistory y0( . ) = {y0(s), - $ \tau$ $ \leq$ s < 0}.

Remark 6.1. For time-invariant system (6.3) we consider the initial time-moment t0 equals to zero, i.e. t0 = 0. $ \Box$

Remark 6.2. The matrix-function G(s) can be described using M-file. $ \Box$

6.3 Linear time-invariant systems with exponential distributed delays

In many cases linear time-invariant systems with delays have the form

$\displaystyle \dot{x}$(t) = A0 x(t) + A$\scriptstyle \tau$ x(t - $\displaystyle \tau$) + D0$\displaystyle \int\limits_{-\tau}^{0}$eD1sD2 x(t + s) ds (6.5)

where A0, A$\scriptstyle \tau$, D0, D1 and D2 are constant n×n matrices.

The conditional representation of (6.5) is

$\displaystyle \dot{x}$ = A0 x + A$\scriptstyle \tau$ y(- $\displaystyle \tau$) + D0$\displaystyle \int\limits_{-\tau}^{0}$eD1sD2 y(s) ds . (6.6)

One can simulate such systems with respect to the basic initial functions using teste.


next up previous
Next: 7 Time-domain analysis of Up: guide Previous: 5 Representation of systems