ece_350_lecture12_fall08
22 Pages

ece_350_lecture12_fall08

Course Number: TH 350, Fall 2008

College/University: CSU Northridge

Word Count: 1344

Rating:

Document Preview

Lecture 12 Lecture Overview End-of-Year Calendar Review Solving State Variable Equations Homogeneous Solution Eigenvalues & Eigenvectors Particular Solution Complete Solution ECE 350 D. van Alphen 1 Repeating: End-of-Year Calendar Tuesday Dec. 2 State Variables Dec. 9 Quiz 3 Thursday Dec. 4 Homework Review ECE 350 D. van Alphen 2 Review Input-output relationship for LTI systems: h(t) x(t) H()...

Unformatted Document Excerpt
Coursehero >> California >> CSU Northridge >> TH 350

Course Hero has millions of student submitted documents similar to the one
below including study guides, practice problems, reference materials, practice exams, textbook help and tutor support.

Course Hero has millions of student submitted documents similar to the one below including study guides, practice problems, reference materials, practice exams, textbook help and tutor support.

12 Lecture Lecture Overview End-of-Year Calendar Review Solving State Variable Equations Homogeneous Solution Eigenvalues & Eigenvectors Particular Solution Complete Solution ECE 350 D. van Alphen 1 Repeating: End-of-Year Calendar Tuesday Dec. 2 State Variables Dec. 9 Quiz 3 Thursday Dec. 4 Homework Review ECE 350 D. van Alphen 2 Review Input-output relationship for LTI systems: h(t) x(t) H() y(t) = x(t) * h(t) Y() = X() H() Y(s) = X(s) H(s) State Variables 2 dt = 1 2 Parsevals Theorem: E x = | x( t) | | X() | d 2 For circuits: assign a state to the voltage across a capacitor, and a state to the current through an inductor. For an Nth order differential equation: assign the states y(t) = q1(t), dy/dt = q2(t), d(N-1y/dtN-1 = qN(t) ECE 350 D. van Alphen 3 Summarizing State Variable Equations for Circuit Problems State variables, continued State equation, vector/matrix form: dq/dt = A q(t) + B x(t) Output equation: y(t) = C q(t) + D x(t) (1) Remaining Problem: How do we solve equation (1) for q(t)? Approach: first find the natural (or homogeneous) solution; i.e., the solution when the input x(t) = 0. Solve dq/dt = A q(t) (2) ECE 350 D. van Alphen 4 Solving the Homogeneous Equation Aside: differential equation review - For the scalar differential equation: & q( t) = a q( t), the solution is: q(t) = q(0) eat I.C. : q(0) Claim: the solution of homogeneous equation (2) is of the same form: q(t) = eAt q(0) where Problem: not a closed-form. k 2 3 At t Ak = 1 + t A + t A 2 + t A3 + K e 2 3! k =0 k! (compare to Maclaurin Series for ex) ECE 350 D. van Alphen 5 Linear Algebra Aside: Eigenvalues & Eigenvectors Multiplying an arbitrary vector x by a conformable square matrix A yields product Ax that is typically a scaled (longer or shorter) and rotated version of the original vector x: Ax 2-dimensional sketch of arbitrary matrix multiplication: x However, for any square matrix A, there will be some vectors, zi, (called the eigenvectors of A) that do not get rotated under the matrix multiplication. They do get scaled, and we call the scalars the eigenvalues, i. ECE 350 D. van Alphen 6 Eigenvalues & Eigenvectors, continued In other words, the eigenvectors of square matrix A are the vectors that do not change direction* when multiplied by A. The eigenvalues are the scalars by which the vectors are shrunk or stretched: 2-dimensional sketch of multiplying an eigenvector of A by A Az = z z Defining equation for eigenvalues () and eigenvectors (z) of square matrix A: Az = z * Here, we do not consider a rotation of 180 a change in direction. ECE 350 D. van Alphen 7 Eigenvalues & Eigenvectors, continued: From the previous equation, we obtain: (A I) z = 0 and det(A - I) = 0 (use to find eigenvalues, s) (called the characteristic eq. of A) 0 I = 0 (use to find eigenvectors, z, once s are known) Example: A = 4 5 , 2 3 1 0 I= , 0 1 4 5 0 A I = 0 = 2 3 ECE 350 _____________ To find A-I: subtract from the diagonal elements of A. 8 D. van Alphen Eigenvalues & Eigenvectors, continued Hence det(A I) = 0 4 2 5 3 = (4 )(3 ) + 10 = 0 2 12 + 10 = 0 2 2 = 0 ( 2) ( + 1) = 0 Now use = 2, -1 Finding the eigenvalues by solving the characteristic equation of A : (A - ) z = 0 to find the eigenvector for each eigenvalue: 5 z1 0 4 2 ( A I)z = z = 0 3 2 2 2 For = 2: 5 2 5 z1 0 = 2z1 5z 2 = 0 z1 = z 2 2 5 z 0 2 2 ECE 350 D. van Alphen 9 Eigenvalues & Eigenvectors, continued 2.5 z= So for = 2: z1 = (5/2) z2 eigenvector is 1 Any multiple of the eigenvector above would also work (eigenvectors are not unique); To standardize our answers, we typically give the unit-length version of the eigenvector ( divide z above by its length): length( z) = (5 / 2)2 + 12 = 7.25 2.69 Unit-length eigenvector: 2.5 / 2.69 .929 z= = .371 1/ 2.69 ECE 350 D. van Alphen 10 Eigenvalues & Eigenvectors, continued Now use : (A - ) z = 0 to find the eigenvector z for = -1: 4 + 1 5 z1 0 ( A I)z = z = 0 3 + 1 2 2 5 5 z1 0 1 = 5z1 5z 2 = 0 z1 = z 2 z = 2 2 z 0 2 1 For standardization, we normalize z above, dividing it by its length: length( z) = (1)2 + 12 = 2 1.414 1/ 1.414 .707 z= = .707 1/ 1.414 ECE 350 (normalized) 11 D. van Alphen MATLAB: Eigenvalues & Eigenvectors To obtain just the eigenvalues of A, type: >> eig(A) or >> e_values = eig(A) To obtain the both eigenvectors and the corresponding eigenvalues, type: >> [V D] = eig(A) Note: V is a matrix whose columns are the eigenvectors of A; D is a diagonal matrix, with eigenvalues on the diagonal. The order of the eigenvalues in D matches the order of the eigenvectors in V. ECE 350 D. van Alphen 12 MATLAB: Eigenvalues & Eigenvectors MATLAB Code for eigenvalues & eigenvectors of A = 2 3 >> A = [4 -5; 2 -3]; >> eig(A) ans = 2 -1 >> [V D] = eig(A) V= 0.9285 0.7071 0.3714 0.7071 D= 20 0 -1 ECE 350 D. van Alphen 13 4 5 Finding Integer Powers of Matrix A Cayley-Hamilton Theorem: Any N-by-N matrix A satisfies its own characteristic equation: det(A I) = 0. Example: 5 4 A= 1 2 5 4 = (5 )(2 ) 4 = 0 2 Characteristic equation: det 1 2 7 + 6 = 0 So, per Cayley-Hamilton: ________________ A2 = _____________ ECE 350 D. van Alphen (3) 14 Finding Integer Powers of Matrix A Repeating equation (3): A2 5 4 A= = 7A 6I, for 1 2 Multiplying (3) by A: A3 = 7A2 6A = 7(7A 6I) 6A = 43A 42I Continuing, we could find any integer power of A as a linear combination of A and I. Claim 1: for any N-by-N matrix A, we can find integer powers of A as a linear combination of Am, m = 0, 1, , N-1. Claim 2: for any N-by-N matrix A, we can evaluate eAt as a timevarying linear combination of Am, m = 0, 1, , N-1. At = N1 ( t) Am e m m=0 ECE 350 D. van Alphen (4) 15 Finding eAt e At = m ( t ) Am m=0 N1 Approach: 1. Find the eigenvalues of A MATLAB function: eig(A) 2. Find the time-varying coefficients, i(t), by writing N equations, replacing A (in eq. (4)) by the (distinct) eigenvalues. t = N1 ( t) m e m m =0 (4) 3. Return to equation (4), and plug in the coefficients, m(t), found in step 2. ECE 350 D. van Alphen 16 eAt Example: 3 / 4 0 A= 1/ 2 1/ 2 Step 1 - Using MATLAB to find the eigenvalues: > A = [3/4 0; 1/2 1/2]; >> eig(A) ans = 0.5000 0.7500 Step 2 - System of equations (from (4)): ( = 1/ 2) : t / 2 = 1 ( t ) m = ( t ) + (1/ 2) ( t ) e m 0 1 m=0 ( = 3 / 4) : e ECE 350 3t / 4 = m ( t) m = 0 ( t) + (3 / 4)1( t) m=0 17 1 D. van Alphen eAt Example: 3 / 4 0 A= 1/ 2 1/ 2 Repeating the system of equations, and subtracting the top equation from the bottom equation: e t / 2 = 0 ( t) + (1/ 2)1( t) 1/ 4 1( t) = e3t / 4 e t / 2 1( t) = 4e3t / 4 4e t / 2 e3t / 4 = 0 ( t) + (3 / 4)1( t) Plugging the answer for 1(t) into the top equation: e t / 2 = 0 ( t) + 2e3t / 4 2e t / 2 ( ) 0 ( t) = 3e t / 2 + 2e3t / 4 Step 3 - Now plug in the values for 0(t), 1(t) into the original equation (4): At = 1 ( t ) Am = ( t ) I + ( t ) A1 e m 0 1 m =0 ECE 350 D. van Alphen 18 eAt Repeating Example: 3 / 4 0 A= 1/ 2 1/ 2 0 t / 2 + 2e3 t / 4 3e t / 2 + 3t / 4 2e At = ( t ) I + ( t ) A1 = 3e e 0 1 0 + ( 4e At = 3e e t / 2 + 3t / 4 2e 3t / 4 4e t / 2 3 / 4 0 ) 1/ 2 1/ 2 0 3e3t / 4 3e t / 2 0 + t / 2 + 2e3 t / 4 2e3 t / 4 2e t / 2 3e 0 t / 2 + 4e 3 t / 4 e D. van Alphen 0 3 t / 4 2e t / 2 2e At = e 5e 3 t / 4 t / 2 + 2e3 t / 4 2e ECE 350 19 Homogeneous Solution for State Variable Problem Recall from p. 5: q(t) = eAt q(0) 3 / 4 0 1 , and q(0) = Suppose, for example, that: A = 1/ 2 1/ 2 2 Then q(t) = eAt q(0) 5e 3 t / 4 = t / 2 + 2e3 t / 4 2e 5e3t / 4 = =5 3t / 4 10e 1 0 t / 2 + 4e3 t / 4 2 e e3 t / 4 3t / 4 2 e ECE 350 D. van Alphen 20 The Particular Solution for State Variable Equations Recall the defining state variable equations: State equation: Output equation: dq/dt = A q(t) + B x(t) y(t) = C q(t) + D x(t) (1) Notation: The complete solution for the states q(t) in equation (1) is the sum of the homogeneous solution (found by letting x(t) = 0) and the particular solution. So far, we can find the homogeneous solution: qh(t) = eAt q(0) We could show that the particular solution is: qp ( t) = e At B x() d 0 ECE 350 D. van Alphen 21 t The Complete Solution for State Variable Equations The complete solution for the state vector is the sum of the homogeneous solution and the particular solution: q( t) = qh ( t) + qp ( t) = e q(0) + e A( t ) B x( ) d At 0 t Another use (in the state variable model) for eAt: Impulse response: h(t) = C eAt B + D (t), t 0 (0 else0 Stability: An LTI system is stable if and only if all of the eigenvalues of matrix A are in the left half of the complex plane; i.e., iff all the eigenvalues have negative real parts. ECE 350 D. van Alphen 22

Find millions of documents on Course Hero - Study Guides, Lecture Notes, Reference Materials, Practice Exams and more. Course Hero has millions of course specific materials providing students with the best way to expand their education.

Below is a small sample set of documents:

CSU Northridge - TH - 350
ECE 350Fall, 2008 1/5Homework Set 7 - Solutions Solve problems 1, 2, and 3 by inspection, knowing that the signals only will pass through the filter if they are in the passband. 1a. y(t) = 200 + 4 sin(200t) 2a. y(t) = 4 sin(250t) 3a. y(t) = 10 +
CSU Northridge - TH - 355
The Theme of Natural Order in "The Tempest" Lawrence E. Bowling College English, Vol. 12, No. 4. (Jan., 1951), pp. 203-209.Stable URL: http:/links.jstor.org/sici?sici=0010-0994%28195101%2912%3A4%3C203%3ATTONOI%3E2.0.CO%3B2-I College English is curre
CSU Northridge - TH - 600
Esther Dabagyan CSUN SED 600 02/20/2007 Review #1 Video I video-taped my 6th period for a 25 minute period during which I reviewed the warm-up with the class and started them working on a laboratory where they made silly putty. During the 25 minute
CSU Northridge - TH - 600
Review#4page1KevinD.McMahon StudentID#:78513 SED600 March27,2007Review#4:Conference"Itisnotenoughtohaveknowledge,onemustalsoapplyit.Itisnotenoughtohave wishes,onemustalsoaccomplish." ~JohannWolfgangvonGoetheConference OnJanuary19th,2007Iatte
CSU Northridge - TH - 600
Afrodita Fuentes SED 600 05/09/07Conference - Review #4 (Instructions from syllabus: Examination of information in reference to my own research and how will my research is influenced by this information in this conference?)I attended the Secondary
CSU Northridge - TH - 630
Manufacturing Systems Engineering and ManagementMSE 630 Electronic Materials [3 units] Course Description 2000-2002: Prerequisites: consent of instructor Electrical behaviors of materials; conductors; semiconductors and insulators; electronic struct
CSU Northridge - TH - 630
Physics 630 Statistical Physics Spring 2005 Lecture Room: 1100 (Science I, 1st floor) Meeting Time: Tuesday and Thursday 6:30 7:45Logistics pm InstructorNicholas Kioussis Office: Science I 1123 Office Hours: Tu: 5:30 6:25 pm; F: 9:00 10:45 pm;
CSU Northridge - TH - 695
Thursday, February 23, 2006 Dr. Flubber Paper Airplane WorkshopOngoing Workshop 9:00 am to 11:30 am Boardwalk BallroomSchedule of Events9:00am Boardwalk Ballroom 9:30 am Limited to 30 teams of 1-3 students (9th 12th) 10:15 am Limited to 10 te
CSU Northridge - W S - 420
420 Ainsworth Lab 9 - More ANCOVA 1. Calculate the adjusted group means for the data set in Lab #8.Presoak Time a1 = 5 min X Y 1 2 5 6 5 6 4 6 3 5 2 3 a2 = 5 min X Y 1 3 8 9 7 9 4 5 5 7 3 5 X 5 9 5 6 2 9 a3 = 5 min Y 8 11 9 9 6 12 a4 = 5 min X Y 4 6
CSU Northridge - HSCI - 237
The Dimensions of Health Education LiteratureAn Evaluation of Published Health Education Research (20002005)Examples of Information SourcesMagazinesProfessional JournalsNewspapersInternet Based InformationTypes of Information Sources Pr
CSU Northridge - HSCI - 303
HSci303 Highfield1SETGOALSYourGOALS w YOU areyourmostimportantasset. w YOU createyourgrade:A,B,C,DorF. w YOUhavealreadyearned5of100pointsinthisclasswithLASSI (Gofor5morewiththe selfcareplan.) WHYIMPORTANT? w Energytosucceed! w Motivationtocometo
CSU Northridge - HSCI - 303
Nuts&Bolts(Week#3) M.Highfield,RN,Phd Nursing,CSUN Today TurninTreasureHunt.Whatdidyoulearn? Quicklibraryliteraturefollowup Healthlinenews:Whatdidyoulearn? Text:Whatdidyoulearn? LASSI Realityshock EvaluatingInformation Peerreviewed Online Typ
CSU Northridge - HSCI - 345
A Public Health ChallengePoor Nutrition and ObesityThe Public Health IssueIt is estimated that 55% of the adult population in the United States are either obese or overweight.* Obesity substantially increases an individuals risk of suffering from
CSU Northridge - HSCI - 345
Ethical Considerations in Scarce Resource Allocation Paper Instructions Students will form small groups of 5-6 people to conduct a public health ethics task force simulation for the influenza vaccine allocation case study. Each group will take on the
CSU Northridge - HSCI - 345
Protecting Childrens HealthPreventing Childrens Exposure to Endocrine Disruptors by Lawrence H. Keith, Ph.D.Version 1.2Introduction3 333There are many chemicals that are present in our environment that have unknown impacts on children One
CSU Northridge - HSCI - 431
CALIFORNIA STATE UNIVERSITY, NORTHRIDGE HEALTH SCIENCES DEPARTMENT HSCI 431 HEALTH BEHAVIOR FALL 2008 SAMPLE MID-TERM EXAM ITEMSINSTURCTIONS: CHOOSE THE BEST ANSWER OUT OF THE POSSIBLE CHOICES. 1. The ecological approach to understanding factors in
CSU Northridge - HSCI - 431
1.2.3.4.5.6.7.8.California State University, Northridge Health Science 431 Fall 2008 Final Exam Study Items Which of the following is reason for conducting health promotion interventions in schools? a. Health and education are intertw
CSU Northridge - HSCI - 431
COMMUNITY MENTAL HEALTHWhat is Mental Health?The emotional and social wellbeing of an individual, including ones psychological resources for dealing with the day-to-day problems of life.Good Mental Health is the Ability to:Function under advers
CSU Northridge - HSCI - 431
Developmental Disabilities: A Public Health School Issue1What are Developmental Disabilities (DD)?Severe and chronic interference of a persons functionality which is attributable to a mental or physical impairment or a combination of mental
CSU Sacramento - ACCY - 1
California State University, SacramentoACCOUNTING 141: ACCOUNTING INFORMATION SYSTEMS DEVELOPMENTFall 2004 - Section 1 Instructor: Dr. Yan Xiong Office: TAH-2082 Office Phone: 278-7143 E-mail: xiongy@csus.edu WWW address: http:/www.csus.edu/indiv/
CSU Sacramento - ACCY - 1
CSU Sacramento - ACCY - 111
ACCY 111 Fall 2004Accounts Receivable Chapter Five skip pages 326 to 331 for now I. A Few Comments A. We will cover the notes payable section when we to Chapter Six at the end of the semester. B. Accounts receivable is a good example of two major
CSU Sacramento - ACCY - 111
ACCY 111 Fall 2004 I Components of CashChapter 7 Part I Cash and Cash EquivalentsA. For HP this is the disclosure in their Form 10-K: Cash and Cash Equivalents HP classifies investments as cash equivalents if the original maturity of an invest
CSU Sacramento - ACCY - 111
Chapter Three: SOLUTIONS TO BRIEF EXERCISESBRIEF EXERCISE 3-1 May 1 Cash . Common Stock . Equipment .. Accounts Payable . Rent Expense . Cash .. Accounts Receivable. Service Revenue . 3,000 3,000 1,100 1,100 400 400 500 50031321BRIEF EXERCI
CSU Sacramento - ACCY - 112
Intermediate Accounting, 10th Edition Kieso, Weygandt, and WarfieldChapter 13: Current Liabilities and ContingenciesPrepared by Krishnan Ranganathan, Angelo State University, San Angelo, TexasLiabilities in Generalw Liabilitiesare: probable
CSU Sacramento - ACCY - 113
ACCY 113 Practice Exam #1 - Solutions posted on Friday 1. The revenue recognition principle provides that revenue is recognized when a. it is realized. b. it is realizable. c. it is realized or realizable and it is earned. d. none of these. When good
CSU Sacramento - ACCY - 113
ACCY 113 Spring 2005 I. Class NotesChapter 20 Accounting Pension Plans - 11. Check class web page for practice midterm after 4 PM on Wednesday. Solution posted as well. 2. Extra office hours: see exam web page note that I will be in Tahoe 1026
CSU Sacramento - ACCY - 121
ACCY 121 Possible Midterm Questions from Horngren Chapter 2 Cost behavior refers to: Costs that are initially recorded as assets and expensed when sold are called: For a manufacturing company, direct labor cost
CSU Sacramento - ACCY - 121
ACCY 121 A Sample of Final Exam Questions from Chapters 10, 11,15, and 21 (Ignore income taxes in this problem.) Which of the following would be used in the calculation of the internal rate of return of an investment in new machinery to replace old
CSU Sacramento - ACCY - 161
ACCY 161Budgetary Accounting and Journal Entry ExercisePost the following transactions to the t-accounts below. Then complete the "Appropriations,.Ledger." Use only the accounts listed below. What is the "balance" in the budgetary accounts at the
CSU Sacramento - ACCY - 161
ACCY 161 Summer 2005Midterm #1 Comments1. The exam is Thursday, June 9th and will start of 8:15 AM and will end at 10:00 AM you should have more than enough time. When you are finished with the exam you are done for the day. No lecture after th