a confidence limit for hilbert spectrum through stoppage criteria
Post on 20-Dec-2015
222 views
TRANSCRIPT
![Page 1: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/1.jpg)
A Confidence Limit for Hilbert Spectrum
Through stoppage criteria
![Page 2: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/2.jpg)
Need a confidence limit
• As we have presented here, EMD could generate infinite many sets of IMFs. In this case, which one of the infinite many sets really represents the true physics? To answer this question, we need a confidence limit on our result.
• Traditional methods also had the similar problem of generating many answers. Take Fourier analysis for example; we have to assume trigonometric series is basis. How about other basis? Why not slightly distorted sinusoidal wave as basis? ….
![Page 3: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/3.jpg)
Confidence Limit for Fourier Spectrum
• The Confidence limit for Fourier Spectral analysis is based on ergodic assumption.
• It is derived by dividing the data into M sections, and substituting temporal (or spatial) average as ensemble average.
• This approach is valid for linear and stationary processes, and the sub-sections have to be statistically independent.
• By dividing the data into subsections, the resolution will suffer.
![Page 4: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/4.jpg)
Statistical Independence
The probability function of x and y jointly, f(x, y), is equal to f(x) times f(y) if x and y are statistically independent. For any number of variables, x1, x2, ...,
xn, if the joint probability is the product of the several
probability functions, then the variables are all statistically independent. Independent variables are noncorrelated, but not necessarily conversely.
James & James : Mathematics Dictionary
![Page 5: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/5.jpg)
LOD Data
![Page 6: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/6.jpg)
Confidence Limit for Fourier Spectrum
Confidence Limit from 7 sections, each 2048 points.
![Page 7: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/7.jpg)
Are the sub-sections statistically independent?
For narrow band signals, most likely they would not be
independent.
![Page 8: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/8.jpg)
Confidence Limit for Hilbert Spectrum
• Any data can be decomposed into infinitely many different constituting component sets.
• EMD is a method to generate infinitely many different IMF representations based on different sifting parameters.
• Some of the IMFs are better than others based on various properties: for example, Orthogonal Index.
• A Confidence Limit for Hilbert Spectral analysis can be based on an ensemble of ‘valid IMF’ resulting from different sifting parameters S covering the parameter space fairly.
• It is valid for nonlinear and nonstationary processes.
![Page 9: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/9.jpg)
Different Kinds of Confidence Limit
The basic idea is to generate various IMFs, treat the mean as the true answer, and obtain the confidence limit based on the STD from the various solutions.
• By stoppage criteria • By Ensemble EMD• By down sampling• Different spline methods
![Page 10: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/10.jpg)
Empirical Mode DecompositionSifting : to get one IMF component
1 1
1 2 2
k 1 k k
k 1
x( t ) m h ,
h m h ,
.....
.....
h m h
h c .
.
![Page 11: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/11.jpg)
None of the above methods depends on ergodic assumption.
We are truly achieving an ensemble mean.
![Page 12: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/12.jpg)
The Stoppage Criteria : S and SDA. The S number : S is defined as the consecutive
number of siftings, in which the numbers of zero-crossing and extrema are the same for these S siftings.
B. If the mean is smaller than a pre-assigned value.
C. Fixed sifting (iterating) time.
D. SD is small than a pre-set value, whereT
2
k 1 kt 0
T2
k 1t 0
h ( t ) h ( t )SD
h ( t )
![Page 13: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/13.jpg)
Critical Parameters for EMD
• The maximum number of sifting allowed to extract an IMF, N.– Note: N is originally set to guarantee convergence
of sifting, but later found to be superfluous.
• The criterion for accepting a sifting component as an IMF, the Stoppage criterion S.
• Therefore, the nomenclature for the IMF are
CE(N, S) : for extrema sifting
CC(N, S) : for curvature sifting
![Page 14: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/14.jpg)
Sifting with Intermittence Test
• To avoid mode mixing, we have to institute a special criterion to separate oscillation of different time scales into different IMF components.
• The criteria is to select time scale so that oscillations with time scale shorter than this pre-selected criterion is not included in the IMF.
![Page 15: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/15.jpg)
Intermittence Sifting : Data
![Page 16: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/16.jpg)
Intermittence Sifting : IMF
![Page 17: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/17.jpg)
Intermittence Sifting : Hilbert Spectra
![Page 18: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/18.jpg)
Intermittence Sifting : Hilbert Spectra (Low)
![Page 19: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/19.jpg)
Intermittence Sifting : Marginal Spectra
![Page 20: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/20.jpg)
Intermittence Sifting : Marginal spectra (Low)
![Page 21: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/21.jpg)
Intermittence Sifting : Marginal spectra (High)
![Page 22: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/22.jpg)
Critical Parameters for Sifting
• Because of the inclusion of intermittence test there will be one set of intermittence criteria.
• Therefore, the Nomenclature for IMF here are
CEI(N, S: n1, n2, …)
CCI(N, S: n1, n2, …)
with n1, n2 as the intermittence test criteria.
Note: N is originally set to guarantee convergence of
sifting, but later found to be superfluous.
![Page 23: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/23.jpg)
Effects of EMD (Sifting)
• To separate data into components of similar scale.
• To eliminate ridding waves.• To make the results symmetric with respect to
the x-axis and the amplitude more even.
– Note: The first two are necessary for valid IMF, the last effect actually cause the IMF to lost its intrinsic properties.
![Page 24: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/24.jpg)
LOD Data
![Page 25: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/25.jpg)
IMF CE(100, 2)
![Page 26: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/26.jpg)
IMF CE(100, 10)
![Page 27: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/27.jpg)
Orthogonal Index as function of N and S Contour
![Page 28: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/28.jpg)
Orthogonality Index as function of N and S
![Page 29: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/29.jpg)
Confidence Limit without Intermittence Criteria
Number of IMF for different siftings may not be the same; therefore, average of IMF is, in general, not possible. However, we can take the mean of the Hilbert Spectra, for we can make all the spectra having the same frequency and time ranges.
![Page 30: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/30.jpg)
Hilbert Spectrum CE(100, 2)
![Page 31: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/31.jpg)
Mean Hilbert Spectrum : All CEs
![Page 32: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/32.jpg)
STD Hilbert Spectrum : All CEs
![Page 33: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/33.jpg)
Marginal Mean & STD Hilbert Spectra : All CEs
![Page 34: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/34.jpg)
Mean and STD of Marginal Hilbert Spectra
![Page 35: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/35.jpg)
Confidence Limit with Intermittence Criteria
In general, the number of IMFs can be controlled to the same; therefore, averages of IMFs and Hilbert Spectra are all possible.
![Page 36: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/36.jpg)
IMF CEI(100,2; 4,-1^3,45^2,-10)
![Page 37: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/37.jpg)
IMF CEI(100,10; 4,-1^3,45^2,-10)
![Page 38: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/38.jpg)
Envelopes of Selected Annual Cycle IMFs
![Page 39: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/39.jpg)
Orthogonal Indices for CEI cases
![Page 40: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/40.jpg)
IMF : Mean CEI 9 cases
![Page 41: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/41.jpg)
IMF : STD CEI 9 cases
![Page 42: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/42.jpg)
Mean Hilbert Spectrum : All CEIs
![Page 43: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/43.jpg)
STD Hilbert Spectrum for All CEIs
![Page 44: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/44.jpg)
Marginal mean & STD Hilbert Spectra : All CEIs
![Page 45: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/45.jpg)
Mean Marginal Hilbert Spectrum & Confidence Limit :All CEIs
![Page 46: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/46.jpg)
Mean Marginal Hilbert Spectrum & Confidence Limit :All CEs
![Page 47: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/47.jpg)
Individual Annual Cycle IMFs : 9 CEI Cases
![Page 48: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/48.jpg)
Details of Individual Annual Cycle : CEIs
![Page 49: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/49.jpg)
Mean Annual Cycle & Envelope: 9 CEI Cases
![Page 50: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/50.jpg)
Individual Envelopes for Annual Cycle IMFs
![Page 51: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/51.jpg)
Mean Envelopes for Annual Cycle IMFs
![Page 52: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/52.jpg)
Optimal Sifting Parameters
• The Maximum sifting number should be set very high to guarantee that the stoppage criterion is always satisfied.
• The Stoppage criterion should be selected by considering the difference between the individual case with the mean to see if there is an optimal range where the difference is minimum.
• The difference can be computed from the Hilbert spectra or IMF components. It turn out that the IMF is a more sensitive way to determine the optimal sifting parameters.
![Page 53: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/53.jpg)
Computation of the Differences
N
jj 1
2
j jt
1V ( t ) V ( t ; S ) .
N
D V( t ) V ( t ; S )
where V(t) can be IMF or Hilbert Spectrum.
![Page 54: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/54.jpg)
IMF for CEI Cases : Annual Cycle
![Page 55: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/55.jpg)
IMF for CEI Cases : Half-monthly tidal Cycle
![Page 56: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/56.jpg)
Hilbert Spectrum : Deviation Individual form the mean CEI
![Page 57: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/57.jpg)
Hilbert Spectrum : Deviation Individual form the mean CE
![Page 58: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/58.jpg)
Another Example using Earthquake Data
Earthquake data has no fixed time scale; therefore, it is not possible to sift with intermittence. The only way to compute the confidence limit is use an ensemble of Hilbert Spectra.
![Page 59: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/59.jpg)
Earthquake Data
![Page 60: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/60.jpg)
Mean Hilbert Spectrum
![Page 61: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/61.jpg)
Another Example using Earthquake Data
![Page 62: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/62.jpg)
Optimal Selection of Stoppage Criterion
• From the above tests, we can see that the Hilbert spectrum difference is less sensitive to the changes of stoppage criterion S than IMFs.
• From the IMF tests, we suggest that the S number should be set in the range of 3 to 10.
• This selection is in agreement with our past experiences; however, additional quantitative tests should be conduct for other data types.
![Page 63: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/63.jpg)
Summary
• The Confidence limit presented here exists only with respect to the EMD method used.
• The Confidence limit presented here is only one of many possibilities. Instead of using OI as criterion, we can also use the STD of different trials to get a feeling of the stability of the analysis.
• Instead of Stoppage criteria, we can us different spline methods, down sampling and study their variations.
• Most interestingly, we could use Ensemble EMD, to be discussed next.
![Page 64: A Confidence Limit for Hilbert Spectrum Through stoppage criteria](https://reader036.vdocument.in/reader036/viewer/2022062421/56649d495503460f94a25a03/html5/thumbnails/64.jpg)
Envelope of IMF : c1