Tag Info

Hot answers tagged ica

6

There a few commercial algorithms that do exactly that (Dolby Prologic, DTS Neo 6, Lexicon Logic 7, Bose Videostage, etc.). If your project has the funds, you can simply try to license one of those. The inner workings of these algorithms are rather complicated and typically a mixture of time domain and frequency domain feature extraction, some steering ...

4

The model ICA uses says that there exist some unknown, statistically independent sources, $s_i$ that are non-normally distributed (their distributions are something other than Gaussian): $$s_i \sim S(\mu_{s_i}, \sigma^2_{s_i})$$ where $S$ is some (possibly) known but non-Gaussian distribution with mean $\mu_{s_i}$ and variance $\sigma^2_{s_i}$. Then it ...

3

Since ICA is a generalization of PCA (Principal Component Analysis) one could use it for intuition about the process. In PCA we create a new coordinate system where we can represent each data sample. So if we have 10 Dimensional data for each sample we have 10 components on the new coordinate system. Classic clustering approach is to use less than 10 ...

3

Some use the Frequency-Time representation as an additional source and try to do what you ask. Some papers on that: Single Channel Source Separation with ICA Based Time-Frequency Decomposition (PDF Source). Musical Source Separation: An Introduction. By the way, at the end, probably Deep Learning will yield the best results for this kind of task. Pay ...

2

In the line C=(mnp*mnp')/3 you're dividing by 3 when you should really be dividing by the number of observations, which is 1000 in your case.

2

A very simple algorithm for stereo signals which are a combination of "panned" mono sources. This allows you to take a stereo source and "listen" to signals panned at a given point, or within a given range of the pan parameter. Caveats: It does not work well if the mix is more than a sum of panned mono sources (stereo reverb applied to the mix, tracks ...

2

In general, the topographic map describes the distribution of electrical activity across the brain as measured in a number of well known "stations" (electrodes whose position is known) on the scalp. This electrical activity is the combined result of millions of neurons firing below the skull. The strength of the electrical signal they create is measured ...

2

The matrices $\hat{\mathbf{M}}$ and $\tilde{\mathbf{M}}$ are constructed in such a way that the relation $\mathbf{M}\mathbf{x}=\mathbf{y}$ implies $\hat{\mathbf{M}}\hat{\mathbf{x}}=\hat{\mathbf{y}}$ and $\tilde{\mathbf{M}}\tilde{\mathbf{x}}=\tilde{\mathbf{y}}$. Consequently, for constructing the matrix $\tilde{\mathbf{M}}$, each element $m_{kl}$ of $\... 2 In the explanation above, the last statement is not the actual reason why in ICA you assume non-Gaussian sources. After all, you could as well assume Gaussian sources and after the mixing process, the observed data will look Gaussian, which is not a problem at all. The reason for assuming non-Gaussian sources was pointed out by Comon (1994), and it is ... 2 The problem with the first point is that if we cannot restrict the set of possible answers, then we will not be able to identify the ground truth components. So we will have a solution that is as useful as the others, but we cannot guarantee that what we found are the components that generated the data. E.g. if you have the mix of two people talking, you ... 2 Your model$\mathbf{x}_n = \mathbf{s}_n + \mathbf{w}_n$seems too simplistic. It basically says that your output is just some input corrupted by noise. (Unless$\mathbf{s}_n $is not really your input but some transformed version of it.) Usually, it's more complicated than that in physical systems. This is why a better model would be$\mathbf{x}= \mathbf{...

2

This appears similar to a classic least squared solution of an overdetermined equation that proceeds as follows: Starting with: $$\mathbf{x} = \mathbf{A} \mathbf{s}$$ $\mathbf{A}$ is not a square matrix if overdetermined (more equations than unknownns) so therefore an inverse does not exist. What you do then is multiply both sides by the transpose of $\... 2 Yes and no. The original ICA methods were not made for the case where there is noise on the signal. However, there exist relatively straightforward extensions that work for noisy signals (e.g. [*]). In this case, you can expect some degree of noise suppression but it won't be perfect. That means that ideally, your desired signal components will stick out ... 1 You might consider this book, it builds up ICA from the basics. It even has a description and some references to papers that tackled EEG data with ICA. 1 If you are interested about the state of the art I would suggest you to read about chimes. This was submitted to the Deep Noise Suppression challenge at interspeech 2020. You can also search for the other submissions if you want, you can see the results of the challenged are here. Links accessed at Aug 11th 2021 1 You have $$E[s_ks_k^*]=E[|s_k|^2]=1$$ because the complex random variables$s_k$have zero mean and unit variance. That means that all elements of the main diagonal of$E[\mathbf{s}\mathbf{s}^H]$equal$1$. Furthermore, with$s_k=x_k+jy_kwe have \begin{align}E[s_ks^*_l]&=E[x_kx_l+y_ky_l+j(x_ly_k-x_ky_l)]\\&=E[x_kx_l]+E[y_ky_l]+jE[x_ly_k]-jE[... 1 In a way, you have answered this already with: FastICA suppose that we have as many sources as channels, but does not check it in any way. If I am going to run the algorithm on my data, it would extract two sources and a 2x2 mixing matrix. ICA will indeed provide a separation, along the lines of section 3 ("What is independence?") which would also ... 1 This is an interesting journey, which is detailed further here. I was intrigued by the attempt to derive a Photoplethysmography (PPG) signal just from the visible spectrum as it is sensed by a common web camera. An indicative publication can be found here (link shared in chat). I have to say, I am still a bit skeptical about this and I am thinking that any ... 1 You already know one of the sources here, so you don't need ICA. There are infinitely many couples k and B that verify: A = B + kC,$but in your case if you assume that$B$and$C$are uncorrelated, we can choose the solution that minimizes the correlation between$B$and$C\$. Using the previous equation and the linearity of the correlation ...

1

Removing the mean (and not the sum) from a linear filter leads to a filter that removes DC. For example: N = 10; h = [0.5 1 0.5]; conv(ones(N,1), h, 'valid') leads to 2s. N=10 h = [0.5 1 0.5]; h_high=h-mean(h); conv(ones(N,1), h_high, 'valid') leads to zeros. Alternatively, you can use any high pass filter.

1

The goal of FastICA is to rotate your data (unitary transform) so that each axis looks as non-Gaussian as possible. Gaussian data still looks Gaussian when you rotate it. If you don't "sphere" the data, all the algorithm can really do is rotate the whole block to one axis. By bringing the mean to zero (centering), and normalizing the variance in all ...

1

ICA can be seen as a dimension reduction technique similar to PCA. Although PCA consists of the eigenvectors of the covariance matrix i.e. second order statistics, ICA is not restricted and can contain higher order stat. The goal of ICA is to span the signal space in the most statistical independent vectors. So you might get a considerable dimension ...

1

If you just want to estimate the height of the square shaped component you can pass the composite signal thorough a low pass filter with a cut-off frequency close to DC or less than the minimum frequency you can expect to see in the composite signal apart from the square component. The resulting signal would have most of the non-DC components filtered out ...

1

This probably should be a comment, but I didn't earn that option yet. You can try measuring instrument impulse response. If I understand correctly, your system is essentially one-dimensional. Then if you use thin plate or plate of known thickness instead of thick object/material, your output signal will be an approximation of impulse response. Using it as ...

1

From your question, I take it that you are attempting to analyze, in a completely unsupervised way, unlabelled speech data so as to unveil two groups that should (hopefully) line up well with the "male" and "female" labels. Some notes: Cross-correlation is useless for speech application. The notion of speech similarity (be it at the higher level of the ...

1

For pre-processing I used EEGLAB, that really help, because EEGLAB have ICA tools for extract Independent Components, that's mean you can check how much of your signals are artifacts (eye moved, blink,muscular moved ), and how much really is brain signal.

1

One way is to synchronize Event related potentials(ERP) generator with your EEG measurement device. Alpha(8-13 Hz) potentials can be generated when the test subject closes the eyes and is in a relaxed state. Another ways is to flash black and white patterns at predetermined intervals. Literature can be easily found where many such techniques and related ERPs ...

1

How do I know I am not extracting noise (they are all very similar)? You should first find or record a well understood set of signals with which you can test your feature extraction code to verify your setup is correct and that the extracted components are as expected. Then you'll have higher confidence that you're not extracting noise due to errors in ...

Only top voted, non community-wiki answers of a minimum length are eligible