The formulas for derivatives are particularly useful because they reduce ODEs to algebraic expressions. Consider the following ODE d 2 dx + p d
|
|
- Emily Mosley
- 5 years ago
- Views:
Transcription
1 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) x where p and q are constants. Transform both sides of the equations FT d 2 f (x) df (x) + p + qf (x) dx 2 dx ( (ik)2 + p(ik) + q) f (k) FT (R(x)) R(k) where ~ denotes the Fourier transform. We then have R(k) f (k) k 2 + ipk + q so that the solution (particular) is f (x) 1 e ikx f (k)dk 1 e ikx R(k) k 2 + ipk + q dk This very formal solution to the problem is called the integral representation of the solution. In general, we need complex integration techniques to evaluate these integrals. We will see how later. Solving ODEs with the Laplace Transform We now have all the necessary tools to solve ODEs using the Laplace transform. We consider an ODE with constant coefficients Now if d 2 y 2 + a dy L(y(t)) Y (s) + by r(t) we then have L d 2 y + a dy 2 + by L(r(t)) R(s) s 2 Y (s) sy() y'() + a sy (s) y() [ ] + by (s) R(s) which is now an algebraic equation. The solution is (s + a)y() + y'() R(s) Y (s) + s 2 + as + b s 2 + as + b Examples: y(t) L 1 (Y (s)) Page 1
2 (1) Consider the homogeneous ODE d 2 y + 4 dy 2 + 8y with initial conditions y() 2, y'(). We have Now or and or Y (s) 2s + 8 2(s + 2) s 2 + 4s + 8 (s + 2) (s + 2) s L(e at f (t)) F(s a) and L(cosωt) s 2 + ω 2 L(e 2t cos2t) L(cos2(t + 2)) (s + 2) (s + 2) ω L(e at f (t)) F(s a) and L(sinωt) s 2 + ω 2 L(e 2t sin2t) L(sin2(t + 2)) 2 (s + 2) Therefore, 2s + 8 2(s + 2) Y (s) s 2 + 4s + 8 (s + 2) (s + 2) L(e 2t cos2t) + 2L(e 2t sin2t) L(y(t)) L(2e 2t (cos2t + sin2t)) y(t) 2e 2t (cos2t + sin2t) (2) Consider the nonhomogeneous ODE 2 d 2 q + 5q 1sinωt 2 with initial conditions q() i() dq(). This ODE arises from a series circuit with an AC voltage source, a capacitor and an inductance. We have Q(s) (s + a)q() + i() s L(5sinωt) s L(5sinωt) s ω (s 2 + ω 2 )(s ) or using the partial fraction rules when a quadratic remains 5ω Q(s) (s 2 + ω 2 )(s ) A s + A 1 2 s 2 + ω + B s + B s Now for this form Page 2
3 or Therefore, P(s) B 1 (a + ib) + B 2 Q(s) / (s a) 2 + b 2 5ω A 1 (iω) + A 2 s ω B 1 (5i) + B 2 s 2 + ω 2 siω s5i sa+ib 5ω ω A, A 5ω ω 2 5ω ω 2 25 B 1, B 2 5ω ω ω Q(s) 25 ω 2 s 2 + ω ω 2 s or 5 q(t) [ sinωt sin5t ] 25 ω 2 Clearly, this form of the solution is valid only if ω 5 since the amplitude becomes unbounded. This corresponds to resonance in the circuit and we have an unbounded amplitude because there is no damping in this circuit (no resistors). If ω 5, we go back to Q(s) 5ω Q(s) (s ) A s + A (s ) + B s + B s (s ) 2 which using the same procedure gives q(t) 25[ sin5t 5t cos5t] unbounded as t gets large Green Function and Convolution Consider the special equation shown below: d 2 dx + p d 2 dx + q G(x x') δ(x x') x In this equation the inhomogeneous function is a delta-function, x' is an arbitrary point, and we denote this special solution as G(x). It is called the Green function. We now show that the solution of the equation d 2 dx + p d 2 dx + q f (x) R(x) x is easily expressed in terms of the Green function. Page 3
4 We now show that the solution of the inhomogeneous ODE can be written as the expression f (x) G(x x')r(x')dx' Substituting into the ODE we get d 2 dx + p d 2 dx + q G(x x')r(x')dx' as it should! d 2 dx + p d 2 dx + q G(x x')r(x')dx' δ(x x')r(x')dx' R(x) So the Green function representation of the solution is valid. This means if we solve for G(x-x') for a particular differential operator, then we have solved for all inhomogeneous solutions of the ODE given by that particular differential operator. As we saw earlier, an integral of the form above involving G(xx') is called a convolution, and is denoted by the symbol (*). f (x) G(x x')r(x')dx' (G * R) x The Fourier transform of a convolution is always the product of the transforms(same as for Laplace transforms) FT ((G * R) x ) 1 1 dxe ikx G(x x')r(x')dx' 1 1 dxe ik(x x ') G(x x')e ikx ' R(x')dx' 1 1 e ikx ' R(x')dx' dxe ik(x x ') G(x x') FT (G)F(R) In Linear Response Theory, the inhomogeneity function R(x) is called the input to, the solution f(x) is called the output from, the system, while the Green function is called the response function, since it describes how the system responds to the input. Examples: Consider Newton's 2nd law, Page 4
5 m d 2 x(t) F(t) force function 2 This is a 2nd-order inhomogeneous ODE with p q and R F. The Green function satisfies this equation with the force function replaced by a delta-function. A delta-function force is called an impulsive force. Thus, the Green function also describes the response of a mechanical system to an impulsive driving force. Now let us consider a damped, driven oscillator x(t) + 2β x(t) + ω 2 x(t) R(t) The solution will be of the form x(t) G(t t ')R(t ') ' where d 2 G(t t ') dg(t t ') + 2β + ω 2 2 G(t t ') δ(t t ') Let the Fourier transform of G(t-t') be given by FT (G(t t ')) G(ω) 1 e iω (t t ') G(t t ') d(t t ') Now take the Fourier transform of the Green function ODE FT d 2 G(t t ') dg(t t ') + 2β + ω 2 2 G(t t ') FT (δ(t t ')) or ω iβω + ω 1 G(ω) G(ω) 1 1 (ω 2 ω 2 ) + 2iβω Finally we compute the Green function itself by taking the inverse Fourier transform G(t) 1 e iωt 1 1 G(ω)dω e iωt (ω 2 ω 2 ) + 2iβω dω Complex integration(later) gives the result where G(t) 1 ω 1 e βt sin(ω 1 t)h (t) Page 5
6 ω 1 ω 2 β 2 and H (t) 1 t> step function t< Finally, the solution of the original equation is x(t) 1 e β (t t ') sin(ω 1 (t t '))H (t t ')R(t ') ' ω 1 t 1 e β (t t ') sin(ω 1 (t t '))R(t ') ' ω 1 Note the interesting feature that the integral over t' takes into account the effects of all driving forces occurring in the past (t' < t) due to the presence of the step function. It contains no effect due to the driving force in the future (t' > t), because these forces have not yet occurred. Hence, the result is explicitly consistent with the physical requirement of causality. Is the answer correct? Let us choose a driving force where we know the answer. In Physics 7 we discussed sinusoidal driving forces; so we choose R(t) e iωt Direct integration of the x(t) equation then gives the result e iωt x(t) (ω 2 Ω 2 ) + 2iβΩ which corresponds to the resonance amplitude formula of the damped, driven oscillator and agrees with the earlier result. More details of these integration techniques later. Examples: Consider M d 2 x + α dx + Kx d(t) 2 Assume the Fourier transform of d(t) is D(ω) and that X(ω) Substituting we get 1 e iωt x(t) x(t) 1 dωe iωt X(ω) Page 6
7 Now defining we have x(t) 1 ( Mω 2 + iαω + K)X(ω) D(ω) X(ω) dωe iωt D(ω) / M ω 2 iαω / M K / M ω ± iα 2M ± K M α 2 4M 2 D(ω) / M ω 2 iαω / M K / M 1 D(ω) / M dω (ω ω + )(ω ω ) eiωt where x(t) is the response of the oscillator to the forcing function d(t). Response to a Delta-Function Impulse Consider This gives d(t) I δ(t) D(ω) I 1 x(t) 1 I dω / M (ω ω + )(ω ω ) eiωt Using complex integration techniques which we learn later we get t < x(t) t > x(t) I KM α M sin e αt K M α 2 4M t 2 The solution looks like Page 7
8 We note that the solution obeys causality, i.e., there is no motion until after the impulse is applied. A delta-function impulse implies that x() and v x () Notice, however, that we have a solution with no unknown constants. We started with a second-order equation. Normally, that means we need to use 2 boundary conditions to completely define the solution. Somehow we have already imposed them. Where? It turns out that the choice of integration path in the complex integration was the place! Effectively, we have done the following: causality x() δ function impulse v x () I / M Fourier methods fail where the response diverges(unstable systems). In this case we must use the Laplace transform. Page 8
9 Example - Unstable Electric Circuit C o + c _ i(t) + + o -R r o - l + L o This circuit has a negative resistance (a theorist s circuit) and is therefore unstable. We let s (t) t < sinω t t > For t <, all voltages and i(t). We have for t > the integral-diffeq s (t) c (t) + r (t) + l (t) 1 di i(t) R C i(t) + L Taking the derivative wrt t we have d 2 i L R di i(t) d s C which is a second-order non-homogeneous diffeq where i(t) dependent variable and d s / is the driving term. Neither s nor d s / have a valid Fourier transform. Thus, we use Laplace transforms. We get s 2 L I(s) sr I(s) + 1 C I(s) s s (s) where I(s) L(i(t) and s (s) L( s (t)) The initial conditions (due to causality) are Page 9
10 Now i() and s (s) di t ω (s + iω )(s iω ) which gives I(s) ω s (s + iω )(s iω )(L s 2 R s + 1 ω s ) (s + iω )(s iω )L (s s + )(s s ) C where s ± R 2L ± i 1 R 2 L C Note that R R negative resistance Real(s ± ) >. Inverting(using complex integration methods) we have i(t) 1 ds i ω se st (s + iω )(s iω )L (s s + )(s s ) L For simplicity we choose ω 1, L 1 so that i(t) 1 ds i se st (s + i)(s i)(s 1+ 2i)(s 1 2i) L For t < we get i(t). For t > we close to the left, we get i(t) cos ( t +.15π ) et cos( 2t +. ) The first term is a result of the driving voltage at ω 1. It persists at a constant amplitude for all t >. The second term is the characteristic response of the circuit (determined solely by R, L,C ). It grows exponentially in time because of R < (an unstable circuit). It looks like(where we have enhanced the initial interference terms by replacing e t with e.1t. 4L 2 Page 1
11 Page 11
( ) 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 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 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 informationMathQuest: Differential Equations
MathQuest: Differential Equations Laplace Tranforms 1. True or False The Laplace transform method is the only way to solve some types of differential equations. (a) True, and I am very confident (b) True,
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 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 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 informationThe RC Circuit: An Approach with Fourier Transforms
The RC Circuit: An Approach with Fourier Transforms In this article we shall mathematically analyse the Resistor- Capacitor RC) circuit with the help of Fourier transforms FT). This very general technique
More informationDamped Harmonic Oscillator
Damped Harmonic Oscillator Wednesday, 23 October 213 A simple harmonic oscillator subject to linear damping may oscillate with exponential decay, or it may decay biexponentially without oscillating, or
More informatione st f (t) dt = e st tf(t) dt = L {t f(t)} s
Additional operational properties How to find the Laplace transform of a function f (t) that is multiplied by a monomial t n, the transform of a special type of integral, and the transform of a periodic
More informationMath 3313: Differential Equations Second-order ordinary differential equations
Math 3313: Differential Equations Second-order ordinary differential equations Thomas W. Carr Department of Mathematics Southern Methodist University Dallas, TX Outline Mass-spring & Newton s 2nd law Properties
More informationThe (Fast) Fourier Transform
The (Fast) Fourier Transform The Fourier transform (FT) is the analog, for non-periodic functions, of the Fourier series for periodic functions can be considered as a Fourier series in the limit that 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 informationModule 25: Outline Resonance & Resonance Driven & LRC Circuits Circuits 2
Module 25: Driven RLC Circuits 1 Module 25: Outline Resonance & Driven LRC Circuits 2 Driven Oscillations: Resonance 3 Mass on a Spring: Simple Harmonic Motion A Second Look 4 Mass on a Spring (1) (2)
More informationHOMEWORK 4: MATH 265: SOLUTIONS. y p = cos(ω 0t) 9 ω 2 0
HOMEWORK 4: MATH 265: SOLUTIONS. Find the solution to the initial value problems y + 9y = cos(ωt) with y(0) = 0, y (0) = 0 (account for all ω > 0). Draw a plot of the solution when ω = and when ω = 3.
More informationGATE EE Topic wise Questions SIGNALS & SYSTEMS
www.gatehelp.com GATE EE Topic wise Questions YEAR 010 ONE MARK Question. 1 For the system /( s + 1), the approximate time taken for a step response to reach 98% of the final value is (A) 1 s (B) s (C)
More information(f g)(t) = Example 4.5.1: Find f g the convolution of the functions f(t) = e t and g(t) = sin(t). Solution: The definition of convolution is,
.5. Convolutions and Solutions Solutions of initial value problems for linear nonhomogeneous differential equations can be decomposed in a nice way. The part of the solution coming from the initial data
More informationFourier series. Complex Fourier series. Positive and negative frequencies. Fourier series demonstration. Notes. Notes. Notes.
Fourier series Fourier series of a periodic function f (t) with period T and corresponding angular frequency ω /T : f (t) a 0 + (a n cos(nωt) + b n sin(nωt)), n1 Fourier series is a linear sum of cosine
More informationHow many initial conditions are required to fully determine the general solution to a 2nd order linear differential equation?
How many initial conditions are required to fully determine the general solution to a 2nd order linear differential equation? (A) 0 (B) 1 (C) 2 (D) more than 2 (E) it depends or don t know How many of
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 informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Arizona State University Lecture 5: Calculating the Laplace Transform of a Signal Introduction In this Lecture, you will learn: Laplace Transform of Simple
More informationENGIN 211, Engineering Math. Laplace Transforms
ENGIN 211, Engineering Math Laplace Transforms 1 Why Laplace Transform? Laplace transform converts a function in the time domain to its frequency domain. It is a powerful, systematic method in solving
More informationNov : Lecture 22: Differential Operators, Harmonic Oscillators
14 Nov. 3 005: Lecture : Differential Operators, Harmonic Oscillators Reading: Kreyszig Sections:.4 (pp:81 83),.5 (pp:83 89),.8 (pp:101 03) Differential Operators The idea of a function as something that
More informationMath Laplace transform
1 Math 371 - Laplace transform Erik Kjær Pedersen November 22, 2005 The functions we have treated with the power series method are called analytical functions y = a k t k The differential equations we
More information20.6. Transfer Functions. Introduction. Prerequisites. Learning Outcomes
Transfer Functions 2.6 Introduction In this Section we introduce the concept of a transfer function and then use this to obtain a Laplace transform model of a linear engineering system. (A linear engineering
More informationTHE FOURIER TRANSFORM (Fourier series for a function whose period is very, very long) Reading: Main 11.3
THE FOURIER TRANSFORM (Fourier series for a function whose period is very, very long) Reading: Main 11.3 Any periodic function f(t) can be written as a Fourier Series a 0 2 + a n cos( nωt) + b n sin n
More informationENGI 9420 Lecture Notes 1 - ODEs Page 1.01
ENGI 940 Lecture Notes - ODEs Page.0. Ordinary Differential Equations An equation involving a function of one independent variable and the derivative(s) of that function is an ordinary differential equation
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 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 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 informationAdvanced Analog Building Blocks. Prof. Dr. Peter Fischer, Dr. Wei Shen, Dr. Albert Comerma, Dr. Johannes Schemmel, etc
Advanced Analog Building Blocks Prof. Dr. Peter Fischer, Dr. Wei Shen, Dr. Albert Comerma, Dr. Johannes Schemmel, etc 1 Topics 1. S domain and Laplace Transform Zeros and Poles 2. Basic and Advanced current
More information3.5 Undetermined Coefficients
3.5. UNDETERMINED COEFFICIENTS 153 11. t 2 y + ty + 4y = 0, y(1) = 3, y (1) = 4 12. t 2 y 4ty + 6y = 0, y(0) = 1, y (0) = 1 3.5 Undetermined Coefficients In this section and the next we consider the nonhomogeneous
More informationFourier transforms. c n e inπx. f (x) = Write same thing in an equivalent form, using n = 1, f (x) = l π
Fourier transforms We can imagine our periodic function having periodicity taken to the limits ± In this case, the function f (x) is not necessarily periodic, but we can still use Fourier transforms (related
More informationSolving a RLC Circuit using Convolution with DERIVE for Windows
Solving a RLC Circuit using Convolution with DERIVE for Windows Michel Beaudin École de technologie supérieure, rue Notre-Dame Ouest Montréal (Québec) Canada, H3C K3 mbeaudin@seg.etsmtl.ca - Introduction
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 informationInductance, Inductors, RL Circuits & RC Circuits, LC, and RLC Circuits
Inductance, Inductors, RL Circuits & RC Circuits, LC, and RLC Circuits Self-inductance A time-varying current in a circuit produces an induced emf opposing the emf that initially set up the timevarying
More information4.2 Homogeneous Linear Equations
4.2 Homogeneous Linear Equations Homogeneous Linear Equations with Constant Coefficients Consider the first-order linear differential equation with constant coefficients a 0 and b. If f(t) = 0 then this
More informationENGI Second Order Linear ODEs Page Second Order Linear Ordinary Differential Equations
ENGI 344 - Second Order Linear ODEs age -01. Second Order Linear Ordinary Differential Equations The general second order linear ordinary differential equation is of the form d y dy x Q x y Rx dx dx Of
More informationTime Response of Systems
Chapter 0 Time Response of Systems 0. Some Standard Time Responses Let us try to get some impulse time responses just by inspection: Poles F (s) f(t) s-plane Time response p =0 s p =0,p 2 =0 s 2 t p =
More informationAPPLICATIONS OF SECOND-ORDER DIFFERENTIAL EQUATIONS
APPLICATIONS OF SECOND-ORDER DIFFERENTIAL EQUATIONS Second-order linear differential equations have a variety of applications in science and engineering. In this section we explore two of them: the vibration
More informationLecture 6: Differential Equations Describing Vibrations
Lecture 6: Differential Equations Describing Vibrations In Chapter 3 of the Benson textbook, we will look at how various types of musical instruments produce sound, focusing on issues like how the construction
More information7. Find the Fourier transform of f (t)=2 cos(2π t)[u (t) u(t 1)]. 8. (a) Show that a periodic signal with exponential Fourier series f (t)= δ (ω nω 0
Fourier Transform Problems 1. Find the Fourier transform of the following signals: a) f 1 (t )=e 3 t sin(10 t)u (t) b) f 1 (t )=e 4 t cos(10 t)u (t) 2. Find the Fourier transform of the following signals:
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 informationMatemáticas II: Segundo del Grado en Ingeniería Aeroespacial
Matemáticas II: Segundo del Grado en Ingeniería Aeroespacial Sergio Blanes, Dolors Roselló Ecuaciones diferenciales y transformadas de Laplace con aplicaciones L.M. Sánchez, M.P. Legua Ref.: 211-798 Capítulo
More informationWork sheet / Things to know. Chapter 3
MATH 251 Work sheet / Things to know 1. Second order linear differential equation Standard form: Chapter 3 What makes it homogeneous? We will, for the most part, work with equations with constant coefficients
More informationLINEAR SYSTEMS. J. Elder PSYC 6256 Principles of Neural Coding
LINEAR SYSTEMS Linear Systems 2 Neural coding and cognitive neuroscience in general concerns input-output relationships. Inputs Light intensity Pre-synaptic action potentials Number of items in display
More informationLaplace Transforms Chapter 3
Laplace Transforms Important analytical method for solving linear ordinary differential equations. - Application to nonlinear ODEs? Must linearize first. Laplace transforms play a key role in important
More informationPractice Problems For Test 3
Practice Problems For Test 3 Power Series Preliminary Material. Find the interval of convergence of the following. Be sure to determine the convergence at the endpoints. (a) ( ) k (x ) k (x 3) k= k (b)
More informationL = 1 2 a(q) q2 V (q).
Physics 3550 Motion near equilibrium - Small Oscillations Relevant Sections in Text: 5.1 5.6, 11.1 11.3 Motion near equilibrium 1 degree of freedom One of the most important situations in physics is motion
More informationChapter 1 Fundamental Concepts
Chapter 1 Fundamental Concepts 1 Signals A signal is a pattern of variation of a physical quantity, often as a function of time (but also space, distance, position, etc). These quantities are usually the
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 informationRaktim Bhattacharya. . AERO 422: Active Controls for Aerospace Vehicles. Dynamic Response
.. AERO 422: Active Controls for Aerospace Vehicles Dynamic Response Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. . Previous Class...........
More informationApplications of Second-Order Differential Equations
Applications of Second-Order Differential Equations ymy/013 Building Intuition Even though there are an infinite number of differential equations, they all share common characteristics that allow intuition
More informationECEN 420 LINEAR CONTROL SYSTEMS. Lecture 2 Laplace Transform I 1/52
1/52 ECEN 420 LINEAR CONTROL SYSTEMS Lecture 2 Laplace Transform I Linear Time Invariant Systems A general LTI system may be described by the linear constant coefficient differential equation: a n d n
More informationf(t)e st dt. (4.1) Note that the integral defining the Laplace transform converges for s s 0 provided f(t) Ke s 0t for some constant K.
4 Laplace transforms 4. Definition and basic properties The Laplace transform is a useful tool for solving differential equations, in particular initial value problems. It also provides an example of integral
More informationDifferential Equations and Linear Algebra Lecture Notes. Simon J.A. Malham. Department of Mathematics, Heriot-Watt University, Edinburgh EH14 4AS
Differential Equations and Linear Algebra Lecture Notes Simon J.A. Malham Department of Mathematics, Heriot-Watt University, Edinburgh EH4 4AS Contents Chapter. Linear second order ODEs 5.. Newton s second
More informationPractice Problems For Test 3
Practice Problems For Test 3 Power Series Preliminary Material. Find the interval of convergence of the following. Be sure to determine the convergence at the endpoints. (a) ( ) k (x ) k (x 3) k= k (b)
More informationNotes on Fourier Series and Integrals Fourier Series
Notes on Fourier Series and Integrals Fourier Series et f(x) be a piecewise linear function on [, ] (This means that f(x) may possess a finite number of finite discontinuities on the interval). Then f(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 information2.161 Signal Processing: Continuous and Discrete Fall 2008
MIT OpenCourseWare http://ocw.mit.edu 2.6 Signal Processing: Continuous and Discrete Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. MASSACHUSETTS
More informationLinear Second Order ODEs
Chapter 3 Linear Second Order ODEs In this chapter we study ODEs of the form (3.1) y + p(t)y + q(t)y = f(t), where p, q, and f are given functions. Since there are two derivatives, we might expect that
More informationSection 6.4 DEs with Discontinuous Forcing Functions
Section 6.4 DEs with Discontinuous Forcing Functions Key terms/ideas: Discontinuous forcing function in nd order linear IVPs Application of Laplace transforms Comparison to viewing the problem s solution
More informationLaplace Transform Problems
AP Calculus BC Name: Laplace Transformation Day 3 2 January 206 Laplace Transform Problems Example problems using the Laplace Transform.. Solve the differential equation y! y = e t, with the initial value
More information+ + LAPLACE TRANSFORM. Differentiation & Integration of Transforms; Convolution; Partial Fraction Formulas; Systems of DEs; Periodic Functions.
COLOR LAYER red LAPLACE TRANSFORM Differentiation & Integration of Transforms; Convolution; Partial Fraction Formulas; Systems of DEs; Periodic Functions. + Differentiation of Transforms. F (s) e st f(t)
More informationComputer Problems for Methods of Solving Ordinary Differential Equations
Computer Problems for Methods of Solving Ordinary Differential Equations 1. Have a computer make a phase portrait for the system dx/dt = x + y, dy/dt = 2y. Clearly indicate critical points and separatrices.
More informationLaplace Transforms and use in Automatic Control
Laplace Transforms and use in Automatic Control P.S. Gandhi Mechanical Engineering IIT Bombay Acknowledgements: P.Santosh Krishna, SYSCON Recap Fourier series Fourier transform: aperiodic Convolution integral
More information1 Notes on Laplace circuit analysis
Physics - David Kleinfeld - Spring 7 Notes on aplace circuit analysis Background We previously learned that we can transform from the time domain to the frequency domain under steady-state conditions thus
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 information2.4 Models of Oscillation
2.4 Models of Oscillation In this section we give three examples of oscillating physical systems that can be modeled by the harmonic oscillator equation. Such models are ubiquitous in physics, but are
More informationThe Laplace Transform and the IVP (Sect. 6.2).
The Laplace Transform and the IVP (Sect..2). Solving differential equations using L ]. Homogeneous IVP. First, second, higher order equations. Non-homogeneous IVP. Recall: Partial fraction decompositions.
More informationSeries FOURIER SERIES. Graham S McDonald. A self-contained Tutorial Module for learning the technique of Fourier series analysis
Series FOURIER SERIES Graham S McDonald A self-contained Tutorial Module for learning the technique of Fourier series analysis Table of contents Begin Tutorial c 24 g.s.mcdonald@salford.ac.uk 1. Theory
More informationElectromagnetic Oscillations and Alternating Current. 1. Electromagnetic oscillations and LC circuit 2. Alternating Current 3.
Electromagnetic Oscillations and Alternating Current 1. Electromagnetic oscillations and LC circuit 2. Alternating Current 3. RLC circuit in AC 1 RL and RC circuits RL RC Charging Discharging I = emf R
More informationMath 221 Topics since the second exam
Laplace Transforms. Math 1 Topics since the second exam There is a whole different set of techniques for solving n-th order linear equations, which are based on the Laplace transform of a function. For
More informationIntroduction ODEs and Linear Systems
BENG 221 Mathematical Methods in Bioengineering ODEs and Linear Systems Gert Cauwenberghs Department of Bioengineering UC San Diego 1.1 Course Objectives 1. Acquire methods for quantitative analysis and
More informationLinear Systems Theory
ME 3253 Linear Systems Theory Review Class Overview and Introduction 1. How to build dynamic system model for physical system? 2. How to analyze the dynamic system? -- Time domain -- Frequency domain (Laplace
More informationExam. Matrikelnummer: Points. Question Bonus. Total. Grade. Information Theory and Signal Reconstruction Summer term 2013
Exam Name: Matrikelnummer: Question 1 2 3 4 5 Bonus Points Total Grade 1/6 Question 1 You are traveling to the beautiful country of Markovia. Your travel guide tells you that the weather w i in Markovia
More informationReview of Linear Time-Invariant Network Analysis
D1 APPENDIX D Review of Linear Time-Invariant Network Analysis Consider a network with input x(t) and output y(t) as shown in Figure D-1. If an input x 1 (t) produces an output y 1 (t), and an input x
More informationPHYSICS 110A : CLASSICAL MECHANICS HW 2 SOLUTIONS. Here is a sketch of the potential with A = 1, R = 1, and S = 1. From the plot we can see
PHYSICS 11A : CLASSICAL MECHANICS HW SOLUTIONS (1) Taylor 5. Here is a sketch of the potential with A = 1, R = 1, and S = 1. From the plot we can see 1.5 1 U(r).5.5 1 4 6 8 1 r Figure 1: Plot for problem
More informationMath 216 Second Midterm 19 March, 2018
Math 26 Second Midterm 9 March, 28 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 that
More informationInductance, RL and RLC Circuits
Inductance, RL and RLC Circuits Inductance Temporarily storage of energy by the magnetic field When the switch is closed, the current does not immediately reach its maximum value. Faraday s law of electromagnetic
More informationSecond order linear equations
Second order linear equations Samy Tindel Purdue University Differential equations - MA 266 Taken from Elementary differential equations by Boyce and DiPrima Samy T. Second order equations Differential
More informationName:... Section:... Physics 208 Quiz 8. April 11, 2008; due April 18, 2008
Name:... Section:... Problem 1 (6 Points) Physics 8 Quiz 8 April 11, 8; due April 18, 8 Consider the AC circuit consisting of an AC voltage in series with a coil of self-inductance,, and a capacitor of
More informationMath 215/255 Final Exam (Dec 2005)
Exam (Dec 2005) Last Student #: First name: Signature: Circle your section #: Burggraf=0, Peterson=02, Khadra=03, Burghelea=04, Li=05 I have read and understood the instructions below: Please sign: Instructions:.
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 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 informationHandout 11: AC circuit. AC generator
Handout : AC circuit AC generator Figure compares the voltage across the directcurrent (DC) generator and that across the alternatingcurrent (AC) generator For DC generator, the voltage is constant For
More informationLecture 7. Please note. Additional tutorial. Please note that there is no lecture on Tuesday, 15 November 2011.
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
More informationMath Fall Linear Filters
Math 658-6 Fall 212 Linear Filters 1. Convolutions and filters. A filter is a black box that takes an input signal, processes it, and then returns an output signal that in some way modifies the input.
More information2.4 Harmonic Oscillator Models
2.4 Harmonic Oscillator Models In this section we give three important examples from physics of harmonic oscillator models. Such models are ubiquitous in physics, but are also used in chemistry, biology,
More informationSome of the different forms of a signal, obtained by transformations, are shown in the figure. jwt e z. jwt z e
Transform methods Some of the different forms of a signal, obtained by transformations, are shown in the figure. X(s) X(t) L - L F - F jw s s jw X(jw) X*(t) F - F X*(jw) jwt e z jwt z e X(nT) Z - Z X(z)
More informationIntroduction to the Phase Plane
Introduction to the Phase Plane June, 06 The Phase Line A single first order differential equation of the form = f(y) () makes no mention of t in the function f. Such a differential equation is called
More informationOscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. Oscillatory solution
Oscillatory Motion Simple pendulum: linear Hooke s Law restoring force for small angular deviations d 2 θ dt 2 = g l θ θ l Oscillatory solution θ(t) =θ 0 sin(ωt + φ) F with characteristic angular frequency
More informationANSWERS Final Exam Math 250b, Section 2 (Professor J. M. Cushing), 15 May 2008 PART 1
ANSWERS Final Exam Math 50b, Section (Professor J. M. Cushing), 5 May 008 PART. (0 points) A bacterial population x grows exponentially according to the equation x 0 = rx, where r>0is the per unit rate
More informationOrdinary Differential Equation Theory
Part I Ordinary Differential Equation Theory 1 Introductory Theory An n th order ODE for y = y(t) has the form Usually it can be written F (t, y, y,.., y (n) ) = y (n) = f(t, y, y,.., y (n 1) ) (Implicit
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 informationOscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. small angle approximation. Oscillatory solution
Oscillatory Motion Simple pendulum: linear Hooke s Law restoring force for small angular deviations d 2 θ dt 2 = g l θ small angle approximation θ l Oscillatory solution θ(t) =θ 0 sin(ωt + φ) F with characteristic
More informationSpectral Broadening Mechanisms
Spectral Broadening Mechanisms Lorentzian broadening (Homogeneous) Gaussian broadening (Inhomogeneous, Inertial) Doppler broadening (special case for gas phase) The Fourier Transform NC State University
More informationHonors Differential Equations
MIT OpenCourseWare http://ocw.mit.edu 18.034 Honors Differential Equations Spring 009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. LECTURE 7. MECHANICAL
More informationEE 3054: Signals, Systems, and Transforms Summer It is observed of some continuous-time LTI system that the input signal.
EE 34: Signals, Systems, and Transforms Summer 7 Test No notes, closed book. Show your work. Simplify your answers. 3. It is observed of some continuous-time LTI system that the input signal = 3 u(t) produces
More informationLinear Filters. L[e iωt ] = 2π ĥ(ω)eiωt. Proof: Let L[e iωt ] = ẽ ω (t). Because L is time-invariant, we have that. L[e iω(t a) ] = ẽ ω (t a).
Linear Filters 1. Convolutions and filters. A filter is a black box that takes an input signal, processes it, and then returns an output signal that in some way modifies the input. For example, if the
More information