Lecture 7. Please note. Additional tutorial. Please note that there is no lecture on Tuesday, 15 November 2011.
|
|
- Chrystal Harrison
- 5 years ago
- Views:
Transcription
1 Lecture 7 3 Ordinary differential equations (ODEs) (continued) 6 Linear equations of second order 7 Systems of differential equations Please note Please note that there is no lecture on Tuesday, 15 November Additional tutorial The additional tutorial will start on Friday, 11 November 2011, , Infotech seminar room. 1/21
2 Linear equations of second order XIV We consider a damped harmonic oscillator driven by a force F(t). The differential equation is mẍ + dẋ + cx = F(t), where m > 0 is the mass and d > 0 is the damping. We divide the equation by m and set 2α = d m, ω2 0 = c m and f (t) = 1 mf(t). Then ẍ + 2αẋ + ω 2 0x = f (t). We first consider the homogeneous equation with f (t) 0 (no external force). The characteristic polynomial is λ 2 + 2αλ + ω0 2 with zeroes λ 1,2 = α ± α 2 ω0 2. There are three cases depending on the sign of α 2 ω /21
3 Linear equations of second order XV 1st case: α 2 > ω0 2. In this case (strong damping) we have two real zeroes λ 1, λ 2, both of which are negative. The general solution is x(t) = c 1 e λ 1t + c 2 e λ 2t. The constants c 1 and c 2 can be determined from initial values for x(t) and ẋ(t). Since both λ 1, λ 2 are negative, we have x(t) 0 for large times t. There is no oscillation. 2nd case: α 2 = ω0 2. In this case (aperiodic limit case) we have one real zero λ 1,2 = α. The general solution of the ODE is x(t) = (c 1 + c 2 t)e αt. Again, there is no oscillation and x(t) 0 for large times t. 3/21
4 Linear equations of second order XVI 3rd case: α 2 < ω0 2. In this case (weak damping) we have two complex zeroes λ 1,2 = α ± iω 1, where ω 1 denotes the so-called eigenfrequency ω0 2 α2 > 0. The general solution is x(t) = e αt (c 1 cos ω 1 t + c 2 sinω 1 t). Again, we have x(t) 0 for large times t but this time there is an oscillation. We now consider the inhomogeneous equation with a harmonic driving force of frequency ω: ẍ + 2αẋ + ω 2 0x = B cos ωt. Since the solutions of the homogeneous equation all satisfy x h (t) 0 for large times t, in the end there will only remain the particular solution x p (t) of the inhomogeneous equation. 4/21
5 Linear equations of second order XVII We consider the equation with complexified right hand side ẍ + 2αẋ + ω 2 0x = Be iωt and use the ansatz x(t) = Ce iωt for the particular solution. This leads to C = The real part of x(t) is therefore x p (t) = B ω0 2 ω2 + 2iαω = B(ω2 0 ω2 2iαω) (ω0 2 ω2 ) 2 + 4α 2 ω 2. B (ω 2 0 ω2 ) 2 + 4ω 2 α 2((ω2 0 ω 2 )cos ωt + 2αω sinωt). Using trigonometric identities this can be written as x p (t) = B 2αω cos(ωt φ), with φ = arctan (ω 20 ω2 ) 2 + 4ω 2 α 2 ω0 2 ω2. 5/21
6 Linear equations of second order XVIII The factor A(ω) = 1 (ω0 2 ω2 ) 2 + 4ω 2 α 2 is called amplification factor. Note that A(ω) 0 for large ω. Hence driving forces with a very high frequency are of no effect. The amplification factor A(ω) has a maximum if the denominator has a minimum. This happens precisely for 4(ω 2 0 ω 2 )ω + 8ωα 2 = 0, hence for the frequency ω = ω r, where ω 2 r = ω 2 0 2α 2 = ω 2 1 α 2 and ω 1 is the eigenfrequency. The frequency ω r is called resonance frequency and exists only for damping constant α < 1 2 ω 0. 6/21
7 Systems of differential equations I It often happens in applications that several scalar functions x i (t) are coupled in the sense that the first derivative of one function depends on the variable t and on the value of all other functions: ẋ i = v i (t, x 1,...,x n ), 1 i n, where v i are certain functions of n + 1 variables. These n equations form a system of ordinary differential equations. It is useful to write this system in the following form: Let x denote the time dependent column vector (x 1,...,x n ) and v the column vector (v 1,...,v n ). Then we can write x = v(t, x). A solution of this system of ODEs is a map x : I R n, where I R is an interval, such that x(t) = v(t, x(t)). If we demand in addition that x(t 0 ) = x 0 we have an initial value problem. 7/21
8 Systems of differential equations II If we have a system of ODEs x = v(t, x), we can think of v as a time dependent vector field on a subset of R n. A solution x(t) is a curve in this subset of R n which is at each instance of time t tangential to this vector field. Example Consider the non-linear system of equations ẋ = 3y 3x ẏ = (1 + a 2 )x y xz ż = xy z, where a is a constant. It is a simplified model of atmospheric turbulence. Depending on the value of a this system admits so-called chaotic solutions related to the Lorenz attractor. 8/21
9 Systems of differential equations III Note that any explicit differential equation of order n for a scalar function x(t) x (n) = f (t, x, ẋ,...,x (n 1) ) can be transformed into an equivalent system of ODEs of first order ẋ 0 = x 1, ẋ 1 = x 2,. ẋ n 2 = x n 1, ẋ n 1 = f (t, x 0, x 1,...,x n 1 ). The scalar function x(t) is a solution of the n-th order explicit ODE if and only if x(t) = (x(t), ẋ(t),...,x (n 1) (t)) T is a solution of the system of ODEs. 9/21
10 Systems of differential equations IV There are certain fundamental theorems which guarantee the existence and uniqueness of solutions to initial value problems of first order ODEs or systems of ODEs. We will only formulate the general version for first order systems of ODEs. Theorem (Peano) Suppose that the time dependent vector field v : G R n, with G R n+1, is continuous. Then the initial value problem x = v(t, x), x(t 0 ) = x 0, with (t 0, x 0 ) G, has at least one solution. Therefore continuity of v guarantees the existence of a solution. However, it is not sufficient to guarantee the uniqueness. 10/21
11 Systems of differential equations V Example Consider the initial value problem ẋ = 2 x, x(0) = 0. This has the constant solution x 1 (t) 0. However, it also has the solution { t 2 if t 0 x 2 (t) = t 2 if t < 0. Hence the solution is not unique even though the function 2 x is continuous. 11/21
12 Systems of differential equations VI Let G R n+1 be some subset and v : G R n a vector field on G. Definition We say that v satisfies a Lipschitz condition on G if there exists a constant L 0 such that for all (t, x) and (t, x ) in G. v(t, x) v(t, x ) L x x, Here the norm of a vector is defined by x = x x2 n. 12/21
13 Systems of differential equations VII We can now formulate a uniqueness theorem for initial value problems: Theorem (Picard-Lindelöf) Suppose that the vector field v(t, x) satisfies a Lipschitz condition on [a, b] R n. Then the initial value problem x = v(t, x), x(t 0 ) = x 0, for each x 0 R n has a unique solution defined on the interval I = [a, b]. 13/21
14 Systems of differential equations VIII We will only consider systems of the following form: x = A x + b(t), where A R n n is a constant matrix. These systems are called linear systems with constant coefficients. If b(t) 0, the system is called homogeneous. Theorem The space of solutions of the homogeneous system with A R n n is an n-dimensional vector space. Let x p (t) denote some particular solution to the inhomogeneous equation. Then any solution of this equation can be written as x(t) = x p (t) + x h (t), where x h (t) is a solution of the homogeneous equation. 14/21
15 Systems of differential equations IX To describe the general solution of a homogeneous linear system of equations we need the concept of the exponential function of a matrix. First, let p(x) = a 0 + a 1 x a n x n denote some polynomial in the variable x. Then for an n n-matrix we define p(a) = a 0 I n + a 1 A a n A n. Similarly we can define a power series in the square matrix A, where convergence means element-wise convergence. In particular we can define e A = k=0 1 k! Ak = I n + A A ! A and hence e ta = k=0 t k k! Ak = I n + ta + t2 2 A2 + t3 3! A /21
16 Systems of differential equations X The exponential function of a square matrix has the following properties: Exponential function of a square matrix Special case: e 0 = I n, where 0 denotes the n n zero-matrix. Functional equation: If AB = BA, then e A e B = e B e A = e A+B. Inverse matrix: The matrix e A is always invertible with inverse (e A ) 1 = e A (this follows from the first two properties). Derivative: d dt eta = Ae ta = e ta A. Note that the second property only holds in general if A and B commute. 16/21
17 Systems of differential equations XI Theorem Every solution of a homogeneous linear system of ODEs with constant coefficients x = A x is of the form x(t) = e ta c, defined for all t R, where c is any column vector in R n. Proof. Note that x(t) = e ta c is a solution of the differential equation, because x(t) = Ae ta c = A x(t). Let c 1,..., c n denote a basis of R n and x 1,..., x n the corresponding solutions. Then the solutions are linearly independent at each time t, hence the space spanned by these solutions is n-dimensional. Since the space of solutions is also n-dimensional, every solution is of this form for some vector c. 17/21
18 Systems of differential equations XII In general the exponential e ta is difficult to calculate if A is any given matrix. Note the following: Lemma If u is an eigenvector of A with eigenvalue λ, then e ta u = e tλ u. The exponential on the right hand side is simply the exponential of a number. Proof. We have e ta u = (I n + ta + t2 2 A2 + t3 3! A3 +...) u = u + tλ u + t2 2 λ2 u + t3 3! λ3 u... = e tλ u. 18/21
19 Systems of differential equations XIII The lemma implies: Theorem Suppose that the matrix A R n n has a basis of (complex) eigenvectors u 1,..., u n with (complex) eigenvalues λ 1,...,λ n. Then the time dependent vectors e tλ 1 u 1,...,e tλn u n form a basis of the solution space of the linear homogeneous system x = A x. Note that a basis of (real) eigenvectors exist in particular if the matrix A is symmetric. Theorem If λ is a non-real complex eigenvalue with complex eigenvector u C n, then both Re(e tλ u) and Im(e tλ u) are linearly independent real solutions of the system x = A x. 19/21
12.1. Exponential shift. The calculation (10.1)
62 12. Resonance and the exponential shift law 12.1. Exponential shift. The calculation (10.1) (1) p(d)e = p(r)e extends to a formula for the effect of the operator p(d) on a product of the form e u, where
More information4. Complex Oscillations
4. Complex Oscillations The most common use of complex numbers in physics is for analyzing oscillations and waves. We will illustrate this with a simple but crucially important model, the damped harmonic
More informationChapter 2: Complex numbers
Chapter 2: Complex numbers Complex numbers are commonplace in physics and engineering. In particular, complex numbers enable us to simplify equations and/or more easily find solutions to equations. We
More informationMath 211. Substitute Lecture. November 20, 2000
1 Math 211 Substitute Lecture November 20, 2000 2 Solutions to y + py + qy =0. Look for exponential solutions y(t) =e λt. Characteristic equation: λ 2 + pλ + q =0. Characteristic polynomial: λ 2 + pλ +
More information10. Operators and the Exponential Response Formula
52 10. Operators and the Exponential Response Formula 10.1. Operators. Operators are to functions as functions are to numbers. An operator takes a function, does something to it, and returns this modified
More informationExercises Lecture 15
AM1 Mathematical Analysis 1 Oct. 011 Feb. 01 Date: January 7 Exercises Lecture 15 Harmonic Oscillators In classical mechanics, a harmonic oscillator is a system that, when displaced from its equilibrium
More informationSolutions to Math 53 Math 53 Practice Final
Solutions to Math 5 Math 5 Practice Final 20 points Consider the initial value problem y t 4yt = te t with y 0 = and y0 = 0 a 8 points Find the Laplace transform of the solution of this IVP b 8 points
More informationAPPM 2360: Final Exam 10:30am 1:00pm, May 6, 2015.
APPM 23: Final Exam :3am :pm, May, 25. ON THE FRONT OF YOUR BLUEBOOK write: ) your name, 2) your student ID number, 3) lecture section, 4) your instructor s name, and 5) a grading table for eight questions.
More informationLinear Differential Equations. Problems
Chapter 1 Linear Differential Equations. Problems 1.1 Introduction 1.1.1 Show that the function ϕ : R R, given by the expression ϕ(t) = 2e 3t for all t R, is a solution of the Initial Value Problem x =
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationMathematics for Engineers II. lectures. Differential Equations
Differential Equations Examples for differential equations Newton s second law for a point mass Consider a particle of mass m subject to net force a F. Newton s second law states that the vector acceleration
More informationMATH 23 Exam 2 Review Solutions
MATH 23 Exam 2 Review Solutions Problem 1. Use the method of reduction of order to find a second solution of the given differential equation x 2 y (x 0.1875)y = 0, x > 0, y 1 (x) = x 1/4 e 2 x Solution
More informationMATH 24 EXAM 3 SOLUTIONS
MATH 4 EXAM 3 S Consider the equation y + ω y = cosω t (a) Find the general solution of the homogeneous equation (b) Find the particular solution of the non-homogeneous equation using the method of Undetermined
More informationTest #2 Math 2250 Summer 2003
Test #2 Math 225 Summer 23 Name: Score: There are six problems on the front and back of the pages. Each subpart is worth 5 points. Show all of your work where appropriate for full credit. ) Show the following
More informationAnalytical Mechanics ( AM )
Analytical Mechanics ( AM ) Olaf Scholten KVI, kamer v8; tel nr 6-55; email: scholten@kvinl Web page: http://wwwkvinl/ scholten Book: Classical Dynamics of Particles and Systems, Stephen T Thornton & Jerry
More informationLinear second-order differential equations with constant coefficients and nonzero right-hand side
Linear second-order differential equations with constant coefficients and nonzero right-hand side We return to the damped, driven simple harmonic oscillator d 2 y dy + 2b dt2 dt + ω2 0y = F sin ωt We note
More informationLinear Algebra and ODEs review
Linear Algebra and ODEs review Ania A Baetica September 9, 015 1 Linear Algebra 11 Eigenvalues and eigenvectors Consider the square matrix A R n n (v, λ are an (eigenvector, eigenvalue pair of matrix A
More informationODE Homework 1. Due Wed. 19 August 2009; At the beginning of the class
ODE Homework Due Wed. 9 August 2009; At the beginning of the class. (a) Solve Lẏ + Ry = E sin(ωt) with y(0) = k () L, R, E, ω are positive constants. (b) What is the limit of the solution as ω 0? (c) Is
More informationP441 Analytical Mechanics - I. RLC Circuits. c Alex R. Dzierba. In this note we discuss electrical oscillating circuits: undamped, damped and driven.
Lecture 10 Monday - September 19, 005 Written or last updated: September 19, 005 P441 Analytical Mechanics - I RLC Circuits c Alex R. Dzierba Introduction In this note we discuss electrical oscillating
More informationMath 216 Second Midterm 16 November, 2017
Math 216 Second Midterm 16 November, 2017 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material
More informationThe Harmonic Oscillator
The Harmonic Oscillator Math 4: Ordinary Differential Equations Chris Meyer May 3, 008 Introduction The harmonic oscillator is a common model used in physics because of the wide range of problems it can
More informationTransform Solutions to LTI Systems Part 3
Transform Solutions to LTI Systems Part 3 Example of second order system solution: Same example with increased damping: k=5 N/m, b=6 Ns/m, F=2 N, m=1 Kg Given x(0) = 0, x (0) = 0, find x(t). The revised
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear
More informationDIAGONALIZABLE LINEAR SYSTEMS AND STABILITY. 1. Algebraic facts. We first recall two descriptions of matrix multiplication.
DIAGONALIZABLE LINEAR SYSTEMS AND STABILITY. 1. Algebraic facts. We first recall two descriptions of matrix multiplication. Let A be n n, P be n r, given by its columns: P = [v 1 v 2... v r ], where the
More informationThe goal of this chapter is to study linear systems of ordinary differential equations: dt,..., dx ) T
1 1 Linear Systems The goal of this chapter is to study linear systems of ordinary differential equations: ẋ = Ax, x(0) = x 0, (1) where x R n, A is an n n matrix and ẋ = dx ( dt = dx1 dt,..., dx ) T n.
More informationMath Assignment 5
Math 2280 - Assignment 5 Dylan Zwick Fall 2013 Section 3.4-1, 5, 18, 21 Section 3.5-1, 11, 23, 28, 35, 47, 56 Section 3.6-1, 2, 9, 17, 24 1 Section 3.4 - Mechanical Vibrations 3.4.1 - Determine the period
More informationClassical Mechanics Phys105A, Winter 2007
Classical Mechanics Phys5A, Winter 7 Wim van Dam Room 59, Harold Frank Hall vandam@cs.ucsb.edu http://www.cs.ucsb.edu/~vandam/ Phys5A, Winter 7, Wim van Dam, UCSB Midterm New homework has been announced
More informationIntroductory Physics. Week 2015/05/29
2015/05/29 Part I Summary of week 6 Summary of week 6 We studied the motion of a projectile under uniform gravity, and constrained rectilinear motion, introducing the concept of constraint force. Then
More informationNo books, no notes, no calculators. You must show work, unless the question is a true/false, yes/no, or fill-in-the-blank question.
Math 304 Final Exam (May 8) Spring 206 No books, no notes, no calculators. You must show work, unless the question is a true/false, yes/no, or fill-in-the-blank question. Name: Section: Question Points
More informationSolutions for homework 5
1 Section 4.3 Solutions for homework 5 17. The following equation has repeated, real, characteristic roots. Find the general solution. y 4y + 4y = 0. The characteristic equation is λ 4λ + 4 = 0 which has
More informationThis work is licensed under a Creative Commons Attribution-Noncommercial-Share Alike 4.0 License.
University of Rhode Island DigitalCommons@URI Classical Dynamics Physics Course Materials 015 14. Oscillations Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative Commons License This wor
More information93 Analytical solution of differential equations
1 93 Analytical solution of differential equations 1. Nonlinear differential equation The only kind of nonlinear differential equations that we solve analytically is the so-called separable differential
More informationFourier transforms, Generalised functions and Greens functions
Fourier transforms, Generalised functions and Greens functions T. Johnson 2015-01-23 Electromagnetic Processes In Dispersive Media, Lecture 2 - T. Johnson 1 Motivation A big part of this course concerns
More informationDifferential Equations Grinshpan Two-Dimensional Homogeneous Linear Systems with Constant Coefficients. Purely Imaginary Eigenvalues. Recall the equation mẍ+k = of a simple harmonic oscillator with frequency
More informationControl Systems. Dynamic response in the time domain. L. Lanari
Control Systems Dynamic response in the time domain L. Lanari outline A diagonalizable - real eigenvalues (aperiodic natural modes) - complex conjugate eigenvalues (pseudoperiodic natural modes) - phase
More informationDamped Oscillation Solution
Lecture 19 (Chapter 7): Energy Damping, s 1 OverDamped Oscillation Solution Damped Oscillation Solution The last case has β 2 ω 2 0 > 0. In this case we define another real frequency ω 2 = β 2 ω 2 0. In
More information18.03SC Practice Problems 14
1.03SC Practice Problems 1 Frequency response Solution suggestions In this problem session we will work with a second order mass-spring-dashpot system driven by a force F ext acting directly on the mass:
More informationLinear and Nonlinear Oscillators (Lecture 2)
Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical
More information3.4.1 Distinct Real Roots
Math 334 3.4. CONSTANT COEFFICIENT EQUATIONS 34 Assume that P(x) = p (const.) and Q(x) = q (const.), so that we have y + py + qy = 0, or L[y] = 0 (where L := d2 dx 2 + p d + q). (3.7) dx Guess a solution
More informationSection 3.4. Second Order Nonhomogeneous. The corresponding homogeneous equation
Section 3.4. Second Order Nonhomogeneous Equations y + p(x)y + q(x)y = f(x) (N) The corresponding homogeneous equation y + p(x)y + q(x)y = 0 (H) is called the reduced equation of (N). 1 General Results
More informationAutonomous system = system without inputs
Autonomous system = system without inputs State space representation B(A,C) = {y there is x, such that σx = Ax, y = Cx } x is the state, n := dim(x) is the state dimension, y is the output Polynomial representation
More informationPhysics 351 Monday, January 22, 2018
Physics 351 Monday, January 22, 2018 Phys 351 Work on this while you wait for your classmates to arrive: Show that the moment of inertia of a uniform solid sphere rotating about a diameter is I = 2 5 MR2.
More informationNotes on the Periodically Forced Harmonic Oscillator
Notes on the Periodically orced Harmonic Oscillator Warren Weckesser Math 38 - Differential Equations 1 The Periodically orced Harmonic Oscillator. By periodically forced harmonic oscillator, we mean the
More informationDifferential Equations
Differential Equations A differential equation (DE) is an equation which involves an unknown function f (x) as well as some of its derivatives. To solve a differential equation means to find the unknown
More informationManifesto on Numerical Integration of Equations of Motion Using Matlab
Manifesto on Numerical Integration of Equations of Motion Using Matlab C. Hall April 11, 2002 This handout is intended to help you understand numerical integration and to put it into practice using Matlab
More informationClasses of Linear Operators Vol. I
Classes of Linear Operators Vol. I Israel Gohberg Seymour Goldberg Marinus A. Kaashoek Birkhäuser Verlag Basel Boston Berlin TABLE OF CONTENTS VOLUME I Preface Table of Contents of Volume I Table of Contents
More informationLab #2 - Two Degrees-of-Freedom Oscillator
Lab #2 - Two Degrees-of-Freedom Oscillator Last Updated: March 0, 2007 INTRODUCTION The system illustrated in Figure has two degrees-of-freedom. This means that two is the minimum number of coordinates
More informationVibrations and waves: revision. Martin Dove Queen Mary University of London
Vibrations and waves: revision Martin Dove Queen Mary University of London Form of the examination Part A = 50%, 10 short questions, no options Part B = 50%, Answer questions from a choice of 4 Total exam
More informationWave Phenomena Physics 15c. Lecture 2 Damped Oscillators Driven Oscillators
Wave Phenomena Physics 15c Lecture Damped Oscillators Driven Oscillators What We Did Last Time Analyzed a simple harmonic oscillator The equation of motion: The general solution: Studied the solution m
More informationSection 3.4. Second Order Nonhomogeneous. The corresponding homogeneous equation. is called the reduced equation of (N).
Section 3.4. Second Order Nonhomogeneous Equations y + p(x)y + q(x)y = f(x) (N) The corresponding homogeneous equation y + p(x)y + q(x)y = 0 (H) is called the reduced equation of (N). 1 General Results
More informationDifferential Equations 2280 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 2015 at 12:50pm
Differential Equations 228 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 25 at 2:5pm Instructions: This in-class exam is 5 minutes. No calculators, notes, tables or books. No answer check is
More informationPhysics III: Final Solutions (Individual and Group)
Massachusetts Institute of Technology MITES 7 Physics III First and Last Name: Physics III: Final Solutions (Individual and Group Instructor: Mobolaji Williams Teaching Assistant: Rene García (Tuesday
More informationStudy guide - Math 220
Study guide - Math 220 November 28, 2012 1 Exam I 1.1 Linear Equations An equation is linear, if in the form y + p(t)y = q(t). Introducing the integrating factor µ(t) = e p(t)dt the solutions is then in
More informationL = 1 2 a(q) q2 V (q).
Physics 3550, Fall 2011 Motion near equilibrium - Small Oscillations Relevant Sections in Text: 5.1 5.6 Motion near equilibrium 1 degree of freedom One of the most important situations in physics is motion
More informationMATH 213 Linear Algebra and Ordinary Differential Equations Spring 2015 Study Sheet for Final Exam. Topics
MATH 213 Linear Algebra and Ordinary Differential Equations Spring 2015 Study Sheet for Final Exam This study sheet will not be allowed during the test. Books and notes will not be allowed during the test.
More informationMATH 53H - Solutions to Problem Set III
MATH 53H - Solutions to Problem Set III. Fix λ not an eigenvalue of L. Then det(λi L) 0 λi L is invertible. We have then p L+N (λ) = det(λi L N) = det(λi L) det(i (λi L) N) = = p L (λ) det(i (λi L) N)
More informationSolution of Linear State-space Systems
Solution of Linear State-space Systems Homogeneous (u=0) LTV systems first Theorem (Peano-Baker series) The unique solution to x(t) = (t, )x 0 where The matrix function is given by is called the state
More informationVibrations: Second Order Systems with One Degree of Freedom, Free Response
Single Degree of Freedom System 1.003J/1.053J Dynamics and Control I, Spring 007 Professor Thomas Peacock 5//007 Lecture 0 Vibrations: Second Order Systems with One Degree of Freedom, Free Response Single
More informationSelected Topics in Physics a lecture course for 1st year students by W.B. von Schlippe Spring Semester 2007
Selected Topics in Physics a lecture course for st year students by W.B. von Schlippe Spring Semester 7 Lecture : Oscillations simple harmonic oscillations; coupled oscillations; beats; damped oscillations;
More informationUnderstand the existence and uniqueness theorems and what they tell you about solutions to initial value problems.
Review Outline To review for the final, look over the following outline and look at problems from the book and on the old exam s and exam reviews to find problems about each of the following topics.. Basics
More informationThe formulas for derivatives are particularly useful because they reduce ODEs to algebraic expressions. Consider the following ODE d 2 dx + p d
Solving ODEs using Fourier Transforms The formulas for derivatives are particularly useful because they reduce ODEs to algebraic expressions. Consider the following ODE d 2 dx + p d 2 dx + q f (x) R(x)
More information1 (2n)! (-1)n (θ) 2n
Complex Numbers and Algebra The real numbers are complete for the operations addition, subtraction, multiplication, and division, or more suggestively, for the operations of addition and multiplication
More informationACM/CMS 107 Linear Analysis & Applications Fall 2016 Assignment 4: Linear ODEs and Control Theory Due: 5th December 2016
ACM/CMS 17 Linear Analysis & Applications Fall 216 Assignment 4: Linear ODEs and Control Theory Due: 5th December 216 Introduction Systems of ordinary differential equations (ODEs) can be used to describe
More informationNormal Modes, Wave Motion and the Wave Equation Hilary Term 2011 Lecturer: F Hautmann
Normal Modes, Wave Motion and the Wave Equation Hilary Term 2011 Lecturer: F Hautmann Part A: Normal modes ( 4 lectures) Part B: Waves ( 8 lectures) printed lecture notes slides will be posted on lecture
More informationProf. Krstic Nonlinear Systems MAE281A Homework set 1 Linearization & phase portrait
Prof. Krstic Nonlinear Systems MAE28A Homework set Linearization & phase portrait. For each of the following systems, find all equilibrium points and determine the type of each isolated equilibrium. Use
More informationOld Math 330 Exams. David M. McClendon. Department of Mathematics Ferris State University
Old Math 330 Exams David M. McClendon Department of Mathematics Ferris State University Last updated to include exams from Fall 07 Contents Contents General information about these exams 3 Exams from Fall
More informationSeminar 6: COUPLED HARMONIC OSCILLATORS
Seminar 6: COUPLED HARMONIC OSCILLATORS 1. Lagrangian Equations of Motion Let consider a system consisting of two harmonic oscillators that are coupled together. As a model, we will use two particles attached
More informationdx n a 1(x) dy
HIGHER ORDER DIFFERENTIAL EQUATIONS Theory of linear equations Initial-value and boundary-value problem nth-order initial value problem is Solve: a n (x) dn y dx n + a n 1(x) dn 1 y dx n 1 +... + a 1(x)
More information1 Controller Optimization according to the Modulus Optimum
Controller Optimization according to the Modulus Optimum w G K (s) F 0 (s) x The goal of applying a control loop usually is to get the control value x equal to the reference value w. x(t) w(t) X(s) W (s)
More informationPhysics 141, Lecture 7. Outline. Course Information. Course information: Homework set # 3 Exam # 1. Quiz. Continuation of the discussion of Chapter 4.
Physics 141, Lecture 7. Frank L. H. Wolfs Department of Physics and Astronomy, University of Rochester, Lecture 07, Page 1 Outline. Course information: Homework set # 3 Exam # 1 Quiz. Continuation of the
More information( ) f (k) = FT (R(x)) = R(k)
Solving ODEs using Fourier Transforms The formulas for derivatives are particularly useful because they reduce ODEs to algebraic expressions. Consider the following ODE d 2 dx + p d 2 dx + q f (x) = R(x)
More informationELECTROMAGNETIC WAVES
Physics 4D ELECTROMAGNETIC WAVE Hans P. Paar 26 January 2006 i Chapter 1 Vector Calculus 1.1 Introduction Vector calculus is a branch of mathematics that allows differentiation and integration of (scalar)
More information221B Lecture Notes on Resonances in Classical Mechanics
1B Lecture Notes on Resonances in Classical Mechanics 1 Harmonic Oscillators Harmonic oscillators appear in many different contexts in classical mechanics. Examples include: spring, pendulum (with a small
More informationThe particular integral :
Second order linear equation with constant coefficients The particular integral : d f df Lf = a + a + af = h( x) Solutions with combinations of driving functions d f df Lf = a + a + af = h( x) + h( x)
More informationThe theory of harmonically driven, clamped silicon micro-beams
The theory of harmonically driven, clamped silicon micro-beams Henrik Bruus Mikroelektronik Centret, bldg. 345, Technical University of Denmark, DK-28 yngby (Dated: June 1, 22) In this note the theory
More information0.1 Rational Canonical Forms
We have already seen that it is useful and simpler to study linear systems using matrices. But matrices are themselves cumbersome, as they are stuffed with many entries, and it turns out that it s best
More informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016
Prof. Dr. Eleni Chatzi Lecture 4-09. March, 2016 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural
More informationORDINARY DIFFERENTIAL EQUATIONS FOR PHYSICISTS. Lecture Notes
ORDINARY DIFFERENTIAL EQUATIONS FOR PHYSICISTS Wolfgang König Lecture Notes University Leipzig summer term 25 Contents 1 Introduction 1 1.1 Motivation and basic notions............................. 1
More informationBackground ODEs (2A) Young Won Lim 3/7/15
Background ODEs (2A) Copyright (c) 2014-2015 Young W. Lim. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or any
More informationINHOMOGENEOUS LINEAR SYSTEMS OF DIFFERENTIAL EQUATIONS
INHOMOGENEOUS LINEAR SYSTEMS OF DIFFERENTIAL EQUATIONS Definitions and a general fact If A is an n n matrix and f(t) is some given vector function, then the system of differential equations () x (t) Ax(t)
More informationStabilization of persistently excited linear systems
Stabilization of persistently excited linear systems Yacine Chitour Laboratoire des signaux et systèmes & Université Paris-Sud, Orsay Exposé LJLL Paris, 28/9/2012 Stabilization & intermittent control Consider
More information1. Find the solution of the following uncontrolled linear system. 2 α 1 1
Appendix B Revision Problems 1. Find the solution of the following uncontrolled linear system 0 1 1 ẋ = x, x(0) =. 2 3 1 Class test, August 1998 2. Given the linear system described by 2 α 1 1 ẋ = x +
More informationProblem 1: Lagrangians and Conserved Quantities. Consider the following action for a particle of mass m moving in one dimension
105A Practice Final Solutions March 13, 01 William Kelly Problem 1: Lagrangians and Conserved Quantities Consider the following action for a particle of mass m moving in one dimension S = dtl = mc dt 1
More informationEntrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4.
Entrance Exam, Differential Equations April, 7 (Solve exactly 6 out of the 8 problems). Consider the following initial value problem: { y + y + y cos(x y) =, y() = y. Find all the values y such that the
More informationAnalysis of Dynamical Systems
2 YFX1520 Nonlinear phase portrait Mathematical Modelling and Nonlinear Dynamics Coursework Assignments 1 ϕ (t) 0-1 -2-6 -4-2 0 2 4 6 ϕ(t) Dmitri Kartofelev, PhD 2018 Variant 1 Part 1: Liénard type equation
More informationComplex Numbers. The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition,
Complex Numbers Complex Algebra The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition, and (ii) complex multiplication, (x 1, y 1 )
More informationFourier Series. Spectral Analysis of Periodic Signals
Fourier Series. Spectral Analysis of Periodic Signals he response of continuous-time linear invariant systems to the complex exponential with unitary magnitude response of a continuous-time LI system at
More informationEcon Slides from Lecture 7
Econ 205 Sobel Econ 205 - Slides from Lecture 7 Joel Sobel August 31, 2010 Linear Algebra: Main Theory A linear combination of a collection of vectors {x 1,..., x k } is a vector of the form k λ ix i for
More informationREAL LINEAR ALGEBRA: PROBLEMS WITH SOLUTIONS
REAL LINEAR ALGEBRA: PROBLEMS WITH SOLUTIONS The problems listed below are intended as review problems to do before the final They are organied in groups according to sections in my notes, but it is not
More informationAPPM 2360: Midterm exam 3 April 19, 2017
APPM 36: Midterm exam 3 April 19, 17 On the front of your Bluebook write: (1) your name, () your instructor s name, (3) your lecture section number and (4) a grading table. Text books, class notes, cell
More informationVANDERBILT UNIVERSITY. MATH 2610 ORDINARY DIFFERENTIAL EQUATIONS Practice for test 1 solutions
VANDERBILT UNIVERSITY MATH 2610 ORDINARY DIFFERENTIAL EQUATIONS Practice for test 1 solutions The first test will cover all material discussed up to (including) section 4.5. Important: The solutions below
More informationOrdinary differential equations
Class 7 Today s topics The nonhomogeneous equation Resonance u + pu + qu = g(t). The nonhomogeneous second order linear equation This is the nonhomogeneous second order linear equation u + pu + qu = g(t).
More information1 Pushing your Friend on a Swing
Massachusetts Institute of Technology MITES 017 Physics III Lecture 05: Driven Oscillations In these notes, we derive the properties of both an undamped and damped harmonic oscillator under the influence
More informationMTH 464: Computational Linear Algebra
MTH 464: Computational Linear Algebra Lecture Outlines Exam 4 Material Prof. M. Beauregard Department of Mathematics & Statistics Stephen F. Austin State University April 15, 2018 Linear Algebra (MTH 464)
More informationSolutions to homework assignment #3 Math 119B UC Davis, Spring = 12x. The Euler-Lagrange equation is. 2q (x) = 12x. q(x) = x 3 + x.
1. Find the stationary points of the following functionals: (a) 1 (q (x) + 1xq(x)) dx, q() =, q(1) = Solution. The functional is 1 L(q, q, x) where L(q, q, x) = q + 1xq. We have q = q, q = 1x. The Euler-Lagrange
More informationLinear Algebra, part 2 Eigenvalues, eigenvectors and least squares solutions
Linear Algebra, part 2 Eigenvalues, eigenvectors and least squares solutions Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Main problem of linear algebra 2: Given
More informationu n 2 4 u n 36 u n 1, n 1.
Exercise 1 Let (u n ) be the sequence defined by Set v n = u n 1 x+ u n and f (x) = 4 x. 1. Solve the equations f (x) = 1 and f (x) =. u 0 = 0, n Z +, u n+1 = u n + 4 u n.. Prove that if u n < 1, then
More informationDynamic interpretation of eigenvectors
EE263 Autumn 2015 S. Boyd and S. Lall Dynamic interpretation of eigenvectors invariant sets complex eigenvectors & invariant planes left eigenvectors modal form discrete-time stability 1 Dynamic interpretation
More informationDIFFERENTIAL EQUATIONS REVIEW. Here are notes to special make-up discussion 35 on November 21, in case you couldn t make it.
DIFFERENTIAL EQUATIONS REVIEW PEYAM RYAN TABRIZIAN Here are notes to special make-up discussion 35 on November 21, in case you couldn t make it. Welcome to the special Friday after-school special of That
More information1 Systems of Differential Equations
March, 20 7- Systems of Differential Equations Let U e an open suset of R n, I e an open interval in R and : I R n R n e a function from I R n to R n The equation ẋ = ft, x is called a first order ordinary
More information