Uncertainty and Spectrogram Geometry
|
|
- Easter Warren
- 6 years ago
- Views:
Transcription
1 From uncertainty to localization Spectrogram geometry CNRS & École Normale Supérieure de Lyon, France Erwin Schrödinger Institute, December 2012 * based on joint work with François Auger and Éric Chassande-Mottin
2 From uncertainty to localization Spectrogram geometry from or to and Heisenberg refined Heisenberg (classical) Second-order (variance-type) measures for the individual and spreadings of a signal x(t) L 2 (R) with spectrum X(ω): 2 t (x) := 1 t 2 x(t) 2 dt E x 2 ω(x) := 1 E x ω 2 X(ω) 2 dω 2π Theorem (Weyl, 27; Gabor, 46,... ) t (x) ω (X) 1 2, with lower bound attained for Gaussians x (t) = C e α t2, α < 0.
3 From uncertainty to localization Spectrogram geometry from or to and Heisenberg refined Heisenberg (-) Second-order (variance-type) measure for the joint - spreading of a signal x(t) L 2 (R) with energy distribution C x (t, ω; ϕ) within Cohen s class: tω (C x ) := 1 E x ( t 2 ) T 2 + T 2 ω 2 C x (t, ω; ϕ) dt dω 2π Theorem (Janssen, 91) Wigner tω (W x ) 1 Spectrogram tω (S h x ) 2, with lower bound attained for Gaussian signals and matched Gaussian windows.
4 From uncertainty to localization Spectrogram geometry from or to and Heisenberg refined - covariance Definition (Cohen, 95) c(x) := t x(t) 2 d arg x(t) dt dt 1 Interpretation: covariance quantifies the coupling between and instantaneous c(x) = t ω x (t) 2 Intuition: covariance is zero in case of no coupling c(x) = t ω x (t) = t ω x (t) = t ω = 0
5 From uncertainty to localization Spectrogram geometry from or to and Heisenberg refined Schrödinger Theorem (Schrödinger, 30) t (x) ω (X) c 2 (x), with lower bound attained for Gaussians x (t) = e α t2 +β t+γ, Re{α} < 0. 1 Terminology: waveforms x attaining the lower bound correspond to squeezed states in quantum mechanics and linear chirps in signal theory 2 Interpretation: possibility of localization in the plane beyond pointwise energy concentration
6 From uncertainty to localization Spectrogram geometry sharp localization reassignment chirps, Hermite and logons chirp localization
7 From uncertainty to localization Spectrogram geometry sharp localization reassignment chirps, Hermite and logons how? 1 Effective localization: in the Wigner case, we have lim W x (t, ω) = δ (ω (β + 2 Im{α}t)) Re{α} 0 2 Generalization: localization on more arbitrary curves of the plane by modifying the symmetry rules underlying the Wigner distribution (Gonçalvès and F., 96) 3 Caveat: global localization holds for monocomponent signals only 4 Way out: reassignment (Kodera et al., 76)
8 From uncertainty to localization Spectrogram geometry sharp localization reassignment chirps, Hermite and logons rationale Starting with the smoothing relationship Sx h (t, ω) = W x (s, ξ) W h (s t, ξ ω) dt dω 2π, the key idea is 1 to replace the geometrical center of the smoothing - domain by the center of mass of the Wigner distribution over this domain, and 2 to reassign the value of the smoothed distribution to this local centroïd: ( ) Ŝx h (t, ω) = Sx h (τ, ξ)δ t ˆt x (τ, ξ), ω ˆω x (τ, ξ) dτ dξ 2π
9 From uncertainty to localization Spectrogram geometry sharp localization reassignment chirps, Hermite and logons spectrogram = smoothed Wigner Wigner-Ville spectrogram
10 From uncertainty to localization Spectrogram geometry sharp localization reassignment chirps, Hermite and logons spreading of auto-terms Wigner-Ville spectrogram
11 From uncertainty to localization Spectrogram geometry sharp localization reassignment chirps, Hermite and logons cancelling of cross-terms Wigner-Ville spectrogram
12 From uncertainty to localization Spectrogram geometry sharp localization reassignment chirps, Hermite and logons reassignment (Kodera et al., 76, Auger & F., 95) Wigner-Ville reassigned spectrogram
13 From uncertainty to localization Spectrogram geometry sharp localization reassignment chirps, Hermite and logons reassignment in action 1 Spectrogram: implicit computation of the local centroïds (Auger & F., 95) via two additional short- Fourier transforms (STFTs) with the companion windows (T h)(t) = t h(t) and (Dh)(t) = (dh/dt)(t) 2 Gaussian spectrogram: alternative computation of the local centroïds via partial derivatives of the magnitude of the STFT (Auger, Chassande-Mottin & F., 12) 3 Beyond spectrograms: possible generalizations to other smoothings (smoothed pseudo-wigner-ville, scalogram, etc.)
14 From uncertainty to localization Spectrogram geometry sharp localization reassignment chirps, Hermite and logons linear chirps Framework: Gaussian window h(t) = π 1/4 e t2 /2 Result: the linear chirp parameterized by α = 1 ( ) 1 2 T 2 i a ; β = 0 is such that lim Ŝx h (t, ω) = 1 δ(ω at) T 2π Interpretation: perfect localization
15 From uncertainty to localization Spectrogram geometry sharp localization reassignment chirps, Hermite and logons Hermite functions Framework: Gaussian window g(t) = π 1/4 e t2 /2 Hermite functions: with h k (t) = C k 1 1 T H k 1 ( 2πt/T )e π(t/t )2, H n (α) = ( 1) n e α2 dα n e α2, n = 0, 1 Radial reassignment: d n ˆω(t, ω) ˆt(t, ω) = ω t
16 From uncertainty to localization Spectrogram geometry sharp localization reassignment chirps, Hermite and logons Hermite functions Result (F., JFAA 12): Ŝ (g) h k (t, ω) = with ˆV δ and: 1 2 k 1 (k 1)! ˆV ) ( t 2 + ω 2 1 D (t, ω), D = {(t, ω) t 2 + ω 2 t 2 m + ω 2 m = 2(k 1)} Interpretation: clearance area and no perfect localization Asymptotic localization: Ŝ (g) h k (t, ω) k δ(t 2 + ω 2 2(k 1))
17 From uncertainty to localization Spectrogram geometry sharp localization reassignment chirps, Hermite and logons Hermite functions M = 2 Wigner spectro. reass. spectro.
18 From uncertainty to localization Spectrogram geometry sharp localization reassignment chirps, Hermite and logons Hermite functions order = 2 (theory) order = 2 (simulation)
19 From uncertainty to localization Spectrogram geometry sharp localization reassignment chirps, Hermite and logons Hermite functions M = 7 Wigner spectro. reass. spectro.
20 From uncertainty to localization Spectrogram geometry sharp localization reassignment chirps, Hermite and logons Hermite functions order = 7 (theory) order = 7 (simulation)
21 From uncertainty to localization Spectrogram geometry sharp localization reassignment chirps, Hermite and logons Hermite functions M = 18 Wigner spectro. reass. spectro.
22 From uncertainty to localization Spectrogram geometry sharp localization reassignment chirps, Hermite and logons Hermite functions order = 18 (theory) order = 18 (simulation)
23 From uncertainty to localization Spectrogram geometry sharp localization reassignment chirps, Hermite and logons logons Framework: Gaussian window h(t) = π 1/4 e t2 /2 Result: the Gabor logon, i.e., the previous linear chirp with a = 0, is such that S h x (t, ω) = e 1 2 (t2 +ω 2) tω (S h x ) = 2 W x (t, ω) = 2 e (t2 +ω 2) tω (W x ) = 1 Ŝ h x (t, ω) = 4 e 2(t2 +ω 2) tω (Ŝh x ) = 1 2 Interpretation: Heisenberg defeated?
24 From uncertainty to localization Spectrogram geometry sharp localization reassignment chirps, Hermite and logons resolving the paradox 1 General remark: as for Fourier, sharp localization not to be confused with resolution (i.e., ability to separate closely spaced components) 2 Complete picture: reassignment = squeezed distribution + vector field r x (t, ω) = (ˆt x (t, ω) t, ˆω x (t, ω) ω) t such that r x (t, ω) = 1 2 log Sh x (t, ω) 3 Interpretation: basins of attraction 4 Open question: geometry of such basins and relationship with Heisenberg?
25 From uncertainty to localization Spectrogram geometry on extrema Voronoi and Delaunay a simplified model factorization Bargmann representation: with a circular Gaussian window, the STFT can be expressed as F h x (z) = F h x (z) e z 2 /4, where z = ω + jt and F h x (z) is an entire function of order at most 2 Weierstrass-Hadamard factorization: it follows that Fx h (z) = e ( ) Q(z) (1 z n ) exp z n + z n/2 2, n where Q(z) is a quadratic polynomial and z n = z/z n Interpretation: complete characterization by zeroes and, by duality, by local maxima
26 From uncertainty to localization Spectrogram geometry on extrema Voronoi and Delaunay a simplified model more on zeroes Result (Balasz et al., 11) Universal singularity of the phase derivative around zeroes of the STFT
27 From uncertainty to localization Spectrogram geometry on extrema Voronoi and Delaunay a simplified model phase-magnitude relationship In the Gaussian case, simple explanation based on: Proposition (Auger et al., 12) Phase and magnitude of a Gaussian STFT form a Cauchy pair, i.e., Φ h x (t, ω) = t ω log Mh x(t, ω) Φ h x ω (t, ω) = t log Mh x(t, ω), with F h x (t, ω) = M h x(t, ω) e jφh x (t,ω) e (t2 +ω 2 )/4
28 From uncertainty to localization Spectrogram geometry on extrema Voronoi and Delaunay a simplified model phase singularity 1 Magnitude: local behaviour in the vicinity of a zero (t 0, ω 0 ) M h x(t, ω) [(ω ω 0 ) 2 + (t t 0 ) 2] 1/2 2 Phase from magnitude: universal hyperbolic divergence of the phase derivative Φ h x (t 0, ω) t (ω ω 0 ) 1 ω ω0 Φ h x ω (t, ω 0) (t 0 t) 1 t t0
29 From uncertainty to localization Spectrogram geometry on extrema Voronoi and Delaunay a simplified model zeroes phase dislocations (Nye & Berry, 74) magnitude phase gradient in phase gradient in phase
30 From uncertainty to localization Spectrogram geometry on extrema Voronoi and Delaunay a simplified model phase dislocation ( domain)
31 From uncertainty to localization Spectrogram geometry on extrema Voronoi and Delaunay a simplified model phase dislocation (complex plane)
32 From uncertainty to localization Spectrogram geometry on extrema Voronoi and Delaunay a simplified model phase dislocation (zoom)
33 From uncertainty to localization Spectrogram geometry on extrema Voronoi and Delaunay a simplified model more on extrema Proposed approach 1 Data: white Gaussian noise 2 Time- representation: Gaussian spectrogram 3 Characterization: identification of local extrema (zeroes and maxima) + Voronoi tessellation and Delaunay triangulation 4 Analysis: distribution of cell areas, lengths between extrema and heights of local maxima
34 From uncertainty to localization Spectrogram geometry on extrema Voronoi and Delaunay a simplified model an example spectrogram (wgn) spectrogram (logon) Voronoi/Delaunay (min) Voronoi/Delaunay (max)
35 From uncertainty to localization Spectrogram geometry on extrema Voronoi and Delaunay a simplified model mean arrangement Voronoi/Delaunay spectrogram (logon) 1 Observation: average connectivity 6 tiling with hexagonal cells 2 Interpretation: maximum packing of circular patches
36 From uncertainty to localization Spectrogram geometry on extrema Voronoi and Delaunay a simplified model predictions D M /d m = 3; N M /N m = 1/3; A M /A m = 3
37 From uncertainty to localization Spectrogram geometry on extrema Voronoi and Delaunay a simplified model simulation results distance : 5.24 number: 24 area : 17.6 maxima 0 10 distance : number: area : 6.3 minima
38 From uncertainty to localization Spectrogram geometry on extrema Voronoi and Delaunay a simplified model comparison with model (d M /d m )/sqrt(3) (A m /A M )/ mode theory median mean number of bins 0.5 mode theory median mean number of bins 1.5 (A M /d M 2 )/(sqrt(3)/2) 1.5 C/ mode theory median mean number of bins 0.5 mode theory median mean number of bins
39 From uncertainty to localization Spectrogram geometry on extrema Voronoi and Delaunay a simplified model analysis 1 Distributions: significant dispersion for lengths and areas 2 Ratios max/min: reasonable agreement between experimental results and theoretical predictions 3 Ranges of values: if we call effective domain of the minimum uncertainty logon the circular domain which encompasses 95% of its energy, its radius and area are equal to 2.6 and 21.8, to be compared to the values d M / 3 3 and 2π/(3 3)A M 21.8 attached to the hexagonal tiling 4 Interpretation: tiling cells minimum uncertainty logons
40 From uncertainty to localization Spectrogram geometry on extrema Voronoi and Delaunay a simplified model local maxima and Voronoi cells areas (1) 1 Heights: well-described by a Gamma distribution 2 Areas: idem and similar to the heights distribution for a proper renormalization 2 (t 2 + ω 2 ) Fx h (t, ω) 2 dt dω ( 2π F R 2 (A) 2 Proposition A The value F of a local maximum of the STFT magnitude and the area A of the associated Voronoi cell satisfy the uncertainty-type inequality A. F 3 6 ) 2
41 From uncertainty to localization Spectrogram geometry on extrema Voronoi and Delaunay a simplified model local maxima and Voronoi cells areas (2) areas of Voronoi cells (red) and local maxima of STFT magnitude (blue) p.d.f o : data + : Gamma fit 0.02 areas of Voronoi cells joint p.d.f local maxima of STFT magnitude
42 From uncertainty to localization Spectrogram geometry on extrema Voronoi and Delaunay a simplified model towards a random Gabor expansion Previous results suggest a possible modeling of Gabor spectrograms as with Sx h (t, ω) = c mn Fh h (t t m, ω ω n ) m n 1 locations (t m, ω n ) distributed on some suitable randomized version of a triangular grid 2 magnitudes of the weights c mn Gamma-distributed 3 locations and weights partly correlated 2
43 From uncertainty to localization Spectrogram geometry on extrema Voronoi and Delaunay a simplified model concluding remarks 1 Spectrogram: geometry constrained by - uncertainty 2 Reassignment: access to simplified properties, both geometrical and statistical 3 Gabor case: complete characterization by zeroes or local maxima? 4 Modeling: new ways of sparsifying energy distributions?
On Time-Frequency Sparsity and Uncertainty
CNRS & E cole Normale Supe rieure de Lyon, France * partially based on joint works with Franc ois Auger, Pierre Borgnat and E ric Chassande-Mottin Heisenberg (classical) from or to and Heisenberg refined
More informationUNCERTAINTY AND SPECTROGRAM GEOMETRY
UNCERTAINTY AND SPECTROGRAM GEOMETRY Patrick Flandrin 1 Éric Chassande-Mottin François Auger 3 1 Laboratoire de Physique de l École Normale Supérieure de Lyon, CNRS and Université de Lyon, France Laboratoire
More informationDrawing sounds, listening to images The art of time-frequency an
Fourier notes localization oscillations Drawing sounds, listening to images The art of time-frequency analysis CNRS & École Normale Supérieure de Lyon, France CIRMMT Montréal (CA), April. 9, Fourier notes
More informationthree observations representations and tools three examples Chirps everywhere Patrick Flandrin CNRS & École Normale Supérieure de Lyon, France
three observations representations and tools three examples CNRS & École Normale Supérieure de Lyon, France three observations representations and tools three examples Euler s disk pendulum Doppler Leonhard
More informationMULTITAPER TIME-FREQUENCY REASSIGNMENT. Jun Xiao and Patrick Flandrin
MULTITAPER TIME-FREQUENCY REASSIGNMENT Jun Xiao and Patrick Flandrin Laboratoire de Physique (UMR CNRS 5672), École Normale Supérieure de Lyon 46, allée d Italie 69364 Lyon Cedex 07, France {Jun.Xiao,flandrin}@ens-lyon.fr
More informationIntroduction to time-frequency analysis. From linear to energy-based representations
Introduction to time-frequency analysis. From linear to energy-based representations Rosario Ceravolo Politecnico di Torino Dep. Structural Engineering UNIVERSITA DI TRENTO Course on «Identification and
More informationA note on the time-frequency analysis of GW150914
A note on the - analysis of GW15914 Patrick Flandrin To cite this version: Patrick Flandrin. A note on the - analysis of GW15914. [Research Report] Ecole normale supérieure de Lyon. 216.
More informationTIME-FREQUENCY ANALYSIS AND HARMONIC GAUSSIAN FUNCTIONS
TIME-FREQUENCY ANALYSIS AND HARMONIC GAUSSIAN FUNCTIONS Tokiniaina Ranaivoson, Raoelina Andriambololona, Rakotoson Hanitriarivo Theoretical Physics Department Institut National des Sciences et Techniques
More informationIntroduction to Time-Frequency Distributions
Introduction to Time-Frequency Distributions Selin Aviyente Department of Electrical and Computer Engineering Michigan State University January 19, 2010 Motivation for time-frequency analysis When you
More informationChirp Rate and Instantaneous Frequency Estimation: Application to Recursive Vertical Synchrosqueezing
Chirp ate and Instantaneous Frequency Estimation: Application to ecursive Vertical Synchrosqueezing Dominique Fourer, François Auger, Krzysztof Czarnecki, Sylvain Meignen, Patrick Flandrin To cite this
More informationTheory and applications of time-frequency analysis
Theory and applications of time-frequency analysis Ville Turunen (ville.turunen@aalto.fi) Abstract: When and how often something happens in a signal? By properly quantizing these questions, we obtain the
More informationShift Covariant Time-Frequency Distributions of Discrete. Signals. Jerey C. O'Neill. of the requirements for the degree of. Doctor of Philosophy
Shift Covariant Time-Frequency Distributions of Discrete Signals by Jerey C. O'Neill A dissertation submitted in partial fulllment of the requirements for the degree of Doctor of Philosophy (Electrical
More informationLecture Hilbert-Huang Transform. An examination of Fourier Analysis. Existing non-stationary data handling method
Lecture 12-13 Hilbert-Huang Transform Background: An examination of Fourier Analysis Existing non-stationary data handling method Instantaneous frequency Intrinsic mode functions(imf) Empirical mode decomposition(emd)
More informationDigital Image Processing Lectures 13 & 14
Lectures 13 & 14, Professor Department of Electrical and Computer Engineering Colorado State University Spring 2013 Properties of KL Transform The KL transform has many desirable properties which makes
More informationTime-Frequency Toolbox For Use with MATLAB
Time-Frequency Toolbox For Use with MATLAB François Auger * Patrick Flandrin * Paulo Gonçalvès Olivier Lemoine * * CNRS (France) Rice University (USA) 1995-1996 2 3 Copyright (C) 1996 CNRS (France) and
More informationTIME-FREQUENCY ANALYSIS: TUTORIAL. Werner Kozek & Götz Pfander
TIME-FREQUENCY ANALYSIS: TUTORIAL Werner Kozek & Götz Pfander Overview TF-Analysis: Spectral Visualization of nonstationary signals (speech, audio,...) Spectrogram (time-varying spectrum estimation) TF-methods
More informationTime-frequency localization from sparsity constraints
Time- localization from sparsity constraints Pierre Borgnat, Patrick Flandrin To cite this version: Pierre Borgnat, Patrick Flandrin. Time- localization from sparsity constraints. 4 pages, 3 figures, 1
More informationInformation geometry for bivariate distribution control
Information geometry for bivariate distribution control C.T.J.Dodson + Hong Wang Mathematics + Control Systems Centre, University of Manchester Institute of Science and Technology Optimal control of stochastic
More informationComputational Harmonic Analysis (Wavelet Tutorial) Part II
Computational Harmonic Analysis (Wavelet Tutorial) Part II Understanding Many Particle Systems with Machine Learning Tutorials Matthew Hirn Michigan State University Department of Computational Mathematics,
More informationExtremal Bounds of Gaussian Gabor Frames and Properties of Jacobi s Theta Functions
Extremal Bounds of Gaussian Gabor Frames and Properties of Jacobi s Theta Functions Supervisor: Univ.Prof. Dr. Karlheinz Gröchenig Public Defense of Doctoral Thesis February 28, 2017 Contents 1 Gabor Systems
More informationAccuracy and Precision in PDV Data Analysis
Accuracy and Precision in PDV Data Analysis Michael R. Furlanetto P-23, LANL Second Annual Heterodyne Workshop, LLNL, August 17, 2007 LA-UR-07-5610 Motivation Given a PDV data set [A(t)], with what accuracy
More informationComparative study of different techniques for Time-Frequency Analysis
Comparative study of different techniques for Time-Frequency Analysis A.Vishwadhar M.Tech Student Malla Reddy Institute Of Technology And Science,Maisammaguda, Dulapally, Secunderabad. Abstract-The paper
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 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 informationTIME-FREQUENCY ANALYSIS AND HARMONIC GAUSSIAN FUNCTIONS
TIME-FREQUENCY ANALYSIS AND HARMONIC GAUSSIAN FUNCTIONS Tokiniaina Ranaivoson *, Raoelina Andriambololona **, Rakotoson Hanitriarivo *** Theoretical Physics Department Institut National des Sciences et
More informationIN NONSTATIONARY contexts, it is well-known [8] that
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 55, NO. 6, JUNE 2007 2851 Multitaper Time-Frequency Reassignment for Nonstationary Spectrum Estimation and Chirp Enhancement Jun Xiao and Patrick Flandrin,
More informationThe ASTRES Toolbox for Mode Extraction of Non-Stationary Multicomponent Signals
The ASTES Toolbo for Mode Etraction of Non-Stationary Multicomponent Signals Dominique Fourer, Jinane Harmouche, Jérémy Schmitt, Thomas Oberlin, Sylvain Meignen, François Auger and Patrick Flandrin dominique@fourer.fr
More informationLecture 3 Dynamics 29
Lecture 3 Dynamics 29 30 LECTURE 3. DYNAMICS 3.1 Introduction Having described the states and the observables of a quantum system, we shall now introduce the rules that determine their time evolution.
More informationSignature Analysis of Mechanical Watch Movements by Reassigned Spectrogram
Proceedings of the 6th WSEAS International Conference on Signal Processing, Robotics and Automation, Corfu Island, Greece, February 16-19, 2007 177 Signature Analysis of Mechanical Watch Movements by Reassigned
More informationVector Derivatives and the Gradient
ECE 275AB Lecture 10 Fall 2008 V1.1 c K. Kreutz-Delgado, UC San Diego p. 1/1 Lecture 10 ECE 275A Vector Derivatives and the Gradient ECE 275AB Lecture 10 Fall 2008 V1.1 c K. Kreutz-Delgado, UC San Diego
More informationLecture 3 Kernel properties and design in Cohen s class time-frequency distributions
Lecture 3 Kernel properties and design in Cohen s class time-frequency distributions Time-frequency analysis, adaptive filtering and source separation José Biurrun Manresa 22.02.2011 Time-Frequency representations
More informationTHE POLE BEHAVIOUR OF THE PHASE DERIVATIVE OF THE SHORT-TIME FOURIER TRANSFORM
THE POLE BEHAVIOUR OF THE PHASE DERIVATIVE OF THE SHORT-TIME FOURIER TRANSFORM PETER BALAZS A), DOMINIK BAYER A), FLORENT JAILLET B) AND PETER SØNDERGAARD C) Abstract. The short-time Fourier transform
More informationMedian Filter Based Realizations of the Robust Time-Frequency Distributions
TIME-FREQUENCY SIGNAL ANALYSIS 547 Median Filter Based Realizations of the Robust Time-Frequency Distributions Igor Djurović, Vladimir Katkovnik, LJubiša Stanković Abstract Recently, somenewefficient tools
More informationELEG 833. Nonlinear Signal Processing
Nonlinear Signal Processing ELEG 833 Gonzalo R. Arce Department of Electrical and Computer Engineering University of Delaware arce@ee.udel.edu February 15, 25 2 NON-GAUSSIAN MODELS 2 Non-Gaussian Models
More informationResearch Article Attenuation Analysis of Lamb Waves Using the Chirplet Transform
Hindawi Publishing Corporation EURASIP Journal on Advances in Signal Processing Volume 00, Article ID 3757, 6 pages doi:0.55/00/3757 Research Article Attenuation Analysis of Lamb Waves Using the Chirplet
More informationIntroduction to time-frequency analysis Centre for Doctoral Training in Healthcare Innovation
Introduction to time-frequency analysis Centre for Doctoral Training in Healthcare Innovation Dr. Gari D. Clifford, University Lecturer & Director, Centre for Doctoral Training in Healthcare Innovation,
More informationTIME-FREQUENCY ANALYSIS EE3528 REPORT. N.Krishnamurthy. Department of ECE University of Pittsburgh Pittsburgh, PA 15261
TIME-FREQUENCY ANALYSIS EE358 REPORT N.Krishnamurthy Department of ECE University of Pittsburgh Pittsburgh, PA 56 ABSTRACT - analysis, is an important ingredient in signal analysis. It has a plethora of
More information(x k ) sequence in F, lim x k = x x F. If F : R n R is a function, level sets and sublevel sets of F are any sets of the form (respectively);
STABILITY OF EQUILIBRIA AND LIAPUNOV FUNCTIONS. By topological properties in general we mean qualitative geometric properties (of subsets of R n or of functions in R n ), that is, those that don t depend
More informationMath 302 Outcome Statements Winter 2013
Math 302 Outcome Statements Winter 2013 1 Rectangular Space Coordinates; Vectors in the Three-Dimensional Space (a) Cartesian coordinates of a point (b) sphere (c) symmetry about a point, a line, and a
More informationReassigned Scalograms and Singularities
Reassigned Scalograms and Singularities Eric Chassande-Mottin and Patrick Flandrin Abstract. Reassignment is a general nonlinear technique aimed at increasing the localization of time-frequency and time-scale
More informationNonlinear Signal Processing ELEG 833
Nonlinear Signal Processing ELEG 833 Gonzalo R. Arce Department of Electrical and Computer Engineering University of Delaware arce@ee.udel.edu May 5, 2005 8 MYRIAD SMOOTHERS 8 Myriad Smoothers 8.1 FLOM
More information04. Random Variables: Concepts
University of Rhode Island DigitalCommons@URI Nonequilibrium Statistical Physics Physics Course Materials 215 4. Random Variables: Concepts Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative
More informationLQR, Kalman Filter, and LQG. Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin
LQR, Kalman Filter, and LQG Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin May 2015 Linear Quadratic Regulator (LQR) Consider a linear system
More informationNotes on uniform convergence
Notes on uniform convergence Erik Wahlén erik.wahlen@math.lu.se January 17, 2012 1 Numerical sequences We begin by recalling some properties of numerical sequences. By a numerical sequence we simply mean
More informationFourier Series Example
Fourier Series Example Let us compute the Fourier series for the function on the interval [ π,π]. f(x) = x f is an odd function, so the a n are zero, and thus the Fourier series will be of the form f(x)
More informationLecture 1 Some Time-Frequency Transformations
Lecture 1 Some Time-Frequency Transformations David Walnut Department of Mathematical Sciences George Mason University Fairfax, VA USA Chapman Lectures, Chapman University, Orange, CA 6-10 November 2017
More informationAdaptive Gabor transforms
Appl. Comput. Harmon. Anal. 13 (2002) 1 21 www.academicpress.com Adaptive Gabor transforms Ingrid Daubechies and Fabrice Planchon,1 Program in Applied and Computational Mathematics, Princeton University,
More informationWavelets and Affine Distributions A Time-Frequency Perspective
Wavelets and Affine Distributions A Time-Frequency Perspective Franz Hlawatsch Institute of Communications and Radio-Frequency Engineering Vienna University of Technology INSTITUT FÜR NACHRICHTENTECHNIK
More informationInstantaneous frequency computation: theory and practice
Instantaneous requency computation: theory and practice Matthew J. Yedlin 1, Gary F. Margrave 2, Yochai Ben Horin 3 1 Department o Electrical and Computer Engineering, UBC 2 Geoscience, University o Calgary
More informationFiltering in Time-Frequency Domain using STFrFT
Filtering in Time-Frequency Domain using STFrFT Pragati Rana P.G Student Vaibhav Mishra P.G Student Rahul Pachauri Sr. Lecturer. ABSTRACT The Fractional Fourier Transform is a generalized form of Fourier
More informationPre-Calculus: Functions and Their Properties (Solving equations algebraically and graphically, matching graphs, tables, and equations, and
Pre-Calculus: 1.1 1.2 Functions and Their Properties (Solving equations algebraically and graphically, matching graphs, tables, and equations, and finding the domain, range, VA, HA, etc.). Name: Date:
More informationIndex. Geostatistics for Environmental Scientists, 2nd Edition R. Webster and M. A. Oliver 2007 John Wiley & Sons, Ltd. ISBN:
Index Akaike information criterion (AIC) 105, 290 analysis of variance 35, 44, 127 132 angular transformation 22 anisotropy 59, 99 affine or geometric 59, 100 101 anisotropy ratio 101 exploring and displaying
More information3F1 Random Processes Examples Paper (for all 6 lectures)
3F Random Processes Examples Paper (for all 6 lectures). Three factories make the same electrical component. Factory A supplies half of the total number of components to the central depot, while factories
More informationLet x(t) be a finite energy signal with Fourier transform X(jω). Let E denote the energy of the signal. Thus E <. We have.
Notes on the Uncertainty principle Let x(t) be a finite energy signal with Fourier transform X(jω). Let E denote the energy of the signal. Thus E
More informationA6523 Linear, Shift-invariant Systems and Fourier Transforms
A6523 Linear, Shift-invariant Systems and Fourier Transforms Linear systems underly much of what happens in nature and are used in instrumentation to make measurements of various kinds. We will define
More informationRice method for the maximum of Gaussian fields
Rice method for the maximum of Gaussian fields IWAP, July 8, 2010 Jean-Marc AZAÏS Institut de Mathématiques, Université de Toulouse Jean-Marc AZAÏS ( Institut de Mathématiques, Université de Toulouse )
More informationPath integrals and the classical approximation 1 D. E. Soper 2 University of Oregon 14 November 2011
Path integrals and the classical approximation D. E. Soper University of Oregon 4 November 0 I offer here some background for Sections.5 and.6 of J. J. Sakurai, Modern Quantum Mechanics. Introduction There
More informationRenormalization of microscopic Hamiltonians. Renormalization Group without Field Theory
Renormalization of microscopic Hamiltonians Renormalization Group without Field Theory Alberto Parola Università dell Insubria (Como - Italy) Renormalization Group Universality Only dimensionality and
More informationFourier Analysis Linear transformations and lters. 3. Fourier Analysis. Alex Sheremet. April 11, 2007
Stochastic processes review 3. Data Analysis Techniques in Oceanography OCP668 April, 27 Stochastic processes review Denition Fixed ζ = ζ : Function X (t) = X (t, ζ). Fixed t = t: Random Variable X (ζ)
More information221A Lecture Notes Convergence of Perturbation Theory
A Lecture Notes Convergence of Perturbation Theory Asymptotic Series An asymptotic series in a parameter ɛ of a function is given in a power series f(ɛ) = f n ɛ n () n=0 where the series actually does
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 information1 Unitary representations of the Virasoro algebra
Week 5 Reading material from the books Polchinski, Chapter 2, 15 Becker, Becker, Schwartz, Chapter 3 Ginspargs lectures, Chapters 3, 4 1 Unitary representations of the Virasoro algebra Now that we have
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 informationProbability and Statistics for Final Year Engineering Students
Probability and Statistics for Final Year Engineering Students By Yoni Nazarathy, Last Updated: May 24, 2011. Lecture 6p: Spectral Density, Passing Random Processes through LTI Systems, Filtering Terms
More informationNonlinear Optimization
Nonlinear Optimization (Com S 477/577 Notes) Yan-Bin Jia Nov 7, 2017 1 Introduction Given a single function f that depends on one or more independent variable, we want to find the values of those variables
More informationPhysics 227 Exam 2. Rutherford said that if you really understand something you should be able to explain it to your grandmother.
Physics 227 Exam 2 Rutherford said that if you really understand something you should be able to explain it to your grandmother. For each of the topics on the next two pages, write clear, concise, physical
More informationMultiresolution analysis
Multiresolution analysis Analisi multirisoluzione G. Menegaz gloria.menegaz@univr.it The Fourier kingdom CTFT Continuous time signals + jωt F( ω) = f( t) e dt + f() t = F( ω) e jωt dt The amplitude F(ω),
More informationAccuracy and transferability of GAP models for tungsten
Accuracy and transferability of GAP models for tungsten Wojciech J. Szlachta Albert P. Bartók Gábor Csányi Engineering Laboratory University of Cambridge 5 November 214 Motivation Number of atoms 1 1 2
More informationLecture Wigner-Ville Distributions
Introduction to Time-Frequency Analysis and Wavelet Transforms Prof. Arun K. Tangirala Department of Chemical Engineering Indian Institute of Technology, Madras Lecture - 6.1 Wigner-Ville Distributions
More informationDigital Image Processing Lectures 15 & 16
Lectures 15 & 16, Professor Department of Electrical and Computer Engineering Colorado State University CWT and Multi-Resolution Signal Analysis Wavelet transform offers multi-resolution by allowing for
More informationReliability Theory of Dynamic Loaded Structures (cont.) Calculation of Out-Crossing Frequencies Approximations to the Failure Probability.
Outline of Reliability Theory of Dynamic Loaded Structures (cont.) Calculation of Out-Crossing Frequencies Approximations to the Failure Probability. Poisson Approximation. Upper Bound Solution. Approximation
More informationLecture 3 Partial Differential Equations
Lecture 3 Partial Differential Equations Prof. Massimo Guidolin Prep Course in Investments August-September 2016 Plan of the lecture Motivation and generalities The heat equation and its applications in
More informationStochastic Spectral Approaches to Bayesian Inference
Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to
More informationKernel Method: Data Analysis with Positive Definite Kernels
Kernel Method: Data Analysis with Positive Definite Kernels 2. Positive Definite Kernel and Reproducing Kernel Hilbert Space Kenji Fukumizu The Institute of Statistical Mathematics. Graduate University
More informationReview (Probability & Linear Algebra)
Review (Probability & Linear Algebra) CE-725 : Statistical Pattern Recognition Sharif University of Technology Spring 2013 M. Soleymani Outline Axioms of probability theory Conditional probability, Joint
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 informationAdaptive Short-Time Fractional Fourier Transform Used in Time-Frequency Analysis
Adaptive Short-Time Fractional Fourier Transform Used in Time-Frequency Analysis 12 School of Electronics and Information,Yili Normal University, Yining, 830054, China E-mail: tianlin20110501@163.com In
More informationELEG 3143 Probability & Stochastic Process Ch. 6 Stochastic Process
Department of Electrical Engineering University of Arkansas ELEG 3143 Probability & Stochastic Process Ch. 6 Stochastic Process Dr. Jingxian Wu wuj@uark.edu OUTLINE 2 Definition of stochastic process (random
More informationPOLARIZATION SPECTROGRAM OF BIVARIATE SIGNALS
POLARIZATION SPECTROGRAM OF BIVARIATE SIGNALS Julien Flamant Pierre Chainais Nicolas Le Bihan Univ. Lille, CNRS, Centrale Lille, UMR 9189 - CRIStAL, 59000 Lille, France CNRS/GIPSA-Lab, Grenoble, France
More informationENSC327 Communications Systems 2: Fourier Representations. Jie Liang School of Engineering Science Simon Fraser University
ENSC327 Communications Systems 2: Fourier Representations Jie Liang School of Engineering Science Simon Fraser University 1 Outline Chap 2.1 2.5: Signal Classifications Fourier Transform Dirac Delta Function
More informationReview Guideline for Final
Review Guideline for Final Here is the outline of the required skills for the final exam. Please read it carefully and find some corresponding homework problems in the corresponding sections to practice.
More informationUNDERSTANDING ENGINEERING MATHEMATICS
UNDERSTANDING ENGINEERING MATHEMATICS JOHN BIRD WORKED SOLUTIONS TO EXERCISES 1 INTRODUCTION In Understanding Engineering Mathematic there are over 750 further problems arranged regularly throughout the
More information2A1H Time-Frequency Analysis II Bugs/queries to HT 2011 For hints and answers visit dwm/courses/2tf
Time-Frequency Analysis II (HT 20) 2AH 2AH Time-Frequency Analysis II Bugs/queries to david.murray@eng.ox.ac.uk HT 20 For hints and answers visit www.robots.ox.ac.uk/ dwm/courses/2tf David Murray. A periodic
More informationRadar Signal Intra-Pulse Feature Extraction Based on Improved Wavelet Transform Algorithm
Int. J. Communications, Network and System Sciences, 017, 10, 118-17 http://www.scirp.org/journal/ijcns ISSN Online: 1913-373 ISSN Print: 1913-3715 Radar Signal Intra-Pulse Feature Extraction Based on
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 14: The Power Function and Native Space Error Estimates Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu
More informationF. Hlawatsch and F. Auger (eds.), Time-Frequency Analysis: Concepts and Methods, London (UK): ISTE and Wiley, 2008, 440 pages ISBN: 9781848210332 http://www.iste.co.uk/index.php?p=a&action=view&id=62 Contents
More informationThe Laplace Transform
The Laplace Transform Introduction There are two common approaches to the developing and understanding the Laplace transform It can be viewed as a generalization of the CTFT to include some signals with
More informationStochastic representation of random positive-definite tensor-valued properties: application to 3D anisotropic permeability random fields
Sous la co-tutelle de : LABORATOIRE DE MODÉLISATION ET SIMULATION MULTI ÉCHELLE CNRS UPEC UNIVERSITÉ PARIS-EST CRÉTEIL UPEM UNIVERSITÉ PARIS-EST MARNE-LA-VALLÉE Stochastic representation of random positive-definite
More informationEAS 305 Random Processes Viewgraph 1 of 10. Random Processes
EAS 305 Random Processes Viewgraph 1 of 10 Definitions: Random Processes A random process is a family of random variables indexed by a parameter t T, where T is called the index set λ i Experiment outcome
More informationMultidimensional partitions of unity and Gaussian terrains
and Gaussian terrains Richard A. Bale, Jeff P. Grossman, Gary F. Margrave, and Michael P. Lamoureu ABSTRACT Partitions of unity play an important rôle as amplitude-preserving windows for nonstationary
More informationMaxima and Minima. (a, b) of R if
Maxima and Minima Definition Let R be any region on the xy-plane, a function f (x, y) attains its absolute or global, maximum value M on R at the point (a, b) of R if (i) f (x, y) M for all points (x,
More informationRaktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Norms for Signals and Systems
. AERO 632: Design of Advance Flight Control System Norms for Signals and. Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. Norms for Signals ...
More informationStochastic Process II Dr.-Ing. Sudchai Boonto
Dr-Ing Sudchai Boonto Department of Control System and Instrumentation Engineering King Mongkuts Unniversity of Technology Thonburi Thailand Random process Consider a random experiment specified by the
More informationCompensation. June 29, Arizona State University. Edge Detection from Nonuniform Data via. Convolutional Gridding with Density.
Arizona State University June 29, 2012 Overview gridding with density compensation is a common method for function reconstruction nonuniform Fourier data. It is computationally inexpensive compared with
More information2. As we shall see, we choose to write in terms of σ x because ( X ) 2 = σ 2 x.
Section 5.1 Simple One-Dimensional Problems: The Free Particle Page 9 The Free Particle Gaussian Wave Packets The Gaussian wave packet initial state is one of the few states for which both the { x } and
More informationSOLUTIONS for Homework #2. 1. The given wave function can be normalized to the total probability equal to 1, ψ(x) = Ne λ x.
SOLUTIONS for Homework #. The given wave function can be normalized to the total probability equal to, ψ(x) = Ne λ x. () To get we choose dx ψ(x) = N dx e λx =, () 0 N = λ. (3) This state corresponds to
More informationSingle Variable Calculus, Early Transcendentals
Single Variable Calculus, Early Transcendentals 978-1-63545-100-9 To learn more about all our offerings Visit Knewtonalta.com Source Author(s) (Text or Video) Title(s) Link (where applicable) OpenStax
More informationFourier Series. Fourier Transform
Math Methods I Lia Vas Fourier Series. Fourier ransform Fourier Series. Recall that a function differentiable any number of times at x = a can be represented as a power series n= a n (x a) n where the
More information2 2 + x =
Lecture 30: Power series A Power Series is a series of the form c n = c 0 + c 1 x + c x + c 3 x 3 +... where x is a variable, the c n s are constants called the coefficients of the series. n = 1 + x +
More informationCommunications and Signal Processing Spring 2017 MSE Exam
Communications and Signal Processing Spring 2017 MSE Exam Please obtain your Test ID from the following table. You must write your Test ID and name on each of the pages of this exam. A page with missing
More information