Page 188 - Numerical Analysis Using MATLAB and Excel
P. 188

Solution of Single State Equations



                               x · 1        0   1   0   0           x 1        0
                           ·
                           x =  x · 2  ,       A =  0  0  1  0  ,       x =  x 2  ,       b =  0  ,       and u =  ut()
                               x · 3        0   0   0   1           x 3        0
                               x · 4         a –  0  a –  1  a –  2  a –  3  x 4  1




               5.8 Solution of Single State Equations

               Let us consider the state equations
                                                    ·
                                                    x =  αx + βu                                      (5.111)
                                                    y =  k x +  k u
                                                               2
                                                          1
               where  , , αβ k 1 , and k 2  are scalar constants, and the initial condition, if non−zero, is denoted as

                                                      x =  xt()                                       (5.112)
                                                              0
                                                       0
               We will now prove that the solution of the first state equation in (5.111) is

                                                   (
                                                 α t –  t )  αt  t  – ατ
                                         xt() =  e    0  x +  e  ∫  e  βu τ() τ                       (5.113)
                                                                          d
                                                        0
                                                               t 0
               Proof:
               First, we must show that (5.113) satisfies the initial condition of (5.112). This is done by substitu-
               tion of t =  t 0   in (5.113). Then,

                                                                t
                                                  (
                                                 α t –  t )  0  αt  0  – ατ
                                                   0
                                                                           d
                                        xt () =  e      x +  e  ∫  e  βu τ() τ                        (5.114)
                                           0
                                                         0
                                                                t 0
               The first term in the right side of (5.114) reduces to x 0  since
                                                 α t –  t )  0  0
                                                  (
                                                   0
                                                e      x =  e x =   x 0                               (5.115)
                                                               0
                                                        0
               The second term of (5.114) is zero since the upper and lower limits of integration are the same.
               Therefore, (5.114) reduces to xt () =  x 0  and thus the initial condition is satisfied.
                                               0
               Next, we must prove that (5.113) satisfies also the first equation in (5.111). To prove this, we dif-
               ferentiate (5.113) with respect to   and we obtain
                                                t
                                                               d ⎧
                                                     (
                                         ·
                                                                -
                                                                 ⎨
                                         x t() =  ----- e (  d  α t –  t )  0  x )  +  ---- e αt ∫  t e – ατ βu τ() τ ⎬  d  ⎫
                                                dt         0  dt ⎩    t 0          ⎭
               Numerical Analysis Using MATLAB® and Excel®, Third Edition                              5−27

               Copyright © Orchard Publications
   183   184   185   186   187   188   189   190   191   192   193