Lecture 4 January 14, 2016
|
|
- Abraham Morgan
- 6 years ago
- Views:
Transcription
1 MATH 262/CME 372: Applied Fourier Analysis and Winer 206 Elemens of Modern Signal Processing Lecure 4 January 4, 206 Prof. Emmanuel Candes Scribe: Carlos A. Sing-Long; Edied by E. Candes & E. Baes Ouline Agenda: Uncerainy Principle. Weyl-Heisenberg Uncerainy Principle 2. Quanum mechanical inerpreaion Las Time: Moivaed by he boxcar funcion, whose Fourier ransform is no summable, we inroduced he space of square-inegrable funcions L 2 (R). We proved he Parseval-Plancherel heorem, which shows he Fourier ransform preserves he L 2 -inner produc and is herefore an isomery (modulo a facor 2π). Using his resul, we were able o exend he Fourier ransform o L 2 (R) and hus make sense of he inverse Fourier ransform of he boxcar funcion as a squareinegrable funcion. The properies we proved also apply o his exension, wih he remark ha some equaliies hold almos everywhere. We also defined he Fourier ransform in higher dimensions by simply ieraing he single-variable ransform. 2 Weyl-Heisenberg s Uncerainy Principle The uncerainy principle is commonly known in physics as saying ha one canno know simulaneously he posiion and he momenum of a paricular wih infinie precision. In fac, his saemen is an implicaion of ha mahemaical observaion ha f and ˆf canno boh be concenraed. This is somehing we can inuiively see from some examples. For insance, our calculaions involving Gaussians (Lecure ) revealed ha increasing or decreasing he sandard deviaion in he ime (or spaial) domain has he opposie effec in he frequency (or Fourier) domain. Fig. illusraes his phenomenon. Taking his example o he exreme cases of 0 or infinie sandard deviaion gives he dela funcion in one domain and he consan funcion in he oher. The former is compleely localized, while he laer is compleely delocalized. Even wihou hese examples, one can iniially undersand he uncerainy principle from he ideniy F () = af(a) ˆF (ω) = ˆf(ω/a). Tha is, dilaing (compressing) ime is equivalen o compressing (dilaing) frequency. Wha he uncerainy principle does is provide a precise, quaniaive saemen of his noion.
2 Time domain: The blue curve corresponds o a signal f concenraed in ime, whereas he red curve corresponds o a signal f 2 spread in ime. Frequency domain: The Fourier ransform ˆf is spread in frequency, whereas he Fourier ransform ˆf 2 concenraes in frequency. Figure : Visual evidence of he uncerainy principle To sae he resul, we need o make a definiion. For f L 2 (R), define he spread of f o be σ 2 (f) = inf 0 R f 2 ( 0 ) 2 f() 2 d. We can inerpre his quaniy in probabilisic erms by noing ha p() = f() 2 / f 2 is a probabiliy densiy funcion, as (i) p() 0, (ii) p() d =. Therefore, one can recognize σ 2 (f) as he variance of a random variable X wih densiy p(), where he minimizing 0 is exacly he mean 0 = EX = p() d. Then σ(f) is, of course, he sandard deviaion of X. Similarly, q(ω) = ˆf(ω) 2 / ˆf 2 is he densiy of random variable Y having mean ω 0 = EY = ωq(ω) dω and sandard deviaion σ( ˆf). Wih his noaion, we can now sae he uncerainy principle: Theorem (Weyl-Heisenberg Uncerainy Principle). If f L 2 (R) is no idenically 0, hen (i) σ(f)σ( ˆf) /2. 2
3 (ii) Equaliy holds if and only if f is a ranslaion and modulaion of a Gaussian funcion. Thus, in erms of simulaneous ime and frequency concenraion, Gaussians are maximal wih respec o mean squared-deviaion. Proof. Le us begin wih some simplifying reducions. Firs, we may assume ha X from above is cenered (EX = 0), since ranslaions and modulaions do no affec spread: F = f( + 0 ) σ 2 (F ) = σ 2 (f) and σ 2 ( ˆF ) = σ 2 (e iω 0 ˆf) = σ 2 ( ˆf). Second, we may assume ha Y is also cenered, for he same reason: F () = e iω 0 f() σ 2 (F ) = σ 2 (f) and σ 2 ( ˆF ) = σ 2 ( ˆf(ω + ω 0 )) = σ 2 ( ˆf). Third, f by f/ f does no change spread, and he densiy funcions p() and q(ω) are lef unchanged. So we may assume f has norm. Numerically, our assumpions ell us (i) f() 2 d = ω ˆf(ω) 2 dω = 0, and (ii) f() 2 d =. Finally, smooh funcions are dense in L 2 (R), and he Fourier ransform preserves L 2 -norm (modulo a facor of 2π). So i suffices o demonsrae he uncerainy principle for smooh, in paricular differeniable, funcions. Now we will express σ( ˆf) in erms of f. Noe ha f () F iω ˆf(ω), and hen he Parseval-Plancherel heorem yields σ 2 ( ˆf) = ˆf ω 2 ˆf(ω) 2 dω = iω 2 2π ˆf(ω) 2 dω = f () 2 d, where we used ha f 2 = implies ˆf 2 = 2π. Consequenly, he uncerainy principle is equivalen o 2 f() 2 d f () 2 d 4, () meaning i is also a saemen abou a rade-off beween concenraion and regulariy. Namely, a funcion canno simulaneously be very concenraed and have small derivaives (in he L 2 sense). To complee he proof by checking (), we will employ he wo mos ubiquious ools of analysis: he Cauchy-Schwarz inequaliy and inegraion by pars. Recall he saemen of Cauchy-Schwarz: For g, h L 2 (R), g()h() d = g, h g h = g() 2 d h() 2 d. Following Weyl s proof (see [2], p. 393), we assume ha f() 2 decays faser han / a infiniy: f() 2 0 as. (2) 3
4 When f L 2 (R), (2) holds; he verificaion is echnical, and he ineresed reader can find i a he end of he proof. When f / L 2 (R), we have σ( ˆf) =, and so he uncerainy inequaliy cerainly holds. Apply Cauchy-Schwarz o g() = f() and h() = f (). Since g = σ(f) and h = σ( ˆf), we have f()f () d σ(f)σ( ˆf). We also have f()f () d σ(f)σ( ˆf) and so f()f () + f()f () d 2σ(f)σ( ˆf) Now inegraion by pars gives f()f () d + f()f () d = d d f() 2 d = f() 2 f() 2 d =. (from (2) and (ii)) We have hus shown σ(f)σ( ˆf) /2. Equaliy holds in Cauchy-Schwarz if and only if g h. So he uncerainy principle is saisfied wih equaliy if and only if f() f (). In his case, one can solve he differenial equaion o see f () = f() (3) σ2 f() = f(0)e 2 /2σ 2. Since we allow ranslaions, modulaions, and rescalings (i.e. undoing our reducions from he beginning), in general we have f() e iω0 e ( 0 ) 2 2σ 2. Tha is, f is a modulaion and ranslaion of a Gaussian funcion, as claimed. Remark : Assume f L 2 (R) is differeniable, wih f L 2 (R). Verificaion of (2) proceeds as follows: Suppose agains he claim ha here exiss ϵ > 0 such ha f() 2 > ϵ for arbirarily large. By possibly replacing f() wih f( ), we may assume his holds for arbirarily large posiive. Since f L 2 (R), we may choose so large ha f (u) 2 du 4. Applying Cauchy-Schwarz, we find ha for any, 2 ( 2 f( ) f() 2 = f (u) du f (u) du) ( ) f (u) 2 du. 4
5 So for [, + ϵ/], Hence for such. I follows ha Bu f L 2 (R), and so u 2 f(u) 2 du Therefore, (4) is a conradicion. f( ) f() 2 ϵ 4. f( ) f() 2 + ϵ lim u 2 f(u) 2 du > ϵ > ϵ 2 + ϵ u 2 f(u) 2 du = ϵ d = ϵ2 4. (4) Remark 2: If he RHS of (3) were given posiive sign, f() would be an exponenial wih posiive exponen and hus no in L 2 (R). Remark 3: We can verify direcly ha Gaussians saisfy he uncerainy principle wih equaliy: and Indeed, f() = e 2 /2 ˆf(ω) = 2πe ω2 /2 = f() 2 = e 2 = σ 2 (f) = 2, = ˆf(ω) 2 = 2πe ω2 = σ 2 ( ˆf) = 2. σ 2 (f)σ 2 ( ˆf) = 4. 3 Quanum mechanical inerpreaion Le us consider he descripion of a simple physical sysem, such as an elecron. From a classical poin of view, he physical sae of he elecron a any insan is characerized by six quaniies, ha is, is hree posiion coordinaes, and is hree momenum coordinaes. In oher words, o each ime insan we associae a poin in he phase space. Physical quaniies of ineres, such as he energy, can be measured from hese quaniies. However, in quanum mechanics he physical descripion is characerized by a sae vecor, Sae vecor: ψ H, where H is an infinie-dimensional Hilber space. The quaniies of ineres are represened by observables, which are symmeric operaors acing on he Hilber space Observable: A L(H), A = A. For insance, he posiion X and momenum P are observables. In classical mechanics, hese only exrac he posiion coordinaes or momenum coordinaes from a poin in he phase space, bu in quanum mechanics hey are symmeric linear operaors acing on H. 5
6 Since hey are symmeric, hey can be diagonalized and have a complee se of eigenvecors and real eigenvalues. We can define x eigenvecor of X wih eigenvalue x, ha is, X x = x x, and p eigenvecor of P wih eigenvalue p, ha is, P p = p p. As usual, we le eigenvecors be normalized. Since { x } x and { p } p are boh complee orhonormal ses, we can decompose he ideniy operaor as I = dx x x = dp p p. Fundamenal posulae of Quanum Mechanics: For an observable A wih eigensaes a and a sae ψ, he oucome of measuring A is a wih probabiliy a ψ 2 and he sysem ransiions o he sae a, ha is, ψ measure A a wih probabiliy a ψ 2. An experimen ha shows his behavior is Sern-Gerlach experimen. For a descripion of his experimen and is significance, you can see hp://en.wikipedia.org/wiki/sern-gerlach_ experimen. The original paper by Oo Sern and Walher Gerlach (in German) can be found here hp://link.springer.com/aricle/0.007%2fbf Also he fundamenal posulae is par of wha is known as he Copenhagen inerpreaion of quanum mechanics (see hp://en.wikipedia.org/wiki/copenhagen_inerpreaion#principles). 3. Uncerainy relaions Since he oucomes of measuring an obsevable are random variables, we can define he expecaion of an observable. Le a denoe he se of eigenvecors of he observable A and assume he sysem under consideraion is on a sae ψ. Then is expeced value is ā = A = da a a, ψ 2 and is variance ( A) 2 = da (a ā) 2 a, ψ 2. Heisenberg s uncerainy relaion saes ha: ( X) 2 ( P ) 2 (ħ/2) 2 (5) wih ( X) 2 = ( P ) 2 = dx (x x) 2 x, ψ 2 dp (p p) 2 p, ψ 2. 6
7 In oher words, he produc of sandard deviaions is a leas ħ/2, where ħ is he reduced Planck consan h/2π. As a remark, x, ψ is called he he wave funcion and ofen denoed by Ψ(x). These are he coefficiens of our sae ψ in he basis of eigenvecors for X. Clearly Ψ(x) 2 dx = and Ψ(x) 2 is usually inerpreed as he probabiliy finding he sysem (say, an elecron) on an (infiniesimal) neighborhood of x. To see why Heisenberg s uncerainy relaion holds, we make use of his: Claim 2. If x is an eigensae of X wih eigenvalue x and, likewise, p is an eigensae of P wih eigenvalue p, hen he do produc obeys x, p = 2πħ e i ħ px. Define he wave funcion Φ(p) in momemum space via Φ(p) = p, ψ. Since ψ = x x, ψ dx, Φ and Ψ are relaed by Φ(p) }{{} wave funcion in momenum space = p ψ = = x ψ p x dx 2πħ Ψ(x) e i ħ px dx. }{{} wave funcion in posiion space Similarly, we can wrie Ψ(x) }{{} wave funcion in posiion space = x ψ = = p ψ x p dx 2πħ Φ(p) e i ħ px dp. }{{} wave funcion in momenum space The Fourier ransform arises naurally, and can be hough as a change of variables beween posiion space and momenum space. Hence, we can wrie uncerainies abou posiion and momemum as ( X) 2 = (x x) Ψ(x) 2 dx ( P ) 2 = (p p) ˆΨ(p/ħ) 2 dp; from here, Theorem proved earlier gives Heisenberg s uncerainy relaion (5). 7
8 3.2 Momemum We need o argue in favor of Claim 2, and in order o do his, we inroduce he momemum operaor. We begin by considering a special class of maps acing on saes, he ranslaion map U(a) for a R defined as U(a) x x + a. This family of maps defines a group. We can consider is behavior for small δa, ha is U(δa) = I i ħ P δa + o(δa2 ). The operaor appearing on he firs-order erm of he above expansion is by definiion he momenum operaor. Roughly speaking, he momenum operaor causes an infiniesimal displacemen on a physical sae. Informally, U(δa) ψ Ψ(x) = Ψ(x) Ψ (x)δa + o(δa 2 ), and consequenly he momenum operaor is such ha Ψ(x) iħ d dxψ(x). Using our noaion, is acion on a sae ψ is given by ( ) d P ψ = iħ dx x ψ x dx. Wih ψ = p and using P p = p p, his gives ( ) d p p = iħ dx x p x dx. Taking he do produc wih he bra x gives ( ) d p x p = iħ dx x p x x dx = iħ d x p. dx This is a differenial equaion wih soluion x p = 2πħ e i ħ px, and consequenly he eigensaes of he momenum represened in he posiion basis are complex exponenials wih (spaial) frequency given by he momenum eigenvalue p. 3.3 Heisenberg s original formulaion In his 927 paper [], Heisenberg gave a slighly differen formulaion of all of his. Firs noe ha he mean oucome when applying A is given by ( ) A = ψ A ψ = ψ A da a a ψ = da ψ A a a ψ = da a a, ψ 2. In linear algebra, we would wrie ψ A ψ as ψ Aψ. 8
9 Se A = A A I, so ha ( A) 2 = A 2 A 2, is jus he variance of he oucome as before. Las, he commuaor or Lie bracke is [A, B] = AB BA, which is a measure, in some sense, of how incoheren he observables A and B are. Theorem 3 (Heisenberg (927)). For a sae ψ and any observables A and B we have ( A) 2 ( B) 2 4 [A, B] 2. This mahemaical inequaliy is also a sor of Cauchy-Schwarz inequaliy and I will omi he proof. Claim 4. We have Hence, [X, P ] 2 = ħ 2, from where we obain [X, P ] = iħ I. ( X) 2 ( P ) 2 ħ2 4. Proof. We compue he acion of XP in he basis of eigensaes of X, and recall ha X x = x x and P ψ = iħ d dx x, ψ x dx. On he one hand, XP ψ = iħ x d ( x, ψ ) x dx dx while on he oher I hus follows ha [X, P ] ψ = iħ P X ψ = iħ d (x x, ψ ) x dx. dx ( x d ) d x, ψ + (x x, ψ ) x dx = iħ dx dx x, ψ x dx = iħ ψ. References [] Werner Heisenberg. Über den anschaulichen inhal der quanenheoreischen kinemaik und mechanik. Zeischrif für Physik, 43(3-4):72 98, 927. [2] Hermann Weyl. The heory of groups and quanum mechanics. Courier Corporaion,
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 informationLinear 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 information4 Sequences of measurable functions
4 Sequences of measurable funcions 1. Le (Ω, A, µ) be a measure space (complee, afer a possible applicaion of he compleion heorem). In his chaper we invesigae relaions beween various (nonequivalen) convergences
More informationFinish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!
MAT 257, Handou 6: Ocober 7-2, 20. I. Assignmen. Finish reading Chaper 2 of Spiva, rereading earlier secions as necessary. handou and fill in some missing deails! II. Higher derivaives. Also, read his
More informationClass 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 informationMATH 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 informationTwo 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 informationAn 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 informationMath 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 informationLet us start with a two dimensional case. We consider a vector ( x,
Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our
More information23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes
Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals
More informationt 2 B F x,t n dsdt t u x,t dxdt
Evoluion Equaions For 0, fixed, le U U0, where U denoes a bounded open se in R n.suppose ha U is filled wih a maerial in which a conaminan is being ranspored by various means including diffusion and convecion.
More informationLecture 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 information1 Review of Zero-Sum Games
COS 5: heoreical Machine Learning Lecurer: Rob Schapire Lecure #23 Scribe: Eugene Brevdo April 30, 2008 Review of Zero-Sum Games Las ime we inroduced a mahemaical model for wo player zero-sum games. Any
More informationLecture 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 informationT L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB
Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal
More informationChapter 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 informationPhysics 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 information10. State Space Methods
. Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he
More informationMatrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality
Marix Versions of Some Refinemens of he Arihmeic-Geomeric Mean Inequaliy Bao Qi Feng and Andrew Tonge Absrac. We esablish marix versions of refinemens due o Alzer ], Carwrigh and Field 4], and Mercer 5]
More informationTraveling 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 information6. 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 informationNotes for Lecture 17-18
U.C. Berkeley CS278: Compuaional Complexiy Handou N7-8 Professor Luca Trevisan April 3-8, 2008 Noes for Lecure 7-8 In hese wo lecures we prove he firs half of he PCP Theorem, he Amplificaion Lemma, up
More informationLecture Notes 2. The Hilbert Space Approach to Time Series
Time Series Seven N. Durlauf Universiy of Wisconsin. Basic ideas Lecure Noes. The Hilber Space Approach o Time Series The Hilber space framework provides a very powerful language for discussing he relaionship
More informationHamilton- 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 informationKEY. Math 334 Midterm III Winter 2008 section 002 Instructor: Scott Glasgow
KEY Mah 334 Miderm III Winer 008 secion 00 Insrucor: Sco Glasgow Please do NOT wrie on his exam. No credi will be given for such work. Raher wrie in a blue book, or on your own paper, preferably engineering
More informationWavelet Methods for Time Series Analysis. What is a Wavelet? Part I: Introduction to Wavelets and Wavelet Transforms. sines & cosines are big waves
Wavele Mehods for Time Series Analysis Par I: Inroducion o Waveles and Wavele Transforms waveles are analysis ools for ime series and images as a subjec, waveles are relaively new (983 o presen) a synhesis
More informationTHE WAVE EQUATION. part hand-in for week 9 b. Any dilation v(x, t) = u(λx, λt) of u(x, t) is also a solution (where λ is constant).
THE WAVE EQUATION 43. (S) Le u(x, ) be a soluion of he wave equaion u u xx = 0. Show ha Q43(a) (c) is a. Any ranslaion v(x, ) = u(x + x 0, + 0 ) of u(x, ) is also a soluion (where x 0, 0 are consans).
More information2. 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 informationDiebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles
Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance
More informationEchocardiography Project and Finite Fourier Series
Echocardiography Projec and Finie Fourier Series 1 U M An echocardiagram is a plo of how a porion of he hear moves as he funcion of ime over he one or more hearbea cycles If he hearbea repeas iself every
More informationOrdinary Differential Equations
Ordinary Differenial Equaions 5. Examples of linear differenial equaions and heir applicaions We consider some examples of sysems of linear differenial equaions wih consan coefficiens y = a y +... + a
More informationConvergence of the Neumann series in higher norms
Convergence of he Neumann series in higher norms Charles L. Epsein Deparmen of Mahemaics, Universiy of Pennsylvania Version 1.0 Augus 1, 003 Absrac Naural condiions on an operaor A are given so ha he Neumann
More informationEXERCISES FOR SECTION 1.5
1.5 Exisence and Uniqueness of Soluions 43 20. 1 v c 21. 1 v c 1 2 4 6 8 10 1 2 2 4 6 8 10 Graph of approximae soluion obained using Euler s mehod wih = 0.1. Graph of approximae soluion obained using Euler
More informationHamilton- 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 informationGuest Lectures for Dr. MacFarlane s EE3350 Part Deux
Gues Lecures for Dr. MacFarlane s EE3350 Par Deux Michael Plane Mon., 08-30-2010 Wrie name in corner. Poin ou his is a review, so I will go faser. Remind hem o go lisen o online lecure abou geing an A
More informationAn random variable is a quantity that assumes different values with certain probabilities.
Probabiliy The probabiliy PrA) of an even A is a number in [, ] ha represens how likely A is o occur. The larger he value of PrA), he more likely he even is o occur. PrA) means he even mus occur. PrA)
More informationThe 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 informationt 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 informationSolutions from Chapter 9.1 and 9.2
Soluions from Chaper 9 and 92 Secion 9 Problem # This basically boils down o an exercise in he chain rule from calculus We are looking for soluions of he form: u( x) = f( k x c) where k x R 3 and k is
More informationACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H.
ACE 56 Fall 005 Lecure 5: he Simple Linear Regression Model: Sampling Properies of he Leas Squares Esimaors by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Inference in he Simple
More information23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes
Half-Range Series 2.5 Inroducion In his Secion we address he following problem: Can we find a Fourier series expansion of a funcion defined over a finie inerval? Of course we recognise ha such a funcion
More informationVehicle Arrival Models : Headway
Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where
More informationLie Derivatives operator vector field flow push back Lie derivative of
Lie Derivaives The Lie derivaive is a mehod of compuing he direcional derivaive of a vecor field wih respec o anoher vecor field We already know how o make sense of a direcional derivaive of real valued
More informationCourse Notes for EE227C (Spring 2018): Convex Optimization and Approximation
Course Noes for EE7C Spring 018: Convex Opimizaion and Approximaion Insrucor: Moriz Hard Email: hard+ee7c@berkeley.edu Graduae Insrucor: Max Simchowiz Email: msimchow+ee7c@berkeley.edu Ocober 15, 018 3
More informationThe Strong Law of Large Numbers
Lecure 9 The Srong Law of Large Numbers Reading: Grimme-Sirzaker 7.2; David Williams Probabiliy wih Maringales 7.2 Furher reading: Grimme-Sirzaker 7.1, 7.3-7.5 Wih he Convergence Theorem (Theorem 54) and
More informationBernoulli numbers. Francesco Chiatti, Matteo Pintonello. December 5, 2016
UNIVERSITÁ DEGLI STUDI DI PADOVA, DIPARTIMENTO DI MATEMATICA TULLIO LEVI-CIVITA Bernoulli numbers Francesco Chiai, Maeo Pinonello December 5, 206 During las lessons we have proved he Las Ferma Theorem
More informationMath 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 informationChallenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k
Challenge Problems DIS 03 and 0 March 6, 05 Choose one of he following problems, and work on i in your group. Your goal is o convince me ha your answer is correc. Even if your answer isn compleely correc,
More informationChapter 6. Systems of First Order Linear Differential Equations
Chaper 6 Sysems of Firs Order Linear Differenial Equaions We will only discuss firs order sysems However higher order sysems may be made ino firs order sysems by a rick shown below We will have a sligh
More informationInventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions
Muli-Period Sochasic Models: Opimali of (s, S) Polic for -Convex Objecive Funcions Consider a seing similar o he N-sage newsvendor problem excep ha now here is a fixed re-ordering cos (> 0) for each (re-)order.
More informationODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004
ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform
More informationLAPLACE TRANSFORM AND TRANSFER FUNCTION
CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dep. of Chemical and Biological Engineering 5-1 Road Map of he Lecure V Laplace Transform and Transfer funcions
More informationPhysics 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 informationChapter Three Systems of Linear Differential Equations
Chaper Three Sysems of Linear Differenial Equaions In his chaper we are going o consier sysems of firs orer orinary ifferenial equaions. These are sysems of he form x a x a x a n x n x a x a x a n x n
More informationHeat kernel and Harnack inequality on Riemannian manifolds
Hea kernel and Harnack inequaliy on Riemannian manifolds Alexander Grigor yan UHK 11/02/2014 onens 1 Laplace operaor and hea kernel 1 2 Uniform Faber-Krahn inequaliy 3 3 Gaussian upper bounds 4 4 ean-value
More informationk B 2 Radiofrequency pulses and hardware
1 Exra MR Problems DC Medical Imaging course April, 214 he problems below are harder, more ime-consuming, and inended for hose wih a more mahemaical background. hey are enirely opional, bu hopefully will
More information15. Vector Valued Functions
1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,
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.
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 informationHomework sheet Exercises done during the lecture of March 12, 2014
EXERCISE SESSION 2A FOR THE COURSE GÉOMÉTRIE EUCLIDIENNE, NON EUCLIDIENNE ET PROJECTIVE MATTEO TOMMASINI Homework shee 3-4 - Exercises done during he lecure of March 2, 204 Exercise 2 Is i rue ha he parameerized
More informationHamilton 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 informationRandom Processes 1/24
Random Processes 1/24 Random Process Oher Names : Random Signal Sochasic Process A Random Process is an exension of he concep of a Random variable (RV) Simples View : A Random Process is a RV ha is a Funcion
More informationSome 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 informationMath 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 informationTHE MYSTERY OF STOCHASTIC MECHANICS. Edward Nelson Department of Mathematics Princeton University
THE MYSTERY OF STOCHASTIC MECHANICS Edward Nelson Deparmen of Mahemaics Princeon Universiy 1 Classical Hamilon-Jacobi heory N paricles of various masses on a Euclidean space. Incorporae he masses in he
More informationSystem 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 informationAnn. Funct. Anal. 2 (2011), no. 2, A nnals of F unctional A nalysis ISSN: (electronic) URL:
Ann. Func. Anal. 2 2011, no. 2, 34 41 A nnals of F uncional A nalysis ISSN: 2008-8752 elecronic URL: www.emis.de/journals/afa/ CLASSIFICAION OF POSIIVE SOLUIONS OF NONLINEAR SYSEMS OF VOLERRA INEGRAL EQUAIONS
More informationHomework 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 informationLecture 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 informationAnnouncements: 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 informationLecture 6: Wiener Process
Lecure 6: Wiener Process Eric Vanden-Eijnden Chapers 6, 7 and 8 offer a (very) brief inroducion o sochasic analysis. These lecures are based in par on a book projec wih Weinan E. A sandard reference for
More informationConcourse Math Spring 2012 Worked Examples: Matrix Methods for Solving Systems of 1st Order Linear Differential Equations
Concourse Mah 80 Spring 0 Worked Examples: Marix Mehods for Solving Sysems of s Order Linear Differenial Equaions The Main Idea: Given a sysem of s order linear differenial equaions d x d Ax wih iniial
More informationApproximation Algorithms for Unique Games via Orthogonal Separators
Approximaion Algorihms for Unique Games via Orhogonal Separaors Lecure noes by Konsanin Makarychev. Lecure noes are based on he papers [CMM06a, CMM06b, LM4]. Unique Games In hese lecure noes, we define
More informationCHAPTER 2 Signals And Spectra
CHAPER Signals And Specra Properies of Signals and Noise In communicaion sysems he received waveform is usually caegorized ino he desired par conaining he informaion, and he undesired par. he desired par
More informationOptimality Conditions for Unconstrained Problems
62 CHAPTER 6 Opimaliy Condiions for Unconsrained Problems 1 Unconsrained Opimizaion 11 Exisence Consider he problem of minimizing he funcion f : R n R where f is coninuous on all of R n : P min f(x) x
More information( ) is the stretch factor, and x the
(Lecures 7-8) Liddle, Chaper 5 Simple cosmological models (i) Hubble s Law revisied Self-similar srech of he universe All universe models have his characerisic v r ; v = Hr since only his conserves homogeneiy
More informationThe 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 informationIMPLICIT AND INVERSE FUNCTION THEOREMS PAUL SCHRIMPF 1 OCTOBER 25, 2013
IMPLICI AND INVERSE FUNCION HEOREMS PAUL SCHRIMPF 1 OCOBER 25, 213 UNIVERSIY OF BRIISH COLUMBIA ECONOMICS 526 We have exensively sudied how o solve sysems of linear equaions. We know how o check wheher
More information( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is
UNIT IMPULSE RESPONSE, UNIT STEP RESPONSE, STABILITY. Uni impulse funcion (Dirac dela funcion, dela funcion) rigorously defined is no sricly a funcion, bu disribuion (or measure), precise reamen requires
More information14 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 informationOn Carlsson type orthogonality and characterization of inner product spaces
Filoma 26:4 (212), 859 87 DOI 1.2298/FIL124859K Published by Faculy of Sciences and Mahemaics, Universiy of Niš, Serbia Available a: hp://www.pmf.ni.ac.rs/filoma On Carlsson ype orhogonaliy and characerizaion
More information5.1 - Logarithms and Their Properties
Chaper 5 Logarihmic Funcions 5.1 - Logarihms and Their Properies Suppose ha a populaion grows according o he formula P 10, where P is he colony size a ime, in hours. When will he populaion be 2500? We
More informationKEY. Math 334 Midterm III Fall 2008 sections 001 and 003 Instructor: Scott Glasgow
KEY Mah 334 Miderm III Fall 28 secions and 3 Insrucor: Sco Glasgow Please do NOT wrie on his exam. No credi will be given for such work. Raher wrie in a blue book, or on your own paper, preferably engineering
More informationSZG Macro 2011 Lecture 3: Dynamic Programming. SZG macro 2011 lecture 3 1
SZG Macro 2011 Lecure 3: Dynamic Programming SZG macro 2011 lecure 3 1 Background Our previous discussion of opimal consumpion over ime and of opimal capial accumulaion sugges sudying he general decision
More information18 Biological models with discrete time
8 Biological models wih discree ime The mos imporan applicaions, however, may be pedagogical. The elegan body of mahemaical heory peraining o linear sysems (Fourier analysis, orhogonal funcions, and so
More informationLECTURE 1: GENERALIZED RAY KNIGHT THEOREM FOR FINITE MARKOV CHAINS
LECTURE : GENERALIZED RAY KNIGHT THEOREM FOR FINITE MARKOV CHAINS We will work wih a coninuous ime reversible Markov chain X on a finie conneced sae space, wih generaor Lf(x = y q x,yf(y. (Recall ha q
More informationMODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE
Topics MODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES 2-6 3. FUNCTION OF A RANDOM VARIABLE 3.2 PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE 3.3 EXPECTATION AND MOMENTS
More information6.003 Homework #9 Solutions
6.00 Homework #9 Soluions Problems. Fourier varieies a. Deermine he Fourier series coefficiens of he following signal, which is periodic in 0. x () 0 0 a 0 5 a k sin πk 5 sin πk 5 πk for k 0 a k 0 πk j
More informationCash 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 information6.2 Transforms of Derivatives and Integrals.
SEC. 6.2 Transforms of Derivaives and Inegrals. ODEs 2 3 33 39 23. Change of scale. If l( f ()) F(s) and c is any 33 45 APPLICATION OF s-shifting posiive consan, show ha l( f (c)) F(s>c)>c (Hin: In Probs.
More informationLecture 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 informationLecture #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 informationL07. 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 informationwhere 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 informationShort Introduction to Fractional Calculus
. Shor Inroducion o Fracional Calculus Mauro Bologna Deparameno de Física, Faculad de Ciencias Universidad de Tarapacá, Arica, Chile email: mbologna@ua.cl Absrac In he pas few years fracional calculus
More informationLecture 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 informationSimulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010
Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid
More informationO Q L N. Discrete-Time Stochastic Dynamic Programming. I. Notation and basic assumptions. ε t : a px1 random vector of disturbances at time t.
Econ. 5b Spring 999 C. Sims Discree-Time Sochasic Dynamic Programming 995, 996 by Chrisopher Sims. This maerial may be freely reproduced for educaional and research purposes, so long as i is no alered,
More informationSection 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients
Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous
More informationRANDOM LAGRANGE MULTIPLIERS AND TRANSVERSALITY
ECO 504 Spring 2006 Chris Sims RANDOM LAGRANGE MULTIPLIERS AND TRANSVERSALITY 1. INTRODUCTION Lagrange muliplier mehods are sandard fare in elemenary calculus courses, and hey play a cenral role in economic
More information