The Hilbert spectrum via wavelet projections
|
|
- Muriel Watts
- 5 years ago
- Views:
Transcription
1 The Hilbert spectrum via wavelet projections By S. Olhede and A.T. Walden Department of Mathematics, Imperial College of Science, Technology and Medicine, 8 Queen s Gate, London SW7 BZ, UK. Non-stationary signals are increasingly analysed in the - domain to determine the variation of components with. It was recently proposed in this journal that such signals could be analysed by projections onto the - plane giving a set of monocomponent signals. These could then be converted to analytic signals using the Hilbert transform and their instantaneous calculated, which when weighted by the energy yields the Hilbert energy spectrum for that projection. Agglomeration over projections yields the complete Hilbert spectrum. We show that superior results can be obtained using waveletbased projections. The maximal-overlap (undecimated/stationary/translation invariant) discrete wavelet transform and wavelet packet transforms are used, with the Fejér-Korovkin class of wavelet filters. These transforms produce decompositions which are conducive to statistical analysis, in particular enabling noise reduction methodology to be developed and easily and successfully applied. Keywords: denoising; empirical mode decomposition; Hilbert spectrum; instantaneous ; wavelet packets; wavelet transform. Introduction Most naturally occurring or man-made signals display non-stationary behaviour, i.e., some structural aspect such as their spectrum, statistical correlation or transfer function vary with. Classical analysis techniques which assume stationarity are unsuitable, and investigators are drawn to using - methods (e.g., Cohen, 995) to analyse such signals. Consider a generally complex-valued signal z(t) which can be written in terms of a -domain amplitude A(t) and phase φ(t) asz(t) =A(t)e iφ(t). Suppose its energy is highly concentrated along some trajectory in the - plane. A general representation would then be in terms of a a curve defined by the average at each the so-called instantaneous. Boashash (99a) remarks that for monocomponent signals which obey certain conditions the instantaneous will have a clear physical meaning. However, multicomponent signals are frequently encountered. Loughlin and Tacer (997) consider the two-component signal z(t) =A e iπft + A e iπft, where A and A are real constants. The instantaneous of the signal is ν(t) = φ (t) π = (f + f ) (f f )(A A + ) {A + A +A A cos[π(f f )t]}. (.) Proc. R. Soc. Lond. A., 46, , 4 TEX Paper
2 S. Olhede and A.T. Walden This is not bounded by f or f unless A = A, in which case the instantaneous is consistent with the interpretation as the average at each of this multicomponent signal. They also consider more general twocomponent signals and find in general that the instantaneous of the signal exhibits large fluctuations and extends beyond the band defined by the instantaneous of the components. Instantaneous is generally a poor measure for multicomponent signals. For this reason the analysis of such multicomponent signals is best done using projections onto the - plane corresponding to monocomponent signals which individually exhibit coherent behaviour at any one. Recently in this journal a - energy spectrum estimator was proposed by Huang et al (998) which projects the signal in the - plane (using the empirical mode decomposition (EMD) ) and then calculates the instantaneous for each projection. Each projection is intended to ensure coherent behaviour at any, or a monocomponent signal, so that the instantaneous is well-behaved. Finally for each projection the appropriate energy is associated with the instantaneous at each, and amalgamation over projections gives the overall Hilbert energy spectrum. The EMD projections are data-dependent and ad-hoc constructions, which vary with each signal this is potentially a strength (adaptivity) and a weakness (resistance to statistical analysis). The approach of Huang et al (998) is aimed at analysis of both non-stationary and non-linear proceses. Our interest is directed at the analysis of non-stationary processes for which we find the EMD method produces oscillatory or poorly-defined Hilbert spectra, often with evidence of notable mode mixing. We show that superior results can be obtained by replacing the EMD projections by wavelet-based ones. In particular we utilise the maximal-overlap (undecimated/stationary/translation invariant) discrete wavelet transform and wavelet packet transforms, denoted by MODWT and MODWPT, respectively. These non-adaptive transforms produce decompositions which are conducive to statistical analysis, enabling noise reduction methodology to be developed and utilised. In and 3 we give a brief summary of the EMD, and the Hilbert spectrum, respectively. The background to wavelet-based - projections via socalled MODWT detail coefficients is described in 4 along with the computation of the Hilbert spectrum. Projections via the wavelet packet transform (MODWPT), a generalization of the MODWT, are developed in 5, where details are given on choosing a suitable set of wavelet packets. The importance of low - leakage for the projections is also addressed in 5, where it is shown that the Fejér- Korovkin class of wavelet filters are particularly appropriate. Examples showing that the wavelet-based approach gives excellent results, and outperforms EMD, are given in 6. Sums of intermittent and complete sinusoids are analysed, as well as a bat echolocation pulse. The theory behind Hilbert spectrum estimation in the presence of noise, via thresholding of wavelet details, is given in 7, followed by an example showing its efficacy. Throughout the paper a t in parentheses indicates continuous, and a subscript t indicates discrete. Proc. R. Soc. Lond. A., 46, , 4
3 The Hilbert spectrum via wavelet projections 3. The Empirical Mode Decomposition (a) Projections in -: continuous Given a real-valued continuous- signal x(t), consider a series of projection operators P Rj which project x(t) into a - region R j. These operators need not be orthogonal but we assume that J + j= (P R j x)(t) =x(t) for some J. For j =,...,J define d j (t) =(P Rj x)(t) and s J (t) =(P RJ x)(t). Then + J x(t) = d j (t)+s J (t). (.) j= Such an additive decomposition of x(t) is achieved by the empirical mode decomposition method. (b) The EMD algorithm Huang et al (998, 999) define an empirical mode decomposition (EMD) of a signal into intrinsic mode functions (IMFs). An IMF has the same number of zero crossings and extrema, and symmetric envelopes defined by the local maxima and minima. Each such IMF should give rise to a well-defined instantaneous. The EMD of x(t) proceeds as follows. Upper, u (t), and lower, l (t), envelopes of x(t) are found by interpolation between extrema. The local mean function is defined as m (t) = [u (t)+l (t)]. Ideally h, (t) =x(t) m (t), should be an IMF, but invariably a sifting refinement is required: h, (t) replaces x(t) and the above step is repeated to give m (t) and h, (t). After k such repeats h,k (t) is obtained and a stopping criterion equates it to the first IMF, d (t) say. Various stopping criteria exist. Huang et al (998, p. 9) test whether the accumulated relative squared differences between h,k (t) and h,k (t) is less than.3; see also Magrin-Chagnolleau and Baraniuk (999). Huang et al (999, p. 45) look for an equal number of extrema and zero crossings. A more detailed combined criterion is used by Rilling et al (3) and this is the scheme used in this paper; for details see also flandrin/software.html. We have also repeated our analyses using the other stopping criteria, and the basic conclusions of the paper are unchanged. Denoting the sum of the mean functions by m tot (t) we note that d (t) =x(t) m tot (t). The second IMF is found by performing the above procedure on the residual s (t) =x(t) d (t), and so on, until s J (t) satisfies a criterion indicating that no more IMFs can be extracted (no more than one extrema). s J (t) will be a trend or constant. Then (.) is obtained. Locally, each IMF contains lower oscillations than that previously extracted. The EMD gives an adaptive subband decomposition of the signal. If the - plane has been tiled ideally, each projection will contain at most one monocomponent signal. 3. The Hilbert spectrum and instantaneous Given a real-valued x(t), and hence real-valued d j (t), we can create the analytic signal d a j (t) = d j(t) +ih{d j (t)} = d j (t) +id (H) j (t), where H{ } denotes Proc. R. Soc. Lond. A., 46, , 4
4 4 S. Olhede and A.T. Walden the Hilbert transform. Then the phase is φ j (t) = tan (I{d j (t)} / R{d j (t)}) = tan (d (H) j (t)/d j (t)), and the amplitude is A j (t) = d a j (t) = {d j (t)+[d(h) j (t)] }. The corresponding instantaneous is ν j (t) =φ j (t)/(π). The Hilbert energy spectrum of d j (t) is given by S j (t, f) =A j(t)δ(f ν j (t)), where δ( ) is the Dirac delta function. The overall Hilbert energy spectrum for x(t) is defined (Huang et al 998) as the amalgamation of the spectra of each of the d j (t) s. We can also define a -varying mean for the entire signal, a useful summary curve. ν (t) = J j= A j (t) ν j (t) J j= A j (t), 4. Discrete wavelet decompositions (a) Wavelet projections Analogous to section (a) we consider projections of a real-valued signal which attempt to produce monocomponent separation, only here we use discrete- methods. Assume we have sampled a continuous- signal at intervals t = to get a vector of observations X = [X,...,X N ] T. Assume also we have a series of projection matrices P Rj, which project a vector X into a particular region R j. These operators need not be orthogonal but we shall assume J+ j= P R j X = X for some particular J.Forj =,...,J define D j = P Rj X and S J = P RJ X. Then + X = J j= D j + S J. (4.) Decomposition (4.) is the discrete equivalent of (.). So we can decompose the sampled signal into a linear combination of the contribution of X in projected spaces. If the - plane has been tiled ideally, each projection will contain at most one monocomponent signal for which the concept of instantaneous is well defined. The maximal overlap discrete wavelet transform (MODWT) is a version of the discrete wavelet transform which does away with the downsampling step, thus ensuring the important property of circular shift equivariance. Transforms that are essentially the same include the undecimated discrete wavelet transform (Shensa, 99) and translation invariant discrete wavelet transform (Liang and Parks, 996). Full details of the MODWT may be found in Percival and Walden (, chapter 5); we adopt their notation here. For t =,...,N the jth level wavelet and scaling coefficients are given by W j,t = L j l= h j,l X (t l) modn and Ṽj,t = Lj l= g j,l X (t l) modn, where { h j,l } and { g j,l } are level-j filters of length L j = ( j )(L ) + derived from the basic DWT wavelet and scaling filters {h l } Proc. R. Soc. Lond. A., 46, , 4
5 The Hilbert spectrum via wavelet projections 5 and {g l } of length L. (Note the tilde notation is used for MODWT filters and coefficients to distinguish them from the DWT equivalents). In matrix form we can write W j = W j X and Ṽj = ṼjX, where W j and Ṽj are N-length column vectors, and W j and Ṽj are N N matrices composed of elements from { h j,l } and { g j,l } respectively, and zeros. Furthermore, the signal can be recovered as in (4.), where D j = W j T W j X and S j = ṼT j ṼjX (Percival and Walden, p. 73). The corresponding projection matrices are P Rj = W j T W j and P = RJ + ṼT J Ṽ J, where P Rj is approximately the identity matrix on - region R j, and zero elsewhere, where { [, (N )] { [ ) ( R j =, j ]} j+, j+ if j j J [, (N )] { [ ], J + J } if j = J +. + This is then an example of a set of projection operators as described above. S j represents the over-all -varying mean trend, while the D j represent local - information. The MODWT of course tiles the - plane in octave bands. The elements of S j and D j are called jth level smooth and detail coefficients. We can write (Percival and Walden, p. 99) D j,t = L j l= (L j ) h j h j,l X (t l) modn, where h j h j,l denotes the cross-correlation of { h j,l } with itself, i.e., the autocorrelation of { h j,l }. The effective filter is symmetric, zero-phase, and considerably smoother than { h j,l } itself. Analogous results hold for S j. So while the utility of wavelet methods in efficiently representing discontinuities is well-known, in fact our projections are well-suited to represent smooth, possibly sinusoidal, sequences. (b) The Hilbert spectrum via wavelets The MODWT detail coefficients for level j, { D j,t,t=,...,n }, will be real for any real-valued wavelet filters. We can create analytic wavelet details via D j,t +i( D j q) t = D j,t +i D (H) j,t, (4.) where {q l } is the Hilbert transform filter and D (H) j,t transformed version of D j,t. The amplitude and phase sequences are given by denotes the discrete Hilbert A j,t = { D j,t +[ D (H) j,t ] }; φ j,t = tan (H) ( D j,t / D j,t ). (4.3) We compute the derivative of the phase at discrete t via a fourth-order generalized phase difference estimator (Boashash, 99b, p. 54), φ j,t = (φ j,t 8φ j,t +8φ j,t+ φ j,t+ ). (4.4) Proc. R. Soc. Lond. A., 46, , 4
6 6 S. Olhede and A.T. Walden G(4 k N ) W, k H(4 N ) G( k N ) W, G( k N ) H( k N ) H(4 k N ) W, k G(4 N ) X = W, G(4 k N ) W, k H(4 N ) k H( N ) W, H( k N ) G( k N ) H(4 k N ) W,3 G(4 k N ) W 3, W 3, W 3, W 3,3 W 3,4 W 3,5 W 3,6 W 3,7 /6 /8 3/6 /4 5/6 3/8 7/6 / f Figure. First three levels of the MODWPT of X. All coefficient vectors together form a (maximal overlap) wavelet packet table. This formula is skew-symmetric about t. At the end-points we cannot use this centred five-point formula but use instead φ j,t = ± ( 5φ j,t +48φ j,t± 36φ j,t± +6φ j,t±3 3φ j,t±4 ). (4.5) Both formulae can be found in Burden and Faires (993, p. 6). The instantaneous is then computed from ν j,t = φ j,t /(π). Now, for a large M, define M frequencies f l = l f, l =,...,M, where f = /[(M )], covering the interval [, /]. Then the Hilbert energy spectrum of { D j,t } may be defined as S j;t,fl = A j,tδ l, νj,t/ f, t =,...,N ; f l =,...,/, (4.6) where δ l,m is the Kronecker delta, and x denotes the integer closest to x; we used M = 5. The overall Hilbert energy spectrum is the amalgamation of the S j;t,fl s over the j. The -varying mean for the entire signal follows as ν t = J j= A j,t ν j,t J. (4.7) j= A j,t 5. Discrete wavelet packet decompositions (a) Wavelet packet projections The MODWT has the disadvantage of imposing a fixed octave band tiling on the - plane. The maximal overlap discrete wavelet packet transform (MODWPT) (Walden & Contreras Cristán, 998) is a generalization of the MODWT which at level j of the transform partitions the axis into j equal width bands, often labelled n =,..., j. This gives better resolution. Let G(f) = L l= g le iπfl be the transfer function of { g l } where Proc. R. Soc. Lond. A., 46, , 4
7 The Hilbert spectrum via wavelet projections 7 level /6 /8 3/6 /4 5/6 3/8 7/6 / Figure. The shaded segments determine a possible disjoint dyadic decomposition a partitioning of the interval [, /] which may be derived from the WPT of Fig.. g l = g l /, and let the transfer function corresponding to { h l } with h l = h l /, be defined similarly. At the first level of the transform X is circularly filtered by the low-pass filter { g l }, with corresponding transfer function G(f), to give an N-length vector of first level coefficients W, and X is also circularly filtered by the high-pass filter { h l }, with transfer function H(f), to give the vector W,. For subsequent levels j of the transform we insert j zeros, j, between the elements of { g l }; the resulting filter has a transfer function given by G( j f). We can do likewise for { h l } to obtain H( j f). The jth level coefficients are obtained by filtering the level j coefficients with the circular filter having discrete Fourier transform { H( j k N )} or { G( j k N )}, as appropriate. The filtering steps up to level j =3 are shown in Fig. and the collection of MODWPT coefficients is called a (maximal overlap) wavelet packet table (WPT). If, as here, the ordering of the filters means that at level j, band n is associated with frequencies I j,n = [ ], j+, j+ for n =, and I j,n = [ ) ( n+ n, j+ n ] j+, n+ j+, for n>, then the transform is j+ said to be sequency ordered. At level j and band n the wavelet packet coefficients for t =,...,N can be written as W j,n,t = L j l= ũ j,n,l X (t l) modn where {ũ j,n,l } is a jth level and nth band MODWPT wavelet packet filter of length L j =( j )(L )+ with a transfer function which may be derived from the filtering steps shown in Fig. ; for example for j =3,n= 3 this transfer function is given by G(f) H(f) G(4f). In matrix form we can write W j,n = W j,n X where W j,n is a N-length column vector, and W j,n is a N N matrix composed of elements from {ũ j,n,l } and zeros. If we let D j,n = W j,n T W j,n X then the signal can be recovered via X = j D n= j,n. For n> the set { D j,n,t,t=,...,n } are level j, band n, MODWPT detail coefficients. However, the decomposition in (4.) makes use of different levels of the transform, and the same can be done for the MODWPT. The doublets (j, n) that form the indices for the MODWPT coefficient vectors W j,n can be collected together to form a set N {(j, n) :j =,...,J ; n =,..., j }. Each (j, n) N Proc. R. Soc. Lond. A., 46, , 4
8 8 S. Olhede and A.T. Walden is associated with the interval I j,n. The WPT can be used to form a large collection of different transforms known as disjoint dyadic decompositions. By definition, each such decomposition is associated with a subset, say C, of doublets from N that satisfies two properties. First, the union of all the intervals I j,n for each (j, n) Cis exactly the interval [, /], and, second, if (j, n) and (j,n ) are any two distinct elements of C, then I j,n and I j,n have no frequencies in common. One such decomposition is shown by the grey shading in Fig. ; here C = {(, ), (, 3), (3, ), (3, ), (3, 4), (3, 5)}). Fundamentally, for any such decomposition C we have X = D j,n + D j,, (j,n> C) for some j J. The last term covers low frequencies and trends. The projection matrices are P Rj,n = W j,n T W j,n where P Rj,n is approximately the identity matrix on - region R j,n, and zero elsewhere, where R j,n =[, (N )] I j,n. This is another example of a set of projection operators. When utilising the wavelet packet transform we must decide on some method of choosing a suitable disjoint dyadic decomposition C which has the effect of producing a - tiling commensurate with separating the original signal into a set of monocomponent signals. (b) Choosing the disjoint dyadic decomposition Consider an algorithm starting at level j = J. If the lowest feature of interest in the signal occurs at a higher than f L say, we need to choose J so that <f J + L. (This is because the Hilbert spectrum discards the lowest components to avoid the zero or infinite scale). Consequently, D J, and D J, are always preserved in any decomposition. Now consider a measure of the bandwidth of the details D j,n,t,t=,...,n at level j, band n>. The discrete Fourier transform (DFT) of the details at frequencies f k = k/n is ξ j,n,k = N t= D j,n,t e iπtk/n, k =,...,N, and the normalised modulus squared, which sums to unity, takes the form / N/ Υ j,n,k = ξ j,n,k ξ j,n,k, k =,..., N/, k= where x denotes the greatest integer x. We denote the bandwidth by B j,n with B j,n = N/ (f k f j,n ) Υ j,n,k where f j,n = k= N/ k= f k Υ j,n,k. The fundamental idea behind the decomposition algorithm is illustrated in Fig. 3. Consider whether a band I j,n should be split into two subbands I j+,n and I j+,n+ ; we note Proc. R. Soc. Lond. A., 46, , 4
9 The Hilbert spectrum via wavelet projections Figure 3. Illustration of points (i) (left) and (ii) (right) of the decomposition algorithm. Left: one chirp signal extending from f =.5 to f =.5 (diagonal line) The bandwidth (long vertical line) is.65. If split into two, we obtain two chirps extending from f =.5 to.875 and one from.875 to.5 The separate chirps both have a bandwidth (short vertical lines) of.35 which is half that of the entire signal. Right: two separate chirps, each with bandwidths of about.66. If they are regarded as one signal, then this has a bandwidth of about.468 > (i) if I j,n contains a monocomponent signal then if it is split, the sum of the ensuing bandwidths is equal to that of the entire signal (Fig. 3, left), while (ii) if the subbands each contain separate monocomponent signals then the sum of their bandwidths will be less than the bandwidth obtained if they are regarded as one signal (Fig. 3, right). The algorithm starts at level j = J and moves up the tree. If B j,n >α(b j+,n + B j+,n+ ) (5.) we keep D j+,n and D j+,n+ (i.e., use subbands I j+,n and I j+,n+ ) and in order to ensure that the split remains as we move up the tree we set B j,n =. Here α is a constant chosen to avoid splitting the band due simply to random variation. It is more important to avoid not splitting a band when we should, rather than splitting when not required, as we seek monocomponent signals in each band if we do not split a band when we should, then a meaningless instantaneous will result. A value α =. is suitable for this requirement. A slight problem arises when dealing with - tiles which are empty of signal, so that the bandwidth is ill-defined. It is advisable to factor in the actual energy of the details at level j and band n and modify the bracketed term in (5.) to read N B j+,n t= A j+,n,t + B j+,n+ N t= A j+,n+,t N t= A j+,n,t + N t= A j+,n+,t (5.) where A j,n,t,t=,...,n are the amplitudes of the analytic MODWPT details, defined next. Proc. R. Soc. Lond. A., 46, , 4
10 S. Olhede and A.T. Walden (c) The Hilbert spectrum via wavelet packets We can create analytic MODWPT details for level j and band n> using D a j,n,t = D j,n,t +i( D j,n q) t = D j,n,t +i D (H) j,n,t. (5.3) Given the analytic details (5.3), the corresponding amplitude, A j,n,t, phase, φ j,n,t and instantaneous, ν j,n,t, are computed as in section 4(b), (an additional subscript n being required throughout). The Hilbert energy spectrum of { D j,n,t } is given by S j,n;t,fl = A j,n,tδ l, νj,n,t/ f, t =,...,N ; f l =,...,/, (5.4) and the -varying mean for the entire signal follows as (j,n> C) ν t = A j,n,t ν j,n,t. (5.5) (j,n> C) A j,n,t (d) Leakage of - projections A major problem with any projection is the unavoidable leakage caused by the fact that no projection can create a signal which is simultaneously and limited. Using finite-length filters, wavelet packets lose some of their resolution at high frequencies (e.g., Hess-Nielsen & Wickerhauser, 996). Nielsen () derived classes of half-band low-pass DWT filters {g l,l=,...,l } such that the L p norm of the difference between the power transfer function G(f), (where G(f) is the Fourier transform of {g l }), and the power transfer function of the ideal half-band filter rectangular on [ /4, /4] is O(/L) for p [, ]. One class of such optimal asymptotic resolution filters are the Fejér-Korovkin filters. Fig. 4 gives the power transfer functions of level j = 3 wavelet packet filters {ũ j,n,l } generated from the L = 8 point Fejér-Korovkin filter {g l } abbreviated to FK(8) and L = point Daubechies least asymmetric filter {g l } abbreviated to LA(). Even though it is shorter by two points than the frequently used Daubechies least asymmetric filter, the Fejér-Korovkin filter gives rise to the better concentrated wavelet packets, the power transfer functions being in particular more rectangular over their pass-bands I j,n. The minimally short two point Haar wavelet filter generates wavelet packet filters with much poorer concentration than the LA(). The efficacy of the FK(8), LA() and Haar wavelet filters for Hilbert spectrum or -varying mean calculation is demonstrated in the examples which follow. 6. Examples In this section we will compare plots of the Hilbert energy spectrum and the varying mean made using the MODWT and EMD approaches. For the EMD approach the same computational scheme is used as set out in Section 4b except that the wavelet detail coefficients are replaced by the IMF sequences; in particular the same phase derivative method is used see (4.4) and (4.5). Proc. R. Soc. Lond. A., 46, , 4
11 The Hilbert spectrum via wavelet projections power power /6 /8 3/6 /4 5/6 3/8 7/6 / Figure 4. Power transfer functions of level j = 3 wavelet filters {ũ j,n,l } for (from top to bottom) bands n =,...,7, generated from the 8-point Fejér-Korovkin filter (solid lines) or -point Daubechies least asymmetric filter (dashed lines). (a) Sum of sinusoids Consider a signal which is the sum of two equal-amplitude sine waves: x (t) = sin(πf t) + sin(πf t), (6.) where f =.8 and f =.6; the theoretical mean is.. Figs. 5(a) (c) show the Hilbert energy spectrum and the -varying mean calculated using the MODWT with J = 5 when this signal is sampled with a unit sample interval, N = 5 s starting at t =. The half-band low-pass filter {g l,l=,...,l } used in the MODWT was (a) Haar, (b) Daubechies LA(), and (c) Fejér-Korovkin FK(8). Leakage-generated oscillations in the Hilbert spectrum are present in all three cases. Leakage of the projections means that in a particular band not only is the appropriate sinusoid present, but also sinusoidal contributions from neighbouring bands, with all the resulting complications such as (.). These oscillations decrease to almost nothing as the concentration of the filters improves; the FK(8) result is strikingly good. Fig. 6 gives the MODWT detail sequences for j =,...,5 for the case of the Fejér-Korovkin FK(8) filter. The five sequences are ordered (top to bottom) according to the top five bands delineated by dotted lines in Figs. 5(c); it is clear that leakage between bands is very small. The representation of the sinusoids is excellent as predicted in the discussion at the end of 4a. Fig. 7 gives the five IMFs produced by EMD. The differences between the components determined by the MODWT in Fig. 6 and EMD in Fig. 7 do not seem significant, yet they are; this can be seen from the Hilbert spectrum given in Fig. 5(d) Proc. R. Soc. Lond. A., 46, , 4
12 S. Olhede and A.T. Walden (a) (b) (c) 5 5 (d) 5 Figure 5. Sum of sines example. f =.8 and f =.6 and theoretical mean is. (a) Left: the Hilbert spectrum using the MODWT, Haar wavelet filter, and J = 5. Right: corresponding -varying mean. (b) Left: the Hilbert spectrum using the MODWT, Daubechies LA() wavelet filter, and J = 5. Right: corresponding -varying mean. (c) Left: the Hilbert spectrum using the MODWT, Fejér-Korovkin FK(8) wavelet filter, and J = 5. Right: corresponding -varying mean. The horizontal dotted lines delineate the MODWT octave bands. (d) Left: the Hilbert spectrum using EMD and five IMFs. Right: corresponding -varying mean. using the five IMFs where relatively poor resolution is achieved and the result is inferior to all but the Haar MODWT. This is somewhat surprising since we would expect the EMD method to produce good representations of sinusoids since theoretically sinusoidal components are IMFs. Proc. R. Soc. Lond. A., 46, , 4
13 The Hilbert spectrum via wavelet projections Figure 6. The five MODWT detail sequences (FK(8) filter) for sum of sines Figure 7. Five IMFs for sum of sines. (b) Intermittent and complete sinusoids Here we consider a signal made up from four sinusoids: x (t) =.sin(πf t)i T (t) + sin(πf t)i T (t)+.sin(πf 3 t) + sin(πf 4 t), (6.) Proc. R. Soc. Lond. A., 46, , 4
14 4 S. Olhede and A.T. Walden (a) (b) 5 Figure 8. Intermittent and complete sines example. f =.938, f =.875, f 3 =.35 and f 4 = (a) Left: the Hilbert spectrum using the MODWPT, Fejér-Korovkin FK(8) wavelet filter, and J = 4. The horizontal dotted lines delineate the MODWPT bands. Additional short tick marks indicating the frequencies f,...,f 4 are not visible under the estimated frequencies. Right: corresponding -varying mean on a magnified scale, with the theoretical mean (dotted) overlain. (b) Left: the Hilbert spectrum using EMD and nine IMFs. The additional short tick marks indicating the frequencies f,...,f 4 are now visible. Right: corresponding -varying mean. where f =.938, f =.875, f 3 =.35, f 4 =.4375, I T (t) =ift T and zero otherwise, T =[,.97] [.4, ) and T =[,.597] [.748, ). The first two sinusoids are intermittent. This signal is sampled with a unit sample interval, N = 4 s starting at t =. The two higher sinusoids have frequencies both falling in the band (/4, /] and so cannot be separated by the MODWT. The more general MODWPT is thus required, and we can allow the algorithm in section 5(b) to select an appropriate disjoint dyadic decomposition Proc. R. Soc. Lond. A., 46, , 4
15 The Hilbert spectrum via wavelet projections 5 the results are shown in Fig. 8(a). The MODWPT method gives a very good and accurate description of the signals, with a few oscillations near the discontinuities in the signal. The corresponding results using EMD are given in Fig. 8(b). This used the nine IMFs found by the algorithm. (Use of only the the first four IMFs there are four monocomponent signals present produced an indistinguishable result.) Huge oscillations are found with the EMD algorithm, and there are also problems with mode-mixing. (c) Application to real data We now analyse the digitized echolocation pulse shown in Fig. 9(a) emitted by the Large Brown Bat, Eptesicus Fuscus. The sample interval is 7µs. Fig. 9(b) shows the MODWPT results, while the corresponding results for EMD are shown in Fig. 9(c). For EMD the seven IMFs found by the algorithm are used, but taking only the first four caused no appreciable difference. The MODWPT clearly outperforms the EMD approach in resolving the hyperbolic chirps present in this signal. Notice that the algorithm of section 5b for choosing the disjoint dyadic decomposition for the MODWPT has worked successfully: at any the different chirps have been separated into different subbands. 7. Denoising (a) Theory Since we are in general dealing with signals contaminated by noise, it makes sense to denoise the signal where possible. We assume x(t) = b(t) + ɛ(t), where b(t) is deterministic signal and ɛ(t) is real-valued noise such that discretization at unit sample interval gives the N-length vector X = b + ɛ, (7.) where {ɛ t,t =,...,N } is a white Gaussian noise sequence with mean zero and variance σɛ. Coifman and Donoho (995) introduced the technique of cycle spinning for denoising of X in a translation invariant way; this is equivalent (Percival and Walden,, p. 49) to applying standard thresholding to the MODWT wavelet coefficients. Here our denoising scheme will determine which analytic wavelet detail magnitudes A j,t should be used in the Hilbert energy spectrum of (4.6) and -varying mean given by (4.7). The MODWT details vector is given by D j = W j T W j X so that the additive noise vector is transformed to the noise details vector Ẽj = W j T W j ɛ. But the covariance of the (Gaussian) elements of Ẽj is given by (Appendix A) cov{ẽj,t, Ẽj,t } { σ ɛ / j, t = t, σɛ {sin[π(t t )/ j ] sin[π(t t )/ j ]} [π(t t )], t t. (7.) Now suppose {η t } is a stationary process with cov{η t,η t } = s t t, say, and define the analytic process by η (a) t = η t +i(η q) t = η t +iη (H) t, Proc. R. Soc. Lond. A., 46, , 4
16 6 S. Olhede and A.T. Walden.. amplitude (s) (a) 7 x 4 7 x 4 7 x 4 7 x (s) (b).. (s).. (s) (c).. (s) Figure 9. (a) Bat signal. (b) Left: the Hilbert spectrum using the MODWPT, Fejér-Korovkin FK(8) wavelet, and J = 3. Right: corresponding -varying mean. (c) Left: the Hilbert spectrum using EMD and seven IMFs. Right: corresponding -varying mean. where, {q l } is the Hilbert transform filter. Then, since the Hilbert transform is an all-pass filter, it does not affect the covariance structure of {η t } so that cov{η (H) t,η (H) t } = s t t, also. Moreover, E{η t η (H) t } = E{η t q l η t l} = τ l s t t +l q l, so that if t = t we get cov{η t,η (H) t } = l s lq l =, since {s l } is even and {q l } is odd. Hence at any point t, the random variables η t and η (H) t are uncorrelated; however this does not mean the processes {η t } and {η (H) t } are uncorrelated, let alone independent. However, if {η t } is Gaussian, then the processes {η t } and {η (H) t } are Proc. R. Soc. Lond. A., 46, , 4
17 The Hilbert spectrum via wavelet projections amplitude (a) (b) (c) 5 5 (d) 5 Figure. (a) Noisy data. (b) Left: the Hilbert spectrum using the MODWT, Fejér-Korovkin FK(8) wavelet, and J = 4, with thresholding of wavelet details. Right: corresponding -varying mean. (c) As in (b) but with no thresholding. (d) Left: the Hilbert spectrum using EMD and six IMFs. Right: corresponding -varying mean. jointly (dependent) Gaussian (Papoulis, 97, p. 475) and so the random variables η t and η (H) t are independent and identically distributed for any specified t. Since from (7.) cov{ẽj,t, Ẽj,t } depends only on t t we let Ẽj,t play the role of the stationary process η t and conclude that, up to inaccuracies introduced by filtering end-effects, the random variables Ẽj,t and Ẽ (H) j,t can be treated as independent and identically distributed for each t =,...,N. In particular Ẽj,t and Ẽ (H) j,t have the same variance given by σ = Ẽ,j σ ɛ / j ; see (7.). To estimate the noise standard deviation σ ɛ we use the MODWT-based MAD estimator (Percival and Proc. R. Soc. Lond. A., 46, , 4
18 8 S. Olhede and A.T. Walden Walden,, p. 49), ˆσ ɛ = / median{ W,, W,,..., W,N }, (7.3).6745 which gives an estimate of ˆσ Ẽ,j =ˆσ ɛ / j. The independence and Gaussianity of Ẽj,t and Ẽ (H) j,t ) ( Ẽ (a) = ( Ẽj,t σẽ,j ) + (Ẽ (H) j,t σẽ,j j,t σẽ,j means that ) χ where χ means is distributed as a chi-square random variable with two degrees of freedom, and Ẽ (a) j,t = {Ẽ j,t (H) +[Ẽ j,t ] }. We thus have a sequence of dependent chi-square random variables ( Ẽ (a) j,t /σ Ẽ,j ),t=,...,n. Thresholding proceeds under the hypothesis of no deterministic signal being present; in this case D j = W j T W j X = Ẽj = W j T W j ɛ and hence A j,t Ẽ (a) j,t. A suitable threshold level for N independent chi-square random variables was developed in Sardy (), and for our case of N dependent chi-square random variables we can show (Olhede and Walden, 3) that the threshold level of Sardy () provides a conservative threshold level. Under the hypothesis of no deterministic signal present we use the hard-thresholding rule, combined with Sardy s threshold level of [ log(n log N)], { A j,t, A (ht) j,t = if A j,t > ˆσẼ,j [ log(n log N)], otherwise. for levels j =,...,J. The Hilbert energy spectrum and -varying mean calculations then proceed as in (4.6) and (4.7) with A (ht) j,t replacing A j,t Notice that our strategy is to threshold the amplitudes only, leaving the phases unchanged, since only when the amplitude is low do we wish to eliminate the contribution to X. Thresholding the real and imaginary parts of Dj,t separately would change the phase of the analytic signal and hence is avoided. (b) Example Here we consider a linear chirp, b(t) = sin( πt 4N ), plus noise. This signal has instantaneous given by t/(4n), i.e., a straight line with changing from to.5 as goes from to N. We use a discretization of the signal at 5 unit sample intervals, plus a white Gaussian noise {ɛ t,t =,...,5} with σ ɛ =.5. A realization is shown in Fig. (a). The MODWT-based Hilbert spectrum and corresponding -varying mean following the thresholding scheme are given in Fig. (b), while the MODWT results omitting thresholding are shown in Fig. (c). The EMD results are shown in Fig. (d) where the number of IMFs used was six, as found by the algorithm; a reduction to the first three IMFs produced worse results. Our denoising method, with thresholding of the amplitudes only, works very well. With regard to the Hilbert spectrum, the standard MODWT-based approach again performs better than the EMD, but the omission of thresholding significantly worsens the results. Proc. R. Soc. Lond. A., 46, , 4
19 The Hilbert spectrum via wavelet projections 9 8. Conclusions The analysis of a possibly multicomponent non-stationary signal via projections onto the - plane followed by computation of the Hilbert spectrum is an appealing way to determine the -variation of components. However, the EMD method seems a problematic step, and we have shown that superior results can be obtained with wavelet-based (MODWT or MODWPT) projections utilising the Fejér-Korovkin class of wavelet filters. These transforms produce decompositions which are conducive to statistical analysis, with the important result that a denoising step can easily be incorporated into this approach. In the bat data analysis of Fig. 9, we note that although, as required in Hilbert spectrum multicomponent analysis, at any the different chirps have been successfully separated into different subbands, each component spans at least two subbands and is noticeably affected by the band-edge imperfections. A paper covering work on increasing the flexibility of our approach to deal with such drawbacks and cope better with highly non-linear chirps, via suitable pre- and post-processing, is nearing completion. Sofia Olhede was Beit Scientific Research Fellow, recipient of an EPSRC (UK) grant, and beneficiary of the Helge Ax:son Johnsons Stiftelse. The authors wish to thank Curtis Condon, Ken White and Al Feng of the Beckman Center at the University of Illinois for the bat data and for permission to use it in this paper. The authors thank the referees for many helpful comments leading to an improved exposition. Appendix A. Here we shall prove result (7.). From Percival and Walden (, p. 7), where H j ( k N But Ẽ j,t = N N k= H ( k ) j N Ξk e iπtk/n, ) Lj = l= h j,l e iπlk/n and Ξ k = N t= ɛ te iπtk/n. Then, E{Ẽj,tẼ j,t } = E { E{Ξ k Ξ l } = N N k= N = σ ɛ N l= N u= v= N u= H ( k ) ( j N l ) } Hj N Ξk Ξ l e iπ(tk t l)/n. E {ɛ u ɛ v} e iπ(uk vl)/n e iπu(k l)/n = { Nσɛ, if k = l,, otherwise. Proc. R. Soc. Lond. A., 46, , 4
20 S. Olhede and A.T. Walden Hence, E{Ẽj,tẼ j,t } = σ ɛ N However, so that / σɛ / = σ ɛ N N k= N k= H j (f) N l= H j (f) 4 e iπf(t t ) df H ( k ) ( j N l ) Hj N δk,l e iπ(tk t l)/n H ( k ) 4 j N e iπk(t t )/N σɛ / / {, if f { [ ) (, j ]} j+, j+ j, otherwise, = H j (f) 4 e iπf(t t ) df. σɛ [ iπ(t t e iπf(t t ) /j+ +e iπf(t t ) ) / j { σ ɛ / j, t = t, σɛ {sin[π(t t )/ j ] sin[π(t t )/ j ]} [π(t t )], t t. /j / j+ and since E{Ẽj,t} =, we have cov{ẽj,t, Ẽj,t } = E{Ẽj,tẼ j,t }, and the result is proven. References Boashash, B. 99a Estimating and interpreting the instantaneous of a signal- Part : Fundamentals, Proc. IEEE, 8, Boashash, B. 99b Estimating and interpreting the instantaneous of a signal- Part : Algorithms and Applications, Proc. IEEE, 8, Burden & Faires. 993 Numerical analysis, 5th Edition. Boston MA: PWS Publishing Company. Coifman, R. R. and Donoho, D. L. 995 Translation-invariant denoising. In Wavelets and Statistics (Lecture Notes in Statistics, Volume 3), (A. Antoniadis and G. Oppenheim, eds.), 5 5, New York: Springer-Verlag. Cohen, L. 995 Time- analysis. Englewood Cliffs NJ: Prentice-Hall. Hess-Nielsen, N. & Wickerhauser, M. V. 996 Wavelets and - analysis, Proc. IEEE, 84, Huang, N. E., Shen, Z., Long, S. R., Wu, M. C., Shih, H. H., Zheng, Q., Yen, N-C., Tung, C. C. & Liu, H. H. 998 The empirical mode decomposition and Hilbert spectrum for nonlinear and non-stationary series analysis, Proc. R. Soc. Lond. A, 454, Huang, N. E., Shen, Z. & Long, S. R. 999 A new view of nonlinear water waves: the Hilbert spectrum, Annu. Rev. Fluid Mech., 3, Liang, J. & Parks, T. W. 996 A translation-invariant wavelet representation algorithm with applications, IEEE Trans. Signal Processing, 44, 5 3. Loughlin, P. J. & Tacer, B. 997 Comments on the interpretation of instantaneous, IEEE Signal Processing Letters, 4, 3 5. Magrin-Chagnolleau, I. & Baraniuk, R. G. 999 Empirical mode decomposition based - attributes. In Proc. 69th SEG Meeting, Houston. Soc. Expl. Geophys. ] Proc. R. Soc. Lond. A., 46, , 4
21 The Hilbert spectrum via wavelet projections Nielsen, M. On the construction and localization of finite orthogonal quadrature filters, J. Approx. Theory, 8, Olhede, S. & Walden, A. T. 3 Analytic wavelet thresholding. Preprint. Papoulis, A. 97 Probability, random variables, and stochastic processes. McGraw-Hill. Percival, D. B. & Walden, A. T. Wavelet methods for series analysis. Cambridge UK: Cambridge University Press. Sardy, S. Minimax threshold for denoising complex signals with Waveshrink, IEEE Trans. Signal Processing, 48, 3 8. Rilling, G., Flandrin, P. & Gonçalvès, P. (3) On empirical mode decomposition and its algorithms. Presented at IEEE-EURASIP Workshop on Nonlinear Signal and Image Processing NSIP-3, Grado, Italy, June 8. Shensa, M. J. 99 The discrete wavelet transform: wedding the à trous and Mallat algorithms, IEEE Transactions on Signal Processing, 4, Walden, A. T. & Contreras Cristán, A. 998 The phase-corrected undecimated discrete wavelet packet transform and its application to interpreting the timing of events. Proc. R. Soc. Lond. A, 454, Proc. R. Soc. Lond. A., 46, , 4
Wavelet Analysis of Discrete Time Series
Wavelet Analysis of Discrete Time Series Andrew T. Walden Abstract. We give a brief review of some of the wavelet-based techniques currently available for the analysis of arbitrary-length discrete time
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 informationStudy of nonlinear phenomena in a tokamak plasma using a novel Hilbert transform technique
Study of nonlinear phenomena in a tokamak plasma using a novel Hilbert transform technique Daniel Raju, R. Jha and A. Sen Institute for Plasma Research, Bhat, Gandhinagar-382428, INDIA Abstract. A new
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 informationDefining the Discrete Wavelet Transform (DWT)
Defining the Discrete Wavelet Transform (DWT) can formulate DWT via elegant pyramid algorithm defines W for non-haar wavelets (consistent with Haar) computes W = WX using O(N) multiplications brute force
More informationA study of the characteristics of white noise using the empirical mode decomposition method
1.198/rspa.3.11 A study of the characteristics of white noise using the empirical mode decomposition method By Zhaohua W u 1 and Norden E. Huang 1 Center for Ocean Land Atmosphere Studies, Suite 3, 441
More informationA new optimization based approach to the. empirical mode decomposition, adaptive
Annals of the University of Bucharest (mathematical series) 4 (LXII) (3), 9 39 A new optimization based approach to the empirical mode decomposition Basarab Matei and Sylvain Meignen Abstract - In this
More informationWavelet Methods for Time Series Analysis
Wavelet Methods for Time Series Analysis Donald B. Percival UNIVERSITY OF WASHINGTON, SEATTLE Andrew T. Walden IMPERIAL COLLEGE OF SCIENCE, TECHNOLOGY AND MEDICINE, LONDON CAMBRIDGE UNIVERSITY PRESS Contents
More informationQuadrature Prefilters for the Discrete Wavelet Transform. Bruce R. Johnson. James L. Kinsey. Abstract
Quadrature Prefilters for the Discrete Wavelet Transform Bruce R. Johnson James L. Kinsey Abstract Discrepancies between the Discrete Wavelet Transform and the coefficients of the Wavelet Series are known
More informationEMD ALGORITHM BASED ON BANDWIDTH AND THE APPLICATION ON ONE ECONOMIC DATA ANALYSIS
15th European Signal Processing Conference (EUSIPCO 7), Poznan, Poland, September 3-7, 7, copyright by EURASIP EMD ALGORITHM BASED ON BANDWIDTH AND THE APPLICATION ON ONE ECONOMIC DATA ANALYSIS Xie Qiwei
More informationQuantitative Finance II Lecture 10
Quantitative Finance II Lecture 10 Wavelets III - Applications Lukas Vacha IES FSV UK May 2016 Outline Discrete wavelet transformations: MODWT Wavelet coherence: daily financial data FTSE, DAX, PX Wavelet
More informationMode Decomposition Analysis Applied to Study the Low-Frequency Embedded in the Vortex Shedding Process. Abstract
Tainan,Taiwan,R.O.C., -3 December 3 Mode Decomposition Analysis Applied to Study the Low-Frequency Embedded in the Vortex Shedding Process Chin-Tsan Wang Department of Electrical Engineering Kau Yuan Institute
More informationEmpirical Mode Decomposition of Financial Data
International Mathematical Forum, 3, 28, no. 25, 1191-122 Empirical Mode Decomposition of Financial Data Konstantinos Drakakis 1 UCD CASL 2 University College Dublin Abstract We propose a novel method
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 informationA Tutorial on Stochastic Models and Statistical Analysis for Frequency Stability Measurements
A Tutorial on Stochastic Models and Statistical Analysis for Frequency Stability Measurements Don Percival Applied Physics Lab, University of Washington, Seattle overheads for talk available at http://staff.washington.edu/dbp/talks.html
More informationWavelet Methods for Time Series Analysis. Part IV: Wavelet-Based Decorrelation of Time Series
Wavelet Methods for Time Series Analysis Part IV: Wavelet-Based Decorrelation of Time Series DWT well-suited for decorrelating certain time series, including ones generated from a fractionally differenced
More informationEnhanced Active Power Filter Control for Nonlinear Non- Stationary Reactive Power Compensation
Enhanced Active Power Filter Control for Nonlinear Non- Stationary Reactive Power Compensation Phen Chiak See, Vin Cent Tai, Marta Molinas, Kjetil Uhlen, and Olav Bjarte Fosso Norwegian University of Science
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 informationON THE FILTERING PROPERTIES OF THE EMPIRICAL MODE DECOMPOSITION
Advances in Adaptive Data Analysis Vol. 2, No. 4 (2010) 397 414 c World Scientific Publishing Company DOI: 10.1142/S1793536910000604 ON THE FILTERING PROPERTIES OF THE EMPIRICAL MODE DECOMPOSITION ZHAOHUA
More informationOne or two Frequencies? The empirical Mode Decomposition Answers.
One or two Frequencies? The empirical Mode Decomposition Answers. Gabriel Rilling, Patrick Flandrin To cite this version: Gabriel Rilling, Patrick Flandrin. One or two Frequencies? The empirical Mode Decomposition
More informationUltrasonic Thickness Inspection of Oil Pipeline Based on Marginal Spectrum. of Hilbert-Huang Transform
17th World Conference on Nondestructive Testing, 25-28 Oct 2008, Shanghai, China Ultrasonic Thickness Inspection of Oil Pipeline Based on Marginal Spectrum of Hilbert-Huang Transform Yimei MAO, Peiwen
More informationDigital Image Processing
Digital Image Processing Wavelets and Multiresolution Processing (Wavelet Transforms) Christophoros Nikou cnikou@cs.uoi.gr University of Ioannina - Department of Computer Science 2 Contents Image pyramids
More informationUsing modern time series analysis techniques to predict ENSO events from the SOI time series
Nonlinear Processes in Geophysics () 9: 4 45 Nonlinear Processes in Geophysics c European Geophysical Society Using modern time series analysis techniques to predict ENSO events from the SOI time series
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 informationMultiscale Characterization of Bathymetric Images by Empirical Mode Decomposition
Multiscale Characterization of Bathymetric Images by Empirical Mode Decomposition El-Hadji Diop & A.O. Boudraa, IRENav, Ecole Navale, Groupe ASM, Lanvéoc Poulmic, BP600, 29240 Brest-Armées, France A. Khenchaf
More informationWavelet Methods for Time Series Analysis. Quick Comparison of the MODWT to the DWT
Wavelet Methods for Time Series Analysis Part IV: MODWT and Examples of DWT/MODWT Analysis MODWT stands for maximal overlap discrete wavelet transform (pronounced mod WT ) transforms very similar to the
More informationWavelet Methods for Time Series Analysis. Quick Comparison of the MODWT to the DWT
Wavelet Methods for Time Series Analysis Part III: MODWT and Examples of DWT/MODWT Analysis MODWT stands for maximal overlap discrete wavelet transform (pronounced mod WT ) transforms very similar to the
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 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 informationComputational Data Analysis!
12.714 Computational Data Analysis! Alan Chave (alan@whoi.edu)! Thomas Herring (tah@mit.edu),! http://geoweb.mit.edu/~tah/12.714! Introduction to Spectral Analysis! Topics Today! Aspects of Time series
More informationAN ALTERNATIVE ALGORITHM FOR EMPIRICAL MODE DECOMPOSITION. 1. Introduction
AN ALTERNATIVE ALGORITHM FOR EMPIRICAL MODE DECOMPOSITION LUAN LIN, YANG WANG, AND HAOMIN ZHOU Abstract. The empirical mode decomposition (EMD) was a method pioneered by Huang et al [8] as an alternative
More informationEnsemble empirical mode decomposition of Australian monthly rainfall and temperature data
19th International Congress on Modelling and Simulation, Perth, Australia, 12 16 December 2011 http://mssanz.org.au/modsim2011 Ensemble empirical mode decomposition of Australian monthly rainfall and temperature
More informationAn Introduction to Wavelets and some Applications
An Introduction to Wavelets and some Applications Milan, May 2003 Anestis Antoniadis Laboratoire IMAG-LMC University Joseph Fourier Grenoble, France An Introduction to Wavelets and some Applications p.1/54
More informationWavelet Methods for Time Series Analysis. Motivating Question
Wavelet Methods for Time Series Analysis Part VII: Wavelet-Based Bootstrapping start with some background on bootstrapping and its rationale describe adjustments to the bootstrap that allow it to work
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 informationDiscrete Simulation of Power Law Noise
Discrete Simulation of Power Law Noise Neil Ashby 1,2 1 University of Colorado, Boulder, CO 80309-0390 USA 2 National Institute of Standards and Technology, Boulder, CO 80305 USA ashby@boulder.nist.gov
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 informationEMD-BASED STOCHASTIC SUBSPACE IDENTIFICATION OF CIVIL ENGINEERING STRUCTURES UNDER OPERATIONAL CONDITIONS
EMD-BASED STOCHASTIC SUBSPACE IDENTIFICATION OF CIVIL ENGINEERING STRUCTURES UNDER OPERATIONAL CONDITIONS Wei-Xin Ren, Department of Civil Engineering, Fuzhou University, P. R. China Dan-Jiang Yu Department
More informationBSc Project Fault Detection & Diagnosis in Control Valve
BSc Project Fault Detection & Diagnosis in Control Valve Supervisor: Dr. Nobakhti Content 2 What is fault? Why we detect fault in a control loop? What is Stiction? Comparing stiction with other faults
More informationComparing the Time-Deformation Method with the Fractional Fourier Transform in Filtering Non-Stationary Processes
Journal of Signal and Information Processing, 1, 3, 491-51 http://dx.doi.org/1.436/jsip.1.346 Published Online November 1 (http://www.scirp.org/journal/jsip) Comparing the Time-Deformation Method with
More informationω 0 = 2π/T 0 is called the fundamental angular frequency and ω 2 = 2ω 0 is called the
he ime-frequency Concept []. Review of Fourier Series Consider the following set of time functions {3A sin t, A sin t}. We can represent these functions in different ways by plotting the amplitude versus
More informationA WAVELET BASED CODING SCHEME VIA ATOMIC APPROXIMATION AND ADAPTIVE SAMPLING OF THE LOWEST FREQUENCY BAND
A WAVELET BASED CODING SCHEME VIA ATOMIC APPROXIMATION AND ADAPTIVE SAMPLING OF THE LOWEST FREQUENCY BAND V. Bruni, D. Vitulano Istituto per le Applicazioni del Calcolo M. Picone, C. N. R. Viale del Policlinico
More informationDigital Image Processing
Digital Image Processing, 2nd ed. Digital Image Processing Chapter 7 Wavelets and Multiresolution Processing Dr. Kai Shuang Department of Electronic Engineering China University of Petroleum shuangkai@cup.edu.cn
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 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 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 informationIMPROVEMENTS IN MODAL PARAMETER EXTRACTION THROUGH POST-PROCESSING FREQUENCY RESPONSE FUNCTION ESTIMATES
IMPROVEMENTS IN MODAL PARAMETER EXTRACTION THROUGH POST-PROCESSING FREQUENCY RESPONSE FUNCTION ESTIMATES Bere M. Gur Prof. Christopher Niezreci Prof. Peter Avitabile Structural Dynamics and Acoustic Systems
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 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 informationLOCALIZED analyses in one dimension(al) (1-D) have
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 55, NO. 3, MARCH 2007 921 Multiple Multidimensional Morse Wavelets Georgios Metikas Sofia C. Olhede Abstract This paper defines a set of operators that localize
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 informationWavelet Methods for Time Series Analysis. Part IX: Wavelet-Based Bootstrapping
Wavelet Methods for Time Series Analysis Part IX: Wavelet-Based Bootstrapping start with some background on bootstrapping and its rationale describe adjustments to the bootstrap that allow it to work with
More informationHilbert-Huang Transform-based Local Regions Descriptors
Hilbert-Huang Transform-based Local Regions Descriptors Dongfeng Han, Wenhui Li, Wu Guo Computer Science and Technology, Key Laboratory of Symbol Computation and Knowledge Engineering of the Ministry of
More informationDerivative-Optimized Empirical Mode Decomposition (DEMD) for the Hilbert- Huang Transform
Derivative-Optimized Empirical Mode Decomposition (DEMD) for the Hilbert- Huang Transform Peter C. Chu 1), Chenwu Fan 1), and Norden Huang 2) 1)Naval Postgraduate School Monterey, California, USA 2)National
More informationA new optimization based approach to the empirical mode decomposition
A new optimization based approach to the empirical mode decomposition Basarab Matei, Sylvain Meignen To cite this version: Basarab Matei, Sylvain Meignen. A new optimization based approach to the empirical
More informationMODWT Based Time Scale Decomposition Analysis. of BSE and NSE Indexes Financial Time Series
Int. Journal of Math. Analysis, Vol. 5, 211, no. 27, 1343-1352 MODWT Based Time Scale Decomposition Analysis of BSE and NSE Indexes Financial Time Series Anu Kumar 1* and Loesh K. Joshi 2 Department of
More informationarxiv:gr-qc/ v1 24 Dec 2004
LIGO-P15133 Multiresolution techniques for the detection of gravitational-wave bursts arxiv:gr-qc/412119v1 24 Dec 24 S Chatterji, L Blackburn, G Martin and E Katsavounidis MIT LIGO Laboratory, NW17-161,
More informationRepresentation of Signals and Systems. Lecturer: David Shiung
Representation of Signals and Systems Lecturer: David Shiung 1 Abstract (1/2) Fourier analysis Properties of the Fourier transform Dirac delta function Fourier transform of periodic signals Fourier-transform
More information5 Analog carrier modulation with noise
5 Analog carrier modulation with noise 5. Noisy receiver model Assume that the modulated signal x(t) is passed through an additive White Gaussian noise channel. A noisy receiver model is illustrated in
More informationIntroduction to Wavelets and Wavelet Transforms
Introduction to Wavelets and Wavelet Transforms A Primer C. Sidney Burrus, Ramesh A. Gopinath, and Haitao Guo with additional material and programs by Jan E. Odegard and Ivan W. Selesnick Electrical and
More informationAn Introduction to Wavelets with Applications in Environmental Science
An Introduction to Wavelets with Applications in Environmental Science Don Percival Applied Physics Lab, University of Washington Data Analysis Products Division, MathSoft overheads for talk available
More informationReview of Concepts from Fourier & Filtering Theory. Fourier theory for finite sequences. convolution/filtering of infinite sequences filter cascades
Review of Concepts from Fourier & Filtering Theory precise definition of DWT requires a few basic concepts from Fourier analysis and theory of linear filters will start with discussion/review of: basic
More informationA Tutorial on Stochastic Models and Statistical Analysis for Frequency Stability Measurements
A Tutorial on Stochastic Models and Statistical Analysis for Frequency Stability Measurements Don Percival Applied Physics Lab, University of Washington, Seattle overheads for talk available at http://staffwashingtonedu/dbp/talkshtml
More informationAn Adaptive Data Analysis Method for nonlinear and Nonstationary Time Series: The Empirical Mode Decomposition and Hilbert Spectral Analysis
An Adaptive Data Analysis Method for nonlinear and Nonstationary Time Series: The Empirical Mode Decomposition and Hilbert Spectral Analysis Norden E Huang and Zhaohua Wu Abstract An adaptive data analysis
More informationTHE GENERALIZATION OF DISCRETE STOCKWELL TRANSFORMS
19th European Signal Processing Conference (EUSIPCO 2011) Barcelona, Spain, August 29 - September 2, 2011 THE GEERALIZATIO OF DISCRETE STOCKWELL TRASFORMS Yusong Yan, Hongmei Zhu Department of Mathematics
More informationNoise reduction of ship-radiated noise based on noise-assisted bivariate empirical mode decomposition
Indian Journal of Geo-Marine Sciences Vol. 45(4), April 26, pp. 469-476 Noise reduction of ship-radiated noise based on noise-assisted bivariate empirical mode decomposition Guohui Li,2,*, Yaan Li 2 &
More informationMultiresolution techniques for the detection of gravitational-wave bursts
INSTITUTE OF PHYSICS PUBLISHING Class. Quantum Grav. 21 (24) S189 S1818 CLASSICAL AND QUANTUM GRAVITY PII: S264-9381(4)8352-6 Multiresolution techniques for the detection of gravitational-wave bursts S
More informationSPECTRAL METHODS ASSESSMENT IN JOURNAL BEARING FAULT DETECTION APPLICATIONS
The 3 rd International Conference on DIAGNOSIS AND PREDICTION IN MECHANICAL ENGINEERING SYSTEMS DIPRE 12 SPECTRAL METHODS ASSESSMENT IN JOURNAL BEARING FAULT DETECTION APPLICATIONS 1) Ioannis TSIAFIS 1),
More informationUsing frequency domain techniques with US real GNP data: A Tour D Horizon
Using frequency domain techniques with US real GNP data: A Tour D Horizon Patrick M. Crowley TAMUCC October 2011 Patrick M. Crowley (TAMUCC) Bank of Finland October 2011 1 / 45 Introduction "The existence
More informationMULTI-SCALE IMAGE DENOISING BASED ON GOODNESS OF FIT (GOF) TESTS
MULTI-SCALE IMAGE DENOISING BASED ON GOODNESS OF FIT (GOF) TESTS Naveed ur Rehman 1, Khuram Naveed 1, Shoaib Ehsan 2, Klaus McDonald-Maier 2 1 Department of Electrical Engineering, COMSATS Institute of
More informationAdaptive Decomposition Into Mono-Components
Adaptive Decomposition Into Mono-Components Tao Qian, Yan-Bo Wang and Pei Dang Abstract The paper reviews some recent progress on adaptive signal decomposition into mono-components that are defined to
More informationApplication of Hilbert-Huang signal processing to ultrasonic non-destructive testing of oil pipelines *
13 Mao et al. / J Zhejiang Univ SCIENCE A 26 7(2):13-134 Journal of Zhejiang University SCIENCE A ISSN 19-395 http://www.zju.edu.cn/jzus E-mail: jzus@zju.edu.cn Application of Hilbert-Huang signal processing
More informationMULTI-ORDERED POSETS. Lisa Bishop Department of Mathematics, Occidental College, Los Angeles, CA 90041, United States.
INTEGERS: ELECTRONIC JOURNAL OF COMBINATORIAL NUMBER THEORY 7 (2007), #A06 MULTI-ORDERED POSETS Lisa Bishop Department of Mathematics, Occidental College, Los Angeles, CA 90041, United States lbishop@oxy.edu
More informationFourier spectra of radially periodic images with a non-symmetric radial period
J. Opt. A: Pure Appl. Opt. 1 (1999) 621 625. Printed in the UK PII: S1464-4258(99)00768-0 Fourier spectra of radially periodic images with a non-symmetric radial period Isaac Amidror Laboratoire de Systèmes
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 informationON EMPIRICAL MODE DECOMPOSITION AND ITS ALGORITHMS
ON EMPIRICAL MODE DECOMPOSITION AND ITS ALGORITHMS Gabriel Rilling, Patrick Flandrin and Paulo Gonçalvès Laboratoire de Physique (UMR CNRS 5672), École Normale Supérieure de Lyon 46, allée d Italie 69364
More informationIntroduction to Wavelet. Based on A. Mukherjee s lecture notes
Introduction to Wavelet Based on A. Mukherjee s lecture notes Contents History of Wavelet Problems of Fourier Transform Uncertainty Principle The Short-time Fourier Transform Continuous Wavelet Transform
More informationThe use of the Hilbert transforra in EMG Analysis
The use of the Hilbert transforra in EMG Analysis Sean Tallier BSc, Dr P.Kyberd Oxford Orthopaedic Engineering Center Oxford University 1.1 INTRODUCTION The Fourier transform has traditionally been used
More informationCommunication Theory II
Communication Theory II Lecture 4: Review on Fourier analysis and probabilty theory Ahmed Elnakib, PhD Assistant Professor, Mansoura University, Egypt Febraury 19 th, 2015 1 Course Website o http://lms.mans.edu.eg/eng/
More informationEXAMINATION OF RELATIONSHIPS BETWEEN TIME SERIES BY WAVELETS: ILLUSTRATION WITH CZECH STOCK AND INDUSTRIAL PRODUCTION DATA
EXAMINATION OF RELATIONSHIPS BETWEEN TIME SERIES BY WAVELETS: ILLUSTRATION WITH CZECH STOCK AND INDUSTRIAL PRODUCTION DATA MILAN BAŠTA University of Economics, Prague, Faculty of Informatics and Statistics,
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 informationIEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 54, NO. 2, FEBRUARY
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL 54, NO 2, FEBRUARY 2006 423 Underdetermined Blind Source Separation Based on Sparse Representation Yuanqing Li, Shun-Ichi Amari, Fellow, IEEE, Andrzej Cichocki,
More informationElec4621 Advanced Digital Signal Processing Chapter 11: Time-Frequency Analysis
Elec461 Advanced Digital Signal Processing Chapter 11: Time-Frequency Analysis Dr. D. S. Taubman May 3, 011 In this last chapter of your notes, we are interested in the problem of nding the instantaneous
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 informationIn the Name of God. Lectures 15&16: Radial Basis Function Networks
1 In the Name of God Lectures 15&16: Radial Basis Function Networks Some Historical Notes Learning is equivalent to finding a surface in a multidimensional space that provides a best fit to the training
More informationREVIEW OF SINGLE CHANNEL SOURCE SEPARATION TECHNIQUES
REVIEW OF SINGLE CHANNEL SOURCE SEPARATION TECHNIQUES Kedar Patki University of Rochester Dept. of Electrical and Computer Engineering kedar.patki@rochester.edu ABSTRACT The paper reviews the problem of
More informationArticle (peer-reviewed)
Title Author(s) Influence of noise intensity on the spectrum of an oscillator Swain, Rabi Sankar; Gleeson, James P.; Kennedy, Michael Peter Publication date 2005-11 Original citation Type of publication
More informationSignal Period Analysis Based on Hilbert-Huang Transform and Its Application to Texture Analysis
Signal Period Analysis Based on Hilbert-Huang Transform and Its Application to Texture Analysis Zhihua Yang 1, Dongxu Qi and Lihua Yang 3 1 School of Information Science and Technology Sun Yat-sen University,
More informationA Spectral Approach for Sifting Process in Empirical Mode Decomposition
5612 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 58, NO. 11, NOVEMBER 2010 A Spectral Approach for Sifting Process in Empirical Mode Decomposition Oumar Niang, Éric Deléchelle, and Jacques Lemoine Abstract
More informationWavelet-Based Bootstrapping for Non-Gaussian Time Series. Don Percival
Wavelet-Based Bootstrapping for Non-Gaussian Time Series Don Percival Applied Physics Laboratory Department of Statistics University of Washington Seattle, Washington, USA overheads for talk available
More informationAnalytic discrete cosine harmonic wavelet transform(adchwt) and its application to signal/image denoising
Analytic discrete cosine harmonic wavelet transform(adchwt) and its application to signal/image denoising M. Shivamurti and S. V. Narasimhan Digital signal processing and Systems Group Aerospace Electronic
More informationANALYSIS OF TEMPORAL VARIATIONS IN TURBIDITY FOR A COASTAL AREA USING THE HILBERT-HUANG-TRANSFORM
ANALYSIS OF TEMPORAL VARIATIONS IN TURBIDITY FOR A COASTAL AREA USING THE HILBERT-HUANG-TRANSFORM Shigeru Kato, Magnus Larson 2, Takumi Okabe 3 and Shin-ichi Aoki 4 Turbidity data obtained by field observations
More informationDigital Image Processing
Digital Image Processing Wavelets and Multiresolution Processing () Christophoros Nikou cnikou@cs.uoi.gr University of Ioannina - Department of Computer Science 2 Contents Image pyramids Subband coding
More informationExamples of DWT & MODWT Analysis: Overview
Examples of DWT & MODWT Analysis: Overview look at DWT analysis of electrocardiogram (ECG) data discuss potential alignment problems with the DWT and how they are alleviated with the MODWT look at MODWT
More informationOPTIMAL PERTURBATION OF UNCERTAIN SYSTEMS
Stochastics and Dynamics c World Scientific Publishing Company OPTIMAL PERTURBATION OF UNCERTAIN SYSTEMS BRIAN F. FARRELL Division of Engineering and Applied Sciences, Harvard University Pierce Hall, 29
More informationOPTIMAL SURE PARAMETERS FOR SIGMOIDAL WAVELET SHRINKAGE
17th European Signal Processing Conference (EUSIPCO 009) Glasgow, Scotland, August 4-8, 009 OPTIMAL SURE PARAMETERS FOR SIGMOIDAL WAVELET SHRINKAGE Abdourrahmane M. Atto 1, Dominique Pastor, Gregoire Mercier
More informationDetection of Anomalous Observations Using Wavelet Analysis in Frequency Domain
International Journal of Scientific and Research Publications, Volume 7, Issue 8, August 2017 76 Detection of Anomalous Observations Using Wavelet Analysis in Frequency Domain Aideyan D.O. Dept of Mathematical
More informationLogarithmic quantisation of wavelet coefficients for improved texture classification performance
Logarithmic quantisation of wavelet coefficients for improved texture classification performance Author Busch, Andrew, W. Boles, Wageeh, Sridharan, Sridha Published 2004 Conference Title 2004 IEEE International
More informationTime Series. Anthony Davison. c
Series Anthony Davison c 2008 http://stat.epfl.ch Periodogram 76 Motivation............................................................ 77 Lutenizing hormone data..................................................
More informationWavelet Footprints: Theory, Algorithms, and Applications
1306 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 51, NO. 5, MAY 2003 Wavelet Footprints: Theory, Algorithms, and Applications Pier Luigi Dragotti, Member, IEEE, and Martin Vetterli, Fellow, IEEE Abstract
More information