Dear Matlab Community,
I am trying to learn how to apply a bessel filter to force data from a vertical jump, as the force data is similar to an inverted square wave when airborne. After reading the matlab pages on the besself commands, I cannot wrap my head around how to apply the filter to the data, similar to how I would with a butterworth filter. I am hoping that you all can help me figure this part out, or explain the inaccuracies of my understanding (if they exist) so that i can get this sorted.
Here's what my code looks like for a typical lowpass Butterworth filter:
sf = 1000; % sampling rate
cf = 50/(sf/2); % normalize cutoff frequency for low pass filter, 50 Hz cutoff
[b,a] = butter(4,cf,'low'); % determine filter coefficients (obtain coefficients)
Fz = filtfilt(b,a,Fz); % overwrite Fz raw data with filtered data
Here's what I have thus far for the Bessel filter:
sf = sampling rate
cf = 50(sf/2); % normalize cutoff frequency for low pass filter, 50 Hz cutoff
[b,a] = besself(4,cf,'low'); % determine filter coefficients (obtain coefficients)
What is minning now is how to apply the filter to the Fz data. I cannot locate any examples for this... All I can find are examples of plotting the frequency response.
Thanks for your assistance!