THE MYSTERY OF STOCHASTIC MECHANICS. Edward Nelson Department of Mathematics Princeton University

Size: px
Start display at page:

Download "THE MYSTERY OF STOCHASTIC MECHANICS. Edward Nelson Department of Mathematics Princeton University"

Transcription

1 THE MYSTERY OF STOCHASTIC MECHANICS Edward Nelson Deparmen of Mahemaics Princeon Universiy 1

2 Classical Hamilon-Jacobi heory N paricles of various masses on a Euclidean space. Incorporae he masses in he fla Riemannian meric m ij, he mass ensor. Then if v i is a velociy, v i = m ij v j is a momenum. Kineic energy: 1 2 vi v i. Poenial energy: V. Lagrangian: L = 1 2 vi v i V. Posiion a ime of he configuraion: ξ(). Iniial ime:. Final ime: 1. Hamilon s principal funcion: S(x, ) = 1 L ( ξ(s) ) ds. Subsanial derivaive (derivaive along rajecories): Then DS = L. D = + vi i. 2

3 Vecor field: v wih dξ d = v. Principle of leas acion in Hamilon-Jacobi heory: v is a criical poin for S, for unconsrained variaions. Tha is, le v be anoher vecor field, le δv = v v, and denoe by a prime quaniies wih v replaced by v. Then D(S S) = D S DS + (D D )S = L L δv i i S = L L δv i i S + o(δv). Now L L = v i δv i + o(δv), S S = 1 (v i i S)δv i ds + o(δv). Since his is rue for all variaions, we have he Hamilon-Jacobi condiion: v i = i S. 3

4 Togeher wih DS = L; i.e., s + vi i S = 1 2 vi v i V, his gives he Hamilon-Jacobi equaion S i i S + V = 0. If we ake he gradien we obain Newon s equaion F = ma. 4

5 Sochasic Hamilon-Jacobi heory Following Guerra and Morao, consruc a conservaive Markovian dynamics for a Markov process ξ. Wiener process (Brownian moion) w wih Edw i ()dw i () = hd. Kinemaics: dξ() = b ( ξ(), ) d + dw(). Dynamics: δe L d = 0 (heurisically). The rajecories are non-differeniable, so wha is he meaning of 1 2 dξ i d dξ i d in he Lagrangian L? 5

6 Le d > 0 and dξ() = ξ( + d) ξ(), so dξ d is a quoien, no a derivaive. Compue E 1 dξ i dξ i up o o(1). 2 d d dξ i = +d Noe: dw is of order d 1/2. b i( ξ(r), r ) dr + dw i. dξ i = where +d + dw i b i( ξ() + r = b i d + k b i W k + dw i + O(d 2 ) b ( ) ξ(s), s)ds + w(r) w(), r dr W k = +d [w k (r) w k ()] dr. Therefore 6

7 1 2 dξi dξ i = 1 2 bi b i d + b i dw i d + k b i W k dw i dwi dw i + o(d 2 ). Now EW k dw i = hδ k i +d (r )d = h 2 δk i, so k b i W k = h 2 ib i. The erm b i dw i d is singular, of order d 3/2, bu is expecaion is 0. Finally, Edw i dw i = h 2 nd. Hence we have he sough-for resul: E 1 2 dξ i d dξ i d = 1 2 bi b i + h 2 ib i + h 2 n d + o(1). 7

8 The singular erm h n is a consan, no depending on he rajecory, and i drops ou of he 2 d variaion. Le and L + = 1 2 bi b i + h 2 ib i V I = E 1 L + ( ξ(s) ) ds. Le δb be a vecor field, b = b + δb, and as before denoe by primes quaniies wih b replaced by b. Le ξ saisfy dξ () = b ( ξ (), ) d + dw() wih ξ () = ξ() for he iniial ime. Definiion: ξ is criical for L in case I I = o(δb). 8

9 Sochasic Hamilon s principal funcion: Then S(x, ) = E x, 1 DS = L + L + ( ξ(s), s ) ds. where now D is he mean forward derivaive: ( ) df() Df = lim E. d 0+ d and As before, D(S S) = D S DS + (D D )S Now = L + L + δb i i S = L + L + δb i i S + o(δb). L + L + = b i δb i + h 2 iδb i + o(δb) ( ) I I = E 1 ( b i i S + h ) 2 i δb i ds + o(δb). 9

10 Digression on Markovian kinemaics: Markov process: given he presen, pas and fuure are independen. The wo direcions of ime are on an equal fooing. In addiion o he forward drif b here is he backward drif b. v = b + b 2 u = b b 2 curren velociy, osmoic velociy. Le ρ be he probabiliy densiy of ξ. Then v = (vρ) u = h ρ 2 ρ curren equaion, osmoic equaion. End of digression. 10

11 The equaion (*), namely I I = E 1 ( b i i S + h ) 2 i δb i ds + o(δb), is awkward because i involves i δb i. Inegrae by pars: E h 2 iδb i( ξ(s), s ) h ( = i δb i) ρ dx R n 2 = R δbi u i ρ dx. n Since b u = v, I I = E 1 (v i i S)δb i ds + o(δb). This is rue for all variaions, so ξ is criical for L if and only if he sochasic Hamilon-Jacobi condiion holds: v i = i S. 11

12 Le R = h 2 log ρ, so i R = u i. If we wrie ou DS = L + and express everyhing in erms of R and S, we obain he sochasic Hamilon-Jacobi equaion: (1) S i S i S + V 1 2 i R i R h R = 0. 2 Expressing he curren equaion v = (vρ) in erms of R and S we obain: (2) R + ir i S + h S = 0. 2 These wo equaions are a sysem of coupled nonlinear parial differenial equaions expressing necessary and sufficien condiions for a Markov process o be criical. How can we solve hem? 12

13 Le ψ = e 1 h (R+iS). Then he sysem is equivalen o he Schrödinger equaion ψ = ī h ] [ h2 + V ψ. 13

14 magneic fields Riemannian manifolds spin Bose-Einsein and Fermi-Dirac saisics exisence of finie-energy Markov processes momenum inerference 14

15 Non-localiy Consider wo paricles on R 1 wih iniial wave funcion ψ(0) = 1 2π e 1 4 σ 1 (0)x x where σ 1 (0) = ( 1 ) and x = ( x 1 x 2 ). The x 1 and x 2 axes are unrelaed; he paricles may be separaed by an arbirarily large amoun a. A ime 0, urn on a linear resoring force (harmonic oscillaor wih circular frequency ω) for he second paricle. Then he paricles are dynamically uncoupled. 15

16 Since he paricles are very widely separaed and dynamically uncoupled, we should expec ha does no depend on ω. Eξ 1 ()ξ 1 (0) In fac i does no o fourh order in, bu neverheless he rajecory of he firs paricle is immediaely affeced by he choice of ω in a far disan place. For me his is unphysical, especially since he effec does no depend on he separaion a. 16

17 A wrong predicion Consider wo harmonic oscillaors, abou wo widely separaed poins a 1 and a 2, wih circular frequency 1, and le X i be he Heisenberg posiion operaor, in unis of disance from a i, wih Heisenberg momenum operaor P i. Then for i = 1, 2. X i () = cos X i (0) + sin P i (0) P i () = sin X i (0) + cos P i (0) Le he Heisenberg sae vecor ψ 0 be a real Gaussian cenered a (0, 0), and wrie A = (ψ 0, Aψ 0 ). Then X i () = P i () = 0 since his is rue for = 0. The operaors X 1 () and X 2 (s) commue. Choose ψ 0 so ha he correlaion X 1 (0)X 2 (0) is 90%. Thus he oscillaors are enangled bu dynamically uncoupled. The quanum mechanical correlaion funcion X 1 ()X 2 (0) is periodic of period 2π since X 1 () is. Hence X 1 ()X 2 (0) =.9 whenever is a muliple of 2π. 17

18 Bu le (ξ 1, ξ 2 ) be he corresponding Markov process of sochasic mechanics. This is a diffusion process and i evenually loses all memory of where i sared. Thus whereas lim E( ξ 1 (2πn)ξ 2 (0) ) = 0 n X 1 (2πn)X 2 (0) =.9. Here we have an empirical difference beween he predicions of quanum mechanics and sochasic mechanics. Measuremens of he posiion of he firs paricle a ime and of he second paricle a ime 0 do no inerfere wih each oher, and he wo heories predic oally differen saisics. Does anyone doub ha quanum mechanics is righ and sochasic mechanics is wrong? 18

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

An Introduction to Malliavin calculus and its applications

An Introduction to Malliavin calculus and its applications An Inroducion o Malliavin calculus and is applicaions Lecure 5: Smoohness of he densiy and Hörmander s heorem David Nualar Deparmen of Mahemaics Kansas Universiy Universiy of Wyoming Summer School 214

More information

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t M ah 5 2 7 Fall 2 0 0 9 L ecure 1 0 O c. 7, 2 0 0 9 Hamilon- J acobi Equaion: Explici Formulas In his lecure we ry o apply he mehod of characerisics o he Hamilon-Jacobi equaion: u + H D u, x = 0 in R n

More information

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and

More information

The expectation value of the field operator.

The expectation value of the field operator. The expecaion value of he field operaor. Dan Solomon Universiy of Illinois Chicago, IL dsolom@uic.edu June, 04 Absrac. Much of he mahemaical developmen of quanum field heory has been in suppor of deermining

More information

Lecture 4 Kinetics of a particle Part 3: Impulse and Momentum

Lecture 4 Kinetics of a particle Part 3: Impulse and Momentum MEE Engineering Mechanics II Lecure 4 Lecure 4 Kineics of a paricle Par 3: Impulse and Momenum Linear impulse and momenum Saring from he equaion of moion for a paricle of mass m which is subjeced o an

More information

Some Basic Information about M-S-D Systems

Some Basic Information about M-S-D Systems Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,

More information

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems. Mah 2250-004 Week 4 April 6-20 secions 7.-7.3 firs order sysems of linear differenial equaions; 7.4 mass-spring sysems. Mon Apr 6 7.-7.2 Sysems of differenial equaions (7.), and he vecor Calculus we need

More information

Lecture #11: Wavepacket Dynamics for Harmonic Oscillator

Lecture #11: Wavepacket Dynamics for Harmonic Oscillator Lecure #11: Wavepacke Dynamics for Harmonic Oscillaor and PIB Las ime: Time Dependen Schrödinger Equaion Ψ HHΨ = iħ Express Ψ in complee basis se of eigenfuncions of ime independen H H {ψ n (x), E n }

More information

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution Physics 7b: Saisical Mechanics Fokker-Planck Equaion The Langevin equaion approach o he evoluion of he velociy disribuion for he Brownian paricle migh leave you uncomforable. A more formal reamen of his

More information

Predator - Prey Model Trajectories and the nonlinear conservation law

Predator - Prey Model Trajectories and the nonlinear conservation law Predaor - Prey Model Trajecories and he nonlinear conservaion law James K. Peerson Deparmen of Biological Sciences and Deparmen of Mahemaical Sciences Clemson Universiy Ocober 28, 213 Ouline Drawing Trajecories

More information

2.3 SCHRÖDINGER AND HEISENBERG REPRESENTATIONS

2.3 SCHRÖDINGER AND HEISENBERG REPRESENTATIONS Andrei Tokmakoff, MIT Deparmen of Chemisry, 2/22/2007 2-17 2.3 SCHRÖDINGER AND HEISENBERG REPRESENTATIONS The mahemaical formulaion of he dynamics of a quanum sysem is no unique. So far we have described

More information

Problem Set on Differential Equations

Problem Set on Differential Equations Problem Se on Differenial Equaions 1. Solve he following differenial equaions (a) x () = e x (), x () = 3/ 4. (b) x () = e x (), x (1) =. (c) xe () = + (1 x ()) e, x () =.. (An asse marke model). Le p()

More information

ψ(t) = V x (0)V x (t)

ψ(t) = V x (0)V x (t) .93 Home Work Se No. (Professor Sow-Hsin Chen Spring Term 5. Due March 7, 5. This problem concerns calculaions of analyical expressions for he self-inermediae scaering funcion (ISF of he es paricle in

More information

Homogenization of random Hamilton Jacobi Bellman Equations

Homogenization of random Hamilton Jacobi Bellman Equations Probabiliy, Geomery and Inegrable Sysems MSRI Publicaions Volume 55, 28 Homogenizaion of random Hamilon Jacobi Bellman Equaions S. R. SRINIVASA VARADHAN ABSTRACT. We consider nonlinear parabolic equaions

More information

14 Autoregressive Moving Average Models

14 Autoregressive Moving Average Models 14 Auoregressive Moving Average Models In his chaper an imporan parameric family of saionary ime series is inroduced, he family of he auoregressive moving average, or ARMA, processes. For a large class

More information

Module 3: The Damped Oscillator-II Lecture 3: The Damped Oscillator-II

Module 3: The Damped Oscillator-II Lecture 3: The Damped Oscillator-II Module 3: The Damped Oscillaor-II Lecure 3: The Damped Oscillaor-II 3. Over-damped Oscillaions. This refers o he siuaion where β > ω (3.) The wo roos are and α = β + α 2 = β β 2 ω 2 = (3.2) β 2 ω 2 = 2

More information

Two Coupled Oscillators / Normal Modes

Two Coupled Oscillators / Normal Modes Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own

More information

IB Physics Kinematics Worksheet

IB Physics Kinematics Worksheet IB Physics Kinemaics Workshee Wrie full soluions and noes for muliple choice answers. Do no use a calculaor for muliple choice answers. 1. Which of he following is a correc definiion of average acceleraion?

More information

Major Concepts. Brownian Motion & More. Chemical Kinetics Master equation & Detailed Balance Relaxation rate & Inverse Phenomenological Rate

Major Concepts. Brownian Motion & More. Chemical Kinetics Master equation & Detailed Balance Relaxation rate & Inverse Phenomenological Rate Major Conceps Brownian Moion & More Langevin Equaion Model for a agged subsysem in a solven Harmonic bah wih emperaure, T Fricion & Correlaed forces (FDR) Markovian/Ohmic vs. Memory Fokker-Planck Equaion

More information

Fractional Method of Characteristics for Fractional Partial Differential Equations

Fractional Method of Characteristics for Fractional Partial Differential Equations Fracional Mehod of Characerisics for Fracional Parial Differenial Equaions Guo-cheng Wu* Modern Teile Insiue, Donghua Universiy, 188 Yan-an ilu Road, Shanghai 51, PR China Absrac The mehod of characerisics

More information

Applications of the Basic Equations Chapter 3. Paul A. Ullrich

Applications of the Basic Equations Chapter 3. Paul A. Ullrich Applicaions of he Basic Equaions Chaper 3 Paul A. Ullrich paullrich@ucdavis.edu Par 1: Naural Coordinaes Naural Coordinaes Quesion: Why do we need anoher coordinae sysem? Our goal is o simplify he equaions

More information

Lecture 10: The Poincaré Inequality in Euclidean space

Lecture 10: The Poincaré Inequality in Euclidean space Deparmens of Mahemaics Monana Sae Universiy Fall 215 Prof. Kevin Wildrick n inroducion o non-smooh analysis and geomery Lecure 1: The Poincaré Inequaliy in Euclidean space 1. Wha is he Poincaré inequaliy?

More information

Lecture 9: Advanced DFT concepts: The Exchange-correlation functional and time-dependent DFT

Lecture 9: Advanced DFT concepts: The Exchange-correlation functional and time-dependent DFT Lecure 9: Advanced DFT conceps: The Exchange-correlaion funcional and ime-dependen DFT Marie Curie Tuorial Series: Modeling Biomolecules December 6-11, 2004 Mark Tuckerman Dep. of Chemisry and Couran Insiue

More information

The Quantum Theory of Atoms and Molecules: The Schrodinger equation. Hilary Term 2008 Dr Grant Ritchie

The Quantum Theory of Atoms and Molecules: The Schrodinger equation. Hilary Term 2008 Dr Grant Ritchie e Quanum eory of Aoms and Molecules: e Scrodinger equaion Hilary erm 008 Dr Gran Ricie An equaion for maer waves? De Broglie posulaed a every paricles as an associaed wave of waveleng: / p Wave naure of

More information

From Particles to Rigid Bodies

From Particles to Rigid Bodies Rigid Body Dynamics From Paricles o Rigid Bodies Paricles No roaions Linear velociy v only Rigid bodies Body roaions Linear velociy v Angular velociy ω Rigid Bodies Rigid bodies have boh a posiion and

More information

Homework Solution Set # 3. Thursday, September 22, Textbook: Claude Cohen Tannoudji, Bernard Diu and Franck Lalo, Second Volume Complement G X

Homework Solution Set # 3. Thursday, September 22, Textbook: Claude Cohen Tannoudji, Bernard Diu and Franck Lalo, Second Volume Complement G X Deparmen of Physics Quanum Mechanics II, 570 Temple Universiy Insrucor: Z.-E. Meziani Homework Soluion Se # 3 Thursday, Sepember, 06 Texbook: Claude Cohen Tannoudji, Bernard Diu and Franck Lalo, Second

More information

Lecture 20: Riccati Equations and Least Squares Feedback Control

Lecture 20: Riccati Equations and Least Squares Feedback Control 34-5 LINEAR SYSTEMS Lecure : Riccai Equaions and Leas Squares Feedback Conrol 5.6.4 Sae Feedback via Riccai Equaions A recursive approach in generaing he marix-valued funcion W ( ) equaion for i for he

More information

Class Meeting # 10: Introduction to the Wave Equation

Class Meeting # 10: Introduction to the Wave Equation MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion

More information

Hamilton Jacobi equations

Hamilton Jacobi equations Hamilon Jacobi equaions Inoducion o PDE The rigorous suff from Evans, mosly. We discuss firs u + H( u = 0, (1 where H(p is convex, and superlinear a infiniy, H(p lim p p = + This by comes by inegraion

More information

where the coordinate X (t) describes the system motion. X has its origin at the system static equilibrium position (SEP).

where the coordinate X (t) describes the system motion. X has its origin at the system static equilibrium position (SEP). Appendix A: Conservaion of Mechanical Energy = Conservaion of Linear Momenum Consider he moion of a nd order mechanical sysem comprised of he fundamenal mechanical elemens: ineria or mass (M), siffness

More information

Announcements: Warm-up Exercise:

Announcements: Warm-up Exercise: Fri Apr 13 7.1 Sysems of differenial equaions - o model muli-componen sysems via comparmenal analysis hp//en.wikipedia.org/wiki/muli-comparmen_model Announcemens Warm-up Exercise Here's a relaively simple

More information

What Ties Return Volatilities to Price Valuations and Fundamentals? On-Line Appendix

What Ties Return Volatilities to Price Valuations and Fundamentals? On-Line Appendix Wha Ties Reurn Volailiies o Price Valuaions and Fundamenals? On-Line Appendix Alexander David Haskayne School of Business, Universiy of Calgary Piero Veronesi Universiy of Chicago Booh School of Business,

More information

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation: M ah 5 7 Fall 9 L ecure O c. 4, 9 ) Hamilon- J acobi Equaion: Weak S oluion We coninue he sudy of he Hamilon-Jacobi equaion: We have shown ha u + H D u) = R n, ) ; u = g R n { = }. ). In general we canno

More information

Optimizing heat exchangers

Optimizing heat exchangers Opimizing hea echangers Jean-Luc Thiffeaul Deparmen of Mahemaics, Universiy of Wisconsin Madison, 48 Lincoln Dr., Madison, WI 5376, USA wih: Florence Marcoe, Charles R. Doering, William R. Young (Daed:

More information

A. Using Newton s second law in one dimension, F net. , write down the differential equation that governs the motion of the block.

A. Using Newton s second law in one dimension, F net. , write down the differential equation that governs the motion of the block. Simple SIMPLE harmonic HARMONIC moion MOTION I. Differenial equaion of moion A block is conneced o a spring, one end of which is aached o a wall. (Neglec he mass of he spring, and assume he surface is

More information

6. Stochastic calculus with jump processes

6. Stochastic calculus with jump processes A) Trading sraegies (1/3) Marke wih d asses S = (S 1,, S d ) A rading sraegy can be modelled wih a vecor φ describing he quaniies invesed in each asse a each insan : φ = (φ 1,, φ d ) The value a of a porfolio

More information

2001 November 15 Exam III Physics 191

2001 November 15 Exam III Physics 191 1 November 15 Eam III Physics 191 Physical Consans: Earh s free-fall acceleraion = g = 9.8 m/s 2 Circle he leer of he single bes answer. quesion is worh 1 poin Each 3. Four differen objecs wih masses:

More information

KINEMATICS IN ONE DIMENSION

KINEMATICS IN ONE DIMENSION KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec

More information

Math 334 Fall 2011 Homework 11 Solutions

Math 334 Fall 2011 Homework 11 Solutions Dec. 2, 2 Mah 334 Fall 2 Homework Soluions Basic Problem. Transform he following iniial value problem ino an iniial value problem for a sysem: u + p()u + q() u g(), u() u, u () v. () Soluion. Le v u. Then

More information

An Introduction to Backward Stochastic Differential Equations (BSDEs) PIMS Summer School 2016 in Mathematical Finance.

An Introduction to Backward Stochastic Differential Equations (BSDEs) PIMS Summer School 2016 in Mathematical Finance. 1 An Inroducion o Backward Sochasic Differenial Equaions (BSDEs) PIMS Summer School 2016 in Mahemaical Finance June 25, 2016 Chrisoph Frei cfrei@ualbera.ca This inroducion is based on Touzi [14], Bouchard

More information

Backward stochastic dynamics on a filtered probability space

Backward stochastic dynamics on a filtered probability space Backward sochasic dynamics on a filered probabiliy space Gechun Liang Oxford-Man Insiue, Universiy of Oxford based on join work wih Terry Lyons and Zhongmin Qian Page 1 of 15 gliang@oxford-man.ox.ac.uk

More information

A quantum method to test the existence of consciousness

A quantum method to test the existence of consciousness A quanum mehod o es he exisence of consciousness Gao Shan The Scieniss Work Team of Elecro-Magneic Wave Velociy, Chinese Insiue of Elecronics -0, NO.0 Building, YueTan XiJie DongLi, XiCheng Disric Beijing

More information

k 1 k 2 x (1) x 2 = k 1 x 1 = k 2 k 1 +k 2 x (2) x k series x (3) k 2 x 2 = k 1 k 2 = k 1+k 2 = 1 k k 2 k series

k 1 k 2 x (1) x 2 = k 1 x 1 = k 2 k 1 +k 2 x (2) x k series x (3) k 2 x 2 = k 1 k 2 = k 1+k 2 = 1 k k 2 k series Final Review A Puzzle... Consider wo massless springs wih spring consans k 1 and k and he same equilibrium lengh. 1. If hese springs ac on a mass m in parallel, hey would be equivalen o a single spring

More information

Math 333 Problem Set #2 Solution 14 February 2003

Math 333 Problem Set #2 Solution 14 February 2003 Mah 333 Problem Se #2 Soluion 14 February 2003 A1. Solve he iniial value problem dy dx = x2 + e 3x ; 2y 4 y(0) = 1. Soluion: This is separable; we wrie 2y 4 dy = x 2 + e x dx and inegrae o ge The iniial

More information

System of Linear Differential Equations

System of Linear Differential Equations Sysem of Linear Differenial Equaions In "Ordinary Differenial Equaions" we've learned how o solve a differenial equaion for a variable, such as: y'k5$e K2$x =0 solve DE yx = K 5 2 ek2 x C_C1 2$y''C7$y

More information

Chapter 12: Velocity, acceleration, and forces

Chapter 12: Velocity, acceleration, and forces To Feel a Force Chaper Spring, Chaper : A. Saes of moion For moion on or near he surface of he earh, i is naural o measure moion wih respec o objecs fixed o he earh. The 4 hr. roaion of he earh has a measurable

More information

Chapter 3 Boundary Value Problem

Chapter 3 Boundary Value Problem Chaper 3 Boundary Value Problem A boundary value problem (BVP) is a problem, ypically an ODE or a PDE, which has values assigned on he physical boundary of he domain in which he problem is specified. Le

More information

Summary of shear rate kinematics (part 1)

Summary of shear rate kinematics (part 1) InroToMaFuncions.pdf 4 CM465 To proceed o beer-designed consiuive equaions, we need o know more abou maerial behavior, i.e. we need more maerial funcions o predic, and we need measuremens of hese maerial

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

Cash Flow Valuation Mode Lin Discrete Time

Cash Flow Valuation Mode Lin Discrete Time IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728,p-ISSN: 2319-765X, 6, Issue 6 (May. - Jun. 2013), PP 35-41 Cash Flow Valuaion Mode Lin Discree Time Olayiwola. M. A. and Oni, N. O. Deparmen of Mahemaics

More information

( ) = b n ( t) n " (2.111) or a system with many states to be considered, solving these equations isn t. = k U I ( t,t 0 )! ( t 0 ) (2.

( ) = b n ( t) n  (2.111) or a system with many states to be considered, solving these equations isn t. = k U I ( t,t 0 )! ( t 0 ) (2. Andrei Tokmakoff, MIT Deparmen of Chemisry, 3/14/007-6.4 PERTURBATION THEORY Given a Hamilonian H = H 0 + V where we know he eigenkes for H 0 : H 0 n = E n n, we can calculae he evoluion of he wavefuncion

More information

Second quantization and gauge invariance.

Second quantization and gauge invariance. 1 Second quanizaion and gauge invariance. Dan Solomon Rauland-Borg Corporaion Moun Prospec, IL Email: dsolom@uic.edu June, 1. Absrac. I is well known ha he single paricle Dirac equaion is gauge invarian.

More information

Solutions of Sample Problems for Third In-Class Exam Math 246, Spring 2011, Professor David Levermore

Solutions of Sample Problems for Third In-Class Exam Math 246, Spring 2011, Professor David Levermore Soluions of Sample Problems for Third In-Class Exam Mah 6, Spring, Professor David Levermore Compue he Laplace ransform of f e from is definiion Soluion The definiion of he Laplace ransform gives L[f]s

More information

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should Cambridge Universiy Press 978--36-60033-7 Cambridge Inernaional AS and A Level Mahemaics: Mechanics Coursebook Excerp More Informaion Chaper The moion of projeciles In his chaper he model of free moion

More information

On Option Pricing by Quantum Mechanics Approach

On Option Pricing by Quantum Mechanics Approach On Opion Pricing by Quanum Mechanics Approach Hiroshi Inoue School of Managemen okyo Universiy of Science Kuki, Saiama 46-851 Japan e-mail:inoue@mskukiusacp Absrac We discuss he pah inegral mehodology

More information

Undetermined coefficients for local fractional differential equations

Undetermined coefficients for local fractional differential equations Available online a www.isr-publicaions.com/jmcs J. Mah. Compuer Sci. 16 (2016), 140 146 Research Aricle Undeermined coefficiens for local fracional differenial equaions Roshdi Khalil a,, Mohammed Al Horani

More information

Optimal Path Planning for Flexible Redundant Robot Manipulators

Optimal Path Planning for Flexible Redundant Robot Manipulators 25 WSEAS In. Conf. on DYNAMICAL SYSEMS and CONROL, Venice, Ialy, November 2-4, 25 (pp363-368) Opimal Pah Planning for Flexible Redundan Robo Manipulaors H. HOMAEI, M. KESHMIRI Deparmen of Mechanical Engineering

More information

PROBLEMS FOR MATH 162 If a problem is starred, all subproblems are due. If only subproblems are starred, only those are due. SLOPES OF TANGENT LINES

PROBLEMS FOR MATH 162 If a problem is starred, all subproblems are due. If only subproblems are starred, only those are due. SLOPES OF TANGENT LINES PROBLEMS FOR MATH 6 If a problem is sarred, all subproblems are due. If onl subproblems are sarred, onl hose are due. 00. Shor answer quesions. SLOPES OF TANGENT LINES (a) A ball is hrown ino he air. Is

More information

Vanishing Viscosity Method. There are another instructive and perhaps more natural discontinuous solutions of the conservation law

Vanishing Viscosity Method. There are another instructive and perhaps more natural discontinuous solutions of the conservation law Vanishing Viscosiy Mehod. There are anoher insrucive and perhaps more naural disconinuous soluions of he conservaion law (1 u +(q(u x 0, he so called vanishing viscosiy mehod. This mehod consiss in viewing

More information

2.1: What is physics? Ch02: Motion along a straight line. 2.2: Motion. 2.3: Position, Displacement, Distance

2.1: What is physics? Ch02: Motion along a straight line. 2.2: Motion. 2.3: Position, Displacement, Distance Ch: Moion along a sraigh line Moion Posiion and Displacemen Average Velociy and Average Speed Insananeous Velociy and Speed Acceleraion Consan Acceleraion: A Special Case Anoher Look a Consan Acceleraion

More information

Math 527 Lecture 6: Hamilton-Jacobi Equation: Explicit Formulas

Math 527 Lecture 6: Hamilton-Jacobi Equation: Explicit Formulas Mah 527 Lecure 6: Hamilon-Jacobi Equaion: Explici Formulas Sep. 23, 2 Mehod of characerisics. We r o appl he mehod of characerisics o he Hamilon-Jacobi equaion: u +Hx, Du = in R n, u = g on R n =. 2 To

More information

Random Walk on Circle Imagine a Markov process governing the random motion of a particle on a circular

Random Walk on Circle Imagine a Markov process governing the random motion of a particle on a circular Random Walk on Circle Imagine a Markov process governing he random moion of a paricle on a circular laice: 1 2 γ γ γ The paricle moves o he righ or lef wih probabiliy γ and says where i is wih probabiliy

More information

arxiv: v1 [math.pr] 19 Feb 2011

arxiv: v1 [math.pr] 19 Feb 2011 A NOTE ON FELLER SEMIGROUPS AND RESOLVENTS VADIM KOSTRYKIN, JÜRGEN POTTHOFF, AND ROBERT SCHRADER ABSTRACT. Various equivalen condiions for a semigroup or a resolven generaed by a Markov process o be of

More information

Electrical Circuits. 1. Circuit Laws. Tools Used in Lab 13 Series Circuits Damped Vibrations: Energy Van der Pol Circuit

Electrical Circuits. 1. Circuit Laws. Tools Used in Lab 13 Series Circuits Damped Vibrations: Energy Van der Pol Circuit V() R L C 513 Elecrical Circuis Tools Used in Lab 13 Series Circuis Damped Vibraions: Energy Van der Pol Circui A series circui wih an inducor, resisor, and capacior can be represened by Lq + Rq + 1, a

More information

Week 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem)

Week 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem) Week 1 Lecure Problems, 5 Wha if somehing oscillaes wih no obvious spring? Wha is ω? (problem se problem) Sar wih Try and ge o SHM form E. Full beer can in lake, oscillaing F = m & = ge rearrange: F =

More information

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still. Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in

More information

Stochastic Model for Cancer Cell Growth through Single Forward Mutation

Stochastic Model for Cancer Cell Growth through Single Forward Mutation Journal of Modern Applied Saisical Mehods Volume 16 Issue 1 Aricle 31 5-1-2017 Sochasic Model for Cancer Cell Growh hrough Single Forward Muaion Jayabharahiraj Jayabalan Pondicherry Universiy, jayabharahi8@gmail.com

More information

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t... Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger

More information

Math 10B: Mock Mid II. April 13, 2016

Math 10B: Mock Mid II. April 13, 2016 Name: Soluions Mah 10B: Mock Mid II April 13, 016 1. ( poins) Sae, wih jusificaion, wheher he following saemens are rue or false. (a) If a 3 3 marix A saisfies A 3 A = 0, hen i canno be inverible. True.

More information

Traveling Waves. Chapter Introduction

Traveling Waves. Chapter Introduction Chaper 4 Traveling Waves 4.1 Inroducion To dae, we have considered oscillaions, i.e., periodic, ofen harmonic, variaions of a physical characerisic of a sysem. The sysem a one ime is indisinguishable from

More information

Estimation of Poses with Particle Filters

Estimation of Poses with Particle Filters Esimaion of Poses wih Paricle Filers Dr.-Ing. Bernd Ludwig Chair for Arificial Inelligence Deparmen of Compuer Science Friedrich-Alexander-Universiä Erlangen-Nürnberg 12/05/2008 Dr.-Ing. Bernd Ludwig (FAU

More information

We here collect a few numerical tests, in order to put into evidence the potentialities of HBVMs [4, 6, 7].

We here collect a few numerical tests, in order to put into evidence the potentialities of HBVMs [4, 6, 7]. Chaper Numerical Tess We here collec a few numerical ess, in order o pu ino evidence he poenialiies of HBVMs [4, 6, 7]. Tes problem Le us consider he problem characerized by he polynomial Hamilonian (4.)

More information

Introduction D P. r = constant discount rate, g = Gordon Model (1962): constant dividend growth rate.

Introduction D P. r = constant discount rate, g = Gordon Model (1962): constant dividend growth rate. Inroducion Gordon Model (1962): D P = r g r = consan discoun rae, g = consan dividend growh rae. If raional expecaions of fuure discoun raes and dividend growh vary over ime, so should he D/P raio. Since

More information

Utility maximization in incomplete markets

Utility maximization in incomplete markets Uiliy maximizaion in incomplee markes Marcel Ladkau 27.1.29 Conens 1 Inroducion and general seings 2 1.1 Marke model....................................... 2 1.2 Trading sraegy.....................................

More information

3, so θ = arccos

3, so θ = arccos Mahemaics 210 Professor Alan H Sein Monday, Ocober 1, 2007 SOLUTIONS This problem se is worh 50 poins 1 Find he angle beween he vecors (2, 7, 3) and (5, 2, 4) Soluion: Le θ be he angle (2, 7, 3) (5, 2,

More information

Today in Physics 218: radiation reaction

Today in Physics 218: radiation reaction Today in Physics 18: radiaion reacion Radiaion reacion The Abraham-Lorenz formula; radiaion reacion force The pah of he elecron in oday s firs example (radial decay grealy exaggeraed) 6 March 004 Physics

More information

4.5 Constant Acceleration

4.5 Constant Acceleration 4.5 Consan Acceleraion v() v() = v 0 + a a() a a() = a v 0 Area = a (a) (b) Figure 4.8 Consan acceleraion: (a) velociy, (b) acceleraion When he x -componen of he velociy is a linear funcion (Figure 4.8(a)),

More information

(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4)

(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4) Phase Plane Analysis of Linear Sysems Adaped from Applied Nonlinear Conrol by Sloine and Li The general form of a linear second-order sysem is a c b d From and b bc d a Differeniaion of and hen subsiuion

More information

EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE

EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Version April 30, 2004.Submied o CTU Repors. EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Per Krysl Universiy of California, San Diego La Jolla, California 92093-0085,

More information

2. Nonlinear Conservation Law Equations

2. Nonlinear Conservation Law Equations . Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear

More information

L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms

L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS NA568 Mobile Roboics: Mehods & Algorihms Today s Topic Quick review on (Linear) Kalman Filer Kalman Filering for Non-Linear Sysems Exended Kalman Filer (EKF)

More information

Ch1: Introduction and Review

Ch1: Introduction and Review //6 Ch: Inroducion and Review. Soli and flui; Coninuum hypohesis; Transpor phenomena (i) Solid vs. Fluid No exernal force : An elemen of solid has a preferred shape; fluid does no. Under he acion of a

More information

1. VELOCITY AND ACCELERATION

1. VELOCITY AND ACCELERATION 1. VELOCITY AND ACCELERATION 1.1 Kinemaics Equaions s = u + 1 a and s = v 1 a s = 1 (u + v) v = u + as 1. Displacemen-Time Graph Gradien = speed 1.3 Velociy-Time Graph Gradien = acceleraion Area under

More information

HOMEWORK # 2: MATH 211, SPRING Note: This is the last solution set where I will describe the MATLAB I used to make my pictures.

HOMEWORK # 2: MATH 211, SPRING Note: This is the last solution set where I will describe the MATLAB I used to make my pictures. HOMEWORK # 2: MATH 2, SPRING 25 TJ HITCHMAN Noe: This is he las soluion se where I will describe he MATLAB I used o make my picures.. Exercises from he ex.. Chaper 2.. Problem 6. We are o show ha y() =

More information

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x WEEK-3 Reciaion PHYS 131 Ch. 3: FOC 1, 3, 4, 6, 14. Problems 9, 37, 41 & 71 and Ch. 4: FOC 1, 3, 5, 8. Problems 3, 5 & 16. Feb 8, 018 Ch. 3: FOC 1, 3, 4, 6, 14. 1. (a) The horizonal componen of he projecile

More information

Section 3.8, Mechanical and Electrical Vibrations

Section 3.8, Mechanical and Electrical Vibrations Secion 3.8, Mechanical and Elecrical Vibraions Mechanical Unis in he U.S. Cusomary and Meric Sysems Disance Mass Time Force g (Earh) Uni U.S. Cusomary MKS Sysem CGS Sysem fee f slugs seconds sec pounds

More information

Solution: b All the terms must have the dimension of acceleration. We see that, indeed, each term has the units of acceleration

Solution: b All the terms must have the dimension of acceleration. We see that, indeed, each term has the units of acceleration PHYS 54 Tes Pracice Soluions Spring 8 Q: [4] Knowing ha in he ne epression a is acceleraion, v is speed, is posiion and is ime, from a dimensional v poin of view, he equaion a is a) incorrec b) correc

More information

III. Direct evolution of the density: The Liouville Operator

III. Direct evolution of the density: The Liouville Operator Cem 564 Lecure 8 3mar From Noes 8 003,005,007, 009 TIME IN QUANTUM MECANICS. I Ouline I. Te ime dependen Scroedinger equaion; ime dependence of energy eigensaes II.. Sae vecor (wave funcion) ime evoluion

More information

Two Popular Bayesian Estimators: Particle and Kalman Filters. McGill COMP 765 Sept 14 th, 2017

Two Popular Bayesian Estimators: Particle and Kalman Filters. McGill COMP 765 Sept 14 th, 2017 Two Popular Bayesian Esimaors: Paricle and Kalman Filers McGill COMP 765 Sep 14 h, 2017 1 1 1, dx x Bel x u x P x z P Recall: Bayes Filers,,,,,,, 1 1 1 1 u z u x P u z u x z P Bayes z = observaion u =

More information

AP Calculus BC Chapter 10 Part 1 AP Exam Problems

AP Calculus BC Chapter 10 Part 1 AP Exam Problems AP Calculus BC Chaper Par AP Eam Problems All problems are NO CALCULATOR unless oherwise indicaed Parameric Curves and Derivaives In he y plane, he graph of he parameric equaions = 5 + and y= for, is a

More information

Zürich. ETH Master Course: L Autonomous Mobile Robots Localization II

Zürich. ETH Master Course: L Autonomous Mobile Robots Localization II Roland Siegwar Margaria Chli Paul Furgale Marco Huer Marin Rufli Davide Scaramuzza ETH Maser Course: 151-0854-00L Auonomous Mobile Robos Localizaion II ACT and SEE For all do, (predicion updae / ACT),

More information

LAB 6: SIMPLE HARMONIC MOTION

LAB 6: SIMPLE HARMONIC MOTION 1 Name Dae Day/Time of Lab Parner(s) Lab TA Objecives LAB 6: SIMPLE HARMONIC MOTION To undersand oscillaion in relaion o equilibrium of conservaive forces To manipulae he independen variables of oscillaion:

More information

Chapter 8 The Complete Response of RL and RC Circuits

Chapter 8 The Complete Response of RL and RC Circuits Chaper 8 The Complee Response of RL and RC Circuis Seoul Naional Universiy Deparmen of Elecrical and Compuer Engineering Wha is Firs Order Circuis? Circuis ha conain only one inducor or only one capacior

More information

Solutions to Assignment 1

Solutions to Assignment 1 MA 2326 Differenial Equaions Insrucor: Peronela Radu Friday, February 8, 203 Soluions o Assignmen. Find he general soluions of he following ODEs: (a) 2 x = an x Soluion: I is a separable equaion as we

More information

MATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence

MATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence MATH 433/533, Fourier Analysis Secion 6, Proof of Fourier s Theorem for Poinwise Convergence Firs, some commens abou inegraing periodic funcions. If g is a periodic funcion, g(x + ) g(x) for all real x,

More information

MOMENTUM CONSERVATION LAW

MOMENTUM CONSERVATION LAW 1 AAST/AEDT AP PHYSICS B: Impulse and Momenum Le us run an experimen: The ball is moving wih a velociy of V o and a force of F is applied on i for he ime inerval of. As he resul he ball s velociy changes

More information

Optima and Equilibria for Traffic Flow on a Network

Optima and Equilibria for Traffic Flow on a Network Opima and Equilibria for Traffic Flow on a Nework Albero Bressan Deparmen of Mahemaics, Penn Sae Universiy bressan@mah.psu.edu Albero Bressan (Penn Sae) Opima and equilibria for raffic flow 1 / 1 A Traffic

More information

u(x) = e x 2 y + 2 ) Integrate and solve for x (1 + x)y + y = cos x Answer: Divide both sides by 1 + x and solve for y. y = x y + cos x

u(x) = e x 2 y + 2 ) Integrate and solve for x (1 + x)y + y = cos x Answer: Divide both sides by 1 + x and solve for y. y = x y + cos x . 1 Mah 211 Homework #3 February 2, 2001 2.4.3. y + (2/x)y = (cos x)/x 2 Answer: Compare y + (2/x) y = (cos x)/x 2 wih y = a(x)x + f(x)and noe ha a(x) = 2/x. Consequenly, an inegraing facor is found wih

More information

4.6 One Dimensional Kinematics and Integration

4.6 One Dimensional Kinematics and Integration 4.6 One Dimensional Kinemaics and Inegraion When he acceleraion a( of an objec is a non-consan funcion of ime, we would like o deermine he ime dependence of he posiion funcion x( and he x -componen of

More information