Hi Pseudoscientist,
I understand that you are attempting to apply a bandstop filter to fMRI data represented as a 4D matrix in MATLAB. You want to filter out specific frequency bands in the temporal dimension but your current code does not seem to achieve the desired filtering effect.
Assuming "timeseries" is your 4D fMRI data matrix and "f_sampling" is the sampling frequency of your temporal data, you are using nested loops to apply the "bandstop" filter to each spatial slice of the data. It's important to ensure that the sampling frequency is set correctly to match the temporal resolution of your fMRI data.
Make sure to replace "f_sampling" with the actual sampling frequency of your data. If the filtering still does not work as expected, you may need to verify the design of your filter, ensure the frequency specifications are correct relative to the Nyquist frequency, and check that the "bandstop" function is receiving the correct parameters.
For more information on how to use the `bandstop` function, refer to the documentation of "bandstop filter" - https://www.mathworks.com/help/signal/ref/bandstop.html
To understand more about digital filtering and design methods, you can look at the Filter Design and Analysis Tool documentation - https://www.mathworks.com/help/signal/filter-design.html
Hope this helps!
