Laplace Transform Integration (ACM 40520)
|
|
- Derrick Singleton
- 5 years ago
- Views:
Transcription
1 Laplace Transform Integration (ACM 40520) Peter Lynch School of Mathematical Sciences
2 Outline Basic Theory Residue Theorem Ordinary Differential Equations General Vector NWP Equation Lagrangian Formulation
3 Outline Basic Theory Residue Theorem Ordinary Differential Equations General Vector NWP Equation Lagrangian Formulation
4 The Laplace Transform: Definition For a function of time f (t), the LT is defined as ˆf (s) = e st f (t) dt. 0 Here, s is complex and ˆf (s) is a complex function of s.
5 The Laplace Transform: Definition For a function of time f (t), the LT is defined as ˆf (s) = e st f (t) dt. 0 Here, s is complex and ˆf (s) is a complex function of s. The inversion from ˆf (s) back to f (t) is f (t) = 1 e st ˆf (s) ds. 2πi C 1 where C 1 is a contour in the s-plane, parallel to the imaginary axis, to the right of all singularities of ˆf (s).
6
7 Integral Transforms The LT is one of a large family of integral transforms They can be defined as F(s) = R K (s, t) f (t) dt where R is the range of f (t) and K (s, t) is called the kernel of the transform.
8 Integral Transforms The LT is one of a large family of integral transforms They can be defined as F(s) = R K (s, t) f (t) dt where R is the range of f (t) and K (s, t) is called the kernel of the transform. For example, the Fourier transform is f (ω) = e iωt f (t) dt f (t) = 1 2π e +iωt f (ω) dω
9 Integral Transforms The LT is one of a large family of integral transforms They can be defined as F(s) = R K (s, t) f (t) dt where R is the range of f (t) and K (s, t) is called the kernel of the transform. For example, the Fourier transform is f (ω) = e iωt f (t) dt f (t) = 1 2π e +iωt f (ω) dω The Hilbert transform is another... and many more.
10 The LT is a linear operator Therefore L{αf (t)} = L{f (t)} = ˆf (s) e st αf (t) dt = α e st f (t) dt. e st f (t) dt = αl{f (t)}.
11 The LT is a linear operator Therefore L{αf (t)} = Also L{f (t)} = ˆf (s) L{f (t)+g(t)} = e st αf (t) dt = α 0 e st f (t) dt. e st f (t) dt = αl{f (t)}. e st [f (t)+g(t)] dt = L{f (t)}+l{g(t)}.
12 The LT is a linear operator Therefore L{αf (t)} = Also L{f (t)} = ˆf (s) L{f (t)+g(t)} = Therefore e st αf (t) dt = α 0 e st f (t) dt. e st f (t) dt = αl{f (t)}. e st [f (t)+g(t)] dt = L{f (t)}+l{g(t)}. L{αf (t) + βg(t)} = αl{f (t)} + βl{g(t)}.
13 Basic Properties of the LT L{a} = a/s L{exp(at)} = 1/(s a)
14 Basic Properties of the LT L{a} = a/s L{exp(at)} = 1/(s a) L{exp(iωt)} = 1/(s iω)
15 Basic Properties of the LT L{a} = a/s L{exp(at)} = 1/(s a) L{exp(iωt)} = 1/(s iω) L{sin at} = a/(s 2 + a 2 ) L{cos at} = s/(s 2 + a 2 )
16 Basic Properties of the LT L{a} = a/s L{exp(at)} = 1/(s a) L{exp(iωt)} = 1/(s iω) L{sin at} = a/(s 2 + a 2 ) L{cos at} = s/(s 2 + a 2 ) L{df /dt} = sˆf (s) f (0)
17 Basic Properties of the LT L{a} = a/s L{exp(at)} = 1/(s a) L{exp(iωt)} = 1/(s iω) L{sin at} = a/(s 2 + a 2 ) L{cos at} = s/(s 2 + a 2 ) L{df /dt} = sˆf (s) f (0) All these results can be demonstrated immediately by using the definition of the Laplace transform L{f (t)}.
18 Outline Basic Theory Residue Theorem Ordinary Differential Equations General Vector NWP Equation Lagrangian Formulation
19 A Simple Oscillation Let f (t) have a single harmonic component f (t) = α exp(iωt)
20 A Simple Oscillation Let f (t) have a single harmonic component f (t) = α exp(iωt) Exercise: Show that the LT of f (t) is ˆf (s) = α s iω, a holomorphic function with a simple pole at s = iω.
21 A Simple Oscillation Let f (t) have a single harmonic component f (t) = α exp(iωt) Exercise: Show that the LT of f (t) is ˆf (s) = α s iω, a holomorphic function with a simple pole at s = iω. A pure (monochrome) oscillation in time transforms to a function with a single pole.
22 A Simple Oscillation Let f (t) have a single harmonic component f (t) = α exp(iωt) Exercise: Show that the LT of f (t) is ˆf (s) = α s iω, a holomorphic function with a simple pole at s = iω. A pure (monochrome) oscillation in time transforms to a function with a single pole. The position of the pole is determined by the frequency of the oscillation.
23 Again ˆf (s) = L{α exp(iωt)} = α s iω.
24 Again ˆf (s) = L{α exp(iωt)} = α s iω. The inverse transform of ˆf (s) is f (t) = 1 e st 1 ˆf (s) ds = 2πi C 1 2πi α exp(st) C 1 s iω ds.
25 Again ˆf (s) = L{α exp(iωt)} = α s iω. The inverse transform of ˆf (s) is f (t) = 1 e st 1 ˆf (s) ds = 2πi C 1 2πi α exp(st) C 1 s iω ds. We augment C 1 by a semi-circular arc C 2 in the left half-plane. Denote the resulting closed contour by C 0 = C 1 C 2.
26 Again ˆf (s) = L{α exp(iωt)} = α s iω. The inverse transform of ˆf (s) is f (t) = 1 e st 1 ˆf (s) ds = 2πi C 1 2πi α exp(st) C 1 s iω ds. We augment C 1 by a semi-circular arc C 2 in the left half-plane. Denote the resulting closed contour by C 0 = C 1 C 2. If it can be shown that this leaves the value of the integral unchanged, the inverse is an integral around a closed contour.
27 Again ˆf (s) = L{α exp(iωt)} = α s iω. The inverse transform of ˆf (s) is f (t) = 1 e st 1 ˆf (s) ds = 2πi C 1 2πi α exp(st) C 1 s iω ds. We augment C 1 by a semi-circular arc C 2 in the left half-plane. Denote the resulting closed contour by C 0 = C 1 C 2. If it can be shown that this leaves the value of the integral unchanged, the inverse is an integral around a closed contour. See textbook of Doetsch
28
29 For an integral around a closed contour, f (t) = 1 α exp(st) 2πi C 0 s iω ds, we can apply the residue theorem:
30 For an integral around a closed contour, f (t) = 1 α exp(st) 2πi C 0 s iω ds, we can apply the residue theorem: f (t) = C 0 [Residues of ( )] α exp(st) s iω so f (t) is the sum of the residues of the integrand within C 0.
31
32 Again f (t) = C 0 [Residues of ( )] α exp(st) s iω
33 Again f (t) = C 0 [Residues of ( )] α exp(st) s iω There is just one pole, at s = iω. The residue is ( ) α exp(st) lim (s iω) = α exp(iωt) s iω s iω
34 Again f (t) = C 0 [Residues of ( )] α exp(st) s iω There is just one pole, at s = iω. The residue is ( ) α exp(st) lim (s iω) = α exp(iωt) s iω s iω So we recover the input function: f (t) = α exp(iωt)
35 A Two-Component Oscillation Let f (t) have two harmonic components f (t) = a exp(iωt) + A exp(iωt) ω Ω
36 A Two-Component Oscillation Let f (t) have two harmonic components f (t) = a exp(iωt) + A exp(iωt) ω Ω The LT is a linear operator, so the transform of f (t) is ˆf (s) = a s iω + A s iω, which has two simple poles, at s = iω and s = iω.
37 A Two-Component Oscillation Let f (t) have two harmonic components f (t) = a exp(iωt) + A exp(iωt) ω Ω The LT is a linear operator, so the transform of f (t) is ˆf (s) = a s iω + A s iω, which has two simple poles, at s = iω and s = iω. The LF pole, at s = iω, is close to the origin. The HF pole, at s = iω, is far from the origin.
38 Again ˆf (s) = a s iω + A s iω.
39 Again ˆf (s) = a s iω + A s iω. The inverse transform of ˆf (s) is f (t) = 1 a exp(st) 2πi C 0 s iω ds + 1 A exp(st) 2πi C 0 s iω ds = a exp(iωt) + A exp(iωt).
40 Again ˆf (s) = a s iω + A s iω. The inverse transform of ˆf (s) is f (t) = 1 a exp(st) 2πi C 0 s iω ds + 1 A exp(st) 2πi C 0 s iω ds = a exp(iωt) + A exp(iωt). We now replace C 0 by a circular contour C centred at the origin, with radius γ such that ω < γ < Ω.
41
42 Again: We replace C 0 by C with ω < γ < Ω.
43 Again: We replace C 0 by C with ω < γ < Ω. We denote the modified operator by L.
44 Again: We replace C 0 by C with ω < γ < Ω. We denote the modified operator by L. Since the pole s = iω falls within the contour C, it contributes to the integral. Since the pole s = iω falls outside the contour C, it makes no contribution.
45 Again: We replace C 0 by C with ω < γ < Ω. We denote the modified operator by L. Since the pole s = iω falls within the contour C, it contributes to the integral. Since the pole s = iω falls outside the contour C, it makes no contribution. Therefore, f (t) L {ˆf (s)} = 1 2πi a exp(st) C s iω ds = a exp(iωt).
46 Again: We replace C 0 by C with ω < γ < Ω. We denote the modified operator by L. Since the pole s = iω falls within the contour C, it contributes to the integral. Since the pole s = iω falls outside the contour C, it makes no contribution. Therefore, f (t) L {ˆf (s)} = 1 2πi a exp(st) C s iω ds = a exp(iωt). We have filtered f (t): the function f (t) is the LF component of f (t). The HF component is gone.
47 Exercise Consider the test function f (t) = α 1 cos(ω 1 t ψ 1 ) + α 2 cos(ω 2 t ψ 2 ) ω 1 < ω 2
48 Exercise Consider the test function f (t) = α 1 cos(ω 1 t ψ 1 ) + α 2 cos(ω 2 t ψ 2 ) ω 1 < ω 2 Show that the LT is [ α 1 e iψ 1 ] ˆf (s) = + eiψ1 2 s iω 1 s + iω 1 + α 2 2 [ e iψ 2 s iω 2 + eiψ2 s + iω 2 ]
49 Exercise Consider the test function f (t) = α 1 cos(ω 1 t ψ 1 ) + α 2 cos(ω 2 t ψ 2 ) ω 1 < ω 2 Show that the LT is [ α 1 e iψ 1 ] ˆf (s) = + eiψ1 2 s iω 1 s + iω 1 + α 2 2 [ e iψ 2 s iω 2 + eiψ2 s + iω 2 ] Show how, by choosing C with ω 1 < γ < ω 2, the HF component can be eliminated.
50
51 Outline Basic Theory Residue Theorem Ordinary Differential Equations General Vector NWP Equation Lagrangian Formulation
52 Applying LT to an ODE We consider a nonlinear ordinary differential equation dw dt + iωw + n(w) = 0 w(0) = w 0
53 Applying LT to an ODE We consider a nonlinear ordinary differential equation dw dt + iωw + n(w) = 0 w(0) = w 0 The LT of the equation is (sŵ w 0 ) + iωŵ + n 0 s = 0. We have frozen n(w) at its initial value n 0 = n(w 0 ).
54 Applying LT to an ODE We consider a nonlinear ordinary differential equation dw dt + iωw + n(w) = 0 w(0) = w 0 The LT of the equation is (sŵ w 0 ) + iωŵ + n 0 s = 0. We have frozen n(w) at its initial value n 0 = n(w 0 ). We can immediately solve for the transform solution: ŵ(s) = 1 [ w 0 n ] 0 s + iω s
55 Using partial fractions, we write the transform as ( ) w0 ŵ(s) = + n ( 0 1 s + iω iω s 1 ) s + iω There are two poles, at s = iω and at s = 0.
56 Using partial fractions, we write the transform as ( ) w0 ŵ(s) = + n ( 0 1 s + iω iω s 1 ) s + iω There are two poles, at s = iω and at s = 0. The pole at s = 0 always falls within the contour C. The pole at s = iω may or may not fall within C.
57 Using partial fractions, we write the transform as ( ) w0 ŵ(s) = + n ( 0 1 s + iω iω s 1 ) s + iω There are two poles, at s = iω and at s = 0. The pole at s = 0 always falls within the contour C. The pole at s = iω may or may not fall within C. Thus, the solution is ( w 0 n ) 0 w (t) = iω exp( iωt) + n 0 iω : ω < γ n 0 iω : ω > γ
58 Again, ( w 0 n 0 w (t) = iω ) exp( iωt) + n 0 iω : ω < γ n 0 : ω > γ iω
59 Again, ( w 0 n 0 w (t) = iω ) exp( iωt) + n 0 iω : ω < γ n 0 : ω > γ iω So we see that, for a LF oscillation ( ω < γ), the solution w (t) is the same as the full solution w(t) of the ODE.
60 Again, ( w 0 n 0 w (t) = iω ) exp( iωt) + n 0 iω : ω < γ n 0 : ω > γ iω So we see that, for a LF oscillation ( ω < γ), the solution w (t) is the same as the full solution w(t) of the ODE. For a HF oscillation ( ω > γ), the solution contains only a constant term.
61 Again, ( w 0 n 0 w (t) = iω ) exp( iωt) + n 0 iω : ω < γ n 0 : ω > γ iω So we see that, for a LF oscillation ( ω < γ), the solution w (t) is the same as the full solution w(t) of the ODE. For a HF oscillation ( ω > γ), the solution contains only a constant term. Thus, high frequencies are filtered out.
62 Outline Basic Theory Residue Theorem Ordinary Differential Equations General Vector NWP Equation Lagrangian Formulation
63 A General Vector Equation We write the general NWP equations symbolically as dx dt + i LX + N(X) = 0 where X(t) is the state vector at time t.
64 A General Vector Equation We write the general NWP equations symbolically as dx dt + i LX + N(X) = 0 where X(t) is the state vector at time t. We apply the Laplace transform to get (s ˆX X 0 ) + i LˆX + 1 s N 0 = 0 where X 0 is the initial value of X and N 0 = N(X 0 ) is held constant at its initial value.
65 A General Vector Equation We write the general NWP equations symbolically as dx dt + i LX + N(X) = 0 where X(t) is the state vector at time t. We apply the Laplace transform to get (s ˆX X 0 ) + i LˆX + 1 s N 0 = 0 where X 0 is the initial value of X and N 0 = N(X 0 ) is held constant at its initial value. The frequencies are entangled. How do we proceed?
66 Eigenanalysis Ẋ + i LX + N(X) = 0
67 Eigenanalysis Ẋ + i LX + N(X) = 0 Assume the eigenanalysis of L is LE = EΛ where Λ = diag(λ 1,..., λ N ) and E = (e 1,..., e N ).
68 Eigenanalysis Ẋ + i LX + N(X) = 0 Assume the eigenanalysis of L is LE = EΛ where Λ = diag(λ 1,..., λ N ) and E = (e 1,..., e N ). More explicitly, assume that the eigenfrequencies split in two: [ ] ΛY 0 Λ = 0 Λ Z Λ Y Λ Z : Frequencies of rotational modes (LF) : Frequencies of gravity-inertia modes (HF)
69 We define a new set of variables: W = E 1 X.
70 We define a new set of variables: W = E 1 X. Multiplying the equation by E 1 we get E 1 Ẋ + i E 1 L(EE 1 )X + E 1 N(X) = 0
71 We define a new set of variables: W = E 1 X. Multiplying the equation by E 1 we get E 1 Ẋ + i E 1 L(EE 1 )X + E 1 N(X) = 0 This is just Ẇ + i ΛW + E 1 N(X) = 0
72 We define a new set of variables: W = E 1 X. Multiplying the equation by E 1 we get E 1 Ẋ + i E 1 L(EE 1 )X + E 1 N(X) = 0 This is just Ẇ + i ΛW + E 1 N(X) = 0 This equation separates into two sub-systems: Ẏ + i Λ Y Y + N Y (Y, Z) = 0 Ż + i Λ Z Z + N Z (Y, Z) = 0 where W = (Y, Z) T.
73 We define a new set of variables: W = E 1 X. Multiplying the equation by E 1 we get E 1 Ẋ + i E 1 L(EE 1 )X + E 1 N(X) = 0 This is just Ẇ + i ΛW + E 1 N(X) = 0 This equation separates into two sub-systems: where W = (Y, Z) T. Ẏ + i Λ Y Y + N Y (Y, Z) = 0 Ż + i Λ Z Z + N Z (Y, Z) = 0 The variables Y and Z are all coupled through the nonlinear terms N Y (Y, Z) and N Z (Y, Z).
74 We first consider a single component w: ẇ + iλw + n(w) = 0 Note that all other components may occur in the nonlinear term.
75 We first consider a single component w: ẇ + iλw + n(w) = 0 Note that all other components may occur in the nonlinear term. Holding the nonlinear term constant, we get or (sŵ w 0 ) + iλŵ + n 0 s = 0 ŵ(s) = 1 [ w 0 n ] 0 s + iλ s
76 We first consider a single component w: ẇ + iλw + n(w) = 0 Note that all other components may occur in the nonlinear term. Holding the nonlinear term constant, we get or (sŵ w 0 ) + iλŵ + n 0 s = 0 ŵ(s) = 1 [ w 0 n ] 0 s + iλ s This has two poles, at s = 0 and s = iλ: ŵ(s) = w 0 + n 0 /iλ s + iλ n 0/iλ s
77 Again, ŵ(s) = w 0 + n 0 /iλ s + iλ n 0/iλ s
78 Again, ŵ(s) = w 0 + n 0 /iλ s + iλ n 0/iλ s If λ is small, both poles are within C, so ( w (t) = L {ŵ} = w 0 + n ) ( 0 e iλt n0 ) iλ iλ an oscillation with frequency λ.,
79 Again, ŵ(s) = w 0 + n 0 /iλ s + iλ n 0/iλ s If λ is small, both poles are within C, so ( w (t) = L {ŵ} = w 0 + n ) ( 0 e iλt n0 ) iλ iλ an oscillation with frequency λ. If λ is large, the pole at s = iλ falls outside C, and ( w (t) = L n0 ) {ŵ} =. iλ,
80 Again, ŵ(s) = w 0 + n 0 /iλ s + iλ n 0/iλ s If λ is small, both poles are within C, so ( w (t) = L {ŵ} = w 0 + n ) ( 0 e iλt n0 ) iλ iλ an oscillation with frequency λ. If λ is large, the pole at s = iλ falls outside C, and ( w (t) = L n0 ) {ŵ} =. iλ This corresponds to putting ẇ = 0 in the equation: iλw + n 0 = 0,
81 General Solution Method We recall that the Laplace transform of the equation is (s ˆX X 0 ) + i LˆX + 1 s N 0 = 0 where X 0 is the initial value of X and N 0 = N(X 0 ) is held constant at its initial value.
82 General Solution Method We recall that the Laplace transform of the equation is (s ˆX X 0 ) + i LˆX + 1 s N 0 = 0 where X 0 is the initial value of X and N 0 = N(X 0 ) is held constant at its initial value. But now we take n t to be the initial time: (s ˆX X n ) + i LˆX + 1 s Nn = 0
83 General Solution Method We recall that the Laplace transform of the equation is (s ˆX X 0 ) + i LˆX + 1 s N 0 = 0 where X 0 is the initial value of X and N 0 = N(X 0 ) is held constant at its initial value. But now we take n t to be the initial time: (s ˆX X n ) + i LˆX + 1 s Nn = 0 The solution can be written formally: ˆX(s) = (s I + i L) 1 [ X n 1 s Nn ]
84 Again, ˆX(s) = (s I + i L) 1 [ X n 1 s Nn ]
85 Again, ˆX(s) = (s I + i L) 1 [ X n 1 s Nn ] We recover the filtered solution by applying L at time (n + 1) t. X = L {ˆX(s)} t= t
86 Again, ˆX(s) = (s I + i L) 1 [ X n 1 s Nn ] We recover the filtered solution by applying L at time (n + 1) t. X = L {ˆX(s)} t= t The procedure may now be iterated to produce a forecast of any length.
87 Again, ˆX(s) = (s I + i L) 1 [ X n 1 s Nn ] We recover the filtered solution by applying L at time (n + 1) t. X = L {ˆX(s)} t= t The procedure may now be iterated to produce a forecast of any length. Further details are given in Clancy and Lynch, 2011a,b
88 Outline Basic Theory Residue Theorem Ordinary Differential Equations General Vector NWP Equation Lagrangian Formulation
89 Lagrangian Formulation We now consider how to combine the Laplace transform approach with Lagrangian advection.
90 Lagrangian Formulation We now consider how to combine the Laplace transform approach with Lagrangian advection. The general form of the equation is DX Dt + i LX + N(X) = 0 where advection is now included in the time derivative.
91 Lagrangian Formulation We now consider how to combine the Laplace transform approach with Lagrangian advection. The general form of the equation is DX Dt + i LX + N(X) = 0 where advection is now included in the time derivative. We re-define the Laplace transform to be the integral in time along the trajectory of a fluid parcel: ˆX(s) e st X(t) dt T
92 We consider parcels that arrive at the gridpoints at time (n + 1) t. They originate at locations not corresponding to gridpoints at time n t.
93 We consider parcels that arrive at the gridpoints at time (n + 1) t. They originate at locations not corresponding to gridpoints at time n t. We denote the value at an arrival point by X n+1 A. The value at the departure point is X n D.
94 We consider parcels that arrive at the gridpoints at time (n + 1) t. They originate at locations not corresponding to gridpoints at time n t. We denote the value at an arrival point by X n+1 A. The value at the departure point is X n D. The initial values when transforming the Lagrangian time derivatives are X n D.
95 We consider parcels that arrive at the gridpoints at time (n + 1) t. They originate at locations not corresponding to gridpoints at time n t. We denote the value at an arrival point by X n+1 A. The value at the departure point is X n D. The initial values when transforming the Lagrangian time derivatives are X n D. The equations thus transform to (s ˆX X n D) + i LˆX + 1 s Nn+ 1 2 M = 0 where we evaluate nonlinear terms at a mid-point, interpolated in space and extrapolated in time.
96
97 The solution can be written formally: ˆX(s) = (s I + i L) 1 [X n D 1 s Nn+ 1 2 M ]
98 The solution can be written formally: ˆX(s) = (s I + i L) 1 [X n D 1 s Nn+ 1 2 M ] The values at the departure point and mid-point are computed by interpolation.
99 The solution can be written formally: ˆX(s) = (s I + i L) 1 [X n D 1 s Nn+ 1 2 M ] The values at the departure point and mid-point are computed by interpolation. We recover the filtered solution by applying L at time (n + 1) t, or t after the initial time. X = L {ˆX(s)} t= t
100 The solution can be written formally: ˆX(s) = (s I + i L) 1 [X n D 1 s Nn+ 1 2 M ] The values at the departure point and mid-point are computed by interpolation. We recover the filtered solution by applying L at time (n + 1) t, or t after the initial time. X = L {ˆX(s)} t= t The procedure may now be iterated to produce a forecast of any length.
101 The solution can be written formally: ˆX(s) = (s I + i L) 1 [X n D 1 s Nn+ 1 2 M ] The values at the departure point and mid-point are computed by interpolation. We recover the filtered solution by applying L at time (n + 1) t, or t after the initial time. X = L {ˆX(s)} t= t The procedure may now be iterated to produce a forecast of any length. Further details are given in Clancy and Lynch, 2011a,b
4. Complex Oscillations
4. Complex Oscillations The most common use of complex numbers in physics is for analyzing oscillations and waves. We will illustrate this with a simple but crucially important model, the damped harmonic
More information10. Operators and the Exponential Response Formula
52 10. Operators and the Exponential Response Formula 10.1. Operators. Operators are to functions as functions are to numbers. An operator takes a function, does something to it, and returns this modified
More informationIntroduction to Path Integrals
Introduction to Path Integrals Consider ordinary quantum mechanics of a single particle in one space dimension. Let s work in the coordinate space and study the evolution kernel Ut B, x B ; T A, x A )
More informationDamped Harmonic Oscillator
Damped Harmonic Oscillator Wednesday, 23 October 213 A simple harmonic oscillator subject to linear damping may oscillate with exponential decay, or it may decay biexponentially without oscillating, or
More informationFourier transforms, Generalised functions and Greens functions
Fourier transforms, Generalised functions and Greens functions T. Johnson 2015-01-23 Electromagnetic Processes In Dispersive Media, Lecture 2 - T. Johnson 1 Motivation A big part of this course concerns
More informationUnderstand the existence and uniqueness theorems and what they tell you about solutions to initial value problems.
Review Outline To review for the final, look over the following outline and look at problems from the book and on the old exam s and exam reviews to find problems about each of the following topics.. Basics
More informationTransient Phenomena in Quantum Bound States Subjected to a Sudden Perturbation
Symmetry, Integrability and Geometry: Methods and Applications Vol. (5), Paper 3, 9 pages Transient Phenomena in Quantum Bound States Subjected to a Sudden Perturbation Marcos MOSHINSKY and Emerson SADURNÍ
More informationProblems 3 (due 27 January)
Problems 3 due 7 January). You previously found singularities of the functions below. Here, identify the poles, order of these poles and residues for these functions: a) tanhz) = ez e z e z + e z Solution:
More informationUnit 2: Modeling in the Frequency Domain Part 2: The Laplace Transform. The Laplace Transform. The need for Laplace
Unit : Modeling in the Frequency Domain Part : Engineering 81: Control Systems I Faculty of Engineering & Applied Science Memorial University of Newfoundland January 1, 010 1 Pair Table Unit, Part : Unit,
More informationHARMONIC FUNCTIONS. x 2 + 2
HARMONIC FUNCTIONS DR. RITU AGARWAL MALAVIYA NATIONAL INSTITUTE OF TECHNOLOGY JAIPUR Contents 1. Harmonic functions 1 1.1. Use of Harmonic mappings 1 1.2. Harmonic functions and holomorphicity 2 1.3. Harmonic
More informationPHYS 3900 Homework Set #03
PHYS 3900 Homework Set #03 Part = HWP 3.0 3.04. Due: Mon. Feb. 2, 208, 4:00pm Part 2 = HWP 3.05, 3.06. Due: Mon. Feb. 9, 208, 4:00pm All textbook problems assigned, unless otherwise stated, are from the
More informationLaplace Theory Examples
Laplace Theory Examples Harmonic oscillator s-differentiation Rule First shifting rule Trigonometric formulas Exponentials Hyperbolic functions s-differentiation Rule First Shifting Rule I and II Damped
More informationDamped harmonic motion
Damped harmonic motion March 3, 016 Harmonic motion is studied in the presence of a damping force proportional to the velocity. The complex method is introduced, and the different cases of under-damping,
More informationA sufficient condition for the existence of the Fourier transform of f : R C is. f(t) dt <. f(t) = 0 otherwise. dt =
Fourier transform Definition.. Let f : R C. F [ft)] = ˆf : R C defined by The Fourier transform of f is the function F [ft)]ω) = ˆfω) := ft)e iωt dt. The inverse Fourier transform of f is the function
More informationf(x)e ikx dx. (1) f(k)e ikx dk. (2) e ikx dx = 2sinak. k f g(x) = g f(x) = f(s)e iks ds
Mathematical Methods: Page 15 2 Fourier transforms 2.1 Integral transforms The Fourier transform is studied in this chapter and the Laplace transform in the next. They are both integral transforms that
More informationLinearization of Differential Equation Models
Linearization of Differential Equation Models 1 Motivation We cannot solve most nonlinear models, so we often instead try to get an overall feel for the way the model behaves: we sometimes talk about looking
More informationSolutions to homework assignment #3 Math 119B UC Davis, Spring = 12x. The Euler-Lagrange equation is. 2q (x) = 12x. q(x) = x 3 + x.
1. Find the stationary points of the following functionals: (a) 1 (q (x) + 1xq(x)) dx, q() =, q(1) = Solution. The functional is 1 L(q, q, x) where L(q, q, x) = q + 1xq. We have q = q, q = 1x. The Euler-Lagrange
More informationThe RC Circuit: An Approach with Fourier Transforms
The RC Circuit: An Approach with Fourier Transforms In this article we shall mathematically analyse the Resistor- Capacitor RC) circuit with the help of Fourier transforms FT). This very general technique
More informationLaplace Transforms Chapter 3
Laplace Transforms Important analytical method for solving linear ordinary differential equations. - Application to nonlinear ODEs? Must linearize first. Laplace transforms play a key role in important
More information221B Lecture Notes on Resonances in Classical Mechanics
1B Lecture Notes on Resonances in Classical Mechanics 1 Harmonic Oscillators Harmonic oscillators appear in many different contexts in classical mechanics. Examples include: spring, pendulum (with a small
More informationPhysics 106a, Caltech 13 November, Lecture 13: Action, Hamilton-Jacobi Theory. Action-Angle Variables
Physics 06a, Caltech 3 November, 08 Lecture 3: Action, Hamilton-Jacobi Theory Starred sections are advanced topics for interest and future reference. The unstarred material will not be tested on the final
More informationThe Principle of Least Action
The Principle of Least Action In their never-ending search for general principles, from which various laws of Physics could be derived, physicists, and most notably theoretical physicists, have often made
More informationMAT389 Fall 2016, Problem Set 11
MAT389 Fall 216, Problem Set 11 Improper integrals 11.1 In each of the following cases, establish the convergence of the given integral and calculate its value. i) x 2 x 2 + 1) 2 ii) x x 2 + 1)x 2 + 2x
More informationEquations of linear stellar oscillations
Chapter 4 Equations of linear stellar oscillations In the present chapter the equations governing small oscillations around a spherical equilibrium state are derived. The general equations were presented
More informationLinear and Nonlinear Oscillators (Lecture 2)
Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical
More informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
More informationBefore you begin read these instructions carefully:
NATURAL SCIENCES TRIPOS Part IB & II (General) Friday, 27 May, 2016 9:00 am to 12:00 pm MATHEMATICS (2) Before you begin read these instructions carefully: You may submit answers to no more than six questions.
More informationFrancis X. Giraldo,
1 Time-Integrators Francis X. Giraldo, giraldo@nrlmry.navy.mil, www.nrlmry.navy.mil/~giraldo/projects/nseam.html 1.1 Introduction Roughly speaking, there are 2 classes of TIs: 1. EulerianMethods(fixed-frame-e.g.,arockatthebottomofaflowing
More informationf (t) K(t, u) d t. f (t) K 1 (t, u) d u. Integral Transform Inverse Fourier Transform
Integral Transforms Massoud Malek An integral transform maps an equation from its original domain into another domain, where it might be manipulated and solved much more easily than in the original domain.
More informationMOMENTS OF HYPERGEOMETRIC HURWITZ ZETA FUNCTIONS
MOMENTS OF HYPERGEOMETRIC HURWITZ ZETA FUNCTIONS ABDUL HASSEN AND HIEU D. NGUYEN Abstract. This paper investigates a generalization the classical Hurwitz zeta function. It is shown that many of the properties
More informationClassical field theory 2012 (NS-364B) Feynman propagator
Classical field theory 212 (NS-364B Feynman propagator 1. Introduction States in quantum mechanics in Schrödinger picture evolve as ( Ψt = Û(t,t Ψt, Û(t,t = T exp ı t dt Ĥ(t, (1 t where Û(t,t denotes the
More informationProblem 1: Lagrangians and Conserved Quantities. Consider the following action for a particle of mass m moving in one dimension
105A Practice Final Solutions March 13, 01 William Kelly Problem 1: Lagrangians and Conserved Quantities Consider the following action for a particle of mass m moving in one dimension S = dtl = mc dt 1
More informationThe spectral transform method
The spectral transform method by Nils Wedi European Centre for Medium-Range Weather Forecasts wedi@ecmwf.int Advanced Numerical Methods for Earth-System Modelling Slide 1 Advanced Numerical Methods for
More informationThe Laplace Transform. Background: Improper Integrals
The Laplace Transform Background: Improper Integrals Recall: Definite Integral: a, b real numbers, a b; f continuous on [a, b] b a f(x) dx 1 Improper integrals: Type I Infinite interval of integration
More informationIntroduction to numerical simulations for Stochastic ODEs
Introduction to numerical simulations for Stochastic ODEs Xingye Kan Illinois Institute of Technology Department of Applied Mathematics Chicago, IL 60616 August 9, 2010 Outline 1 Preliminaries 2 Numerical
More informationProperties of the S-matrix
Properties of the S-matrix In this chapter we specify the kinematics, define the normalisation of amplitudes and cross sections and establish the basic formalism used throughout. All mathematical functions
More informationX(t)e 2πi nt t dt + 1 T
HOMEWORK 31 I) Use the Fourier-Euler formulae to show that, if X(t) is T -periodic function which admits a Fourier series decomposition X(t) = n= c n exp (πi n ) T t, then (1) if X(t) is even c n are all
More informationApplied Mathematics Masters Examination Fall 2016, August 18, 1 4 pm.
Applied Mathematics Masters Examination Fall 16, August 18, 1 4 pm. Each of the fifteen numbered questions is worth points. All questions will be graded, but your score for the examination will be the
More informationChapter 31. The Laplace Transform The Laplace Transform. The Laplace transform of the function f(t) is defined. e st f(t) dt, L[f(t)] =
Chapter 3 The Laplace Transform 3. The Laplace Transform The Laplace transform of the function f(t) is defined L[f(t)] = e st f(t) dt, for all values of s for which the integral exists. The Laplace transform
More information13 Definite integrals
3 Definite integrals Read: Boas h. 4. 3. Laurent series: Def.: Laurent series (LS). If f() is analytic in a region R, then converges in R, with a n = πi f() = a n ( ) n + n= n= f() ; b ( ) n+ n = πi b
More information1 Equal-time and Time-ordered Green Functions
1 Equal-time and Time-ordered Green Functions Predictions for observables in quantum field theories are made by computing expectation values of products of field operators, which are called Green functions
More informationZeta Functions and Regularized Determinants for Elliptic Operators. Elmar Schrohe Institut für Analysis
Zeta Functions and Regularized Determinants for Elliptic Operators Elmar Schrohe Institut für Analysis PDE: The Sound of Drums How Things Started If you heard, in a dark room, two drums playing, a large
More informationTransform Solutions to LTI Systems Part 3
Transform Solutions to LTI Systems Part 3 Example of second order system solution: Same example with increased damping: k=5 N/m, b=6 Ns/m, F=2 N, m=1 Kg Given x(0) = 0, x (0) = 0, find x(t). The revised
More informationComplex Analysis MATH 6300 Fall 2013 Homework 4
Complex Analysis MATH 6300 Fall 2013 Homework 4 Due Wednesday, December 11 at 5 PM Note that to get full credit on any problem in this class, you must solve the problems in an efficient and elegant manner,
More informationNew variables in spherical geometry. David G. Dritschel. Mathematical Institute University of St Andrews.
New variables in spherical geometry David G Dritschel Mathematical Institute University of St Andrews http://www-vortexmcsst-andacuk Collaborators: Ali Mohebalhojeh (Tehran St Andrews) Jemma Shipton &
More informationElectric and Magnetic Forces in Lagrangian and Hamiltonian Formalism
Electric and Magnetic Forces in Lagrangian and Hamiltonian Formalism Benjamin Hornberger 1/26/1 Phy 55, Classical Electrodynamics, Prof. Goldhaber Lecture notes from Oct. 26, 21 Lecture held by Prof. Weisberger
More informationLecture 5: Single Step Methods
Lecture 5: Single Step Methods J.K. Ryan@tudelft.nl WI3097TU Delft Institute of Applied Mathematics Delft University of Technology 1 October 2012 () Single Step Methods 1 October 2012 1 / 44 Outline 1
More informationMORE CONSEQUENCES OF CAUCHY S THEOREM
MOE CONSEQUENCES OF CAUCHY S THEOEM Contents. The Mean Value Property and the Maximum-Modulus Principle 2. Morera s Theorem and some applications 3 3. The Schwarz eflection Principle 6 We have stated Cauchy
More informationThe Laplace Transform. Background: Improper Integrals
The Laplace Transform Background: Improper Integrals Recall: Definite Integral: a, b real numbers, a b; f continuous on [a, b] b a f(x) dx 1 Improper integrals: Type I Infinite interval of integration
More informationBTL What is the value of m if the vector is solenoidal. BTL What is the value of a, b, c if the vector may be irrotational.
VALLIAMMAI ENGINEERING OLLEGE SRM NAGAR, KATTANDKULATHUR Department of Mathematics MA65 - MATHEMATIS II QUESTION BANK - 6 UNIT - I VETOR ALULUS Part - A. Find, if at (, -, ). BTL-. Find the Directional
More informationPHY6426/Fall 07: CLASSICAL MECHANICS HOMEWORK ASSIGNMENT #1 due by 9:35 a.m. Wed 09/05 Instructor: D. L. Maslov Rm.
PHY646/Fall 07: CLASSICAL MECHANICS HOMEWORK ASSIGNMENT # due by 9:35 a.m. Wed 09/05 Instructor: D. L. Maslov maslov@phys.ufl.edu 39-053 Rm. 4 Please help your instructor by doing your work neatly.. Goldstein,
More informationA Filtering Laplace Transform Integration Scheme For Numerical Weather Prediction
A Filtering Laplace Transform Integration Scheme For Numerical Weather Prediction by Colm Clancy B.Sc. A dissertation presented to University College Dublin in partial fulfillment of the requirements for
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Linearization and Characteristic Relations 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
More informationLecture 12. AO Control Theory
Lecture 12 AO Control Theory Claire Max with many thanks to Don Gavel and Don Wiberg UC Santa Cruz February 18, 2016 Page 1 What are control systems? Control is the process of making a system variable
More informationTaylor and Laurent Series
Chapter 4 Taylor and Laurent Series 4.. Taylor Series 4... Taylor Series for Holomorphic Functions. In Real Analysis, the Taylor series of a given function f : R R is given by: f (x + f (x (x x + f (x
More information( ) f (k) = FT (R(x)) = R(k)
Solving ODEs using Fourier Transforms The formulas for derivatives are particularly useful because they reduce ODEs to algebraic expressions. Consider the following ODE d 2 dx + p d 2 dx + q f (x) = R(x)
More informationODE Homework 1. Due Wed. 19 August 2009; At the beginning of the class
ODE Homework Due Wed. 9 August 2009; At the beginning of the class. (a) Solve Lẏ + Ry = E sin(ωt) with y(0) = k () L, R, E, ω are positive constants. (b) What is the limit of the solution as ω 0? (c) Is
More informationMath Fall Linear Filters
Math 658-6 Fall 212 Linear Filters 1. Convolutions and filters. A filter is a black box that takes an input signal, processes it, and then returns an output signal that in some way modifies the input.
More informationA Note on the Differential Equations with Distributional Coefficients
MATEMATIKA, 24, Jilid 2, Bil. 2, hlm. 115 124 c Jabatan Matematik, UTM. A Note on the Differential Equations with Distributional Coefficients Adem Kilicman Department of Mathematics, Institute for Mathematical
More informationThe 3-wave PDEs for resonantly interac7ng triads
The 3-wave PDEs for resonantly interac7ng triads Ruth Mar7n & Harvey Segur Waves and singulari7es in incompressible fluids ICERM, April 28, 2017 What are the 3-wave equa.ons? What is a resonant triad?
More informationKrylov-Subspace Based Model Reduction of Nonlinear Circuit Models Using Bilinear and Quadratic-Linear Approximations
Krylov-Subspace Based Model Reduction of Nonlinear Circuit Models Using Bilinear and Quadratic-Linear Approximations Peter Benner and Tobias Breiten Abstract We discuss Krylov-subspace based model reduction
More informationA Lecture on Selective RF-pulses in MRI
A Lecture on Selective RF-pulses in MRI Charles L. Epstein February, 23 Introduction We describe the notion of selective excitation and explain how the Bloch equation is used to design selective RF-pulses.
More informationChapter 1. Harmonic Oscillator. 1.1 Energy Analysis
Chapter 1 Harmonic Oscillator Figure 1.1 illustrates the prototypical harmonic oscillator, the mass-spring system. A mass is attached to one end of a spring. The other end of the spring is attached to
More informationSolutions to practice problems for the final
Solutions to practice problems for the final Holomorphicity, Cauchy-Riemann equations, and Cauchy-Goursat theorem 1. (a) Show that there is a holomorphic function on Ω = {z z > 2} whose derivative is z
More information22.2. Applications of Eigenvalues and Eigenvectors. Introduction. Prerequisites. Learning Outcomes
Applications of Eigenvalues and Eigenvectors 22.2 Introduction Many applications of matrices in both engineering and science utilize eigenvalues and, sometimes, eigenvectors. Control theory, vibration
More informationDSP-I DSP-I DSP-I DSP-I
DSP-I DSP-I DSP-I DSP-I Digital Signal Processing I (8-79) Fall Semester, 005 OTES FOR 8-79 LECTURE 9: PROPERTIES AD EXAPLES OF Z-TRASFORS Distributed: September 7, 005 otes: This handout contains in outline
More informationChapter 2. The continuous equations
Chapter. The continuous equations Fig. 1.: Schematic of a forecast with slowly varying weather-related variations and superimposed high frequency Lamb waves. Note that even though the forecast of the slow
More informationResponse to Periodic and Non-periodic Loadings. Giacomo Boffi.
Periodic and Non-periodic http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano March 7, 218 Outline Undamped SDOF systems
More informationComplex Numbers. The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition,
Complex Numbers Complex Algebra The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition, and (ii) complex multiplication, (x 1, y 1 )
More informationω (rad/s)
1. (a) From the figure we see that the signal has energy content in frequencies up to about 15rad/s. According to the sampling theorem, we must therefore sample with at least twice that frequency: 3rad/s
More informationThe Laplace Transform
C H A P T E R 6 The Laplace Transform Many practical engineering problems involve mechanical or electrical systems acted on by discontinuous or impulsive forcing terms. For such problems the methods described
More informationP321(b), Assignement 1
P31(b), Assignement 1 1 Exercise 3.1 (Fetter and Walecka) a) The problem is that of a point mass rotating along a circle of radius a, rotating with a constant angular velocity Ω. Generally, 3 coordinates
More information16. Local theory of regular singular points and applications
16. Local theory of regular singular points and applications 265 16. Local theory of regular singular points and applications In this section we consider linear systems defined by the germs of meromorphic
More informationAnalytical Mechanics ( AM )
Analytical Mechanics ( AM ) Olaf Scholten KVI, kamer v8; tel nr 6-55; email: scholten@kvinl Web page: http://wwwkvinl/ scholten Book: Classical Dynamics of Particles and Systems, Stephen T Thornton & Jerry
More informationAspects of Continuous- and Discrete-Time Signals and Systems
Aspects of Continuous- and Discrete-Time Signals and Systems C.S. Ramalingam Department of Electrical Engineering IIT Madras C.S. Ramalingam (EE Dept., IIT Madras) Networks and Systems 1 / 45 Scaling the
More informationHeisenberg's inequality for Fourier transform
Heisenberg's inequality for Fourier transform Riccardo Pascuzzo Abstract In this paper, we prove the Heisenberg's inequality using the Fourier transform. Then we show that the equality holds for the Gaussian
More informationLecture 4: Birkhoff normal forms
Lecture 4: Birkhoff normal forms Walter Craig Department of Mathematics & Statistics Waves in Flows - Lecture 4 Prague Summer School 2018 Czech Academy of Science August 31 2018 Outline Two ODEs Water
More informationScattered Interpolation Survey
Scattered Interpolation Survey j.p.lewis u. southern california Scattered Interpolation Survey p.1/53 Scattered vs. Regular domain Scattered Interpolation Survey p.2/53 Motivation modeling animated character
More information1 Constant Pressure Inner Boundary Condition
Jim Lambers ENERGY 281 Spring Quarter 2007-08 Lecture 8 Notes These notes are based on Rosalind Archer s PE281 lecture notes, with some revisions by Jim Lambers. 1 Constant Pressure Inner Boundary Condition
More informationComplex Variables. Instructions Solve any eight of the following ten problems. Explain your reasoning in complete sentences to maximize credit.
Instructions Solve any eight of the following ten problems. Explain your reasoning in complete sentences to maximize credit. 1. The TI-89 calculator says, reasonably enough, that x 1) 1/3 1 ) 3 = 8. lim
More informationSome Fun with Divergent Series
Some Fun with Divergent Series 1. Preliminary Results We begin by examining the (divergent) infinite series S 1 = 1 + 2 + 3 + 4 + 5 + 6 + = k=1 k S 2 = 1 2 + 2 2 + 3 2 + 4 2 + 5 2 + 6 2 + = k=1 k 2 (i)
More informationCausality. but that does not mean it is local in time, for = 1. Let us write ɛ(ω) = ɛ 0 [1 + χ e (ω)] in terms of the electric susceptibility.
We have seen that the issue of how ɛ, µ n depend on ω raises questions about causality: Can signals travel faster than c, or even backwards in time? It is very often useful to assume that polarization
More informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
More informationThe Hopf equation. The Hopf equation A toy model of fluid mechanics
The Hopf equation A toy model of fluid mechanics 1. Main physical features Mathematical description of a continuous medium At the microscopic level, a fluid is a collection of interacting particles (Van
More information2 Write down the range of values of α (real) or β (complex) for which the following integrals converge. (i) e z2 dz where {γ : z = se iα, < s < }
Mathematical Tripos Part II Michaelmas term 2007 Further Complex Methods, Examples sheet Dr S.T.C. Siklos Comments and corrections: e-mail to stcs@cam. Sheet with commentary available for supervisors.
More informationPATH INTEGRAL SOLUTION FOR A PARTICLE WITH POSITION DEPENDENT MASS
Vol. 9 ACTA PHYSICA POLONICA B No PATH INTEGRAL SOLUTION FOR A PARTICLE WITH POSITION DEPENDENT MASS N. Bouchemla, L. Chetouani Département de Physique, Faculté de Sciences Exactes Université Mentouri,
More informationExercises for Part 1
MATH200 Complex Analysis. Exercises for Part Exercises for Part The following exercises are provided for you to revise complex numbers. Exercise. Write the following expressions in the form x + iy, x,y
More informationSTABILITY. Phase portraits and local stability
MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),
More information27. Euler s Equations
27 Euler s Equations Michael Fowler Introduction We ve just seen that by specifying the rotational direction and the angular phase of a rotating body using Euler s angles, we can write the Lagrangian in
More informationBalance in the Atmosphere: Implications for Numerical Weather Prediction
Balance in the Atmosphere: Implications for Numerical Weather Prediction Peter Lynch School of Mathematical Sciences University College Dublin BIRS Summer School, Banff, 10 15 July, 2011 Outline Introduction
More informationMATH 311: COMPLEX ANALYSIS CONTOUR INTEGRALS LECTURE
MATH 3: COMPLEX ANALYSIS CONTOUR INTEGRALS LECTURE Recall the Residue Theorem: Let be a simple closed loop, traversed counterclockwise. Let f be a function that is analytic on and meromorphic inside. Then
More informationControl Systems. Frequency domain analysis. L. Lanari
Control Systems m i l e r p r a in r e v y n is o Frequency domain analysis L. Lanari outline introduce the Laplace unilateral transform define its properties show its advantages in turning ODEs to algebraic
More informationLecture 9 February 2, 2016
MATH 262/CME 372: Applied Fourier Analysis and Winter 26 Elements of Modern Signal Processing Lecture 9 February 2, 26 Prof. Emmanuel Candes Scribe: Carlos A. Sing-Long, Edited by E. Bates Outline Agenda:
More informationPoincaré Map, Floquet Theory, and Stability of Periodic Orbits
Poincaré Map, Floquet Theory, and Stability of Periodic Orbits CDS140A Lecturer: W.S. Koon Fall, 2006 1 Poincaré Maps Definition (Poincaré Map): Consider ẋ = f(x) with periodic solution x(t). Construct
More informationChapter 2: Complex numbers
Chapter 2: Complex numbers Complex numbers are commonplace in physics and engineering. In particular, complex numbers enable us to simplify equations and/or more easily find solutions to equations. We
More information10. Zwanzig-Mori Formalism
University of Rhode Island DigitalCommons@URI Nonequilibrium Statistical Physics Physics Course Materials 0-9-205 0. Zwanzig-Mori Formalism Gerhard Müller University of Rhode Island, gmuller@uri.edu Follow
More informationSpiral Wave Chimeras
Spiral Wave Chimeras Carlo R. Laing IIMS, Massey University, Auckland, New Zealand Erik Martens MPI of Dynamics and Self-Organization, Göttingen, Germany Steve Strogatz Dept of Maths, Cornell University,
More informationGFD 2012 Lecture 1: Dynamics of Coherent Structures and their Impact on Transport and Predictability
GFD 2012 Lecture 1: Dynamics of Coherent Structures and their Impact on Transport and Predictability Jeffrey B. Weiss; notes by Duncan Hewitt and Pedram Hassanzadeh June 18, 2012 1 Introduction 1.1 What
More informationλ n = L φ n = π L eınπx/l, for n Z
Chapter 32 The Fourier Transform 32. Derivation from a Fourier Series Consider the eigenvalue problem y + λy =, y( L = y(l, y ( L = y (L. The eigenvalues and eigenfunctions are ( nπ λ n = L 2 for n Z +
More informationPhysics 250 Green s functions for ordinary differential equations
Physics 25 Green s functions for ordinary differential equations Peter Young November 25, 27 Homogeneous Equations We have already discussed second order linear homogeneous differential equations, which
More information