14.02.2013 Views

Mathematics in Independent Component Analysis

Mathematics in Independent Component Analysis

Mathematics in Independent Component Analysis

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

Chapter 15. Neurocomput<strong>in</strong>g, 69:1485-1501, 2006 221<br />

Signal [a.u.] Signal [a.u.]<br />

10 8 6 4 2 0 -2 -4<br />

10 8 6 4 2 0<br />

14<br />

12<br />

10<br />

8<br />

6<br />

4<br />

2<br />

0<br />

-2<br />

-2 -4<br />

14<br />

12<br />

10<br />

8<br />

6<br />

4<br />

2<br />

0<br />

-2<br />

δ [ppm]<br />

10 8 6 4 2 0 -2 0<br />

10 8 6 4 2 0 -2 0<br />

δ [ppm]<br />

Fig. 8. The graph uncovers the differences of the LICA and KPCA denois<strong>in</strong>g algorithms.<br />

As a reference the correspond<strong>in</strong>g 1D slice of the orig<strong>in</strong>al P11 spectrum is<br />

displayed on top. From top to bottom the three curves show: The difference of the<br />

orig<strong>in</strong>al and the spectrum with the GEVD-MP algorithm applied, the difference between<br />

the orig<strong>in</strong>al and the LICA denoised spectrum and the difference between the<br />

orig<strong>in</strong>al and the KPCA denoised spectrum. To compare the graphs <strong>in</strong> one diagram<br />

the three graphs are translated vertically by 2, 4 and 6 respectively.<br />

where<br />

2σ 2 =<br />

1<br />

400 ∗ 399<br />

is the width parameter σ, was chosen.<br />

400<br />

�<br />

i,j=1<br />

� xi − xj � 2<br />

8<br />

7<br />

6<br />

5<br />

4<br />

3<br />

2<br />

1<br />

(24)<br />

F<strong>in</strong>ally the kernel matrix K was expressed <strong>in</strong> terms of its EVD (equation 17)<br />

which lead to the expansion parameters α necessary to determ<strong>in</strong>e the pr<strong>in</strong>cipal<br />

axes of the correspond<strong>in</strong>g feature space Ω (m) :<br />

400<br />

�<br />

ω = αiΦ(xi). (25)<br />

i=1<br />

Similar to the orig<strong>in</strong>al data, the noisy data of the reconstructed spectra were<br />

used to form six 400 × 1024 dimensional pattern matrices P (m) , m = 1, . . . , 6.<br />

24

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

Saved successfully!

Ooh no, something went wrong!