Model Reduction for Dynamical Systems Lecture 6
|
|
- Cecil Davis
- 5 years ago
- Views:
Transcription
1 Oo-von-Guericke Universiä Magdeburg Faculy of Mahemaics Summer erm 07 Model Reducion for Dynamical Sysems ecure 6 v eer enner and ihong Feng Max lanck Insiue for Dynamics of Complex echnical Sysems Compuaional Mehods in Sysems and Conrol heory Magdeburg Germany benner@mpi-magdeburg.mpg.de feng@mpi-magdeburg.mpg.de hp://
2 Ouline MOR: alanced runcaion
3 Overlook alanced runcaion: firs balancing hen runcae. Given a I sysem: / x y u x d dx For convenience of discussion we denoe he sysem as a block form: alancing runcae reduced model he unimporan par is runcaed
4 Eigenspaces of and make he wo measures pracically compuable ill now i seems we could do he runcaion by finding subspace spanned by he eigenvecors corresponding o he small eigenvalues of or. However i could happen ha saes which are difficul o reach produce he maximal energy of observaion; saes which produce he smalles energy of observaion are neverheless he easies o reach! For such sysem we do no know which saes o runcae!
5 Example: Consider he following I sysem / x y u x d dx he wo Gramians are: heir eigenvalues and eigenvecors are: λ ξ λ ξ λ ξ λ ξ Eigenspaces of and make he wo measures pracically compuable
6 Eigenspaces of and make he wo measures pracically compuable ξ λ ξ λ ξ he angle beween ξ and is very small. his means if S is he subspace spanned by hen he easily observable saes ξ ξ ξ migh also be in S. ξ αξ α α x α >> - I ells us if we runcae he saes which are difficul o reach he saes locae in S we risk runcaing he saes which are easy o observe produce he maximal energy of observaion because hey migh also be in S.
7 Eigenspaces of and make he wo measures pracically compuable However if and have he same eigenvalues and eigenvecors hen he problems is solved. he saes in he subspace spanned by he eigenvecors of corresponding o he small eigenvalues always in he subspace spanned by he eigenvecors of corresponding o he small eigenvalues because he eigenvalues are he same and eigenvecors are he same herefore he subspaces are he same. Can we achieve his? Yes. e can achieve i by balancing.
8 MOR: alanced runcaion alancing asic idea of balancing ransformaion: x x Use sae space ransformaion o ge anoher realizaion of he same sysem so ha he ransformed Gramians are diagonal marices. Definiion of alancing ransformaion: Finding a nonsingular marix such ha and. Definiion of alanced sysem: he reachable observable and sable I sysem is balanced if is wo Gramians are equal i is principal-axis balanced if. diag σ σ n.
9 MOR: alanced runcaion alancing asic idea of balancing ransformaion: x x Use sae space ransformaion o ge anoher realizaion of he same sysem so ha he ransformed Gramians are equal and are diagonal marices. I.e. How o consruc? Recall ha. Since we have which means. should be he inverse of he marix Y of eigenvecors of.
10 MOR: alanced runcaion Check :? How o make alancing? If UU hen UU U UU U Here we mus have he relaion U. / / If furher U hen U UU How o compue? / Subsiue U he lef hand side / / U I ino we ge / U U / / U UU U U if /. U I. /. ook a he righ hand side we ge / / U U i. e. U U. herefore is he inverse of he marix of eigenvecors of U U. Furunaely he ranspose of U U is a s.p.d. marix. herefore he inverse of he marix of he marix iself.so ha we do no have o compue he inverse. eigenvecors is exacly
11 MOR: alanced runcaion alancing he above analysis clearly shows ha: Exisence of balancing ransformaion: dx / d y x Given a reachable observable and sable I sysem and he corresponding Gramians and a principal axis balancing ransformaion is given as follows: x u K U and UK / / Here UU is he Cholesky facorizaion of. U U K is he eigen-decomposiion of U U. Symmeric posiive semi-definie marix has real non-negive eigenvalues and orhogonal eigenvecors. Here he Eigenvecors in K are aken as orhonormal K
12 MOR: alanced runcaion alancing ha is he corresponding balanced sysem? pply he sae space ansformaion: x x o he original realizaion: dx / d y x x u x x dx / d y x x u
13 alancing : MOR: alanced runcaion alancing Given dx / d x y x u Compue. Compue UU U U K K he eigenvalues are ordered from he larges o he smalles dx / d x u / K U dx / d x u y x / UK y x
14 balanced sysem: dx / d y x x u he uni vecors e are he eigenvecors of : e σ e ssume ha he elemens on he herefore e e easily observable saes. r i i i i i n. diagonal of is already ordered as : σ σ σ. span he subspace conaining he firs r easily conrollable and runcae he difficul - o - observe and difficul - o - conrol saes means : x e x e x e x e x x x 0 0. MOR: alanced runcaion n n I.e. x x 0 0 :. Replace xr x x wih x dx / d x u z : x xr y x r r runcae r in he balanced sysem : z z u d 0 0 u y z 0 n
15 MOR: alanced runcaion runcae z y u u z z d noher realizaion of he same sysem is: is a non-minimal realizaion of a sysem. ˆ z y u z dz he reduced-order model ROM herefore we have he following simple seps for runcaion:
16 alancing: MOR: alanced runcaion runcae dx / d y x x u K U / / UK dx / d y x x u runcae: Small par and Separaed according o he separaion of. dz / d z u yˆ z ROM!
17 MOR: alanced runcaion alancing: dx / d y x x u / K U / UK dx / d y x x u Does i make sense if we do model reducion on he balanced sysem raher han he original sysem? Yes. s a sae ransformaion balancing does no change he ransfer funcion and he HSVs. he balanced sysem is only a differen realizaion of he sysem.
18 Observe: / x y u x d dx / x y u x d dx ˆ / z y u z d dz Y Y Y Yx x x :. produc marix he of eigenvecors he are : in columns he Here.. hen and as as separae If x Y Y x Yx x x x x x Y Y Y Y Y Y Y Y Y Y Y Y Y Y Y Y MOR: alanced runcaion
19 MOR: alanced runcaion herefore he wo ROMs are he same: dz / d yˆ z z u dx / d Y x y Y x u Conclusion: balanced runcaion is erov-galerkin projecion as below: e x Y x dx / d y x u x erov - Galerkin Y x x dx / d y Y x using Y x u herefore balanced runcaion is equivalen o: finding he invarian subspace rangey of and remaining only he par Y which corresponds o he larges HSVs square roo of he eigenvalues of.
20 MOR: alanced runcaion lgorihm alancing:. Compue. Given dx / d x u y x. Compue UU 3. K U / U U K / UK K 4. alancing and separaing according o he separaion of : runcae: 5. Form he reduced model: dxˆ / d yˆ xˆ xˆ u
21 MOR: alanced runcaion compuaional deails Recall: In M use command: lyap ' lyap '
22 MOR: alanced runcaion Numerical issues he balancing marix is: / K U UU. Compuaion of U may cause numerical insabiliy because U is usually near singular. U is usually near singular because he marix has numerically low-rank i.e. near singular. is near singular because in may cases is eigenvalues decay rapidly o zero some eigenvalues are very close o zero e.g.. and behaves similarly as. λ 0 0 i However in algorihm we need o compue: / / K U UK Can we avoid compuing U?
23 MOR: alanced runcaion p p p p compuaional deails If using Choelsky facorizaion of boh Observe Use SVD insead of eigen- decomposiion U V Comparing wih defined in lgorihm we immediaely ge / U p o avoid compuing he inverse of U V V p U p we have: / V
24 MOR: alanced runcaion Numerical issues he balanced sysem which is balanced by and is: V / / U V V V / / / / / / U U U Furhermore if Define:. V V / V / U V we have: hen V V V U U U
25 MOR: alanced runcaion Numerical issues lgorihm SR mehod Geing he reduced model wihou compuing U :. Do Cholesky facorizaion of he wo Gramians: are lower riangular marices.. Do Singular value decomposiion SVD of marix i.e. here are wo orhonormal marices U V U U I V V I such ha V U V U U. V / 3. e / V V U. 4. e ˆ ˆ ˆ V V. 5. he reduced model is dxˆ / d x ˆ ˆ u ˆ yˆ ˆ xˆ
26 MOR: alanced runcaion Numerical issues lgorihm someimes canno coninue eiher because he Cholesky facorizaion of canno be done. his is because ha in some cases and include oo small eigenvalues like: λ 0 0 which is considered by he algorihm as a singular marix herefore Cholesky facorizaion canno be coninued. aper [enner 05] provides an algorihm compuing he numerically full rank facors n n n n of and which are in he forms ˆ ˆ R R n ˆ << n he full rank facors numerically saisfy:. [enner 05]. enner E.S. uiana-ori Model reducion based on specral projecion mehods. In:. enner V.. Mehrmann D.C. Sorensen eds. "Dimension Reduion of arge-scale Sysems" vol. 45 of ecure Noes in Compuaional Science and Engineering pp Springer-Verlag erlin/heidelberg 005. lgorihm 4 in he paper
27 MOR: alanced runcaion Numerical issues lgorihm 3 Geing he reduced model using full-rank facors [enner 05]:. Compue full-rank facors of he Gramians:. Compue SVD. V V U U V U / V 3. e. / U V 5. he reduced model is ˆ ˆ ˆ ˆ ˆ ˆ / ˆ x y u x d dx 4. e. ˆ ˆ ˆ V V. ˆ ˆ ˆ n n R R n n n n <<
28 MOR: alanced runcaion Error bound heorem [enner 05] If he original I sysem is sable hen he reduced model obained by lgorihm lgorihm lgorihm 3 saisfies: he reduced model is balanced minimal and sable. I s Gramians are equal o he same diagonal marix. he absolue error bound proof in [noulas 05] Chaper 7 holds. H s Hˆ s σ H n k r k
Then. 1 The eigenvalues of A are inside R = n i=1 R i. 2 Union of any k circles not intersecting the other (n k)
Ger sgorin Circle Chaper 9 Approimaing Eigenvalues Per-Olof Persson persson@berkeley.edu Deparmen of Mahemaics Universiy of California, Berkeley Mah 128B Numerical Analysis (Ger sgorin Circle) Le A be
More information10. State Space Methods
. Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he
More informationt is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...
Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger
More informationChapter 2. First Order Scalar Equations
Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.
More informationMath 334 Fall 2011 Homework 11 Solutions
Dec. 2, 2 Mah 334 Fall 2 Homework Soluions Basic Problem. Transform he following iniial value problem ino an iniial value problem for a sysem: u + p()u + q() u g(), u() u, u () v. () Soluion. Le v u. Then
More information2.160 System Identification, Estimation, and Learning. Lecture Notes No. 8. March 6, 2006
2.160 Sysem Idenificaion, Esimaion, and Learning Lecure Noes No. 8 March 6, 2006 4.9 Eended Kalman Filer In many pracical problems, he process dynamics are nonlinear. w Process Dynamics v y u Model (Linearized)
More informationSolutions of Sample Problems for Third In-Class Exam Math 246, Spring 2011, Professor David Levermore
Soluions of Sample Problems for Third In-Class Exam Mah 6, Spring, Professor David Levermore Compue he Laplace ransform of f e from is definiion Soluion The definiion of he Laplace ransform gives L[f]s
More informationLinear Response Theory: The connection between QFT and experiments
Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and
More informationd 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3
and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or
More informationCHAPTER 2 Signals And Spectra
CHAPER Signals And Specra Properies of Signals and Noise In communicaion sysems he received waveform is usually caegorized ino he desired par conaining he informaion, and he undesired par. he desired par
More informationHamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t
M ah 5 2 7 Fall 2 0 0 9 L ecure 1 0 O c. 7, 2 0 0 9 Hamilon- J acobi Equaion: Explici Formulas In his lecure we ry o apply he mehod of characerisics o he Hamilon-Jacobi equaion: u + H D u, x = 0 in R n
More informationOptimality Conditions for Unconstrained Problems
62 CHAPTER 6 Opimaliy Condiions for Unconsrained Problems 1 Unconsrained Opimizaion 11 Exisence Consider he problem of minimizing he funcion f : R n R where f is coninuous on all of R n : P min f(x) x
More informationInstitute for Mathematical Methods in Economics. University of Technology Vienna. Singapore, May Manfred Deistler
MULTIVARIATE TIME SERIES ANALYSIS AND FORECASTING Manfred Deisler E O S Economerics and Sysems Theory Insiue for Mahemaical Mehods in Economics Universiy of Technology Vienna Singapore, May 2004 Inroducion
More informationHomework sheet Exercises done during the lecture of March 12, 2014
EXERCISE SESSION 2A FOR THE COURSE GÉOMÉTRIE EUCLIDIENNE, NON EUCLIDIENNE ET PROJECTIVE MATTEO TOMMASINI Homework shee 3-4 - Exercises done during he lecure of March 2, 204 Exercise 2 Is i rue ha he parameerized
More information(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4)
Phase Plane Analysis of Linear Sysems Adaped from Applied Nonlinear Conrol by Sloine and Li The general form of a linear second-order sysem is a c b d From and b bc d a Differeniaion of and hen subsiuion
More informationChapter 6. Systems of First Order Linear Differential Equations
Chaper 6 Sysems of Firs Order Linear Differenial Equaions We will only discuss firs order sysems However higher order sysems may be made ino firs order sysems by a rick shown below We will have a sligh
More informationCourse Notes for EE227C (Spring 2018): Convex Optimization and Approximation
Course Noes for EE7C Spring 018: Convex Opimizaion and Approximaion Insrucor: Moriz Hard Email: hard+ee7c@berkeley.edu Graduae Insrucor: Max Simchowiz Email: msimchow+ee7c@berkeley.edu Ocober 15, 018 3
More informationSystem of Linear Differential Equations
Sysem of Linear Differenial Equaions In "Ordinary Differenial Equaions" we've learned how o solve a differenial equaion for a variable, such as: y'k5$e K2$x =0 solve DE yx = K 5 2 ek2 x C_C1 2$y''C7$y
More informationGuest Lectures for Dr. MacFarlane s EE3350 Part Deux
Gues Lecures for Dr. MacFarlane s EE3350 Par Deux Michael Plane Mon., 08-30-2010 Wrie name in corner. Poin ou his is a review, so I will go faser. Remind hem o go lisen o online lecure abou geing an A
More informationMatrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality
Marix Versions of Some Refinemens of he Arihmeic-Geomeric Mean Inequaliy Bao Qi Feng and Andrew Tonge Absrac. We esablish marix versions of refinemens due o Alzer ], Carwrigh and Field 4], and Mercer 5]
More informationOrdinary Differential Equations
Ordinary Differenial Equaions 5. Examples of linear differenial equaions and heir applicaions We consider some examples of sysems of linear differenial equaions wih consan coefficiens y = a y +... + a
More informationSPECTRAL EVOLUTION OF A ONE PARAMETER EXTENSION OF A REAL SYMMETRIC TOEPLITZ MATRIX* William F. Trench. SIAM J. Matrix Anal. Appl. 11 (1990),
SPECTRAL EVOLUTION OF A ONE PARAMETER EXTENSION OF A REAL SYMMETRIC TOEPLITZ MATRIX* William F Trench SIAM J Marix Anal Appl 11 (1990), 601-611 Absrac Le T n = ( i j ) n i,j=1 (n 3) be a real symmeric
More informationRecursive Least-Squares Fixed-Interval Smoother Using Covariance Information based on Innovation Approach in Linear Continuous Stochastic Systems
8 Froniers in Signal Processing, Vol. 1, No. 1, July 217 hps://dx.doi.org/1.2266/fsp.217.112 Recursive Leas-Squares Fixed-Inerval Smooher Using Covariance Informaion based on Innovaion Approach in Linear
More informationTechnical Report Doc ID: TR March-2013 (Last revision: 23-February-2016) On formulating quadratic functions in optimization models.
Technical Repor Doc ID: TR--203 06-March-203 (Las revision: 23-Februar-206) On formulaing quadraic funcions in opimizaion models. Auhor: Erling D. Andersen Convex quadraic consrains quie frequenl appear
More informationMATH 128A, SUMMER 2009, FINAL EXAM SOLUTION
MATH 28A, SUMME 2009, FINAL EXAM SOLUTION BENJAMIN JOHNSON () (8 poins) [Lagrange Inerpolaion] (a) (4 poins) Le f be a funcion defined a some real numbers x 0,..., x n. Give a defining equaion for he Lagrange
More information1 Review of Zero-Sum Games
COS 5: heoreical Machine Learning Lecurer: Rob Schapire Lecure #23 Scribe: Eugene Brevdo April 30, 2008 Review of Zero-Sum Games Las ime we inroduced a mahemaical model for wo player zero-sum games. Any
More informationLecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.
Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in
More informationHamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:
M ah 5 7 Fall 9 L ecure O c. 4, 9 ) Hamilon- J acobi Equaion: Weak S oluion We coninue he sudy of he Hamilon-Jacobi equaion: We have shown ha u + H D u) = R n, ) ; u = g R n { = }. ). In general we canno
More informationA problem related to Bárány Grünbaum conjecture
Filoma 27:1 (2013), 109 113 DOI 10.2298/FIL1301109B Published by Faculy of Sciences and Mahemaics, Universiy of Niš, Serbia Available a: hp://www.pmf.ni.ac.rs/filoma A problem relaed o Bárány Grünbaum
More informationSUFFICIENT CONDITIONS FOR EXISTENCE SOLUTION OF LINEAR TWO-POINT BOUNDARY PROBLEM IN MINIMIZATION OF QUADRATIC FUNCTIONAL
HE PUBLISHING HOUSE PROCEEDINGS OF HE ROMANIAN ACADEMY, Series A, OF HE ROMANIAN ACADEMY Volume, Number 4/200, pp 287 293 SUFFICIEN CONDIIONS FOR EXISENCE SOLUION OF LINEAR WO-POIN BOUNDARY PROBLEM IN
More informationTwo Coupled Oscillators / Normal Modes
Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own
More informationSTABILITY OF RESET SWITCHING SYSTEMS
STABILITY OF RESET SWITCHING SYSTEMS J.P. Paxman,G.Vinnicombe Cambridge Universiy Engineering Deparmen, UK, jpp7@cam.ac.uk Keywords: Sabiliy, swiching sysems, bumpless ransfer, conroller realizaion. Absrac
More information4. Advanced Stability Theory
Applied Nonlinear Conrol Nguyen an ien - 4 4 Advanced Sabiliy heory he objecive of his chaper is o presen sabiliy analysis for non-auonomous sysems 41 Conceps of Sabiliy for Non-Auonomous Sysems Equilibrium
More information( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is
UNIT IMPULSE RESPONSE, UNIT STEP RESPONSE, STABILITY. Uni impulse funcion (Dirac dela funcion, dela funcion) rigorously defined is no sricly a funcion, bu disribuion (or measure), precise reamen requires
More informationOptimal Paired Choice Block Designs. Supplementary Material
Saisica Sinica: Supplemen Opimal Paired Choice Block Designs Rakhi Singh 1, Ashish Das 2 and Feng-Shun Chai 3 1 IITB-Monash Research Academy, Mumbai, India 2 Indian Insiue of Technology Bombay, Mumbai,
More informationADDITIONAL PROBLEMS (a) Find the Fourier transform of the half-cosine pulse shown in Fig. 2.40(a). Additional Problems 91
ddiional Problems 9 n inverse relaionship exiss beween he ime-domain and freuency-domain descripions of a signal. Whenever an operaion is performed on he waveform of a signal in he ime domain, a corresponding
More informationExercises: Similarity Transformation
Exercises: Similariy Transformaion Problem. Diagonalize he following marix: A [ 2 4 Soluion. Marix A has wo eigenvalues λ 3 and λ 2 2. Since (i) A is a 2 2 marix and (ii) i has 2 disinc eigenvalues, we
More informationKEY. Math 334 Midterm III Winter 2008 section 002 Instructor: Scott Glasgow
KEY Mah 334 Miderm III Winer 008 secion 00 Insrucor: Sco Glasgow Please do NOT wrie on his exam. No credi will be given for such work. Raher wrie in a blue book, or on your own paper, preferably engineering
More informationL07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms
L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS NA568 Mobile Roboics: Mehods & Algorihms Today s Topic Quick review on (Linear) Kalman Filer Kalman Filering for Non-Linear Sysems Exended Kalman Filer (EKF)
More informationSTATE-SPACE MODELLING. A mass balance across the tank gives:
B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing
More informationOperators related to the Jacobi setting, for all admissible parameter values
Operaors relaed o he Jacobi seing, for all admissible parameer values Peer Sjögren Universiy of Gohenburg Join work wih A. Nowak and T. Szarek Alba, June 2013 () 1 / 18 Le Pn α,β be he classical Jacobi
More informationSecond Order Linear Differential Equations
Second Order Linear Differenial Equaions Second order linear equaions wih consan coefficiens; Fundamenal soluions; Wronskian; Exisence and Uniqueness of soluions; he characerisic equaion; soluions of homogeneous
More informationContinuous Time. Time-Domain System Analysis. Impulse Response. Impulse Response. Impulse Response. Impulse Response. ( t) + b 0.
Time-Domain Sysem Analysis Coninuous Time. J. Robers - All Righs Reserved. Edied by Dr. Rober Akl 1. J. Robers - All Righs Reserved. Edied by Dr. Rober Akl 2 Le a sysem be described by a 2 y ( ) + a 1
More informationLaplace Transform and its Relation to Fourier Transform
Chaper 6 Laplace Transform and is Relaion o Fourier Transform (A Brief Summary) Gis of he Maer 2 Domains of Represenaion Represenaion of signals and sysems Time Domain Coninuous Discree Time Time () [n]
More informationSignal and System (Chapter 3. Continuous-Time Systems)
Signal and Sysem (Chaper 3. Coninuous-Time Sysems) Prof. Kwang-Chun Ho kwangho@hansung.ac.kr Tel: 0-760-453 Fax:0-760-4435 1 Dep. Elecronics and Informaion Eng. 1 Nodes, Branches, Loops A nework wih b
More informationTHE MATRIX-TREE THEOREM
THE MATRIX-TREE THEOREM 1 The Marix-Tree Theorem. The Marix-Tree Theorem is a formula for he number of spanning rees of a graph in erms of he deerminan of a cerain marix. We begin wih he necessary graph-heoreical
More informationMath 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:
Mah 4 Eam Review Problems Problem. Calculae he 3rd Taylor polynomial for arcsin a =. Soluion. Le f() = arcsin. For his problem, we use he formula f() + f () + f ()! + f () 3! for he 3rd Taylor polynomial
More informationLinear Dynamic Models
Linear Dnamic Models and Forecasing Reference aricle: Ineracions beween he muliplier analsis and he principle of acceleraion Ouline. The sae space ssem as an approach o working wih ssems of difference
More informationt + t sin t t cos t sin t. t cos t sin t dt t 2 = exp 2 log t log(t cos t sin t) = Multiplying by this factor and then integrating, we conclude that
ODEs, Homework #4 Soluions. Check ha y ( = is a soluion of he second-order ODE ( cos sin y + y sin y sin = 0 and hen use his fac o find all soluions of he ODE. When y =, we have y = and also y = 0, so
More informationLAPLACE TRANSFORM AND TRANSFER FUNCTION
CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dep. of Chemical and Biological Engineering 5-1 Road Map of he Lecure V Laplace Transform and Transfer funcions
More informationEE363 homework 1 solutions
EE363 Prof. S. Boyd EE363 homework 1 soluions 1. LQR for a riple accumulaor. We consider he sysem x +1 = Ax + Bu, y = Cx, wih 1 1 A = 1 1, B =, C = [ 1 ]. 1 1 This sysem has ransfer funcion H(z) = (z 1)
More informationRemoving Useless Productions of a Context Free Grammar through Petri Net
Journal of Compuer Science 3 (7): 494-498, 2007 ISSN 1549-3636 2007 Science Publicaions Removing Useless Producions of a Conex Free Grammar hrough Peri Ne Mansoor Al-A'ali and Ali A Khan Deparmen of Compuer
More informationBernoulli numbers. Francesco Chiatti, Matteo Pintonello. December 5, 2016
UNIVERSITÁ DEGLI STUDI DI PADOVA, DIPARTIMENTO DI MATEMATICA TULLIO LEVI-CIVITA Bernoulli numbers Francesco Chiai, Maeo Pinonello December 5, 206 During las lessons we have proved he Las Ferma Theorem
More informationFor example, the comb filter generated from. ( ) has a transfer function. e ) has L notches at ω = (2k+1)π/L and L peaks at ω = 2π k/l,
Comb Filers The simple filers discussed so far are characeried eiher by a single passband and/or a single sopband There are applicaions where filers wih muliple passbands and sopbands are required The
More informationRobotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1.
Roboics I April 11, 017 Exercise 1 he kinemaics of a 3R spaial robo is specified by he Denavi-Harenberg parameers in ab 1 i α i d i a i θ i 1 π/ L 1 0 1 0 0 L 3 0 0 L 3 3 able 1: able of DH parameers of
More informationElements of Computer Graphics
CS580: Compuer Graphics Min H. Kim KAIST School of Compuing Elemens of Compuer Graphics Geomery Maerial model Ligh Rendering Virual phoography 2 Foundaions of Compuer Graphics A PINHOLE CAMERA IN 3D 3
More information3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon
3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of
More informationProperties Of Solutions To A Generalized Liénard Equation With Forcing Term
Applied Mahemaics E-Noes, 8(28), 4-44 c ISSN 67-25 Available free a mirror sies of hp://www.mah.nhu.edu.w/ amen/ Properies Of Soluions To A Generalized Liénard Equaion Wih Forcing Term Allan Kroopnick
More informationHomework 6 AERE331 Spring 2019 Due 4/24(W) Name Sec. 1 / 2
Homework 6 AERE33 Spring 9 Due 4/4(W) Name Sec / PROBLEM (5p In PROBLEM 4 of HW4 we used he frequency domain o design a yaw/rudder feedback conrol sysem for a plan wih ransfer funcion 46 Gp () s The conroller
More informationLecture 9: September 25
0-725: Opimizaion Fall 202 Lecure 9: Sepember 25 Lecurer: Geoff Gordon/Ryan Tibshirani Scribes: Xuezhi Wang, Subhodeep Moira, Abhimanu Kumar Noe: LaTeX emplae couresy of UC Berkeley EECS dep. Disclaimer:
More informationDifferential Equations
Mah 21 (Fall 29) Differenial Equaions Soluion #3 1. Find he paricular soluion of he following differenial equaion by variaion of parameer (a) y + y = csc (b) 2 y + y y = ln, > Soluion: (a) The corresponding
More informationComputer-Aided Analysis of Electronic Circuits Course Notes 3
Gheorghe Asachi Technical Universiy of Iasi Faculy of Elecronics, Telecommunicaions and Informaion Technologies Compuer-Aided Analysis of Elecronic Circuis Course Noes 3 Bachelor: Telecommunicaion Technologies
More informationSection 4.4 Logarithmic Properties
Secion. Logarihmic Properies 5 Secion. Logarihmic Properies In he previous secion, we derived wo imporan properies of arihms, which allowed us o solve some asic eponenial and arihmic equaions. Properies
More informationSMT 2014 Calculus Test Solutions February 15, 2014 = 3 5 = 15.
SMT Calculus Tes Soluions February 5,. Le f() = and le g() =. Compue f ()g (). Answer: 5 Soluion: We noe ha f () = and g () = 6. Then f ()g () =. Plugging in = we ge f ()g () = 6 = 3 5 = 5.. There is a
More informationKEY. Math 334 Midterm III Fall 2008 sections 001 and 003 Instructor: Scott Glasgow
KEY Mah 334 Miderm III Fall 28 secions and 3 Insrucor: Sco Glasgow Please do NOT wrie on his exam. No credi will be given for such work. Raher wrie in a blue book, or on your own paper, preferably engineering
More informationLECTURE 1: GENERALIZED RAY KNIGHT THEOREM FOR FINITE MARKOV CHAINS
LECTURE : GENERALIZED RAY KNIGHT THEOREM FOR FINITE MARKOV CHAINS We will work wih a coninuous ime reversible Markov chain X on a finie conneced sae space, wih generaor Lf(x = y q x,yf(y. (Recall ha q
More informationUnit Root Time Series. Univariate random walk
Uni Roo ime Series Univariae random walk Consider he regression y y where ~ iid N 0, he leas squares esimae of is: ˆ yy y y yy Now wha if = If y y hen le y 0 =0 so ha y j j If ~ iid N 0, hen y ~ N 0, he
More informationMA 214 Calculus IV (Spring 2016) Section 2. Homework Assignment 1 Solutions
MA 14 Calculus IV (Spring 016) Secion Homework Assignmen 1 Soluions 1 Boyce and DiPrima, p 40, Problem 10 (c) Soluion: In sandard form he given firs-order linear ODE is: An inegraing facor is given by
More information4 Sequences of measurable functions
4 Sequences of measurable funcions 1. Le (Ω, A, µ) be a measure space (complee, afer a possible applicaion of he compleion heorem). In his chaper we invesigae relaions beween various (nonequivalen) convergences
More informationGeneralized Least Squares
Generalized Leas Squares Augus 006 1 Modified Model Original assumpions: 1 Specificaion: y = Xβ + ε (1) Eε =0 3 EX 0 ε =0 4 Eεε 0 = σ I In his secion, we consider relaxing assumpion (4) Insead, assume
More informationPhysics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle
Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,
More informationVanishing Viscosity Method. There are another instructive and perhaps more natural discontinuous solutions of the conservation law
Vanishing Viscosiy Mehod. There are anoher insrucive and perhaps more naural disconinuous soluions of he conservaion law (1 u +(q(u x 0, he so called vanishing viscosiy mehod. This mehod consiss in viewing
More informationSecond quantization and gauge invariance.
1 Second quanizaion and gauge invariance. Dan Solomon Rauland-Borg Corporaion Moun Prospec, IL Email: dsolom@uic.edu June, 1. Absrac. I is well known ha he single paricle Dirac equaion is gauge invarian.
More informationODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004
ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform
More informationMath 10B: Mock Mid II. April 13, 2016
Name: Soluions Mah 10B: Mock Mid II April 13, 016 1. ( poins) Sae, wih jusificaion, wheher he following saemens are rue or false. (a) If a 3 3 marix A saisfies A 3 A = 0, hen i canno be inverible. True.
More informationIntroduction to Simulation - Lecture 16. Methods for Computing Periodic Steady-State - Part II Jacob White
Inroducion o Simuaion - ecure 6 Mehods for Compuing Periodic Seady-Sae - Par II Jacob Whie hanks o Deepak Ramaswamy, Micha Rewienski, and Karen Veroy Ouine hree Mehods so far ime inegraion uni seady-sae
More informationSOLUTIONS TO ECE 3084
SOLUTIONS TO ECE 384 PROBLEM 2.. For each sysem below, specify wheher or no i is: (i) memoryless; (ii) causal; (iii) inverible; (iv) linear; (v) ime invarian; Explain your reasoning. If he propery is no
More informationPade and Laguerre Approximations Applied. to the Active Queue Management Model. of Internet Protocol
Applied Mahemaical Sciences, Vol. 7, 013, no. 16, 663-673 HIKARI Ld, www.m-hikari.com hp://dx.doi.org/10.1988/ams.013.39499 Pade and Laguerre Approximaions Applied o he Acive Queue Managemen Model of Inerne
More informationOrthogonal Rational Functions, Associated Rational Functions And Functions Of The Second Kind
Proceedings of he World Congress on Engineering 2008 Vol II Orhogonal Raional Funcions, Associaed Raional Funcions And Funcions Of The Second Kind Karl Deckers and Adhemar Bulheel Absrac Consider he sequence
More informationTHE GENERALIZED PASCAL MATRIX VIA THE GENERALIZED FIBONACCI MATRIX AND THE GENERALIZED PELL MATRIX
J Korean Mah Soc 45 008, No, pp 479 49 THE GENERALIZED PASCAL MATRIX VIA THE GENERALIZED FIBONACCI MATRIX AND THE GENERALIZED PELL MATRIX Gwang-yeon Lee and Seong-Hoon Cho Reprined from he Journal of he
More informationLatent Spaces and Matrix Factorization
Compuaional Linguisics Laen Spaces and Marix Facorizaion Dierich Klakow FR 4.7 Allgemeine Linguisik (Compuerlinguisik) Universiä des Saarlandes Summer 0 Goal Goal: rea documen clusering and word clusering
More informationTHE BERNOULLI NUMBERS. t k. = lim. = lim = 1, d t B 1 = lim. 1+e t te t = lim t 0 (e t 1) 2. = lim = 1 2.
THE BERNOULLI NUMBERS The Bernoulli numbers are defined here by he exponenial generaing funcion ( e The firs one is easy o compue: (2 and (3 B 0 lim 0 e lim, 0 e ( d B lim 0 d e +e e lim 0 (e 2 lim 0 2(e
More information= ( ) ) or a system of differential equations with continuous parametrization (T = R
XIII. DIFFERENCE AND DIFFERENTIAL EQUATIONS Ofen funcions, or a sysem of funcion, are paramerized in erms of some variable, usually denoed as and inerpreed as ime. The variable is wrien as a funcion of
More informationMATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence
MATH 433/533, Fourier Analysis Secion 6, Proof of Fourier s Theorem for Poinwise Convergence Firs, some commens abou inegraing periodic funcions. If g is a periodic funcion, g(x + ) g(x) for all real x,
More informationWe just finished the Erdős-Stone Theorem, and ex(n, F ) (1 1/(χ(F ) 1)) ( n
Lecure 3 - Kövari-Sós-Turán Theorem Jacques Versraëe jacques@ucsd.edu We jus finished he Erdős-Sone Theorem, and ex(n, F ) ( /(χ(f ) )) ( n 2). So we have asympoics when χ(f ) 3 bu no when χ(f ) = 2 i.e.
More informationdt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.
Mah Rahman Exam Review Soluions () Consider he IVP: ( 4)y 3y + 4y = ; y(3) = 0, y (3) =. (a) Please deermine he longes inerval for which he IVP is guaraneed o have a unique soluion. Soluion: The disconinuiies
More informationChapter 3 Boundary Value Problem
Chaper 3 Boundary Value Problem A boundary value problem (BVP) is a problem, ypically an ODE or a PDE, which has values assigned on he physical boundary of he domain in which he problem is specified. Le
More informationIntroduction D P. r = constant discount rate, g = Gordon Model (1962): constant dividend growth rate.
Inroducion Gordon Model (1962): D P = r g r = consan discoun rae, g = consan dividend growh rae. If raional expecaions of fuure discoun raes and dividend growh vary over ime, so should he D/P raio. Since
More informationLet us start with a two dimensional case. We consider a vector ( x,
Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our
More informationSection 4.4 Logarithmic Properties
Secion. Logarihmic Properies 59 Secion. Logarihmic Properies In he previous secion, we derived wo imporan properies of arihms, which allowed us o solve some asic eponenial and arihmic equaions. Properies
More informationBEng (Hons) Telecommunications. Examinations for / Semester 2
BEng (Hons) Telecommunicaions Cohor: BTEL/14/FT Examinaions for 2015-2016 / Semeser 2 MODULE: ELECTROMAGNETIC THEORY MODULE CODE: ASE2103 Duraion: 2 ½ Hours Insrucions o Candidaes: 1. Answer ALL 4 (FOUR)
More information4/9/2012. Signals and Systems KX5BQY EE235. Today s menu. System properties
Signals and Sysems hp://www.youube.com/v/iv6fo KX5BQY EE35 oday s menu Good weeend? Sysem properies iy Superposiion! Sysem properies iy: A Sysem is if i mees he following wo crieria: If { x( )} y( ) and
More informationMath Week 15: Section 7.4, mass-spring systems. These are notes for Monday. There will also be course review notes for Tuesday, posted later.
Mah 50-004 Week 5: Secion 7.4, mass-spring sysems. These are noes for Monday. There will also be course review noes for Tuesday, posed laer. Mon Apr 3 7.4 mass-spring sysems. Announcemens: Warm up exercise:
More informationMATH 2050 Assignment 9 Winter Do not need to hand in. 1. Find the determinant by reducing to triangular form for the following matrices.
MATH 2050 Assignmen 9 Winer 206 Do no need o hand in Noe ha he final exam also covers maerial afer HW8, including, for insance, calculaing deerminan by row operaions, eigenvalues and eigenvecors, similariy
More information11!Hí MATHEMATICS : ERDŐS AND ULAM PROC. N. A. S. of decomposiion, properly speaking) conradics he possibiliy of defining a counably addiive real-valu
ON EQUATIONS WITH SETS AS UNKNOWNS BY PAUL ERDŐS AND S. ULAM DEPARTMENT OF MATHEMATICS, UNIVERSITY OF COLORADO, BOULDER Communicaed May 27, 1968 We shall presen here a number of resuls in se heory concerning
More informationPredator - Prey Model Trajectories and the nonlinear conservation law
Predaor - Prey Model Trajecories and he nonlinear conservaion law James K. Peerson Deparmen of Biological Sciences and Deparmen of Mahemaical Sciences Clemson Universiy Ocober 28, 213 Ouline Drawing Trajecories
More informationSome Ramsey results for the n-cube
Some Ramsey resuls for he n-cube Ron Graham Universiy of California, San Diego Jozsef Solymosi Universiy of Briish Columbia, Vancouver, Canada Absrac In his noe we esablish a Ramsey-ype resul for cerain
More informationLatent Spaces and Matrix Factorization
Compuaional Linguisics Laen Spaces and Marix Facorizaion Sefan Thaer & Dierich Klakow FR 4.7 Allgemeine Linguisik (Compuerlinguisik) Universiä des Saarlandes Summer 2013 Goal Goal: rea documen clusering
More informationChapter Three Systems of Linear Differential Equations
Chaper Three Sysems of Linear Differenial Equaions In his chaper we are going o consier sysems of firs orer orinary ifferenial equaions. These are sysems of he form x a x a x a n x n x a x a x a n x n
More informationSolutions from Chapter 9.1 and 9.2
Soluions from Chaper 9 and 92 Secion 9 Problem # This basically boils down o an exercise in he chain rule from calculus We are looking for soluions of he form: u( x) = f( k x c) where k x R 3 and k is
More information