6.3: Connecting Continuous - Time Models with DiscreteTime Models
- Page ID
- 7798
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)Continuous-time models and discrete-time models are different mathematical models with different mathematical properties. But it is still possible to develop a “similar” continuous-time model from a discrete-time model, and vice versa. Here we discuss how you can jump across the border back and forth between the two time treatments.
Assume you already have an autonomous first-order discrete-time model
\[x_{t}=F(x_{t-1}), \label{6.16} \]
and you want to develop a continuous-time model analogous to it. You set up the following “container” differential equation
\[\dfrac{dx}{dt} =G(x),\label{6.17} \]
and try to find out how \(F\) and \(G\) are related to each other.
Here, let me introduce a very simple yet useful analogy between continuous- and discrete-time models:
\[\dfrac{dx}{dt} \approx \dfrac{\Delta{x}}{\Delta{t}} \label{(6.18)} \]
This may look almost tautological. But the left hand side is a ratio between two infinitesimally small quantities, while the right hand side is a ratio between two quantities that are small yet have definite non-zero sizes. \(∆x\) is the difference between \(x(t + ∆t)\) and \(x(t)\), and \(∆t\) is the finite time interval between two consecutive discrete time points. Using this analogy, you can rewrite Equation \ref{6.17} as
\[\dfrac{\Delta{x}}{\Delta{t}} = \dfrac{x(t+\Delta{t}) -x(t)}{\Delta{t}} \approx G(x(t)), \label{6.19} \]
\[x(t+\Delta{t}) \approx x(t) +G(x(t))\Delta{t}.\label{(6.20)} \]
By comparing this with Eq.(6.16), we notice the following analogous relationship between \(F\) and \(G\):
\[F(x) \leftrightarrow \dfrac{F(x) -x}{\Delta{t}} \label{6.22} \]
For linear systems in particular, \(F(x)\) and \(G(x)\) are just the product of a coefficient matrix and a state vector. If \(F(x) = Ax\) and \(G(x) = Bx\), then the analogous relationships become
\[Ax \leftrightarrow x+Bx\Delta{t}, \label{6.23} \]
\[A \leftrightarrow I +B\Delta{t} \label{6.24} \]
or
\[B \leftrightarrow \dfrac{A-I}{\Delta{t}}. \label{6.25} \]
I should emphasize that these analogous relationships between discrete-time and continuous-time models do not mean they are mathematically equivalent. They simply mean that the models are constructed according to similar assumptions and thus they may have similar properties. In fact, analogous models often share many identical mathematical properties, yet there are certain fundamental differences between them. For example, one- or two-dimensional discrete-time iterative maps can show chaotic behaviors, but their continuous-time counterparts never show chaos. We will discuss this issue in more detail later.
Nonetheless, knowing these analogous relationships between discrete-time and continuous-time models is helpful in several ways. First, they can provide convenient pathways when you develop your own mathematical models. Some natural phenomena may be conceived more easily as discrete-time, stepwise processes, while others may be better conceived as continuous-time, smooth processes. You can start building your model in either way, and when needed, convert the model from discrete-time to continuous-time or vice versa. Second, Equation \ref{(6.20)} offers a simple method to numerically simulate continuous-time models. While this method is rather crude and prone to accumulating numerical errors, the meaning of the formula is quite straight forward, and the implementation of simulation is very easy, so we will use this method in the following section. Third, the relationship between the coefficient matrices given in Equation \ref{6.25} is very helpful for understanding mathematical differences of stability criteria between discrete-time and continuous-time models. This will be detailed in the next chapter.
Consider the dynamics of a system made of three parts, A, B, and C. Each takes a real-valued state whose range is [−1,1]. The system behaves according to the following state transitions:
- A adopts B’s current state as its next state.
- B adopts C’s current state as its next state.
- C adopts the average of the current states of A and B as its next state.
First, create a discrete-time model of this system, and then convert it into a continuous-time model using Equation \ref{6.25}.