Note that this reference documentation is identical to the help that is displayed in MATLAB when you type “help ft_preproc_bandstopfilter”.
FT_PREPROC_BANDSTOPFILTER applies a band-stop filter to the data and thereby removes the spectral components in the specified frequency band Use as [filt] = ft_preproc_bandstopfilter(dat, Fsample, Fbp, N, type, dir, instabilityfix) where dat data matrix (Nchans X Ntime) Fsample sampling frequency in Hz Fbp frequency band, specified as [Fhp Flp] N optional filter order, default is 4 (but) or dependent upon frequency band and data length (fir/firls) type optional filter type, can be 'but' Butterworth IIR filter (default) 'firws' windowed sinc FIR filter 'fir' FIR filter using MATLAB fir1 function 'firls' FIR filter using MATLAB firls function (requires MATLAB Signal Processing Toolbox) 'brickwall' Frequency-domain filter using MATLAB FFT and iFFT function dir optional filter direction, can be 'onepass' forward filter only 'onepass-reverse' reverse filter only, i.e. backward in time 'twopass' zero-phase forward and reverse filter (default except for firws) 'twopass-reverse' zero-phase reverse and forward filter 'twopass-average' average of the twopass and the twopass-reverse 'onepass-zerophase' zero-phase forward filter with delay compensation (default for firws, linear-phase symmetric FIR only) 'onepass-minphase' minimum-phase converted forward filter (non-linear!, firws only) instabilityfix optional method to deal with filter instabilities 'no' only detect and give error (default) 'reduce' reduce the filter order 'split' split the filter in two lower-order filters, apply sequentially df optional transition width (firws) wintype optional window type (firws), can be 'hann' (max passband deviation 0.0063 [0.63%], stopband attenuation -44dB) 'hamming' (default) (max passband deviation 0.0022 [0.22%], stopband attenuation -53dB) 'blackman' (max passband deviation 0.0002 [0.02%], stopband attenuation -74dB) 'kaiser' dev optional max passband deviation/stopband attenuation (firws with kaiser window, default = 0.001 [0.1%, -60 dB]) plotfiltresp optional, 'yes' or 'no', plot filter responses (firws, default = 'no') usefftfilt optional, 'yes' or 'no', use fftfilt instead of filter (firws, default = 'no') Note that a one- or two-pass filter has consequences for the strength of the filter, i.e. a two-pass filter with the same filter order will attenuate the signal twice as strong. Further note that the filter type 'brickwall' filters in the frequency domain, but may have severe issues. For instance, it has the implication that the time domain signal is periodic. Another issue pertains to that frequencies are not well defined over short time intervals; particularly for low frequencies. See also PREPROC