STATE-SPACE MODELLING. A mass balance across the tank gives:
|
|
- Audrey Turner
- 6 years ago
- Views:
Transcription
1 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 use of sae-space modelling in conrol books and applicaions. Alhough sae-space is ofen considered o be a complex represenaion, i is in he end simply a mehod for describing he behaviour of dynamic sysems. The growing populariy of sae-space modelling can be aribued o he fac ha i has a very general form and is able o represen linear, nonlinear, ime-invarian and ime-varying dynamics in a relaively compac form. Furhermore, in many applicaion areas he sae-space approach is very inuiive, allowing rapid developmen of dynamic models from firs principles. This aricle will illusrae he ideas ha underpin sae-space analysis and illusrae hem wih a sraighforward example. The Basics, and an Example: In he process indusries, conrol heory has ypically been based on dynamic models ha are represened as Ordinary Differenial Equaions (ODEs) and/or ransfer funcions. As an example, consider he sysem illusraed in Figure. This sysem is used in many ex-books and concerns he relaionship beween he flow of liquid ino he ank, qin () and he level of liquid, h. From his model a feedback conrol sysem can be designed o mainain he level of liquid in he ank. q in h cross secional area = q = kh( ) ou Figure : Liquid Level Sysem The flow of liquid ou of he ank is proporional o level of liquid h if he flow is laminar: qou () = kh() where k is a consan parameer. A mass balance across he ank gives: dv() qin() qou () = i.e. he rae of change of liquid volume, V, is equal o he volumeric flow rae of liquid in minus he flow rae of liquid ou. The mass balance can be re-arranged o give: d h() qin() qou () = and since A T is consan: Adh T () qin() qou () =. Replacing qou () by kh() gives: Adh T () qin() kh() = and more re-arrangemen gives: dh() k + h () = q in() () Equaion () describes a linear, firs order relaionship beween he flow ino he ank and he liquid level. If zero iniial condiions are assumed, hen his equaion can be ransformed ino he Laplace domain as follows: k sh( s) + H( s) = Qin( s) where H( s ) and Qin ( s ) are he Laplace ransforms of he ime rends h () and qin ( ). For insance, if qin () is a sep change of magniude q, hen Qin ( s) = q s. Afer some algebra, he ransfer funcion for he ank is: Hs ( ) k = = () Qin( s) s+ k ( k) s+ The generic form of a firs order ransfer funcion K p is where K p is he sysem gain and τ is sτ + he ime consan, hence by comparison, he gain of he ank sysem is k and is ime consan is k. So far, he analysis has involved an ODE wih one dependen variable h(), and one independen variable qin (). Applicaion of he Laplace ransform o he ODE gave a singleinpu-single-oupu ransfer funcion showing he gain and ime consan. Nex, in anicipaion of more complicaed sysems wih more han one dependen variable, le s examine a more general formulaion of he same equaion. Equaion () is re-considered and a new
2 B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer variable, x () is now defined, such ha x () = h (). Equaion () can be re-wrien as: () k = x () + q in() which in a more general form can be expressed a dynamic sae equaion: = x() = Ax() + Bu () (3) where: k x () = x(), u () = qin(), A = and B = A T Equaion (3) is he general, linear sae equaion, where he vecor x represens he sae variables of he process. I anicipaes more complex sysems by use of vecors and marices. In he ank example, however, here is a single sae variable h() and a single inpu variable qin ( ) so he A and B marices are scalars. Sae variables are he smalles subse of sysem variables ha describe he enire dynamic characerisics of he sysem. They can be hough of as inernal elemens of he sysem ha are relaed o, or in some case are acually equal o, he oupu variables. Alhough i can someimes be beneficial, i is no necessary for he sae variables o be measured or even have physical meaning. In he even more general non-linear case he sae equaion is defined as: () = f ( xu,, ) where f is a non-linear funcion. q in h h ( ) q = k h h Figure : A Two Tank Sysem q = k h 3 3 A Two Sae Sysem: The coupled wo-ank sysem shown in Figure is an example of a sysem wih more han one sae. The coupled wo ank sysem has wo physical saes, h and h and, if he anks each have uni cross secional area (i.e. A T = ), he sae equaion is: dh k k h () = + qin dh k ( k k + 3) h () Expanding ou he marices and rearranging shows ha he sae equaion is nohing more han a volumeric balance for each ank. dh = k( h h) + qin dh = k( h h) k3h Thus, here is no new physics in a sae equaion. The reason why conrol engineers use i is ha is a convenien way of sudying he mahemaical properies of he physical equaions. Using sae-space leads o some very powerful conroller designs ha would be oo cumbersome o consider wihou he compac marix formulaion. Referring back o he general sae equaion: = x() = Ax() + Bu () In he wo-ank example, he sae vecor x () is h, he inpu u () is q h in and he consancoefficien A and B marices are: k k A = and B = k ( k + k 3) Observer Equaion: The sae equaion (3) describes how he saes vary wih ime. In he case of he examples above, he oupu variables were he liquid levels, which were also defined as he sae variables. However, in many examples, he sae variables will be differen from he oupu variables. I is herefore necessary o have anoher equaion which describes he relaionship beween he process oupus and sae variables. For general sysems, his equaion is defined as follows: y( ) = g( x, u, ) where g is a non-linear funcion. For a linear, ime-invarian sysem his equaion always akes he form: y( ) = Cx( ) + Du ( ) (4) where y() is a vecor of oupu variables. In he wo-ank example: C = and D =
3 B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer The fac C is a uni marix shows he oupu variables are equal o he sae variables in he wo ank example. The above equaion for y () is referred o as he observer equaion. The combinaion of he sae equaion (3) and observer equaion (4) is called a sae-space model. Anoher sae-space form: For any dynamic sysem, here is no unique sae-space model. In fac here is an infinie number of sae-space equaions ha can be developed, depending on how he saes are defined. As an example, he physical balance equaions for he wo ank sysem naurally yielded wo simulaneous firs order differenial equaions. Bu anoher way o rearrange he equaions is o eliminae eiher h or h o give a single second order differenial equaion. For insance, h may be eliminaed o give an equaion for h as follows: Adding he wo mass balances: dh = k( h h) + qin dh = k( h h) k3h gives dh dh + = qin k3h and, by rearrangemen: dh dh = qin k3h The second mass balance can be differeniaed o give: and now = k k3 d h dh dh dh dh can be replaced o give: = k qin k3h k3 d h dh dh dh and finally: dh ( k3 k) kkh 3 kq in d h = (5) Again, here is no new physics, bu now he focus is on h which obeys a second order differenial equaion. Now, he saes can be represened no as h and h, bu as h and dh. Therefore he sae vecor would be: h () x () = = () () x dh x Clearly (by definiion) one of he sae equaions is = x () and he second comes from he rearrangemen of equaion (5): dh ( k3 k) k k3h kqin d h = + + (6) dh Using h( ) = x( ) and = x () gives: = ( k3 + k) x() kk3x() + kqin Hence he sae equaions are: = x () = kkx 3 () ( k3 + k) x() + kq in and hey can be expressed in marix form as: x () k = + q kk 3 ( k3 k + ) x () This way of describing a sae-space is referred o as he conrollable canonical form. There are however many oher ways in which he A, B, C and D marices can be re-arranged o produce, for example, he observable canonical form. Relaionship Beween Sae-Space, ODEs and Transfer Funcions: Process conrol engineers end o be familiar wih ODEs and he ransfer funcion derived from he Laplace ransform of he ODE. As shown in he previous examples sae-space does no describe anyhing new, i is simply a re-arrangemen of he dynamic equaions ino a form ha is more compac and convenien for conrol sysem developmen and analysis. Therefore i is useful o see how he coefficiens in he ODE relae o he parameers of he ransfer funcion and he elemens in he sae-space A, B, C and D marices. The ODE in he general case is a single inpu, single oupu, high order dynamic sysem described by he following differenial equaion: n n () () d y d y + a n n + n m d u() + ay () = bm + m m d u b () m m () + + b u in
4 B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer h h qin The ransfer funcion for his sysem, derived by applying Laplace ransforms wih zero iniial condiion is: m m Y( s) s + bm s bs + b Gs ( ) = = U n n ( s ) s + an s as + a and he equivalen sae-space equaion in conrollable canonical form is as follows: A =, B = a a a a n C = [ b b bm ] and D = As can be seen, he same coefficiens appear in all hree forms indicaing ha hey are hree differen ways of looking a he same equaions. Relaionship wih sep-response models: The dynamics for elecrical and mechanical sysems are ofen reasonably well defined and can be derived from a physical undersanding of he sysem. In conras, process sysems may no be undersood wih sufficien accuracy o develop a model from firs principles and i is ypical for sep or oher ess o be performed o beer undersand he dynamics. This secion will demonsrae ha he resuls of sep response ess can be used o provide boh ransfer funcion models and a sae-space model Time (secs) Figure 3. Sep responses of he wo-ank sysem. Figure 3 shows sep ess for a wo-ank sysem in which he inle flow rae was doubled a ime sec. The levels h and h boh changed. Level h shows a response wih a seep iniial gradien ha is characerisic of a firs order sysem. Figure 3 indicaes ha he gain is because he uni sep in q in produces a change in level of on he h scale. The ime consan is approximaely seconds, and hence: H ( s ) in ( ) Q s s+ This ransfer funcion is approximae, because alhough he sep response looks firs order here is in pracice is a back-pressure on he firs ank from he second ank. The dynamics are acually of a higher order, however he effec is small and was no deeced by he sep es. The level in he second ank, h, has he S- shaped response ha is characerisic of a second order sysem. I requires a numerical mehod o fi a ransfer funcion model o such sep response daa, for insance by using he ARX command in Malab. The resuls from sysem idenificaion (o one s.f.) were: H( s). in ( ) = Q s s +.7s+.5 If boh levels are now considered, he complee ransfer funcion model for his sysem is: H( s) H( s). Qin( s) s + s +.7s+.5 An equivalen sae-space equaion is as follows: x() = Ax() + Bu() y Cx Du () = () + () where: h() h() x() = h (), y() =, () = qin() h () u dh () and / A = B =.5.7. C =, = D The ransfer funcions and sae-space models have he same inpu-oupu behaviour, again illusraing ha he sae-space form does no include any new physical informaion. However, as is beginning o become apparen from his example, wih more inpu and oupu variables he ransfer funcion model quickly becomes unwieldy. In conras, he complexiy of he saespace represenaion does no really change. The sae-space model always consiss of wo compac equaions, one is he dynamic sae equaion and he oher is he algebraic observer
5 B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer equaion. The only change as he sysem becomes more complex is ha he dimensions of he marices increase. There are hree saes in he above model, which effecively combine he saes h () and h () from he model based on volumeric balancing wih he saes h () and dh ha were presen in he conrollable canonical form. I can ake some judgemen o selec he bes saes for a given purpose, however he abiliy o make good choices grows wih experience. There are also some useful algorihms such as sub-space idenificaion which can help.. Commens and Conclusions: The aricle has shown how o explore and invesigae he physical equaions of a sysem by means of is inernal saes. The saes may have a direc physical meaning, or hey may be absrac ones such as when he conrollable canonical form is used. There are many reasons for making he saes of a sysem explici in he mahemaical formulaion. For insance, if he saes can be esimaed from process measuremens, hen here is an opporuniy o use sae feedback. Rae feedback, where he rae of change of he conrolled variable is used as a feedback signal, is an example of sae feedback. Here is a shor lis of some of he oher hings ha working wih sae-space makes possible:. Provides a model srucure from which feedback conrol sysems can be designed wih relaive ease.. Enables closed-loop characerisics, such as robusness and sabiliy, o be analysed and considered during conrol design. The complexiy of he mahemaics involved for his is such ha i would no be pracical o apply i o oher modelling formas. 3. In siuaions where he oupu variables are no direcly measured, he observer equaion can be used o esimae hese measuremens from he sae variables. Despie heir advanages, here is a very serious weakness o sae-space models. When using ransfer models, paricularly firs and second order, i is sraigh forward o visualize he sepresponse of he sysem. Unforunaely, his is no he case wih sae-space models which are more complicaed o inerpre and ypically require he use of a compuer o simulae sep responses and evaluae he sysem ime consans. On he oher hand, i is quie easy o obain a sae-space model from a sep-response, as shown in his aricle. The example considered was sraighforward, however more realisic cases require care and judgemen in selecion of he saes, and opimal sae-space models may no resul. Mehods for obaining good saespace models for complex sysems using pseudo random binary inpu sequences or sepresponse daa will be considered in a fuure aricle on sub-space idenificaion. Furher reading. Those who are ineresed o dig deeper migh like o visi he URL below: hp://en.wikipedia.org/wiki/sae_space_(conrols)
10. 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 informationModal identification of structures from roving input data by means of maximum likelihood estimation of the state space model
Modal idenificaion of srucures from roving inpu daa by means of maximum likelihood esimaion of he sae space model J. Cara, J. Juan, E. Alarcón Absrac The usual way o perform a forced vibraion es is o fix
More informationDesigning Information Devices and Systems I Spring 2019 Lecture Notes Note 17
EES 16A Designing Informaion Devices and Sysems I Spring 019 Lecure Noes Noe 17 17.1 apaciive ouchscreen In he las noe, we saw ha a capacior consiss of wo pieces on conducive maerial separaed by a nonconducive
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 informationSimulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010
Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid
More informationLecture 1 Overview. course mechanics. outline & topics. what is a linear dynamical system? why study linear systems? some examples
EE263 Auumn 27-8 Sephen Boyd Lecure 1 Overview course mechanics ouline & opics wha is a linear dynamical sysem? why sudy linear sysems? some examples 1 1 Course mechanics all class info, lecures, homeworks,
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 informationLecture 20: Riccati Equations and Least Squares Feedback Control
34-5 LINEAR SYSTEMS Lecure : Riccai Equaions and Leas Squares Feedback Conrol 5.6.4 Sae Feedback via Riccai Equaions A recursive approach in generaing he marix-valued funcion W ( ) equaion for i for he
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 informationR t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t
Exercise 7 C P = α + β R P + u C = αp + βr + v (a) (b) C R = α P R + β + w (c) Assumpions abou he disurbances u, v, w : Classical assumions on he disurbance of one of he equaions, eg. on (b): E(v v s P,
More informationCHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK
175 CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK 10.1 INTRODUCTION Amongs he research work performed, he bes resuls of experimenal work are validaed wih Arificial Neural Nework. From he
More informationNotes on Kalman Filtering
Noes on Kalman Filering Brian Borchers and Rick Aser November 7, Inroducion Daa Assimilaion is he problem of merging model predicions wih acual measuremens of a sysem o produce an opimal esimae of he curren
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 informationA DELAY-DEPENDENT STABILITY CRITERIA FOR T-S FUZZY SYSTEM WITH TIME-DELAYS
A DELAY-DEPENDENT STABILITY CRITERIA FOR T-S FUZZY SYSTEM WITH TIME-DELAYS Xinping Guan ;1 Fenglei Li Cailian Chen Insiue of Elecrical Engineering, Yanshan Universiy, Qinhuangdao, 066004, China. Deparmen
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 informationChapter 7 Response of First-order RL and RC Circuits
Chaper 7 Response of Firs-order RL and RC Circuis 7.- The Naural Response of RL and RC Circuis 7.3 The Sep Response of RL and RC Circuis 7.4 A General Soluion for Sep and Naural Responses 7.5 Sequenial
More informationDesign of a control system
SE3 Prof. Davide Manca Poliecnico di Milano Dynamics and Conrol of Chemical Processes Soluion o Lab #3 Design of a conrol sysem Davide Manca Dynamics and Conrol of Chemical Processes Maser Degree in ChemEng
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 informationINVERSE RESPONSE COMPENSATION BY ESTIMATING PARAMETERS OF A PROCESS COMPRISING OF TWO FIRST ORDER SYSTEMS
Inernaional Journal of Informaion Technology and nowledge Managemen July-December 0, Volume 5, No., pp. 433-438 INVERSE RESPONSE COMPENSATION BY ESTIMATING PARAMETERS OF A PROCESS COMPRISING OF TWO FIRST
More informationnot to be republished NCERT MATHEMATICAL MODELLING Appendix 2 A.2.1 Introduction A.2.2 Why Mathematical Modelling?
256 MATHEMATICS A.2.1 Inroducion In class XI, we have learn abou mahemaical modelling as an aemp o sudy some par (or form) of some real-life problems in mahemaical erms, i.e., he conversion of a physical
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 informationChapter 1 Fundamental Concepts
Chaper 1 Fundamenal Conceps 1 Signals A signal is a paern of variaion of a physical quaniy, ofen as a funcion of ime (bu also space, disance, posiion, ec). These quaniies are usually he independen variables
More informationEcon107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)
I. Definiions and Problems A. Perfec Mulicollineariy Econ7 Applied Economerics Topic 7: Mulicollineariy (Sudenmund, Chaper 8) Definiion: Perfec mulicollineariy exiss in a following K-variable regression
More informationClass Meeting # 10: Introduction to the Wave Equation
MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion
More informationAN OPTIMAL CONTROL PROBLEM FOR SPLINES ASSOCIATED TO LINEAR DIFFERENTIAL OPERATORS
CONTROLO 6 7h Poruguese Conference on Auomaic Conrol Insiuo Superior Técnico, Lisboa, Porugal Sepember -3, 6 AN OPTIMAL CONTROL PROBLEM FOR SPLINES ASSOCIATED TO LINEAR DIFFERENTIAL OPERATORS Rui C. Rodrigues,
More informationBasic Circuit Elements Professor J R Lucas November 2001
Basic Circui Elemens - J ucas An elecrical circui is an inerconnecion of circui elemens. These circui elemens can be caegorised ino wo ypes, namely acive and passive elemens. Some Definiions/explanaions
More informationElectrical and current self-induction
Elecrical and curren self-inducion F. F. Mende hp://fmnauka.narod.ru/works.hml mende_fedor@mail.ru Absrac The aricle considers he self-inducance of reacive elemens. Elecrical self-inducion To he laws of
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 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 informationMechanical Fatigue and Load-Induced Aging of Loudspeaker Suspension. Wolfgang Klippel,
Mechanical Faigue and Load-Induced Aging of Loudspeaker Suspension Wolfgang Klippel, Insiue of Acousics and Speech Communicaion Dresden Universiy of Technology presened a he ALMA Symposium 2012, Las Vegas
More informationA Specification Test for Linear Dynamic Stochastic General Equilibrium Models
Journal of Saisical and Economeric Mehods, vol.1, no.2, 2012, 65-70 ISSN: 2241-0384 (prin), 2241-0376 (online) Scienpress Ld, 2012 A Specificaion Tes for Linear Dynamic Sochasic General Equilibrium Models
More informationSome Basic Information about M-S-D Systems
Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,
More informationINDEX. Transient analysis 1 Initial Conditions 1
INDEX Secion Page Transien analysis 1 Iniial Condiions 1 Please inform me of your opinion of he relaive emphasis of he review maerial by simply making commens on his page and sending i o me a: Frank Mera
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 informationBias in Conditional and Unconditional Fixed Effects Logit Estimation: a Correction * Tom Coupé
Bias in Condiional and Uncondiional Fixed Effecs Logi Esimaion: a Correcion * Tom Coupé Economics Educaion and Research Consorium, Naional Universiy of Kyiv Mohyla Academy Address: Vul Voloska 10, 04070
More informationEECE 301 Signals & Systems Prof. Mark Fowler
EECE 31 Signals & Sysems Prof. Mark Fowler Noe Se #1 C-T Sysems: Convoluion Represenaion Reading Assignmen: Secion 2.6 of Kamen and Heck 1/11 Course Flow Diagram The arrows here show concepual flow beween
More informationAn introduction to the theory of SDDP algorithm
An inroducion o he heory of SDDP algorihm V. Leclère (ENPC) Augus 1, 2014 V. Leclère Inroducion o SDDP Augus 1, 2014 1 / 21 Inroducion Large scale sochasic problem are hard o solve. Two ways of aacking
More informationThe expectation value of the field operator.
The expecaion value of he field operaor. Dan Solomon Universiy of Illinois Chicago, IL dsolom@uic.edu June, 04 Absrac. Much of he mahemaical developmen of quanum field heory has been in suppor of deermining
More informationNavneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi
Creep in Viscoelasic Subsances Numerical mehods o calculae he coefficiens of he Prony equaion using creep es daa and Herediary Inegrals Mehod Navnee Saini, Mayank Goyal, Vishal Bansal (23); Term Projec
More informationVehicle Arrival Models : Headway
Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where
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 informationOrdinary Differential Equations
Lecure 22 Ordinary Differenial Equaions Course Coordinaor: Dr. Suresh A. Karha, Associae Professor, Deparmen of Civil Engineering, IIT Guwahai. In naure, mos of he phenomena ha can be mahemaically described
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 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 informationSection 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients
Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous
More informationSliding Mode Controller for Unstable Systems
S. SIVARAMAKRISHNAN e al., Sliding Mode Conroller for Unsable Sysems, Chem. Biochem. Eng. Q. 22 (1) 41 47 (28) 41 Sliding Mode Conroller for Unsable Sysems S. Sivaramakrishnan, A. K. Tangirala, and M.
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION DOI: 0.038/NCLIMATE893 Temporal resoluion and DICE * Supplemenal Informaion Alex L. Maren and Sephen C. Newbold Naional Cener for Environmenal Economics, US Environmenal Proecion
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 informationAnalytical Solutions of an Economic Model by the Homotopy Analysis Method
Applied Mahemaical Sciences, Vol., 26, no. 5, 2483-249 HIKARI Ld, www.m-hikari.com hp://dx.doi.org/.2988/ams.26.6688 Analyical Soluions of an Economic Model by he Homoopy Analysis Mehod Jorge Duare ISEL-Engineering
More information0.1 MAXIMUM LIKELIHOOD ESTIMATION EXPLAINED
0.1 MAXIMUM LIKELIHOOD ESTIMATIO EXPLAIED Maximum likelihood esimaion is a bes-fi saisical mehod for he esimaion of he values of he parameers of a sysem, based on a se of observaions of a random variable
More informationConcourse Math Spring 2012 Worked Examples: Matrix Methods for Solving Systems of 1st Order Linear Differential Equations
Concourse Mah 80 Spring 0 Worked Examples: Marix Mehods for Solving Sysems of s Order Linear Differenial Equaions The Main Idea: Given a sysem of s order linear differenial equaions d x d Ax wih iniial
More informationSub Module 2.6. Measurement of transient temperature
Mechanical Measuremens Prof. S.P.Venkaeshan Sub Module 2.6 Measuremen of ransien emperaure Many processes of engineering relevance involve variaions wih respec o ime. The sysem properies like emperaure,
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 informationMath Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.
Mah 2250-004 Week 4 April 6-20 secions 7.-7.3 firs order sysems of linear differenial equaions; 7.4 mass-spring sysems. Mon Apr 6 7.-7.2 Sysems of differenial equaions (7.), and he vecor Calculus we need
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 informationGeorey E. Hinton. University oftoronto. Technical Report CRG-TR February 22, Abstract
Parameer Esimaion for Linear Dynamical Sysems Zoubin Ghahramani Georey E. Hinon Deparmen of Compuer Science Universiy oftorono 6 King's College Road Torono, Canada M5S A4 Email: zoubin@cs.orono.edu Technical
More informationRC, RL and RLC circuits
Name Dae Time o Complee h m Parner Course/ Secion / Grade RC, RL and RLC circuis Inroducion In his experimen we will invesigae he behavior of circuis conaining combinaions of resisors, capaciors, and inducors.
More information2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes
Some common engineering funcions 2.7 Inroducion This secion provides a caalogue of some common funcions ofen used in Science and Engineering. These include polynomials, raional funcions, he modulus funcion
More informationReading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4.
PHY1 Elecriciy Topic 7 (Lecures 1 & 11) Elecric Circuis n his opic, we will cover: 1) Elecromoive Force (EMF) ) Series and parallel resisor combinaions 3) Kirchhoff s rules for circuis 4) Time dependence
More informationDiebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles
Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance
More information6.2 Transforms of Derivatives and Integrals.
SEC. 6.2 Transforms of Derivaives and Inegrals. ODEs 2 3 33 39 23. Change of scale. If l( f ()) F(s) and c is any 33 45 APPLICATION OF s-shifting posiive consan, show ha l( f (c)) F(s>c)>c (Hin: In Probs.
More informationA Dynamic Model of Economic Fluctuations
CHAPTER 15 A Dynamic Model of Economic Flucuaions Modified for ECON 2204 by Bob Murphy 2016 Worh Publishers, all righs reserved IN THIS CHAPTER, OU WILL LEARN: how o incorporae dynamics ino he AD-AS model
More information18 Biological models with discrete time
8 Biological models wih discree ime The mos imporan applicaions, however, may be pedagogical. The elegan body of mahemaical heory peraining o linear sysems (Fourier analysis, orhogonal funcions, and so
More informationExponential Weighted Moving Average (EWMA) Chart Under The Assumption of Moderateness And Its 3 Control Limits
DOI: 0.545/mjis.07.5009 Exponenial Weighed Moving Average (EWMA) Char Under The Assumpion of Moderaeness And Is 3 Conrol Limis KALPESH S TAILOR Assisan Professor, Deparmen of Saisics, M. K. Bhavnagar Universiy,
More informationCE 395 Special Topics in Machine Learning
CE 395 Special Topics in Machine Learning Assoc. Prof. Dr. Yuriy Mishchenko Fall 2017 DIGITAL FILTERS AND FILTERING Why filers? Digial filering is he workhorse of digial signal processing Filering is a
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 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 informationInventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions
Muli-Period Sochasic Models: Opimali of (s, S) Polic for -Convex Objecive Funcions Consider a seing similar o he N-sage newsvendor problem excep ha now here is a fixed re-ordering cos (> 0) for each (re-)order.
More informationLinear Control System EE 711. Design. Lecture 8 Dr. Mostafa Abdel-geliel
Linear Conrol Sysem EE 7 MIMO Sae Space Analysis and Design Lecure 8 Dr. Mosafa Abdel-geliel Course Conens Review Sae Space SS modeling and analysis Sae feed back design Oupu feedback design Observer design
More informationAnnouncements: Warm-up Exercise:
Fri Apr 13 7.1 Sysems of differenial equaions - o model muli-componen sysems via comparmenal analysis hp//en.wikipedia.org/wiki/muli-comparmen_model Announcemens Warm-up Exercise Here's a relaively simple
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 informationModule 2 F c i k c s la l w a s o s f dif di fusi s o i n
Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms
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 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 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 information23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes
Half-Range Series 2.5 Inroducion In his Secion we address he following problem: Can we find a Fourier series expansion of a funcion defined over a finie inerval? Of course we recognise ha such a funcion
More informationLongest Common Prefixes
Longes Common Prefixes The sandard ordering for srings is he lexicographical order. I is induced by an order over he alphabe. We will use he same symbols (,
More informationA First Course on Kinetics and Reaction Engineering. Class 19 on Unit 18
A Firs ourse on Kineics and Reacion Engineering lass 19 on Uni 18 Par I - hemical Reacions Par II - hemical Reacion Kineics Where We re Going Par III - hemical Reacion Engineering A. Ideal Reacors B. Perfecly
More informationAssignment 6. Tyler Shendruk December 6, 2010
Assignmen 6 Tyler Shendruk December 6, 1 1 Harden Problem 1 Le K be he coupling and h he exernal field in a 1D Ising model. From he lecures hese can be ransformed ino effecive coupling and fields K and
More information6.01: Introduction to EECS I Lecture 8 March 29, 2011
6.01: Inroducion o EES I Lecure 8 March 29, 2011 6.01: Inroducion o EES I Op-Amps Las Time: The ircui Absracion ircuis represen sysems as connecions of elemens hrough which currens (hrough variables) flow
More informationSingle and Double Pendulum Models
Single and Double Pendulum Models Mah 596 Projec Summary Spring 2016 Jarod Har 1 Overview Differen ypes of pendulums are used o model many phenomena in various disciplines. In paricular, single and double
More informationNCSS Statistical Software. , contains a periodic (cyclic) component. A natural model of the periodic component would be
NCSS Saisical Sofware Chaper 468 Specral Analysis Inroducion This program calculaes and displays he periodogram and specrum of a ime series. This is someimes nown as harmonic analysis or he frequency approach
More informationProblemas das Aulas Práticas
Mesrado Inegrado em Engenharia Elecroécnica e de Compuadores Conrolo em Espaço de Esados Problemas das Aulas Práicas J. Miranda Lemos Fevereiro de 3 Translaed o English by José Gaspar, 6 J. M. Lemos, IST
More informationSolution of Integro-Differential Equations by Using ELzaki Transform
Global Journal of Mahemaical Sciences: Theory and Pracical. Volume, Number (), pp. - Inernaional Research Publicaion House hp://www.irphouse.com Soluion of Inegro-Differenial Equaions by Using ELzaki Transform
More informationEECE251. Circuit Analysis I. Set 4: Capacitors, Inductors, and First-Order Linear Circuits
EEE25 ircui Analysis I Se 4: apaciors, Inducors, and Firs-Order inear ircuis Shahriar Mirabbasi Deparmen of Elecrical and ompuer Engineering Universiy of Briish olumbia shahriar@ece.ubc.ca Overview Passive
More informationEXERCISES FOR SECTION 1.5
1.5 Exisence and Uniqueness of Soluions 43 20. 1 v c 21. 1 v c 1 2 4 6 8 10 1 2 2 4 6 8 10 Graph of approximae soluion obained using Euler s mehod wih = 0.1. Graph of approximae soluion obained using Euler
More informationWeek 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem)
Week 1 Lecure Problems, 5 Wha if somehing oscillaes wih no obvious spring? Wha is ω? (problem se problem) Sar wih Try and ge o SHM form E. Full beer can in lake, oscillaing F = m & = ge rearrange: F =
More informationOnline Convex Optimization Example And Follow-The-Leader
CSE599s, Spring 2014, Online Learning Lecure 2-04/03/2014 Online Convex Opimizaion Example And Follow-The-Leader Lecurer: Brendan McMahan Scribe: Sephen Joe Jonany 1 Review of Online Convex Opimizaion
More information) were both constant and we brought them from under the integral.
YIELD-PER-RECRUIT (coninued The yield-per-recrui model applies o a cohor, bu we saw in he Age Disribuions lecure ha he properies of a cohor do no apply in general o a collecion of cohors, which is wha
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 informationKINEMATICS IN ONE DIMENSION
KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec
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 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 informationThe fundamental mass balance equation is ( 1 ) where: I = inputs P = production O = outputs L = losses A = accumulation
Hea (iffusion) Equaion erivaion of iffusion Equaion The fundamenal mass balance equaion is I P O L A ( 1 ) where: I inpus P producion O oupus L losses A accumulaion Assume ha no chemical is produced or
More informationTime Domain Transfer Function of the Induction Motor
Sudies in Engineering and Technology Vol., No. ; Augus 0 ISSN 008 EISSN 006 Published by Redfame Publishing URL: hp://se.redfame.com Time Domain Transfer Funcion of he Inducion Moor N N arsoum Correspondence:
More informationOn-line Adaptive Optimal Timing Control of Switched Systems
On-line Adapive Opimal Timing Conrol of Swiched Sysems X.C. Ding, Y. Wardi and M. Egersed Absrac In his paper we consider he problem of opimizing over he swiching imes for a muli-modal dynamic sysem when
More informationAn recursive analytical technique to estimate time dependent physical parameters in the presence of noise processes
WHAT IS A KALMAN FILTER An recursive analyical echnique o esimae ime dependen physical parameers in he presence of noise processes Example of a ime and frequency applicaion: Offse beween wo clocks PREDICTORS,
More informationEE100 Lab 3 Experiment Guide: RC Circuits
I. Inroducion EE100 Lab 3 Experimen Guide: A. apaciors A capacior is a passive elecronic componen ha sores energy in he form of an elecrosaic field. The uni of capaciance is he farad (coulomb/vol). Pracical
More information23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes
Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals
More informationModeling the Dynamics of an Ice Tank Carriage
Modeling he Dynamics of an Ice Tank Carriage The challenge: To model he dynamics of an Ice Tank Carriage and idenify a mechanism o alleviae he backlash inheren in he design of he gearbox. Maplesof, a division
More information