appendix j

4
E1BAPP10 11/02/2010 12:18:0 Page 1 Appendix J: Derivation of the Time Domain Solution of State Equations J.1 Derivation Rather than using the Laplace transformation, we can solve the equations directly in the time domain using a method closely allied to the classical solution of differential equations. We will find that the final solution consists of two parts that are different from the forced and natural responses First, assume a homogeneous state equation of the form _ xðtÞ¼ AxðtÞ ðJ:1Þ Since we want to solve for x, we assume a series solution, just as we did in elementary scalar differential equations. Thus, xðtÞ¼ b 0 þ b 1 t þ b 2 t 2 þþ b k t k þ b kþ1 t kþ1 þ ðJ:2Þ Substituting Eq. (J.2) into (J.1) we get b 1 þ 2b 2 t þþ kb k t k1 þðk þ 1Þb kþ1 t k þ ¼ Aðb 0 þ b 1 t þ b 2 t 2 þþ b k t k þ b kþ1 t kþ1 þ Þ ðJ:3Þ Equating like coefficients yields b 1 ¼ Ab 0 ðJ:4aÞ b 2 ¼ 1 2 Ab 1 ¼ 1 2 A 2 b 0 ðJ:4bÞ . . . b k ¼ 1 k! A k b 0 ðJ:4cÞ b kþ1 ¼ 1 ðk þ 1Þ! A kþ1 b 0 . . . ðJ:4dÞ 1

Upload: abdullah1345

Post on 15-Sep-2015

227 views

Category:

Documents


0 download

DESCRIPTION

Appendix j

TRANSCRIPT

  • E1BAPP10 11/02/2010 12:18:0 Page 1

    Appendix J: Derivationof the Time Domain Solutionof State Equations

    J.1 Derivation

    Rather than using the Laplace transformation, we can solve the equations directly inthe time domain using a method closely allied to the classical solution of differentialequations. We will find that the final solution consists of two parts that are differentfrom the forced and natural responses

    First, assume a homogeneous state equation of the form

    _xt Axt J:1Since we want to solve for x, we assume a series solution, just as we did in elementaryscalar differential equations. Thus,

    xt b0 b1t b2t2 bktk bk1tk1 J:2Substituting Eq. (J.2) into (J.1) we get

    b1 2b2t kbktk1 k 1bk1tk Ab0 b1t b2t2 bktk bk1tk1 J:3

    Equating like coefficients yields

    b1 Ab0 J:4a

    b2 12Ab1 1

    2A2b0 J:4b

    ..

    .

    bk 1k!Akb0 J:4c

    bk1 1k 1!Ak1b0

    ..

    .

    J:4d

    1

  • E1BAPP10 11/02/2010 12:18:0 Page 2

    Substituting these values into Eq. (J.2) yields

    xt b0 Ab0t 12A2b0t

    2 1k!Akb0t

    k 1k 1!Ak1b0tk1

    IAt 12A2t2 1

    k!Aktk 1k 1!A

    k1tk1

    b0

    J:5

    But, from Eq. (J.2),

    x0 b0 J:61

    Therefore,

    xt IAt 12A2t2 1

    k!Aktk 1k 1!A

    k1tk1

    x0 J:7

    Let

    eAt IAt 12A2t2 1

    k!Aktk 1k 1!A

    k1tk1

    J:8

    where eAt is simply a notation for the matrix formed by the right-hand side of Eq.(J.8). We use this definition because the right-hand side of Eq. (J.8) resembles apower series expansion of eat, or

    eat 1 at 12a2t2 1

    k!aktk 1k 1! a

    k1tk1

    J:9

    Using Eq. (J.7), we have

    xt eAtx0 J:10We give a special name to eAt: it is called the state-transition matrix2, since it performsa transformation on x(0), taking x from the initial state, x(0), to the state x(t) at anytime, t. The symbol, F(t), is used to denote eAt. Thus,

    Ft eAt J:11and

    xt Ftx0 J:12There are some properties ofF(t) that we will use later when we solve for x(t)

    in the text. From Eq. (J.12),

    x0 F0x0 J:13Hence, the first property of F(t) is

    F0 I J:14

    1 In this development we consider the initial time, t0, to be 0. More generally, t0 6 0. After completing thisdevelopment, the interested reader should consult Appendix K on www.wiley.com/college/nise for themore general solution in terms of initial time t0 6 0.2 The state-transition matrix here is for the initial time t0 0. The derivation in Appendix Kon www.wiley.com/college/nise for t0 6 0 yields xt eAtt0xt0.

    2 Appendix J: Derivation of the Time Domain Solution of State Equations

  • E1BAPP10 11/02/2010 12:18:1 Page 3

    Bibliography 3

    where I is the identity matrix. Also, differentiating Eq. (J.12) and setting this equal toEq. (J.1) yields

    _xt _Ftx0 Axt J:15which, at t = 0, yields

    _F0x0 Ax0 J:16Thus, the second property of F(t) follows from Eq. (J.16):

    _F0 A J:17In summary, the solution to the homogeneous, or unforced, system is

    xt Ftx0 J:18where

    F0 I J:19and

    _F0 A J:20Let us now solve the forced, or nonhomogeneous, problem. Given the forced

    state equation

    _xtAxt But J:21rearrange and multiply both sides by eAt:

    eAt _xt Axt eAtBut J:22Realizing that the left-hand side is equal to the derivative of the product eAtx(t), weobtain

    d

    dteAtxt eAtBut J:23

    Integrating both sides yields

    eAtxt t0 eAtxt x0 Z t0eAtButdt J:24

    since eAt evaluated at t 0 is the identity matrix (from Eq. (J.8)). Solving for x(t) inEq. (J.24) we obtain

    xt eAtx0 Z t0eAttButdt

    Ftx0 Z t0Ft tButdt

    J:25

    where F(t) eAt by definition.

    BibliographyTimothy, L. K., and Bona, B. E. State Space Analysis: An Introduction. McGraw-Hill,New York, 1968.

  • E1BAPP10 11/02/2010 12:18:1 Page 4

    Copyright# 2011 by John Wiley & Sons, Inc.

    Founded in 1807, John Wiley & Sons, Inc. has been a valued source of knowledge and understanding for morethan 200 years, helping people around the world meet their needs and fulfill their aspirations. Our company isbuilt on a foundation of principles that include responsibility to the communities we serve and where we liveand work. In 2008, we launched a Corporate Citizenship Initiative, a global effort to address the environmental,social, economic, and ethical challenges we face in our business. Among the issues we are addressing are carbonimpact, paper specifications and procurement, ethical conduct within our business and among our vendors, andcommunity and charitable support. For more information, please visit our website: www.wiley.com/go/citizenship.

    No part of this publication may be reproduced, stored in a retrieval system or transmitted in any form or by anymeans, electronic, mechanical, photocopying recording, scanning or otherwise, except as permitted underSections 107 or 108 of the 1976 United States Copyright Act, without either the prior written permission of thePublisher or authorization through payment of the appropriate per-copy fee to the Copyright Clearance Center,222 Rosewood Drive, Danvers, MA 01923, (978) 750-8400, fax (978) 646-8600. Requests to the Publisher forpermission should be addressed to the Permissions Department, John Wiley & Sons, Inc., 111 River Street,Hoboken, NJ 07030-5774, (201) 748-6011, fax (201) 748-6008 or online at http://www.wiley.com/go/permissions.

    Evaluation copies are provided to qualified academics and professionals for review purposes only, for use intheir courses during the next academic year. These copies are licensed and may not be sold or transferred to athird party. Upon completion of the review period, please return the evaluation copy to Wiley. Returninstructions and a free of charge return shipping label are available at www.wiley.com/go/returnlabel. Outsideof the United States, please contact your local representative.

    ISBN 13 978-0470-54756-4

    Printed in the United States of America.

    10 9 8 7 6 5 4 3 2 1