Choosing dimension for filter

19 vues (au cours des 30 derniers jours)
Aboltabol
Aboltabol le 28 Mar 2023
Commenté : Star Strider le 28 Mar 2023
I have voltage recordings from x sensors for y time-frames, each spanning z seconds. This is stored in a x*y*z matrix. Now, I designed an iirnotch to remove the line noise (60 Hz). When I apply the filter to the matrix using filtfilt, am I correct in assuming that the filter will run along the x dimension, and is this choice of dimension appropriate?

Réponses (1)

Star Strider
Star Strider le 28 Mar 2023
I believe so.
From the filtfilt documentation for x:
The function operates along the first array dimension of x unless x is a row vector. If x is a row vector, then the function operates along the second dimension.
So I believe it will operate along each column (down the rows) of the matrix and across all dimensions. (I have only filtered 2D matrices with it, not 3D.)
.
  2 commentaires
Aboltabol
Aboltabol le 28 Mar 2023
Modifié(e) : Aboltabol le 28 Mar 2023
My matrix has sensors as rows, events as columns, and time as pages. Now, ought not the filter run per sensor per event along the time dimension i.e. the third dimension (page)?
Star Strider
Star Strider le 28 Mar 2023
Yes, it should, however that information did not originally appear.
Use the permute function to put them in the order you want them, for example —
A = randn(12, 5, 30); % Sensor x Event x Time
Ap = permute(A, [3 1 2]);
PermutedSize = size(Ap) % Time x Sensor x Event
PermutedSize = 1×3
30 12 5
That appears to me to be appropriate for the filtering you want to do. If it is not, changing it is straightforward.
.

Connectez-vous pour commenter.

Catégories

En savoir plus sur EEG/MEG/ECoG dans Help Center et File Exchange

Tags

Produits


Version

R2022b

Community Treasure Hunt

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

Start Hunting!

Translated by