CSE 245: Computer Aided Circuit Simulation and Verification

Similar documents
Transfer function and the Laplace transformation

REPETITION before the exam PART 2, Transform Methods. Laplace transforms: τ dτ. L1. Derive the formulas : L2. Find the Laplace transform F(s) if.

Boyce/DiPrima 9 th ed, Ch 2.1: Linear Equations; Method of Integrating Factors

Chapter 3: Fourier Representation of Signals and LTI Systems. Chih-Wei Liu

Physics 160 Lecture 3. R. Johnson April 6, 2015

Charging of capacitor through inductor and resistor

Midterm exam 2, April 7, 2009 (solutions)

Review Lecture 5. The source-free R-C/R-L circuit Step response of an RC/RL circuit. The time constant = RC The final capacitor voltage v( )

Lecture 2: Current in RC circuit D.K.Pandey

Boyce/DiPrima 9 th ed, Ch 7.8: Repeated Eigenvalues

Elementary Differential Equations and Boundary Value Problems

C From Faraday's Law, the induced voltage is, C The effect of electromagnetic induction in the coil itself is called selfinduction.

General Article Application of differential equation in L-R and C-R circuit analysis by classical method. Abstract

Lecture 1: Growth and decay of current in RL circuit. Growth of current in LR Circuit. D.K.Pandey

Wave Equation (2 Week)

fiziks Institute for NET/JRF, GATE, IIT JAM, JEST, TIFR and GRE in PHYSICAL SCIENCES MATEMATICAL PHYSICS SOLUTIONS are

LaPlace Transform in Circuit Analysis

Control System Engineering (EE301T) Assignment: 2

Voltage v(z) ~ E(z)D. We can actually get to this wave behavior by using circuit theory, w/o going into details of the EM fields!

where: u: input y: output x: state vector A, B, C, D are const matrices

Inverse Fourier Transform. Properties of Continuous time Fourier Transform. Review. Linearity. Reading Assignment Oppenheim Sec pp.289.

Chapter 5 The Laplace Transform. x(t) input y(t) output Dynamic System

whereby we can express the phase by any one of the formulas cos ( 3 whereby we can express the phase by any one of the formulas

Lecture 1: Numerical Integration The Trapezoidal and Simpson s Rule

H is equal to the surface current J S

On the Speed of Heat Wave. Mihály Makai

S.Y. B.Sc. (IT) : Sem. III. Applied Mathematics. Q.1 Attempt the following (any THREE) [15]

Poisson process Markov process

Revisiting what you have learned in Advanced Mathematical Analysis

10. The Discrete-Time Fourier Transform (DTFT)

EXERCISE - 01 CHECK YOUR GRASP

1. Inverse Matrix 4[(3 7) (02)] 1[(0 7) (3 2)] Recall that the inverse of A is equal to:

DEPARTMENT OF ELECTRICAL &ELECTRONICS ENGINEERING SIGNALS AND SYSTEMS. Assoc. Prof. Dr. Burak Kelleci. Spring 2018

Sinusoidal Response Notes

Logistic equation of Human population growth (generalization to the case of reactive environment).

Let s look again at the first order linear differential equation we are attempting to solve, in its standard form:

The Matrix Exponential

symmetric/hermitian matrices, and similarity transformations

AC STEADY-STATE ANALYSIS

UNIT #5 EXPONENTIAL AND LOGARITHMIC FUNCTIONS

Lagrangian for RLC circuits using analogy with the classical mechanics concepts

SECTION where P (cos θ, sin θ) and Q(cos θ, sin θ) are polynomials in cos θ and sin θ, provided Q is never equal to zero.

EE 434 Lecture 22. Bipolar Device Models

The Matrix Exponential

2.1. Differential Equations and Solutions #3, 4, 17, 20, 24, 35

Circuits and Systems I

MEM 355 Performance Enhancement of Dynamical Systems A First Control Problem - Cruise Control

CHAPTER. Linear Systems of Differential Equations. 6.1 Theory of Linear DE Systems. ! Nullcline Sketching. Equilibrium (unstable) at (0, 0)

EE 350 Signals and Systems Spring 2005 Sample Exam #2 - Solutions

An Indian Journal FULL PAPER. Trade Science Inc. A stage-structured model of a single-species with density-dependent and birth pulses ABSTRACT

A THREE COMPARTMENT MATHEMATICAL MODEL OF LIVER

Linear-Phase FIR Transfer Functions. Functions. Functions. Functions. Functions. Functions. Let

ES 250 Practice Final Exam

Spring 2006 Process Dynamics, Operations, and Control Lesson 2: Mathematics Review

Lecture 4: Laplace Transforms

The transition:transversion rate ratio vs. the T-ratio.

Why Laplace transforms?

EE202 Circuit Theory II , Spring. Dr. Yılmaz KALKAN & Dr. Atilla DÖNÜK

On the Derivatives of Bessel and Modified Bessel Functions with Respect to the Order and the Argument

AR(1) Process. The first-order autoregressive process, AR(1) is. where e t is WN(0, σ 2 )

MATH 308: Diff Eqs, BDP10 EXAMPLES [Belmonte, 2019] 1 Introduction 1.1 Basic Mathematical Models; Direction Fields

I) Title: Rational Expectations and Adaptive Learning. II) Contents: Introduction to Adaptive Learning

Ma/CS 6a Class 15: Flows and Bipartite Graphs

Chapter 12 Introduction To The Laplace Transform

SOLUTIONS. 1. Consider two continuous random variables X and Y with joint p.d.f. f ( x, y ) = = = 15. Stepanov Dalpiaz

FIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS I: Introduction and Linear Systems

Voltage v(z) ~ E(z)D. We can actually get to this wave behavior by using circuit theory, w/o going into details of the EM fields!

AN INTRODUCTION TO FOURIER ANALYSIS PROF. VEDAT TAVSANOĞLU

ANALYSIS IN THE FREQUENCY DOMAIN

7.4 QUANTUM MECHANICAL TREATMENT OF FLUCTUATIONS *

Lecture 26: Leapers and Creepers

4.3 Design of Sections for Flexure (Part II)

University of Kansas, Department of Economics Economics 911: Applied Macroeconomics. Problem Set 2: Multivariate Time Series Analysis

Quasi-Classical States of the Simple Harmonic Oscillator

Instructors Solution for Assignment 3 Chapter 3: Time Domain Analysis of LTIC Systems

Finite element discretization of Laplace and Poisson equations

Part I: Short Answer [50 points] For each of the following, give a short answer (2-3 sentences, or a formula). [5 points each]

On Ψ-Conditional Asymptotic Stability of First Order Non-Linear Matrix Lyapunov Systems

Bifurcation Theory. , a stationary point, depends on the value of α. At certain values

ECE 145A / 218 C, notes set 1: Transmission Line Properties and Analysis

Lecture 13 RC/RL Circuits, Time Dependent Op Amp Circuits

Laplace Transforms recap for ccts

Institute of Actuaries of India

CHAPTER CHAPTER14. Expectations: The Basic Tools. Prepared by: Fernando Quijano and Yvonn Quijano

Fourier Series and Parseval s Relation Çağatay Candan Dec. 22, 2013

Final Exam : Solutions

Chapter 7 Response of First-order RL and RC Circuits

EHBSIM: MATLAB-BASED NONLINEAR CIRCUIT SIMULATION PROGRAM (HARMONIC BALANCE AND NONLINEAR ENVELOPE METHODS)

CPSC 211 Data Structures & Implementations (c) Texas A&M University [ 259] B-Trees

Double Slits in Space and Time

Types of Transfer Functions. Types of Transfer Functions. Types of Transfer Functions. Ideal Filters. Ideal Filters

Module 1-2: LTI Systems. Prof. Ali M. Niknejad

State Observer Design

Complex Dynamic Models of Star and Delta Connected Multi-phase Asynchronous Motors

Types of Transfer Functions. Types of Transfer Functions. Ideal Filters. Ideal Filters. Ideal Filters

Laplace Transform. National Chiao Tung University Chun-Jen Tsai 10/19/2011

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is

a dt a dt a dt dt If 1, then the poles in the transfer function are complex conjugates. Let s look at f t H t f s / s. So, for a 2 nd order system:

Chapter 8 The Complete Response of RL and RC Circuits

2.3 Matrix Formulation

Transcription:

CSE 45: Compur Aidd Circui Simulaion and Vrificaion Fall 4, Sp 8 Lcur : Dynamic Linar Sysm

Oulin Tim Domain Analysis Sa Equaions RLC Nwork Analysis by Taylor Expansion Impuls Rspons in im domain Frquncy Domain Analysis From im domain o Frquncy domain Corrspondnc bwn im domain and frquncy domain Srial xpansion of (si-a) - Lcur.

Oulin (Con ) Modl Ordr Rducion Momns Passiviy, Sabiliy and Ralizabiliy Symbolic Analysis Y-Dla Transformaion BDD Analysis Lcur.3

Lcur.4 Sa of a sysm Th sa of a sysm is a s of daa, h valu of which a any im, oghr wih h inpu o h sysm a im, drmin uniquly h valu of any nwork variabl a im. W can xprss h sa in vcor form x Whr x i () is h sa variabls of h sysm ) (... ) ( ) ( x x x k

Sa Variabl How o Choos Sa Variabl? Th knowldg of h insananous valus of all branch currns and volags drmins his insananous sa Bu NOT ALL hs valus ar rquird in ordr o drmin h insananous sa, som can b drivd from ohrs. choos capacior volags and inducor currns as h sa variabls! Bu no all of hm ar chosn Lcur.5

Dgnra Nwork A nwork ha has a cu-s composd only of inducors and/or currn sourcs or a loop ha conains only of capaciors and/or volag sourcs is calld a dgnra nwork Exampl: Th following nwork is a dgnra nwork sinc C,C and C 5 form a dgnra capacior loop Lcur.6

Dgnra Nwork In a dgnrad nwork, no all h capaciors and inducors can b chosn as sa variabls sinc hr ar som rdundancy On h ohr hand, w choos all h capacior volags and inducors currns as sa variabl in a nondgnra nwork W will giv an xampl of how o choos sa variabl in h following scion Lcur.7

Ordr of Circui n b LC n C - n L n h ordr of circui, oal numbr of indpndn sa variabls b LC oal numbr of capaciors and inducors in h nwork n C numbr of dgnra loops (C-E loops) n L numbr of dgnra cu-ss (L-J cu-ss) n 4 3 In a nondgnra nwork, n quals o h oal numbr of nrgy sorag lmns Lcur.8

Sa Equaions Linar sysm of ordinary diffrnial quaions dx d Sa Ax( ) + Bu( ) Inpu Oupu y( ) Qx( ) + Du( ) Lcur.9

Sa Equaion for RLC Circuis Th sa quaion is of h form C v& E v L & YT il E R il Or M x&() Gx() + Pu() - + Pu v : volag in h runk, capacior volag i l : currn in h loop, inducor currn. Y and R ar h admianc marix and impdanc marix of cu-s and msh E covrs h co-r branchs in h cu-s E T covrs h r runks in h msh analysis Lcur.

Sa Equaions Mx&() Gx() + Pu() If w shif h marix M o h righ hand sid, w hav x& () M - Ax() + M - Bu() L A M - G and B M - P, w hav h sa quaion x&() Ax() + Bu() Toghr wih h oupu quaion y () Qx() + Du() ar calld h Sa Equaions of h linar sysm Lcur.

RLC Nwork Analysis A givn RLC nwork L 6 g 3 C Vs C C 5 g 4 Dgnra Nwork, Choos only volags of C and C 5, currn of L 6 as our sa variabl Lcur.

Tr Srucur Tak ino r as many capaciors as possibl and, as lss inducors as possibl Rsisors can b chosn as ihr r branchs or co-r branchs L 6 g 3 C C /L 6 g 3 C C 5 g4 Vs C C 5 g 4 Vs Lcur.3

Lcur.4 Linar Sa Equaion By a mixd cu-s and msh analysis, considr capacior cu-ss and inducor loops only. w can wri h linar sa quaion as follows M Gx() + Pu() x&() Cu-s KCL Loop KVL Cu-s KCL + + 6 5 L C C C C C C 6 i v v & & & 4 3 g g 6 i v v g 3 - + Vs

Oulin Tim Domain Analysis Sa Equaions RLC Nwork Analysis by Taylor Expansion Impuls Rspons in im domain Frquncy Domain Analysis From im domain o Frquncy domain Corrspondnc bwn im domain and frquncy domain Srial xpansion of (si-a) - Lcur.5

Solving RCL Equaion by Taylor Expansion() Gnral Circui Equaion X AX + BU Considr homognous form firs X AX X A X and A A I + +! A! +... + k A k! k +... Q: How o Compu A k? Lcur.6

Solving RCL Equaion by Taylor Expansion () Assum A has non-dgnra ignvalus,,..., k and corrsponding linarly indpndn ignvcors Χ, Χ,..., Χ k, hn A can b dcomposd as AΧΛΧ L whr L M and Λ M L O M L L k [ Χ Χ ] Χ,...,, Χ k Lcur.7

Solving RCL Equaion by Taylor Expansion (3) Wha s h implicaion hn? AΧΛΧ A Χ A ΧΛ Χ Λ Χ whr Λ To compu h ignvalus: M L L L L O L M M k d( I A) n + c n n +... + c ( p )( + p + p )...(...) ral ignvalu Conjugaiv Complx ignvalu Lcur.8

Lcur.9 In h prvious xampl () () / / / i v X i v l r l c A + Χ ΛΧ X X A 3 3 j j + + whr Χ 3 3 3 j j j + Χ 3 3 j j hnc Χ Χ + A L crl, w hav Solving RCL Equaion by Taylor Expansion (4)

Solving RCL Equaion by Taylor Expansion (5) Wha if marix A has dgnrad ignvalus? Jordan dcomposiion! A ΧJΧ J is in h Jordan Canonical form And sill A Χ J Χ Lcur.

Lcur. Jordan Dcomposiion J + + J L J + + J! L similarly

Oulin Tim Domain Analysis Sa Equaions RLC Nwork Analysis by Taylor Expansion Rspons in im domain Frquncy Domain Analysis From im domain o Frquncy domain Corrspondnc bwn im domain and frquncy domain Srial xpansion of (si-a) - Lcur.

Rspons in im domain W can solv h sa quaion and g h closd form xprssion Th oupu quaion can b xprssd as No: * dnos convoluion Lcur.3

Impuls Rspons Th Impuls Rspons of a sysm is dfind as h Zro Sa Rspons rsuling from an impuls xciaion Thus, in h oupu quaion, rplac u() by h impuls funcion δ(), and l x( ) w hav h() y() Q A B Lcur.4

Oulin Tim Domain Analysis Sa Equaions RLC Nwork Analysis by Taylor Expansion Impuls Rspons in im domain Frquncy Domain Analysis From im domain o Frquncy domain Corrspondnc bwn im domain and frquncy domain Srial xpansion of (si-a) - Lcur.5

Soluions in S domain By solving h sa quaion in s domain, w hav x(s) (si-a) - x( )+ (si-a) - Bu(s) y(s) Qx(s) +Du(s) Q(sI-A) - (x( ) + Bu(s)) +Du(s) Suppos h nwork has zro sa and h oupu vcor dpnds only on h sa vcor x, ha is, x( ) and D, w can driv h ransfr funcion of h nwork y( s) H(s) Q(sI-A) - B u( s) Lcur.6

Lcur.7 Ky Ky Transform Transform Propry: Propry: Bilaral Bilaral Laplac Laplac Transform: Transform: ) ( ) ( ) ( ) ( ) ( s Qx s y s Bu s Ax s sx + ) ( ) ( ) ( ) ( ) ( Qx y Bu Ax d dx + d dx sx(s) d x s x s ) ( ) ( x() Frquncy Domain Rprsnaion

Sysm Transfr Funcion sx( s) Ax( s) y( s) Qx( s) + Bu( s) Exprss y(s) as a funcion of u(s) y( s) Q( si A) Bu( s) Transfr Funcion: H (s) Lcur.8

Transfr Funcion Tim Domain Impuls Rspons Frquncy domain rprsnaion u(s) H(s) Linar sysm y(s) H(s) u(s) Tim domain rprsnaion u() y( ) h( τ ) u( τ ) dτ h() Linar sysm Th ransfr funcion H(s) is h Laplac Transform of h impuls rspons h() Lcur.9

Oulin Tim Domain Analysis Sa Equaions RLC Nwork Analysis by Taylor Expansion Impuls Rspons in im domain Frquncy Domain Analysis From im domain o Frquncy domain Corrspondnc bwn im domain and frquncy domain Srial xpansion of (si-a) - Lcur.3

Corrspondnc bwn im domain and frquncy domain W can driv h im domain soluions of h nwork from h s domain soluions by invrs Laplac Transformaion of h s domain soluions. Sa Equaions in S domain sx(s) x( ) Ax(s) +Bu(s) y(s) Qx(s) +Du(s) Invrs Laplac Transform Sa Equaions in im Domain x() L - [(si-a) - x( ) + (si-a) - Bu(s)] L - [(si-a) - ]x( ) + L - [(si-a) - ]B*u() y() L - [Q(sI-A) - (x( ) + Bu(s)) +Du(s)] Q L - [(si-a) - ] x( ) + {QL - [(si-a) - ]B +Dδ()}* u(s) Lcur.3

Corrspondnc bwn im domain and frquncy domain Soluion from im domain analysis Soluion by invrs Laplac ransform x() L - [(si-a) - x( ) + (si-a) - Bu(s)] L - [(si-a) - ]x( ) + L - [(si-a) - ]B*u() y() L - [Q(sI-A) - (x( ) + Bu(s)) +Du(s)] Q L - [(si-a) - ] x( ) + {QL - [(si-a) - ]B +Dδ()}* u(s) (si-a) - A muliplicaion of u(s) in s domain corrsponds o h convoluion in im domain Lcur.3

Oulin Tim Domain Analysis Sa Equaions RLC Nwork Analysis by Taylor Expansion Impuls Rspons in im domain Frquncy Domain Analysis From im domain o Frquncy domain Corrspondnc bwn im domain and frquncy domain Srial xpansion of (si-a) - Lcur.33

Srial xpansion of (si-a) - Whn s w can wri (si-a) - as (si-a) - -A - (I sa - ) -A - (I + sa - + s A - + + s k A -k + ) Thus, h ransfr funcion can b wrin as H(s) Q(sI-A) - B -QA - (I + sa - + s A - + + s k A -k + )B Whn s w can wri (si-a) - as (si-a) - s - (I s - A) - s - (I + s - A + s - A + + s -k A k + ) Th ransfr funcion can b wrin as H(s) Q(sI-A) - B s - (I + s - A + s - A + + s -k A k + )B Lcur.34

Marix Dcomposiion Assum A has non-dgnra ignvalus,,..., k and corrsponding linarly indpndn ignvcors Χ, Χ,..., Χ k, hn A can b dcomposd as A ΧΛΧ L whr L M and Λ M L O M L L k [ Χ Χ ] Χ,...,, Χ k Lcur.35

Lcur.36 Marix Dcomposiion Thn w can wri (si-a) - in h following form (si-a) - in s domain corrsponds o h xponnial funcion A in im domain, w can wri A as (si-a) - (SI XΛX - ) - X - (si Λ) - X X - n s s s.. X A X - n.. X