To view this presentation, you'll need to enable Flash.

Show me how

After you enable Flash, refresh this webpage and the presentation should play.

Loading...

PPT – The Hilbert Transform and Empirical Mode Decomposition: PowerPoint presentation | free to view - id: ab4f1-OTdhN

The Adobe Flash plugin is needed to view this content

View by Category

Presentations

Products
Sold on our sister site CrystalGraphics.com

About This Presentation

Write a Comment

User Comments (0)

Transcript and Presenter's Notes

The Hilbert Transform and Empirical Mode

Decomposition

Powerful Tools for Data Analysis

- Suz Tolwinski
- University of Arizona
- Program in Applied Mathematics
- Spring 2007 RTG

Deriving the Hilbert Transform

For f(z) u(x,y) iv(x,y) analytic on the upper

half-plane, and decays at infinity

Cauchys Integral Formula Decay of function

Clever rearrangement of terms

Relationship between u(x,y) and v(x,y) on R

Define the Hilbert transform in similar spirit

Another View of the Hilbert Transform.

Looks like minus the convolution of f(t) with

1/pt. Apply Convolution Theorem for something

more manageable in frequency space?

Take away message The Hilbert transform rotates

a signal counter-clockwise by p/2 at every point

in its positive frequency spectrum, and clockwise

at all its negative frequencies.

Real Signals.

A signal is any time-varying quantity containing

information.

Signals in nature are real.

Real signals have even frequency spectra.

This makes them difficult to analyze. We would

like to know, how is the signal energy

distributed in time and/or frequency space?

For a signal with even frequency spectrum,

-Mean frequency? -Spread of

signal in frequency space?

(Energy of s(t) ? s(t)2 dt ? S(?)2 d?)

Electrocardiographic signal in time and frequency

domain

Analytic Signals.

Have positive frequency spectrum only, so lt?gt,

spread in ? are meaningful quantities.

We can construct an analytic signal

corresponding to any real signal (takes only

the positive frequencies)

Frequency spectrum of a real signal.

FANALYT.(?) 1/2(F(?) sgn(?)F(?))

Then in the time domain, the analytic signal is

given by

Frequency spectrum of the corresponding analytic

signal.

fANALYT.(t) 1/2(f(t) iHf(t))

Analytic signal can be represented as

time-varying frequency and amplitude

fANALYT.(t) A(t)ei?(t)

A(t) (1/2f(t))21/2 (Hf(t))2)1/2

?(t) tan-1(Hf(t)/f(t))

Empirical Mode Decomposition of Real Signals.

(Method due to N. Huang, 1998.)

- Creates an adaptive decomposition of signal
- Result is a generalized Fourier series
- (Modes with time-varying amplitude and phase)
- Components are called Intrinsic Mode Functions

(IMFs) - IMFs satisfy two criteria by designs
- -Have only one zero between successive extrema
- -Have zero local mean
- EMD separates phenomena occurring on different

time scales. - Residue shows overall trend in data.

Huangs Sifting Process.

Residue s(t) I1(t) Residue i 1 k 1

while Residue not equal zero or not monotone

while Ii has non-negligible local

mean U(t) spline through local maxima of

Ii L(t) spline through local minima of

Ii Av(t) 1/2 (U(t) L(t)) Ii(t) Ii(t) -

Av(t) i i 1 end

IMFk(t) Ii(t) Residue Residue -

IMFk k k1 end

Huangs Sifting Process.

Residue s(t) I1(t) Residue i 1 k 1

while Residue not equal zero or not monotone

while Ii has non-negligible local

mean U(t) spline through local maxima of

Ii L(t) spline through local minima of

Ii Av(t) 1/2 (U(t) L(t)) Ii(t) Ii(t) -

Av(t) i i 1 end

IMFk(t) Ii(t) Residue Residue -

IMFk k k1 end

Huangs Sifting Process.

Residue s(t) I1(t) Residue i 1 k 1

while Residue not equal zero or not monotone

while Ii has non-negligible local

mean U(t) spline through local maxima of

Ii L(t) spline through local minima of

Ii Av(t) 1/2 (U(t) L(t)) Ii(t) Ii(t) -

Av(t) i i 1 end

IMFk(t) Ii(t) Residue Residue -

IMFk k k1 end

Huangs Sifting Process.

Residue s(t) I1(t) Residue i 1 k 1

while Residue not equal zero or not monotone

while Ii has non-negligible local

mean U(t) spline through local maxima of

Ii L(t) spline through local minima of

Ii Av(t) 1/2 (U(t) L(t)) Ii(t) Ii(t) -

Av(t) i i 1 end

IMFk(t) Ii(t) Residue Residue -

IMFk k k1 end

Huangs Sifting Process.

Residue s(t) I1(t) Residue i 1 k 1

while Residue not equal zero or not monotone

while Ii has non-negligible local

mean U(t) spline through local maxima of

Ii L(t) spline through local minima of

Ii Av(t) 1/2 (U(t) L(t)) Ii(t) Ii(t) -

Av(t) i i 1 end

IMFk(t) Ii(t) Residue Residue -

IMFk k k1 end

Huangs Sifting Process.

Residue s(t) I1(t) Residue i 1 k 1

while Residue not equal zero or not monotone

while Ii has non-negligible local

mean U(t) spline through local maxima of

Ii L(t) spline through local minima of

Ii Av(t) 1/2 (U(t) L(t)) Ii(t) Ii(t) -

Av(t) i i 1 end

IMFk(t) Ii(t) Residue Residue -

IMFk k k1 end

Huangs Sifting Process.

Residue s(t) I1(t) Residue i 1 k 1

while Residue not equal zero or not monotone

while Ii has non-negligible local

mean U(t) spline through local maxima of

Ii L(t) spline through local minima of

Ii Av(t) 1/2 (U(t) L(t)) Ii(t) Ii(t) -

Av(t) (residue--gt) i i 1 end

IMFk(t) Ii(t) Residue

Residue - IMFk k k1

end

Huangs Sifting Process.

Residue s(t) I1(t) Residue i 1 k 1

while Residue not equal zero or not monotone

while Ii has non-negligible local

mean U(t) spline through local maxima of

Ii L(t) spline through local minima of

Ii Av(t) 1/2 (U(t) L(t)) Ii(t) Ii(t) -

Av(t) i i 1 end

IMFk(t) Ii(t) Residue Residue -

IMFk k k1 end

Huangs Sifting Process.

Residue s(t) I1(t) Residue i 1 k 1

while Residue not equal zero or not monotone

while Ii has non-negligible local

mean U(t) spline through local maxima of

Ii L(t) spline through local minima of

Ii Av(t) 1/2 (U(t) L(t)) Ii(t) Ii(t) -

Av(t) i i 1 end

IMFk(t) Ii(t) Residue Residue -

IMFk k k1 end

Huangs Sifting Process.

Residue s(t) I1(t) Residue i 1 k 1

while Residue not equal zero or not monotone

while Ii has non-negligible local

mean U(t) spline through local maxima of

Ii L(t) spline through local minima of

Ii Av(t) 1/2 (U(t) L(t)) Ii(t) Ii(t) -

Av(t) i i 1 end

IMFk(t) Ii(t) Residue Residue -

IMFk k k1 end

Huangs Sifting Process.

Residue s(t) I1(t) Residue i 1 k 1

while Residue not equal zero or not monotone

while Ii has non-negligible local

mean U(t) spline through local maxima of

Ii L(t) spline through local minima of

Ii Av(t) 1/2 (U(t) L(t)) Ii(t) Ii(t) -

Av(t) i i 1 end

IMFk(t) Ii(t) Residue Residue -

IMFk k k1 end

Huangs Sifting Process.

Residue s(t) I1(t) Residue i 1 k 1

while Residue not equal zero or not monotone

while Ii has non-negligible local

mean U(t) spline through local maxima of

Ii L(t) spline through local minima of

Ii Av(t) 1/2 (U(t) L(t)) Ii(t) Ii(t) -

Av(t) i i 1 end

IMFk(t) Ii(t) Residue Residue -

IMFk k k1 end

Huangs Sifting Process.

Residue s(t) I1(t) Residue i 1 k 1

while Residue not equal zero or not monotone

while Ii has non-negligible local

mean U(t) spline through local maxima of

Ii L(t) spline through local minima of

Ii Av(t) 1/2 (U(t) L(t)) Ii(t) Ii(t) -

Av(t) i i 1 end

IMFk(t) Ii(t) Residue Residue -

IMFk k k1 end

Huangs Sifting Process.

Residue s(t) I1(t) Residue i 1 k 1

while Residue not equal zero or not monotone

while Ii has non-negligible local

mean U(t) spline through local maxima of

Ii L(t) spline through local minima of

Ii Av(t) 1/2 (U(t) L(t)) Ii(t) Ii(t) -

Av(t) (residue--gt) i i 1 end

IMFk(t) Ii(t) Residue

Residue - IMFk k k1

end

Huangs Sifting Process.

Residue s(t) I1(t) Residue i 1 k 1

while Residue not equal zero or not monotone

while Ii has non-negligible local

mean U(t) spline through local maxima of

Ii L(t) spline through local minima of

Ii Av(t) 1/2 (U(t) L(t)) Ii(t) Ii(t) -

Av(t) i i 1 end

IMFk(t) Ii(t) Residue Residue -

IMFk k k1 end

Huangs Sifting Process.

Residue s(t) I1(t) Residue i 1 k 1

while Residue not equal zero or not monotone

while Ii has non-negligible local

mean U(t) spline through local maxima of

Ii L(t) spline through local minima of

Ii Av(t) 1/2 (U(t) L(t)) Ii(t) Ii(t) -

Av(t) i i 1 end

IMFk(t) Ii(t) Residue Residue -

IMFk k k1 end

Huangs Sifting Process.

Residue s(t) I1(t) Residue i 1 k 1

while Residue not equal zero or not monotone

while Ii has non-negligible local

mean U(t) spline through local maxima of

Ii L(t) spline through local minima of

Ii Av(t) 1/2 (U(t) L(t)) Ii(t) Ii(t) -

Av(t) i i 1 end

IMFk(t) Ii(t) Residue Residue -

IMFk k k1 end

Huangs Sifting Process.

Residue s(t) I1(t) Residue i 1 k 1

while Residue not equal zero or not monotone

while Ii has non-negligible local

mean U(t) spline through local maxima of

Ii L(t) spline through local minima of

Ii Av(t) 1/2 (U(t) L(t)) Ii(t) Ii(t) -

Av(t) i i 1 end

IMFk(t) Ii(t) Residue Residue -

IMFk k k1 end

Huangs Sifting Process.

Residue s(t) I1(t) Residue i 1 k 1

while Residue not equal zero or not monotone

while Ii has non-negligible local

mean U(t) spline through local maxima of

Ii L(t) spline through local minima of

Ii Av(t) 1/2 (U(t) L(t)) Ii(t) Ii(t) -

Av(t) i i 1 end

IMFk(t) Ii(t) Residue Residue -

IMFk k k1 end

Huangs Sifting Process.

Residue s(t) I1(t) Residue i 1 k 1

while Residue not equal zero or not monotone

while Ii has non-negligible local

mean U(t) spline through local maxima of

Ii L(t) spline through local minima of

Ii Av(t) 1/2 (U(t) L(t)) Ii(t) Ii(t) -

Av(t) i i 1 end

IMFk(t) Ii(t) Residue Residue -

IMFk k k1 end

Huangs Sifting Process.

Residue s(t) I1(t) Residue i 1 k 1

while Residue not equal zero or not monotone

while Ii has non-negligible local

mean U(t) spline through local maxima of

Ii L(t) spline through local minima of

Ii Av(t) 1/2 (U(t) L(t)) Ii(t) Ii(t) -

Av(t) (residue--gt) i i 1 end

IMFk(t) Ii(t) Residue

Residue - IMFk k k1

end

Huangs Sifting Process.

Residue s(t) I1(t) Residue i 1 k 1

while Residue not equal zero or not monotone

while Ii has non-negligible local

mean U(t) spline through local maxima of

Ii L(t) spline through local minima of

Ii Av(t) 1/2 (U(t) L(t)) Ii(t) Ii(t) -

Av(t) i i 1 end

IMFk(t) Ii(t) Residue Residue -

IMFk k k1 end

Huangs Sifting Process.

Residue s(t) I1(t) Residue i 1 k 1

while Residue not equal zero or not monotone

while Ii has non-negligible local

mean U(t) spline through local maxima of

Ii L(t) spline through local minima of

Ii Av(t) 1/2 (U(t) L(t)) Ii(t) Ii(t) -

Av(t) i i 1 end

IMFk(t) Ii(t) Residue Residue -

IMFk k k1 end

Huangs Sifting Process.

Residue s(t) I1(t) Residue i 1 k 1

while Residue not equal zero or not monotone

while Ii has non-negligible local

mean U(t) spline through local maxima of

Ii L(t) spline through local minima of

Ii Av(t) 1/2 (U(t) L(t)) Ii(t) Ii(t) -

Av(t) i i 1 end

IMFk(t) Ii(t) Residue Residue -

IMFk k k1 end

Huangs Sifting Process.

Residue s(t) I1(t) Residue i 1 k 1

while Residue not equal zero or not monotone

while Ii has non-negligible local

mean U(t) spline through local maxima of

Ii L(t) spline through local minima of

Ii Av(t) 1/2 (U(t) L(t)) Ii(t) Ii(t) -

Av(t) i i 1 end

IMFk(t) Ii(t) Residue Residue -

IMFk k k1 end

Huangs Sifting Process.

Residue s(t) I1(t) Residue i 1 k 1

while Residue not equal zero or not monotone

while Ii has non-negligible local

mean U(t) spline through local maxima of

Ii L(t) spline through local minima of

Ii Av(t) 1/2 (U(t) L(t)) Ii(t) Ii(t) -

Av(t) i i 1 end

IMFk(t) Ii(t) Residue Residue -

IMFk k k1 end

Huangs Sifting Process.

Residue s(t) I1(t) Residue i 1 k 1

while Residue not equal zero or not monotone

while Ii has non-negligible local

mean U(t) spline through local maxima of

Ii L(t) spline through local minima of

Ii Av(t) 1/2 (U(t) L(t)) Ii(t) Ii(t) -

Av(t) i i 1 end

IMFk(t) Ii(t) Residue Residue -

IMFk k k1 end

Huangs Sifting Process.

Residue s(t) I1(t) Residue i 1 k 1

while Residue not equal zero or not monotone

while Ii has non-negligible local

mean U(t) spline through local maxima of

Ii L(t) spline through local minima of

Ii Av(t) 1/2 (U(t) L(t)) Ii(t) Ii(t) -

Av(t) i i 1 end

IMFk(t) Ii(t) Residue Residue -

IMFk k k1 end

The Hilbert Spectrum

Now that we have a set of IMFs construct their

analytic counterparts using the Hilbert transform.

Now we can look at f(t) in time and frequency

space simultaneously. Instantaneous frequency is

given by the derivative of the phase angle. The

Hilbert spectrum H(?,t) gives the instantaneous

amplitude (energy) as a function of frequency. A

plot of H(?,t) provides an intuitive, visual

representation of the signal in time and

frequency.

Example EMD for a Test Signal with Known

Components.

Hilbert Spectrum for EMD of Example Signal.

Real-World Example Temperature Data from

Amherst, MA from 1988 - 2005.

(No Transcript)

Physical Interpretation of IMFs for Amherst

Temperature Data.

IMF 8

IMF 5

Huang-Hilbert Spectrum for EMD of Amherst

Temperature Data

-Most of the signal information contained in the

bottom 10 of all frequencies represented -Can we

make a correspondence between what we see here

and the IMFs of the decomposition? -This would

help identify which modes are important in

determining the character of the data.

Hilbert Spectrum for IMF 8

Hilbert Spectrum for IMF 7

Increasing Frequencies Decreasing Information!

IMF 4

IMF 5

IMF 6

(No Transcript)

Interpretation of Residue (Shows Overall Trend).

Bibliography.

- Langton, C. Hilbert Transform, Anayltic Signal

and the Complex Envelope Signal Processing and

Simulation Newsletter, 1999.

http//www.complextoreal.com/tcomplex.htm - Electrocardiograph signal graphic

http//www.neurotraces.com/scilab/scilab2/node61.h

tml - MATLAB code used for all computations, and

sifting graphics - Rilling, G. MATLAB code for computation of

EMD and illustrative slides. - http//perso.ens lyon.fr/patrick.flandrin/emd.h

tml - Temperature data from Amherst, MA
- Williams, C.N., Jr., M.J. Menne, R.S. Vose, and

D.R. Easterling. 2006. United States - Historical Climatology Network Daily

Temperature, Precipitation, and Snow Data. - ORNL/CDIAC-118, NDP-070. Online at

http//cdiac.ornl.gov/epubs/ndp/ushcn/usa.html

site MA190120 - Burning Globe graphic http//river2sea72.wordp

ress.com/2007/03/24/its-not-all-about-carbon-or-is

-it/

About PowerShow.com

You can use PowerShow.com to find and download example online PowerPoint ppt presentations on just about any topic you can imagine so you can learn how to improve your own slides and presentations for free. Or use it to find and download high-quality how-to PowerPoint ppt presentations with illustrated or animated slides that will teach you how to do something new, also for free. Or use it to upload your own PowerPoint slides so you can share them with your teachers, class, students, bosses, employees, customers, potential investors or the world. Or use it to create really cool photo slideshows - with 2D and 3D transitions, animation, and your choice of music - that you can share with your Facebook friends or Google+ circles. That's all free as well!

For a small fee you can get the industry's best online privacy or publicly promote your presentations and slide shows with top rankings. But aside from that it's free. We'll even convert your presentations and slide shows into the universal Flash format with all their original multimedia glory, including animation, 2D and 3D transition effects, embedded music or other audio, or even video embedded in slides. All for free. Most of the presentations and slideshows on PowerShow.com are free to view, many are even free to download. (You can choose whether to allow people to download your original PowerPoint presentations and photo slideshows for a fee or free or not at all.) Check out PowerShow.com today - for FREE. There is truly something for everyone!

presentations for free. Or use it to find and download high-quality how-to PowerPoint ppt presentations with illustrated or animated slides that will teach you how to do something new, also for free. Or use it to upload your own PowerPoint slides so you can share them with your teachers, class, students, bosses, employees, customers, potential investors or the world. Or use it to create really cool photo slideshows - with 2D and 3D transitions, animation, and your choice of music - that you can share with your Facebook friends or Google+ circles. That's all free as well!

For a small fee you can get the industry's best online privacy or publicly promote your presentations and slide shows with top rankings. But aside from that it's free. We'll even convert your presentations and slide shows into the universal Flash format with all their original multimedia glory, including animation, 2D and 3D transition effects, embedded music or other audio, or even video embedded in slides. All for free. Most of the presentations and slideshows on PowerShow.com are free to view, many are even free to download. (You can choose whether to allow people to download your original PowerPoint presentations and photo slideshows for a fee or free or not at all.) Check out PowerShow.com today - for FREE. There is truly something for everyone!

For a small fee you can get the industry's best online privacy or publicly promote your presentations and slide shows with top rankings. But aside from that it's free. We'll even convert your presentations and slide shows into the universal Flash format with all their original multimedia glory, including animation, 2D and 3D transition effects, embedded music or other audio, or even video embedded in slides. All for free. Most of the presentations and slideshows on PowerShow.com are free to view, many are even free to download. (You can choose whether to allow people to download your original PowerPoint presentations and photo slideshows for a fee or free or not at all.) Check out PowerShow.com today - for FREE. There is truly something for everyone!

Recommended

«

/ »

Page of

«

/ »

Promoted Presentations

Related Presentations

Page of

Home About Us Terms and Conditions Privacy Policy Presentation Removal Request Contact Us Send Us Feedback

Copyright 2018 CrystalGraphics, Inc. — All rights Reserved. PowerShow.com is a trademark of CrystalGraphics, Inc.

Copyright 2018 CrystalGraphics, Inc. — All rights Reserved. PowerShow.com is a trademark of CrystalGraphics, Inc.

The PowerPoint PPT presentation: "The Hilbert Transform and Empirical Mode Decomposition:" is the property of its rightful owner.

Do you have PowerPoint slides to share? If so, share your PPT presentation slides online with PowerShow.com. It's FREE!