Borrar filtros
Borrar filtros

Signal auto-correlation in matlab (xcorr function)

6 visualizaciones (últimos 30 días)
Panos Kerezoudis
Panos Kerezoudis el 20 de Abr. de 2023
Comentada: Panos Kerezoudis el 21 de Abr. de 2023
Hi!
I am working on an EEG-electrophysiology dataset ('csd'). It is structured as a 3D array (channels x time x trials). The dimensions are 16x1527x200.
I am trying to calculate trial-averaged autocorrelation for each channel and I am getting the error message below. I think it has to do with the last part of my loop? The dimensions of the ac0 output (dot products) are 1x201 for some reason, I think the way the function is built, it adds one more column?
Any help would be great!
  2 comentarios
Mathieu NOE
Mathieu NOE el 20 de Abr. de 2023
if you expect someone to help you , it's good practice to share some code and data
Panos Kerezoudis
Panos Kerezoudis el 20 de Abr. de 2023
Absolutely. The code is below. I cannot share data unfortunately (confidentiality reasons).
% Calculate trial-averaged auto-correlation
nchan = size(csd, 1); % channels
ntrials = size(csd, 3); % trials
nlags = 100; % n of lags
ac = zeros(1, 2*nlags + 1); % initialize auto-correlation vector
for i = 1:nchan % loop over channels
for k = 1:ntrials % loop over trials
x = csd(i, :, k) - mean(csd(i, :, k)); % mean-normalize data
[acO, lags] = xcorr(x, nlags, 'biased'); % biased estimate
ac(i, :) = ac(i, :) + double(acO(i, :)/ntrials); % add to total scaled by 1/k
end
end

Iniciar sesión para comentar.

Respuesta aceptada

Chunru
Chunru el 21 de Abr. de 2023
csd = randn(4, 100, 20); % dummy data
% Calculate trial-averaged auto-correlation
nchan = size(csd, 1); % channels
ntrials = size(csd, 3); % trials
nlags = 100; % n of lags
ac = zeros(nchan, 2*nlags + 1); % initialize auto-correlation vector
% ^^^^
for i = 1:nchan % loop over channels
for k = 1:ntrials % loop over trials
x = csd(i, :, k) - mean(csd(i, :, k)); % mean-normalize data
[acO, lags] = xcorr(x, nlags, 'biased'); % biased estimate
ac(i, :) = ac(i, :) + double(acO/ntrials); % add to total scaled by 1/k
% ^^^
end
end

Más respuestas (0)

Categorías

Más información sobre EEG/MEG/ECoG en Help Center y File Exchange.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by