14.02.2013 Views

Mathematics in Independent Component Analysis

Mathematics in Independent Component Analysis

Mathematics in Independent Component Analysis

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

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

Table 1<br />

Number of output signals correlated with noise or source signals after step 1 and<br />

step 2 of the algorithm dAMUSE.<br />

1st step 2nd step<br />

SNR NM Sources Noise Sources Noise Total<br />

20 dB 12 6 0 6 0 6<br />

15 dB 12 5 2 6 1 7<br />

10 dB 12 6 2 7 1 8<br />

5 dB 12 6 3 7 2 9<br />

0 dB 12 7 4 8 3 11<br />

matrix C also <strong>in</strong>creases after the application of the first step (last column of<br />

table 1). As the noise <strong>in</strong>creases, an <strong>in</strong>creas<strong>in</strong>g number of ICs will be available<br />

at the output of the two steps. Comput<strong>in</strong>g, <strong>in</strong> the frequency doma<strong>in</strong>, the correlation<br />

coefficients between the output signals of each step of the algorithm<br />

and noise or source signals we confirm that some are related with the sources<br />

and others with noise. Table 1 (columns 3-6) shows that the maximal correlation<br />

coefficients are distributed between noise and source signals to a vary<strong>in</strong>g<br />

degree. We can see that the number of signals correlated with noise is always<br />

higher <strong>in</strong> the first level. Results show that for low noise levels the first step<br />

(which is ma<strong>in</strong>ly a pr<strong>in</strong>cipal component analysis <strong>in</strong> a space of dimension NM)<br />

achieves good solutions already. However, we can also see (for narrow-band<br />

signals and/or M low) that the time doma<strong>in</strong> characteristics of the signals resemble<br />

the orig<strong>in</strong>al source signals only after a GEVD, i.e. at the output of the<br />

second step rather than with a PCA, i.e. at the output of first step. Figure 5<br />

shows examples of signals that have been obta<strong>in</strong>ed <strong>in</strong> the two steps of the<br />

algorithm for SNR = 10 dB. At the output of the first level the 3 signals with<br />

highest frequency correlation were chosen among the 8 output signals. Us<strong>in</strong>g a<br />

similar criterion to choose 3 signals at the output of the 2nd step (last column<br />

of figure 5), we can see that their time course is more similar to the source<br />

signals than after the first step (middle column of figure 5)<br />

4.3 Denois<strong>in</strong>g of prote<strong>in</strong> NMR spectra<br />

In biophysics the determ<strong>in</strong>ation of the 3D structure of biomolecules like prote<strong>in</strong>s<br />

is of utmost importance. Nuclear magnetic resonance techniques provide<br />

<strong>in</strong>dispensable tools to reach this goal. As hydrogen nuclei are the most abundant<br />

and most sensitive nuclei <strong>in</strong> prote<strong>in</strong>s, proton NMR spectra of prote<strong>in</strong>s<br />

dissolved <strong>in</strong> water are recorded mostly. S<strong>in</strong>ce the concentration of the solvent<br />

is by magnitudes larger than the prote<strong>in</strong> concentration, there is always a large<br />

19

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

Saved successfully!

Ooh no, something went wrong!