I have a discusion with someone about noise covariance matrix, he said:
You can specify the noise level in brainstorm by specifying noise covariance matrix (S). If you want to just consider sensor noise (with uniform noise across sensors) then you can specify the noise covariance matrix to be S = sigma^2 I, where I is the Nsensors x Nsensors identity matrix and sigma^2 is the noise variance. For example, for OPMs sigma^2 would be (6 fT)^2, if you consider just unit bandwidth (the noise variance is sigma^2 = (6fT/sqrt(Hz))^2 * 1 Hz).
So as I understand to do it in brainstorm I should do:
- Noise covariane ---> compute from recordings
- Export it to matlab (name it noise_cov)
- Our sensor have a noise variance of: s2 = (6e-15)^2;
Then noise_cov.NoiseCov = s2*eye(size(gain,1)); Is this true?
The noise_cov matrix contain several matrix, I just need to update the NoiseCov?
The other part of his sentence I did not understand it:
, if you consider just unit bandwidth (the noise variance is sigma^2 = (6fT/sqrt(Hz))^2 * 1 Hz).
do you have explanation for it? does what I did above is enough? and what is the application for if we want to consider a unit bandwidth, and how ?
Thanks in advance!