Hello everyone, I wrote a simple code to analyze an accelerometric signal in the frequency domain and I saw that the magnitudes over the 300 Hz are quite negligible. I used a Chebyshev type 1 lowpass filter and it works:
%%Design a Lowpass IIR Filter
N=7;Fp=300;Ap=1;h=fdesign.lowpass('N,Fp,Ap', N, Fp, Ap, Fs);d=design(h, 'cheby1');%%Apply the filter to to Smooth out the Signal
xfilter = filter(d,x2);%%Overlay the filtered signal on the original signal.
% Filtered signal is delayed
figure;plot(t2,x2,'b',t2,xfilter,'r');grid on;legend({'Original Signal','Filtered Signal'});%set(gcf,'NumberTitle','Off', 'Name','Filtered Signal vs. Actual Signal');
%%Compare the original signal and delay compensated filtered signal
figure;xfiltfilt = filtfilt(d.sosMatrix,d.ScaleValues,x2);plot(t2,x2,t2,xfiltfilt);grid onlegend({'Original Signal','Actual (filtered and delayed signal)'});
Now, I'd like to use a Butterworth lowpass filter: how can I modify my script? Thank you very much!
Best Answer