Solution Methods

Total Page:16

File Type:pdf, Size:1020Kb

Solution Methods Solution Methods Jesús Fernández-Villaverde University of Pennsylvania March 16, 2016 Jesús Fernández-Villaverde (PENN) Solution Methods March 16, 2016 1 / 36 Functional equations A large class of problems in macroeconomics search for a function d that solves a functional equation: (d) = 0 H More formally: 1 Let J1 and J2 be two functional spaces and let : J1 J2 be an operator between these two spaces. H ! 2 Let W Rl . 3 Then, we need to find a function d : W Rm such that (d) = 0. ! H Notes: 1 Regular equations are particular examples of functional equations. 2 0 is the space zero, different in general that the zero in the reals. Jesús Fernández-Villaverde (PENN) Solution Methods March 16, 2016 2 / 36 Example I Let’sgo back to our basic stochastic neoclassical growth model: ¥ t max E0 ∑ b u (ct , lt ) t=0 zt a 1 a ct + kt+1 = e Ak l + (1 d) kt , t > 0 t t 8 zt = rzt 1 + s#t , #t (0, 1) N 1 2 2 log st = (1 r ) log s + r log st 1 + 1 r hut s s s The first order condition: zt+1 a 1 1 a u0 (ct , lt ) = bEt u0 (ct+1, lt+1) 1 + ae Ak l d t+1 t+1 Where is the stochastic volatility? Jesús Fernández-Villaverde (PENN) Solution Methods March 16, 2016 3 / 36 Example II Define: xt = (kt , zt 1, log st 1, #t , ut ) There is a decision rule (a.k.a. policy function) that gives the optimal choice of consumption and capital tomorrow given the states today: 1 d (xt ) = lt d = 2 d (xt ) = kt+1 From these two choices, we can find: zt a 1 1 a 2 ct = e Ak d (xt ) + (1 d) kt d (xt ) t Jesús Fernández-Villaverde (PENN) Solution Methods March 16, 2016 4 / 36 Example III Then: 1 = u0 ct , d (xt ) H 1 u0 ct+1, d (xt+1) bEt z 2 a 1 1 1 a = 0 1 + ae t+1 Ad (xt ) d (xt+1) d ( ) If we find g, and a transversality condition is satisfied, we are done! Jesús Fernández-Villaverde (PENN) Solution Methods March 16, 2016 5 / 36 Example IV There is a recursive problem associated with the previous sequential problem: V (xt ) = max u (ct , lt ) + bEt V (xt+1) kt+1,lt f g zt a 1 a ct + kt+1 = e Ak l + (1 d) kt , t > 0 t t 8 zt = rzt 1 + s#t , #t (0, 1) N 1 2 2 log st = (1 r ) log s + r log st 1 + 1 r hut s s s Then: d (xt ) = V (xt ) and (d) = d (xt ) max u (ct , lt ) + bEt d (xt+1) = 0 H kt+1,lt f g Jesús Fernández-Villaverdee (PENN) Solution Methods March 16, 2016 6 / 36 How do we solve functional equations? General idea: substitute d (x) by dn (x, q) where q is an n dim vector of coeffi cients to be determined. Two Main Approaches: 1 Perturbation methods: n n i d (x, q) = ∑ qi (x x0) i=0 We use implicit-function theorems to find qi . 2 Projection methods: n n d (x, q) = ∑ qi Yi (x) i=0 ¥ We pick a basis Yi (x) and “project” ( ) against that basis. f gi=0 H · Jesús Fernández-Villaverde (PENN) Solution Methods March 16, 2016 7 / 36 Comparison with traditional solution methods Linearization (or loglinearization): equivalent to a first-order perturbation. Linear-quadratic approximation: equivalent (under certain conditions) to a first-order perturbation. Parameterized expectations: a particular example of projection. Value function iteration:it can be interpreted as an iterative procedure to solve a particular projection method. Nevertheless, I prefer to think about it as a different family of problems. Policy function iteration: similar to VFI. Jesús Fernández-Villaverde (PENN) Solution Methods March 16, 2016 8 / 36 Advantages of the functional equation approach Generality: abstract framework highlights commonalities across problems. Large set of existing theoretical and numerical results in applied math. It allows us to identify more clearly issue and challenges specific to economic problems (for example, importance of expectations). It allows us to deal effi ciently with nonlinearities. Jesús Fernández-Villaverde (PENN) Solution Methods March 16, 2016 9 / 36 Perturbation: motivation Perturbation builds a Taylor-series approximation of the exact solution. Very accurate around the point where the approximation is undertakes. Often, surprisingly good global properties. Only approach that handle models with dozens of state variables. Relation between uncertainty shocks and the curse of dimensionality. Jesús Fernández-Villaverde (PENN) Solution Methods March 16, 2016 10 / 36 References General: 1 A First Look at Perturbation Theory by James G. Simmonds and James E. Mann Jr. 2 Advanced Mathematical Methods for Scientists and Engineers: Asymptotic Methods and Perturbation Theory by Carl M. Bender, Steven A. Orszag. Economics: 1 Perturbation Methods for General Dynamic Stochastic Models” by Hehui Jin and Kenneth Judd. 2 Perturbation Methods with Nonlinear Changes of Variables” by Kenneth Judd. 3 A gentle introduction: “Solving Dynamic General Equilibrium Models Using a Second-Order Approximation to the Policy Function” by Martín Uribe and Stephanie Schmitt-Grohe. Jesús Fernández-Villaverde (PENN) Solution Methods March 16, 2016 11 / 36 Our RBC-SV model Let me come back to our RBC-SV model. Three changes: 1 Eliminate labor supply and have a log utility function for consumption. 2 Full depreciation. 3 A = 1. Why? Jesús Fernández-Villaverde (PENN) Solution Methods March 16, 2016 12 / 36 Environment ¥ t max E0 ∑ b log ct t=0 zt a s.t. ct + kt+1 = e k , t > 0 t 8 zt = rzt 1 + s#t 1 2 2 log st = (1 r ) log s + r log st 1 + 1 r hut s s s Equilibrium conditions: 1 1 zt+1 a 1 = bEt ae kt+1 ct ct+1 zt a ct + kt+1 = e kt zt = rzt 1 + s#t 1 2 2 log st = (1 r ) log s + r log st 1 + 1 r hut s s s Jesús Fernández-Villaverde (PENN) Solution Methods March 16, 2016 13 / 36 Solution Exact solution (found by “guess and verify”): zt a ct = (1 ab) e k t zt a kt+1 = abe kt Note that this solution is the same than in the model without stochastic volatility. Intuition. However, the dynamics of the model is affected by the law of motion for zt . Jesús Fernández-Villaverde (PENN) Solution Methods March 16, 2016 14 / 36 Steady state Steady state is also easy to find: 1 1 a k = (ab) a 1 c = (ab) 1 a (ab) 1 a z = 0 log s = log s Steady state in more general models. Jesús Fernández-Villaverde (PENN) Solution Methods March 16, 2016 15 / 36 The goal Define xt = (kt , zt 1, log st 1, #t , ut ; l). Role of l: perturbation parameter. We are searching for decision rules: c = c (x ) d = t t kt+1 = k (xt ) Then, we have a system of functional equations: 1 1 zt+1 a 1 = bEt ae k (xt ) c (xt ) c (xt+1) zt a c (xt ) + k (xt ) = e kt zt = rzt 1 + st l#t 1 2 2 log st = (1 r ) log s + r log st 1 + 1 r hlut s s s Jesús Fernández-Villaverde (PENN) Solution Methods March 16, 2016 16 / 36 Taylor’stheorem We will build a local approximation around x = (k, 0, log s, 0, 0; l). Given equilibrium conditions: 1 1 zt+1 a 1 Et b ae k (xt ) = 0 c (xt ) c (xt+1) zt a c (xt ) + k (xt ) e k = 0 t We will take derivatives with respect to (k, z, log s, #, u; l) and evaluate them around x = (k, 0, log s, 0, 0; l). Why? Apply Taylor’stheorem and a version of the implicit-function theorem. Jesús Fernández-Villaverde (PENN) Solution Methods March 16, 2016 17 / 36 Taylor series expansion I ct = c (kt , zt 1, log st 1, #t , ut ; 1) k,0,0 = c (x) j +ck (x)(kt k) + cz (x) zt 1 + clog s (x) log st 1 + +c# (x) #t + cu (x) ut + cl (x) 1 2 1 + ckk (x)(kt k) + ckz (x)(kt k) zt 1 + 2 2 1 + ck log s (x)(kt k) log st 1 + ... 2 Jesús Fernández-Villaverde (PENN) Solution Methods March 16, 2016 18 / 36 Taylor series expansion II kt = k (kt , zt 1, log st 1, #t , ut ; 1) k,0,0 = k (x) j +kk (x)(kt k) + kz (x) zt 1 + klog s (x) log st 1 + +k# (x) #t + ku (x) ut + kl (x) 1 2 1 + kkk (x)(kt k) + kkz (x)(kt k) zt 1 + 2 2 1 + kk log s (x)(kt k) log st 1 + ... 2 Jesús Fernández-Villaverde (PENN) Solution Methods March 16, 2016 19 / 36 Comment on notation From now on, to save on notation, I will write 1 1 z a 1 b ae t+1 k (xt ) 0 F (x ) = E c(xt ) c(xt+1 ) = t t z a c (xt ) + k (xt ) e t k 0 " t # I will take partial derivatives of F (xt ) and evaluate them at the steady state Do these derivatives exist? Jesús Fernández-Villaverde (PENN) Solution Methods March 16, 2016 20 / 36 First-order approximation We take first-order derivatives of F (xt ) around x. Thus, we get: DF (x) = 0 A matrix quadratic system. Why quadratic? Stable and unstable manifold. However, it has a nice recursive structure that we can and should exploit.
Recommended publications
  • 1 Probelms on Implicit Differentiation 2 Problems on Local Linearization
    Math-124 Calculus Chapter 3 Review 1 Probelms on Implicit Di®erentiation #1. A function is de¯ned implicitly by x3y ¡ 3xy3 = 3x + 4y + 5: Find y0 in terms of x and y. In problems 2-6, ¯nd the equation of the tangent line to the curve at the given point. x3 + 1 #2. + 2y2 = 1 ¡ 2x + 4y at the point (2; ¡1). y 1 #3. 4ey + 3x = + (y + 1)2 + 5x at the point (1; 0). x #4. (3x ¡ 2y)2 + x3 = y3 ¡ 2x ¡ 4 at the point (1; 2). p #5. xy + x3 = y3=2 ¡ y ¡ x at the point (1; 4). 2 #6. x sin(y ¡ 3) + 2y = 4x3 + at the point (1; 3). x #7. Find y00 for the curve xy + 2y3 = x3 ¡ 22y at the point (3; 1): #8. Find the points at which the curve x3y3 = x + y has a horizontal tangent. 2 Problems on Local Linearization #1. Let f(x) = (x + 2)ex. Find the value of f(0). Use this to approximate f(¡:2). #2. f(2) = 4 and f 0(2) = 7. Use linear approximation to approximate f(2:03). 6x4 #3. f(1) = 9 and f 0(x) = : Use a linear approximation to approximate f(1:02). x2 + 1 #4. A linear approximation is used to approximate y = f(x) at the point (3; 1). When ¢x = :06 and ¢y = :72. Find the equation of the tangent line. 3 Problems on Absolute Maxima and Minima 1 #1. For the function f(x) = x3 ¡ x2 ¡ 8x + 1, ¯nd the x-coordinates of the absolute max and 3 absolute min on the interval ² a) ¡3 · x · 5 ² b) 0 · x · 5 3 #2.
    [Show full text]
  • DYNAMICAL SYSTEMS Contents 1. Introduction 1 2. Linear Systems 5 3
    DYNAMICAL SYSTEMS WILLY HU Contents 1. Introduction 1 2. Linear Systems 5 3. Non-linear systems in the plane 8 3.1. The Linearization Theorem 11 3.2. Stability 11 4. Applications 13 4.1. A Model of Animal Conflict 13 4.2. Bifurcations 14 Acknowledgments 15 References 15 Abstract. This paper seeks to establish the foundation for examining dy- namical systems. Dynamical systems are, very broadly, systems that can be modelled by systems of differential equations. In this paper, we will see how to examine the qualitative structure of a system of differential equations and how to model it geometrically, and what information can be gained from such an analysis. We will see what it means for focal points to be stable and unstable, and how we can apply this to examining population growth and evolution, bifurcations, and other applications. 1. Introduction This paper is based on Arrowsmith and Place's book, Dynamical Systems.I have included corresponding references for propositions, theorems, and definitions. The images included in this paper are also from their book. Definition 1.1. (Arrowsmith and Place 1.1.1) Let X(t; x) be a real-valued function of the real variables t and x, with domain D ⊆ R2. A function x(t), with t in some open interval I ⊆ R, which satisfies dx (1.2) x0(t) = = X(t; x(t)) dt is said to be a solution satisfying x0. In other words, x(t) is only a solution if (t; x(t)) ⊆ D for each t 2 I. We take I to be the largest interval for which x(t) satisfies (1.2).
    [Show full text]
  • Appendix a Short Course in Taylor Series
    Appendix A Short Course in Taylor Series The Taylor series is mainly used for approximating functions when one can identify a small parameter. Expansion techniques are useful for many applications in physics, sometimes in unexpected ways. A.1 Taylor Series Expansions and Approximations In mathematics, the Taylor series is a representation of a function as an infinite sum of terms calculated from the values of its derivatives at a single point. It is named after the English mathematician Brook Taylor. If the series is centered at zero, the series is also called a Maclaurin series, named after the Scottish mathematician Colin Maclaurin. It is common practice to use a finite number of terms of the series to approximate a function. The Taylor series may be regarded as the limit of the Taylor polynomials. A.2 Definition A Taylor series is a series expansion of a function about a point. A one-dimensional Taylor series is an expansion of a real function f(x) about a point x ¼ a is given by; f 00ðÞa f 3ðÞa fxðÞ¼faðÞþf 0ðÞa ðÞþx À a ðÞx À a 2 þ ðÞx À a 3 þÁÁÁ 2! 3! f ðÞn ðÞa þ ðÞx À a n þÁÁÁ ðA:1Þ n! © Springer International Publishing Switzerland 2016 415 B. Zohuri, Directed Energy Weapons, DOI 10.1007/978-3-319-31289-7 416 Appendix A: Short Course in Taylor Series If a ¼ 0, the expansion is known as a Maclaurin Series. Equation A.1 can be written in the more compact sigma notation as follows: X1 f ðÞn ðÞa ðÞx À a n ðA:2Þ n! n¼0 where n ! is mathematical notation for factorial n and f(n)(a) denotes the n th derivation of function f evaluated at the point a.
    [Show full text]
  • Linearization of Nonlinear Differential Equation by Taylor's Series
    International Journal of Theoretical and Applied Science 4(1): 36-38(2011) ISSN No. (Print) : 0975-1718 International Journal of Theoretical & Applied Sciences, 1(1): 25-31(2009) ISSN No. (Online) : 2249-3247 Linearization of Nonlinear Differential Equation by Taylor’s Series Expansion and Use of Jacobian Linearization Process M. Ravi Tailor* and P.H. Bhathawala** *Department of Mathematics, Vidhyadeep Institute of Management and Technology, Anita, Kim, India **S.S. Agrawal Institute of Management and Technology, Navsari, India (Received 11 March, 2012, Accepted 12 May, 2012) ABSTRACT : In this paper, we show how to perform linearization of systems described by nonlinear differential equations. The procedure introduced is based on the Taylor's series expansion and on knowledge of Jacobian linearization process. We develop linear differential equation by a specific point, called an equilibrium point. Keywords : Nonlinear differential equation, Equilibrium Points, Jacobian Linearization, Taylor's Series Expansion. I. INTRODUCTION δx = a δ x In order to linearize general nonlinear systems, we will This linear model is valid only near the equilibrium point. use the Taylor Series expansion of functions. Consider a function f(x) of a single variable x, and suppose that x is a II. EQUILIBRIUM POINTS point such that f( x ) = 0. In this case, the point x is called Consider a nonlinear differential equation an equilibrium point of the system x = f( x ), since we have x( t )= f [ x ( t ), u ( t )] ... (1) x = 0 when x= x (i.e., the system reaches an equilibrium n m n at x ). Recall that the Taylor Series expansion of f(x) around where f:.
    [Show full text]
  • Student Understanding of Taylor Series Expansions in Statistical Mechanics
    PHYSICAL REVIEW SPECIAL TOPICS - PHYSICS EDUCATION RESEARCH 9, 020110 (2013) Student understanding of Taylor series expansions in statistical mechanics Trevor I. Smith,1 John R. Thompson,2,3 and Donald B. Mountcastle2 1Department of Physics and Astronomy, Dickinson College, Carlisle, Pennsylvania 17013, USA 2Department of Physics and Astronomy, University of Maine, Orono, Maine 04469, USA 3Maine Center for Research in STEM Education, University of Maine, Orono, Maine 04469, USA (Received 2 March 2012; revised manuscript received 23 May 2013; published 8 August 2013) One goal of physics instruction is to have students learn to make physical meaning of specific mathematical expressions, concepts, and procedures in different physical settings. As part of research investigating student learning in statistical physics, we are developing curriculum materials that guide students through a derivation of the Boltzmann factor using a Taylor series expansion of entropy. Using results from written surveys, classroom observations, and both individual think-aloud and teaching interviews, we present evidence that many students can recognize and interpret series expansions, but they often lack fluency in creating and using a Taylor series appropriately, despite previous exposures in both calculus and physics courses. DOI: 10.1103/PhysRevSTPER.9.020110 PACS numbers: 01.40.Fk, 05.20.Ày, 02.30.Lt, 02.30.Mv students’ use of Taylor series expansion to derive and I. INTRODUCTION understand the Boltzmann factor as a component of proba- Over the past several decades, much work has been done bility in the canonical ensemble. While this is a narrow investigating ways in which students understand and learn topic, we feel that it serves as an excellent example of physics (see Refs.
    [Show full text]
  • Linearization Extreme Values
    Math 31A Discussion Notes Week 6 November 3 and 5, 2015 This week we'll review two of last week's lecture topics in preparation for the quiz. Linearization One immediate use we have for derivatives is local linear approximation. On small neighborhoods around a point, a differentiable function behaves linearly. That is, if we zoom in enough on a point on a curve, the curve will eventually look like a straight line. We can use this fact to approximate functions by their tangent lines. You've seen all of this in lecture, so we'll jump straight to the formula for local linear approximation. If f is differentiable at x = a and x is \close" to a, then f(x) ≈ L(x) = f(a) + f 0(a)(x − a): Example. Use local linear approximation to estimate the value of sin(47◦). (Solution) We know that f(x) := sin(x) is differentiable everywhere, and we know the value of sin(45◦). Since 47◦ is reasonably close to 45◦, this problem is ripe for local linear approximation. We know that f 0(x) = π cos(x), so f 0(45◦) = πp . Then 180 180 2 1 π 90 + π sin(47◦) ≈ sin(45◦) + f 0(45◦)(47 − 45) = p + 2 p = p ≈ 0:7318: 2 180 2 90 2 For comparison, Google says that sin(47◦) = 0:7314, so our estimate is pretty good. Considering the fact that most folks now have (extremely powerful) calculators in their pockets, the above example is a very inefficient way to compute sin(47◦). Local linear approximation is no longer especially useful for estimating particular values of functions, but it can still be a very useful tool.
    [Show full text]
  • Nonlinear Optics
    Nonlinear Optics: A very brief introduction to the phenomena Paul Tandy ADVANCED OPTICS (PHYS 545), Spring 2010 Dr. Mendes INTRO So far we have considered only the case where the polarization density is proportionate to the electric field. The electric field itself is generally supplied by a laser source. PE= ε0 χ Where P is the polarization density, E is the electric field, χ is the susceptibility, and ε 0 is the permittivity of free space. We can imagine that the electric field somehow induces a polarization in the material. As the field increases the polarization also in turn increases. But we can also imagine that this process does not continue indefinitely. As the field grows to be very large the field itself starts to become on the order of the electric fields between the atoms in the material (medium). As this occurs the material itself becomes slightly altered and the behavior is no longer linear in nature but instead starts to become weakly nonlinear. We can represent this as a Taylor series expansion model in the electric field. ∞ n 23 PAEAEAEAE==++∑ n 12 3+... n=1 This model has become so popular that the coefficient scaling has become standardized using the letters d and χ (3) as parameters corresponding the quadratic and cubic terms. Higher order terms would play a role too if the fields were extremely high. 2(3)3 PEdEE=++()εχ0 (2) ( 4χ ) + ... To get an idea of how large a typical field would have to be before we see a nonlinear effect, let us scale the above equation by the first order term.
    [Show full text]
  • Linearization Via the Lie Derivative ∗
    Electron. J. Diff. Eqns., Monograph 02, 2000 http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu or ejde.math.unt.edu (login: ftp) Linearization via the Lie Derivative ∗ Carmen Chicone & Richard Swanson Abstract The standard proof of the Grobman–Hartman linearization theorem for a flow at a hyperbolic rest point proceeds by first establishing the analogous result for hyperbolic fixed points of local diffeomorphisms. In this exposition we present a simple direct proof that avoids the discrete case altogether. We give new proofs for Hartman’s smoothness results: A 2 flow is 1 linearizable at a hyperbolic sink, and a 2 flow in the C C C plane is 1 linearizable at a hyperbolic rest point. Also, we formulate C and prove some new results on smooth linearization for special classes of quasi-linear vector fields where either the nonlinear part is restricted or additional conditions on the spectrum of the linear part (not related to resonance conditions) are imposed. Contents 1 Introduction 2 2 Continuous Conjugacy 4 3 Smooth Conjugacy 7 3.1 Hyperbolic Sinks . 10 3.1.1 Smooth Linearization on the Line . 32 3.2 Hyperbolic Saddles . 34 4 Linearization of Special Vector Fields 45 4.1 Special Vector Fields . 46 4.2 Saddles . 50 4.3 Infinitesimal Conjugacy and Fiber Contractions . 50 4.4 Sources and Sinks . 51 ∗Mathematics Subject Classifications: 34-02, 34C20, 37D05, 37G10. Key words: Smooth linearization, Lie derivative, Hartman, Grobman, hyperbolic rest point, fiber contraction, Dorroh smoothing. c 2000 Southwest Texas State University. Submitted November 14, 2000.
    [Show full text]
  • Linearization and Stability Analysis of Nonlinear Problems
    Rose-Hulman Undergraduate Mathematics Journal Volume 16 Issue 2 Article 5 Linearization and Stability Analysis of Nonlinear Problems Robert Morgan Wayne State University Follow this and additional works at: https://scholar.rose-hulman.edu/rhumj Recommended Citation Morgan, Robert (2015) "Linearization and Stability Analysis of Nonlinear Problems," Rose-Hulman Undergraduate Mathematics Journal: Vol. 16 : Iss. 2 , Article 5. Available at: https://scholar.rose-hulman.edu/rhumj/vol16/iss2/5 Rose- Hulman Undergraduate Mathematics Journal Linearization and Stability Analysis of Nonlinear Problems Robert Morgana Volume 16, No. 2, Fall 2015 Sponsored by Rose-Hulman Institute of Technology Department of Mathematics Terre Haute, IN 47803 Email: [email protected] a http://www.rose-hulman.edu/mathjournal Wayne State University, Detroit, MI Rose-Hulman Undergraduate Mathematics Journal Volume 16, No. 2, Fall 2015 Linearization and Stability Analysis of Nonlinear Problems Robert Morgan Abstract. The focus of this paper is on the use of linearization techniques and lin- ear differential equation theory to analyze nonlinear differential equations. Often, mathematical models of real-world phenomena are formulated in terms of systems of nonlinear differential equations, which can be difficult to solve explicitly. To overcome this barrier, we take a qualitative approach to the analysis of solutions to nonlinear systems by making phase portraits and using stability analysis. We demonstrate these techniques in the analysis of two systems of nonlinear differential equations. Both of these models are originally motivated by population models in biology when solutions are required to be non-negative, but the ODEs can be un- derstood outside of this traditional scope of population models.
    [Show full text]
  • 1 Introduction 2 Linearization
    I.2 Quadratic Eigenvalue Problems 1 Introduction The quadratic eigenvalue problem (QEP) is to find scalars λ and nonzero vectors u satisfying Q(λ)x = 0, (1.1) where Q(λ)= λ2M + λD + K, M, D and K are given n × n matrices. Sometimes, we are also interested in finding the left eigenvectors y: yH Q(λ) = 0. Note that Q(λ) has 2n eigenvalues λ. They are the roots of det[Q(λ)] = 0. 2 Linearization A common way to solve the QEP is to first linearize it to a linear eigenvalue problem. For example, let λu z = , u Then the QEP (1.1) is equivalent to the generalized eigenvalue problem Lc(λ)z = 0 (2.2) where M 0 D K L (λ)= λ + ≡ λG + C. c 0 I −I 0 Lc(λ) is called a companion form or a linearization of Q(λ). Definition 2.1. A matrix pencil L(λ)= λG + C is called a linearization of Q(λ) if Q(λ) 0 E(λ)L(λ)F (λ)= (2.3) 0 I for some unimodular matrices E(λ) and F (λ).1 For the pencil Lc(λ) in (2.2), the identity (2.3) holds with I λM + D λI I E(λ)= , F (λ)= . 0 −I I 0 There are various ways to linearize a quadratic eigenvalue problem. Some are preferred than others. For example if M, D and K are symmetric and K is nonsingular, then we can preserve the symmetry property and use the following linearization: M 0 D K L (λ)= λ + .
    [Show full text]
  • Fundamentals of Duct Acoustics
    Fundamentals of Duct Acoustics Sjoerd W. Rienstra Technische Universiteit Eindhoven 16 November 2015 Contents 1 Introduction 3 2 General Formulation 4 3 The Equations 8 3.1 AHierarchyofEquations ........................... 8 3.2 BoundaryConditions. Impedance.. 13 3.3 Non-dimensionalisation . 15 4 Uniform Medium, No Mean Flow 16 4.1 Hard-walled Cylindrical Ducts . 16 4.2 RectangularDucts ............................... 21 4.3 SoftWallModes ................................ 21 4.4 AttenuationofSound.............................. 24 5 Uniform Medium with Mean Flow 26 5.1 Hard-walled Cylindrical Ducts . 26 5.2 SoftWallandUniformMeanFlow . 29 6 Source Expansion 32 6.1 ModalAmplitudes ............................... 32 6.2 RotatingFan .................................. 32 6.3 Tyler and Sofrin Rule for Rotor-Stator Interaction . ..... 33 6.4 PointSourceinaLinedFlowDuct . 35 6.5 PointSourceinaDuctWall .......................... 38 7 Reflection and Transmission 40 7.1 A Discontinuity in Diameter . 40 7.2 OpenEndReflection .............................. 43 VKI - 1 - CONTENTS CONTENTS A Appendix 49 A.1 BesselFunctions ................................ 49 A.2 AnImportantComplexSquareRoot . 51 A.3 Myers’EnergyCorollary ............................ 52 VKI - 2 - 1. INTRODUCTION CONTENTS 1 Introduction In a duct of constant cross section, with a medium and boundary conditions independent of the axial position, the wave equation for time-harmonic perturbations may be solved by means of a series expansion in a particular family of self-similar solutions, called modes. They are related to the eigensolutions of a two-dimensional operator, that results from the wave equation, on a cross section of the duct. For the common situation of a uniform medium without flow, this operator is the well-known Laplace operator 2. For a non- uniform medium, and in particular with mean flow, the details become mo∇re complicated, but the concept of duct modes remains by and large the same1.
    [Show full text]
  • Fourier Series
    Chapter 10 Fourier Series 10.1 Periodic Functions and Orthogonality Relations The differential equation ′′ y + 2y = F cos !t models a mass-spring system with natural frequency with a pure cosine forcing function of frequency !. If 2 = !2 a particular solution is easily found by undetermined coefficients (or by∕ using Laplace transforms) to be F y = cos !t. p 2 !2 − If the forcing function is a linear combination of simple cosine functions, so that the differential equation is N ′′ 2 y + y = Fn cos !nt n=1 X 2 2 where = !n for any n, then, by linearity, a particular solution is obtained as a sum ∕ N F y (t)= n cos ! t. p 2 !2 n n=1 n X − This simple procedure can be extended to any function that can be repre- sented as a sum of cosine (and sine) functions, even if that summation is not a finite sum. It turns out that the functions that can be represented as sums in this form are very general, and include most of the periodic functions that are usually encountered in applications. 723 724 10 Fourier Series Periodic Functions A function f is said to be periodic with period p> 0 if f(t + p)= f(t) for all t in the domain of f. This means that the graph of f repeats in successive intervals of length p, as can be seen in the graph in Figure 10.1. y p 2p 3p 4p 5p Fig. 10.1 An example of a periodic function with period p.
    [Show full text]