Using Mathematica to Teach Linear Differential Operators and the Method of Undetermined Coefficients
|
|
- Ashley Casey
- 5 years ago
- Views:
Transcription
1 Using Mathematica to Teach Linear Differential Operators and the Method of Undetermined oefficients Itai Seggev Knox ollege / Wolfram Research Joint Mathematics Meetings January 9, 2011 New Orleans, Louisiana
2 2 JMM2011.nb Problem Statement Given a linear differential operator PD and a polynomial-exponential forcing function f, we wish to find the general solution to the homogeneous equation PDy 0 one solution to the driven equation PDy f by means of an ansatz We can then solve any IVP for PD, but we need to enter PD into Mathematica
3 JMM2011.nb Ingredients In[11]:= Solve[] oefficientlist[] omplexexpand[], Re[], Im[] (if dealing with complex and or trigonometric equations) Some magic code: leardpluss DPlusS : PowerDPlusSt_, s_:0, n_ : Functionf,NestD, t s &, f, n DPlusSt_, s_:0^n_f ^: DPlusSt, s n f DPlusSt_,s_:0f_: Df,ts f UnprotectPlus, Times; Plus : Plus a_:1 o1_function, cruft f_ : a o1f Pluscruftf Plus : Plus a_:1dplusst_, s, cruft f_ : a DPlusSt, sf Pluscruftf Times : Timesa_, DPlusSt_,s f_ : a DPlusSt,sf Times : Timesa_, o1_function f_ : a o1f ProtectPlus, Times;
4 4 JMM2011.nb Playing with DPlusSt, s D s The basic operator In[21]:= Out[21]= DPlusSt, syt s yt y t Monomials Positive powers: In[22]:= DPlusSt, 0^2yt DPlusSt^2yt Out[22]= Out[23]= y t y t Positive powers with shift: In[24]:= Out[24]= DPlusSt, s^2yt Expand s 2 yt 2 s y t y t The identity operator In[25]:= Out[25]= DPlusSt, s^0yt yt A polynomial In[26]:= Out[28]= In[29]:= learpofd PofDt_, s_ : 0 : L DPlusSt, s^2 R DPlusSt, s PofDt, 0yt yt R y t L y t PofDt, Ωyt Simplify DPlusSt, s0 Out[29]= 1 Ω R L Ω yt R 2 L Ω y t L y t
5 JMM2011.nb Step 1: Finding the basic solution set Simple roots Assuming a constant coefficient operator, we know the solution is given by the characteristic roots In[30]:= PofDt_, s_ : 0 : DPlusSt, s 2 5 DPlusSt, s 6 DPlusSt, s 0 The characteristic polynomial is obtain with D 0, i.e., computing PD1 In[31]:= PofDt, s1 Solve 0, s Exps t. Out[31]= 6 5 s s 2 Out[32]= Out[33]= s 3, s 2 3 t, 2 t So, we have our basic solution set. Repeated roots In[34]:= PofDt_, s_ : 0 : DPlusSt, s 2 4 DPlusSt, s 4 DPlusSt, s 0 PofDt, s1 Solve 0, s Exps t. Out[35]= 4 4 s s 2 Out[36]= Out[37]= s 2, s 2 2 t, 2 t Our set is linearly dependent, we need to replace one vector with t e 2 t
6 6 JMM2011.nb Step 2: Find the particular solution (real case) In[38]:= learpofd PofDt_, s_ : 0 : DPlusSt, s 2 5 DPlusSt, s 6 DPlusSt, s 0 A couple of forcing functions In[40]:= f1t_ : t 2 3 f2t_ : Exp7 t t 2 f 1 We know that the particular solution has the same form as the driving polynomial in the case of a simple polynomial driving term. In[42]:= y1t_ : a2 t 2 a1 t a0 We must substitute in and solve for the coefficients. oefficientlist[] to the rescue! In[43]:= Out[43]= oefficientlistf1t, t oefficientlistpofdty1t, t Solve, a0, a1, a2 y1t. First 3, 0, 1 Out[44]= 6 a0 5 a1 2 a2, 6 a1 10 a2, 6 a2 Out[45]= a , a1 5 18, a2 1 6 Out[46]= t t f 2 When there is an exponential term, it is conventional to separate the solution into its polynomial and exponential parts. In[47]:= leary2, h2 h2t_ : a1 t a0 y2t_ : h2t Exp7 t We then equate the polynomial part of the driving terms with the shifted operator acting on the polynomial ansatz:
7 JMM2011.nb In[50]:= Out[50]= oefficientlistf2t Exp7 t, t oefficientlistpofdt, 7h2t, t Solve, a0, a1 y2t. First 2, 1 Out[51]= 90 a0 19 a1, 90 a1 Out[52]= a0 161 Out[53]= 7 t t 8100, a Alternate solution While it is conventional to do the separation, it is not strictly required. oefficientlist[] will consider the exponential part of the "coefficient", allowing its cancellation. In[54]:= oefficientlistf2t, t oefficientlistpofdty2t, t Solve, a0, a1 y2t. First Out[54]= Out[55]= 2 7 t, 7 t 90 a0 7 t 19 a1 7 t, 90 a1 7 t Out[56]= a0 161 Out[57]= 7 t t 8100, a
8 8 JMM2011.nb Step 2: Find the particular solution (trig/complex cases) Sinusoidal driving term In[58]:= learpofd, f PofDt_, s_ : 0 : DPlusSt, s 2 5 DPlusSt, s 6 DPlusSt, s 0 ft_ : Sin7 t t 2 Since it we have a sinusoid, we need to complexify the problem, solving with a driving term f c t t 2 7 t. In[61]:= learyc, hc hct_ : a1 t a0 fct_ : Exp7 t t 2 yct_ : hct Exp7 t Finding the complex solution is as before. In[65]:= Out[65]= oefficientlistfct Exp7 t, t oefficientlistpofdt, 7 hct, t Solve, a0, a1 yct_ yct. First 2, 1 Out[66]= a a1, a Out[67]= a Out[68]= 7 t , a t We can then find the solution by taking the imaginary part of our complex solution. In[69]:= Out[69]= yt_ omplexexpandimyct Simplify t os7 t t Sin7 t Notice that the omplexexpand[] is essential, since the value above is not valid if t is complex. In[70]:= Out[70]= Imyct Simplify Im t t For free, we also get the solution for a cosine driving term. In[71]:= Out[71]= omplexexpandreyct Simplify t os7 t t Sin7 t Alternate solution As before, we are not obligated to separate out the polynomial if we don't want to.
9 JMM2011.nb In[72]:= yct_ : hct Exp7 t oefficientlistfct, t oefficientlistpofdtyct, t Solve, a0, a1 Out[73]= Out[74]= 2 7 t, 7 t a0 7 t 5 14 a1 7 t, a1 7 t Out[75]= a , a
10 10 JMM2011.nb Bonus slide: additional cases Degeneracy A degeneracy occurs when the characteristic polynomial has a repeated root, and/or the shift equals one of the roots. In[76]:= In[79]:= learpofd, f PofDt_, s_ : 0 : DPlusSt, s 2 4 DPlusSt, s 4 DPlusSt, s 0 ft_ : t^2 2 t 3 Exp2 t leary, h ht_ : a2 t^2 a1 t a0 yt_ : ht Exp2 t Now, what happens if we proceed with our normal ansatz without paying attention to the degeneracy? In[82]:= Out[82]= oefficientlistft Exp2 t, t oefficientlistpofdt, 2ht, t Solve, a0, a1, a2 3, 2, 1 Out[83]= Out[84]= 2 a2 We get no solution: something has gone wrong. But, there is more information. The two lowest terms have dropped out entirely: we have a double degeneracy (repeated root equal to shift) Hence, we use a polynomial of degree two higher (but possessing the same number of terms). In[85]:= Out[86]= ht_ : a4 t^4 a3 t^3 a2 t^2 oefficientlistft Exp2 t, t oefficientlistpofdt, 2ht, t Solve, a2, a3, a4 yt. First 3, 2, 1 Out[87]= 2 a2, 6 a3, 12 a4 Out[88]= a2 3 2, a3 1 1, a Out[89]= 2 t 3 t 2 2 t3 3 t4 12 Higher order operators In[90]:= learpofd, f PofDt_, s_ : 0 : DPlusSt, s 4 4 DPlusSt, s 3 3 DPlusSt, s 2 4 DPlusSt, s 4 DPlusSt, s 0 ft_ : t^2 2 t 3 Expt haracteristic roots are -2, -2, 1, and -1.
11 JMM2011.nb In[93]:= PofDt, s1 Simplify Out[93]= 2 s 2 1 s 2 We have single degenerac: shift equal to a simple root. So we need a polynomial of one degree higher. Rest of process is unchanged. In[94]:= Out[97]= leary, h ht_ : a3 t^3 a2 t^2 a1 t yt_ : ht Expt oefficientlistft Expt, t oefficientlistpofdt, 1ht, t Solve, a1, a2, a3 yt. First 3, 2, 1 Out[98]= 2 a1 6 a2, 4 a2 18 a3, 6 a3 Out[99]= a1 9 4, a2 1 4, a3 1 6 Out[100]= t 9 t 4 t2 t3 4 6 RL circuit (complete example from lecture notes) An RL circuit obeys L q.. R q q t (1) If we wanted to express the relationship in terms of current, we would need to differentiate the previous equation to get: L İ. R I I t. (2) Getting equation (2) is helpful because expressing charge in terms of current involves a definite integral: qt qt 0 t0 t Is s. If we substitute this equation into equation (1), we'd get a so-called integro-differential equation. Differentiating both sides with respect to t gives us a plain old differential equation. For convenience, we put a bunch of assumptions into the global variable $Assumptions which will shorten many Simplify[] commands. In[101]:= Out[101]= $Assumptions a Reals, R 0, Ω 0, L 0, 0, t Reals a Reals, R 0, Ω 0, L 0, 0, t Reals We define our differential operator In[102]:= learpofd PofDt_, s_ : 0 : L DPlusSt, s^2 R DPlusSt, s DPlusSt, s 0 We want to solve equation (2) assuming t a cosωt. First, we must find the characteristic roots:
12 12 JMM2011.nb In[104]:= SolvePofDt, s1 0, s Simplify Out[104]= s R 4 L 2 L R2, s R 4 L 2 L R2 As long as R is non-zero (i.e., there actually is a resistor in the circuit), the characteristic roots have non-zero real part and hence cannot equal s Ω. This eliminates degenerate cases of undetermined coefficients. We use Mathematica to find the solution, but as in example 3 this is unnecessary. Since the polynomial multiplying the cosine is a constant, the polynomial ht in the driven solution is a ht P Ω We ask Mathematica to verify: In[105]:= Out[107]= learh, y, z ht_ : a0 SolveoefficientListPofDt, Ωht, t a, a0 zt_ ht Ω t. 1 a a0 1 R Ω L Ω2 Out[108]= 1 a t Ω R Ω L Ω2 Since we're focused on a cosine driving term, we need to take the real part of z. The following would be the steps corresponding to finding the real part by hand. We make the denominator real by multiplying by its complex conjugate: In[109]:= Out[109]= 1 R Ω L 1 Ω2 R Ω L Ω2 1 2 L Ω2 2 R 2 Ω 2 L 2 Ω 4 Simplify We do the same thing to the numerator and use Euler's identity to convert to sines and cosines. 1 In[110]:= ExpToTriga t Ω R Ω L Ω2 Expand Out[110]= a ost Ω a Sint Ω a R Ω ost Ω a L Ω 2 ost Ω a R Ω Sint Ω a L Ω 2 Sint Ω We then take combine the numerator and denominator. Since the denominator is now real, this simply means take the real part of the numerator and keep the denominator unchanged. In[111]:= Out[111]= yprovisional FullSimplify Re a 1 L Ω 2 ost Ω R Ω Sint Ω 1 Ω 2 R 2 L 2 L Ω 2 Notice that we could have gotten to y provisional in one step using omplexexpand[]:
13 JMM2011.nb In[112]:= yprovisional omplexexpandrezt Out[112]= a ost Ω a L Ω2 ost Ω a R Ω Sint Ω R 2 Ω 2 1 L Ω2 2 R 2 Ω 2 1 L Ω2 2 R 2 Ω 2 1 L Ω2 2 Also, here is the imaginary part. It looks similar, but the roles of sine and cosine have been reversed: In[113]:= Out[113]= omplexexpandimzt a R Ω ost Ω R 2 Ω 2 1 L Ω2 2 a Sint Ω a L Ω2 Sint Ω R 2 Ω 2 1 L Ω2 2 R 2 Ω 2 1 L Ω2 2 Let Θ arctan 1 L Ω2, R Ω. In[114]:= φ ArcTan 1 L Ω2, R Ω Out[114]= ArcTan 1 L Ω2, R Ω The two-variable arctangent gives an angle between Π and Π, taking into account the quadrant of the point. Then the expression in the Evaluate[] in the cell below is equivalent to y provisional. Notice that they have the same denominator, and the equality of the numerators follows from the angle addition formulas and then basic facts cosarctanx, y x, sinarctanx, y y. Just to be sure, we ask Mathematica to verify: In[115]:= Out[115]= a osω t φ FullSimplifyEvaluate TrigExpand yprovisional, 1 L Ω2 2 R Ω 2 TransformationFunctions Automatic,. osarctand_, n_ d, SinArcTand_, n_ n & True So, the general solution is one of the following The only remaining issue is what is the size of 4 L versus R 2, which affects the homogeneous solution. When 4 L is smaller, we have two decaying exponentials. When 4 L is larger, we have two oscillating, decaying exponentials. When they are equal, we have an exponential and a polynomial exponential. This is the electrical equivalent of a critically damped harmonic oscillator, in which the system returns to equilibrium as quickly as possible. learly it is a very interesting case from an engineering standpoint. So here are the three possible solutions: In[116]:= yoverdampedt_ : a osω t φ 1 L Ω2 2 R Ω 2 c1 R 4 L R2 2 L t c2 R 4 L R2 2 L t yriticallydampedt_ : yunderdampedt_ : a osω t φ 1 L Ω2 2 R Ω 2 c1 R t R 2 L c2 t t 2 L a osω t φ 1 L Ω2 2 R Ω 2 c1 R t 2 L os 4 L 2 L R2 t c2 R 2 L t Sin 4 L 2 L R2 t
14 14 JMM2011.nb All in three cases, however, the homogeneous solution decays to zero, leaving just the driven solution asymptotically. This piece has its largest amplitude when the denominator is minimized, which for fixed R occurs at 1 L Ω2, or Ω 1 L. Recall that 1 L is what we found for the natural frequency of an L circuit, so the amplitude is maximized at resonance. Indeed, if there is no resistance, the amplitude will grow without bound, as we found on problem Being finished with our assumptions, we reset them by assigning True to $Assumptions In[119]:= Out[119]= $Assumptions True True
15 JMM2011.nb Is that a parameter which I see varying before me? Who said we need constant coefficients? In[120]:= learpofd, y PofDt_ : DPlusSt DPlusSt t t^2 DPlusSt0 Basic solution set We can construct a basic solution set in the standard manner In[122]:= DSolvePofDtyt 0, y1 1, y'1 0, yt, t DSolvePofDtyt 0, y1 0, y'1 1, yt, t Wronskianyt. Join,, t Out[122]= Out[123]= yt 2 t t 2 yt t t 2 Out[124]= t 2 By adding and subtracting the first to/from the second, we would get t, t 2 is a basic solution set In[125]:= Out[128]= 0 Out[129]= 0 leary1, y2 y1t_ : t y2t_ : t^2 PofDty1t PofDty2t Wronskiany1t, y2t, t Out[130]= t 2 Variation of Parameters onsider a logarithmic forcing function In[131]:= ft_ : Logt 2 We can construct the
16 16 JMM2011.nb In[132]:= y2s fs ygent_ SimplifyIntegrate, s, 1, t y1t Wronskiany1s, y2s, s y1s fs Integrate, s, 1, t y2t Wronskiany1s, y2s, s c1 y1t c2 y2t, Assumptions t 0 Expand Out[132]= 2 t c1 t 2 t 2 c2 t 2 2 t 2 Logt t 2 Logt t2 Logt 3 In[133]:= PofDtyGent Simplify ygen1 ygen'1 Out[133]= Logt 2 Out[134]= Out[135]= In[136]:= c1 c2 c1 2 c2 PlotEvaluateyGent. c1 0, c2 0, t, 0, Out[136]=
17 JMM2011.nb Use and Reception ode given as black box to students advanced S students can analyze/appreciate it Students assigned these problems, including higher order equations Students had little difficulty mastering the use these operators and showed they could self correct if started with wrong polynomial Emphasize operator/linear algebra view There is a single theory of linear equations!
2.3 Oscillation. The harmonic oscillator equation is the differential equation. d 2 y dt 2 r y (r > 0). Its solutions have the form
2. Oscillation So far, we have used differential equations to describe functions that grow or decay over time. The next most common behavior for a function is to oscillate, meaning that it increases and
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 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 information1. (10 points) Find the general solution to the following second-order differential equation:
Math 307A, Winter 014 Midterm Solutions Page 1 of 8 1. (10 points) Find the general solution to the following second-order differential equation: 4y 1y + 9y = 9t. To find the general solution to this nonhomogeneous
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 information/ \ ( )-----/\/\/\/ \ / In Lecture 3 we offered this as an example of a first order LTI system.
18.03 Class 17, March 12, 2010 Linearity and time invariance [1] RLC [2] Superposition III [3] Time invariance [4] Review of solution methods [1] We've spent a lot of time with mx" + bx' + cx = q(t). There
More informationChapter 10: Sinusoidal Steady-State Analysis
Chapter 10: Sinusoidal Steady-State Analysis 1 Objectives : sinusoidal functions Impedance use phasors to determine the forced response of a circuit subjected to sinusoidal excitation Apply techniques
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 information12.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 informationSecond In-Class Exam Solutions Math 246, Professor David Levermore Thursday, 31 March 2011
Second In-Class Exam Solutions Math 246, Professor David Levermore Thursday, 31 March 211 (1) [6] Give the interval of definition for the solution of the initial-value problem d 4 y dt 4 + 7 1 t 2 dy dt
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 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 informationForced Mechanical Vibrations
Forced Mechanical Vibrations Today we use methods for solving nonhomogeneous second order linear differential equations to study the behavior of mechanical systems.. Forcing: Transient and Steady State
More informationTable of contents. d 2 y dx 2, As the equation is linear, these quantities can only be involved in the following manner:
M ath 0 1 E S 1 W inter 0 1 0 Last Updated: January, 01 0 Solving Second Order Linear ODEs Disclaimer: This lecture note tries to provide an alternative approach to the material in Sections 4. 4. 7 and
More informationLecture Notes for Math 251: ODE and PDE. Lecture 16: 3.8 Forced Vibrations Without Damping
Lecture Notes for Math 25: ODE and PDE. Lecture 6:.8 Forced Vibrations Without Damping Shawn D. Ryan Spring 202 Forced Vibrations Last Time: We studied non-forced vibrations with and without damping. We
More informationAPPM 2360: Midterm 3 July 12, 2013.
APPM 2360: Midterm 3 July 12, 2013. ON THE FRONT OF YOUR BLUEBOOK write: (1) your name, (2) your instructor s name, (3) your recitation section number and (4) a grading table. Text books, class notes,
More information37. f(t) sin 2t cos 2t 38. f(t) cos 2 t. 39. f(t) sin(4t 5) 40.
28 CHAPTER 7 THE LAPLACE TRANSFORM EXERCISES 7 In Problems 8 use Definition 7 to find {f(t)} 2 3 4 5 6 7 8 9 f (t),, f (t) 4,, f (t) t,, f (t) 2t,, f (t) sin t,, f (t), cos t, t t t 2 t 2 t t t t t t t
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 informationLinear Second-Order Differential Equations LINEAR SECOND-ORDER DIFFERENTIAL EQUATIONS
11.11 LINEAR SECOND-ORDER DIFFERENTIAL EQUATIONS A Spring with Friction: Damped Oscillations The differential equation, which we used to describe the motion of a spring, disregards friction. But there
More informationDynamic circuits: Frequency domain analysis
Electronic Circuits 1 Dynamic circuits: Contents Free oscillation and natural frequency Transfer functions Frequency response Bode plots 1 System behaviour: overview 2 System behaviour : review solution
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 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 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 informationA: Brief Review of Ordinary Differential Equations
A: Brief Review of Ordinary Differential Equations Because of Principle # 1 mentioned in the Opening Remarks section, you should review your notes from your ordinary differential equations (odes) course
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 informationDamped harmonic motion
Damped harmonic motion March 3, 016 Harmonic motion is studied in the presence of a damping force proportional to the velocity. The complex method is introduced, and the different cases of under-damping,
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 informationSection 3: Complex numbers
Essentially: Section 3: Complex numbers C (set of complex numbers) up to different notation: the same as R 2 (euclidean plane), (i) Write the real 1 instead of the first elementary unit vector e 1 = (1,
More informationMATH 251 Week 6 Not collected, however you are encouraged to approach all problems to prepare for exam
MATH 51 Week 6 Not collected, however you are encouraged to approach all problems to prepare for exam A collection of previous exams could be found at the coordinator s web: http://www.math.psu.edu/tseng/class/m51samples.html
More informationspring mass equilibrium position +v max
Lecture 20 Oscillations (Chapter 11) Review of Simple Harmonic Motion Parameters Graphical Representation of SHM Review of mass-spring pendulum periods Let s review Simple Harmonic Motion. Recall we used
More informationPHYSICS. Chapter 15 Lecture FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E RANDALL D. KNIGHT Pearson Education, Inc.
PHYSICS FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E Chapter 15 Lecture RANDALL D. KNIGHT Chapter 15 Oscillations IN THIS CHAPTER, you will learn about systems that oscillate in simple harmonic
More informationCHAPTER 4 FOURIER SERIES S A B A R I N A I S M A I L
CHAPTER 4 FOURIER SERIES 1 S A B A R I N A I S M A I L Outline Introduction of the Fourier series. The properties of the Fourier series. Symmetry consideration Application of the Fourier series to circuit
More informationHIGHER-ORDER LINEAR ORDINARY DIFFERENTIAL EQUATIONS II: Nonhomogeneous Equations. David Levermore Department of Mathematics University of Maryland
HIGHER-ORDER LINEAR ORDINARY DIFFERENTIAL EQUATIONS II: Nonhomogeneous Equations David Levermore Department of Mathematics University of Maryland 14 March 2012 Because the presentation of this material
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 informationSupplemental Notes on Complex Numbers, Complex Impedance, RLC Circuits, and Resonance
Supplemental Notes on Complex Numbers, Complex Impedance, RLC Circuits, and Resonance Complex numbers Complex numbers are expressions of the form z = a + ib, where both a and b are real numbers, and i
More informationCHEE 319 Tutorial 3 Solutions. 1. Using partial fraction expansions, find the causal function f whose Laplace transform. F (s) F (s) = C 1 s + C 2
CHEE 39 Tutorial 3 Solutions. Using partial fraction expansions, find the causal function f whose Laplace transform is given by: F (s) 0 f(t)e st dt (.) F (s) = s(s+) ; Solution: Note that the polynomial
More information2. Algebraic functions, power functions, exponential functions, trig functions
Math, Prep: Familiar Functions (.,.,.5, Appendix D) Name: Names of collaborators: Main Points to Review:. Functions, models, graphs, tables, domain and range. Algebraic functions, power functions, exponential
More informationLab 1: Damped, Driven Harmonic Oscillator
1 Introduction Lab 1: Damped, Driven Harmonic Oscillator The purpose of this experiment is to study the resonant properties of a driven, damped harmonic oscillator. This type of motion is characteristic
More informationCircuits Advanced Topics by Dr. Colton (Fall 2016)
ircuits Advanced Topics by Dr. olton (Fall 06). Time dependence of general and L problems General and L problems can always be cast into first order ODEs. You can solve these via the particular solution
More informationCMPT 318: Lecture 5 Complex Exponentials, Spectrum Representation
CMPT 318: Lecture 5 Complex Exponentials, Spectrum Representation Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University January 23, 2006 1 Exponentials The exponential is
More informationLab 1: damped, driven harmonic oscillator
Lab 1: damped, driven harmonic oscillator 1 Introduction The purpose of this experiment is to study the resonant properties of a driven, damped harmonic oscillator. This type of motion is characteristic
More informationSinusoids and Phasors
CHAPTER 9 Sinusoids and Phasors We now begins the analysis of circuits in which the voltage or current sources are time-varying. In this chapter, we are particularly interested in sinusoidally time-varying
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 informationPHY217: Vibrations and Waves
Assessed Problem set 1 Issued: 5 November 01 PHY17: Vibrations and Waves Deadline for submission: 5 pm Thursday 15th November, to the V&W pigeon hole in the Physics reception on the 1st floor of the GO
More informationPhysics 4 Spring 1989 Lab 5 - AC Circuits
Physics 4 Spring 1989 Lab 5 - AC Circuits Theory Consider the series inductor-resistor-capacitor circuit shown in figure 1. When an alternating voltage is applied to this circuit, the current and voltage
More informationExam II Review: Selected Solutions and Answers
November 9, 2011 Exam II Review: Selected Solutions and Answers NOTE: For additional worked problems see last year s review sheet and answers, the notes from class, and your text. Answers to problems from
More informationDr. Ian R. Manchester
Dr Ian R. Manchester Week Content Notes 1 Introduction 2 Frequency Domain Modelling 3 Transient Performance and the s-plane 4 Block Diagrams 5 Feedback System Characteristics Assign 1 Due 6 Root Locus
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 information1.5 Inverse Trigonometric Functions
1.5 Inverse Trigonometric Functions Remember that only one-to-one functions have inverses. So, in order to find the inverse functions for sine, cosine, and tangent, we must restrict their domains to intervals
More informationChapter 14 Oscillations
Chapter 14 Oscillations Chapter Goal: To understand systems that oscillate with simple harmonic motion. Slide 14-2 Chapter 14 Preview Slide 14-3 Chapter 14 Preview Slide 14-4 Chapter 14 Preview Slide 14-5
More informationTrig Identities. or (x + y)2 = x2 + 2xy + y 2. Dr. Ken W. Smith Other examples of identities are: (x + 3)2 = x2 + 6x + 9 and
Trig Identities An identity is an equation that is true for all values of the variables. Examples of identities might be obvious results like Part 4, Trigonometry Lecture 4.8a, Trig Identities and Equations
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 informationCM2202: Scientific Computing and Multimedia Applications General Maths: 3. Complex Numbers
CM2202: Scientific Computing and Multimedia Applications General Maths: 3. Complex Numbers Prof. David Marshall School of Computer Science & Informatics A problem when solving some equations There are
More informationMATHia Unit MATHia Workspace Overview TEKS
1 Function Overview Searching for Patterns Exploring and Analyzing Patterns Comparing Familiar Function Representations Students watch a video about a well-known mathematician creating an expression for
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 informationChapter 7: Trigonometric Equations and Identities
Chapter 7: Trigonometric Equations and Identities In the last two chapters we have used basic definitions and relationships to simplify trigonometric expressions and equations. In this chapter we will
More informationSolutions to the Homework Replaces Section 3.7, 3.8
Solutions to the Homework Replaces Section 3.7, 3.8 1. Our text (p. 198) states that µ ω 0 = ( 1 γ2 4km ) 1/2 1 1 2 γ 2 4km How was this approximation made? (Hint: Linearize 1 x) SOLUTION: We linearize
More informationThursday, August 4, 2011
Chapter 16 Thursday, August 4, 2011 16.1 Springs in Motion: Hooke s Law and the Second-Order ODE We have seen alrealdy that differential equations are powerful tools for understanding mechanics and electro-magnetism.
More informationAC analysis. EE 201 AC analysis 1
AC analysis Now we turn to circuits with sinusoidal sources. Earlier, we had a brief look at sinusoids, but now we will add in capacitors and inductors, making the story much more interesting. What are
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 informationSection 9.3 Phase Plane Portraits (for Planar Systems)
Section 9.3 Phase Plane Portraits (for Planar Systems) Key Terms: Equilibrium point of planer system yꞌ = Ay o Equilibrium solution Exponential solutions o Half-line solutions Unstable solution Stable
More informationMAT292 - Calculus III - Fall Solution for Term Test 2 - November 6, 2014 DO NOT WRITE ON THE QR CODE AT THE TOP OF THE PAGES.
MAT9 - Calculus III - Fall 4 Solution for Term Test - November 6, 4 Time allotted: 9 minutes. Aids permitted: None. Full Name: Last First Student ID: Email: @mail.utoronto.ca Instructions DO NOT WRITE
More informationChapter 7: Trigonometric Equations and Identities
Chapter 7: Trigonometric Equations and Identities In the last two chapters we have used basic definitions and relationships to simplify trigonometric expressions and equations. In this chapter we will
More informationThis is the number of cycles per unit time, and its units are, for example,
16 4. Sinusoidal solutions Many things in nature are periodic, even sinusoidal. We will begin by reviewing terms surrounding periodic functions. If an LTI system is fed a periodic input signal, we have
More informationName (print): Lab (circle): W8 Th8 Th11 Th2 F8. θ (radians) θ (degrees) cos θ sin θ π/ /2 1/2 π/4 45 2/2 2/2 π/3 60 1/2 3/2 π/
Name (print): Lab (circle): W8 Th8 Th11 Th2 F8 Trigonometric Identities ( cos(θ) = cos(θ) sin(θ) = sin(θ) sin(θ) = cos θ π ) 2 Cosines and Sines of common angles Euler s Formula θ (radians) θ (degrees)
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 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 informationCircuit Analysis-III. Circuit Analysis-II Lecture # 3 Friday 06 th April, 18
Circuit Analysis-III Sinusoids Example #1 ü Find the amplitude, phase, period and frequency of the sinusoid: v (t ) =12cos(50t +10 ) Signal Conversion ü From sine to cosine and vice versa. ü sin (A ± B)
More information2. Higher-order Linear ODE s
2. Higher-order Linear ODE s 2A. Second-order Linear ODE s: General Properties 2A-1. On the right below is an abbreviated form of the ODE on the left: (*) y + p()y + q()y = r() Ly = r() ; where L is the
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 informationCh 3.7: Mechanical & Electrical Vibrations
Ch 3.7: Mechanical & Electrical Vibrations Two important areas of application for second order linear equations with constant coefficients are in modeling mechanical and electrical oscillations. We will
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 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 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 information1MA1 Introduction to the Maths Course
1MA1/-1 1MA1 Introduction to the Maths Course Preamble Throughout your time as an engineering student at Oxford you will receive lectures and tuition in the range of applied mathematical tools that today
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 informationM A : Ordinary Differential Equations
M A 2 0 5 1: Ordinary Differential Equations Essential Class Notes & Graphics D 19 * 2018-2019 Sections D07 D11 & D14 1 1. INTRODUCTION CLASS 1 ODE: Course s Overarching Functions An introduction to the
More informationMATHEMATICS Lecture. 4 Chapter.8 TECHNIQUES OF INTEGRATION By Dr. Mohammed Ramidh
MATHEMATICS Lecture. 4 Chapter.8 TECHNIQUES OF INTEGRATION By TECHNIQUES OF INTEGRATION OVERVIEW The Fundamental Theorem connects antiderivatives and the definite integral. Evaluating the indefinite integral,
More informationFourier Series and Fourier Transforms
Fourier Series and Fourier Transforms EECS2 (6.082), MIT Fall 2006 Lectures 2 and 3 Fourier Series From your differential equations course, 18.03, you know Fourier s expression representing a T -periodic
More informationComplex Numbers Review
Complex Numbers view ference: Mary L. Boas, Mathematical Methods in the Physical Sciences Chapter 2 & 14 George Arfken, Mathematical Methods for Physicists Chapter 6 The real numbers (denoted R) are incomplete
More informationWave Phenomena Physics 15c
Wave Phenomena Physics 15c Lecture Harmonic Oscillators (H&L Sections 1.4 1.6, Chapter 3) Administravia! Problem Set #1! Due on Thursday next week! Lab schedule has been set! See Course Web " Laboratory
More informationMath 216 Second Midterm 28 March, 2013
Math 26 Second Midterm 28 March, 23 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 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 informationFinal Exam Review Part 2
Final Exam Review Part 2 Exponential & Logarithmic Functions and Equations Polynomial & Sinusoidal Functions Rational Expressions and Equations Exponential Functions To describe, orally and in written
More informationPortable Assisted Study Sequence ALGEBRA IIB
SCOPE This course is divided into two semesters of study (A & B) comprised of five units each. Each unit teaches concepts and strategies recommended for intermediate algebra students. The second half of
More informationThe Rayleigh Pulse Forming Network
The Rayleigh Pulse Forming Network The ideal power supply for many high current applications is one which supplies a square voltage wave over a short, but predetermined, time. Since these applications
More informationOn Exponential Decay and the Riemann Hypothesis
On Exponential Decay and the Riemann Hypothesis JEFFREY N. COOK ABSTRACT. A Riemann operator is constructed in which sequential elements are removed from a decaying set by means of prime factorization,
More information8. Introduction and Chapter Objectives
Real Analog - Circuits Chapter 8: Second Order Circuits 8. Introduction and Chapter Objectives Second order systems are, by definition, systems whose input-output relationship is a second order differential
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 informationWave Phenomena Physics 15c
Wave Phenomena Physics 15c Lecture Harmonic Oscillators (H&L Sections 1.4 1.6, Chapter 3) Administravia! Problem Set #1! Due on Thursday next week! Lab schedule has been set! See Course Web " Laboratory
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 informationChapter 7: Trigonometric Equations and Identities
Section 7. Solving Trigonometric Equations and Identities 5 Chapter 7: Trigonometric Equations and Identities In the last two chapters we have used basic definitions and relationships to simplify trigonometric
More informationSecond-Order Linear ODEs
Chap. 2 Second-Order Linear ODEs Sec. 2.1 Homogeneous Linear ODEs of Second Order On pp. 45-46 we extend concepts defined in Chap. 1, notably solution and homogeneous and nonhomogeneous, to second-order
More informationSTABILITY. Have looked at modeling dynamic systems using differential equations. and used the Laplace transform to help find step and impulse
SIGNALS AND SYSTEMS: PAPER 3C1 HANDOUT 4. Dr David Corrigan 1. Electronic and Electrical Engineering Dept. corrigad@tcd.ie www.sigmedia.tv STABILITY Have looked at modeling dynamic systems using differential
More information11. The Series RLC Resonance Circuit
Electronicsab.nb. The Series RC Resonance Circuit Introduction Thus far we have studied a circuit involving a () series resistor R and capacitor C circuit as well as a () series resistor R and inductor
More informationExperiment 3: Resonance in LRC Circuits Driven by Alternating Current
Experiment 3: Resonance in LRC Circuits Driven by Alternating Current Introduction In last week s laboratory you examined the LRC circuit when constant voltage was applied to it. During this laboratory
More informationCHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION. Professor Dae Ryook Yang
CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dept. of Chemical and Biological Engineering 5-1 Road Map of the Lecture V Laplace Transform and Transfer
More informationPHYS 3900 Homework Set #02
PHYS 3900 Homework Set #02 Part = HWP 2.0, 2.02, 2.03. Due: Mon. Jan. 22, 208, 4:00pm Part 2 = HWP 2.04, 2.05, 2.06. Due: Fri. Jan. 26, 208, 4:00pm All textbook problems assigned, unless otherwise stated,
More informationElectric Circuit Theory
Electric Circuit Theory Nam Ki Min nkmin@korea.ac.kr 010-9419-2320 Chapter 11 Sinusoidal Steady-State Analysis Nam Ki Min nkmin@korea.ac.kr 010-9419-2320 Contents and Objectives 3 Chapter Contents 11.1
More informationAMATH 351 Mar 15, 2013 FINAL REVIEW. Instructor: Jiri Najemnik
AMATH 351 Mar 15, 013 FINAL REVIEW Instructor: Jiri Najemni ABOUT GRADES Scores I have so far will be posted on the website today sorted by the student number HW4 & Exam will be added early next wee Let
More information