Hamiltonian systems. Marc R. Roussel. October 25, 2005
|
|
- Dorothy Jocelyn Hamilton
- 5 years ago
- Views:
Transcription
1 Hamiltoia systems Marc R. Roussel October 25, 25 1 Itroductio Today s otes will deviate somewhat from the mai lie of lectures to itroduce a importat class of dyamical systems which were first studied i mechaics, amely Hamiltoia systems. There is a large literature o Hamiltoia systems. The itetio here is ot to comprehesively survey this literature, which would be quite impossible eve if we devoted a whole course to the subject, but to discuss some of the key ideas uderlyig the treatmet of this class of systems. We begi with some defiitios: Defiitio 1 Suppose that H(x,p) is a smooth fuctio of its argumets for x ad p R. The the dyamical system ẋ i = H p i, (1a) ṗ i = H x i (i = 1,2,...,) is called a Hamiltoia system ad H is the Hamiltoia fuctio (or just the Hamiltoia) of the system. Equatios 1 are called Hamilto s equatios. Defiitio 2 The umber of degrees of freedom of a Hamiltoia system is the umber of (x i, p i ) pairs i Hamilto s equatios, i.e. the value of. The phase space is therefore 2-dimesioal. I mechaics, the vector x represets the geeralized coordiates of the compoets of the system (positios, agles, etc.), while p is a set of geeralized mometa. There is a elaborate theory for costructig geeralized mometa ad Hamiltoias which we wo t go ito here. Note that the geeralized mometa are ot always ordiary liear mometa (mv i ), although that is ofte the case. Note that the Hamiltoia fuctio is a cotat of the motio: Ḣ = = H ẋ i + x i H H x i p i + H p i ṗ i 1 H p i ( H ) =. x i (1b)
2 E (J) p (kg m/s) x (m) -1 Figure 1: Harmoic oscillator orbits for k = 2 N/m, m =.1 kg, ad various values of the total eergy E. Maple code which ca be used to draw such a figure is give i Appedix A. (I actually used guplot because I thik the output looks a little icer, but the priciple is the same.) Note that these closed orbits are ot limit cycles sice differet iitial coditios (which i geeral lead to differet values of E) will produce differet orbits. I fact, the Hamiltoia is ofte just the total eergy i mechaical systems, although this is t always the case. Let us for the momet specialize the discussio to plaar systems, i.e. systems for which = 1. The fact that H is costat is meas that the motio is costraied to the curve H(x, p) = h, where h is the value of the Hamiltoia fuctio implied by the iitial coditios. This curve is of course just the orbit of the system i phase space. If we chage the total eergy, we get a differet orbit. Example 1.1 A harmoic oscillator is a mass-sprig system with potetial eergy 1 2 kx2, where x is the displacemet of the sprig from equilibrium. For simple systems like this oe i which the potetial eergy simply depeds o the positio, the Hamiltoia is just the total eergy: H(x, p) = 1 2 kx2 + p2 2m, (2) where p is the mometum. Because H is a costat, the orbits are just the family of ellipses 1 2 kx2 + p2 2m = E. The value of E is fixed by the iitial coditios. Differet values of E correspod to ellipses of differet sizes. Figure 1 shows a example. 2
3 If we are iterested i the equatios of motio, we ca recover them from Hamilto s equatios: ẋ = H p = p/m. ṗ = H x = kx. Sice v = p/m, the first equatio just says that ẋ is the velocity of the particle. Recall that F = ṗ. (We usually write this as F = ma i elemetary physics courses.) Thus, the secod equatio gives us the force law correspodig to the potetial eergy 1 2 kx2, which you may recogize from your earlier studies of physics. As you well kow from mechaics, mechaical systems are subject to other coservatio laws: mometum, agular mometum, ad so o. If we ca fid eough coservatio laws for a system the, as i the example above, the busiess of determiig the behavior of the system becomes rather simple. O the other had, if we have more degrees of freedom tha costats, the behavior ca be complex. These issues are discussed briefly i the followig sectios. (3a) (3b) 2 Itegrable systems We start with some specialized machiery ad defiitios pertaiig to Hamiltoia systems. Defiitio 3 Let H(x,p) ad L(x,p) be differetiable fuctios of their argumets for x ad p R. The Poisso bracket of H with L, {H,L}, is defied by {H,L} = ( H L H ) L. p i x i x i p i Defiitio 4 A quatity L is called a first itegral of a Hamiltoia system if it is a costat of the motio, i.e. if L = uder the flow implied by Hamilto s equatios. Note that the Hamiltoia itself is a first itegral accordig to this defiitio. Theorem 1 The quatity L is a first itegral of a Hamiltoia system with Hamiltoia H if {H,L} =. The proof of this propositio is similar to the proof that H is a costat, ad is left as a exercise to the reader. Defiitio 5 A Hamiltoia system is said to be completely itegrable if it has first itegrals (icludig the Hamiltoia itself), where is the umber of degrees of freedom. I mechaical systems, the first itegrals are ofte familiar quatities. 3
4 Example 2.1 Cosider the two-body orbital problem with m 1 m 2. I this case, we ca treat the body of mass m 1 as beig essetially motioless. 1 Moreover, with just two bodies, the coordiate system ca always be chose i such a way that m 1 is at the origi, ad that z =, i.e. the motio is cofied to the (x,y) plae. The Hamiltoia is H = p2 x + p2 y 2m 2 Gm 1m 2 x 2 + y 2. Based o our kowledge of mechaics, we expect the agular mometum L = r p = xp y yp x to be a costat of the motio. If this is so, the Poisso bracket of H with L should be zero: {H,L} = H L p x x H L + H L x p x p y y H L y p y = p x p y Gm 1m 2 x m 2 (x 2 + y 2 ) 3/2 ( y) + p y ( p x ) Gm 1m 2 y m 2 (x 2 + y 2 ) 3/2 x =. The agular mometum is therefore, as expected, a costat of the motio. This meas that the system is cofied to the itersectio of the hyperplaes H = E ad L = l, where E ad l are costats determied by the iitial coditios. This itersectio is a ivariat maifold of the differetial equatios sice all possible trajectories for a give E ad l lie withi it. What does this maifold look like? We have two equatios coectig our four variables, so these equatios implicitly defie a two-dimesioal surface i the four-dimesioal phase space. I this case, this surface is easy to plot usig Maple. See Appedix B for sample code. We start by solvig for y (say) from the agular mometum: y = (xp y l)/p x. We the substitute this expressio ito the eergy equatio. Depedig o the total eergy, we get differet kids of results: 1. The surface is a torus if E < (Fig. 2). The trajectories o the torus ca either be periodic (which they are for this simple problem) or quasiperiodic, i.e. they wid aroud ad aroud the torus without repeatig themselves. I the latter case, a trajectory evetually covers the whole torus ad is said to be ergodic. I the former case, slightly differet iitial coditios give differet orbits. 2. The surface is essetially a cylider if E > (Fig. 3) or if E = (Fig. 4). This correspods to a case where the lighter particle has eough eergy to escape the gravitatioal attractio of its heavier eighbor. The trajectories wid aroud this cylider, either toward x = or toward x =. 1 A more careful treatmet would trasform to cetre-of-mass ad relative coordiates. We would obtai essetially the same Hamiltoia, with a small adjustmet i the costats appearig i the equatio. 4
5 4 2 x px.5.5 py Figure 2: Ivariat torus arisig from the sigle-plaet orbital problem. This case correspods to m 1 = 1kg, m 2 = 1kg, E = J, ad l = 1kgm 2 /s. The foregoig example brigs out importat features of itegrable Hamiltoia problems: Ivariat tori or cyliders are geerically see. The study of these ivariat maifolds ad how they chage whe we chage the Hamiltoia is a importat part of the theory of Hamiltoia dyamical systems. The study of trajectories o these maifolds is of course also iterestig ad importat. I the case of itegrable Hamiltoias, the behavior ca t be much more complicated tha a quasiperiodic trajectory. 3 Numerical itegratio We wo t go ito great detail, but there is a importat issue which we must discuss i relatio to the umerical itegratio of Hamiltoia systems: Most umerical methods do ot preserve the costacy of H. May umerical methods do a good job of holdig H reasoably costat for simple problems, but ot for more complex problems, especially if log trajectories are wated. Numerical itegratio methods which boud the variatio of H are said to be symplectic. We will discuss this problem ad the costructio of symplectic itegrators i this sectio. However, we will ot udertake a full implemetatio, the details beig rather messy. As a itroductio to this topic, let us discuss the simplest umerical itegratio method, amely the forward Euler method. The idea behid this algorithm is very simple: The differetial 5
6 1 5 x py px Figure 3: Ivariat cylider arisig from the sigle-plaet orbital problem. All parameters are as i Fig. 2, except E = J ad l = 12kgm 2 /s. equatio dz = f(z,t), (4) dt where z is the vector of phase-space coordiates, ca be iterpreted as follows, accordig to the basic defiitio of the derivative: z(t + ) z(t) lim = f(z(t),t). Euler s method simply cosists i takig small values of ad doig away with the limit. Sice is fixed, we will be calculatig z at followig sequece of times:,2,3,..., j,... We therefore label the values of z by j, the umber of multiples of the step size which have passed. Euler s method ca thus be summarized by the formula z j+1 = z j + f(z j, j). (5) Equatio 5 defies a dyamical system with discrete time (idexed by j) i which the state at the ext time icremet is calculated from a simple formula. This is called a map. I the limit, the Euler map coverges o the solutios of the differetial equatio 4. Hopefully, for small values of, Euler s method will give results which are at least reasoably represetative of those of the origial differetial equatio. 6
7 1 5 x px py Figure 4: Ivariat surface arisig from the sigle-plaet orbital problem with E = ad l = 115kgm 2 /s, all other parameters beig set as i Fig. 2. As you ca imagie, Euler s method is a very aïve techique, ad teds ot to be very accurate. Nevertheless, it has its uses, amog them pedagogical applicatios. We will use it here to illustrate the problem you ru ito with ordiary umerical methods ad Hamiltoia systems. Example 3.1 Suppose that we wat to solve the differetial equatios 3 which arise from the Hamiltoia 2 by Euler s method. The we have x j+1 = x j + p j /m, p j+1 = p j kx j. The value of the Hamiltoia at time idex j + 1 is H(x j+1, p j+1 ) = 1 2 kx2 j+1 + p2 j+1 2m = k ( x j + p j /m ) 2 1 ( ) 2 + p j kx j 2 2m ( ) = 1 2 kx2 j + p2 j 2m + k 1 m ()2 2 kx2 j + p2 j 2m = H(x j, p j ) (1 + km ) ()2. 7
8 The value of the Hamiltoia therefore grows by a factor of 1 + k m ()2 at every step. The problem is less serious the smaller we make, but it remais that it is ofte just ot acceptable to have H grow (or decrease) with time. Most of the commo umerical methods you may have heard of (Ruge-Kutta, Gear, etc.) are ot symplectic. For some problems, some of these methods give reasoable results. However, if you rely o these methods for Hamiltoia systems, you will sooer or later get bured. Euler s method is a explicit method: Give the curret poit i phase space, we ca calculate the ext poit just by pluggig umbers ito a formula. Ufortuately, all geeral-purpose symplectic itegratio methods are implicit: We ca t i geeral write dow a formula for the ext poit. Rather, we have to solve a equatio (usually umerically) to calculate (x j+1,p j+1 ). There are however some explicit methods for special forms of the Hamiltoia which we will ot discuss here. The key to maitaiig the costacy of the Hamiltoia is to derive a umerical method directly from Hamilto s equatios 1. Let s start with equatio 1a. We ca estimate the time derivative o the left-had side as we did i Euler s method: dx i dt x i, j+1 x i, j. I this equatio, x i, j+1 is the value of x i at time step j. The right-had side of equatio 1a is also a derivative. I each time step, we will icremet the variables, for istace icremetig p i from p i, j to p i, j+1. Sice we oly take small time steps, p i should oly chage by a small amout. This gives us a opportuity to evaluate H/ p i by usig aother fiite-differece approximatio. You ca probably imagie what we eed to do, but it s a bit of a mess to write dow without some special otatio. Let p j p i, j+1 be the vector obtaied by replacig p i, j by p i, j+1 i the vector of mometa at time step k, p j. The, Puttig the two sides together, we have H H(x j,p j p i, j+1 ) H(x j,p j ). p i p i, j+1 p i, j x i, j+1 x i, j = H(x j,p j p i, j+1 ) H(x j,p j ) p i, j+1 p i, j. (6a) We ca repeat this procedure for equatio 1b, with oe mior variatio: p i, j+1 p i, j = H(x j x i, j+1,p j p i, j+1 ) H(x j,p j p i, j+1 ) x i, j+1 x i, j. (6b) The reaso that we used the updated mometum i equatio 6b is that this decisio makes the scheme symplectic. To see this, rearrage equatios 6a ad 6b to ad 1 ( )( ) xi, j+1 x i, j pi, j+1 p i, j = H(x j,p j p i, j+1 ) H(x j,p j ), 1 ( )( ) [ xi, j+1 x i, j pi, j+1 p i, j = H(x j x i, j+1,p j p i, j+1 ) H(x j,p j p i, j+1 ) ]. 8
9 Now subtract these two equatios: or = H(x j x i, j+1,p j p i, j+1 ) H(x j,p j ), H(x j x i, j+1,p j p i, j+1 ) = H(x j,p j ). I other words, the update rule cosistig of equatios 6 keeps the value of H costat for each pair of cojugate variables (x i, p i ) to which it is applied. Thus, the method is symplectic. The tradeoff is that we have a implicit method: We eed values of the Hamiltoia with some of the coordiates evaluated at the ukow ext time poit. We therefore have 2 equatios i 2 ukows. Uless the Hamiltoia is of a particularly simple form, we wo t be able to reduce this problem to a simple mappig. Not all symplectic itegratio methods hold the value of the Hamiltoia exactly costat. Rather, their effect is to limit the variatio i the value of the Hamiltoia. The drift i the Hamiltoia ecoutered with ordiary umerical methods ca be disastrous. Holdig the Hamiltoia withi tight bouds is ofte eough to give good results, ad oe is ofte cotet to trade off exact costacy of the Hamiltoia agaist step size i order to accelerate log itegratios. Appedices A Maple code for drawig the orbits of the harmoic oscillator k := 2; m :=.1; H := (x,p) -> pˆ2/(2*m) + k*xˆ2/2; with(plots): cotourplot(h(x,p),x= ,p= ,cotours=[$1..6]); B Maple code for drawig the ivariat tori or cyliders of Hamiltoias with two degrees of freedom y := (x*py - ell)/px; # Calculate y from agular mometum H := (x,y,px,py) -> (pxˆ2+pyˆ2)/(2*m2) - G*m1*m2/sqrt(xˆ2+yˆ2); m1 := 1e11; # A slightly differet set of masses tha i the otes m2 := 1/1e6; G := 6.67e-11; x := 1; # Iitial poit (x,y,px,py) = (x,,,p) p :=.1/sqrt(1e6); E := H(x,,,p); # Set the eergy ell := x*p; # Set the agular mometum with(plots): implicitplot3d(h(x,y,px,py)=e,px=-p..p,py=-p..p,x=-5*x..5*x, axes=boxed,grid=[1,4,4],orietatio=[-35,5],shadig=none); 9
Physics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001.
Physics 324, Fall 2002 Dirac Notatio These otes were produced by David Kapla for Phys. 324 i Autum 2001. 1 Vectors 1.1 Ier product Recall from liear algebra: we ca represet a vector V as a colum vector;
More information1. Linearization of a nonlinear system given in the form of a system of ordinary differential equations
. Liearizatio of a oliear system give i the form of a system of ordiary differetial equatios We ow show how to determie a liear model which approximates the behavior of a time-ivariat oliear system i a
More informationx a x a Lecture 2 Series (See Chapter 1 in Boas)
Lecture Series (See Chapter i Boas) A basic ad very powerful (if pedestria, recall we are lazy AD smart) way to solve ay differetial (or itegral) equatio is via a series expasio of the correspodig solutio
More informationInfinite Sequences and Series
Chapter 6 Ifiite Sequeces ad Series 6.1 Ifiite Sequeces 6.1.1 Elemetary Cocepts Simply speakig, a sequece is a ordered list of umbers writte: {a 1, a 2, a 3,...a, a +1,...} where the elemets a i represet
More informationOnce we have a sequence of numbers, the next thing to do is to sum them up. Given a sequence (a n ) n=1
. Ifiite Series Oce we have a sequece of umbers, the ext thig to do is to sum them up. Give a sequece a be a sequece: ca we give a sesible meaig to the followig expressio? a = a a a a While summig ifiitely
More information6.3 Testing Series With Positive Terms
6.3. TESTING SERIES WITH POSITIVE TERMS 307 6.3 Testig Series With Positive Terms 6.3. Review of what is kow up to ow I theory, testig a series a i for covergece amouts to fidig the i= sequece of partial
More informationmx bx kx F t. dt IR I LI V t, Q LQ RQ V t,
Lecture 5 omplex Variables II (Applicatios i Physics) (See hapter i Boas) To see why complex variables are so useful cosider first the (liear) mechaics of a sigle particle described by Newto s equatio
More informationChapter 10: Power Series
Chapter : Power Series 57 Chapter Overview: Power Series The reaso series are part of a Calculus course is that there are fuctios which caot be itegrated. All power series, though, ca be itegrated because
More informationMath 2784 (or 2794W) University of Connecticut
ORDERS OF GROWTH PAT SMITH Math 2784 (or 2794W) Uiversity of Coecticut Date: Mar. 2, 22. ORDERS OF GROWTH. Itroductio Gaiig a ituitive feel for the relative growth of fuctios is importat if you really
More informationCHAPTER 10 INFINITE SEQUENCES AND SERIES
CHAPTER 10 INFINITE SEQUENCES AND SERIES 10.1 Sequeces 10.2 Ifiite Series 10.3 The Itegral Tests 10.4 Compariso Tests 10.5 The Ratio ad Root Tests 10.6 Alteratig Series: Absolute ad Coditioal Covergece
More informationSequences, Mathematical Induction, and Recursion. CSE 2353 Discrete Computational Structures Spring 2018
CSE 353 Discrete Computatioal Structures Sprig 08 Sequeces, Mathematical Iductio, ad Recursio (Chapter 5, Epp) Note: some course slides adopted from publisher-provided material Overview May mathematical
More information1 Approximating Integrals using Taylor Polynomials
Seughee Ye Ma 8: Week 7 Nov Week 7 Summary This week, we will lear how we ca approximate itegrals usig Taylor series ad umerical methods. Topics Page Approximatig Itegrals usig Taylor Polyomials. Defiitios................................................
More informationMath 257: Finite difference methods
Math 257: Fiite differece methods 1 Fiite Differeces Remember the defiitio of a derivative f f(x + ) f(x) (x) = lim 0 Also recall Taylor s formula: (1) f(x + ) = f(x) + f (x) + 2 f (x) + 3 f (3) (x) +...
More information1 Adiabatic and diabatic representations
1 Adiabatic ad diabatic represetatios 1.1 Bor-Oppeheimer approximatio The time-idepedet Schrödiger equatio for both electroic ad uclear degrees of freedom is Ĥ Ψ(r, R) = E Ψ(r, R), (1) where the full molecular
More informationChapter 5 Vibrational Motion
Fall 4 Chapter 5 Vibratioal Motio... 65 Potetial Eergy Surfaces, Rotatios ad Vibratios... 65 Harmoic Oscillator... 67 Geeral Solutio for H.O.: Operator Techique... 68 Vibratioal Selectio Rules... 7 Polyatomic
More informationChapter 6 Infinite Series
Chapter 6 Ifiite Series I the previous chapter we cosidered itegrals which were improper i the sese that the iterval of itegratio was ubouded. I this chapter we are goig to discuss a topic which is somewhat
More informationKinetics of Complex Reactions
Kietics of Complex Reactios by Flick Colema Departmet of Chemistry Wellesley College Wellesley MA 28 wcolema@wellesley.edu Copyright Flick Colema 996. All rights reserved. You are welcome to use this documet
More informationsubcaptionfont+=small,labelformat=parens,labelsep=space,skip=6pt,list=0,hypcap=0 subcaption ALGEBRAIC COMBINATORICS LECTURE 8 TUESDAY, 2/16/2016
subcaptiofot+=small,labelformat=pares,labelsep=space,skip=6pt,list=0,hypcap=0 subcaptio ALGEBRAIC COMBINATORICS LECTURE 8 TUESDAY, /6/06. Self-cojugate Partitios Recall that, give a partitio λ, we may
More informationMath 312 Lecture Notes One Dimensional Maps
Math 312 Lecture Notes Oe Dimesioal Maps Warre Weckesser Departmet of Mathematics Colgate Uiversity 21-23 February 25 A Example We begi with the simplest model of populatio growth. Suppose, for example,
More informationMath 113 Exam 3 Practice
Math Exam Practice Exam will cover.-.9. This sheet has three sectios. The first sectio will remid you about techiques ad formulas that you should kow. The secod gives a umber of practice questios for you
More information62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +
62. Power series Defiitio 16. (Power series) Give a sequece {c }, the series c x = c 0 + c 1 x + c 2 x 2 + c 3 x 3 + is called a power series i the variable x. The umbers c are called the coefficiets of
More informationChapter 4. Fourier Series
Chapter 4. Fourier Series At this poit we are ready to ow cosider the caoical equatios. Cosider, for eample the heat equatio u t = u, < (4.) subject to u(, ) = si, u(, t) = u(, t) =. (4.) Here,
More informationINFINITE SEQUENCES AND SERIES
11 INFINITE SEQUENCES AND SERIES INFINITE SEQUENCES AND SERIES 11.4 The Compariso Tests I this sectio, we will lear: How to fid the value of a series by comparig it with a kow series. COMPARISON TESTS
More informationLecture 6: Integration and the Mean Value Theorem. slope =
Math 8 Istructor: Padraic Bartlett Lecture 6: Itegratio ad the Mea Value Theorem Week 6 Caltech 202 The Mea Value Theorem The Mea Value Theorem abbreviated MVT is the followig result: Theorem. Suppose
More informationChapter 9: Numerical Differentiation
178 Chapter 9: Numerical Differetiatio Numerical Differetiatio Formulatio of equatios for physical problems ofte ivolve derivatives (rate-of-chage quatities, such as velocity ad acceleratio). Numerical
More informationSECTION 2 Electrostatics
SECTION Electrostatics This sectio, based o Chapter of Griffiths, covers effects of electric fields ad forces i static (timeidepedet) situatios. The topics are: Electric field Gauss s Law Electric potetial
More informationStatistics 511 Additional Materials
Cofidece Itervals o mu Statistics 511 Additioal Materials This topic officially moves us from probability to statistics. We begi to discuss makig ifereces about the populatio. Oe way to differetiate probability
More informationThe Born-Oppenheimer approximation
The Bor-Oppeheimer approximatio 1 Re-writig the Schrödiger equatio We will begi from the full time-idepedet Schrödiger equatio for the eigestates of a molecular system: [ P 2 + ( Pm 2 + e2 1 1 2m 2m m
More informationRecursive Algorithms. Recurrences. Recursive Algorithms Analysis
Recursive Algorithms Recurreces Computer Sciece & Egieerig 35: Discrete Mathematics Christopher M Bourke cbourke@cseuledu A recursive algorithm is oe i which objects are defied i terms of other objects
More informationMA131 - Analysis 1. Workbook 3 Sequences II
MA3 - Aalysis Workbook 3 Sequeces II Autum 2004 Cotets 2.8 Coverget Sequeces........................ 2.9 Algebra of Limits......................... 2 2.0 Further Useful Results........................
More informationMa 530 Infinite Series I
Ma 50 Ifiite Series I Please ote that i additio to the material below this lecture icorporated material from the Visual Calculus web site. The material o sequeces is at Visual Sequeces. (To use this li
More informationDiscrete Mathematics for CS Spring 2008 David Wagner Note 22
CS 70 Discrete Mathematics for CS Sprig 2008 David Wager Note 22 I.I.D. Radom Variables Estimatig the bias of a coi Questio: We wat to estimate the proportio p of Democrats i the US populatio, by takig
More information6 Integers Modulo n. integer k can be written as k = qn + r, with q,r, 0 r b. So any integer.
6 Itegers Modulo I Example 2.3(e), we have defied the cogruece of two itegers a,b with respect to a modulus. Let us recall that a b (mod ) meas a b. We have proved that cogruece is a equivalece relatio
More information3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,
3. Z Trasform Referece: Etire Chapter 3 of text. Recall that the Fourier trasform (FT) of a DT sigal x [ ] is ω ( ) [ ] X e = j jω k = xe I order for the FT to exist i the fiite magitude sese, S = x [
More informationThe picture in figure 1.1 helps us to see that the area represents the distance traveled. Figure 1: Area represents distance travelled
1 Lecture : Area Area ad distace traveled Approximatig area by rectagles Summatio The area uder a parabola 1.1 Area ad distace Suppose we have the followig iformatio about the velocity of a particle, how
More informationAppendix: The Laplace Transform
Appedix: The Laplace Trasform The Laplace trasform is a powerful method that ca be used to solve differetial equatio, ad other mathematical problems. Its stregth lies i the fact that it allows the trasformatio
More informationSequences. Notation. Convergence of a Sequence
Sequeces A sequece is essetially just a list. Defiitio (Sequece of Real Numbers). A sequece of real umbers is a fuctio Z (, ) R for some real umber. Do t let the descriptio of the domai cofuse you; it
More informationCOURSE INTRODUCTION: WHAT HAPPENS TO A QUANTUM PARTICLE ON A PENDULUM π 2 SECONDS AFTER IT IS TOSSED IN?
COURSE INTRODUCTION: WHAT HAPPENS TO A QUANTUM PARTICLE ON A PENDULUM π SECONDS AFTER IT IS TOSSED IN? DROR BAR-NATAN Follows a lecture give by the author i the trivial otios semiar i Harvard o April 9,
More informationConvergence of random variables. (telegram style notes) P.J.C. Spreij
Covergece of radom variables (telegram style otes).j.c. Spreij this versio: September 6, 2005 Itroductio As we kow, radom variables are by defiitio measurable fuctios o some uderlyig measurable space
More information2.4 - Sequences and Series
2.4 - Sequeces ad Series Sequeces A sequece is a ordered list of elemets. Defiitio 1 A sequece is a fuctio from a subset of the set of itegers (usually either the set 80, 1, 2, 3,... < or the set 81, 2,
More informationCHAPTER I: Vector Spaces
CHAPTER I: Vector Spaces Sectio 1: Itroductio ad Examples This first chapter is largely a review of topics you probably saw i your liear algebra course. So why cover it? (1) Not everyoe remembers everythig
More information1.3 Convergence Theorems of Fourier Series. k k k k. N N k 1. With this in mind, we state (without proof) the convergence of Fourier series.
.3 Covergece Theorems of Fourier Series I this sectio, we preset the covergece of Fourier series. A ifiite sum is, by defiitio, a limit of partial sums, that is, a cos( kx) b si( kx) lim a cos( kx) b si(
More informationLecture 2: Monte Carlo Simulation
STAT/Q SCI 43: Itroductio to Resamplig ethods Sprig 27 Istructor: Ye-Chi Che Lecture 2: ote Carlo Simulatio 2 ote Carlo Itegratio Assume we wat to evaluate the followig itegratio: e x3 dx What ca we do?
More informationTrue Nature of Potential Energy of a Hydrogen Atom
True Nature of Potetial Eergy of a Hydroge Atom Koshu Suto Key words: Bohr Radius, Potetial Eergy, Rest Mass Eergy, Classical Electro Radius PACS codes: 365Sq, 365-w, 33+p Abstract I cosiderig the potetial
More informationINFINITE SEQUENCES AND SERIES
INFINITE SEQUENCES AND SERIES INFINITE SEQUENCES AND SERIES I geeral, it is difficult to fid the exact sum of a series. We were able to accomplish this for geometric series ad the series /[(+)]. This is
More informationCALCULATION OF FIBONACCI VECTORS
CALCULATION OF FIBONACCI VECTORS Stuart D. Aderso Departmet of Physics, Ithaca College 953 Daby Road, Ithaca NY 14850, USA email: saderso@ithaca.edu ad Dai Novak Departmet of Mathematics, Ithaca College
More informationSimilarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall
Iteratioal Mathematical Forum, Vol. 9, 04, o. 3, 465-475 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/imf.04.48 Similarity Solutios to Usteady Pseudoplastic Flow Near a Movig Wall W. Robi Egieerig
More informationPHYC - 505: Statistical Mechanics Homework Assignment 4 Solutions
PHYC - 55: Statistical Mechaics Homewor Assigmet 4 Solutios Due February 5, 14 1. Cosider a ifiite classical chai of idetical masses coupled by earest eighbor sprigs with idetical sprig costats. a Write
More informationSeptember 2012 C1 Note. C1 Notes (Edexcel) Copyright - For AS, A2 notes and IGCSE / GCSE worksheets 1
September 0 s (Edecel) Copyright www.pgmaths.co.uk - For AS, A otes ad IGCSE / GCSE worksheets September 0 Copyright www.pgmaths.co.uk - For AS, A otes ad IGCSE / GCSE worksheets September 0 Copyright
More informationProduct measures, Tonelli s and Fubini s theorems For use in MAT3400/4400, autumn 2014 Nadia S. Larsen. Version of 13 October 2014.
Product measures, Toelli s ad Fubii s theorems For use i MAT3400/4400, autum 2014 Nadia S. Larse Versio of 13 October 2014. 1. Costructio of the product measure The purpose of these otes is to preset the
More information(A sequence also can be thought of as the list of function values attained for a function f :ℵ X, where f (n) = x n for n 1.) x 1 x N +k x N +4 x 3
MATH 337 Sequeces Dr. Neal, WKU Let X be a metric space with distace fuctio d. We shall defie the geeral cocept of sequece ad limit i a metric space, the apply the results i particular to some special
More informationNumerical Methods for Ordinary Differential Equations
Numerical Methods for Ordiary Differetial Equatios Braislav K. Nikolić Departmet of Physics ad Astroomy, Uiversity of Delaware, U.S.A. PHYS 460/660: Computatioal Methods of Physics http://www.physics.udel.edu/~bikolic/teachig/phys660/phys660.html
More informationEDEXCEL NATIONAL CERTIFICATE UNIT 4 MATHEMATICS FOR TECHNICIANS OUTCOME 4 - CALCULUS
EDEXCEL NATIONAL CERTIFICATE UNIT 4 MATHEMATICS FOR TECHNICIANS OUTCOME 4 - CALCULUS TUTORIAL 1 - DIFFERENTIATION Use the elemetary rules of calculus arithmetic to solve problems that ivolve differetiatio
More informationWe are mainly going to be concerned with power series in x, such as. (x)} converges - that is, lims N n
Review of Power Series, Power Series Solutios A power series i x - a is a ifiite series of the form c (x a) =c +c (x a)+(x a) +... We also call this a power series cetered at a. Ex. (x+) is cetered at
More informationAnalytical mechanics
Divisio of Mechaics Lud Uiversity Aalytical mechaics 7059 Name (write i block letters):. Id.-umber: Writte examiatio with five tasks. Please check that all tasks are icluded. A clea copy of the solutios
More informationSeries III. Chapter Alternating Series
Chapter 9 Series III With the exceptio of the Null Sequece Test, all the tests for series covergece ad divergece that we have cosidered so far have dealt oly with series of oegative terms. Series with
More informationCastiel, Supernatural, Season 6, Episode 18
13 Differetial Equatios the aswer to your questio ca best be epressed as a series of partial differetial equatios... Castiel, Superatural, Seaso 6, Episode 18 A differetial equatio is a mathematical equatio
More informationMicroscopic Theory of Transport (Fall 2003) Lecture 6 (9/19/03) Static and Short Time Properties of Time Correlation Functions
.03 Microscopic Theory of Trasport (Fall 003) Lecture 6 (9/9/03) Static ad Short Time Properties of Time Correlatio Fuctios Refereces -- Boo ad Yip, Chap There are a umber of properties of time correlatio
More informationSeunghee Ye Ma 8: Week 5 Oct 28
Week 5 Summary I Sectio, we go over the Mea Value Theorem ad its applicatios. I Sectio 2, we will recap what we have covered so far this term. Topics Page Mea Value Theorem. Applicatios of the Mea Value
More informationDiscrete-Time Systems, LTI Systems, and Discrete-Time Convolution
EEL5: Discrete-Time Sigals ad Systems. Itroductio I this set of otes, we begi our mathematical treatmet of discrete-time s. As show i Figure, a discrete-time operates or trasforms some iput sequece x [
More informationNotes on iteration and Newton s method. Iteration
Notes o iteratio ad Newto s method Iteratio Iteratio meas doig somethig over ad over. I our cotet, a iteratio is a sequece of umbers, vectors, fuctios, etc. geerated by a iteratio rule of the type 1 f
More informationA sequence of numbers is a function whose domain is the positive integers. We can see that the sequence
Sequeces A sequece of umbers is a fuctio whose domai is the positive itegers. We ca see that the sequece,, 2, 2, 3, 3,... is a fuctio from the positive itegers whe we write the first sequece elemet as
More informationLecture 7: Density Estimation: k-nearest Neighbor and Basis Approach
STAT 425: Itroductio to Noparametric Statistics Witer 28 Lecture 7: Desity Estimatio: k-nearest Neighbor ad Basis Approach Istructor: Ye-Chi Che Referece: Sectio 8.4 of All of Noparametric Statistics.
More information4.3 Growth Rates of Solutions to Recurrences
4.3. GROWTH RATES OF SOLUTIONS TO RECURRENCES 81 4.3 Growth Rates of Solutios to Recurreces 4.3.1 Divide ad Coquer Algorithms Oe of the most basic ad powerful algorithmic techiques is divide ad coquer.
More informationMechanics Physics 151
Mechaics Physics 151 Lecture 4 Cotiuous Systems ad Fields (Chapter 13) What We Did Last Time Built Lagragia formalism for cotiuous system Lagragia L = L dxdydz d L L Lagrage s equatio = dx η, η Derived
More information10.6 ALTERNATING SERIES
0.6 Alteratig Series Cotemporary Calculus 0.6 ALTERNATING SERIES I the last two sectios we cosidered tests for the covergece of series whose terms were all positive. I this sectio we examie series whose
More information1 6 = 1 6 = + Factorials and Euler s Gamma function
Royal Holloway Uiversity of Lodo Departmet of Physics Factorials ad Euler s Gamma fuctio Itroductio The is a self-cotaied part of the course dealig, essetially, with the factorial fuctio ad its geeralizatio
More informationNumerical Methods in Fourier Series Applications
Numerical Methods i Fourier Series Applicatios Recall that the basic relatios i usig the Trigoometric Fourier Series represetatio were give by f ( x) a o ( a x cos b x si ) () where the Fourier coefficiets
More informationSequences A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence
Sequeces A sequece of umbers is a fuctio whose domai is the positive itegers. We ca see that the sequece 1, 1, 2, 2, 3, 3,... is a fuctio from the positive itegers whe we write the first sequece elemet
More informationChapter 2: Numerical Methods
Chapter : Numerical Methods. Some Numerical Methods for st Order ODEs I this sectio, a summar of essetial features of umerical methods related to solutios of ordiar differetial equatios is give. I geeral,
More informationAssignment 2 Solutions SOLUTION. ϕ 1 Â = 3 ϕ 1 4i ϕ 2. The other case can be dealt with in a similar way. { ϕ 2 Â} χ = { 4i ϕ 1 3 ϕ 2 } χ.
PHYSICS 34 QUANTUM PHYSICS II (25) Assigmet 2 Solutios 1. With respect to a pair of orthoormal vectors ϕ 1 ad ϕ 2 that spa the Hilbert space H of a certai system, the operator  is defied by its actio
More informationMath 525: Lecture 5. January 18, 2018
Math 525: Lecture 5 Jauary 18, 2018 1 Series (review) Defiitio 1.1. A sequece (a ) R coverges to a poit L R (writte a L or lim a = L) if for each ǫ > 0, we ca fid N such that a L < ǫ for all N. If the
More informationThe Growth of Functions. Theoretical Supplement
The Growth of Fuctios Theoretical Supplemet The Triagle Iequality The triagle iequality is a algebraic tool that is ofte useful i maipulatig absolute values of fuctios. The triagle iequality says that
More informationThe Random Walk For Dummies
The Radom Walk For Dummies Richard A Mote Abstract We look at the priciples goverig the oe-dimesioal discrete radom walk First we review five basic cocepts of probability theory The we cosider the Beroulli
More informationACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER / Statistics
ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER 1 018/019 DR. ANTHONY BROWN 8. Statistics 8.1. Measures of Cetre: Mea, Media ad Mode. If we have a series of umbers the
More informationNUMERICAL METHODS COURSEWORK INFORMAL NOTES ON NUMERICAL INTEGRATION COURSEWORK
NUMERICAL METHODS COURSEWORK INFORMAL NOTES ON NUMERICAL INTEGRATION COURSEWORK For this piece of coursework studets must use the methods for umerical itegratio they meet i the Numerical Methods module
More informationChapter 6 Overview: Sequences and Numerical Series. For the purposes of AP, this topic is broken into four basic subtopics:
Chapter 6 Overview: Sequeces ad Numerical Series I most texts, the topic of sequeces ad series appears, at first, to be a side topic. There are almost o derivatives or itegrals (which is what most studets
More informationMAT1026 Calculus II Basic Convergence Tests for Series
MAT026 Calculus II Basic Covergece Tests for Series Egi MERMUT 202.03.08 Dokuz Eylül Uiversity Faculty of Sciece Departmet of Mathematics İzmir/TURKEY Cotets Mootoe Covergece Theorem 2 2 Series of Real
More informationPHYS 321 Solutions to Practice Final (December 2002).
PHYS Solutios to Practice Fial (December ) Two masses, m ad m are coected by a sprig of costat k, leadig to the potetial V( r) = k( r ) r a) What is the Lagragia for this system? (Assume -dimesioal motio)
More information17 Phonons and conduction electrons in solids (Hiroshi Matsuoka)
7 Phoos ad coductio electros i solids Hiroshi Matsuoa I this chapter we will discuss a miimal microscopic model for phoos i a solid ad a miimal microscopic model for coductio electros i a simple metal.
More informationMATH 10550, EXAM 3 SOLUTIONS
MATH 155, EXAM 3 SOLUTIONS 1. I fidig a approximate solutio to the equatio x 3 +x 4 = usig Newto s method with iitial approximatio x 1 = 1, what is x? Solutio. Recall that x +1 = x f(x ) f (x ). Hece,
More information18.01 Calculus Jason Starr Fall 2005
Lecture 18. October 5, 005 Homework. Problem Set 5 Part I: (c). Practice Problems. Course Reader: 3G 1, 3G, 3G 4, 3G 5. 1. Approximatig Riema itegrals. Ofte, there is o simpler expressio for the atiderivative
More informationMATH 21 SECTION NOTES
MATH SECTION NOTES EVAN WARNER. March 9.. Admiistrative miscellay. These weekly sectios will be for some review ad may example problems, i geeral. Attedace will be take as per class policy. We will be
More informationR is a scalar defined as follows:
Math 8. Notes o Dot Product, Cross Product, Plaes, Area, ad Volumes This lecture focuses primarily o the dot product ad its may applicatios, especially i the measuremet of agles ad scalar projectio ad
More informationLecture 6: Integration and the Mean Value Theorem
Math 8 Istructor: Padraic Bartlett Lecture 6: Itegratio ad the Mea Value Theorem Week 6 Caltech - Fall, 2011 1 Radom Questios Questio 1.1. Show that ay positive ratioal umber ca be writte as the sum of
More information(3) If you replace row i of A by its sum with a multiple of another row, then the determinant is unchanged! Expand across the i th row:
Math 5-4 Tue Feb 4 Cotiue with sectio 36 Determiats The effective way to compute determiats for larger-sized matrices without lots of zeroes is to ot use the defiitio, but rather to use the followig facts,
More informationLinear Differential Equations of Higher Order Basic Theory: Initial-Value Problems d y d y dy
Liear Differetial Equatios of Higher Order Basic Theory: Iitial-Value Problems d y d y dy Solve: a( ) + a ( )... a ( ) a0( ) y g( ) + + + = d d d ( ) Subject to: y( 0) = y0, y ( 0) = y,..., y ( 0) = y
More informationThe Binomial Theorem
The Biomial Theorem Robert Marti Itroductio The Biomial Theorem is used to expad biomials, that is, brackets cosistig of two distict terms The formula for the Biomial Theorem is as follows: (a + b ( k
More informationOptimization Methods MIT 2.098/6.255/ Final exam
Optimizatio Methods MIT 2.098/6.255/15.093 Fial exam Date Give: December 19th, 2006 P1. [30 pts] Classify the followig statemets as true or false. All aswers must be well-justified, either through a short
More information10-701/ Machine Learning Mid-term Exam Solution
0-70/5-78 Machie Learig Mid-term Exam Solutio Your Name: Your Adrew ID: True or False (Give oe setece explaatio) (20%). (F) For a cotiuous radom variable x ad its probability distributio fuctio p(x), it
More informationMost text will write ordinary derivatives using either Leibniz notation 2 3. y + 5y= e and y y. xx tt t
Itroductio to Differetial Equatios Defiitios ad Termiolog Differetial Equatio: A equatio cotaiig the derivatives of oe or more depedet variables, with respect to oe or more idepedet variables, is said
More informationWHAT IS THE PROBABILITY FUNCTION FOR LARGE TSUNAMI WAVES? ABSTRACT
WHAT IS THE PROBABILITY FUNCTION FOR LARGE TSUNAMI WAVES? Harold G. Loomis Hoolulu, HI ABSTRACT Most coastal locatios have few if ay records of tsuami wave heights obtaied over various time periods. Still
More informationRoberto s Notes on Series Chapter 2: Convergence tests Section 7. Alternating series
Roberto s Notes o Series Chapter 2: Covergece tests Sectio 7 Alteratig series What you eed to kow already: All basic covergece tests for evetually positive series. What you ca lear here: A test for series
More informationMAT 271 Project: Partial Fractions for certain rational functions
MAT 7 Project: Partial Fractios for certai ratioal fuctios Prerequisite kowledge: partial fractios from MAT 7, a very good commad of factorig ad complex umbers from Precalculus. To complete this project,
More informationCHAPTER 8 SYSTEMS OF PARTICLES
CHAPTER 8 SYSTES OF PARTICLES CHAPTER 8 COLLISIONS 45 8. CENTER OF ASS The ceter of mass of a system of particles or a rigid body is the poit at which all of the mass are cosidered to be cocetrated there
More informationInfinite Series. Definition. An infinite series is an expression of the form. Where the numbers u k are called the terms of the series.
Ifiite Series Defiitio. A ifiite series is a expressio of the form uk = u + u + u + + u + () 2 3 k Where the umbers u k are called the terms of the series. Such a expressio is meat to be the result of
More informationf(x) dx as we do. 2x dx x also diverges. Solution: We compute 2x dx lim
Math 3, Sectio 2. (25 poits) Why we defie f(x) dx as we do. (a) Show that the improper itegral diverges. Hece the improper itegral x 2 + x 2 + b also diverges. Solutio: We compute x 2 + = lim b x 2 + =
More information9.3 The INTEGRAL TEST; p-series
Lecture 9.3 & 9.4 Math 0B Nguye of 6 Istructor s Versio 9.3 The INTEGRAL TEST; p-series I this ad the followig sectio, you will study several covergece tests that apply to series with positive terms. Note
More information7 Sequences of real numbers
40 7 Sequeces of real umbers 7. Defiitios ad examples Defiitio 7... A sequece of real umbers is a real fuctio whose domai is the set N of atural umbers. Let s : N R be a sequece. The the values of s are
More informationChapter 7: Numerical Series
Chapter 7: Numerical Series Chapter 7 Overview: Sequeces ad Numerical Series I most texts, the topic of sequeces ad series appears, at first, to be a side topic. There are almost o derivatives or itegrals
More information