CDS 101: Lecture 5.1 Reachability and State Space Feedback

Size: px
Start display at page:

Download "CDS 101: Lecture 5.1 Reachability and State Space Feedback"

Transcription

1 CDS, Lctur 5. CDS : Lctur 5. Rachability ad Stat Spac Fdback Richard M. Murray ad Hido Mabuchi 5 Octobr 4 Goals: Di rachability o a cotrol systm Giv tsts or rachability o liar systms ad apply to ampls Dscrib th dsig o stat dback cotrollrs or liar systms Radig: Åström ad Murray, Aalysis ad Dsig o Fdback Systms, Ch 5 Lctur Rviw 4.: rom Liar Last Systms Wk u = A + Bu ( = y t ( = At A( t τ ( + ( τ τ + ( τ = yt C C Bu d t Proprtis o liar systms Liarity with rspct to iitial coditio ad iputs Stability charactrizd by igvalus May applicatios ad tools availabl Provid local dscriptio or oliar systms 5 Oct 4 H.Mabuchi, Caltch CDS 5 Octobr 4

2 CDS, Lctur 5. Cotrol Dsig Cocpts Systm dscriptio: sigl iput, sigl output systm (MIMO also OK = (, u R, ( giv y = h(, u u, y Stability: stabiliz th systm aroud a quilibrium poit Giv quilibrium poit R, id cotrol law u=α( such that lim t ( = or all ( Rachability: str th systm btw two poits Giv, R, id a iput u(t such that = u t t = T = (, ( taks ( ( Trackig: track a giv output trajctory Giv y d (t, id u=α(,t such that ( lim yt ( y( t = or all ( d y(t y d (t t 5 Oct 4 H.Mabuchi, Caltch CDS 3 Rachability o Iput/Output Systms = (, u y = h(, u D A iput/output systm is rachabl i or ay, R ad ay tim T > thr ists a iput u:[,t] R such that th solutio o th dyamics startig rom (= ad applyig iput u(t givs (T=. Rmarks I th diitio, ad do ot hav to b quilibrium poits w do t cssarily stay at atr tim T. Rachability is did i trms o stats dos t dpd o output For liar systms, ca charactriz rachability by lookig at th gral solutio: T = A + Bu AT A( T τ ( T = + Bu( τ dτ τ = I itgral is surjctiv (as a liar oprator, th w ca id a iput to achiv ay dsird ial stat. 5 Oct 4 H.Mabuchi, Caltch CDS 4 5 Octobr 4

3 CDS, Lctur 5. = A + Bu Tsts or Rachability T & AT A( T τ = + τ = ( T Bu( τ dτ Thm A liar systm is rachabl i ad oly i th rachability matri is ull rak. Rmarks Vry simpl tst to apply. I MATLAB, us ctrb(a,b ad chck rak w/ dt( I this tst is satisid, w say th pair (A,B is rachabl Som isight ito th proo ca b s by padig th matri potial 5 Oct 4 H.Mabuchi, Caltch CDS 5 Eampl #: Liarizd pdulum o a cart θ m Qustio: ca w locally cotrol th positio o th cart by propr choic o iput? F M Approach: look at th liarizatio aroud th upright positio (good approimatio to th ull dyamics i θ rmais small 5 Oct 4 H.Mabuchi, Caltch CDS 6 5 Octobr 4 3

4 CDS, Lctur 5. Eampl #, co t: Liarizd pdulum o a cart m θ u M Rachability matri B AB A B A 3 B Simpliy by sttig b = Full rak as log as costats ar such that colums ad 3 ar ot multipls o ach othr rachabl as log as g(m+m ca str liarizatio btw poits by propr choic o iput 5 Oct 4 H.Mabuchi, Caltch CDS 7 Cotrol Dsig Cocpts Systm dscriptio: sigl iput, sigl output systm (MIMO also OK = (, u y = h(, u Stability: stabiliz th systm aroud a quilibrium poit Giv quilibrium poit R, id cotrol law u=α( such that lim t ( = or all ( Rachability: str th systm btw two poits Giv, R, id a iput u(t such that = (, u( t taks ( t = ( t = Trackig: track a giv output trajctory Giv y d (t, id u=α(,t such that ( lim yt ( y( t = or all ( d y(t y d (t t 5 Oct 4 H.Mabuchi, Caltch CDS 8 5 Octobr 4 4

5 CDS, Lctur 5. Stat spac cotrollr dsig or liar systms = A + Bu T ( = AT A( T τ + ( τ τ = T Bu dτ Goal: id a liar cotrol law u=k such that th closd loop systm = A + BK = ( A + BK is stabl at =. Rmarks Stability basd o igvalus us K to mak igvalus o (A+BK stabl Ca also lik igvalus to prormac (g, iitial coditio rspos Qustio: wh ca w plac th igvalus ayplac that w wat? Thorm Th igvalus o (A+BK ca b st to arbitrary valus i ad oly i th pair (A,B is rachabl. MATLAB: K = plac(a, B, igs 5 Oct 4 H.Mabuchi, Caltch CDS 9 Natural dyamics = br a = a d Eampl #: Prdator pry Cotrolld dyamics: modulat ood supply = br ( + u a = a d Q: ca w mov rom som iitial populatio o os ad rabbits to a spciid o i tim T by modulatio i th ood supply? Q: ca w stabiliz th populatio aroud th dsird quilibrum poit Approach: try to aswr this qustio locally, aroud th atural quilibrium poit Oct 4 H.Mabuchi, Caltch CDS ustabl stabl 5 Octobr 4 5

6 CDS, Lctur 5. Eampl #: Problm stup Equilibrium poit calculatio = br ( + u a = a d = (5, 35 Liarizatio Comput liarizatio aroud quil. poit, : A= B = u (, u (, u % Comput th quil poit % prdpry.m cotais dyamics = ili('prdpry(,'; q = solv(, [5,5]; % Comput liarizatio A = [ br - a*q( - a*q(; a*q(, -d + a*q( ]; B = [br*q(; ]; Rdi local variabls: z=-, v=u-u d z br a, a, z b r, v dt z = a, d a, z + + Rachabl? YES, i b r,a (chck [B AB] ca locally str to ay poit 5 Oct 4 H.Mabuchi, Caltch CDS Eampl #: Stabilizatio via igvalu assigmt d z br a, a, z b r, v dt z = + a, d + a, z Cotrol dsig: v = Kz = K( u = u + v = u K( Plac pols at stabl valus Choos λ=-, - K = plac(a, B, [-; -]; Modiy dyamics to iclud cotrol = b ( K( a = a d r stabl ustabl Oct 4 H.Mabuchi, Caltch CDS 5 Octobr 4 6

7 CDS, Lctur 5. Implmtatio Dtails Eigvalus dtrmi prormac For ach igvalu λ i =σ i + jω i, gt cotributio o th orm ( ω ω σt y ( t = asi( t + bcos( t i Imag Ais - Pol-zro map ω p Amplitud To: Y( Stp Rspos From: U( T π/ω p Rpatd igvalus ca giv additioal trms o th orm t k σ +jω Ral Ais Tim (sc. Us stimator to dtrmi th currt stat i you ca t masur it u = A+ Bu y Estimator looks at iputs ad outputs o plat ad stimats th currt stat Ca show that i a systm is obsrvabl th you ca costruct ad stimator Us th stimatd stat as th dback Estimator ˆ u = Kˆ Kalma iltr is a ampl o a stimator 5 Oct 4 H.Mabuchi, Caltch CDS 3 Summary: Rachability ad Stat Spac Fdback = A + Bu u = u K( Ky cocpts Rachability: id u s.t. Rachability rak tst or liar systms Stat dback to assig igvalus Oct 4 H.Mabuchi, Caltch CDS 4 5 Octobr 4 7

CDS 101: Lecture 5.1 Reachability and State Space Feedback

CDS 101: Lecture 5.1 Reachability and State Space Feedback CDS, Lctur 5. CDS : Lctur 5. Rachability ad Stat Spac Fdback Richard M. Murray 7 Octobr 3 Goals: Di rachability o a cotrol systm Giv tsts or rachability o liar systms ad apply to ampls Dscrib th dsig o

More information

CDS 101: Lecture 5.1 Controllability and State Space Feedback

CDS 101: Lecture 5.1 Controllability and State Space Feedback CDS, Lecture 5. CDS : Lecture 5. Cotrollability ad State Space Feedback Richard M. Murray 8 October Goals: Deie cotrollability o a cotrol system Give tests or cotrollability o liear systems ad apply to

More information

1985 AP Calculus BC: Section I

1985 AP Calculus BC: Section I 985 AP Calculus BC: Sctio I 9 Miuts No Calculator Nots: () I this amiatio, l dots th atural logarithm of (that is, logarithm to th bas ). () Ulss othrwis spcifid, th domai of a fuctio f is assumd to b

More information

Chapter Taylor Theorem Revisited

Chapter Taylor Theorem Revisited Captr 0.07 Taylor Torm Rvisitd Atr radig tis captr, you sould b abl to. udrstad t basics o Taylor s torm,. writ trascdtal ad trigoomtric uctios as Taylor s polyomial,. us Taylor s torm to id t valus o

More information

CDS 101: Lecture 5-1 Reachability and State Space Feedback

CDS 101: Lecture 5-1 Reachability and State Space Feedback CDS 11: Lecture 5-1 Reachability and State Space Feedback Richard M. Murray 23 October 26 Goals: Define reachability of a control system Give tests for reachability of linear systems and apply to examples

More information

Worksheet: Taylor Series, Lagrange Error Bound ilearnmath.net

Worksheet: Taylor Series, Lagrange Error Bound ilearnmath.net Taylor s Thorm & Lagrag Error Bouds Actual Error This is th ral amout o rror, ot th rror boud (worst cas scario). It is th dirc btw th actual () ad th polyomial. Stps:. Plug -valu ito () to gt a valu.

More information

MONTGOMERY COLLEGE Department of Mathematics Rockville Campus. 6x dx a. b. cos 2x dx ( ) 7. arctan x dx e. cos 2x dx. 2 cos3x dx

MONTGOMERY COLLEGE Department of Mathematics Rockville Campus. 6x dx a. b. cos 2x dx ( ) 7. arctan x dx e. cos 2x dx. 2 cos3x dx MONTGOMERY COLLEGE Dpartmt of Mathmatics Rockvill Campus MATH 8 - REVIEW PROBLEMS. Stat whthr ach of th followig ca b itgratd by partial fractios (PF), itgratio by parts (PI), u-substitutio (U), or o of

More information

PURE MATHEMATICS A-LEVEL PAPER 1

PURE MATHEMATICS A-LEVEL PAPER 1 -AL P MATH PAPER HONG KONG EXAMINATIONS AUTHORITY HONG KONG ADVANCED LEVEL EXAMINATION PURE MATHEMATICS A-LEVEL PAPER 8 am am ( hours) This papr must b aswrd i Eglish This papr cosists of Sctio A ad Sctio

More information

Option 3. b) xe dx = and therefore the series is convergent. 12 a) Divergent b) Convergent Proof 15 For. p = 1 1so the series diverges.

Option 3. b) xe dx = and therefore the series is convergent. 12 a) Divergent b) Convergent Proof 15 For. p = 1 1so the series diverges. Optio Chaptr Ercis. Covrgs to Covrgs to Covrgs to Divrgs Covrgs to Covrgs to Divrgs 8 Divrgs Covrgs to Covrgs to Divrgs Covrgs to Covrgs to Covrgs to Covrgs to 8 Proof Covrgs to π l 8 l a b Divrgt π Divrgt

More information

Lectures 9 IIR Systems: First Order System

Lectures 9 IIR Systems: First Order System EE3054 Sigals ad Systms Lcturs 9 IIR Systms: First Ordr Systm Yao Wag Polytchic Uivrsity Som slids icludd ar xtractd from lctur prstatios prpard by McCllla ad Schafr Lics Ifo for SPFirst Slids This work

More information

Chapter 2 Infinite Series Page 1 of 11. Chapter 2 : Infinite Series

Chapter 2 Infinite Series Page 1 of 11. Chapter 2 : Infinite Series Chatr Ifiit Sris Pag of Sctio F Itgral Tst Chatr : Ifiit Sris By th d of this sctio you will b abl to valuat imror itgrals tst a sris for covrgc by alyig th itgral tst aly th itgral tst to rov th -sris

More information

Review Exercises. 1. Evaluate using the definition of the definite integral as a Riemann Sum. Does the answer represent an area? 2

Review Exercises. 1. Evaluate using the definition of the definite integral as a Riemann Sum. Does the answer represent an area? 2 MATHEMATIS --RE Itgral alculus Marti Huard Witr 9 Rviw Erciss. Evaluat usig th dfiitio of th dfiit itgral as a Rima Sum. Dos th aswr rprst a ara? a ( d b ( d c ( ( d d ( d. Fid f ( usig th Fudamtal Thorm

More information

ECE602 Exam 1 April 5, You must show ALL of your work for full credit.

ECE602 Exam 1 April 5, You must show ALL of your work for full credit. ECE62 Exam April 5, 27 Nam: Solution Scor: / This xam is closd-book. You must show ALL of your work for full crdit. Plas rad th qustions carfully. Plas chck your answrs carfully. Calculators may NOT b

More information

Problem Value Score Earned No/Wrong Rec -3 Total

Problem Value Score Earned No/Wrong Rec -3 Total GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL & COMPUTER ENGINEERING ECE6 Fall Quiz # Writt Eam Novmr, NAME: Solutio Kys GT Usram: LAST FIRST.g., gtiit Rcitatio Sctio: Circl t dat & tim w your Rcitatio

More information

Washington State University

Washington State University he 3 Ktics ad Ractor Dsig Sprg, 00 Washgto Stat Uivrsity Dpartmt of hmical Egrg Richard L. Zollars Exam # You will hav o hour (60 muts) to complt this xam which cosists of four (4) problms. You may us

More information

Exercises for lectures 23 Discrete systems

Exercises for lectures 23 Discrete systems Exrciss for lcturs 3 Discrt systms Michal Šbk Automatické říí 06 30-4-7 Stat-Spac a Iput-Output scriptios Automatické říí - Kybrtika a robotika Mols a trasfrs i CSTbx >> F=[ ; 3 4]; G=[ ;]; H=[ ]; J=0;

More information

Addition of angular momentum

Addition of angular momentum Addition of angular momntum April, 0 Oftn w nd to combin diffrnt sourcs of angular momntum to charactriz th total angular momntum of a systm, or to divid th total angular momntum into parts to valuat th

More information

First derivative analysis

First derivative analysis Robrto s Nots on Dirntial Calculus Chaptr 8: Graphical analysis Sction First drivativ analysis What you nd to know alrady: How to us drivativs to idntiy th critical valus o a unction and its trm points

More information

Periodic Structures. Filter Design by the Image Parameter Method

Periodic Structures. Filter Design by the Image Parameter Method Prioic Structurs a Filtr sig y th mag Paramtr Mtho ECE53: Microwav Circuit sig Pozar Chaptr 8, Sctios 8. & 8. Josh Ottos /4/ Microwav Filtrs (Chaptr Eight) microwav filtr is a two-port twork us to cotrol

More information

CDS 101: Lecture 8.2 Tools for PID & Loop Shaping

CDS 101: Lecture 8.2 Tools for PID & Loop Shaping CDS : Lecture 8. Tools for PID & Loop Shapig Richard M. Murray 7 November 4 Goals: Show how to use loop shapig to achieve a performace specificatio Itroduce ew tools for loop shapig desig: Ziegler-Nichols,

More information

Frequency Measurement in Noise

Frequency Measurement in Noise Frqucy Masurmt i ois Porat Sctio 6.5 /4 Frqucy Mas. i ois Problm Wat to o look at th ct o ois o usig th DFT to masur th rqucy o a siusoid. Cosidr sigl complx siusoid cas: j y +, ssum Complx Whit ois Gaussia,

More information

15/03/1439. Lectures on Signals & systems Engineering

15/03/1439. Lectures on Signals & systems Engineering Lcturs o Sigals & syms Egirig Dsigd ad Prd by Dr. Ayma Elshawy Elsfy Dpt. of Syms & Computr Eg. Al-Azhar Uivrsity Email : aymalshawy@yahoo.com A sigal ca b rprd as a liar combiatio of basic sigals. Th

More information

CS537. Numerical Analysis and Computing

CS537. Numerical Analysis and Computing CS57 Numerical Aalysis ad Computig Lecture Locatig Roots o Equatios Proessor Ju Zhag Departmet o Computer Sciece Uiversity o Ketucky Leigto KY 456-6 Jauary 9 9 What is the Root May physical system ca be

More information

CDS 101: Lecture 4.1 Linear Systems

CDS 101: Lecture 4.1 Linear Systems CDS : Lecture 4. Linear Systems Richard M. Murray 8 October 4 Goals: Describe linear system models: properties, eamples, and tools Characterize stability and performance of linear systems in terms of eigenvalues

More information

Theory of Control: I. Overview. Types of Structural Control Active Control Systems. Types of Structural Control

Theory of Control: I. Overview. Types of Structural Control Active Control Systems. Types of Structural Control Thory of Cotrol: I Asia Pacific Summr School o Smart s Tchology Richard Christso Uivrsity of Cocticut Ovrviw Itroductio to structural cotrol Cotrol thory Basic fdback cotrol Optimal cotrol stat fdback

More information

Frequency Response & Digital Filters

Frequency Response & Digital Filters Frquy Rspos & Digital Filtrs S Wogsa Dpt. of Cotrol Systms ad Istrumtatio Egirig, KUTT Today s goals Frquy rspos aalysis of digital filtrs LTI Digital Filtrs Digital filtr rprstatios ad struturs Idal filtrs

More information

Probability & Statistics,

Probability & Statistics, Probability & Statistics, BITS Pilai K K Birla Goa Campus Dr. Jajati Kshari Sahoo Dpartmt of Mathmatics BITS Pilai, K K Birla Goa Campus Poisso Distributio Poisso Distributio: A radom variabl X is said

More information

2617 Mark Scheme June 2005 Mark Scheme 2617 June 2005

2617 Mark Scheme June 2005 Mark Scheme 2617 June 2005 Mark Schm 67 Ju 5 GENERAL INSTRUCTIONS Marks i th mark schm ar plicitly dsigatd as M, A, B, E or G. M marks ("mthod" ar for a attmpt to us a corrct mthod (ot mrly for statig th mthod. A marks ("accuracy"

More information

Introduction to Quantum Information Processing. Overview. A classical randomised algorithm. q 3,3 00 0,0. p 0,0. Lecture 10.

Introduction to Quantum Information Processing. Overview. A classical randomised algorithm. q 3,3 00 0,0. p 0,0. Lecture 10. Itroductio to Quatum Iformatio Procssig Lctur Michl Mosca Ovrviw! Classical Radomizd vs. Quatum Computig! Dutsch-Jozsa ad Brsti- Vazirai algorithms! Th quatum Fourir trasform ad phas stimatio A classical

More information

Control systems

Control systems Last tim,.5 Cotrol sstms Cotrollabilit ad obsrvabilit (Chaptr ) Two approahs to stat fdbak dsig (Chaptr 8) Usig otrollabl aoial form B solvig matri quatios Toda, w otiu to work o fdbak dsig (Chaptr 8)

More information

Solution to 1223 The Evil Warden.

Solution to 1223 The Evil Warden. Solutio to 1 Th Evil Ward. This is o of thos vry rar PoWs (I caot thik of aothr cas) that o o solvd. About 10 of you submittd th basic approach, which givs a probability of 47%. I was shockd wh I foud

More information

CS321. Numerical Analysis and Computing

CS321. Numerical Analysis and Computing CS Numerical Aalysis ad Computig Lecture Locatig Roots o Equatios Proessor Ju Zhag Departmet o Computer Sciece Uiversity o Ketucky Leigto KY 456-6 September 8 5 What is the Root May physical system ca

More information

Time : 1 hr. Test Paper 08 Date 04/01/15 Batch - R Marks : 120

Time : 1 hr. Test Paper 08 Date 04/01/15 Batch - R Marks : 120 Tim : hr. Tst Papr 8 D 4//5 Bch - R Marks : SINGLE CORRECT CHOICE TYPE [4, ]. If th compl umbr z sisfis th coditio z 3, th th last valu of z is qual to : z (A) 5/3 (B) 8/3 (C) /3 (D) o of ths 5 4. Th itgral,

More information

Data Assimilation 1. Alan O Neill National Centre for Earth Observation UK

Data Assimilation 1. Alan O Neill National Centre for Earth Observation UK Data Assimilation 1 Alan O Nill National Cntr for Earth Obsrvation UK Plan Motivation & basic idas Univariat (scalar) data assimilation Multivariat (vctor) data assimilation 3d-Variational Mthod (& optimal

More information

Last time: Ground rules for filtering and control system design

Last time: Ground rules for filtering and control system design 6.3 Stochatic Etimatio ad Cotrol, Fall 004 Lecture 7 Lat time: Groud rule for filterig ad cotrol ytem deig Gral ytem Sytem parameter are cotaied i w( t ad w ( t. Deired output i grated by takig the igal

More information

Ordinary Differential Equations

Ordinary Differential Equations Ordiary Diffrtial Equatio Aftr radig thi chaptr, you hould b abl to:. dfi a ordiary diffrtial quatio,. diffrtiat btw a ordiary ad partial diffrtial quatio, ad. Solv liar ordiary diffrtial quatio with fid

More information

Continuous Models for Eigenvalue Problems

Continuous Models for Eigenvalue Problems Cotiuous Models for Eigevalue Problems Li-Zhi Liao (http://www.math.hkbu.edu.hk/~liliao) Departmet of Mathematics Hog Kog Baptist Uiversity Dedicate this talk to the memory of Gee Golub Program for Gee

More information

Technical Support Document Bias of the Minimum Statistic

Technical Support Document Bias of the Minimum Statistic Tchical Support Documt Bias o th Miimum Stattic Itroductio Th papr pla how to driv th bias o th miimum stattic i a radom sampl o siz rom dtributios with a shit paramtr (also kow as thrshold paramtr. Ths

More information

ln x = n e = 20 (nearest integer)

ln x = n e = 20 (nearest integer) H JC Prlim Solutios 6 a + b y a + b / / dy a b 3/ d dy a b at, d Giv quatio of ormal at is y dy ad y wh. d a b () (,) is o th curv a+ b () y.9958 Qustio Solvig () ad (), w hav a, b. Qustio d.77 d d d.77

More information

Addition of angular momentum

Addition of angular momentum Addition of angular momntum April, 07 Oftn w nd to combin diffrnt sourcs of angular momntum to charactriz th total angular momntum of a systm, or to divid th total angular momntum into parts to valuat

More information

APPENDIX: STATISTICAL TOOLS

APPENDIX: STATISTICAL TOOLS I. Nots o radom samplig Why do you d to sampl radomly? APPENDI: STATISTICAL TOOLS I ordr to masur som valu o a populatio of orgaisms, you usually caot masur all orgaisms, so you sampl a subst of th populatio.

More information

Discrete Fourier Transform. Nuno Vasconcelos UCSD

Discrete Fourier Transform. Nuno Vasconcelos UCSD Discrt Fourir Trasform uo Vascoclos UCSD Liar Shift Ivariat (LSI) systms o of th most importat cocpts i liar systms thory is that of a LSI systm Dfiitio: a systm T that maps [ ito y[ is LSI if ad oly if

More information

Mixed Mode Oscillations as a Mechanism for Pseudo-Plateau Bursting

Mixed Mode Oscillations as a Mechanism for Pseudo-Plateau Bursting Mixd Mod Oscillatios as a Mchaism for Psudo-Platau Burstig Richard Brtram Dpartmt of Mathmatics Florida Stat Uivrsity Tallahass, FL Collaborators ad Support Thodor Vo Marti Wchslbrgr Joël Tabak Uivrsity

More information

NET/JRF, GATE, IIT JAM, JEST, TIFR

NET/JRF, GATE, IIT JAM, JEST, TIFR Istitut for NET/JRF, GATE, IIT JAM, JEST, TIFR ad GRE i PHYSICAL SCIENCES Mathmatical Physics JEST-6 Q. Giv th coditio φ, th solutio of th quatio ψ φ φ is giv by k. kφ kφ lφ kφ lφ (a) ψ (b) ψ kφ (c) ψ

More information

Consider serial transmission. In Proakis notation, we receive

Consider serial transmission. In Proakis notation, we receive 5..3 Dciio-Dirctd Pha Trackig [P 6..4] 5.-1 Trackr commoly work o radom data igal (plu oi), o th kow-igal modl do ot apply. W till kow much about th tructur o th igal, though, ad w ca xploit it. Coidr

More information

LECTURE 13 Filling the bands. Occupancy of Available Energy Levels

LECTURE 13 Filling the bands. Occupancy of Available Energy Levels LUR 3 illig th bads Occupacy o Availabl rgy Lvls W hav dtrmid ad a dsity o stats. W also d a way o dtrmiig i a stat is illd or ot at a giv tmpratur. h distributio o th rgis o a larg umbr o particls ad

More information

Lyapunov Stability Analysis for Feedback Control Design

Lyapunov Stability Analysis for Feedback Control Design Copyright F.L. Lewis 008 All rights reserved Updated: uesday, November, 008 Lyapuov Stability Aalysis for Feedbac Cotrol Desig Lyapuov heorems Lyapuov Aalysis allows oe to aalyze the stability of cotiuous-time

More information

KISS: A Bit Too Simple. Greg Rose

KISS: A Bit Too Simple. Greg Rose KI: A Bit Too impl Grg Ros ggr@qualcomm.com Outli KI radom umbr grator ubgrators Efficit attack N KI ad attack oclusio PAGE 2 O approach to PRNG scurity "A radom umbr grator is lik sx: Wh it's good, its

More information

Ideal crystal : Regulary ordered point masses connected via harmonic springs

Ideal crystal : Regulary ordered point masses connected via harmonic springs Statistical thrmodyamics of crystals Mooatomic crystal Idal crystal : Rgulary ordrd poit masss coctd via harmoic sprigs Itratomic itractios Rprstd by th lattic forc-costat quivalt atom positios miima o

More information

Partial Derivatives: Suppose that z = f(x, y) is a function of two variables.

Partial Derivatives: Suppose that z = f(x, y) is a function of two variables. Chaptr Functions o Two Variabls Applid Calculus 61 Sction : Calculus o Functions o Two Variabls Now that ou hav som amiliarit with unctions o two variabls it s tim to start appling calculus to hlp us solv

More information

Self-Adjointness and Its Relationship to Quantum Mechanics. Ronald I. Frank 2016

Self-Adjointness and Its Relationship to Quantum Mechanics. Ronald I. Frank 2016 Ronald I. Frank 06 Adjoint https://n.wikipdia.org/wiki/adjoint In gnral thr is an oprator and a procss that dfin its adjoint *. It is thn slf-adjoint if *. Innr product spac https://n.wikipdia.org/wiki/innr_product_spac

More information

Physics 302 Exam Find the curve that passes through endpoints (0,0) and (1,1) and minimizes 1

Physics 302 Exam Find the curve that passes through endpoints (0,0) and (1,1) and minimizes 1 Physis Exam 6. Fid th urv that passs through dpoits (, ad (, ad miimizs J [ y' y ]dx Solutio: Si th itgrad f dos ot dpd upo th variabl of itgratio x, w will us th sod form of Eulr s quatio: f f y' y' y

More information

2.29 Numerical Fluid Mechanics Spring 2015 Lecture 12

2.29 Numerical Fluid Mechanics Spring 2015 Lecture 12 REVIEW Lctur 11: Numrical Fluid Mchaics Sprig 2015 Lctur 12 Fiit Diffrcs basd Polyomial approximatios Obtai polyomial (i gral u-qually spacd), th diffrtiat as dd Nwto s itrpolatig polyomial formulas Triagular

More information

[ ] Review. For a discrete-time periodic signal xn with period N, the Fourier series representation is

[ ] Review. For a discrete-time periodic signal xn with period N, the Fourier series representation is Discrt-tim ourir Trsform Rviw or discrt-tim priodic sigl x with priod, th ourir sris rprsttio is x + < > < > x, Rviw or discrt-tim LTI systm with priodic iput sigl, y H ( ) < > < > x H r rfrrd to s th

More information

NEW APPLICATIONS OF THE ABEL-LIOUVILLE FORMULA

NEW APPLICATIONS OF THE ABEL-LIOUVILLE FORMULA NE APPLICATIONS OF THE ABEL-LIOUVILLE FORMULA Mirca I CÎRNU Ph Dp o Mathmatics III Faculty o Applid Scincs Univrsity Polithnica o Bucharst Cirnumirca @yahoocom Abstract In a rcnt papr [] 5 th indinit intgrals

More information

COHORT MBA. Exponential function. MATH review (part2) by Lucian Mitroiu. The LOG and EXP functions. Properties: e e. lim.

COHORT MBA. Exponential function. MATH review (part2) by Lucian Mitroiu. The LOG and EXP functions. Properties: e e. lim. MTH rviw part b Lucian Mitroiu Th LOG and EXP functions Th ponntial function p : R, dfind as Proprtis: lim > lim p Eponntial function Y 8 6 - -8-6 - - X Th natural logarithm function ln in US- log: function

More information

On a problem of J. de Graaf connected with algebras of unbounded operators de Bruijn, N.G.

On a problem of J. de Graaf connected with algebras of unbounded operators de Bruijn, N.G. O a problm of J. d Graaf coctd with algbras of uboudd oprators d Bruij, N.G. Publishd: 01/01/1984 Documt Vrsio Publishr s PDF, also kow as Vrsio of Rcord (icluds fial pag, issu ad volum umbrs) Plas chck

More information

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

Bifurcation Theory. , a stationary point, depends on the value of α. At certain values Dnamic Macroconomic Thor Prof. Thomas Lux Bifurcation Thor Bifurcation: qualitativ chang in th natur of th solution occurs if a paramtr passs through a critical point bifurcation or branch valu. Local

More information

Quasi-Classical States of the Simple Harmonic Oscillator

Quasi-Classical States of the Simple Harmonic Oscillator Quasi-Classical Stats of th Simpl Harmonic Oscillator (Draft Vrsion) Introduction: Why Look for Eignstats of th Annihilation Oprator? Excpt for th ground stat, th corrspondnc btwn th quantum nrgy ignstats

More information

Cramér-Rao Inequality: Let f(x; θ) be a probability density function with continuous parameter

Cramér-Rao Inequality: Let f(x; θ) be a probability density function with continuous parameter WHEN THE CRAMÉR-RAO INEQUALITY PROVIDES NO INFORMATION STEVEN J. MILLER Abstract. W invstigat a on-paramtr family of probability dnsitis (rlatd to th Parto distribution, which dscribs many natural phnomna)

More information

Discrete Fourier Transform (DFT)

Discrete Fourier Transform (DFT) Discrt Fourir Trasorm DFT Major: All Egirig Majors Authors: Duc guy http://umricalmthods.g.us.du umrical Mthods or STEM udrgraduats 8/3/29 http://umricalmthods.g.us.du Discrt Fourir Trasorm Rcalld th xpotial

More information

y = 2xe x + x 2 e x at (0, 3). solution: Since y is implicitly related to x we have to use implicit differentiation: 3 6y = 0 y = 1 2 x ln(b) ln(b)

y = 2xe x + x 2 e x at (0, 3). solution: Since y is implicitly related to x we have to use implicit differentiation: 3 6y = 0 y = 1 2 x ln(b) ln(b) 4. y = y = + 5. Find th quation of th tangnt lin for th function y = ( + ) 3 whn = 0. solution: First not that whn = 0, y = (1 + 1) 3 = 8, so th lin gos through (0, 8) and thrfor its y-intrcpt is 8. y

More information

10. EXTENDING TRACTABILITY

10. EXTENDING TRACTABILITY Coping with NP-compltnss 0. EXTENDING TRACTABILITY ining small vrtx covrs solving NP-har problms on trs circular arc covrings vrtx covr in bipartit graphs Q. Suppos I n to solv an NP-complt problm. What

More information

H2 Mathematics Arithmetic & Geometric Series ( )

H2 Mathematics Arithmetic & Geometric Series ( ) H Mathmatics Arithmtic & Gomtric Sris (08 09) Basic Mastry Qustios Arithmtic Progrssio ad Sris. Th rth trm of a squc is 4r 7. (i) Stat th first four trms ad th 0th trm. (ii) Show that th squc is a arithmtic

More information

A Simple Proof that e is Irrational

A Simple Proof that e is Irrational Two of th most bautiful ad sigificat umbrs i mathmatics ar π ad. π (approximatly qual to 3.459) rprsts th ratio of th circumfrc of a circl to its diamtr. (approximatly qual to.788) is th bas of th atural

More information

DSP-First, 2/e. LECTURE # CH2-3 Complex Exponentials & Complex Numbers TLH MODIFIED. Aug , JH McClellan & RW Schafer

DSP-First, 2/e. LECTURE # CH2-3 Complex Exponentials & Complex Numbers TLH MODIFIED. Aug , JH McClellan & RW Schafer DSP-First, / TLH MODIFIED LECTURE # CH-3 Complx Exponntials & Complx Numbrs Aug 016 1 READING ASSIGNMENTS This Lctur: Chaptr, Scts. -3 to -5 Appndix A: Complx Numbrs Complx Exponntials Aug 016 LECTURE

More information

Chemical Physics II. More Stat. Thermo Kinetics Protein Folding...

Chemical Physics II. More Stat. Thermo Kinetics Protein Folding... Chmical Physics II Mor Stat. Thrmo Kintics Protin Folding... http://www.nmc.ctc.com/imags/projct/proj15thumb.jpg http://nuclarwaponarchiv.org/usa/tsts/ukgrabl2.jpg http://www.photolib.noaa.gov/corps/imags/big/corp1417.jpg

More information

Digital Signal Processing, Fall 2006

Digital Signal Processing, Fall 2006 Digital Sigal Procssig, Fall 6 Lctur 9: Th Discrt Fourir Trasfor Zhg-Hua Ta Dpartt of Elctroic Systs Aalborg Uivrsity, Dar zt@o.aau.d Digital Sigal Procssig, I, Zhg-Hua Ta, 6 Cours at a glac MM Discrt-ti

More information

Basic Polyhedral theory

Basic Polyhedral theory Basic Polyhdral thory Th st P = { A b} is calld a polyhdron. Lmma 1. Eithr th systm A = b, b 0, 0 has a solution or thr is a vctorπ such that π A 0, πb < 0 Thr cass, if solution in top row dos not ist

More information

ENGG 1203 Tutorial. Difference Equations. Find the Pole(s) Finding Equations and Poles

ENGG 1203 Tutorial. Difference Equations. Find the Pole(s) Finding Equations and Poles ENGG 03 Tutoial Systms ad Cotol 9 Apil Laig Obctivs Z tasfom Complx pols Fdbac cotol systms Ac: MIT OCW 60, 6003 Diffc Equatios Cosid th systm pstd by th followig diffc quatio y[ ] x[ ] (5y[ ] 3y[ ]) wh

More information

INTRODUCTION TO SAMPLING DISTRIBUTIONS

INTRODUCTION TO SAMPLING DISTRIBUTIONS http://wiki.stat.ucla.du/socr/id.php/socr_courss_2008_thomso_econ261 INTRODUCTION TO SAMPLING DISTRIBUTIONS By Grac Thomso INTRODUCTION TO SAMPLING DISTRIBUTIONS Itro to Samplig 2 I this chaptr w will

More information

Statistics 3858 : Likelihood Ratio for Exponential Distribution

Statistics 3858 : Likelihood Ratio for Exponential Distribution Statistics 3858 : Liklihood Ratio for Expotial Distributio I ths two xampl th rjctio rjctio rgio is of th form {x : 2 log (Λ(x)) > c} for a appropriat costat c. For a siz α tst, usig Thorm 9.5A w obtai

More information

Chapter (8) Estimation and Confedence Intervals Examples

Chapter (8) Estimation and Confedence Intervals Examples Chaptr (8) Estimatio ad Cofdc Itrvals Exampls Typs of stimatio: i. Poit stimatio: Exampl (1): Cosidr th sampl obsrvatios, 17,3,5,1,18,6,16,10 8 X i i1 17 3 5 118 6 16 10 116 X 14.5 8 8 8 14.5 is a poit

More information

DTFT Properties. Example - Determine the DTFT Y ( e ) of n. Let. We can therefore write. From Table 3.1, the DTFT of x[n] is given by 1

DTFT Properties. Example - Determine the DTFT Y ( e ) of n. Let. We can therefore write. From Table 3.1, the DTFT of x[n] is given by 1 DTFT Proprtis Exampl - Dtrmi th DTFT Y of y α µ, α < Lt x α µ, α < W ca thrfor writ y x x From Tabl 3., th DTFT of x is giv by ω X ω α ω Copyright, S. K. Mitra Copyright, S. K. Mitra DTFT Proprtis DTFT

More information

Slide 1. Slide 2. Slide 3 DIGITAL SIGNAL PROCESSING CLASSIFICATION OF SIGNALS

Slide 1. Slide 2. Slide 3 DIGITAL SIGNAL PROCESSING CLASSIFICATION OF SIGNALS Slid DIGITAL SIGAL PROCESSIG UIT I DISCRETE TIME SIGALS AD SYSTEM Slid Rviw of discrt-tim signals & systms Signal:- A signal is dfind as any physical quantity that varis with tim, spac or any othr indpndnt

More information

Fooling Newton s Method a) Find a formula for the Newton sequence, and verify that it converges to a nonzero of f. A Stirling-like Inequality

Fooling Newton s Method a) Find a formula for the Newton sequence, and verify that it converges to a nonzero of f. A Stirling-like Inequality Foolig Nwto s Mthod a Fid a formla for th Nwto sqc, ad vrify that it covrgs to a ozro of f. ( si si + cos 4 4 3 4 8 8 bt f. b Fid a formla for f ( ad dtrmi its bhavior as. f ( cos si + as A Stirlig-li

More information

Ordinary Differential Equations

Ordinary Differential Equations Basi Nomlatur MAE 0 all 005 Egirig Aalsis Ltur Nots o: Ordiar Diffrtial Equatios Author: Profssor Albrt Y. Tog Tpist: Sakurako Takahashi Cosidr a gral O. D. E. with t as th idpdt variabl, ad th dpdt variabl.

More information

Fourier Transforms and the Wave Equation. Key Mathematics: More Fourier transform theory, especially as applied to solving the wave equation.

Fourier Transforms and the Wave Equation. Key Mathematics: More Fourier transform theory, especially as applied to solving the wave equation. Lur 7 Fourir Transforms and th Wav Euation Ovrviw and Motivation: W first discuss a fw faturs of th Fourir transform (FT), and thn w solv th initial-valu problm for th wav uation using th Fourir transform

More information

Schematic of a mixed flow reactor (both advection and dispersion must be accounted for)

Schematic of a mixed flow reactor (both advection and dispersion must be accounted for) Cas stuy 6.1, R: Chapra an Canal, p. 769. Th quation scribin th concntration o any tracr in an lonat ractor is known as th avction-isprsion quation an may b writtn as: Schmatic o a mi low ractor (both

More information

Castiel, Supernatural, Season 6, Episode 18

Castiel, Supernatural, Season 6, Episode 18 13 Differetial Equatios the aswer to your questio ca best be epressed as a series of partial differetial equatios... Castiel, Superatural, Seaso 6, Episode 18 A differetial equatio is a mathematical equatio

More information

CALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems

CALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems CDS 101 1. For each of the following linear systems, determine whether the origin is asymptotically stable and, if so, plot the step response and frequency response for the system. If there are multiple

More information

Molecular Orbitals in Inorganic Chemistry

Molecular Orbitals in Inorganic Chemistry Outlin olcular Orbitals in Inorganic Chmistry Dr. P. Hunt p.hunt@imprial.ac.uk Rm 167 (Chmistry) http://www.ch.ic.ac.uk/hunt/ octahdral complxs forming th O diagram for Oh colour, slction ruls Δoct, spctrochmical

More information

Folding of Hyperbolic Manifolds

Folding of Hyperbolic Manifolds It. J. Cotmp. Math. Scics, Vol. 7, 0, o. 6, 79-799 Foldig of Hyprbolic Maifolds H. I. Attiya Basic Scic Dpartmt, Collg of Idustrial Educatio BANE - SUEF Uivrsity, Egypt hala_attiya005@yahoo.com Abstract

More information

Section 11.6: Directional Derivatives and the Gradient Vector

Section 11.6: Directional Derivatives and the Gradient Vector Sction.6: Dirctional Drivativs and th Gradint Vctor Practic HW rom Stwart Ttbook not to hand in p. 778 # -4 p. 799 # 4-5 7 9 9 35 37 odd Th Dirctional Drivativ Rcall that a b Slop o th tangnt lin to th

More information

Where k is either given or determined from the data and c is an arbitrary constant.

Where k is either given or determined from the data and c is an arbitrary constant. Exponntial growth and dcay applications W wish to solv an quation that has a drivativ. dy ky k > dx This quation says that th rat of chang of th function is proportional to th function. Th solution is

More information

Introduction to Arithmetic Geometry Fall 2013 Lecture #20 11/14/2013

Introduction to Arithmetic Geometry Fall 2013 Lecture #20 11/14/2013 18.782 Introduction to Arithmtic Gomtry Fall 2013 Lctur #20 11/14/2013 20.1 Dgr thorm for morphisms of curvs Lt us rstat th thorm givn at th nd of th last lctur, which w will now prov. Thorm 20.1. Lt φ:

More information

[ ] 1+ lim G( s) 1+ s + s G s s G s Kacc SYSTEM PERFORMANCE. Since. Lecture 10: Steady-state Errors. Steady-state Errors. Then

[ ] 1+ lim G( s) 1+ s + s G s s G s Kacc SYSTEM PERFORMANCE. Since. Lecture 10: Steady-state Errors. Steady-state Errors. Then SYSTEM PERFORMANCE Lctur 0: Stady-tat Error Stady-tat Error Lctur 0: Stady-tat Error Dr.alyana Vluvolu Stady-tat rror can b found by applying th final valu thorm and i givn by lim ( t) lim E ( ) t 0 providd

More information

Notes on iteration and Newton s method. Iteration

Notes on iteration and Newton s method. Iteration Notes o iteratio ad Newto s method Iteratio Iteratio meas doig somethig over ad over. I our cotet, a iteratio is a sequece of umbers, vectors, fuctios, etc. geerated by a iteratio rule of the type 1 f

More information

Iterative Methods of Order Four for Solving Nonlinear Equations

Iterative Methods of Order Four for Solving Nonlinear Equations Itrativ Mods of Ordr Four for Solvig Noliar Equatios V.B. Kumar,Vatti, Shouri Domii ad Mouia,V Dpartmt of Egirig Mamatis, Formr Studt of Chmial Egirig Adhra Uivrsity Collg of Egirig A, Adhra Uivrsity Visakhapatam

More information

EEO 401 Digital Signal Processing Prof. Mark Fowler

EEO 401 Digital Signal Processing Prof. Mark Fowler EEO 401 Digital Signal Procssing Prof. Mark Fowlr Dtails of th ot St #19 Rading Assignmnt: Sct. 7.1.2, 7.1.3, & 7.2 of Proakis & Manolakis Dfinition of th So Givn signal data points x[n] for n = 0,, -1

More information

Chapter 3 Fourier Series Representation of Periodic Signals

Chapter 3 Fourier Series Representation of Periodic Signals Chptr Fourir Sris Rprsttio of Priodic Sigls If ritrry sigl x(t or x[] is xprssd s lir comitio of som sic sigls th rspos of LI systm coms th sum of th idividul rsposs of thos sic sigls Such sic sigl must:

More information

On the approximation of the constant of Napier

On the approximation of the constant of Napier Stud. Uiv. Babş-Bolyai Math. 560, No., 609 64 O th approximatio of th costat of Napir Adri Vrscu Abstract. Startig from som oldr idas of [] ad [6], w show w facts cocrig th approximatio of th costat of

More information

Chapter Five. More Dimensions. is simply the set of all ordered n-tuples of real numbers x = ( x 1

Chapter Five. More Dimensions. is simply the set of all ordered n-tuples of real numbers x = ( x 1 Chatr Fiv Mor Dimsios 51 Th Sac R W ar ow rard to mov o to sacs of dimsio gratr tha thr Ths sacs ar a straightforward gralizatio of our Euclida sac of thr dimsios Lt b a ositiv itgr Th -dimsioal Euclida

More information

TuLiP: A Software Toolbox for Receding Horizon Temporal Logic Planning & Computer Lab 2

TuLiP: A Software Toolbox for Receding Horizon Temporal Logic Planning & Computer Lab 2 TuLiP: A Softwar Toolbox for Rcding Horizon Tmporal Logic Planning & Computr Lab 2 Nok Wongpiromsarn Richard M. Murray Ufuk Topcu EECI, 21 March 2013 Outlin Ky Faturs of TuLiP Embddd control softwar synthsis

More information

Taylor and Maclaurin Series

Taylor and Maclaurin Series Taylor ad Maclauri Sris Taylor ad Maclauri Sris Thory sctio which dals with th followig topics: - Th Sigma otatio for summatio. - Dfiitio of Taylor sris. - Commo Maclauri sris. - Taylor sris ad Itrval

More information

ELG3150 Assignment 3

ELG3150 Assignment 3 ELG350 Aigmt 3 Aigmt 3: E5.7; P5.6; P5.6; P5.9; AP5.; DP5.4 E5.7 A cotrol ytm for poitioig th had of a floppy dik driv ha th clodloop trafr fuctio 0.33( + 0.8) T ( ) ( + 0.6)( + 4 + 5) Plot th pol ad zro

More information

Dynamic Response of Linear Systems

Dynamic Response of Linear Systems Dyamic Respose of Liear Systems Liear System Respose Superpositio Priciple Resposes to Specific Iputs Dyamic Respose of st Order Systems Characteristic Equatio - Free Respose Stable st Order System Respose

More information

How much air is required by the people in this lecture theatre during this lecture?

How much air is required by the people in this lecture theatre during this lecture? 3 NTEGRATON tgrtio is us to swr qustios rltig to Ar Volum Totl qutity such s: Wht is th wig r of Boig 747? How much will this yr projct cost? How much wtr os this rsrvoir hol? How much ir is rquir y th

More information

Mathematics. Complex Number rectangular form. Quadratic equation. Quadratic equation. Complex number Functions: sinusoids. Differentiation Integration

Mathematics. Complex Number rectangular form. Quadratic equation. Quadratic equation. Complex number Functions: sinusoids. Differentiation Integration Mathmatics Compl numbr Functions: sinusoids Sin function, cosin function Diffrntiation Intgration Quadratic quation Quadratic quations: a b c 0 Solution: b b 4ac a Eampl: 1 0 a= b=- c=1 4 1 1or 1 1 Quadratic

More information

3.4 Properties of the Stress Tensor

3.4 Properties of the Stress Tensor cto.4.4 Proprts of th trss sor.4. trss rasformato Lt th compots of th Cauchy strss tsor a coordat systm wth bas vctors b. h compots a scod coordat systm wth bas vctors j,, ar gv by th tsor trasformato

More information