Synchrosqueezed Curvelet Transforms for 2D mode decomposition
|
|
- Merry Greene
- 5 years ago
- Views:
Transcription
1 Synchrosqueezed Curvelet Transforms for 2D mode decomposition Haizhao Yang Department of Mathematics, Stanford University Collaborators: Lexing Ying and Jianfeng Lu Department of Mathematics and ICME, Stanford University Department of Mathematics, Duke University SIAM Conference on Imaging Science, May 24
2 Geophysics A superposition of several wave-like components. Wave field separation and ground roll removal problems
3 Materials science Atomic crystal analysis Observation: an assemblage of wave-like components; Goal: Crystal segmentations, crystal rotations, crystal defects, crystal deformations. Figure : Left: An atomic crystal image. Middle: Estimated crystal rotation. Right: Identified crystal defects. Takes about s for a 24*24 image, while other methods using GPU computing take at least 4s.
4 D mode decomposition Known: A superposition of wave-like components K K f (t) = f k (t) = α k (t)e 2πiN k φ k (t). k= k= Unknown: Number K, components f k (t), smooth instantaneous amplitudes α k (t), smooth instantaneous frequencies N k φ k (t). Existing methods: Empirical mode decomposition methods (Huang et al. 98, 9); Synchrosqueezed wavelet transform (Daubechies et al. 9, ); Synchrosqueezed wave packet transform (Y. 3); Data-driven time-frequency analysis (Hou et al., 2, 3); Regularized nonstationary autoregression (Fomel 3);
5 D synchrosqueezed wavelet transform (SSWT) Continuous wavelet transform of a signal s(t): W s (a, b) = s(t), φ ab (t) = s(t)a /2 φ( t b ) dt. a EX: s(t) = A cos(ωt). W s (a, b) = ŝ(ξ)a /2 φ(aξ)e ibξ dξ = A 2π 4π a/2 φ(aω)e ibω. Synchrosqueezing for better readability Figure : Numerical examples by Daubechies et al, signal f (t) = sin(8t).
6 Definitions of SSWT EX: s(t) = A cos(ωt). W s (a, b) = A 4π a/2 φ(aω)e ibω bw s (a,b) iw s (a,b) = ω Definition: Instantaneous frequency estimate ω s (a, b) = bw s (a, b) iw s (a, b). Definition: Synchrosqueezed wavelet transform T s (ω, b) = W s (a, b)a 3/2 δ(ω s (a, b) ω) da. {a:w s (a,b) }
7 Theory of SSWT Theorem: (Daubechies, Lu, Wu ACHA) If K K f (t) = f k (t) = α k (t)e 2πiN k φ k (t) k= and f k (t) are well-separated, then k= T f (a, b) has well-separated supports Z k concentrating (N k φ k (b), b); f k (t) can be accurately recovered by applying an inverse transform on I Zk (a, b)t f (a, b). where I Zk (a, b) is an indication function.
8 2D mode decomposition Known: A superposition of wave-like components K K f (x) = f k (x) = α k (x)e 2πiN k φ k (x). k= Unknown: Number K, components f k (x), local amplitudes α k (x), local wave vectors N k φ k (x). Existing methods: k=. 2D EMD methods (Huang et al.); 2. 2D Empirical wavelet, curvelet transforms (Gilles, Tran, Osher 3); 3. 2D SS wave packet and curvelet transforms (Y. and Ying 3, 4).
9 2D wavelet transform or wave packet transform? Consider two plane waves e ip x and e iq x again. ξ ξ ξ ξ ξ ξ Left: Supports of continuous wavelets and plane waves with different wave numbers in the Fourier domain. Radial separation. Middle: Supports of wavelets and plane waves with the same wave number. No angular separation. Right: Supports of wave packets and plane waves with the same wave number. Angular separation.
10 Definition: Wave Packets Given the mother wave packet w(x) and s (/2, ), the family of wave packets {w pb (x), p, b R 2 } is defined as or equivalently in Fourier domain w pb (x) = p s w( p s (x b))e 2πi(x b) p, ŵ pb (ξ) = p s e 2πib ξ ŵ( p s (ξ p)). Definition: The 2D wave packet transform of a function f (x) is a function of p, b R 2 W f (p, b) = w pb (x)f (x) dx.
11 Definition: The local wavevector estimation of a function f (x) at (p, b) is v f (p, b) = bw f (p, b) 2πiW f (p, b) for p, b R 2 such that W f (p, b). Definition: Given f (x), for v, b R 2, W f (p, b), and v f (p, b), the synchrosqueezed energy distribution T f (v, b) is T f (v, b) = W f (p, b) 2 δ(v f (p, b) v) dp. {p:w f (p,b) } Remark: The support of T f (v, b) is concentrating around local wavevectors.
12 p 2 v 2 v 2 v 2 Sketch of 2D mode decomposition Two components: f (x) = α (x)e 2πiNφ(x) + α 2 (x)e 2πiNφ2(x) p.5 b v.5 b v.5 b v.5 b 2 Fast algorithms for forward and inverse transforms. O(L 2 log L) for a L L image.
13 Theory of 2D SS wave packet transform (SSWPT) Theorem [Y.,Ying 3, SIAM Imaging Science] For a well-separated superposition f (x) = K k= α k(x)e 2πiNφ k (x) and ɛ >, we define and R f,ɛ = {(p, b) : W f (p, b) p s ɛ} Z f,k = {(p, b) : p N φ k (b) p s } for k K. There exists a constant ɛ (K) > such that for any ɛ (, ɛ ) there exists a constant N (K, ɛ) > such that for any N > N (K, ɛ) the following statements hold. (i) {Z f,k : k K} are disjoint and R f,ɛ k K Z f,k; (ii) For any (p, b) R f,ɛ Z f,k, v f (p, b) N φ k (b) N φ k (b) ɛ.
14 Banded wave-like components x x x x x x x x.2. Figure : Top-left: A banded deformed plane wave. Top-right: Number of wave packet coefficients W f (a, b) > δ. Bottom-left: Relative error of local wave-vector estimates using SSWPT. Bottom-right: Relative error of local wave-vector estimates using SSCT.
15 General curvelet transform: Some notations. The scaling matrix A a = ( ) a t a s. 2. The rotation angle θ and rotation matrix ( cos θ sin θ R θ = sin θ cos θ ). 3. The unit vector e θ = (cos θ, sin θ) T of rotation angle θ. 4. θ α represents the argument of given vector α. Definition: 2D general curvelets For 2 < s < t <, define w aθb (x) = a t+s 2 e 2πia(x b) e θ w(a a R (x b)). θ A family of curvelets {w aθb (x), a [, ), θ [, 2π), b R 2 }.
16 Definition: 2D general curvelet transform W f (a, θ, b) = w aθb (x)f (x) dx R 2 for a [, ), θ [, 2π), b R 2. Definition: local wave vector estimation Local wave-vector estimation at (a, θ, b) is v f (a, θ, b) = bw f (a, θ, b) 2πiW f (a, θ, b) for a [, ), θ [, 2π), b R 2 such that W f (a, θ, b). Definition: synchrosqueezed energy distribution Synchrosqueezed energy distribution T f (v, b) is T f (v, b) = W f (a, θ, b) 2 δ(rv f (a, θ, b) v)a da dθ for v R 2, b R 2. {(a,θ):w f (a,θ,b) }
17 Theory of the SS curvelet transform (SSCT) Theorem [Y., Ying SIAM Mathematical Analysis 4] Suppose 2 < s < η < t are fixed, and a well-separated superposition K f (x) = e (φ k (x) c k ) 2 /σ 2 k αk (x)e 2πiNφ k (x), k= where σ k N η. For any ɛ >, define { R f,ɛ = (a, θ, b) : W f (a, θ, b) a } s+t 2 ɛ and Z f,k = { (a, θ, b) : A a R θ (a e θ N φ k (b)) } for k K. For any ɛ, there exists N (ɛ) > such that for any N > N (ɛ) the following statements hold. (i) {Z f,k : k K} are disjoint and R f,ɛ k K Z f,k; (ii) For any (a, θ, b) R f,ɛ Z f,k, v f (a, θ, b) N φ k (b) N φ k (b) ɛ.
18 Applications in Geophysics Noisy Example for SS curvelet transform Figure : Top: A noisy superposition of two components. Left: First recovered component. Right: Second recovered component.
19 Noisy Example 2 for SS curvelet transform Figure : Left: A noisy superposition of two components. Right: The recovered main component.
20 A real example comment Figure : A superposition of several components.
21 Results of the real example Figure : Four recovered components.
22 Toy example of an atomic crystal image x x Energy 3.5 x Frequency Figure : Left: Crystal image. Crystal rotations, 7.5 degrees on the left and 45 on the right. Middle: 2D Fourier power spectrum. Right: Radially average Fourier power spectrum.
23 Results of this toy example x x Figure : Left: Crystal image. Middle: The angle estimates of crystal rotations. Right: The boundary indicator function.
24 Phase field crystal image Figure : Left: A crystal image and its zoomed-in image. Middle: The crystal rotation and its zoomed-in result. Right: The crystal defects and its zoomed-in result. Takes about s, while other methods take at least 4s.
25 Conclusion Developed 2D synchrosqueezed transforms to analyze images. A method with rigorous mathematical analysis. A method with both radial and angular separations. The SS curvelet transform can adapt different data structure: s and t. Fast algorithms to compute the forward and inverse transforms. Succesful applications in geophysics and materials science.
26 Results of 2D EEMD x x x x x x x x x x x x x Figure : A superposition of two components with similar local wave number.
27 Results of 2D Empirical curvelet transforms 2D Fourier power spetrum of a superposition of two components and Fourier domain partitions using empirical curvelet transforms.
Synchrosqueezed Transforms and Applications
Synchrosqueezed Transforms and Applications Haizhao Yang Department of Mathematics, Stanford University Collaborators: Ingrid Daubechies,JianfengLu ] and Lexing Ying Department of Mathematics, Duke University
More informationDownloaded 01/15/16 to Redistribution subject to SIAM license or copyright; see
MULTISCALE MODEL. SIMUL. Vol. 13, No. 4, pp. 1542 1572 c 215 Society for Industrial and Applied Mathematics Downloaded 1/15/16 to 171.64.38.112. Redistribution subject to SIAM license or copyright; see
More informationStatistical Analysis of Synchrosqueezed Transforms
Statistical Analysis of Synchrosqueezed Transforms Haizhao Yang Department of Mathematics, Due University October 4, Revised: August 6 Abstract Synchrosqueezed transforms are non-linear processes for a
More informationTime-frequency analysis of seismic data using synchrosqueezing wavelet transform a
Time-frequency analysis of seismic data using synchrosqueezing wavelet transform a a Published in Journal of Seismic Exploration, 23, no. 4, 303-312, (2014) Yangkang Chen, Tingting Liu, Xiaohong Chen,
More informationRobustness Analysis of Synchrosqueezed Transforms
Robustness Analysis of Synchrosqueezed Transforms Haizhao Yang Department of Mathematics, Due University October 4, Revised: July 5 Abstract Synchrosqueezed transforms are non-linear processes for a sharpened
More informationEmpirical Wavelet Transform
Jérôme Gilles Department of Mathematics, UCLA jegilles@math.ucla.edu Adaptive Data Analysis and Sparsity Workshop January 31th, 013 Outline Introduction - EMD 1D Empirical Wavelets Definition Experiments
More informationONE OR TWO FREQUENCIES? THE SYNCHROSQUEEZING ANSWERS
Advances in Adaptive Data Analysis Vol. 3, Nos. & 2 (2) 29 39 c World Scientific Publishing Company DOI:.42/S793536974X ONE OR TWO FREQUENCIES? THE SYNCHROSQUEEZING ANSWERS HAU-TIENG WU,, PATRICK FLANDRIN,
More informationEnergy during a burst of deceleration
Problem 1. Energy during a burst of deceleration A particle of charge e moves at constant velocity, βc, for t < 0. During the short time interval, 0 < t < t its velocity remains in the same direction but
More informationIterative Matching Pursuit and its Applications in Adaptive Time-Frequency Analysis
Iterative Matching Pursuit and its Applications in Adaptive Time-Frequency Analysis Zuoqiang Shi Mathematical Sciences Center, Tsinghua University Joint wor with Prof. Thomas Y. Hou and Sparsity, Jan 9,
More informationHAIZHAO YANG CURRENT POSITION
HAIZHAO YANG Department of Mathematics National Universify of Singapore Singapore Phone: 65162798 Email: matyh@nus.edu.sg Homepage: http://www.math.nus.edu.sg/ matyh/ Assistant Professor Department of
More informationNWC Distinguished Lecture. Blind-source signal decomposition according to a general mathematical model
NWC Distinguished Lecture Blind-source signal decomposition according to a general mathematical model Charles K. Chui* Hong Kong Baptist university and Stanford University Norbert Wiener Center, UMD February
More informationDifferential Equations: Homework 8
Differential Equations: Homework 8 Alvin Lin January 08 - May 08 Section.6 Exercise Find a general solution to the differential equation using the method of variation of parameters. y + y = tan(t) r +
More informationLinear second-order differential equations with constant coefficients and nonzero right-hand side
Linear second-order differential equations with constant coefficients and nonzero right-hand side We return to the damped, driven simple harmonic oscillator d 2 y dy + 2b dt2 dt + ω2 0y = F sin ωt We note
More informationInstantaneous Frequency Estimation Based on Synchrosqueezing Wavelet Transform
Instantaneous Frequency Estimation Based on Synchrosqueezing Wavelet Transform Qingtang Jiang and Bruce W. Suter February 216 1st revision in July 216 2nd revision in February 217 Abstract Recently, the
More informationHilbert-Huang and Morlet wavelet transformation
Hilbert-Huang and Morlet wavelet transformation Sonny Lion (sonny.lion@obspm.fr) LESIA, Observatoire de Paris The Hilbert-Huang Transform The main objective of this talk is to serve as a guide for understanding,
More informationApplied and Computational Harmonic Analysis
JID:YACHA AID:892 /FLA [m3g; v 1.85; Prn:31/10/2012; 14:04] P.1 (1-25) Appl. Comput. Harmon. Anal. ( ) Contents lists available at SciVerse ScienceDirect Applied and Computational Harmonic Analysis www.elsevier.com/locate/acha
More informationJim Lambers MAT 280 Summer Semester Practice Final Exam Solution. dy + xz dz = x(t)y(t) dt. t 3 (4t 3 ) + e t2 (2t) + t 7 (3t 2 ) dt
Jim Lambers MAT 28 ummer emester 212-1 Practice Final Exam olution 1. Evaluate the line integral xy dx + e y dy + xz dz, where is given by r(t) t 4, t 2, t, t 1. olution From r (t) 4t, 2t, t 2, we obtain
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 informationJean Morlet and the Continuous Wavelet Transform (CWT)
Jean Morlet and the Continuous Wavelet Transform (CWT) Brian Russell 1 and Jiajun Han 1 CREWES Adjunct Professor CGG GeoSoftware Calgary Alberta. www.crewes.org Introduction In 198 Jean Morlet a geophysicist
More informationSolutions of the wave equation as superpositions of its localized solutions
Solutions of the wave equation as superpositions of its localized solutions Department of Mathematical Physics, Physics Faculty St Petersburg University, Russian Federation July 26, 2011 Abstract of the
More informationDifferential Equations: Homework 2
Differential Equations: Homework Alvin Lin January 08 - May 08 Section.3 Exercise The direction field for provided x 0. dx = 4x y is shown. Verify that the straight lines y = ±x are solution curves, y
More informationIntegration by Substitution
November 22, 2013 Introduction 7x 2 cos(3x 3 )dx =? 2xe x2 +5 dx =? Chain rule The chain rule: d dx (f (g(x))) = f (g(x)) g (x). Use the chain rule to find f (x) and then write the corresponding anti-differentiation
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 informationJean Morlet and the Continuous Wavelet Transform
Jean Brian Russell and Jiajun Han Hampson-Russell, A CGG GeoSoftware Company, Calgary, Alberta, brian.russell@cgg.com ABSTRACT Jean Morlet was a French geophysicist who used an intuitive approach, based
More informationAn Introduction to HILBERT-HUANG TRANSFORM and EMPIRICAL MODE DECOMPOSITION (HHT-EMD) Advanced Structural Dynamics (CE 20162)
An Introduction to HILBERT-HUANG TRANSFORM and EMPIRICAL MODE DECOMPOSITION (HHT-EMD) Advanced Structural Dynamics (CE 20162) M. Ahmadizadeh, PhD, PE O. Hemmati 1 Contents Scope and Goals Review on transformations
More informationEdge preserved denoising and singularity extraction from angles gathers
Edge preserved denoising and singularity extraction from angles gathers Felix Herrmann, EOS-UBC Martijn de Hoop, CSM Joint work Geophysical inversion theory using fractional spline wavelets: ffl Jonathan
More informationCh 4: The Continuous-Time Fourier Transform
Ch 4: The Continuous-Time Fourier Transform Fourier Transform of x(t) Inverse Fourier Transform jt X ( j) x ( t ) e dt jt x ( t ) X ( j) e d 2 Ghulam Muhammad, King Saud University Continuous-time aperiodic
More informationDetermination of thin elastic inclusions from boundary measurements.
Determination of thin elastic inclusions from boundary measurements. Elena Beretta in collaboration with E. Francini, S. Vessella, E. Kim and J. Lee September 7, 2010 E. Beretta (Università di Roma La
More informationGaussian beams for high frequency waves: some recent developments
Gaussian beams for high frequency waves: some recent developments Jianliang Qian Department of Mathematics, Michigan State University, Michigan International Symposium on Geophysical Imaging with Localized
More informationSolving physics pde s using Gabor multipliers
Solving physics pde s using Michael P. Lamoureux, Gary F. Margrave, Safa Ismail ABSTRACT We develop the mathematical properties of, which are a nonstationary version of Fourier multipliers. Some difficulties
More informationShort-time Fourier transform for quaternionic signals
Short-time Fourier transform for quaternionic signals Joint work with Y. Fu and U. Kähler P. Cerejeiras Departamento de Matemática Universidade de Aveiro pceres@ua.pt New Trends and Directions in Harmonic
More informationComplex Analysis, Stein and Shakarchi The Fourier Transform
Complex Analysis, Stein and Shakarchi Chapter 4 The Fourier Transform Yung-Hsiang Huang 2017.11.05 1 Exercises 1. Suppose f L 1 (), and f 0. Show that f 0. emark 1. This proof is observed by Newmann (published
More informationHölder regularity estimation by Hart Smith and Curvelet transforms
Hölder regularity estimation by Hart Smith and Curvelet transforms Jouni Sampo Lappeenranta University Of Technology Department of Mathematics and Physics Finland 18th September 2007 This research is done
More information2D HILBERT-HUANG TRANSFORM. Jérémy Schmitt, Nelly Pustelnik, Pierre Borgnat, Patrick Flandrin
2D HILBERT-HUANG TRANSFORM Jérémy Schmitt, Nelly Pustelnik, Pierre Borgnat, Patrick Flandrin Laboratoire de Physique de l Ecole Normale Suprieure de Lyon, CNRS and Université de Lyon, France first.last@ens-lyon.fr
More informationFinal: Solutions Math 118A, Fall 2013
Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or
More informationSPARSE TIME FREQUENCY REPRESENTATIONS AND DYNAMICAL SYSTEMS
COMMUN. MATH. SCI. Vol. 3, No. 3, pp. 673 694 c 5 International Press SPARSE TIME FREQUENCY REPRESENTATIONS AND DYNAMICAL SYSTEMS THOMAS Y. HOU, ZUOQIANG SHI, AND PEYMAN TAVALLALI Dedicated to George Papanicolaou
More informationPROBLEM SET 6. E [w] = 1 2 D. is the smallest for w = u, where u is the solution of the Neumann problem. u = 0 in D u = h (x) on D,
PROBLEM SET 6 UE ATE: - APR 25 Chap 7, 9. Questions are either directl from the text or a small variation of a problem in the text. Collaboration is oka, but final submission must be written individuall.
More informationConstruction of Orthonormal Quasi-Shearlets based on quincunx dilation subsampling
Construction of Orthonormal Quasi-Shearlets based on quincunx dilation subsampling Rujie Yin Department of Mathematics Duke University USA Email: rujie.yin@duke.edu arxiv:1602.04882v1 [math.fa] 16 Feb
More informationFigure 1: Surface waves
4 Surface Waves on Liquids 1 4 Surface Waves on Liquids 4.1 Introduction We consider waves on the surface of liquids, e.g. waves on the sea or a lake or a river. These can be generated by the wind, by
More informationMathematics (Course B) Lent Term 2005 Examples Sheet 2
N12d Natural Sciences, Part IA Dr M. G. Worster Mathematics (Course B) Lent Term 2005 Examples Sheet 2 Please communicate any errors in this sheet to Dr Worster at M.G.Worster@damtp.cam.ac.uk. Note that
More informationCLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE
CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE 1. Linear Partial Differential Equations A partial differential equation (PDE) is an equation, for an unknown function u, that
More information2D Wavelets. Hints on advanced Concepts
2D Wavelets Hints on advanced Concepts 1 Advanced concepts Wavelet packets Laplacian pyramid Overcomplete bases Discrete wavelet frames (DWF) Algorithme à trous Discrete dyadic wavelet frames (DDWF) Overview
More informationLecture Notes 5: Multiresolution Analysis
Optimization-based data analysis Fall 2017 Lecture Notes 5: Multiresolution Analysis 1 Frames A frame is a generalization of an orthonormal basis. The inner products between the vectors in a frame and
More informationMassachusetts Institute of Technology Physics 8.03 Fall 2004 Final Exam Thursday, December 16, 2004
You have 3 hours Do all eight problems You may use calculators Massachusetts Institute of Technology Physics 8.03 Fall 004 Final Exam Thursday, December 16, 004 This is a closed-book exam; no notes are
More informationPHYSICS 110A : CLASSICAL MECHANICS HW 2 SOLUTIONS. Here is a sketch of the potential with A = 1, R = 1, and S = 1. From the plot we can see
PHYSICS 11A : CLASSICAL MECHANICS HW SOLUTIONS (1) Taylor 5. Here is a sketch of the potential with A = 1, R = 1, and S = 1. From the plot we can see 1.5 1 U(r).5.5 1 4 6 8 1 r Figure 1: Plot for problem
More informationL = 1 2 a(q) q2 V (q).
Physics 3550, Fall 2011 Motion near equilibrium - Small Oscillations Relevant Sections in Text: 5.1 5.6 Motion near equilibrium 1 degree of freedom One of the most important situations in physics is motion
More information1 Introduction. 2 Shannon Wavelet. Jaime Hernandez Jr. Wavelets and Fourier Analysis SEMESTER PROJECT December 12, 2007
Jaime Hernandez Jr. Wavelets and Fourier Analysis SEMESTER PROJECT December 12, 27 1 Introduction This report will address the construction of objects called Curvelets, which have been pioneered by Candes
More informationIntroduction to Biomedical Engineering
Introduction to Biomedical Engineering Biosignal processing Kung-Bin Sung 6/11/2007 1 Outline Chapter 10: Biosignal processing Characteristics of biosignals Frequency domain representation and analysis
More informationFundamentals of the gravitational wave data analysis V
Fundamentals of the gravitational wave data analysis V - Hilbert-Huang Transform - Ken-ichi Oohara Niigata University Introduction The Hilbert-Huang transform (HHT) l It is novel, adaptive approach to
More informationUNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH Faculty of Science and Engineering Department of Mathematics and Statistics END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MA4006 SEMESTER: Spring 2011 MODULE TITLE:
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 informationLegendre s Equation. PHYS Southern Illinois University. October 13, 2016
PHYS 500 - Southern Illinois University October 13, 2016 PHYS 500 - Southern Illinois University Legendre s Equation October 13, 2016 1 / 10 The Laplacian in Spherical Coordinates The Laplacian is given
More informationON NUMERICAL ANALYSIS AND EXPERIMENT VERIFICATION OF CHARACTERISTIC FREQUENCY OF ANGULAR CONTACT BALL-BEARING IN HIGH SPEED SPINDLE SYSTEM
ON NUMERICAL ANALYSIS AND EXPERIMENT VERIFICATION OF CHARACTERISTIC FREQUENCY OF ANGULAR CONTACT BALL-BEARING IN HIGH SPEED SPINDLE SYSTEM Tian-Yau Wu and Chun-Che Sun Department of Mechanical Engineering,
More informationkg meter ii) Note the dimensions of ρ τ are kg 2 velocity 2 meter = 1 sec 2 We will interpret this velocity in upcoming slides.
II. Generalizing the 1-dimensional wave equation First generalize the notation. i) "q" has meant transverse deflection of the string. Replace q Ψ, where Ψ may indicate other properties of the medium that
More informationIntroduction to Seismology Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 12.510 Introduction to Seismology Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 12.510 Introduction
More informationExtended isochron rays in prestack depth (map) migration
Extended isochron rays in prestack depth (map) migration A.A. Duchkov and M.V. de Hoop Purdue University, 150 N.University st., West Lafayette, IN, 47907 e-mail: aduchkov@purdue.edu (December 15, 2008)
More informationMATH 391 Test 1 Fall, (1) (12 points each)compute the general solution of each of the following differential equations: = 4x 2y.
MATH 391 Test 1 Fall, 2018 (1) (12 points each)compute the general solution of each of the following differential equations: (a) (b) x dy dx + xy = x2 + y. (x + y) dy dx = 4x 2y. (c) yy + (y ) 2 = 0 (y
More informationNew Multiscale Methods for 2D and 3D Astronomical Data Set
New Multiscale Methods for 2D and 3D Astronomical Data Set Jean Luc Starck Dapnia/SEDI SAP, CEA Saclay, France. jstarck@cea.fr http://jstarck.free.fr Collaborators: D.L. Donoho, O. Levi, Department of
More informationCalculus of Variations and Computer Vision
Calculus of Variations and Computer Vision Sharat Chandran Page 1 of 23 Computer Science & Engineering Department, Indian Institute of Technology, Bombay. http://www.cse.iitb.ernet.in/ sharat January 8,
More informationTraveling Harmonic Waves
Traveling Harmonic Waves 6 January 2016 PHYC 1290 Department of Physics and Atmospheric Science Functional Form for Traveling Waves We can show that traveling waves whose shape does not change with time
More informationVIBRATING BASE PENDULUM. Math 485 Project team Thomas Bello Emily Huang Fabian Lopez Kellin Rumsey Tao Tao
VIBRATING BASE PENDULUM Math 485 Project team Thomas Bello Emily Huang Fabian Lopez Kellin Rumsey Tao Tao Agenda Midterm Recap Equation of Motion & Energy Modeling Effective Potential Stability Analysis
More information2 Infinite products and existence of compactly supported φ
415 Wavelets 1 Infinite products and existence of compactly supported φ Infinite products.1 Infinite products a n need to be defined via limits. But we do not simply say that a n = lim a n N whenever the
More informationOscillatory Motion. Solutions of Selected Problems
Chapter 15 Oscillatory Motion. Solutions of Selected Problems 15.1 Problem 15.18 (In the text book) A block-spring system oscillates with an amplitude of 3.50 cm. If the spring constant is 250 N/m and
More informationHomework 3. 1 Goldstein Part (a) Theoretical Dynamics September 24, The Hamiltonian is given by
Theoretical Dynamics September 4, 010 Instructor: Dr. Thomas Cohen Homework 3 Submitted by: Vivek Saxena 1 Goldstein 8.1 1.1 Part (a) The Hamiltonian is given by H(q i, p i, t) = p i q i L(q i, q i, t)
More informationEilenberg-Steenrod properties. (Hatcher, 2.1, 2.3, 3.1; Conlon, 2.6, 8.1, )
II.3 : Eilenberg-Steenrod properties (Hatcher, 2.1, 2.3, 3.1; Conlon, 2.6, 8.1, 8.3 8.5 Definition. Let U be an open subset of R n for some n. The de Rham cohomology groups (U are the cohomology groups
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 informationTheoretical Estimation of Diffusion Coefficients for Binary Gas Mixtures
52 Appendix A Theoretical Estimation of Diffusion Coefficients for Binary Gas Mixtures The diffusion coefficient D 2 for the isothermal diffusion of species through constantpressure binary mixture of species
More informationMATH 131P: PRACTICE FINAL SOLUTIONS DECEMBER 12, 2012
MATH 3P: PRACTICE FINAL SOLUTIONS DECEMBER, This is a closed ook, closed notes, no calculators/computers exam. There are 6 prolems. Write your solutions to Prolems -3 in lue ook #, and your solutions to
More informationWave Phenomena Physics 15c. Lecture 10 Fourier Transform
Wave Phenomena Physics 15c Lecture 10 Fourier ransform What We Did Last ime Reflection of mechanical waves Similar to reflection of electromagnetic waves Mechanical impedance is defined by For transverse/longitudinal
More informationRotational Kinematics
Rotational Kinematics Rotational Coordinates Ridged objects require six numbers to describe their position and orientation: 3 coordinates 3 axes of rotation Rotational Coordinates Use an angle θ to describe
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 informationContinuous-Time Fourier Transform
Signals and Systems Continuous-Time Fourier Transform Chang-Su Kim continuous time discrete time periodic (series) CTFS DTFS aperiodic (transform) CTFT DTFT Lowpass Filtering Blurring or Smoothing Original
More informationSparse linear models
Sparse linear models Optimization-Based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_spring16 Carlos Fernandez-Granda 2/22/2016 Introduction Linear transforms Frequency representation Short-time
More informationPhysics 200 Lecture 4. Integration. Lecture 4. Physics 200 Laboratory
Physics 2 Lecture 4 Integration Lecture 4 Physics 2 Laboratory Monday, February 21st, 211 Integration is the flip-side of differentiation in fact, it is often possible to write a differential equation
More informationLINOEP vectors, spiral of Theodorus, and nonlinear time-invariant system models of mode decomposition
LINOEP vectors, spiral of Theodorus, and nonlinear time-invariant system models of mode decomposition Pushpendra Singh 1,2, 1 Department of EE, Indian Institute of Technology Delhi, India 2 Jaypee Institute
More informationFaculty of Computers and Information. Basic Science Department
18--018 FCI 1 Faculty of Computers and Information Basic Science Department 017-018 Prof. Nabila.M.Hassan 18--018 FCI Aims of Course: The graduates have to know the nature of vibration wave motions with
More informationSymmetries 2 - Rotations in Space
Symmetries 2 - Rotations in Space This symmetry is about the isotropy of space, i.e. space is the same in all orientations. Thus, if we continuously rotated an entire system in space, we expect the system
More informationHHT: the theory, implementation and application. Yetmen Wang AnCAD, Inc. 2008/5/24
HHT: the theory, implementation and application Yetmen Wang AnCAD, Inc. 2008/5/24 What is frequency? Frequency definition Fourier glass Instantaneous frequency Signal composition: trend, periodical, stochastic,
More informationPlane waves and spatial frequency. A plane wave
Plane waves and spatial frequency A plane wave Complex representation E(,) z t = E cos( ωt kz) = E cos( ωt kz) o Ezt (,) = Ee = Ee j( ωt kz) j( ωt kz) o = 1 2 A B t + + + [ cos(2 ω α β ) cos( α β )] {
More informationMultiresolution analysis & wavelets (quick tutorial)
Multiresolution analysis & wavelets (quick tutorial) Application : image modeling André Jalobeanu Multiresolution analysis Set of closed nested subspaces of j = scale, resolution = 2 -j (dyadic wavelets)
More informationAmplitude preserving AMO from true amplitude DMO and inverse DMO
Stanford Exploration Project, Report 84, May 9, 001, pages 1 167 Amplitude preserving AMO from true amplitude DMO and inverse DMO Nizar Chemingui and Biondo Biondi 1 ABSTRACT Starting from the definition
More informationIMAGE RECONSTRUCTION FROM UNDERSAMPLED FOURIER DATA USING THE POLYNOMIAL ANNIHILATION TRANSFORM
IMAGE RECONSTRUCTION FROM UNDERSAMPLED FOURIER DATA USING THE POLYNOMIAL ANNIHILATION TRANSFORM Anne Gelb, Rodrigo B. Platte, Rosie Renaut Development and Analysis of Non-Classical Numerical Approximation
More informationA SECOND-ORDER SYNCHROSQUEEZING S-TRANSFORM AND ITS APPLICATION IN SEISMIC SPECTRAL DECOMPOSITION
CHINESE JOURNAL OF GEOPHYSICS Vol.60, No.6, 2017, pp: 627 639 DOI: 10.6038/cjg20170728 A SECOND-ORDER SYNCHROSQUEEZING S-TRANSFORM AND ITS APPLICATION IN SEISMIC SPECTRAL DECOMPOSITION HUANG Zhong-Lai
More informationQuantum Mechanics in One Dimension. Solutions of Selected Problems
Chapter 6 Quantum Mechanics in One Dimension. Solutions of Selected Problems 6.1 Problem 6.13 (In the text book) A proton is confined to moving in a one-dimensional box of width.2 nm. (a) Find the lowest
More informationFirst Year Physics: Prelims CP1. Classical Mechanics: Prof. Neville Harnew. Problem Set III : Projectiles, rocket motion and motion in E & B fields
HT017 First Year Physics: Prelims CP1 Classical Mechanics: Prof Neville Harnew Problem Set III : Projectiles, rocket motion and motion in E & B fields Questions 1-10 are standard examples Questions 11-1
More informationJoint Distributions: Part Two 1
Joint Distributions: Part Two 1 STA 256: Fall 2018 1 This slide show is an open-source document. See last slide for copyright information. 1 / 30 Overview 1 Independence 2 Conditional Distributions 3 Transformations
More informationWavepackets. Outline. - Review: Reflection & Refraction - Superposition of Plane Waves - Wavepackets - ΔΔk Δx Relations
Wavepackets Outline - Review: Reflection & Refraction - Superposition of Plane Waves - Wavepackets - ΔΔk Δx Relations 1 Sample Midterm (one of these would be Student X s Problem) Q1: Midterm 1 re-mix (Ex:
More informationBernstein s inequality and Nikolsky s inequality for R d
Bernstein s inequality and Nikolsky s inequality for d Jordan Bell jordan.bell@gmail.com Department of athematics University of Toronto February 6 25 Complex Borel measures and the Fourier transform Let
More informationSection B. Electromagnetism
Prelims EM Spring 2014 1 Section B. Electromagnetism Problem 0, Page 1. An infinite cylinder of radius R oriented parallel to the z-axis has uniform magnetization parallel to the x-axis, M = m 0ˆx. Calculate
More information1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point
Solving a Linear System τ = trace(a) = a + d = λ 1 + λ 2 λ 1,2 = τ± = det(a) = ad bc = λ 1 λ 2 Classification of Fixed Points τ 2 4 1. < 0: the eigenvalues are real and have opposite signs; the fixed point
More informationWave Phenomena Physics 15c
Wave Phenomena Physics 5c Lecture Fourier Analysis (H&L Sections 3. 4) (Georgi Chapter ) What We Did Last ime Studied reflection of mechanical waves Similar to reflection of electromagnetic waves Mechanical
More informationWave equation techniques for attenuating multiple reflections
Wave equation techniques for attenuating multiple reflections Fons ten Kroode a.tenkroode@shell.com Shell Research, Rijswijk, The Netherlands Wave equation techniques for attenuating multiple reflections
More informationHow to extract the oscillating components of a signal? A wavelet-based approach compared to the Empirical Mode Decomposition
How to extract the oscillating components of a signal? A wavelet-based approach compared to the Empirical Mode Decomposition Adrien DELIÈGE University of Liège, Belgium Louvain-La-Neuve, February 7, 27
More informationMicrostates of AdS black holes and supersymmetric localization
Microstates of AdS black holes and supersymmetric localization Seyed Morteza Hosseini Università di Milano-Bicocca IPM, Tehran, May 8-11, 2017 Recent Trends in String Theory and Related Topics in collaboration
More informationPractice problems for Exam 1. a b = (2) 2 + (4) 2 + ( 3) 2 = 29
Practice problems for Exam.. Given a = and b =. Find the area of the parallelogram with adjacent sides a and b. A = a b a ı j k b = = ı j + k = ı + 4 j 3 k Thus, A = 9. a b = () + (4) + ( 3)
More information3: THE SHANNON SAMPLING THEOREM
3: THE SHANNON SAMPLING THEOEM STEVEN HEILMAN 1. Introduction In this final set of notes, we would like to end where we began. In the first set of notes, we mentioned that mathematics could be used to
More informationarxiv: v1 [cs.it] 18 Apr 2016
Time-Frequency analysis via the Fourier Representation Pushpendra Singh 1,2 1 Department of Electrical Engineering, IIT Delhi 2 Department of Electronics & Communication Engineering, JIIT Noida April 19,
More informationTensor intro 1. SIAM Rev., 51(3), Tensor Decompositions and Applications, Kolda, T.G. and Bader, B.W.,
Overview 1. Brief tensor introduction 2. Stein s lemma 3. Score and score matching for fitting models 4. Bringing it all together for supervised deep learning Tensor intro 1 Tensors are multidimensional
More informationTides in Higher-Dimensional Newtonian Gravity
Tides in Higher-Dimensional Newtonian Gravity Philippe Landry Department of Physics University of Guelph 23 rd Midwest Relativity Meeting October 25, 2013 Tides: A Familiar Example Gravitational interactions
More information12 th Class Mathematics Paper
th Class Mathematics Paper Maimum Time: hours Maimum Marks: 00 General Instructions: (i) All questions are compulsory. (ii) The question paper consists of 9 questions divided into four sections A, B, C
More information