Interpolation on same number of samples after FIR filtration

2 views (last 30 days)
Hi,
I have made FIR filtration of my signal (sampled 25 Hz, 4500 samples) and after filtration that signal has about 200 samples less. The original signal is synced with another signal. And I need these two signals to work together and have same number of samples. What is best solution? How to interpolate filtrated signal to the same number of samples like before?

Accepted Answer

Star Strider
Star Strider on 3 May 2022
Filtering should not change the signal length.
I am not certain exactly what you are doing, how your designed your FIR filter, or what your signal is.
I always suggest using the filtfilt function to do the actual filtering.
  1 Comment
Star Strider
Star Strider on 3 May 2022
That is my filter. Before filtration the signal has 4538 samples and after it has just 4477 samples.
I cannot reproduce that behaviour —
y2 = randn(1,4538);
fs = 25;
fc = 0.65; % cut off frequency
Nfir = 124; % order of filter
N = length(y2);
wc=fc/(fs/2);
fil=fir1(Nfir,wc); % Filter design
filt_y_fir = filter(fil,1,y2); % filtering
filtered_signal_size = size(filt_y_fir)
filtered_signal_size = 1×2
1 4538
filt_y_fir = filtfilt(fil,1,y2); % filtering
filtered_signal_size = size(filt_y_fir)
filtered_signal_size = 1×2
1 4538
figure
freqz(fil,1,2^16,fs)
.

Sign in to comment.

More Answers (1)

Lukas Poviser
Lukas Poviser on 3 May 2022
% y2 is vector with original data
fs = 25;
fc = 0.65; % cut off frequency
Nfir = 124; % order of filter
N = length(y2);
wc=fc/(fs/2);
fil=fir1(Nfir,wc); % Filter design
filt_y_fir = filter(fil,1,y2); % filtering
delay = mean(grpdelay(fil)); % removing time shift
tt = tt(1:end-delay);
sn = y2(1:end-delay);
sf = filt_y_fir;
sf(1:delay) = [];
That is my filter. Before filtration the signal has 4538 samples and after it has just 4477 samples.

Products


Release

R2021b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!