ANALYSIS OF MRS SIGNALS WITH METABOLITE-BASED AUTOCORRELATION WAVELETS
|
|
- Camron Neal
- 5 years ago
- Views:
Transcription
1 ANALYSIS OF MRS SIGNALS WITH METABOLITE-BASED AUTOCORRELATION WAVELETS A. Schuck Jr. A. Suvichakorn C. Lemke J.-P. Antoine Unité de Physique Théorique et Mathematique, FYMA Université Catholique de Louvain, UCL Louvain-la-Neuve, BELGIUM Wavelets and Fractals, Esneux, April 26-28, 2 BELGIUM Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets /37
2 The traditional analysis The Autocorrelation Function Objectives and outline THE ANALYSIS SO FAR Time-domain and frequency-domain are traditional methods to quantify metabolites; Also Time-Frequency (Scale) methods have been used to analyzing MRS signals; The CWT using Morlet was presented in [4]. An example of pure Creatine using the Morlet wavelet: Fig. Freq. response of CRE at 4.7 T and its Morlet CWT (w = rad/s, σ =, Fs = 46, 4 s ). Extracted from A.-Suvichakorn and J.-P. Antoine, The Continuous Wavelet Transform in MRS, e-book FAST website, 29 Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 2/37
3 The traditional analysis The Autocorrelation Function Objectives and outline QUESTIONS MADE BY THE AUTHORS Questions made by the authors: Can one make a Wavelet function with its spectrum "Matched" with some MRS metabolite signal? 2 Can one estimate this metabolite s parameters using the CWT with such a wavelet function? Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 3/37
4 The traditional analysis The Autocorrelation Function Objectives and outline THE AUTOCORRELATION FUNCTION One answer can be the use of the concept of Autocorrelation Functions. The autocorrelation function estimator R xx (τ) of an ergodic process time series x(t) is ([3]): R xx (t) = x(τ)x(τ t) dt = where x(t) means the complex conjugate of x(t). x(τ)x(τ + t) dτ, () In the frequency domain, the Fourier transform of R xx, called S xx (ω) can be evaluated by the Wiener-Khintchine relation: S xx (ω) = F {R xx (τ)} = X(ω) 2. (2) Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 4/37
5 OBJECTIVES AND OUTLINE The traditional analysis The Autocorrelation Function Objectives and outline Objectives: Create a wavelet from a MRS pure metabolite signal autocorrelation function; 2 Perform the CWT of a complex MRS signal using this wavelet; 3 Estimate the metabolite parameters. Outline: Find an Analytical Solution for CWT using classic MRS models and autocorrelation wavelets form this model; 2 Create the discrete versions of signals and wavelets presented in the analytical part and analyze them with YAWtb Matlab Toolbox(Lorentzian Models); 3 Create discrete versions of signals and wavelets now based on the metabolite database and analyze them with Matlab (Metabolite-based Models); Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 5/37
6 One peak model Generalizing for N-peaks What happens with more realistic signals? LORENTZIAN LINESHAPE A simple classical FID MRS signal model x(t): The Lorentzian lineshape. x(t) = A e Dt e i(ωs t+φ) A >, D >. (3) In the frequency domain this signal is defined by X(ω) = 2πA e iφ δ(ω (ω s + id )), (4) where δ is the Dirac delta function. Fig. Real part (blue line) and imaginary part (red dash) of x(t) for A =, D = s, ω s = 32 rad/s, φ = rad. What is the autocorrelation function of this signal? NOT a CLUE!!! Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 6/37
7 MODIFIED LORENTZIAN LINESHAPE One peak model Generalizing for N-peaks What happens with more realistic signals? A slightly modified FID MRS signal model x (t) can be: x (t) = A e D t e iωs t θ(t), D >, (5) where θ(t) is the Heaviside function (or step function). In the frequency domain this signal is defined by X (ω) = A [D + i(ω ω s )]. (6) Now, the autocorrelation function (on freq. domain) is easily computed: S xx (ω) = A D + i(ω ω s ) 2 = A 2 D 2 + (ω ω s) 2. (7) Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 7/37
8 One peak model Generalizing for N-peaks What happens with more realistic signals? THE ADMISSIBLE WAVELET FUNCTION Is this function admissible? PROBABLY NOT. What to do, then? As with Morlet, use a correction term. Then the admissible wavelet function Ψ adm (ω) becomes: Ψ adm (ω) = Is this term really necessary? A 2 D 2 + (ω ω s ) 2 A 2 D 2 + (ω2 + ω 2 s ). (8) As with Morlet, in practice the appropriate choice of D and ω s make its value numerically negligible, so that the correcting term can indeed be omitted. Ψ(ω) = A 2 D 2 + (ω ω s ) 2. (9) Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 8/37
9 THE WAVELET FUNCTION: One peak model Generalizing for N-peaks What happens with more realistic signals? Fig. Ψ adm (ω) (blue line) and Ψ(ω) (red dot) for A =, D = and ωs = 32 rad/s. Fig. ψ(t), for A =, D = and ωs = 32 rad/s: real part (blue line) and imaginary part (red dash). Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 9/37
10 One peak model Generalizing for N-peaks What happens with more realistic signals? THE CWT OF A LORENTZIAN MRS SIGNAL Then, the CWT of the the Lorentzian lineshape signal x(t) using Ψ(ω) is: S(b, a) = a 2πA e i φ A 2 δ(ω (ω s + i D )) 2π D 2 + (aω ω s) 2 eiωb dω = a x(b) A 2 D 2 + [ωs (a ) + i a D ] 2. () S(b, a) diverges for a = cause signal and wavelet have the same damping factors" D. It cannot be used to estimate parameters, as in the Morlet case, but... Using different damping factors" D for the wavelet function and D for the signal, then S(b, a), for a = will become: S(b, ) = x(b) A 2 D 2 D2 = A e D b So the signal s damping factor can be calculated by: A 2 D 2 D2. () D = d ln S(b, ) (2) db Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets /37
11 THE (BEAUTIFUL) GRAPHICS One peak model Generalizing for N-peaks What happens with more realistic signals? Fig. S(b, a) for A =, D =, ω s = 32 rad/s,b = [, 2] s and a = [.5,.5], a. Fig. arg S(b, a) for A =, D =, ω s = 32 rad/s,b = [, 2] s and a = [.5,.5], a. Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets /37
12 One peak model Generalizing for N-peaks What happens with more realistic signals? GENERALIZING FOR N-PEAKS SIGNALS As many metabolites are multipeaked, we will generalize for N-peaks Lorentzian signals. The signal x N (t), a weighted sum of N Lorentzian components: N x N (t) = A ne Dnt e i(ωnt+φn) ; D n >. (3) The Fourier transform of x N (t): n= X N (ω) = 2π N A ne iφn δ(ω (ω sn + id n)). (4) n= The left truncated version of x N (t): x N (t) = n= where θ(t) is the Heaviside function. Now the Fourier transform of x N (t) N A ne Dnt e i(ωsnt) θ(t), D n >, (5) X N (ω) = N n= A n [D n + i(ω ω sn)]. (6) Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 2/37
13 One peak model Generalizing for N-peaks What happens with more realistic signals? THE N-PEAK WAVELET FUNCTION AND CWT The autocorrelation of (5) in the frequency domain is: S x N x N 2 N A n N (ω) = D n + i(ω ω n= sn) n= A 2 n D 2 n + (ω ω sn) 2. (7) Under the considerations made, the autocorrelation wavelet function is: Ψ N (ω) = N n= A 2 n D 2 n + (ω ω sn) 2 (8) The CWT of the signal (3) using (8) as wavelet function is given by: S N (b, a) = a N N A k e iφ k δ(ω (ω sk + id k )) k= n= A 2 n D 2 n + (aω ω sn ) 2 eiωb dω = N N A 2 n a x k (b) D 2 (9) k= n= n + [a(ω sk + id k ) ω sn] 2 Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 3/37
14 One peak model Generalizing for N-peaks What happens with more realistic signals? (ALSO BEAUTIFUL)GRAPHS OF 2 AND 3 PEAKS CWT Fig. S 2 (b, a) for A = A 2 =, D = D2 =, ω s = 32 and ω s2 = 64 rad/s, b = [, 2] s and a = [., 3], a.; Fig. S 3 (b, a) for A = A 2 = A 3 =, D = D 2 = D 3 =, ω s = 3, ω s2 = 6 and ω s2 = 9 rad/s, b = [, 2] s and a = [., 4], a.; Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 4/37
15 One peak model Generalizing for N-peaks What happens with more realistic signals? THE N-PEAK DAMPING FACTORS ESTIMATION An oscillating" local maxima at a = ; Also, S N (b, a) have horizontal ridges at a = ω sn /ω sk, k n; k, n =, 2,..., N. Estimation of the damping factors: Consider the wavelet function have damping factors D n and the signal have D nn, D n D nn for any n. 2 Choosing one of the other local maxima a = ω sn /ω sk, and considering that the other factors are small enough in this scale, so only one term of (9) will be significative and we can estimate the damping factor as: d db ln S(b, ωsn ω sk ) d db ln e D kk b + d db ln A k A 2 n d db ln D2 n + [(ωs k + id kk ) ω sn ] 2 D kk. (2) Summarizing: if (a) peaks are far enough from each other and;(b) damping factors provide sharp peaks, the damping factors can estimated by D kk d ωsn ln S(b, db ω sk ). Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 5/37
16 One peak model Generalizing for N-peaks What happens with more realistic signals? WHAT HAPPENS WITH MORE REALISTIC SIGNALS? What happens when one analyzes limited, discretized, noisy and more realistic signals? Two procedures proposed: Create Discrete versions of signals and wavelets presented before and analyze them with Matlab and YAWtb Toolbox (Lorentzian Models); 2 Finally, create Discrete versions of signals and wavelets based on the metabolite database and analyze them with Matlab (Metabolite-based Models); Next, in the Numerical Analysis section. Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 6/37
17 NUMERICAL ANALYSIS WITH PEAK LORENTZIAN FUNCTION A MATLAB c function which implements Ψ(ω) = N A n n= D n+i(ω ω sn ) 2 were created (called LorentzNd.m") and added to YAWtb toolbox ([5]); A N component discrete exponential signal was defined by:, n ( N 4 ), x[n] = N n= A n e D (n N 4 )ts e iωsn (n N 4 )ts, ( N 4 ) n ( 3N 4 ),, ( 3N 4 ) n N, where t s = /f s is the sampling period in seconds. The CWT of this signal using the LorentzNd" function, with the same frequency and damping factors as the signals, for n =, 2, 3 was performed; Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 7/37
18 THE CWT RESULTS FOR COMPONENT Fig. CWT x [n] for A =, D =, ω s = 32 rad/s, n = [, 496] and a = [.5,.5] using the Lorentzd" wavelet. Fig. Skeleton of the CWT, for A =, D =, ω s = 32 rad/s, n = [, 496] and a = [.5,.5] using the Lorentzd" wavelet. Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 8/37
19 THE CWT RESULTS FOR 2 COMPONENTS Fig. CWT x 2 [n] for A = A 2 =, D = D 2 = /s, ω s = 32 and ω s2 = 64 rad/s, t s = /256 s and N = 496 for Lorentz2Pkd" wavelet. Fig. Skeleton of CWT for A = A 2 =, D = D 2 = /s, ω s = 32 and ω s2 = 64 rad/s, t s = /256 s and N = 496 for Lorentz2Pkd" wavelet. Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 9/37
20 ...AND THE CWT RESULTS FOR 3 COMPONENTS: Fig. CWT x 3 [n] for A = A 2 = A 3 =, D = D 2 = D 3 = /s, ω s = 3, ω s2 = 6 and ω s2 = 9 rad/s, t s = /256 s, N = 496 and a = [.5, 5] for Lorentz3Pkd" wavelet. Fig. Skeleton of CWT for A = A 2 = A 3 =, D = D 2 = D 3 = /s, ω s = 3, ω s2 = 6 and ω s2 = 9 rad/s, t s = /256 s, N = 496 and a = [.5, 5] for Lorentz3Pkd" wavelet. Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 2/37
21 NUMERICAL ANALYSIS FOR SIGNALS WITH NOISE Here the goal was to find the limit of signal detection in presence of noise by means of autocorrelation wavelets ; White gaussian noise was added to signals at different SNR levels; 2 The CWT was performed (N peak was analyzed by its related N autocorrelation wavelet; 3 The results are in the limit of signal detection by its wavelet (still have the horizontal ridge); 4 The SNR of this limit were in average 2dB Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 2/37
22 SOME NOISE RESULTS Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 22/37
23 CREATING AND USING METABOLITE-BASED WAVELETS Numerical analysis with Metabolite-based wavelets: Take signals from one in vitro metabolite database; 2 Create its related Autocorrelation Wavelets; 3 Analyze a mixture of metabolites by one related autocorrelation wavelet with Matlab; Metabolite profile (Cre) 3 Cre The Algorithm: 25 Calculate autocorrelation function R of metabolite profile φ: Φ Cre (f) 2 5 R[n] = +N k= N φ[k] φ[n k] (2) f [khz].5 using Matlab function xcorr; 4 2 Subtract its mean value: ψ[n] = R[n] E{R[n]} (22) Re{φ Cre (t)} t [ms] Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 23/37
24 WAVELET CONSTRUCTION RESULT Metabolite profile (Cre) 3 Cre Cre-based wavelet x 4 Cre 25 8 Φ Cre (f) Ψ Cre (f) f [khz] f [khz] Re{φ Cre (t)} 2 Re{ψ Cre (t)} t [ms] t [ms] Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 24/37
25 HOW TO DILATE DISCRETE WAVELETS? PROBLEM: Mother wavelet is Discrete, no analytical expression; SOLUTION: An discrete upsampler/downsampler system was used. Fig. Block diagram of Upsampler/downsampler system. L, M (Z ). Wavelet will be expanded by an integer factor of L and contracted by an factor of M. a = L/M. Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 25/37
26 EXAMPLE OF DILATED WAVELET Cre-based wavelet (real part) at different :.5 scale a=.5 scale a=.2 Re{ψ Cre (t)}.5 Re{ψ Cre (t)} t [ms] 5 5 t [ms] Fig. scale a = Fig. scale a =.2 Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 26/37
27 SIGNAL : COMBINATION NAA + CRE WITHOUT LAC Sum of pure Naa and Cre signals, no noise: 2 x 7 Naa 2 x 7 Cre 2 x Y Naa (f) 6 Y Cre (f) 6 S(f) f [khz].5.5 f [khz].5.5 f [khz].5 Fig. Naa part Fig. Cre part Fig. Naa + Cre Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 27/37
28 SIGNAL 2: COMBINATION NAA + CRE + LAC Sum of pure Naa and Cre signals PLUS Lac, no noise: 2 x 7 2 x 7 Lac 2 x S(f) 6 Y Lac (f) 6 S(f) f [khz].5.5 f [khz].5.5 f [khz].5 Fig. Naa + Cre Fig. Lac part Fig. Naa + Cre + Lac Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 28/37
29 CWT ANALYSIS OF SIGNAL WITH NAA-BASED WAVELET CWT Naa: Max= 75.2 at scale=.99 and τ=47ms CWT Naa: Max= at scale= and τ=27ms horizontal ridges, thresh= horizontal ridges, thresh= Fig. Naa-wavelet CWT of Naa reference signal Fig. Naa-wavelet CWT of Naa+Cre composed signal Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 29/37
30 CWT ANALYSIS OF SIGNAL 2 WITH NAA-BASED WAVELET CWT Naa: Max= 75.2 at scale=.99 and τ=47ms CWT Naa: Max= at scale= and τ=27ms horizontal ridges, thresh= horizontal ridges, thresh= Fig. Naa-wavelet CWT of Naa reference signal Fig. Naa-wavelet CWT of Naa+Cre+Lac composed signal Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 3/37
31 CWT ANALYSIS OF SIGNAL WITH CRE-BASED WAVELET CWT Cre: Max= 63.4 at scale=.98 and τ=52ms CWT Cre: Max= at scale= and τ=6ms horizontal ridges, thresh= horizontal ridges, thresh= Fig. Cre-wavelet CWT of Cre reference signal Fig. Cre-wavelet CWT of Naa+Cre composed signal Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 3/37
32 CWT ANALYSIS OF SIGNAL 2 WITH CRE-BASED WAVELET CWT Cre: Max= 63.4 at scale=.98 and τ=52ms CWT Cre: Max= at scale= and τ=6ms horizontal ridges, thresh= horizontal ridges, thresh= Fig. Cre-wavelet CWT of Cre reference signal Fig. Cre-wavelet CWT of Naa+Cre+Lac composed signal Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 32/37
33 CWT ANALYSIS OF SIGNAL WITH LAC-BASED WAVELET CWT Lac: Max= 64.7 at scale=.99 and τ=5ms CWT Lac: Max= at scale=.2 and τ=6ms horizontal ridges, thresh= horizontal ridges, thresh= Fig. Lac-wavelet CWT of Lac reference signal Fig. Lac-wavelet CWT of Naa+Cre composed signal Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 33/37
34 CWT ANALYSIS OF SIGNAL 2 WITH LAC-BASED WAVELET CWT Lac: Max= 64.7 at scale=.99 and τ=5ms CWT Lac: Max= at scale=.2 and τ=3ms horizontal ridges, thresh= horizontal ridges, thresh= Fig. Lac-wavelet CWT of Lac reference signal Fig. Lac-wavelet CWT of Naa+Cre+Lac composed signal Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 34/37
35 THE CONCLUSIONS SO FAR Analytical analysis: Analytical expressions for Wavelet function and CWT, using Lorentzian models were made; 2 Horizontal ridges at a = means Presence of Metabolite which generated the autocorr.wavelet; 3 For peak, the Damping factor can be found. For more than one peak, the Damping factors can be approximated from horizontal ridges at a ; Numerical analysis with Lorentzian signals: Discretization and time limitation of signals changed a little the results ( lateral" Ridges are not smooth); 2 With multi-peaked signal/autocorr. wavelet, the damping estimation is harder; 3 Single signals could be detected at 2dB SNR; Numerical analysis with metabolite-based signals: Algorithms for autocorrelation wavelet creation, dilation and CWT were made; 2 The metabolites presence could be detected in the mixture by its related wavelet; 3 Parameters estimation is almost done (so wait a little bit more). Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 35/37
36 ACKNOWLEDGEMENTS The authors would like to thank Diana Sima, Maria Isabel Osorio Garcia and Sabine van Huffel from KUL; Adalberto Schuck Jr. would like to thank Conselho Nacional de Pesquisa e Desenvolvimento - CNPq, Brasil, for its financial support. Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 36/37
37 REFERENCES Introduction J-P. Antoine and R. Murenzi and P. Vandergheynst and S.T. Ali. Two-Dimensional Wavelets and their Relatives. Cambridge University Press, 24. A. Oppenheim and R.W. Schaffer with J. R. Buck. Discrete-time signal processing, 2nd Ed.. Prentice-Hall, Inc., 999. A. Papoulis and S. U. Pillai. Probability, random variables and stochastic processes, 4th ed.. McGraw-Hill, Inc., 22. A. Suvichakorn and H. Ratiney and A. Bucur and S. Cavassila and J-P. Antoine. Toward a quantitative analysis of in vivo proton magnetic resonance spectroscopic signals using the continuous Morlet wavelet transform. Meas. Sci. Technol., 2:429 44, 29. L. Jacques and A. Coron and P. Vandergheynst and A. Rivoldini. YAWTb toolbox : Yet Another Wavelet Toolbox Adalberto Schuck Jr. Analysis of MRS Signals with Autocorr. Wavelets 37/37
Joint ICTP-TWAS School on Coherent State Transforms, Time- Frequency and Time-Scale Analysis, Applications.
585- Joint ICTP-TWAS School on Coherent State Transforms, Time- Frequency and Time-Scale Analysis, Applications - June Wavelets and multiresolution : from NMR spectroscopy to motion analysis J-P. Antoine
More informationarxiv: v1 [math.ca] 6 Feb 2015
The Fourier-Like and Hartley-Like Wavelet Analysis Based on Hilbert Transforms L. R. Soares H. M. de Oliveira R. J. Cintra Abstract arxiv:150.0049v1 [math.ca] 6 Feb 015 In continuous-time wavelet analysis,
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 informationFigure 3.1 Effect on frequency spectrum of increasing period T 0. Consider the amplitude spectrum of a periodic waveform as shown in Figure 3.2.
3. Fourier ransorm From Fourier Series to Fourier ransorm [, 2] In communication systems, we oten deal with non-periodic signals. An extension o the time-requency relationship to a non-periodic signal
More informationCorrelator I. Basics. Chapter Introduction. 8.2 Digitization Sampling. D. Anish Roshi
Chapter 8 Correlator I. Basics D. Anish Roshi 8.1 Introduction A radio interferometer measures the mutual coherence function of the electric field due to a given source brightness distribution in the sky.
More informationCast of Characters. Some Symbols, Functions, and Variables Used in the Book
Page 1 of 6 Cast of Characters Some s, Functions, and Variables Used in the Book Digital Signal Processing and the Microcontroller by Dale Grover and John R. Deller ISBN 0-13-081348-6 Prentice Hall, 1998
More informationRandom signals II. ÚPGM FIT VUT Brno,
Random signals II. Jan Černocký ÚPGM FIT VUT Brno, cernocky@fit.vutbr.cz 1 Temporal estimate of autocorrelation coefficients for ergodic discrete-time random process. ˆR[k] = 1 N N 1 n=0 x[n]x[n + k],
More informationEEM 409. Random Signals. Problem Set-2: (Power Spectral Density, LTI Systems with Random Inputs) Problem 1: Problem 2:
EEM 409 Random Signals Problem Set-2: (Power Spectral Density, LTI Systems with Random Inputs) Problem 1: Consider a random process of the form = + Problem 2: X(t) = b cos(2π t + ), where b is a constant,
More informationFilter Banks For "Intensity Analysis"
morlet.sdw 4/6/3 /3 Filter Banks For "Intensity Analysis" Introduction In connection with our participation in two conferences in Canada (August 22 we met von Tscharner several times discussing his "intensity
More informationLecture 2: Haar Multiresolution analysis
WAVELES AND MULIRAE DIGIAL SIGNAL PROCESSING Lecture 2: Haar Multiresolution analysis Prof.V. M. Gadre, EE, II Bombay 1 Introduction HAAR was a mathematician, who has given an idea that any continuous
More informationOPTIMIZATION OF MORLET WAVELET FOR MECHANICAL FAULT DIAGNOSIS
Twelfth International Congress on Sound and Vibration OPTIMIZATION OF MORLET WAVELET FOR MECHANICAL FAULT DIAGNOSIS Jiří Vass* and Cristina Cristalli** * Dept. of Circuit Theory, Faculty of Electrical
More informationWavelets and multiresolution representations. Time meets frequency
Wavelets and multiresolution representations Time meets frequency Time-Frequency resolution Depends on the time-frequency spread of the wavelet atoms Assuming that ψ is centred in t=0 Signal domain + t
More informationSpectral Analysis of Random Processes
Spectral Analysis of Random Processes Spectral Analysis of Random Processes Generally, all properties of a random process should be defined by averaging over the ensemble of realizations. Generally, all
More informationWavelet entropy as a measure of solar cycle complexity
Astron. Astrophys. 363, 3 35 () Wavelet entropy as a measure of solar cycle complexity S. Sello Mathematical and Physical Models, Enel Research, Via Andrea Pisano, 56 Pisa, Italy (sello@pte.enel.it) Received
More informationPARAMETERS OPTIMIZATION OF CONTINUOUS WAVELET TRANSFORM AND ITS APPLICATION IN ACOUSTIC EMISSION SIGNAL ANALYSIS OF ROLLING BEARING *
104 CHINESE JOURNAL OF MECHANICAL ENGINEERING Vol. 20,aNo. 2,a2007 ZHANG Xinming HE Yongyong HAO Ruiang CHU Fulei State Key Laboratory of Tribology, Tsinghua University, Being 100084, China PARAMETERS
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 informationSimulation of the Effects of Timing Jitter in Track-and-Hold and Sample-and-Hold Circuits
5. Simulation of the Effects of Timing Jitter in Track-and-Hold and Sample-and-Hold Circuits V. Vasudevan Department of Electrical Engineering Indian Institute of Technology-Madras Chennai-636, India Email:
More informationSignal interactions Cross correlation, cross spectral coupling and significance testing Centre for Doctoral Training in Healthcare Innovation
Signal interactions Cross correlation, cross spectral coupling and significance testing Centre for Doctoral Training in Healthcare Innovation Dr. Gari D. Clifford, University Lecturer & Director, Centre
More informationPhoton Physics. Week 4 26/02/2013
Photon Physics Week 4 6//13 1 Classical atom-field interaction Lorentz oscillator: Classical electron oscillator with frequency ω and damping constant γ Eqn of motion: Final result: Classical atom-field
More information5.74 Introductory Quantum Mechanics II
MIT OpenCourseWare http://ocw.mit.edu 5.74 Introductory Quantum Mechanics II Spring 9 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Andrei Tokmakoff,
More informationA.1 THE SAMPLED TIME DOMAIN AND THE Z TRANSFORM. 0 δ(t)dt = 1, (A.1) δ(t)dt =
APPENDIX A THE Z TRANSFORM One of the most useful techniques in engineering or scientific analysis is transforming a problem from the time domain to the frequency domain ( 3). Using a Fourier or Laplace
More information1. Fourier Transform (Continuous time) A finite energy signal is a signal f(t) for which. f(t) 2 dt < Scalar product: f(t)g(t)dt
1. Fourier Transform (Continuous time) 1.1. Signals with finite energy A finite energy signal is a signal f(t) for which Scalar product: f(t) 2 dt < f(t), g(t) = 1 2π f(t)g(t)dt The Hilbert space of all
More informationSpectra (2A) Young Won Lim 11/8/12
Spectra (A) /8/ Copyright (c) 0 Young W. Lim. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version. or any later version
More informationHEATED WIND PARTICLE S BEHAVIOURAL STUDY BY THE CONTINUOUS WAVELET TRANSFORM AS WELL AS THE FRACTAL ANALYSIS
HEATED WIND PARTICLE S BEHAVIOURAL STUDY BY THE CONTINUOUS WAVELET TRANSFORM AS WELL AS THE FRACTAL ANALYSIS ABSTRACT Sabyasachi Mukhopadhyay 1, Sanmoy Mandal,, P.K.Panigrahi 3, Asish Mitra 4 1,,3 Physical
More informationA6523 Modeling, Inference, and Mining Jim Cordes, Cornell University. False Positives in Fourier Spectra. For N = DFT length: Lecture 5 Reading
A6523 Modeling, Inference, and Mining Jim Cordes, Cornell University Lecture 5 Reading Notes on web page Stochas
More information3 Chemical exchange and the McConnell Equations
3 Chemical exchange and the McConnell Equations NMR is a technique which is well suited to study dynamic processes, such as the rates of chemical reactions. The time window which can be investigated in
More informationPhys460.nb Back to our example. on the same quantum state. i.e., if we have initial condition (5.241) ψ(t = 0) = χ n (t = 0)
Phys46.nb 89 on the same quantum state. i.e., if we have initial condition ψ(t ) χ n (t ) (5.41) then at later time ψ(t) e i ϕ(t) χ n (t) (5.4) This phase ϕ contains two parts ϕ(t) - E n(t) t + ϕ B (t)
More informationMedical Image Processing Using Transforms
Medical Image Processing Using Transforms Hongmei Zhu, Ph.D Department of Mathematics & Statistics York University hmzhu@yorku.ca MRcenter.ca Outline Image Quality Gray value transforms Histogram processing
More informationChapter Three Theoretical Description Of Stochastic Resonance 24
Table of Contents List of Abbreviations and Symbols 5 Chapter One Introduction 8 1.1 The Phenomenon of the Stochastic Resonance 8 1.2 The Purpose of the Study 10 Chapter Two The Experimental Set-up 12
More informationWavelets, wavelet networks and the conformal group
Wavelets, wavelet networks and the conformal group R. Vilela Mendes CMAF, University of Lisbon http://label2.ist.utl.pt/vilela/ April 2016 () April 2016 1 / 32 Contents Wavelets: Continuous and discrete
More informationMassachusetts Institute of Technology
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.011: Introduction to Communication, Control and Signal Processing QUIZ, April 1, 010 QUESTION BOOKLET Your
More informationPhotonic Excess Noise in a p-i-n Photodetector 1
Photonic Excess Noise in a p-i-n Photodetector 1 F. Javier Fraile-Peláez, David J. Santos and José Capmany )+( Dept. de Teoría de la Señal y Comunicaciones, Universidad de Vigo ETSI Telecomunicación, Campus
More informationDIGITAL SIGNAL PROCESSING LECTURE 1
DIGITAL SIGNAL PROCESSING LECTURE 1 Fall 2010 2K8-5 th Semester Tahir Muhammad tmuhammad_07@yahoo.com Content and Figures are from Discrete-Time Signal Processing, 2e by Oppenheim, Shafer, and Buck, 1999-2000
More informationTopic 3: Fourier Series (FS)
ELEC264: Signals And Systems Topic 3: Fourier Series (FS) o o o o Introduction to frequency analysis of signals CT FS Fourier series of CT periodic signals Signal Symmetry and CT Fourier Series Properties
More informationDirectional Wavelets Revisited: Cauchy Wavelets and Symmetry Detection in Patterns
Applied and Computational Harmonic Analysis 6, 314 345 (1999) Article ID acha.1998.0255, available online at http://www.idealibrary.com on Directional Wavelets Revisited: Cauchy Wavelets and Symmetry Detection
More informationReview of Fourier Transform
Review of Fourier Transform Fourier series works for periodic signals only. What s about aperiodic signals? This is very large & important class of signals Aperiodic signal can be considered as periodic
More informationLinear and Nonlinear Oscillators (Lecture 2)
Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical
More informationEE531 (Semester II, 2010) 6. Spectral analysis. power spectral density. periodogram analysis. window functions 6-1
6. Spectral analysis EE531 (Semester II, 2010) power spectral density periodogram analysis window functions 6-1 Wiener-Khinchin theorem: Power Spectral density if a process is wide-sense stationary, the
More informationQuestion: Total. Points:
MATH 308 May 23, 2011 Final Exam Name: ID: Question: 1 2 3 4 5 6 7 8 9 Total Points: 0 20 20 20 20 20 20 20 20 160 Score: There are 9 problems on 9 pages in this exam (not counting the cover sheet). Make
More informationWavelet Transform. Figure 1: Non stationary signal f(t) = sin(100 t 2 ).
Wavelet Transform Andreas Wichert Department of Informatics INESC-ID / IST - University of Lisboa Portugal andreas.wichert@tecnico.ulisboa.pt September 3, 0 Short Term Fourier Transform Signals whose frequency
More informationStochastic Processes. M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno
Stochastic Processes M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno 1 Outline Stochastic (random) processes. Autocorrelation. Crosscorrelation. Spectral density function.
More informationOrdinary Least Squares and its applications
Ordinary Least Squares and its applications Dr. Mauro Zucchelli University Of Verona December 5, 2016 Dr. Mauro Zucchelli Ordinary Least Squares and its applications December 5, 2016 1 / 48 Contents 1
More informationRecursive Gaussian filters
CWP-546 Recursive Gaussian filters Dave Hale Center for Wave Phenomena, Colorado School of Mines, Golden CO 80401, USA ABSTRACT Gaussian or Gaussian derivative filtering is in several ways optimal for
More informationSpectral Broadening Mechanisms
Spectral Broadening Mechanisms Lorentzian broadening (Homogeneous) Gaussian broadening (Inhomogeneous, Inertial) Doppler broadening (special case for gas phase) The Fourier Transform NC State University
More informationAnalyzing the Effect of Moving Resonance on Seismic Response of Structures Using Wavelet Transforms
Analyzing the Effect of Moving Resonance on Seismic Response of Structures Using Wavelet Transforms M.R. Eatherton Virginia Tech P. Naga WSP Cantor Seinuk, New York, NY SUMMARY: When the dominant natural
More information221B Lecture Notes on Resonances in Classical Mechanics
1B Lecture Notes on Resonances in Classical Mechanics 1 Harmonic Oscillators Harmonic oscillators appear in many different contexts in classical mechanics. Examples include: spring, pendulum (with a small
More informationCOMPRESSIVE OPTICAL DEFLECTOMETRIC TOMOGRAPHY
COMPRESSIVE OPTICAL DEFLECTOMETRIC TOMOGRAPHY Adriana González ICTEAM/UCL March 26th, 214 1 ISPGroup - ICTEAM - UCL http://sites.uclouvain.be/ispgroup Université catholique de Louvain, Louvain-la-Neuve,
More information08. Brownian Motion. University of Rhode Island. Gerhard Müller University of Rhode Island,
University of Rhode Island DigitalCommons@URI Nonequilibrium Statistical Physics Physics Course Materials 1-19-215 8. Brownian Motion Gerhard Müller University of Rhode Island, gmuller@uri.edu Follow this
More informationCOMPLEX WAVELET TRANSFORM IN SIGNAL AND IMAGE ANALYSIS
COMPLEX WAVELET TRANSFORM IN SIGNAL AND IMAGE ANALYSIS MUSOKO VICTOR, PROCHÁZKA ALEŠ Institute of Chemical Technology, Department of Computing and Control Engineering Technická 905, 66 8 Prague 6, Cech
More informationChemical Exchange. Spin-interactions External interactions Magnetic field Bo, RF field B1
Chemical Exchange Spin-interactions External interactions Magnetic field Bo, RF field B1 Internal Interactions Molecular motions Chemical shifts J-coupling Chemical Exchange 1 Outline Motional time scales
More informationEvaluation of the modified group delay feature for isolated word recognition
Evaluation of the modified group delay feature for isolated word recognition Author Alsteris, Leigh, Paliwal, Kuldip Published 25 Conference Title The 8th International Symposium on Signal Processing and
More informationLecture #6 Chemical Exchange
Lecture #6 Chemical Exchange Topics Introduction Effects on longitudinal magnetization Effects on transverse magnetization Examples Handouts and Reading assignments Kowalewski, Chapter 13 Levitt, sections
More informationEli Barkai. Wrocklaw (2015)
1/f Noise and the Low-Frequency Cutoff Paradox Eli Barkai Bar-Ilan University Niemann, Kantz, EB PRL (213) Theory Sadegh, EB, Krapf NJP (214) Experiment Stefani, Hoogenboom, EB Physics Today (29) Wrocklaw
More informationTime Response Analysis (Part II)
Time Response Analysis (Part II). A critically damped, continuous-time, second order system, when sampled, will have (in Z domain) (a) A simple pole (b) Double pole on real axis (c) Double pole on imaginary
More informationStochastic Processes. A stochastic process is a function of two variables:
Stochastic Processes Stochastic: from Greek stochastikos, proceeding by guesswork, literally, skillful in aiming. A stochastic process is simply a collection of random variables labelled by some parameter:
More information1 Introduction to Wavelet Analysis
Jim Lambers ENERGY 281 Spring Quarter 2007-08 Lecture 9 Notes 1 Introduction to Wavelet Analysis Wavelets were developed in the 80 s and 90 s as an alternative to Fourier analysis of signals. Some of the
More informationSimple Identification of Nonlinear Modal Parameters Using Wavelet Transform
Proceedings of the 9 th ISSM achen, 7 th -9 th October 4 1 Simple Identification of Nonlinear Modal Parameters Using Wavelet Transform Tegoeh Tjahjowidodo, Farid l-bender, Hendrik Van Brussel Mechanical
More informationHow many initial conditions are required to fully determine the general solution to a 2nd order linear differential equation?
How many initial conditions are required to fully determine the general solution to a 2nd order linear differential equation? (A) 0 (B) 1 (C) 2 (D) more than 2 (E) it depends or don t know How many of
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 informationDevelopment of PC-Based Leak Detection System Using Acoustic Emission Technique
Key Engineering Materials Online: 004-08-5 ISSN: 66-9795, Vols. 70-7, pp 55-50 doi:0.408/www.scientific.net/kem.70-7.55 004 Trans Tech Publications, Switzerland Citation & Copyright (to be inserted by
More informationROBUST VIRTUAL DYNAMIC STRAIN SENSORS FROM ACCELERATION MEASUREMENTS
7th European Workshop on Structural Health Monitoring July 8-11, 2014. La Cité, Nantes, France More Info at Open Access Database www.ndt.net/?id=17230 ROBUST VIRTUAL DYNAMIC STRAIN SENSORS FROM ACCELERATION
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 informationA6523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring
A653 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring 15 http://www.astro.cornell.edu/~cordes/a653 Lecture 3 Power spectrum issues Frequentist approach Bayesian approach (some
More information5.74 Introductory Quantum Mechanics II
MIT OpenCourseWare http://ocw.mit.edu 5.74 Introductory Quantum Mechanics II Spring 009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Andrei Tokmakoff,
More informationRadar Systems Engineering Lecture 3 Review of Signals, Systems and Digital Signal Processing
Radar Systems Engineering Lecture Review of Signals, Systems and Digital Signal Processing Dr. Robert M. O Donnell Guest Lecturer Radar Systems Course Review Signals, Systems & DSP // Block Diagram of
More information2A1H Time-Frequency Analysis II
2AH Time-Frequency Analysis II Bugs/queries to david.murray@eng.ox.ac.uk HT 209 For any corrections see the course page DW Murray at www.robots.ox.ac.uk/ dwm/courses/2tf. (a) A signal g(t) with period
More information742. Time-varying systems identification using continuous wavelet analysis of free decay response signals
74. Time-varying systems identification using continuous wavelet analysis of free decay response signals X. Xu, Z. Y. Shi, S. L. Long State Key Laboratory of Mechanics and Control of Mechanical Structures
More informationWavelet Methods for Time Series Analysis. What is a Wavelet? Part I: Introduction to Wavelets and Wavelet Transforms. sines & cosines are big waves
Wavelet Methods for Time Series Analysis Part I: Introduction to Wavelets and Wavelet Transforms wavelets are analysis tools for time series and images as a subject, wavelets are relatively new (1983 to
More informationLecture 15: Time and Frequency Joint Perspective
WAVELETS AND MULTIRATE DIGITAL SIGNAL PROCESSING Lecture 15: Time and Frequency Joint Perspective Prof.V.M.Gadre, EE, IIT Bombay Introduction In lecture 14, we studied steps required to design conjugate
More informationSIO 211B, Rudnick. We start with a definition of the Fourier transform! ĝ f of a time series! ( )
SIO B, Rudnick! XVIII.Wavelets The goal o a wavelet transorm is a description o a time series that is both requency and time selective. The wavelet transorm can be contrasted with the well-known and very
More informationReview of Discrete-Time System
Review of Discrete-Time System Electrical & Computer Engineering University of Maryland, College Park Acknowledgment: ENEE630 slides were based on class notes developed by Profs. K.J. Ray Liu and Min Wu.
More informationSignals and Spectra - Review
Signals and Spectra - Review SIGNALS DETERMINISTIC No uncertainty w.r.t. the value of a signal at any time Modeled by mathematical epressions RANDOM some degree of uncertainty before the signal occurs
More informationComputer Assignment (CA) No. 11: Autocorrelation And Power Spectral Density
Computer Assignment (CA) No. 11: Autocorrelation And Power Spectral Density Problem Statement Recall the autocorrelation function is defined as: N 1 R(τ) = n=0 x[n]x[n τ],τ = 0, 1, 2,..., M Compute and
More informationECE 3800 Probabilistic Methods of Signal and System Analysis
ECE 3800 Probabilistic Methods of Signal and System Analysis Dr. Bradley J. Bazuin Western Michigan University College of Engineering and Applied Sciences Department of Electrical and Computer Engineering
More informationEstimation of Variance and Skewness of Non-Gaussian Zero mean Color Noise from Measurements of the Atomic Transition Probabilities
International Journal of Electronic and Electrical Engineering. ISSN 974-2174, Volume 7, Number 4 (214), pp. 365-372 International Research Publication House http://www.irphouse.com Estimation of Variance
More information4. Complex Oscillations
4. Complex Oscillations The most common use of complex numbers in physics is for analyzing oscillations and waves. We will illustrate this with a simple but crucially important model, the damped harmonic
More informationFeature extraction 2
Centre for Vision Speech & Signal Processing University of Surrey, Guildford GU2 7XH. Feature extraction 2 Dr Philip Jackson Linear prediction Perceptual linear prediction Comparison of feature methods
More informationADAPTIVE ANTENNAS. SPATIAL BF
ADAPTIVE ANTENNAS SPATIAL BF 1 1-Spatial reference BF -Spatial reference beamforming may not use of embedded training sequences. Instead, the directions of arrival (DoA) of the impinging waves are used
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 informationScattering of Electromagnetic Radiation. References:
Scattering of Electromagnetic Radiation References: Plasma Diagnostics: Chapter by Kunze Methods of experimental physics, 9a, chapter by Alan Desilva and George Goldenbaum, Edited by Loveberg and Griem.
More informationSpectral Broadening Mechanisms. Broadening mechanisms. Lineshape functions. Spectral lifetime broadening
Spectral Broadening echanisms Lorentzian broadening (Homogeneous) Gaussian broadening (Inhomogeneous, Inertial) Doppler broadening (special case for gas phase) The Fourier Transform NC State University
More information1 The Continuous Wavelet Transform The continuous wavelet transform (CWT) Discretisation of the CWT... 2
Contents 1 The Continuous Wavelet Transform 1 1.1 The continuous wavelet transform (CWT)............. 1 1. Discretisation of the CWT...................... Stationary wavelet transform or redundant wavelet
More informationReview of Frequency Domain Fourier Series: Continuous periodic frequency components
Today we will review: Review of Frequency Domain Fourier series why we use it trig form & exponential form how to get coefficients for each form Eigenfunctions what they are how they relate to LTI systems
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 informationSignal Analysis. David Ozog. May 11, Abstract
Signal Analysis David Ozog May 11, 2007 Abstract Signal processing is the analysis, interpretation, and manipulation of any time varying quantity [1]. Signals of interest include sound files, images, radar,
More informationx[n] = x a (nt ) x a (t)e jωt dt while the discrete time signal x[n] has the discrete-time Fourier transform x[n]e jωn
Sampling Let x a (t) be a continuous time signal. The signal is sampled by taking the signal value at intervals of time T to get The signal x(t) has a Fourier transform x[n] = x a (nt ) X a (Ω) = x a (t)e
More informationSubband Coding and Wavelets. National Chiao Tung University Chun-Jen Tsai 12/04/2014
Subband Coding and Wavelets National Chiao Tung Universit Chun-Jen Tsai /4/4 Concept of Subband Coding In transform coding, we use N (or N N) samples as the data transform unit Transform coefficients are
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 informationFourier Series and Fourier Transforms
Fourier Series and Fourier Transforms EECS2 (6.082), MIT Fall 2006 Lectures 2 and 3 Fourier Series From your differential equations course, 18.03, you know Fourier s expression representing a T -periodic
More informationINTRODUCTION TO NMR and NMR QIP
Books (NMR): Spin dynamics: basics of nuclear magnetic resonance, M. H. Levitt, Wiley, 2001. The principles of nuclear magnetism, A. Abragam, Oxford, 1961. Principles of magnetic resonance, C. P. Slichter,
More informationAnalog Digital Sampling & Discrete Time Discrete Values & Noise Digital-to-Analog Conversion Analog-to-Digital Conversion
Analog Digital Sampling & Discrete Time Discrete Values & Noise Digital-to-Analog Conversion Analog-to-Digital Conversion 6.082 Fall 2006 Analog Digital, Slide Plan: Mixed Signal Architecture volts bits
More informationOn the Frequency-Domain Properties of Savitzky-Golay Filters
On the Frequency-Domain Properties of Savitzky-Golay Filters Ronald W Schafer HP Laboratories HPL-2-9 Keyword(s): Savitzky-Golay filter, least-squares polynomial approximation, smoothing Abstract: This
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 informationContinuous wavelet transform Python module
Continuous wavelet transform Python module Erwin Verwichte University of Warwick erwin.verwichte@warwick.ac.uk 31 Oct 017 1 Introduction A continuous wavelet is a well-known fundamental tool that allows
More informationMGA Tutorial, September 08, 2004 Construction of Wavelets. Jan-Olov Strömberg
MGA Tutorial, September 08, 2004 Construction of Wavelets Jan-Olov Strömberg Department of Mathematics Royal Institute of Technology (KTH) Stockholm, Sweden Department of Numerical Analysis and Computer
More informationWavelets and Multiresolution Processing
Wavelets and Multiresolution Processing Wavelets Fourier transform has it basis functions in sinusoids Wavelets based on small waves of varying frequency and limited duration In addition to frequency,
More informationThe BTE with a High B-field
ECE 656: Electronic Transport in Semiconductors Fall 2017 The BTE with a High B-field Mark Lundstrom Electrical and Computer Engineering Purdue University West Lafayette, IN USA 10/11/17 Outline 1) Introduction
More informationECE 636: Systems identification
ECE 636: Systems identification Lectures 3 4 Random variables/signals (continued) Random/stochastic vectors Random signals and linear systems Random signals in the frequency domain υ ε x S z + y Experimental
More informationFile name: Supplementary Information Description: Supplementary Figures, Supplementary Notes and Supplementary References
File name: Supplementary Information Description: Supplementary Figures, Supplementary Notes and Supplementary References File name: Peer Review File Description: Optical frequency (THz) 05. 0 05. 5 05.7
More informationELEN E4810: Digital Signal Processing Topic 11: Continuous Signals. 1. Sampling and Reconstruction 2. Quantization
ELEN E4810: Digital Signal Processing Topic 11: Continuous Signals 1. Sampling and Reconstruction 2. Quantization 1 1. Sampling & Reconstruction DSP must interact with an analog world: A to D D to A x(t)
More information