Prashant Kumar answered . 2024-05-17 21:55:50

signal = ...;                                           % Signal Vector
tv = ...;                                               % Time Vector
Ts = mean(diff(tv));                                    % Sampling Interval
Fs = 1/Ts;                                              % Sampling Frequency (Hz)
Fn = Fs/2;                                              % Nyquist Frequency (Hz)
Wp = 250/Fn;                                            % Passband Frequency (Normalised)
Ws = 252/Fn;                                            % Stopband Frequency (Normalised)
Rp =   1;                                               % Passband Ripple (dB)
Rs =  50;                                               % Stopband Ripple (dB)
[n,Ws] = cheb2ord(Wp,Ws,Rp,Rs);                         % Filter Order
[z,p,k] = cheby2(n,Rs,Ws);                              % Filter Design
[soslp,glp] = zp2sos(z,p,k);                            % Convert To Second-Order-Section For Stability

figure(3)
freqz(soslp, 2^16, Fs)                                  % Filter Bode Plot

filtered_signal = filtfilt(soslp, glp, signal);         % Filter Signal

Note — Your sampling frequency ‘Fs’ must be at least 550 Hz for this filter to work. It has a passband ripple of 1 dB (that is irrelevant here with a Chebyshev Type II filter) and a stopband attenuation of 50 dB. Change those as necessary.

Not satisfied with the answer ?? ASK NOW