Part I: Numerical Fourier analysis of quasi periodic functions
|
|
- Eustacia Marsh
- 5 years ago
- Views:
Transcription
1 Part I: Numerical Fourier analysis of quasi periodic functions I.1 The procedure I.2 Error estimation f(t) Q f (t; A, ω) DFT(f) = DFT(Q f ) Applying the above procedure, the system to be solved can be written as g(y + y) = b + b. The exact frequencies and amplitudes would be obtained if b = 0. A bound for the error is given by y Dg(y) 1 b. We will bound Dg(y) 1 and b. I.3 Applications Academic example, to show the procedure and the accuracy of the error estimation Solar System models: SS = RTBP + Perturbations (Fourier analysis) SS = RTBP + Perturbations(ω 1, ω 2,...) 3
2 Part II: The neighborhood of the collinear equilibrium points in the RTBP Final goal: To get Poincaré section representations of the flow around L 1,2,3,. 0.3 Energy Tools: Parallel algorithms for the computation of different kinds of families of periodic orbits and invariant tori. 4
3 I.1 Refined Fourier analysis procedure Given a real analytic quasi periodic function, f(t), at N equally spaced points on an interval [0, T], {f(t l )} N 1 l=0, t l = l T N, the goal is to compute a trigonometric approximation N f Q f (t) = A c 0 + (A c l cos(2πν l T t) + As l sin(2πν l T t). The Discrete Fourier Transform (DFT) of f is defined through f(t l ) = 1 2 (c f,t,n(0) + c f,t,n ( N 2 )cos(2πn 2 t l T ))+ N/2 1 (c f,t,n (j)cos( 2πjt l T ) + s f,t,n(j)sin( 2πjt l T )). where j=1 c f,t,n (j) = 2 N s f,t,n (j) = 2 N N 1 l=0 N 1 l=0 f(t l )cos(2π j N l), j = 0,..., N 2, f(t l )sin(2π j N l), j = 1,..., N
4 Algorithm 1. Set an starting thresold for collecting peaks of the modulus of the DFT of f(t). 2. Find initial approximations of the frequencies, starting from the peaks of the DFT greater than the thresold. 3. Find the amplitudes of the frequencies found in the previous step (DFT(Q f ) = DFT(f)). 4. Simultaneously refine ALL the frequencies and amplitudes of the current quasi periodic approximation of f. 5. Perform a DFT of the input signal minus the current quasi periodic approximation obtained in step 4, decrease the thresold and go back to step 2. 6-a
5 Equations of step 3 We ask DFT(Q f ) = DFT(f), being N f Q f (t) = A c 0 + (A c l cos(2πν l T t) + As l sin(2πν l T t). The system of equations to be solved is linear and (1 + 2N f ) (1 + 2N f ): Nf A c 0 cn h 1,T,N (0) + (A c l cn h ν l,n (0) + As l c n h ν l,n (0)) = cn h Nf A c 0 cn h 1,T,N (j) + (A c l cn h ν l,n (j) + As l c n h ν l,n (j)) = cn h N f (A c l sn h ν l,t (j) + Ac l s n h ν l,t (j)) = sn h where j = [ν l + 0.5], l = 1 N f, and c n h 1 (j) = cn h 1,T,N (j), c n h ν l,n (j) = cn h cos( 2πν l ),T,N(j), sn h ν l,n (j) = sn h T c n h ν l,n (j) = cn h sin( 2πν l ),T,N(j), sn h ν l,n (j) = sn h T cos( 2πν l T sin( 2πν l T f,t,n (0) f,t,n (j) f,t,n (j) ),T,N(j), ),T,N(j). 6-b
6 Equations for step 4 The system of equations is a (1 + 3N f ) (1 + 3N f ) non linear system which is solved by Newton s method: A c 0 cn h 1,T,N (0) + Nf A c 0 cn h 1,T,N (j i) + A c 0 csn h 1,T,N (j+ i ) + N f N f N f (A c l cn h ν l,n (0) + As l c n h ν l,n (0)) = cn h (A c l cn h ν l,n (j i) + A s l c n h ν l,n (j i)) = c n h (A c l sn h ν l,n (j i) + A s l s n h ν l,n (j i)) = s n h (A c l csn h ν l,n (j+ i ) + A s l cs n h ν l,n (j+ i )) = cs n h f,t,n (0) f,t,n (j i) f,t,n (j i) f,t,n (j+ i ) being j i = [ν i + 0.5], j + i j i, j + i j i = 1. 6-c
7 An example: f(t) = cos(2π0.23t) 1 sin(2π0.27t) + sin(2π0.37t) Starting thresold: 0.8 modulus of the DFT of the input data: peaks j = 61, j = Approximation of frequencies (Laskar s method): peak 61 frequency peak 189 frequency Computation of amplitudes from known frequencies: Frequency Cosine amplitude Sine amplitude modulus of the DFT of the residual Iterative refinement: Frequency Cosine amplitude Sine amplitude a
8 5. modulus of the DFT of input signal minus step 4: New threshold: Approximation of frequencies (Laskar s method): peak 138 frequency Amplitudes from known frequencies: Frequency Cosine amplitude Sine amplitude Iterative refinement: Frequency Cosine amplitude Sine amplitude modulus of the DFT of the residual: 5e-14 4e-14 3e-14 2e-14 1e b
9 I.2 Error estimation We assume f is real analytic and quasi periodic, f(t) = A c 0 + k Z m kω>0 (A c k cos(2πkωt) + As k sin(2πkωt)), and its Fourier coefficients a k satisfy the Cauchy estimates (A c k )2 + (A s k )2 Ce δ k k Z m. The frequency vector ω = (ω 1,... ω m ) satisfies a Diophantine condition with D, τ > 0. kω D k τ, We determine the frequencies {ν l } N f of order r 0 1 (ν l Tkω, 1 k r 0 1) and the corresponding amplitudes {A c l }N f l=0, {As l }N f. 8
10 Strategy Let us denote, f r0 (t) = A c 0 + (A c k cos(2πkωt) + As k sin(2πkωt)). k r 0 1 kω>0 The system of equations used for simultaneous improvement of frequencies and amplitudes can be written as A c 0 + Nf A c 0 cn h 1 (j i) + A c 0 csn h 1 (j+ i ) + N f N f N f (A c l cn h ν l,n (0) + As l c n h ν l,n (0)) = cn h f r0,t,n (0) + cn h f f r0,t,n (0) (A c l cn h ν l,n (j i) + A s l c n h ν l,n (j i)) = c n h f r0,t,n (j i) + c n h f f r0,t,n (j i) (A c l sn h ν l,n (j i) + A s l s n h ν l,n (j i)) = s n h f r0,t,n (j i) + s n h f f r0,t,n (j i) (A c l csn h ν l,n (j+ i ) + A s l cs n h ν l,n (j+ i )) }{{} g(y + y) = cs n h f r0,t,n (j+ i ) + cs n h f f r0,t,n (j+ i ). } {{ } b } {{ } b We would get the exact frequencies and amplitudes if b = 0. The error in frequencies and amplitudes is given, in the first order approximation, by y Dg(y) 1 b. We will obtain bounds for Dg(y) 1 and b. 9
11 Bound for Dg(y) 1 Preliminaries We can write Dg(y) =: M = 2 B 0,1... B 0,Nf 0 B 1,1... B 1,Nf B Nf,1... B Nf,N f. We split M = M D + M O, M = 0 B 1, B Nf,N f + 0 B 0,1... B 0,Nf B 1,Nf B Nf, The idea is to obtain bounds for M 1 D, M O and use (M D + M O ) 1 M 1 D 1 M 1 D M O. The components B i,l of M (i 0) have the following general aspect A c l cn h ν l,n (j i) + A s l cn h ν l,n (j i) c n h ν l,n (j i) c n h ν l,n (j i) A c l sn h ν l,n (j i) + A s l sn h ν l,n (j i) s n h ν l,n (j i) A c l csn h ν l,n (j+ i ) + A s l csn h ν l,n (j+ i ) cs n h ν l,n (j+ i ) cs n h s n h ν l,n (j i) ν l,n (j+ i ). 10
12 First simplification: M M The first simplification consists in the use of the Truncated Continuous Fourier Tansform (C, S), 1 T T 0 H n h T (t)f(t)e i2π j T t =: 1 2 (Cn h f,t (j) + isn h f,t (j)) instead of the DFT (c, s). We define CS n h ν (j) as in the discrete case. In this way, M M = M D + M O (substitute DFT by TCFT) = M D + M O Remark 1: The difference between cs and CS can be bounded from a discrete version of Poisson s summation formula: c n h f,t,n (j) = l= C n + ln h f,t (j ),... T Remark 2: Explicit formulae for the components of M can be obtained. 11
13 Formulae for the components of M The components of the matrix M are given by ( C n h ν (j) = K sin(2π(ν j)) n h ψ nh (ν j).. ( C n h ν (j) = K h r (ν j) n h ψ nh (ν j).. + sin(2π( ν j)) ψ nh ( ν j). h r ( ν j) ψ nh ( ν j). ), ), where K nh = ( 1)n h (n h!) 2 n h 2π ψ nh (x) = (x + l) l= n h h r (x) = 2π cos(2πx) r nh (x)sin(2πx), h i (x) = 2π sin(2πx) r nh (x)(1 cos(2πx)), r nh (x) = n h l= n h 1 x + l = ψ n h (x) ψ nh (x). 12
14 Second simplification: M M The second simplification consist in eliminating the second summands of the expressions of the components of M as given in the previous lemma. In this way, M = M D + M O (remove the second summands) M = M D + M O Remark 1: The difference between the components of M and M ( CS and cs) can be bounded. Remark 2: the expression for cs only depends on the difference ν j. 13
15 Bound for M 1 D Since M D is diagonal, to compute M 1 D to invert its diagonal blocks B i,i A c i cn h ν i,n (j i) + A s i cn h ν i,n (j i) c n h ν i,n (j i) A c i sn h ν i,n (j i) + A s i sn h ν i,n (j i) s n h ν i,n (j i) A c i csn h ν i,n (j+ i ) + A s i csn h ν i,n (j+ i ) cs n h ν i,n (j+ i ) cs n h we only need c n h ν i,n (j i) s n h ν i,n (j i) ν i,n (j+ i ) For that, we first show that, if (A c i, As i ) (0,0), B i,i is invertible either setting cs = c or cs = s.. The actual bounds of B i,i are computed numerically for each n h, by first minimizing w.r.t. cs = c, s and maximizing w.r.t. θ [0,2π], ν j 1/2 the supremum norm of c n h ν i,n (j i)cos θ + c n h ν i,n (j i)sin θ ν i,n (j i)sin θ c n h s n h ν i,n (j i) ν i,n (j i) s n h ν i,n (j i)cos θ + s n h cs n h ν i,n (j+ i )cos θ + cs n h ν i,n(j + i )sin θ cs n h ν i,n (j+ i ) cs n h We call G nh In this way, we get being c n h ν i,n (j i) s n h ν i,n (j i) ν i,n(j + i ) the obtained quantity. Some values are n h G nh (M D ) 1 max(a 1 min,1)g n h, A i = ((A c i) 2 + (A s i) 2 ) 1/2 A min = min i=1 N f {A 1,..., A Nf } 1. 14
16 Bounds for M 1 The computation of the bound for Dg(y) 1 is completed by following the following scheme: M 1 D max(a 1 min,1)g n h M 1 D M 1 D 1 M 1 D M D M D M 1 D M 1 D 1 M 1 D M D M D M O M O + M O M O + M O M O M 1 M 1 D 1 M 1 D M O 15
17 We have b 2C max j J Bound for b k =r 0 e δ k hnh N (Tkω j), where hn h N is the envelope displayed below (N = 16, n h = 0) The Diophantine condition gives a lower bound for Tkω j : Tkω j TD ( k + k j ) τ 1. For k small, hn h N (Tkω j) 1. After some order r, hn h N (Tkω j) may approach a
18 Therefore, b 2C(max j J where: r 1 k =r 0 e δ k hnh N (Tkω j) +max j J k =r e δ k ), The first term is bounded by replacing the DFT by the TCFT. This introduces an additional error term due to this approximation. All the sums are reduced to sums of the form j jα e δj, which are bounded by incomplete Gamma functions. 16-b
19 Final theorem Under the stated hypothesis, the error in frequencies and amplitudes can be bounded as where and M O (n h!) and M 1 D being ε 1 = ε 2 = π 2( N f 2( N f ( N ( 2( f 4 y M 1 b, M 1 TD M 1 D 1 M 1 D M O A l )(π + ln( 1 + n (2r 0 2) τ h ) ln( 2 n (2r 0 2) τ h )) + 2N f TD ( 1 n (2r 0 2) τ h ) 1+2n h A l )(π + ln([ν min ] + n h ) ln([ν min ] 1 n h )) + 2N f ([ν min ] n h ) 1+2n h A l )(π+ln(n Ω 0 +n h ) ln(n Ω 0 1 n h ))+2N f )(1+ 1 M 1 D 1 M 1 ε, M 1 D D 1 (N Ω 0 n h ) 1+2n h TD M 1 D 1 M 1 ε, M 1 D D 2 2n h ) G nh min(1, A min ), ( 2Amax ) 4(n h!) 2 (π+ln(n Ω 0 +n h ) ln(n Ω 0 1 n h ))+2 (1+ 1 π(n Ω 0 n h ) 1+2n h (n h!) 2 ( 2Amax (π+ln(2[ν min ]+n h ) ln(2[ν min ] n h 1)) + 2 Ω 0 = T(2r 0 2) ω +1, π(2[ν min ] n h ) 1+2n h ), 17 ) 2n h ),
20 As for b, ( b χ {r >r0 } 2m+1 C (m 1)! (n h!) 2 e δ(r 0 1) + χ {r >r0 } Error estimates: Final theorem (cont.) m 1 ( m 1 l l=0 ) ( m 2 r 0+1) m 1 l G f (2r 0 1, r 0 +r 2, l+τ(1+2n h ), δ) E π(td) 1+2n h 2(n h!) 2 e δ m 2 (1 + 1 )G 2n f (r 0 + m, r h 2 1+ m, m 1, δ) 2 π(n Ω n h ) 1+2n h + e δ m 2 G (r + m 2, m 1, δ) ), where Ω = T(r + r 0 2) ω + 1 r = max ( r 0,min E = (z 1 n h ) 1+2n h, z = z 1+2n h TD (r + r 0 2) τ, ([( TD max(( (n h!)2 ) 1 1+2n h n π h,2(1 + n h )) [ N 1 nh ]) ) r T ω ) 1 τ r0 + 2], equals 1 if condi- In the above formulas, χ {condition} tion is true and 0 otherwise. 18
21 I.3 Applications Academic example We consider the quasi periodic function f 0.9 (t) = sin(2πω 1 t + ϕ 1 ) sin(2πω 2 t + ϕ 2 ) 1 0.9cos(2πω 1 t + ϕ 1 ) 1 0.9cos(2πω 2 t + ϕ 2 ). Explicit formulae for frequencies and amplitudes can be obtained, as well as the Cauchy estimates and the Diophantine condition. We have performed Fourier analysis of this function for several T, N. µ=0.9 log10(error) log2(t) log2(t/n) µ= log10(error) log2(t/n) Solid line: actual error. Dashed line: estimated error. 19
22 Numerical test of the bounds obtained: actual error vs. predicted error In the preceeding slide, the difference between the error predicted and the actual error is big because of the Diophantine condition, which is reached for for very few orders k. The points in the plot below represent min k =const. kω for k = The curve represents the values of the Diophantine condition / k min kω If we substitute the first term of b by what is obtained just before the use of the Diophantine condition, the following figure is obtained. k log10(error) µ= log2(t/n) 13 log2(t) log10(error) µ= log2(t/n) Solid line: actual error. Dashed line: estimated error. 20
Lecture 2: Modeling Accelerators Calculation of lattice functions and parameters. X. Huang USPAS, January 2015 Hampton, Virginia
Lecture 2: Modeling Accelerators Calculation of lattice functions and parameters X. Huang USPAS, January 2015 Hampton, Virginia 1 Outline Closed orbit Transfer matrix, tunes, Optics functions Chromatic
More informationPart II. The neighborhood of the collinear equilibrium points in the RTBP
Part II The neighborhood of the collinear equilibrium points in the RTBP Part II 9 This part is devoted to the study of the vicinity of the collinear equilibrium points of the RTBP for the Earth Moon
More information1 Orthonormal sets in Hilbert space
Math 857 Fall 15 1 Orthonormal sets in Hilbert space Let S H. We denote by [S] the span of S, i.e., the set of all linear combinations of elements from S. A set {u α : α A} is called orthonormal, if u
More informationAssignment 3 Solutions
Assignment Solutions Networks and systems August 8, 7. Consider an LTI system with transfer function H(jw) = input is sin(t + π 4 ), what is the output? +jw. If the Solution : C For an LTI system with
More information4. Sinusoidal solutions
16 4. Sinusoidal solutions Many things in nature are periodic, even sinusoidal. We will begin by reviewing terms surrounding periodic functions. If an LTI system is fed a periodic input signal, we have
More informationECE 301 Fall 2011 Division 1. Homework 1 Solutions.
ECE 3 Fall 2 Division. Homework Solutions. Reading: Course information handout on the course website; textbook sections.,.,.2,.3,.4; online review notes on complex numbers. Problem. For each discrete-time
More informationDistributed detection of topological changes in communication networks. Riccardo Lucchese, Damiano Varagnolo, Karl H. Johansson
1 Distributed detection of topological changes in communication networks Riccardo Lucchese, Damiano Varagnolo, Karl H. Johansson Thanks to... 2 The need: detecting changes in topological networks 3 The
More informationMathematical methods and its applications Dr. S. K. Gupta Department of Mathematics Indian Institute of Technology, Roorkee
Mathematical methods and its applications Dr. S. K. Gupta Department of Mathematics Indian Institute of Technology, Roorkee Lecture - 56 Fourier sine and cosine transforms Welcome to lecture series on
More informationLecture 3f Polar Form (pages )
Lecture 3f Polar Form (pages 399-402) In the previous lecture, we saw that we can visualize a complex number as a point in the complex plane. This turns out to be remarkable useful, but we need to think
More informationHomework #2 Due Monday, April 18, 2012
12.540 Homework #2 Due Monday, April 18, 2012 Matlab solution codes are given in HW02_2012.m This code uses cells and contains the solutions to all the questions. Question 1: Non-linear estimation problem
More informationMathematical Methods: Fourier Series. Fourier Series: The Basics
1 Mathematical Methods: Fourier Series Fourier Series: The Basics Fourier series are a method of representing periodic functions. It is a very useful and powerful tool in many situations. It is sufficiently
More informationImproving Boundary Conditions in Time- Domain Self-Force Calculations
Improving Boundary Conditions in Time- Domain Self-Force Calculations Carlos F. Sopuerta Institute of Space Sciences National Spanish Research Council Work in Collaboration with Anil Zenginoğlu (Caltech)
More informationCotangent and the Herglotz trick
Cotangent and the Herglotz trick Yang Han Christian Wude May 0, 0 The following script introduces the partial fraction expression of the cotangent function and provides an elegant proof, using the Herglotz
More informationMath 251 December 14, 2005 Final Exam. 1 18pt 2 16pt 3 12pt 4 14pt 5 12pt 6 14pt 7 14pt 8 16pt 9 20pt 10 14pt Total 150pt
Math 251 December 14, 2005 Final Exam Name Section There are 10 questions on this exam. Many of them have multiple parts. The point value of each question is indicated either at the beginning of each question
More informationPeriodic orbits high above the ecliptic plane in the solar sail 3-body problem
Periodic orbits high above the ecliptic plane in the solar sail 3-body problem Thomas Waters Department of Mechanical Engineering, University of Strathclyde, Glasgow. 4th October 2006 In conjunction with
More informationJUHA KINNUNEN. Partial Differential Equations
JUHA KINNUNEN Partial Differential Equations Department of Mathematics and Systems Analysis, Aalto University 207 Contents INTRODUCTION 2 FOURIER SERIES AND PDES 5 2. Periodic functions*.............................
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 informationFirst Year Examination Department of Statistics, University of Florida
First Year Examination Department of Statistics, University of Florida August 20, 2009, 8:00 am - 2:00 noon Instructions:. You have four hours to answer questions in this examination. 2. You must show
More informationSignal Functions (0B)
Signal Functions (0B) Signal Functions Copyright (c) 2009-203 Young W. Lim. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License,
More informationReview of Fundamentals of Digital Signal Processing
Chapter 2 Review of Fundamentals of Digital Signal Processing 2.1 (a) This system is not linear (the constant term makes it non linear) but is shift-invariant (b) This system is linear but not shift-invariant
More informationA fast quantum mechanical algorithm for estimating the median
A fast quantum mechanical algorithm for estimating the median Lov K Grover 3C-404A Bell Labs 600 Mountain Avenue Murray Hill J 07974 lkg@mhcnetattcom Summary Consider the problem of estimating the median
More informationTHE METHOD OF SEPARATION OF VARIABLES
THE METHOD OF SEPARATION OF VARIABES To solve the BVPs that we have encountered so far, we will use separation of variables on the homogeneous part of the BVP. This separation of variables leads to problems
More informationON THE BOUNDEDNESS BEHAVIOR OF THE SPECTRAL FACTORIZATION IN THE WIENER ALGEBRA FOR FIR DATA
ON THE BOUNDEDNESS BEHAVIOR OF THE SPECTRAL FACTORIZATION IN THE WIENER ALGEBRA FOR FIR DATA Holger Boche and Volker Pohl Technische Universität Berlin, Heinrich Hertz Chair for Mobile Communications Werner-von-Siemens
More informationEE 3054: Signals, Systems, and Transforms Summer It is observed of some continuous-time LTI system that the input signal.
EE 34: Signals, Systems, and Transforms Summer 7 Test No notes, closed book. Show your work. Simplify your answers. 3. It is observed of some continuous-time LTI system that the input signal = 3 u(t) produces
More informationEllipticTheta2. Notations. Primary definition. Specific values. Traditional name. Traditional notation. Mathematica StandardForm notation
EllipticTheta Notations Traditional name Jacobi theta function ϑ Traditional notation ϑ z, Mathematica StandardForm notation EllipticTheta, z, Primary definition 09.0.0.0001.01 ϑ z, cos k 1 z ; 1 Specific
More information2.161 Signal Processing: Continuous and Discrete
MIT OpenCourseWare http://ocw.mit.edu.6 Signal Processing: Continuous and Discrete Fall 008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. M MASSACHUSETTS
More informationLecture 34. Fourier Transforms
Lecture 34 Fourier Transforms In this section, we introduce the Fourier transform, a method of analyzing the frequency content of functions that are no longer τ-periodic, but which are defined over the
More informationIntroduction to Time Series Analysis. Lecture 18.
Introduction to Time Series Analysis. Lecture 18. 1. Review: Spectral density, rational spectra, linear filters. 2. Frequency response of linear filters. 3. Spectral estimation 4. Sample autocovariance
More informationTowards stability results for planetary problems with more than three bodies
Towards stability results for planetary problems with more than three bodies Ugo Locatelli [a] and Marco Sansottera [b] [a] Math. Dep. of Università degli Studi di Roma Tor Vergata [b] Math. Dep. of Università
More informationIntroduction to Fourier Analysis Part 2. CS 510 Lecture #7 January 31, 2018
Introduction to Fourier Analysis Part 2 CS 510 Lecture #7 January 31, 2018 OpenCV on CS Dept. Machines 2/4/18 CSU CS 510, Ross Beveridge & Bruce Draper 2 In the extreme, a square wave Graphic from http://www.mechatronics.colostate.edu/figures/4-4.jpg
More informationA LOCALIZATION PROPERTY AT THE BOUNDARY FOR MONGE-AMPERE EQUATION
A LOCALIZATION PROPERTY AT THE BOUNDARY FOR MONGE-AMPERE EQUATION O. SAVIN. Introduction In this paper we study the geometry of the sections for solutions to the Monge- Ampere equation det D 2 u = f, u
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 informationIntroduction to Fourier Analysis. CS 510 Lecture #5 February 2 nd 2015
Introduction to Fourier Analysis CS 510 Lecture 5 February 2 nd 2015 Why? Get used to changing representa6ons! (and this par6cular transforma6on is ubiquitous and important) 2/9/15 CSU CS 510, Ross Beveridge
More informationData Structures and Algorithms CMPSC 465
Data Structures and Algorithms CMPSC 465 LECTURE 9 Solving recurrences Substitution method Adam Smith S. Raskhodnikova and A. Smith; based on slides by E. Demaine and C. Leiserson Review question Draw
More informationElectrostatic Imaging via Conformal Mapping. R. Kress. joint work with I. Akduman, Istanbul and H. Haddar, Paris
Electrostatic Imaging via Conformal Mapping R. Kress Göttingen joint work with I. Akduman, Istanbul and H. Haddar, Paris Or: A new solution method for inverse boundary value problems for the Laplace equation
More informationImage Analysis. 3. Fourier Transform
Image Analysis Image Analysis 3. Fourier Transform Lars Schmidt-Thieme Information Systems and Machine Learning Lab (ISMLL) Institute for Business Economics and Information Systems & Institute for Computer
More informationMath 216 Second Midterm 19 March, 2018
Math 26 Second Midterm 9 March, 28 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material that
More informationExaminer: D. Burbulla. Aids permitted: Formula Sheet, and Casio FX-991 or Sharp EL-520 calculator.
University of Toronto Faculty of Applied Science and Engineering Solutions to Final Examination, June 017 Duration: and 1/ hrs First Year - CHE, CIV, CPE, ELE, ENG, IND, LME, MEC, MMS MAT187H1F - Calculus
More informationTravelling and Standing Waves
Travelling and Standing Waves Many biological phenomena are cyclic. Heartbeats Circadian rhythms, Estrus cycles Many more e.g. s Such events are best described as waves. Therefore the study of waves is
More informationNotes 14 The Steepest-Descent Method
ECE 638 Fall 17 David R. Jackson Notes 14 The Steepest-Descent Method Notes are adapted from ECE 6341 1 Steepest-Descent Method Complex Integral: Ω g z I f z e dz Ω = C The method was published by Peter
More informationWe start with a simple result from Fourier analysis. Given a function f : [0, 1] C, we define the Fourier coefficients of f by
Chapter 9 The functional equation for the Riemann zeta function We will eventually deduce a functional equation, relating ζ(s to ζ( s. There are various methods to derive this functional equation, see
More informationSolar Sailing near a collinear point
Solar Sailing near a collinear point 6th AIMS International Conference on Dynamical Systems and Differential Equations 26 Ariadna Farrès & Àngel Jorba Departament de Matemàtica Aplicada i Anàlisi Universitat
More informationFinal Exam SOLUTIONS MAT 131 Fall 2011
1. Compute the following its. (a) Final Exam SOLUTIONS MAT 131 Fall 11 x + 1 x 1 x 1 The numerator is always positive, whereas the denominator is negative for numbers slightly smaller than 1. Also, as
More informationAlgorithms of Scientific Computing
Algorithms of Scientific Computing Discrete Fourier Transform (DFT) Michael Bader Technical University of Munich Summer 2018 Fast Fourier Transform Outline Discrete Fourier transform Fast Fourier transform
More informationAnalytic Trigonometry. Copyright Cengage Learning. All rights reserved.
Analytic Trigonometry Copyright Cengage Learning. All rights reserved. 7.4 Basic Trigonometric Equations Copyright Cengage Learning. All rights reserved. Objectives Basic Trigonometric Equations Solving
More informationChange point in trending regression
Charles University, Prague ROBUST 2010 joint work with A. Aue, L.Horváth, J. Picek (Charles University, Prague) 1 (Charles University, Prague) 1 2 model-formulation (Charles University, Prague) 1 2 model-formulation
More informationPart 3.3 Differentiation Taylor Polynomials
Part 3.3 Differentiation 3..3.1 Taylor Polynomials Definition 3.3.1 Taylor 1715 and Maclaurin 1742) If a is a fixed number, and f is a function whose first n derivatives exist at a then the Taylor polynomial
More informationReview of Concepts from Fourier & Filtering Theory. Fourier theory for finite sequences. convolution/filtering of infinite sequences filter cascades
Review of Concepts from Fourier & Filtering Theory precise definition of DWT requires a few basic concepts from Fourier analysis and theory of linear filters will start with discussion/review of: basic
More informationKey Intuition: invertibility
Introduction to Fourier Analysis CS 510 Lecture #6 January 30, 2017 In the extreme, a square wave Graphic from http://www.mechatronics.colostate.edu/figures/4-4.jpg 2 Fourier Transform Formally, the Fourier
More informationLecture # 06. Image Processing in Frequency Domain
Digital Image Processing CP-7008 Lecture # 06 Image Processing in Frequency Domain Fall 2011 Outline Fourier Transform Relationship with Image Processing CP-7008: Digital Image Processing Lecture # 6 2
More information2.3 Oscillation. The harmonic oscillator equation is the differential equation. d 2 y dt 2 r y (r > 0). Its solutions have the form
2. Oscillation So far, we have used differential equations to describe functions that grow or decay over time. The next most common behavior for a function is to oscillate, meaning that it increases and
More informationPhysics 70007, Fall 2009 Answers to Final Exam
Physics 70007, Fall 009 Answers to Final Exam December 17, 009 1. Quantum mechanical pictures a Demonstrate that if the commutation relation [A, B] ic is valid in any of the three Schrodinger, Heisenberg,
More informationElliptic Fourier Transform
NCUR, ART GRIGORYAN Frequency Analysis of the signals on ellipses: Elliptic Fourier Transform Joao Donacien N. Nsingui Department of Electrical and Computer Engineering The University of Texas at San Antonio
More informationPuzzle 1 Puzzle 2 Puzzle 3 Puzzle 4 Puzzle 5 /10 /10 /10 /10 /10
MATH-65 Puzzle Collection 6 Nov 8 :pm-:pm Name:... 3 :pm Wumaier :pm Njus 5 :pm Wumaier 6 :pm Njus 7 :pm Wumaier 8 :pm Njus This puzzle collection is closed book and closed notes. NO calculators are allowed
More informationin Electromagnetics Numerical Method Introduction to Electromagnetics I Lecturer: Charusluk Viphavakit, PhD
2141418 Numerical Method in Electromagnetics Introduction to Electromagnetics I Lecturer: Charusluk Viphavakit, PhD ISE, Chulalongkorn University, 2 nd /2018 Email: charusluk.v@chula.ac.th Website: Light
More informationChapter 10: Sinusoidal Steady-State Analysis
Chapter 10: Sinusoidal Steady-State Analysis 1 Objectives : sinusoidal functions Impedance use phasors to determine the forced response of a circuit subjected to sinusoidal excitation Apply techniques
More informationDiscrete-Time Signals: Time-Domain Representation
Discrete-Time Signals: Time-Domain Representation 1 Signals represented as sequences of numbers, called samples Sample value of a typical signal or sequence denoted as x[n] with n being an integer in the
More informationIAC-07-C SOLAR SAIL SURFING ALONG FAMILIES OF EQUILIBRIUM POINTS
58th International Astronautical Congress, Hyderabad, India, 24-28 September 2007. Copyright IAF/IAA. All rights reserved. IAC-07-C1.4.04 SOLAR SAIL SURFING ALONG FAMILIES OF EQUILIBRIUM POINTS Ariadna
More informationFinite Difference Methods for Boundary Value Problems
Finite Difference Methods for Boundary Value Problems October 2, 2013 () Finite Differences October 2, 2013 1 / 52 Goals Learn steps to approximate BVPs using the Finite Difference Method Start with two-point
More informationInterpolation via weighted l 1 -minimization
Interpolation via weighted l 1 -minimization Holger Rauhut RWTH Aachen University Lehrstuhl C für Mathematik (Analysis) Matheon Workshop Compressive Sensing and Its Applications TU Berlin, December 11,
More informationBifurcations thresholds of halo orbits
10 th AstroNet-II Final Meeting, 15 th -19 th June 2015, Tossa Del Mar 1/23 spazio Bifurcations thresholds of halo orbits Dr. Ceccaroni Marta ceccaron@mat.uniroma2.it University of Roma Tor Vergata Work
More informationarxiv:math/ v1 [math.mg] 31 May 2006
Covering spheres with spheres arxiv:math/060600v1 [math.mg] 31 May 006 Ilya Dumer College of Engineering, University of California at Riverside, Riverside, CA 951, USA dumer@ee.ucr.edu Abstract Given a
More informationH = ( H(x) m,n. Ω = T d T x = x + ω (d frequency shift) Ω = T 2 T x = (x 1 + x 2, x 2 + ω) (skewshift)
Chapter One Introduction We will consider infinite matrices indexed by Z (or Z b ) associated to a dynamical system in the sense that satisfies H = ( H(x) m,n )m,n Z H(x) m+1,n+1 = H(T x) m,n where x Ω,
More informationA remarkable representation of the Clifford group
A remarkable representation of the Clifford group Steve Brierley University of Bristol March 2012 Work with Marcus Appleby, Ingemar Bengtsson, Markus Grassl, David Gross and Jan-Ake Larsson Outline Useful
More informationGARCH Models Estimation and Inference
GARCH Models Estimation and Inference Eduardo Rossi University of Pavia December 013 Rossi GARCH Financial Econometrics - 013 1 / 1 Likelihood function The procedure most often used in estimating θ 0 in
More informationDecouplings and applications
April 27, 2018 Let Ξ be a collection of frequency points ξ on some curved, compact manifold S of diameter 1 in R n (e.g. the unit sphere S n 1 ) Let B R = B(c, R) be a ball with radius R 1. Let also a
More informationMathematics for Chemists 2 Lecture 14: Fourier analysis. Fourier series, Fourier transform, DFT/FFT
Mathematics for Chemists 2 Lecture 14: Fourier analysis Fourier series, Fourier transform, DFT/FFT Johannes Kepler University Summer semester 2012 Lecturer: David Sevilla Fourier analysis 1/25 Remembering
More informationFigure 10: Tangent vectors approximating a path.
3 Curvature 3.1 Curvature Now that we re parametrizing curves, it makes sense to wonder how we might measure the extent to which a curve actually curves. That is, how much does our path deviate from being
More informationCentral to this is two linear transformations: the Fourier Transform and the Laplace Transform. Both will be considered in later lectures.
In this second lecture, I will be considering signals from the frequency perspective. This is a complementary view of signals, that in the frequency domain, and is fundamental to the subject of signal
More informationLectures on Dynamical Systems. Anatoly Neishtadt
Lectures on Dynamical Systems Anatoly Neishtadt Lectures for Mathematics Access Grid Instruction and Collaboration (MAGIC) consortium, Loughborough University, 2007 Part 3 LECTURE 14 NORMAL FORMS Resonances
More informationEE301 Signals and Systems In-Class Exam Exam 3 Thursday, Apr. 19, Cover Sheet
EE301 Signals and Systems In-Class Exam Exam 3 Thursday, Apr. 19, 2012 Cover Sheet Test Duration: 75 minutes. Coverage: Chaps. 5,7 Open Book but Closed Notes. One 8.5 in. x 11 in. crib sheet Calculators
More informationGATE EE Topic wise Questions SIGNALS & SYSTEMS
www.gatehelp.com GATE EE Topic wise Questions YEAR 010 ONE MARK Question. 1 For the system /( s + 1), the approximate time taken for a step response to reach 98% of the final value is (A) 1 s (B) s (C)
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 12 Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted for noncommercial,
More informationQuasi-periodic solutions with Sobolev regularity of NLS on T d with a multiplicative potential
Quasi-periodic solutions with Sobolev regularity of NLS on T d with a multiplicative potential Massimiliano Berti, Philippe Bolle Abstract: We prove the existence of quasi-periodic solutions for Schrödinger
More informationGeorge Mason University Signals and Systems I Spring 2016
George Mason University Signals and Systems I Spring 206 Problem Set #6 Assigned: March, 206 Due Date: March 5, 206 Reading: This problem set is on Fourier series representations of periodic signals. The
More informationWeek 6 notes. dz = 0, it follows by taking limit that f has no zeros in B δ/2, contradicting the hypothesis that f(z 0 ) = 0.
Week 6 notes. Proof of Riemann Mapping Theorem Remark.. As far as the Riemann-mapping theorem and properties of the map near the boundary, I am following the approach in O. Costin notes, though I have
More informationCS711008Z Algorithm Design and Analysis
CS711008Z Algorithm Design and Analysis Lecture 5 FFT and Divide and Conquer Dongbo Bu Institute of Computing Technology Chinese Academy of Sciences, Beijing, China 1 / 56 Outline DFT: evaluate a polynomial
More informationThis is the number of cycles per unit time, and its units are, for example,
16 4. Sinusoidal solutions Many things in nature are periodic, even sinusoidal. We will begin by reviewing terms surrounding periodic functions. If an LTI system is fed a periodic input signal, we have
More informationIntroduction to the FFT
Introduction to the FFT 1 Introduction Assume that we have a signal fx where x denotes time. We would like to represent fx as a linear combination of functions e πiax or, equivalently, sinπax and cosπax
More informationInner products. Theorem (basic properties): Given vectors u, v, w in an inner product space V, and a scalar k, the following properties hold:
Inner products Definition: An inner product on a real vector space V is an operation (function) that assigns to each pair of vectors ( u, v) in V a scalar u, v satisfying the following axioms: 1. u, v
More informationFourier Transform with Rotations on Circles and Ellipses in Signal and Image Processing
Fourier Transform with Rotations on Circles and Ellipses in Signal and Image Processing Artyom M. Grigoryan Department of Electrical and Computer Engineering University of Texas at San Antonio One UTSA
More informationGARCH Models Estimation and Inference. Eduardo Rossi University of Pavia
GARCH Models Estimation and Inference Eduardo Rossi University of Pavia Likelihood function The procedure most often used in estimating θ 0 in ARCH models involves the maximization of a likelihood function
More informationLecture 4: Analysis of MIMO Systems
Lecture 4: Analysis of MIMO Systems Norms The concept of norm will be extremely useful for evaluating signals and systems quantitatively during this course In the following, we will present vector norms
More informationNotes for Expansions/Series and Differential Equations
Notes for Expansions/Series and Differential Equations In the last discussion, we considered perturbation methods for constructing solutions/roots of algebraic equations. Three types of problems were illustrated
More informationFourier transform. Stefano Ferrari. Università degli Studi di Milano Methods for Image Processing. academic year
Fourier transform Stefano Ferrari Università degli Studi di Milano stefano.ferrari@unimi.it Methods for Image Processing academic year 27 28 Function transforms Sometimes, operating on a class of functions
More informationPROBABILITY THEORY LECTURE 3
PROBABILITY THEORY LECTURE 3 Per Sidén Division of Statistics Dept. of Computer and Information Science Linköping University PER SIDÉN (STATISTICS, LIU) PROBABILITY THEORY - L3 1 / 15 OVERVIEW LECTURE
More informationAVERAGING AND RECONSTRUCTION IN HAMILTONIAN SYSTEMS
AVERAGING AND RECONSTRUCTION IN HAMILTONIAN SYSTEMS Kenneth R. Meyer 1 Jesús F. Palacián 2 Patricia Yanguas 2 1 Department of Mathematical Sciences University of Cincinnati, Cincinnati, Ohio (USA) 2 Departamento
More informationFigure Circuit for Question 1. Figure Circuit for Question 2
Exercises 10.7 Exercises Multiple Choice 1. For the circuit of Figure 10.44 the time constant is A. 0.5 ms 71.43 µs 2, 000 s D. 0.2 ms 4 Ω 2 Ω 12 Ω 1 mh 12u 0 () t V Figure 10.44. Circuit for Question
More informationCHAPTER 4 FOURIER SERIES S A B A R I N A I S M A I L
CHAPTER 4 FOURIER SERIES 1 S A B A R I N A I S M A I L Outline Introduction of the Fourier series. The properties of the Fourier series. Symmetry consideration Application of the Fourier series to circuit
More informationMATH 251 Week 6 Not collected, however you are encouraged to approach all problems to prepare for exam
MATH 51 Week 6 Not collected, however you are encouraged to approach all problems to prepare for exam A collection of previous exams could be found at the coordinator s web: http://www.math.psu.edu/tseng/class/m51samples.html
More informationMATH 127 SAMPLE FINAL EXAM I II III TOTAL
MATH 17 SAMPLE FINAL EXAM Name: Section: Do not write on this page below this line Part I II III TOTAL Score Part I. Multiple choice answer exercises with exactly one correct answer. Each correct answer
More informationAlexander Ostrowski
Ostrowski p. 1/3 Alexander Ostrowski 1893 1986 Walter Gautschi wxg@cs.purdue.edu Purdue University Ostrowski p. 2/3 Collected Mathematical Papers Volume 1 Determinants Linear Algebra Algebraic Equations
More informationNumerical Algorithms for Visual Computing II 2010/11 Example Solutions for Assignment 6
Numerical Algorithms for Visual Computing II 00/ Example Solutions for Assignment 6 Problem (Matrix Stability Infusion). The matrix A of the arising matrix notation U n+ = AU n takes the following form,
More informationMultifractality in simple systems. Eugene Bogomolny
Multifractality in simple systems Eugene Bogomolny Univ. Paris-Sud, Laboratoire de Physique Théorique et Modèles Statistiques, Orsay, France In collaboration with Yasar Atas Outlook 1 Multifractality Critical
More informationSome Collision solutions of the rectilinear periodically forced Kepler problem
Advanced Nonlinear Studies 1 (2001), xxx xxx Some Collision solutions of the rectilinear periodically forced Kepler problem Lei Zhao Johann Bernoulli Institute for Mathematics and Computer Science University
More informationPerturbation theory, KAM theory and Celestial Mechanics 7. KAM theory
Perturbation theory, KAM theory and Celestial Mechanics 7. KAM theory Alessandra Celletti Department of Mathematics University of Roma Tor Vergata Sevilla, 25-27 January 2016 Outline 1. Introduction 2.
More informationContinuous Wave Data Analysis: Fully Coherent Methods
Continuous Wave Data Analysis: Fully Coherent Methods John T. Whelan School of Gravitational Waves, Warsaw, 3 July 5 Contents Signal Model. GWs from rotating neutron star............ Exercise: JKS decomposition............
More informationSolution of the radiative transfer equation in NLTE stellar atmospheres
Solution of the radiative transfer equation in NLTE stellar atmospheres Jiří Kubát kubat@sunstel.asu.cas.cz Astronomický ústav AV ČR Ondřejov Non-LTE Line Formation for Trace Elements in Stellar Atmospheres,
More informationResummations of self energy graphs and KAM theory G. Gentile (Roma3), G.G. (I.H.E.S.) Communications in Mathematical Physics, 227, , 2002.
Resummations of self energy graphs and KAM theory G. Gentile (Roma3), G.G. (I.H.E.S.) Communications in Mathematical Physics, 227, 421 460, 2002. 1 Hamiltonian for a rotator system H = 1 2 (A2 + B 2 )+εf(α,
More informationDiscrete-Time Signals: Time-Domain Representation
Discrete-Time Signals: Time-Domain Representation 1 Signals represented as sequences of numbers, called samples Sample value of a typical signal or sequence denoted as x[n] with n being an integer in the
More information