13.02.2013 Views

Development of a wavelet-based algorithm to detect and determine ...

Development of a wavelet-based algorithm to detect and determine ...

Development of a wavelet-based algorithm to detect and determine ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

7.1. THEORY 54<br />

bance.<br />

Step 2: Evaluation <strong>of</strong> the square <strong>of</strong> the <strong>wavelet</strong> coefficients found at step 1.<br />

Step 3: Calculation <strong>of</strong> the dis<strong>to</strong>rted signal energy, in each <strong>wavelet</strong> coefficient level.<br />

The "energy" mentioned above is <strong>based</strong> on the Parseval’s theorem: "the energy that<br />

a time domain function contains is equal <strong>to</strong> the sum <strong>of</strong> all energy concentrated in<br />

the different resolution levels <strong>of</strong> the corresponding <strong>wavelet</strong> transformed signal".<br />

This can be mathematically expressed as:<br />

Where:<br />

N�<br />

|f(n)| 2 =<br />

n=1<br />

N�<br />

|aJ(n)| 2 +<br />

n=1<br />

f(n): Time domain signal in study<br />

J�<br />

j=1 n=1<br />

N: Total number <strong>of</strong> samples <strong>of</strong> the signal<br />

N�<br />

|f(n)| 2 : Total energy <strong>of</strong> the f(n) signal<br />

n=1<br />

N�<br />

|dj(n)| 2<br />

N�<br />

|aJ(n)| 2 : Total energy concentrated in the level "j" <strong>of</strong> the approximated<br />

n=1<br />

J�<br />

j=1 n=1<br />

version <strong>of</strong> the signal.<br />

N�<br />

|dj(n)| 2 :Total energy concentrated in the detailed version <strong>of</strong> the signal,<br />

from levels "1" <strong>to</strong> "j".<br />

(7.1)<br />

Step 4: In this stage the steps 1, 2 <strong>and</strong> 3 are repeated for the corresponding "pure sinu-<br />

soidal version" <strong>of</strong> the signal in study.<br />

Step 5: The <strong>to</strong>tal dis<strong>to</strong>rted signal energy <strong>of</strong> the signal in study (found in step 3) is com-<br />

pared <strong>to</strong> the corresponding one <strong>of</strong> the pure sinusoidal signal version (evaluated<br />

in step 4). The result <strong>of</strong> this comparison is a deviation that can be evaluated by<br />

equation 7.2.<br />

� �<br />

energy(j) − reference_energy(j)<br />

deviation_energy(j)(%) =<br />

· 100% (7.2)<br />

reference_energy(7)<br />

where:

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!