Singular perturbation theory
|
|
- Charleen Sullivan
- 5 years ago
- Views:
Transcription
1 Singular perturbation theory Marc R. Rouel June 21, Introduction When we apply the teady-tate approximation (SSA) in chemical kinetic, we typically argue that ome of the intermediate are highly reactive, o that they are removed a fat a they are made. We then et the correponding rate of change to zero. What we are aying i not that thee rate are identically zero, of coure, but that they are much maller than the other rate of reaction. The teady-tate approximation i often urpriingly accurate, but the claical argument lack rigor and don t lead to any atifying mathematical generalization which can be ued to improve on the approximation. There are two mathematical theorie which build on the SSA and provide the required generalization: 1. Slow manifold theory: The curve or urface which arie from the SSA can in fact be undertood a approximation to low invariant manifol of the differential equation. We have touched on the relevant idea in the previou lecture. 2. Singular pertubation theory: Perturbation metho generally try to write olution of equation a power erie in a mall parameter. To undertand the SSA, we will find that we are led to problem in which the mall parameter multiplie a derivative. Thee problem have an intereting tructure, which i both practically and theoretically important, and to which we turn our attention in thi lecture. The two approache are not unrelated, a we hall ee in ection 6. 2 Scaling and balancing One of the problem with the claical SSA argument i that they lead u to a quandary regarding the etablihment of quantitative criteria for their validity. Conider the Michaeli-Menten mechanim: k 1 k 2 E S C E P MM k 1 In thi mechanim, E i the enzyme, S the ubtrate, P the product, and C the intermediate (enzymeubtrate/product) complex. 1
2 The claical argument run omething like thi: 1 Enzyme are highly effective catalyt. We therefore expect the degradation of C to enzyme and product to be a fat proce. Fat compared to what? Compared to formation of the enzyme-ubtrate complex, run the uual anwer. Thi argument i fine, a far a it goe, but what do we mean when we ay that one chemical proce i fater than another? What number are we comparing? It can t be rate, becaue in the teady tate, the rate of formation and of degradation are approximately the ame. It can t be rate contant, becaue k 1 ha different unit from k 1 and k 2. There i no imple way out of thi difficulty until we realize that we re going to get in trouble a long a we are trying to work with value which have unit. If we could get rid of the unit, we might be able to find thing to compare uch that we can make rigorou tatement about what it mean for one proce to be fat compared to another. Let tart by writing down the rate equation for thi mechanim: de k 1 ES k 1 C k 2 C (1a) ds k 1 ES k 1 C (1b) dc k 1ES k 1 C k 2 C (1c) dp k 2C (1d) Ma conervation give u the two equation E 0 E C S 0 S C P Note that thee conervation relation are conequence of the differential equation ince de dc ds dc dp 0. Uing thee algebraic equation therefore doen t change anything about the olution. It only reduce the et of differential equation 1 to the exactly equivalent planar ytem ds k 1 S E 0 C k1 C (2a) dc k 1S E 0 C k1 k 2 C (2b) The inight behind the SSA i that the concentration of highly reactive intermediate like C eventually reach a rough balance between production and detruction, uch that the net rate of change of C i mall. We want to make that balance evident. To do thi, we firt have to chooe new meaurement cale for each of the variable (S, C and t) uch that each i of unit magnitude. In other wor, we want to define new variable, c and τ by (2c) S S (3a) c C C (3b) and τ t t (3c) 1 The claical argument can be put a little better than thi, but the underlying problem remain. 2
3 where the tilde indicate a meaurement cale whoe unit are the ame a thoe of the correponding variable and uch that the new variable are, inofar a thi i poible, value between 0 and 1. It eay to chooe S: In mot experiment, we will tart with S 0 S 0, o S S 0 (4) hould be a good choice. Chooing C i trickier, but till not all that difficult. If we tart at the initial point SC S 0 0, which again i what we typically do in experiment, then C will initially rie to ome maximum, and then fall toward zero (the equilibrium concentration). From calculu, we know that the function C t will reach a maximum when dc 0, i.e. when C k 1 E 0 S k 1 S k 1 k 2 (5) Note that thi i the claical SSA. If we knew the value of S when C reache it maximum, we could plug it in and ue that a C. However, we have no eay way of etimating thi value of S. That being aid, we don t need to know the maximum value of C exactly. We really only need an order-of-magnitude etimate to be ued a a cale factor. We expect the maximum in C to be reached early in the reaction, before much S ha been conumed. We therefore ubtitute S S 0 into equation 5, and get k 1 E 0 S 0 E 0 S 0 C (6) k 1 S 0 k 1 k 2 S 0 K S where K S k 1 k 2 k 1 i the Michaeli contant, well known from claical enzyme kinetic. It not at all clear what we hould ue for t, o for now we ll jut ubtitute equation 3 into the rate equation 2, uing the caling factor 4 and 6, and ee what we get. S S 0 C c C c E 0S 0 S 0 K S t τ t ds d S 0 d τ t S 0 t dτ k1 S 0 E 0 E 0 S 0 c S 0 K S k E 0 S 0 1c S 0 K S k 2 dτ t k 1 E 0 S 0 1 c S 0 K S k 1c E 0 S 0 K S Similarly, for C: Let dτ tk 1 S 0 K S 1 c α S 0 S 0 K S S 0 S 0 K S c K S S 0 K S 3
4 Thi i a parameter between 0 and 1. Our two rate equation become dτ t k 1 E 0 1 E 0 αc k 1 c S 0 K S dτ tk 1 S 0 K S 1 αc c 1 α The term in the brace of the econd of thee equation are now clearly balanced: Each i made up of a product of quantitie which are of unit magnitude. There only remain to chooe t. For many enzyme, the back reaction C E S in t particularly ignificant. Under thee condition, dτ hould be dominated by the firt term, k 1 E 0 1 αc t. Thi ugget that we hould take t 1 k 1 E 0 in order to bring out the balance between the derivative and thi term. To put it another way, we expect k 1 E 0 1 to repreent a low time cale controlling the rate at which the ubtrate i ued up. Obviouly, if the above hypothei i fale for ome particular enzyme, we hould make a different choice for t. In any event, our choice lea to Define dτ 1 k 1 αc c 1 k 1 k α 2 S 0 K S 1 dτ E αc c 1 α 0 β and µ k 1 k 1 k 2 E 0 S 0 K S Then we have and dτ 1 αc βc 1 α (7a) µ dτ 1 αc c 1 α (7b) In claical enzymology, it i often aid that the SSA i valid when the enzyme concentration i mall. Thi correpon to a mall value of µ. The above caling provide a jutification for thi tatement, a well a a precie condition for it validity: If µ i very mall, then the differential equation 7b i cloe to the algebraic equation i.e. the SSA. A mall value of µ implie that 0 1 αc c 1 α (8) E 0 S 0 K S 4
5 i.e. the SSA will be valid under thi condition. Note that thi i a ufficient, but not a neceary condition. If we caled time differently, we might find other condition which lead to the validity of the SSA. Note that there omething funny about what we re doing here. We tarted out with two differential equation, and ended up aying that one of the differential equation degenerate to an algebraic equation. A ytem of two differential equation i quite different from a ingle differential equation with an algebraic contraint: The firt require the value of both and c to pecify the initial condition, while the latter require only a value for, the value of c being computed from equation 8. We have changed the nature of the problem by conidering the limit µ 0, a limit which can never be realized in practice. Problem like thi are called ingular perturbation problem. Becaue the problem we tarted with and the problem we obtain in the limit in which our mall parameter reache zero are different in nature, ingular perturbation problem can be very tricky. However, they how up all over the cience, o it well worth thinking about them a bit more. Depite the oddity of ingular perturbation problem, the ue of the limit µ 0 to approximate the low manifold (which i what we get from equation 8 if we olve for c) i upported by a theorem: Conider a ytem of differential equation of the form In the ingular limit µ dx µ dz 0, we obtain the ytem f xzt (9a) g xzt (9b) dx f xzt z φ xt where the econd of thee equation i the olution of g xzt ytem. Finally, we define the adjoined ytem, which i jut 0. Thi i called the degenerate dz g xzt where now we treat x and t a contant. 2 Theorem 1 (Tikhonov theorem) When µ 0, the olution of the ytem 9 ten to the olution of the degenerate ytem if z φ xt i a table root of the adjoined ytem. (The condition regarding the tability of the adjoined ytem i rarely checked.) The uphot of Tikhonov theorem i that for ufficiently mall µ, there i alway a low manifold whoe formula i approximately given by the SSA. Further detail can be obtained from Klonowki excellent review paper [1]. 2 If you re wondering where the µ went in the adjoined ytem, the time variable can be recaled in uch a way a to make thi coefficient diappear. Since tability doen t depend on how we cale time, there no point carrying thi parameter around. 5
6 Now imagine that we tart an experiment at the point c 10 for a ytem with a mall value of µ. Initially, becaue c 0, we will have ċ dτ 1 µ. Thi i a large (becaue µ i mall), poitive quantity. The intermediate complex concentration will therefore increae rapidly at firt, with a characteritic rie time of µ. However, a c increae, ċ will become maller. Becaue µ i mall, the right-hand ide of equation 7b hould alo remain mall, i.e. equation 8 will become valid. In other wor, we can eparate the time coure of the reaction into two part: 1. The rapid rie in the intermediate concentration: In chemical kinetic, we call thi the induction period, or ometime the tranient. In ingular perturbation theory, thi i called the inner olution A lower decay toward equilibrium during which the SSA i at leat approximately valid. Thi i the outer olution. 3 3 The outer olution In thi ection, we obtain the outer olution to lowet order in µ. Thi ue elementary technique which you hould have een in your undergraduate coure. Higher order approximation can be obtained. See Heineken, Tuchiya and Ari [2] for detail. Solving equation 8 for c, we get c α 1 α (10) Subtituting thi equation into the differential equation 7a, we get dτ 1 α α 1 α β 1 α α 1 α 1 β α 1 α Thi equation i eaily olved by eparation of variable: τ 1 α α 1 α 0 1 α 1 β dτ τ 0 1 α 1 β τ τ0 α 1 α 0 α 0 1 α ln 0 The contant 0 and t 0 repreent the tate of the ytem at ome time after it ha reached the manifold. Thee are therefore arbitrary quantitie not related (at leat not in any imple way) to 3 The term inner olution and outer olution are uch for hitorical reaon: Singular perturbation problem were firt conidered in fluid dynamic, where the inner olution referred to the behavior near a wall, while the outer olution referred to the behavior away from the wall. 6
7 the true initial condition of the ytem. We can therefore combine thee contant into a ingle, arbitrary contant A: 1 α 1 β τ A α 1 α ln (11) We will determine the contant A later. 4 The inner olution The inner olution i a bit trickier to obtain. The problem i that the equation 7 are caled wrong. Given our caling of time, the rie time for the intermediate complex i µ, which i a mall quantity. In order to tudy the inner olution, we have to recale the equation uch that the rie time i not mall. The eaiet way to do thi i a follow: τ µθ The new time variable θ i tretched by a factor of 1 µ relative to τ. Our rate equation become Setting µ and dθ µ 1 αc βc 1 α (12a) dθ 1 αc c 1 α (12b) 0, i.e. to the ame level of approximation a our outer olution, we get dθ 0 with equation 12b remaining a i. It follow that we can, to the lowet level of approximation, aume that i approximately contant during the induction period. To ditinguih thi olution from the outer olution, we define in to be the inner olution. Since S S 0, and S 0 i the initial concentration of S, in 1 to lowet order in µ. With in 1, equation 12b reduce to dθ 1 c Thi i a imple linear differential equation which can be integrated from the initial condition c 0: c 0 1 c θ θ 0 dθ ln 1 c or c 1 e θ (13) 5 Matching the inner and outer olution Since they jut repreent two piece of the ame trajectory, the inner and outer olution have to match up. Thi will provide u with condition from which we will determine the parameter A. 7
8 The outer olution i uppoed to be valid for value of τ which aren t too mall, while the inner olution i valid when τ in t too large. To match the two olution, we conider mall value of τ and large value of θ. At mall value of τ, the outer olution hould repreent early moment after the decay of the tranient, while large value of θ hould correpond to the very late tage of the induction period. If we have done everything right, we hould be able to match up the two olution in thi intermediate regime. If we let τ 0 and in in the outer olution 11, we get A α in 1 α ln in α Thu, atifie 1 α 1 β τ α α 1 α ln (14) in the outer olution region. Matching the value of c yiel no extra information in thi cae: If we put in 1 in equation 10, we get c 1. If we take the limit τ in equation 13, we get c 1. The c component of the olution therefore match up automatically. The overall procedure for doing ingular perturbation work i alway a hown here, at leat in outline: 1. Develop an appropriate caling for the equation. Identify the mall parameter(). 2. Find the outer olution. 3. Find the inner olution. 4. Match up the two olution in an intermediate time range. The only wrinkle are that we may want the olution developed to higher order than we have done here, and that the matching procedure i ometime difficult to carry out properly. 6 Geometric ingular perturbation theory and the outer olution Our work in the previou few ection focued on the temporal evolution of the ytem. Equation 14 for intance relate the time τ to in the outer olution. However, we have hinted that thee conideration are connected to the low manifold, a tructure which i bet undertood in phae pace. We now make thi connection explicit by working out a ingular perturbation expanion of the low manifold. The theoretical framework for doing thi wa firt provided by Fenichel [3] and i known a geometric ingular perturbation theory. Note that equation 10 ugget that the manifold can be written in the form c c. The manifold equation for our ytem would therefore be dτ dτ 8
9 We have dτ and dτ from equation 7. We know that, for mall µ, the manifold reduce to the form 10. One imple way to ue thi knowledge i to write the manifold a a power erie in µ: c γ i µ i (15) We hould find (and will confirm) that γ 0 i given by equation 10. In order to find the unknown function γ i (written imply a γ i from here on), we ubtitute the power erie 15 into the manifold equation. I will carry out thi procedure in tep: dτ dγ i µi 1 1 µ α i µ γ i 1 α γ i µ i µ α i µ γ i 1 1 α γ i µ i 1 µ α 1 α γ i µ i 1 dτ 1 α i µ γ i β 1 α γ i µ i α β 1 α µ α 1 α i µ γ i 1 dγ i µi γ i µ i α β 1 α γ i µ i (16) We now want to collect term in power of µ. We could jut read thee term directly from the equation. For intance, the leading term are the µ 1 term: Since thi equation mut be true for any µ 1 µ α 1 α γ 0 0 γ 0 0, the term in brace mut be zero, i.e. α 1 α (17) which i exactly what we aid thi term hould be. It not too difficult to do thi provided we re only intereted in mall power. However, if we want to develop a general equation which would allow u to obtain the erie expanion of the manifold to arbitrary order, it i convenient to rewrite equation 16 o that each ide appear a a power erie in µ, rather than a a complicated expreion involving, among other thing, a product of erie. Furthermore, it i convenient to write thee erie o that they are um involving µ i rather than, a on the left-hand ide, µ i 1 or other uch expreion. There are rule for rewriting um which you may have learned in ome of your 9
10 mathematic coure. If not, the following example will hopefully let you ee the logic behind thee rule, which in t too hard once you ve done a few example. The firt um on the left-hand ide of equation 16 i fairly eay to rewrite in the deired form: γ i µ i 1 i 1 γ i 1µ i If you can t ee thi immediately, try writing out the firt few term of each ide of thi identity. The trickiet bit i the product of um on the right-hand ide of equation 16. Here i the tranformation: j 0 dγ j µj k µ k 0γ k µ i i dγ j j 0 γ i j The renaming of the ummation indicie on the left-hand ide i a convenience introduced for clarity. The key obervation i that multiplying the two um produce term in µ i when j k i. We re now ready to go back to equation 16: µ α 1 α γ i 1µ i i 1 dγ i µi α β 1 α µ i i dγ j j 0 γ i j Reading off term in each power of µ i now eay: µ 1 : γ 0 α 1 α 0 µ i i 0 : γ i 1 α 1 α dγ i α β 1 α i dγ j j 0 γ i j We have een the firt of thee equation, and even written down it olution (equation 17). The econd equation can be rewritten in the form γ i 1 α 1 1 α dγ i α β 1 α i dγ j j 0 γ i j Note that we only need to know the formula of the γ up to i to determine the i erie. For intance, γ 1 i calculated by 1 t term in the 1 γ 1 α 1 α dγ 0 1 α 2 1 β α 1 α 4 α β 1 α dγ 0 γ 0 A you can imagine, thee calculation get tediou pretty quickly. In the old day, we ued to fill up notebook with perturbation erie, carefully implified by hand. It wa very error-prone, and we would often pend more time verifying our anwer than we did doing the initial calculation. Nowaday, we jut ue Maple. It much fater, but it till a good idea to calculate a few term by hand to check that Maple i doing the calculation correctly. 10
11 Reference [1] W. Klonowki, Biophy. Chem. 18, 73 (1983). [2] F. Heineken, H. Tuchiya, and R. Ari, Math. Bioci. 1, 95 (1967). [3] N. Fenichel, J. Diff. Eq. 31, 53 (1979). 11
Singular perturbation theory
Singular perturbation theory Marc R Roussel October 19, 2005 1 Introduction When we apply the steady-state approximation (SSA) in chemical kinetics, we typically argue that some of the intermediates are
More informationBogoliubov Transformation in Classical Mechanics
Bogoliubov Tranformation in Claical Mechanic Canonical Tranformation Suppoe we have a et of complex canonical variable, {a j }, and would like to conider another et of variable, {b }, b b ({a j }). How
More informationLaplace Transformation
Univerity of Technology Electromechanical Department Energy Branch Advance Mathematic Laplace Tranformation nd Cla Lecture 6 Page of 7 Laplace Tranformation Definition Suppoe that f(t) i a piecewie continuou
More informationCorrection for Simple System Example and Notes on Laplace Transforms / Deviation Variables ECHE 550 Fall 2002
Correction for Simple Sytem Example and Note on Laplace Tranform / Deviation Variable ECHE 55 Fall 22 Conider a tank draining from an initial height of h o at time t =. With no flow into the tank (F in
More information7.2 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 281
72 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 28 and i 2 Show how Euler formula (page 33) can then be ued to deduce the reult a ( a) 2 b 2 {e at co bt} {e at in bt} b ( a) 2 b 2 5 Under what condition
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS Laplace Tranform Paul Dawkin Table of Content Preface... Laplace Tranform... Introduction... The Definition... 5 Laplace Tranform... 9 Invere Laplace Tranform... Step Function...4
More informationNCAAPMT Calculus Challenge Challenge #3 Due: October 26, 2011
NCAAPMT Calculu Challenge 011 01 Challenge #3 Due: October 6, 011 A Model of Traffic Flow Everyone ha at ome time been on a multi-lane highway and encountered road contruction that required the traffic
More informationSocial Studies 201 Notes for November 14, 2003
1 Social Studie 201 Note for November 14, 2003 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the
More informationDIFFERENTIAL EQUATIONS Laplace Transforms. Paul Dawkins
DIFFERENTIAL EQUATIONS Laplace Tranform Paul Dawkin Table of Content Preface... Laplace Tranform... Introduction... The Definition... 5 Laplace Tranform... 9 Invere Laplace Tranform... Step Function...
More informationSocial Studies 201 Notes for March 18, 2005
1 Social Studie 201 Note for March 18, 2005 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the
More informationNonlinear Single-Particle Dynamics in High Energy Accelerators
Nonlinear Single-Particle Dynamic in High Energy Accelerator Part 6: Canonical Perturbation Theory Nonlinear Single-Particle Dynamic in High Energy Accelerator Thi coure conit of eight lecture: 1. Introduction
More informationμ + = σ = D 4 σ = D 3 σ = σ = All units in parts (a) and (b) are in V. (1) x chart: Center = μ = 0.75 UCL =
Our online Tutor are available 4*7 to provide Help with Proce control ytem Homework/Aignment or a long term Graduate/Undergraduate Proce control ytem Project. Our Tutor being experienced and proficient
More informationMarch 18, 2014 Academic Year 2013/14
POLITONG - SHANGHAI BASIC AUTOMATIC CONTROL Exam grade March 8, 4 Academic Year 3/4 NAME (Pinyin/Italian)... STUDENT ID Ue only thee page (including the back) for anwer. Do not ue additional heet. Ue of
More informationMATEMATIK Datum: Tid: eftermiddag. A.Heintz Telefonvakt: Anders Martinsson Tel.:
MATEMATIK Datum: 20-08-25 Tid: eftermiddag GU, Chalmer Hjälpmedel: inga A.Heintz Telefonvakt: Ander Martinon Tel.: 073-07926. Löningar till tenta i ODE och matematik modellering, MMG5, MVE6. Define what
More informationLecture 9: Shor s Algorithm
Quantum Computation (CMU 8-859BB, Fall 05) Lecture 9: Shor Algorithm October 7, 05 Lecturer: Ryan O Donnell Scribe: Sidhanth Mohanty Overview Let u recall the period finding problem that wa et up a a function
More informationLecture 8: Period Finding: Simon s Problem over Z N
Quantum Computation (CMU 8-859BB, Fall 205) Lecture 8: Period Finding: Simon Problem over Z October 5, 205 Lecturer: John Wright Scribe: icola Rech Problem A mentioned previouly, period finding i a rephraing
More informationIntroduction to Laplace Transform Techniques in Circuit Analysis
Unit 6 Introduction to Laplace Tranform Technique in Circuit Analyi In thi unit we conider the application of Laplace Tranform to circuit analyi. A relevant dicuion of the one-ided Laplace tranform i found
More informationConstant Force: Projectile Motion
Contant Force: Projectile Motion Abtract In thi lab, you will launch an object with a pecific initial velocity (magnitude and direction) and determine the angle at which the range i a maximum. Other tak,
More informationChapter 4. The Laplace Transform Method
Chapter 4. The Laplace Tranform Method The Laplace Tranform i a tranformation, meaning that it change a function into a new function. Actually, it i a linear tranformation, becaue it convert a linear combination
More informationDimensional Analysis A Tool for Guiding Mathematical Calculations
Dimenional Analyi A Tool for Guiding Mathematical Calculation Dougla A. Kerr Iue 1 February 6, 2010 ABSTRACT AND INTRODUCTION In converting quantitie from one unit to another, we may know the applicable
More informationGiven the following circuit with unknown initial capacitor voltage v(0): X(s) Immediately, we know that the transfer function H(s) is
EE 4G Note: Chapter 6 Intructor: Cheung More about ZSR and ZIR. Finding unknown initial condition: Given the following circuit with unknown initial capacitor voltage v0: F v0/ / Input xt 0Ω Output yt -
More informationMoment of Inertia of an Equilateral Triangle with Pivot at one Vertex
oment of nertia of an Equilateral Triangle with Pivot at one Vertex There are two wa (at leat) to derive the expreion f an equilateral triangle that i rotated about one vertex, and ll how ou both here.
More informationLecture 3. January 9, 2018
Lecture 3 January 9, 208 Some complex analyi Although you might have never taken a complex analyi coure, you perhap till know what a complex number i. It i a number of the form z = x + iy, where x and
More informationFermi Distribution Function. n(e) T = 0 T > 0 E F
LECTURE 3 Maxwell{Boltzmann, Fermi, and Boe Statitic Suppoe we have a ga of N identical point particle in a box ofvolume V. When we ay \ga", we mean that the particle are not interacting with one another.
More informationDepartment of Mechanical Engineering Massachusetts Institute of Technology Modeling, Dynamics and Control III Spring 2002
Department of Mechanical Engineering Maachuett Intitute of Technology 2.010 Modeling, Dynamic and Control III Spring 2002 SOLUTIONS: Problem Set # 10 Problem 1 Etimating tranfer function from Bode Plot.
More informationChapter 13. Root Locus Introduction
Chapter 13 Root Locu 13.1 Introduction In the previou chapter we had a glimpe of controller deign iue through ome imple example. Obviouly when we have higher order ytem, uch imple deign technique will
More informationComputers and Mathematics with Applications. Sharp algebraic periodicity conditions for linear higher order
Computer and Mathematic with Application 64 (2012) 2262 2274 Content lit available at SciVere ScienceDirect Computer and Mathematic with Application journal homepage: wwweleviercom/locate/camwa Sharp algebraic
More informationSolving Differential Equations by the Laplace Transform and by Numerical Methods
36CH_PHCalter_TechMath_95099 3//007 :8 PM Page Solving Differential Equation by the Laplace Tranform and by Numerical Method OBJECTIVES When you have completed thi chapter, you hould be able to: Find the
More informationPhysics 741 Graduate Quantum Mechanics 1 Solutions to Final Exam, Fall 2014
Phyic 7 Graduate Quantum Mechanic Solution to inal Eam all 0 Each quetion i worth 5 point with point for each part marked eparately Some poibly ueful formula appear at the end of the tet In four dimenion
More informationMidterm Review - Part 1
Honor Phyic Fall, 2016 Midterm Review - Part 1 Name: Mr. Leonard Intruction: Complete the following workheet. SHOW ALL OF YOUR WORK. 1. Determine whether each tatement i True or Fale. If the tatement i
More informationME 375 FINAL EXAM SOLUTIONS Friday December 17, 2004
ME 375 FINAL EXAM SOLUTIONS Friday December 7, 004 Diviion Adam 0:30 / Yao :30 (circle one) Name Intruction () Thi i a cloed book eamination, but you are allowed three 8.5 crib heet. () You have two hour
More informationLecture 17: Analytic Functions and Integrals (See Chapter 14 in Boas)
Lecture 7: Analytic Function and Integral (See Chapter 4 in Boa) Thi i a good point to take a brief detour and expand on our previou dicuion of complex variable and complex function of complex variable.
More informationp. (The electron is a point particle with radius r = 0.)
- pin ½ Recall that in the H-atom olution, we howed that the fact that the wavefunction Ψ(r) i ingle-valued require that the angular momentum quantum nbr be integer: l = 0,,.. However, operator algebra
More informationNAME (pinyin/italian)... MATRICULATION NUMBER... SIGNATURE
POLITONG SHANGHAI BASIC AUTOMATIC CONTROL June Academic Year / Exam grade NAME (pinyin/italian)... MATRICULATION NUMBER... SIGNATURE Ue only thee page (including the bac) for anwer. Do not ue additional
More informationLecture 15 - Current. A Puzzle... Advanced Section: Image Charge for Spheres. Image Charge for a Grounded Spherical Shell
Lecture 15 - Current Puzzle... Suppoe an infinite grounded conducting plane lie at z = 0. charge q i located at a height h above the conducting plane. Show in three different way that the potential below
More informationCHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS
CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.1 INTRODUCTION 8.2 REDUCED ORDER MODEL DESIGN FOR LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.3
More informationA Simplified Methodology for the Synthesis of Adaptive Flight Control Systems
A Simplified Methodology for the Synthei of Adaptive Flight Control Sytem J.ROUSHANIAN, F.NADJAFI Department of Mechanical Engineering KNT Univerity of Technology 3Mirdamad St. Tehran IRAN Abtract- A implified
More informationPHYS 110B - HW #6 Spring 2004, Solutions by David Pace Any referenced equations are from Griffiths Problem statements are paraphrased
PHYS B - HW #6 Spring 4, Solution by David Pace Any referenced equation are from Griffith Problem tatement are paraphraed. Problem. from Griffith Show that the following, A µo ɛ o A V + A ρ ɛ o Eq..4 A
More informationMolecular Dynamics Simulations of Nonequilibrium Effects Associated with Thermally Activated Exothermic Reactions
Original Paper orma, 5, 9 7, Molecular Dynamic Simulation of Nonequilibrium Effect ociated with Thermally ctivated Exothermic Reaction Jerzy GORECKI and Joanna Natalia GORECK Intitute of Phyical Chemitry,
More informationMath Skills. Scientific Notation. Uncertainty in Measurements. Appendix A5 SKILLS HANDBOOK
ppendix 5 Scientific Notation It i difficult to work with very large or very mall number when they are written in common decimal notation. Uually it i poible to accommodate uch number by changing the SI
More informationThe Laplace Transform , Haynes Miller and Jeremy Orloff
The Laplace Tranform 8.3, Hayne Miller and Jeremy Orloff Laplace tranform baic: introduction An operator take a function a input and output another function. A tranform doe the ame thing with the added
More informationinto a discrete time function. Recall that the table of Laplace/z-transforms is constructed by (i) selecting to get
Lecture 25 Introduction to Some Matlab c2d Code in Relation to Sampled Sytem here are many way to convert a continuou time function, { h( t) ; t [0, )} into a dicrete time function { h ( k) ; k {0,,, }}
More informationON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION. Xiaoqun Wang
Proceeding of the 2008 Winter Simulation Conference S. J. Maon, R. R. Hill, L. Mönch, O. Roe, T. Jefferon, J. W. Fowler ed. ON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION Xiaoqun Wang
More informationLecture 21. The Lovasz splitting-off lemma Topics in Combinatorial Optimization April 29th, 2004
18.997 Topic in Combinatorial Optimization April 29th, 2004 Lecture 21 Lecturer: Michel X. Goeman Scribe: Mohammad Mahdian 1 The Lovaz plitting-off lemma Lovaz plitting-off lemma tate the following. Theorem
More informationProblem Set 8 Solutions
Deign and Analyi of Algorithm April 29, 2015 Maachuett Intitute of Technology 6.046J/18.410J Prof. Erik Demaine, Srini Devada, and Nancy Lynch Problem Set 8 Solution Problem Set 8 Solution Thi problem
More informationThe continuous time random walk (CTRW) was introduced by Montroll and Weiss 1.
1 I. CONTINUOUS TIME RANDOM WALK The continuou time random walk (CTRW) wa introduced by Montroll and Wei 1. Unlike dicrete time random walk treated o far, in the CTRW the number of jump n made by the walker
More informationHomework #7 Solution. Solutions: ΔP L Δω. Fig. 1
Homework #7 Solution Aignment:. through.6 Bergen & Vittal. M Solution: Modified Equation.6 becaue gen. peed not fed back * M (.0rad / MW ec)(00mw) rad /ec peed ( ) (60) 9.55r. p. m. 3600 ( 9.55) 3590.45r.
More informationRiemann s Functional Equation is Not a Valid Function and Its Implication on the Riemann Hypothesis. Armando M. Evangelista Jr.
Riemann Functional Equation i Not a Valid Function and It Implication on the Riemann Hypothei By Armando M. Evangelita Jr. armando78973@gmail.com On Augut 28, 28 ABSTRACT Riemann functional equation wa
More informationAssignment for Mathematics for Economists Fall 2016
Due date: Mon. Nov. 1. Reading: CSZ, Ch. 5, Ch. 8.1 Aignment for Mathematic for Economit Fall 016 We now turn to finihing our coverage of concavity/convexity. There are two part: Jenen inequality for concave/convex
More informationV = 4 3 πr3. d dt V = d ( 4 dv dt. = 4 3 π d dt r3 dv π 3r2 dv. dt = 4πr 2 dr
0.1 Related Rate In many phyical ituation we have a relationhip between multiple quantitie, and we know the rate at which one of the quantitie i changing. Oftentime we can ue thi relationhip a a convenient
More informationSolutions to homework #10
Solution to homework #0 Problem 7..3 Compute 6 e 3 t t t 8. The firt tep i to ue the linearity of the Laplace tranform to ditribute the tranform over the um and pull the contant factor outide the tranform.
More informationonline learning Unit Workbook 4 RLC Transients
online learning Pearon BTC Higher National in lectrical and lectronic ngineering (QCF) Unit 5: lectrical & lectronic Principle Unit Workbook 4 in a erie of 4 for thi unit Learning Outcome: RLC Tranient
More informationHalliday/Resnick/Walker 7e Chapter 6
HRW 7e Chapter 6 Page of Halliday/Renick/Walker 7e Chapter 6 3. We do not conider the poibility that the bureau might tip, and treat thi a a purely horizontal motion problem (with the peron puh F in the
More informationA Constraint Propagation Algorithm for Determining the Stability Margin. The paper addresses the stability margin assessment for linear systems
A Contraint Propagation Algorithm for Determining the Stability Margin of Linear Parameter Circuit and Sytem Lubomir Kolev and Simona Filipova-Petrakieva Abtract The paper addree the tability margin aement
More informationEigenvalues and eigenvectors
Eigenvalue and eigenvector Defining and computing uggeted problem olution For each matri give below, find eigenvalue and eigenvector. Give a bai and the dimenion of the eigenpace for each eigenvalue. P:
More informationIEOR 3106: Fall 2013, Professor Whitt Topics for Discussion: Tuesday, November 19 Alternating Renewal Processes and The Renewal Equation
IEOR 316: Fall 213, Profeor Whitt Topic for Dicuion: Tueday, November 19 Alternating Renewal Procee and The Renewal Equation 1 Alternating Renewal Procee An alternating renewal proce alternate between
More informationNotes on Phase Space Fall 2007, Physics 233B, Hitoshi Murayama
Note on Phae Space Fall 007, Phyic 33B, Hitohi Murayama Two-Body Phae Space The two-body phae i the bai of computing higher body phae pace. We compute it in the ret frame of the two-body ytem, P p + p
More informationManprit Kaur and Arun Kumar
CUBIC X-SPLINE INTERPOLATORY FUNCTIONS Manprit Kaur and Arun Kumar manpreet2410@gmail.com, arun04@rediffmail.com Department of Mathematic and Computer Science, R. D. Univerity, Jabalpur, INDIA. Abtract:
More informationRiemann s Functional Equation is Not Valid and its Implication on the Riemann Hypothesis. Armando M. Evangelista Jr.
Riemann Functional Equation i Not Valid and it Implication on the Riemann Hypothei By Armando M. Evangelita Jr. On November 4, 28 ABSTRACT Riemann functional equation wa formulated by Riemann that uppoedly
More informationHELICAL TUBES TOUCHING ONE ANOTHER OR THEMSELVES
15 TH INTERNATIONAL CONFERENCE ON GEOMETRY AND GRAPHICS 0 ISGG 1-5 AUGUST, 0, MONTREAL, CANADA HELICAL TUBES TOUCHING ONE ANOTHER OR THEMSELVES Peter MAYRHOFER and Dominic WALTER The Univerity of Innbruck,
More information1. The F-test for Equality of Two Variances
. The F-tet for Equality of Two Variance Previouly we've learned how to tet whether two population mean are equal, uing data from two independent ample. We can alo tet whether two population variance are
More informationJump condition at the boundary between a porous catalyst and a homogeneous fluid
From the SelectedWork of Francico J. Valde-Parada 2005 Jump condition at the boundary between a porou catalyt and a homogeneou fluid Francico J. Valde-Parada J. Alberto Ochoa-Tapia Available at: http://work.bepre.com/francico_j_valde_parada/12/
More informationAutomatic Control Systems. Part III: Root Locus Technique
www.pdhcenter.com PDH Coure E40 www.pdhonline.org Automatic Control Sytem Part III: Root Locu Technique By Shih-Min Hu, Ph.D., P.E. Page of 30 www.pdhcenter.com PDH Coure E40 www.pdhonline.org VI. Root
More informationLecture 10 Filtering: Applied Concepts
Lecture Filtering: Applied Concept In the previou two lecture, you have learned about finite-impule-repone (FIR) and infinite-impule-repone (IIR) filter. In thee lecture, we introduced the concept of filtering
More informationME 375 FINAL EXAM Wednesday, May 6, 2009
ME 375 FINAL EXAM Wedneday, May 6, 9 Diviion Meckl :3 / Adam :3 (circle one) Name_ Intruction () Thi i a cloed book examination, but you are allowed three ingle-ided 8.5 crib heet. A calculator i NOT allowed.
More informationLecture 4. Chapter 11 Nise. Controller Design via Frequency Response. G. Hovland 2004
METR4200 Advanced Control Lecture 4 Chapter Nie Controller Deign via Frequency Repone G. Hovland 2004 Deign Goal Tranient repone via imple gain adjutment Cacade compenator to improve teady-tate error Cacade
More informationGain and Phase Margins Based Delay Dependent Stability Analysis of Two- Area LFC System with Communication Delays
Gain and Phae Margin Baed Delay Dependent Stability Analyi of Two- Area LFC Sytem with Communication Delay Şahin Sönmez and Saffet Ayaun Department of Electrical Engineering, Niğde Ömer Halidemir Univerity,
More informationRoot Locus Diagram. Root loci: The portion of root locus when k assume positive values: that is 0
Objective Root Locu Diagram Upon completion of thi chapter you will be able to: Plot the Root Locu for a given Tranfer Function by varying gain of the ytem, Analye the tability of the ytem from the root
More informationNOTE: The items d) and e) of Question 4 gave you bonus marks.
MAE 40 Linear ircuit Summer 2007 Final Solution NOTE: The item d) and e) of Quetion 4 gave you bonu mark. Quetion [Equivalent irciut] [4 mark] Find the equivalent impedance between terminal A and B in
More informationThe machines in the exercise work as follows:
Tik-79.148 Spring 2001 Introduction to Theoretical Computer Science Tutorial 9 Solution to Demontration Exercie 4. Contructing a complex Turing machine can be very laboriou. With the help of machine chema
More informationCHAPTER 4 DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL
98 CHAPTER DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL INTRODUCTION The deign of ytem uing tate pace model for the deign i called a modern control deign and it i
More informationLecture 13. Thermodynamic Potentials (Ch. 5)
Lecture 13. hermodynamic Potential (Ch. 5) So far we have been uing the total internal energy U and ometime the enthalpy H to characterize variou macrocopic ytem. hee function are called the thermodynamic
More informationLecture 7: Testing Distributions
CSE 5: Sublinear (and Streaming) Algorithm Spring 014 Lecture 7: Teting Ditribution April 1, 014 Lecturer: Paul Beame Scribe: Paul Beame 1 Teting Uniformity of Ditribution We return today to property teting
More informationCodes Correcting Two Deletions
1 Code Correcting Two Deletion Ryan Gabry and Frederic Sala Spawar Sytem Center Univerity of California, Lo Angele ryan.gabry@navy.mil fredala@ucla.edu Abtract In thi work, we invetigate the problem of
More informationCalculation of the temperature of boundary layer beside wall with time-dependent heat transfer coefficient
Ŕ periodica polytechnica Mechanical Engineering 54/1 21 15 2 doi: 1.3311/pp.me.21-1.3 web: http:// www.pp.bme.hu/ me c Periodica Polytechnica 21 RESERCH RTICLE Calculation of the temperature of boundary
More informationMath 273 Solutions to Review Problems for Exam 1
Math 7 Solution to Review Problem for Exam True or Fale? Circle ONE anwer for each Hint: For effective tudy, explain why if true and give a counterexample if fale (a) T or F : If a b and b c, then a c
More informationFigure 1 Siemens PSSE Web Site
Stability Analyi of Dynamic Sytem. In the lat few lecture we have een how mall ignal Lalace domain model may be contructed of the dynamic erformance of ower ytem. The tability of uch ytem i a matter of
More informationSolving Radical Equations
10. Solving Radical Equation Eential Quetion How can you olve an equation that contain quare root? Analyzing a Free-Falling Object MODELING WITH MATHEMATICS To be proficient in math, you need to routinely
More informationSERIES COMPENSATION: VOLTAGE COMPENSATION USING DVR (Lectures 41-48)
Chapter 5 SERIES COMPENSATION: VOLTAGE COMPENSATION USING DVR (Lecture 41-48) 5.1 Introduction Power ytem hould enure good quality of electric power upply, which mean voltage and current waveform hould
More informationReading assignment: In this chapter we will cover Sections Definition and the Laplace transform of simple functions
Chapter 4 Laplace Tranform 4 Introduction Reading aignment: In thi chapter we will cover Section 4 45 4 Definition and the Laplace tranform of imple function Given f, a function of time, with value f(t
More informationPreemptive scheduling on a small number of hierarchical machines
Available online at www.ciencedirect.com Information and Computation 06 (008) 60 619 www.elevier.com/locate/ic Preemptive cheduling on a mall number of hierarchical machine György Dóa a, Leah Eptein b,
More informationON THE DETERMINATION OF NUMBERS BY THEIR SUMS OF A FIXED ORDER J. L. SELFRIDGE AND E. G. STRAUS
ON THE DETERMINATION OF NUMBERS BY THEIR SUMS OF A FIXED ORDER J. L. SELFRIDGE AND E. G. STRAUS 1. Introduction. We wih to treat the following problem (uggeted by a problem of L. Moer [2]): Let {x} = {x
More informationChapter 2 Sampling and Quantization. In order to investigate sampling and quantization, the difference between analog
Chapter Sampling and Quantization.1 Analog and Digital Signal In order to invetigate ampling and quantization, the difference between analog and digital ignal mut be undertood. Analog ignal conit of continuou
More informationNew bounds for Morse clusters
New bound for More cluter Tamá Vinkó Advanced Concept Team, European Space Agency, ESTEC Keplerlaan 1, 2201 AZ Noordwijk, The Netherland Tama.Vinko@ea.int and Arnold Neumaier Fakultät für Mathematik, Univerität
More informationPythagorean Triple Updated 08--5 Drlnoordzij@leennoordzijnl wwwleennoordzijme Content A Roadmap for generating Pythagorean Triple Pythagorean Triple 3 Dicuion Concluion 5 A Roadmap for generating Pythagorean
More information1. Preliminaries. In [8] the following odd looking integral evaluation is obtained.
June, 5. Revied Augut 8th, 5. VA DER POL EXPASIOS OF L-SERIES David Borwein* and Jonathan Borwein Abtract. We provide concie erie repreentation for variou L-erie integral. Different technique are needed
More informationarxiv: v4 [math.co] 21 Sep 2014
ASYMPTOTIC IMPROVEMENT OF THE SUNFLOWER BOUND arxiv:408.367v4 [math.co] 2 Sep 204 JUNICHIRO FUKUYAMA Abtract. A unflower with a core Y i a family B of et uch that U U Y for each two different element U
More informationDesign By Emulation (Indirect Method)
Deign By Emulation (Indirect Method he baic trategy here i, that Given a continuou tranfer function, it i required to find the bet dicrete equivalent uch that the ignal produced by paing an input ignal
More informationConvex Hulls of Curves Sam Burton
Convex Hull of Curve Sam Burton 1 Introduction Thi paper will primarily be concerned with determining the face of convex hull of curve of the form C = {(t, t a, t b ) t [ 1, 1]}, a < b N in R 3. We hall
More informationLecture 7 Grain boundary grooving
Lecture 7 Grain oundary grooving The phenomenon. A polihed polycrytal ha a flat urface. At room temperature, the urface remain flat for a long time. At an elevated temperature atom move. The urface grow
More informationCHEAP CONTROL PERFORMANCE LIMITATIONS OF INPUT CONSTRAINED LINEAR SYSTEMS
Copyright 22 IFAC 5th Triennial World Congre, Barcelona, Spain CHEAP CONTROL PERFORMANCE LIMITATIONS OF INPUT CONSTRAINED LINEAR SYSTEMS Tritan Pérez Graham C. Goodwin Maria M. Serón Department of Electrical
More informationAn Inequality for Nonnegative Matrices and the Inverse Eigenvalue Problem
An Inequality for Nonnegative Matrice and the Invere Eigenvalue Problem Robert Ream Program in Mathematical Science The Univerity of Texa at Dalla Box 83688, Richardon, Texa 7583-688 Abtract We preent
More informationEE Control Systems LECTURE 6
Copyright FL Lewi 999 All right reerved EE - Control Sytem LECTURE 6 Updated: Sunday, February, 999 BLOCK DIAGRAM AND MASON'S FORMULA A linear time-invariant (LTI) ytem can be repreented in many way, including:
More information6.302 Feedback Systems Recitation 6: Steady-State Errors Prof. Joel L. Dawson S -
6302 Feedback ytem Recitation 6: teadytate Error Prof Joel L Dawon A valid performance metric for any control ytem center around the final error when the ytem reache teadytate That i, after all initial
More information(b) Is the game below solvable by iterated strict dominance? Does it have a unique Nash equilibrium?
14.1 Final Exam Anwer all quetion. You have 3 hour in which to complete the exam. 1. (60 Minute 40 Point) Anwer each of the following ubquetion briefly. Pleae how your calculation and provide rough explanation
More informationExercises for lectures 19 Polynomial methods
Exercie for lecture 19 Polynomial method Michael Šebek Automatic control 016 15-4-17 Diviion of polynomial with and without remainder Polynomial form a circle, but not a body. (Circle alo form integer,
More informationUSPAS Course on Recirculated and Energy Recovered Linear Accelerators
USPAS Coure on Recirculated and Energy Recovered Linear Accelerator G. A. Krafft and L. Merminga Jefferon Lab I. Bazarov Cornell Lecture 6 7 March 005 Lecture Outline. Invariant Ellipe Generated by a Unimodular
More informationLecture 2 Phys 798S Spring 2016 Steven Anlage. The heart and soul of superconductivity is the Meissner Effect. This feature uniquely distinguishes
ecture Phy 798S Spring 6 Steven Anlage The heart and oul of uperconductivity i the Meiner Effect. Thi feature uniquely ditinguihe uperconductivity fro any other tate of atter. Here we dicu oe iple phenoenological
More informationHyperbolic Partial Differential Equations
Hyperbolic Partial Differential Equation Evolution equation aociated with irreverible phyical procee like diffuion heat conduction lead to parabolic partial differential equation. When the equation i a
More informationQuestion 1 Equivalent Circuits
MAE 40 inear ircuit Fall 2007 Final Intruction ) Thi exam i open book You may ue whatever written material you chooe, including your cla note and textbook You may ue a hand calculator with no communication
More information